[postgis-devel] [PostGIS] #2006: st_area returns NaN
PostGIS
trac at osgeo.org
Mon Oct 1 00:50:41 PDT 2012
#2006: st_area returns NaN
----------------------+-----------------------------------------------------
Reporter: evberghe | Owner: pramsey
Type: defect | Status: new
Priority: medium | Milestone:
Component: postgis | Version: 1.5.X
Keywords: |
----------------------+-----------------------------------------------------
Description changed by strk:
Old description:
> The following statement returns a NaN for one of the polygons in my table
> (the one going round the Antarctic):
>
> select st_area(geom::geography)/10^6 from hexgrid1;
>
> output of st_astext for the offending geom is
>
> MULTIPOLYGON(((-180 -88.4923488723258,-180 -90,180 -90,180
> -88.4923488723258,162.118373 -88.171111,149.065274 -86.934953,143.650482
> -85.628588,140.735071 -84.296646,138.904131 -82.949353,137.621854
> -81.588793,136.636021 -80.213815,135.80367 -78.821101,135.024934
> -77.405011,130.614732 -76.788369,126.617235 -76.10412,123.007648
> -75.362082,119.754016 -74.570872,116.821712 -73.737902,114.176337
> -72.869454,111.785425 -71.970806,109.619294 -71.046369,107.651354
> -70.099821,105.858101 -69.134226,104.218929 -68.15214,102.71586
> -67.155693,101.33325 -66.14667,100.057503 -65.126564,98.876793
> -64.09663,97.780829 -63.057924,95.269777 -62.928846,92.784259
> -62.755583,90.330873 -62.53894,87.915631 -62.279883,85.543889
> -61.979524,83.220291 -61.639094,80.948747 -61.259921,78.73243
> -60.843404,76.573798 -60.390995,74.474625 -59.904172,72.436054
> -59.384424,70.458654 -58.833232,68.542479 -58.252054,66.687139
> -57.642314,64.891861 -57.005391,63.155551 -56.342614,60.689007
> -56.770884,58.224295 -57.16813,55.747207 -57.526964,53.250484
> -57.842369,50.731408 -58.110721,48.190396 -58.329286,45.630125
> -58.495949,43.054941 -58.609076,40.465984 -58.63646,37.870999
> -58.609778,35.275127 -58.530124,32.682206 -58.399161,30.09412
> -58.219304,27.50986 -57.994024,24.924061 -57.728423,22.324536
> -57.430319,20.658302 -58.179398,18.924185 -58.904744,17.119621
> -59.604908,15.242268 -60.27835,13.290081 -60.923436,11.261404
> -61.538451,9.15508 -62.1216,6.970555 -62.671026,4.708 -63.184822,2.368426
> -63.661056,-0.046198 -64.097787,-2.532853 -64.493106,-5.087385
> -64.84516,-7.704457 -65.15219,-10.377546 -65.412574,-13.098976
> -65.624861,-14.118509 -66.695816,-15.230334 -67.759689,-16.447764
> -68.815497,-17.786591 -69.862034,-19.265623 -70.897816,-20.907348
> -71.921007,-22.738728 -72.929325,-24.792135 -73.919922,-27.10642
> -74.889224,-29.728014 -75.832727,-32.711907 -76.74474,-36.122086
> -77.618066,-40.030693 -78.443626,-44.514624 -79.210078,-49.647612
> -79.903494,-55.485512 -80.507309,-56.821032 -81.917838,-58.460449
> -83.302906,-60.747466 -84.666001,-64.420025 -86.00515,-71.55536
> -87.304963,-90.525593 -88.486268,-151.608241 -89.002398,-180
> -88.4923488723258)))
>
> PostGIS Version: 1.5 USE_GEOS=1 USE_PROJ=1 USE_STATS=1
> PostgreSQL version: PostgreSQL 8.4.13 on x86_64-pc-linux-gnu, compiled by
> GCC gcc-4.4.real (Ubuntu 4.4.3-4ubuntu5.1) 4.4.3, 64-bit
>
> Any attention to this is much appreciated, as is, in fact the whole of
> PostgreSQL/PostGIS.
New description:
The following statement returns a NaN for one of the polygons in my table
(the one going round the Antarctic):
{{{ select st_area(geom::geography)/10^6 from hexgrid1; }}}
output of st_astext for the offending geom is
MULTIPOLYGON(((-180 -88.4923488723258,-180 -90,180 -90,180
-88.4923488723258,162.118373 -88.171111,149.065274 -86.934953,143.650482
-85.628588,140.735071 -84.296646,138.904131 -82.949353,137.621854
-81.588793,136.636021 -80.213815,135.80367 -78.821101,135.024934
-77.405011,130.614732 -76.788369,126.617235 -76.10412,123.007648
-75.362082,119.754016 -74.570872,116.821712 -73.737902,114.176337
-72.869454,111.785425 -71.970806,109.619294 -71.046369,107.651354
-70.099821,105.858101 -69.134226,104.218929 -68.15214,102.71586
-67.155693,101.33325 -66.14667,100.057503 -65.126564,98.876793
-64.09663,97.780829 -63.057924,95.269777 -62.928846,92.784259
-62.755583,90.330873 -62.53894,87.915631 -62.279883,85.543889
-61.979524,83.220291 -61.639094,80.948747 -61.259921,78.73243
-60.843404,76.573798 -60.390995,74.474625 -59.904172,72.436054
-59.384424,70.458654 -58.833232,68.542479 -58.252054,66.687139
-57.642314,64.891861 -57.005391,63.155551 -56.342614,60.689007
-56.770884,58.224295 -57.16813,55.747207 -57.526964,53.250484
-57.842369,50.731408 -58.110721,48.190396 -58.329286,45.630125
-58.495949,43.054941 -58.609076,40.465984 -58.63646,37.870999
-58.609778,35.275127 -58.530124,32.682206 -58.399161,30.09412
-58.219304,27.50986 -57.994024,24.924061 -57.728423,22.324536
-57.430319,20.658302 -58.179398,18.924185 -58.904744,17.119621
-59.604908,15.242268 -60.27835,13.290081 -60.923436,11.261404
-61.538451,9.15508 -62.1216,6.970555 -62.671026,4.708 -63.184822,2.368426
-63.661056,-0.046198 -64.097787,-2.532853 -64.493106,-5.087385
-64.84516,-7.704457 -65.15219,-10.377546 -65.412574,-13.098976
-65.624861,-14.118509 -66.695816,-15.230334 -67.759689,-16.447764
-68.815497,-17.786591 -69.862034,-19.265623 -70.897816,-20.907348
-71.921007,-22.738728 -72.929325,-24.792135 -73.919922,-27.10642
-74.889224,-29.728014 -75.832727,-32.711907 -76.74474,-36.122086
-77.618066,-40.030693 -78.443626,-44.514624 -79.210078,-49.647612
-79.903494,-55.485512 -80.507309,-56.821032 -81.917838,-58.460449
-83.302906,-60.747466 -84.666001,-64.420025 -86.00515,-71.55536
-87.304963,-90.525593 -88.486268,-151.608241 -89.002398,-180
-88.4923488723258)))
PostGIS Version: 1.5 USE_GEOS=1 USE_PROJ=1 USE_STATS=1
PostgreSQL version: PostgreSQL 8.4.13 on x86_64-pc-linux-gnu, compiled by
GCC gcc-4.4.real (Ubuntu 4.4.3-4ubuntu5.1) 4.4.3, 64-bit
Any attention to this is much appreciated, as is, in fact the whole of
PostgreSQL/PostGIS.
--
--
Ticket URL: <http://trac.osgeo.org/postgis/ticket/2006#comment:1>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.
More information about the postgis-devel
mailing list