[gdal-dev] Set max value of "1.#IO" in Geotiff to nodata

[Even Rouault]

Le mercredi 02 avril 2014 02:20:54, Moskovitz, Bob at DOC a écrit :
> Hello gdal,
> 
> 
> 
> I have a set of geotiffs that have a strange max (infinity?) values that I
> would like to set to nodata.  These originally were netcdf files.  With
> the current of gdal in OSGEO4W, gdalinfo reports:
> 
> 
> 
> gdalinfo -mm pola_cascadiad_sift2_maxspeed.tif
> 
> Driver: GTiff/GeoTIFF
> 
> Files: pola_cascadiad_sift2_maxspeed.tif
> 
> Size is 1810, 1100
> 
> Coordinate System is:
> 
> GEOGCS["WGS 84",
> 
>     DATUM["WGS_1984",
> 
>         SPHEROID["WGS 84",6378137,298.257223563,
> 
>             AUTHORITY["EPSG","7030"]],
> 
>         AUTHORITY["EPSG","6326"]],
> 
>     PRIMEM["Greenwich",0],
> 
>     UNIT["degree",0.0174532925199433],
> 
>     AUTHORITY["EPSG","4326"]]
> 
> Origin = (-118.302500000118700,33.779999999966194)
> 
> Pixel Size = (0.000092592592580,-0.000092592592587)
> 
> Metadata:
> 
>   AREA_OR_POINT=Area
> 
> Image Structure Metadata:
> 
>   INTERLEAVE=BAND
> 
> Corner Coordinates:
> 
> Upper Left  (-118.3025000,  33.7800000) (118d18' 9.00"W, 33d46'48.00"N)
> 
> Lower Left  (-118.3025000,  33.6781481) (118d18' 9.00"W, 33d40'41.33"N)
> 
> Upper Right (-118.1349074,  33.7800000) (118d 8' 5.67"W, 33d46'48.00"N)
> 
> Lower Right (-118.1349074,  33.6781481) (118d 8' 5.67"W, 33d40'41.33"N)
> 
> Center      (-118.2187037,  33.7290741) (118d13' 7.33"W, 33d43'44.67"N)
> 
> Band 1 Block=1810x1 Type=Float64, ColorInterp=Gray
> 
>     Computed Min/Max=0.000,1.#IO
> 
> 
> 
> I've been converting the geotiff to AAIGrids and replace all "1.#INF" to my
> nodata value, but I'm hoping that there is an easier way to do this with
> GDAL.

With GDAL 1.10, you should be able to specify "-a_nodata inf".  But I'm not 
completely sure if 1.#IO is the represnetation of infinity. Might be a not-a-
number ?

> 
> 
> 
> Thanks in advance,
> 
> Bob
> 
> 
> 
> Robert Moskovitz
> California Geological Survey
> Seismic Hazards Zonation Program
> 
> 
> CONFIDENTIALITY NOTICE: This communication is intended only for the use of
> the individual or entity to which it is addressed. This message contains
> information from the State of California, California Geological Survey,
> which may be privileged, confidential and exempt from disclosure under
> applicable law, including the Electronic Communications Privacy Act. If
> the reader of this communication is not the intended recipient, you are
> hereby notified that any dissemination, distribution, or copying of this
> communication is strictly prohibited.

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