[gdal-dev] Using GDAL with gridded XYZ that has no 'nodata' values?

[Even Rouault]

Le mardi 27 mai 2014 00:06:33, Sam Franklin a écrit :
> Hey Even - great tip.  The sort command you suggested worked well with that
> xyz dataset.
> 
> However, I've found a further dataset that does not work with gdal even
> after the sort. This is the error:
> 
> mbp-sf:data08 samfranklin$ gdal_translate data08-02-sorted.xyz
> data08-02.tif
> Input file size is 4371, 11210
> 0...10...20...30...40.ERROR 1: At line 1100207, found 48162.000000 instead
> of 48164.000000 for nBlockYOff = 4951ERROR 1: data08-02-sorted.xyz, band 1:
> IReadBlock failed at X offset 0, Y offset 4951
> 
> ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 4951
> 
> See screenshot below, the data is arranged into three 'strips'. See dropbox
> link to the data |
> https://dl.dropboxusercontent.com/u/4086367/data08-02.xyz.zip
> 
> Any workaround you can think of?

Probably an issue with just sorting on column 2. I didn't have the same error 
as you, but anyway you should try :

sort -nk 2 -nk 1 data08-02.xyz > data08-02-sorted.xyz

with that I manage to translate the file.

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html