<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">On 02/04/11 07:02, Paragon Corporation wrote:
<blockquote cite="mid:310476FAE7A54932A4FCFC69C72E63C8@J" type="cite">
  <pre wrap="">Marcello,<br> <br>So you get an error about ptarray_area_spheroid.  You know off hand which<br>version of PostGIS 2.0.0 you are running.<br> <br>We are running a 2011-01-29 build and your queries return answers for us.<br>Whether the answer is right is a different question.<br> <br> polygon |    gsareageo<br>---------+------------------<br>  161109 | 473946981.913427<br> <br> <br> polygon |    gsareageo<br>---------+------------------<br>  161109 | 772770383.761807<br> <br>We are running on 32-bit Windows:<br>PostgreSQL 9.0.2, compiled by Visual C++ build 1500, 32-bit<br>POSTGIS="2.0.0SVN" GEOS="3.3.0-CAPI-1.7.0" PROJ="Rel. 4.6.1, 21 August 2008"<br>LIBXML="2.7.6" USE_STATS<br> <br>I thought maybe it would be a difference in Proj, but geography doesn't use<br>proj for area.  It could be a difference in 32-bit vs. 64-bit.  Perhaps<br>someone else running PostGIS 2.0.0 or 1.5 can test (preferably on 64-bit<br>Linux).<br> <br>Works fine on
 our 1.5 build as well.<br><br>Marcello:<br><br>Ok now i recompiled postgresql and postgis to last snapshot available:<br><br>"PostgreSQL 9.1alpha3 on x86_64-unknown-linux-gnu, compiled by GCC gcc (SUSE Linux) 4.5.0 20100604 [gcc-4_5-branch revision 160292], 64-bit"<br><br>"POSTGIS="2.0.0SVN" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009" LIBXML="2.7.7" USE_STATS"<br><br>Postgis build is 6780 01/29/11<br><br>So i submit this:<br><br>SELECT 161109 polygon ,ST_AREA(ST_GEOGFROMTEXT('POLYGON((15.8327132977611 0.25,16 0.25,16 0,<br>15.8566682147637 0,15.8566112518311 0.00356000009924173,15.8485898971558 0.124770000576973,15.8327132977611 0.25))')) gsareageo<br><br>and i receive the same error:<br><br>ptarray_area_spheroid: cannot handle ptarray that crosses equator<br><br>If i change the points with 0 lat coordinates with a very small quantity as 0.0000000000000000000000000000000001  the select goes ok...<br><br>SELECT 161109 polygon
 ,ST_AREA(ST_GEOGFROMTEXT('POLYGON((15.8327132977611 0.25,16 0.25,16 0.0000000000000000000000000000000001,<br>15.8566682147637 0.0000000000000000000000000000000001,15.8566112518311 0.00356000009924173,15.8485898971558 0.124770000576973,15.8327132977611 0.25))')) gsareageo<br><br>161109;471828022.374555<br><br>but there is a reliable difference between 0 lat coordinate points : <br>473946981.913427 vs.  471828022.374555   <br><br>I installed postgis in separe schema postgis2, now i can create another schema with postgis 1.5 last build but maybe the problem can be in the last geos library or proj ?<br><br><br></pre>
</blockquote></td></tr></table><br>