def test_point_on_surface(self):
# Reference values.
if oracle:
# SELECT SDO_UTIL.TO_WKTGEOMETRY(SDO_GEOM.SDO_POINTONSURFACE(GEOAPP_COUNTRY.MPOLY, 0.05))
# FROM GEOAPP_COUNTRY;
ref = {
'New Zealand': fromstr('POINT (174.616364 -36.100861)',
srid=4326),
'Texas': fromstr('POINT (-103.002434 36.500397)', srid=4326),
}
else:
# Using GEOSGeometry to compute the reference point on surface values
# -- since PostGIS also uses GEOS these should be the same.
ref = {
'New Zealand':
Country.objects.get(name='New Zealand').mpoly.point_on_surface,
'Texas':
Country.objects.get(name='Texas').mpoly.point_on_surface
}
qs = Country.objects.annotate(
point_on_surface=functions.PointOnSurface('mpoly'))
for country in qs:
tol = 0.00001 # SpatiaLite might have WKT-translation-related precision issues
self.assertTrue(ref[country.name].equals_exact(
country.point_on_surface, tol))
def test_point_on_surface(self):
qs = Country.objects.annotate(
point_on_surface=functions.PointOnSurface("mpoly"))
for country in qs:
self.assertTrue(
country.mpoly.intersection(country.point_on_surface))