[gdal-dev] gridding binary data

[Even Rouault]

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, nhv at 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 at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev