Hi Norman,
Yes this sounds like it should work. I'll try it when I get a chance.
If I made a geotiff with no navigation(WKT) and had X,Y,Z as individual
bands would that work? How do I tell gdal_grid which dimension or band
is X, which is Y, and which is Z? Thanks.
Could you please CC me in any replies, I get the digest and it's easier
to reply if I get the original email.
-Dave
On 1/23/13 10:30 AM, n...@cape.com wrote:
Dave
Seems to me this is a point not a raster process you are looking for
e.g. How do I make a surface from a bunch of discrete points ?
So instead of having 3 bands of raster you want a stream of XYZ triplets
you can then submit these to gdal_grid or other tools that work with discrete
data
pseudo python assuming you have read your files into numpy arrays
XYZ = array( zip( bandX.flatten(), bandY.flatten(), bandZ.flatten() ) )
convert XYZ to any OGR format
call gdal_grid
HTH
Norman
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
