<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=us-ascii" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.7600.16700"></HEAD>
<BODY>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>Marcello,</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>So you get an error about ptarray_area_spheroid. You
know off hand which version of PostGIS 2.0.0 you are
running.</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>We are running a 2011-01-29 build and your queries return
answers for us. Whether the answer is right is a different
question.</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial> polygon |
gsareageo<BR>---------+------------------<BR> 161109 |
473946981.913427</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial> polygon |
gsareageo<BR>---------+------------------<BR> 161109 |
772770383.761807</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>We are running on 32-bit Windows:</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>PostgreSQL 9.0.2, compiled by Visual C++ build 1500, 32-bit
POSTGIS="2.0.0SVN" GEOS="3.3.0-CAPI-1.7.0" PROJ="Rel. 4.6.1, 21 August 2008"
LIBXML="2.7.6" USE_STATS</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>I thought maybe it would be a difference in Proj, but
geography doesn't use proj for area. It could be a difference in 32-bit
vs. 64-bit. Perhaps someone else running PostGIS 2.0.0 or 1.5 can
test (preferably on 64-bit Linux).</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>Works fine on our 1.5 build as well.</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>FWIW: Your code is more verbose than it needs to
be. You can go straight to geography with
ST_GeogFromText:</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial><A
href="http://www.postgis.org/documentation/manual-svn/ST_GeogFromText.html">http://www.postgis.org/documentation/manual-svn/ST_GeogFromText.html</A></FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>SELECT 161109 polygon
,ST_AREA(ST_GEOGFROMTEXT('POLYGON((15.8327132977611 0.25,16 0.25,16
0,15.8566682147637 0,15.8566112518311 0.00356000009924173,15.8485898971558
0.124770000576973,15.8327132977611 0.25))')) gsareageo;</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>Hope that helps,</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial>Regina and Leo</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial><A
href="http://www.postgis.us">http://www.postgis.us</A></FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV>
<DIV dir=ltr align=left><SPAN class=739464905-04022011><FONT color=#0000ff
size=2 face=Arial></FONT></SPAN> </DIV><BR>
<DIV dir=ltr lang=en-us class=OutlookMessageHeader align=left>
<HR tabIndex=-1>
<FONT size=2 face=Tahoma><B>From:</B>
postgis-users-bounces@postgis.refractions.net
[mailto:postgis-users-bounces@postgis.refractions.net] <B>On Behalf Of
</B>marcello marcellone<BR><B>Sent:</B> Thursday, February 03, 2011 3:52
PM<BR><B>To:</B> postgis-users@postgis.refractions.net<BR><B>Subject:</B>
[postgis-users] ERROR: ptarray_area_spheroid: cannot handle ptarraythat crosses
equator<BR></FONT><BR></DIV>
<DIV></DIV>
<TABLE border=0 cellSpacing=0 cellPadding=0>
<TBODY>
<TR>
<TD vAlign=top>i try to calculate the area for a geographic grid 0.25
degree but for some grids they touche equator i have the error in
object.<BR><BR>"PostgreSQL 9.0.0 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><BR>SELECT
161109 polygon
,ST_AREA(GEOGRAPHY(ST_GEOMETRYFROMTEXT('POLYGON((15.8327132977611 0.25,16
0.25,16 0,15.8566682147637 0,15.8566112518311
0.00356000009924173,15.8485898971558 0.124770000576973,15.8327132977611
0.25))'))) gsareageo<BR><BR>SELECT 161109 polygon
,ST_AREA(GEOGRAPHY(ST_GEOMETRYFROMTEXT('POLYGON((15.75 0,15.75 0.25,16
0.25,16 0,15.75 0))'))) gsareageo;<BR><BR>Any idea about it
?<BR><BR></TD></TR></TBODY></TABLE><BR></BODY></HTML>