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

[Moskovitz, Bob@DOC]

I've never seen "1.#IO" for max, but it is interesting that when this geotif was converted to an AAIGrids, it contained "1.#INF".  So this led me to conclude that I'm dealing with infinity.  I'll try what you using "-a_nodata inf" when I get the chance.

-Bob

-----Original Message-----
From: Even Rouault [mailto:even.rouault at mines-paris.org] 
Sent: Wednesday, April 02, 2014 11:33 AM
To: gdal-dev at lists.osgeo.org
Cc: Moskovitz, Bob at DOC
Subject: Re: [gdal-dev] Set max value of "1.#IO" in Geotiff to nodata

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