Opened 13 years ago
Closed 12 years ago
#1864 closed defect (invalid)
ST_Intersection returns wrong values with high latitudes
| Reported by: | mhiper3 | Owned by: | pramsey |
|---|---|---|---|
| Priority: | medium | Milestone: | PostGIS 1.5.4 |
| Component: | postgis | Version: | 1.5.X |
| Keywords: | Cc: |
Description
When using latitudes far from the equator the ST_Intersects() function returns unexpected values. These are some sample queries that show the inconsistent behavior:
testdb=# select st_intersects(st_geographyfromtext('SRID=4326;POINT(-80 65)'), st_geographyfromtext('SRID=4326;POLYGON((-100 67, -30 67, -30 62, -100 62, -100 67))'));
st_intersects
f
(1 row)
testdb=# select st_intersects(st_geographyfromtext('SRID=4326;POINT(-80 65)'), st_geographyfromtext('SRID=4326;POLYGON((-100 67, -40 67, -40 62, -100 62, -100 67))'));
st_intersects
t
(1 row)
testdb=# select st_intersects(st_geographyfromtext('SRID=4326;POINT(-80 15)'), st_geographyfromtext('SRID=4326;POLYGON((-100 17, -30 17, -30 12, -100 12, -100 17))'));
st_intersects
t
(1 row)
Please review in light of this.
http://blog.opengeo.org/2012/04/30/the-earth-is-not-flat-volume-2/