Hi Everyone,
I was wondering if anyone can shed any light on what I could be doing wrong
here. I am trying to convert a wkt from SRID 900913 to a wkt in 28350 (UTM
50). Everything works fine if I go from 900913 -> 4326 but going from
900913 -> 28350 throws an OGR Exception.
I'd be grateful if anyone could shed any light on what I could be doing
wrong.
Thanks
Python 2.7.3 (default, Sep 26 2012, 21:51:14)
[GCC 4.7.2] on linux2
Type "help", "copyright", "credits" or "license" for more information.
(InteractiveConsole)
>>> from osgeo import ogr
>>> srs = ogr.osr.SpatialReference()
>>> srs.ImportFromEPSG(28350)
0
>>> print srs.ExportToPrettyWkt()
PROJCS["GDA94 / MGA zone 50",
GEOGCS["GDA94",
DATUM["Geocentric_Datum_of_Australia_1994",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6283"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4283"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",117],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",10000000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","28350"]]
>>>
>>> from django.contrib.gis.geos import GEOSGeometry
>>> wkt = 'GEOMETRYCOLLECTION (POLYGON ((1496942.7617285000160336
440277.2828613300225697, 146759.0942870005965233 -371789.7055273400037549,
1086017.2977249994874001 -1291480.0297266000416130,
2514472.4821190014481544 -9783.9396191388004809, 1496942.7617285000160336
440277.2828613300225697)))'
>>> wkt_srid = 900913
>>>
>>> g = GEOSGeometry(wkt, srid=wkt_srid)
>>>
>>> g.transform(4326)
>>> print g.srid
4326
>>>
>>> g = GEOSGeometry(wkt, srid=wkt_srid)
>>> g.transform(28350)
Traceback (most recent call last):
File "<console>", line 1, in <module>
File
"/home/benk/workspace/swcc/VirtualEnv/local/lib/python2.7/site-packages/django/contrib/gis/geos/geometry.py",
line 513, in transform
g.transform(ct)
File
"/home/benk/workspace/swcc/VirtualEnv/local/lib/python2.7/site-packages/django/contrib/gis/gdal/geometries.py",
line 414, in transform
capi.geom_transform_to(self.ptr, sr.ptr)
File
"/home/benk/workspace/swcc/VirtualEnv/local/lib/python2.7/site-packages/django/contrib/gis/gdal/prototypes/errcheck.py",
line 108, in check_errcode
check_err(result)
File
"/home/benk/workspace/swcc/VirtualEnv/local/lib/python2.7/site-packages/django/contrib/gis/gdal/error.py",
line 39, in check_err
raise e(msg)
OGRException: OGR failure.
>>>
>>>
--
You received this message because you are subscribed to the Google Groups
"Django users" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to django-users+unsubscr...@googlegroups.com.
To post to this group, send email to django-users@googlegroups.com.
Visit this group at http://groups.google.com/group/django-users?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
