[gdal-dev] gridding binary data
David Hoese
dhoese at gmail.com
Wed Jan 23 10:43:08 PST 2013
Ah sorry about that. I thought OGR formats and GDAL Raster Formats were
the same thing. Thank you.
On 1/23/13 12:12 PM, Etienne Tourigny wrote:
> 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
More information about the gdal-dev
mailing list