[gdal-dev] gridding binary data

[Etienne Tourigny]

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