On Fri, Aug 8, 2008 at 9:00 PM, Michael ODonnell <[EMAIL PROTECTED]> wrote:
> I am trying to run a simple script and the class transformto does not
> appear to be working for me. I have used GDAL before for raster analysis,
> but I have little experience with OSR. Does anyone have any ideas as to why
> this is not working? I have also run a script I found online called
> test_transformation_osr.py and it also produces errors when transforming
> between spatial references. In other words, there seems to be a problem with
> the installation or something was not installed that I should have
> installed. I installed gdal1.5.2 binaries on win32 xp OS system.
>
> Thanks for the help.
> Mike
>
>
> import osr, ogr
>
> # CREATE WKT STRING OF POINT COORDS
> sample = (' 776041.000', ' 3386618.000')
> x = sample[0]
> y = sample[1]
> wkt = 'POINT(%s %s)' % (x, y)
>
> # CREATE PROJECTION OBJECTS
> utm16N = osr.SpatialReference ()
> utm16N.ImportFromEPSG(26916)
> wgs84 = osr.SpatialReference()
> wgs84.ImportFromEPSG(4326)
>
> # CREATE OGR POINT OBJECT, ASSIGN PROJECTION, REPROJECT
> point = ogr.CreateGeometryFromWkt(wkt)
> point.AssignSpatialReference(utm16N)
> point.TransformTo(wgs84)
>
> print wkt
> print point
>
>
>
>
>
>
> _______________________________________________
> gdal-dev mailing list
> [email protected]
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
Hi Mike,
Have you setup GDAL_DATA in the environment variable?
Please go through the following link
http://trac.osgeo.org/gdal/wiki/GdalOgrInPython
Rgds.
Abhay.
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
