Le mercredi 23 janvier 2013 17:57:37, David Hoese a écrit : > 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.
gdal_grid takes a OGR layer as input, so the more simple solution is to output a CSV for example. Or use the OGR API to create the datasource in any format you prefer (PostGIS, Shapefile, CSV, etc etc....) > > 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 _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
