Thanks for fixing the precision loss so quickly. "you now can set the configuration option AAIGRID_DATATYPE to Float64 to force Float64 being used to avoid the loss of precision."
I am not exactly sure whether this option relates to the gdal_translate command line utility or the Python binding (ReadAsArray) - or both. Could you perhaps post the updated script that produced the double precision values correctly? Thanks Ole Nielsen ----- Ole Nielsen Australia - Indonesia Facility for Disaster Reduction -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/gdal-dev-How-to-retain-double-precision-in-GDAL-tp5379507p5387617.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev