Hogendoorn, Rene wrote:
x = numpy.array( [p[0] for p in polygon] )
y = numpy.array( [p[1] for p in polygon] )
z = numpy.array( [p[2] for p in polygon] )
print proj.TransformPoints( len(polygon), x, y, z)
returns:
Traceback (most recent call last):
File "<stdin>", line 37, in <module>
File
"/var/tmp/libgdal0-1.6.0-build/usr/lib64/python2.6/site-packages/osgeo/osr.py",
line 622, in TransformPoints
TypeError: in method 'CoordinateTransformation_TransformPoints',
argument 3 of type 'double *'
I also tried straight python lists of doubles.
So what type of arguments does TransformPoints expect (as there is
obviously nothing like double* in python)?
Rene,
To the best of my knowledge TransformPoints() is not currently working smoothly
with Python. I have just been using TransformPoint(). You might want to file
a ticket on this.
Best regards,
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, [email protected]
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev