Opened 14 years ago
Closed 14 years ago
#450 closed defect (fixed)
ST_Area(geography) returns wrong result for objects crossing dateline
| Reported by: | pramsey | Owned by: | pramsey |
|---|---|---|---|
| Priority: | high | Milestone: | PostGIS 1.5.1 |
| Component: | postgis | Version: | 1.5.X |
| Keywords: | Cc: | noss@… |
Description (last modified by )
http://postgis.refractions.net/pipermail/postgis-users/2010-February/025984.html
A box that doesn't cross the dateline returns a good result:
g4wd=# select st_area(g)/1000000, st_area(transform(geometry(g), 4326)),
astext(geometry(g)) from (select 'POLYGON((178.5 2, 178.5 1, 179.5 1, 179.5
2, 178.5 2))'::geography) as foo(g);
?column? | st_area | astext
-----------------+---------+----------------------------------------------------
12304.814950073 | 1 | POLYGON((178.5 2,178.5 1,179.5 1,179.5 2,178.5
2))
(1 row)
But the same box s hifted over the dateline returns a bad result.
g4wd=# select st_area(g)/1000000, st_area(transform(geometry(g), 4326)),
astext(geometry(g)) from (select 'POLYGON((179.5 2, 179.5 1, -179.5 1,
-179.5 2, 179.5 2))'::geography) as foo(g);
?column? | st_area | astext
------------------+---------+------------------------------------------------------
55910.8880691158 | 359 | POLYGON((179.5 2,179.5 1,-179.5 1,-179.5
2,179.5 2))
(1 row)
Attachments (1)
Change History (5)
comment:1 by , 14 years ago
| Description: | modified (diff) |
|---|
by , 14 years ago
| Attachment: | 450-1.patch added |
|---|
comment:2 by , 14 years ago
| Cc: | added |
|---|
The attached patch fixes the spheroid case. I have found the sphere case ( st_area(geography, false) )also has slightly different results on shifted polygons, even with this fix.
comment:4 by , 14 years ago
| Resolution: | → fixed |
|---|---|
| Status: | new → closed |
Working the sphere area problem in a new ticket (#451) and closing this one out.
Note:
See TracTickets
for help on using tickets.

Fix for spheroid area