gdal_grid takes a vector file (e.g. shapefile) as input, it cannot take a raster (gtiff) as input
please read the documentation before asking other questions http://www.gdal.org/gdal_grid.html Etienne On Wed, Jan 23, 2013 at 2:57 PM, David Hoese <dho...@gmail.com> wrote: > 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 _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
