[gdal-dev] gridding binary data

[David Hoese]

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, 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