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:[email protected]] Sent: Wednesday, April 02, 2014 11:33 AM To: [email protected] Cc: Moskovitz, Bob@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@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 _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
