<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>
<div> </div>
<div>I can't tell you how much thought there was behind the change from "%.<DECIMAL_PRECISION>f" to "%.<DECIMAL_PRECISION>g". At least noone bothered to change the documentation which continues to tell us: "You can alternatively specify the number of decimal places with the DECIMAL_PRECISION creation option." [http://www.gdal.org/frmt_various.html#AAIGrid]</div>
<div> </div>
<div>"DECIMAL_PRECISION" actually became a misnomer now and should perhaps be called "SIGNIFICANT_DIGITS" -- or the format should be changed back to "%<precision>f". Maybe some of the GDAL devs has an opinion about the issue.</div>
<div> </div>
<div>Hermann</div>
<div>
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b> Montag, 02. September 2013 um 13:30 Uhr<br/>
<b>Von:</b> "Casper Børgesen (CABO)" <CABO@NIRAS.DK><br/>
<b>An:</b> peifer <peifer@gmx.eu>, "gdal-dev@lists.osgeo.org" <gdal-dev@lists.osgeo.org><br/>
<b>Betreff:</b> RE: [gdal-dev] Change of DECIMAL_PRECISION in AAIGrid</div>
<div name="quoted-content">Hi Peifer.<br/>
<br/>
Thank you for your comments. They do explain more in depth what is happening. But as seen in your first example, the value 333.123456 is reduced to 300 while the value 0.00012345 is reduced to 0.0001. I think this is a big change in the relative values of each number.<br/>
<br/>
I am curious to know if the "new" implementation in GDAL 1.10 is actually what is wanted/expected?<br/>
<br/>
I use the DECIMAL_PRECISION for rounding purposes to control how many digits I want to the right of the comma, not how many digits I want to use for the entire number.<br/>
I know I get a lot of zeros if I use DECIMAL_PRECISION=10 in GDAL 1.9.2 but I retain my precision.<br/>
<br/>
In GDAL 1.10 I have to determine the number of digits to the left of the comma used by the largest number in my data set, add the number of digits I want to have to the right of the comma and use that value as DECIMAL_PRECISION. The result must be that I get a lot of (waste) digits to the right of the comma for small numbers and the expected result for large numbers. But isn't this similar to what the patch (mentioned ealier) wanted to avoid? :<br/>
<br/>
> The AAIGrid float format string is "%6.20g"<br/>
> Firstly, the 6 doesn't really add anything, and when writing small integers it causes them to be space-padded, which is a waste of bytes.<br/>
> At present AAIGrid supports "-co DECIMAL_PRECISION=10" which changes the formatstring to "%.10f".<br/>
> This is a little inconsistent, and the %f causes integers to be written as "9.000000000" which is a horrible waste of bytes.<br/>
<br/>
I could write my own GeoTIFF to AAIGrid exporter, but I don't believe that should be necessary :)<br/>
<br/>
<br/>
Kind regards, Casper<br/>
<br/>
<br/>
-----Original Message-----<br/>
From: gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev-bounces@lists.osgeo.org] On Behalf Of peifer<br/>
Sent: 2. september 2013 12:21<br/>
To: gdal-dev@lists.osgeo.org<br/>
Subject: Re: [gdal-dev] Change of DECIMAL_PRECISION in AAIGrid<br/>
<br/>
Casper Børgesen wrote<br/>
> This is my syntax:<br/>
><br/>
> gdal_translate -of AAIGrid -ot Float32 -co DECIMAL_PRECISION=2<br/>
> my_source.vrt my_target.asc<br/>
><br/>
> I have tried the above syntax using gdal 1.9.2 and 1.10 and these are<br/>
> examples from my results:<br/>
><br/>
> 1.9.2:<br/>
> ...<br/>
> 3.77 3.83 3.89 3.87 3.79 3.49 3.03 2.50 1.96 1.56 1.29 0.95 0.77 0.54<br/>
> 0.53<br/>
> 0.41<br/>
> ...<br/>
><br/>
> 1.10:<br/>
> ...<br/>
> 3.8 3.8 3.9 3.9 3.8 3.5 3 2.5 2 1.6 1.3 0.95 0.77 0.54 0.53 0.41 ...<br/>
><br/>
> Has the meaning of DECIMAL_PRECISION changed from 1.9.2 to 1.10?<br/>
><br/>
> It seems like in 1.9.2 the meaning is X characters to the right of the<br/>
> comma and in 1.10 it means X describing characters?<br/>
><br/>
> I haven't tried the 1.10.1 release because I can't find the binaries<br/>
> on Tamas gisinternals.<br/>
<br/>
Hmm. After a look into the source code [0] and some testing [1], it looks to me that DECIMAL_PRECISION is used as precision of a "%.<precision>g" format.<br/>
<br/>
Format letter g means: "Print a number in either scientific notation or in floating-point notation, whichever uses fewer characters."<br/>
<br/>
For format letter g, the precision value defines: "Maximum number of<br/>
*significant* digits."<br/>
<br/>
I see there has been this change recently:<br/>
<a href="http://trac.osgeo.org/gdal/changeset/25032/trunk/gdal/frmts/aaigrid/aaigriddataset.cpp" target="_blank">http://trac.osgeo.org/gdal/changeset/25032/trunk/gdal/frmts/aaigrid/aaigriddataset.cpp</a><br/>
<br/>
Hermann<br/>
<br/>
[0]<br/>
<a href="http://trac.osgeo.org/gdal/browser/trunk/gdal/frmts/aaigrid/aaigriddataset.cpp#L1051" target="_blank">http://trac.osgeo.org/gdal/browser/trunk/gdal/frmts/aaigrid/aaigriddataset.cpp#L1051</a><br/>
<br/>
[1]<br/>
<br/>
# Input file<br/>
<br/>
$ cat in.asc<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1<br/>
yllcorner 1<br/>
cellsize 1<br/>
0.00012345<br/>
1.12345678<br/>
22.1234567<br/>
333.123456<br/>
<br/>
# Test run<br/>
<br/>
$ for i in 1 2 3 4 5 ; do<br/>
echo DECIMAL_PRECISION=$i &&<br/>
gdal_translate -of AAIGrid -ot Float32 -co DECIMAL_PRECISION=$i in.asc out.asc &&<br/>
cat out.asc<br/>
done<br/>
<br/>
# Output<br/>
<br/>
DECIMAL_PRECISION=1<br/>
<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1.000000000000<br/>
yllcorner 1.000000000000<br/>
cellsize 1.000000000000<br/>
0.0001<br/>
1<br/>
2e+01<br/>
3e+02<br/>
<br/>
DECIMAL_PRECISION=2<br/>
<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1.000000000000<br/>
yllcorner 1.000000000000<br/>
cellsize 1.000000000000<br/>
0.00012<br/>
1.1<br/>
22<br/>
3.3e+02<br/>
<br/>
DECIMAL_PRECISION=3<br/>
<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1.000000000000<br/>
yllcorner 1.000000000000<br/>
cellsize 1.000000000000<br/>
0.000123<br/>
1.12<br/>
22.1<br/>
333<br/>
<br/>
DECIMAL_PRECISION=4<br/>
<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1.000000000000<br/>
yllcorner 1.000000000000<br/>
cellsize 1.000000000000<br/>
0.0001234<br/>
1.123<br/>
22.12<br/>
333.1<br/>
<br/>
DECIMAL_PRECISION=5<br/>
<br/>
ncols 1<br/>
nrows 4<br/>
xllcorner 1.000000000000<br/>
yllcorner 1.000000000000<br/>
cellsize 1.000000000000<br/>
0.00012345<br/>
1.1235<br/>
22.123<br/>
333.12<br/>
<br/>
<br/>
<br/>
<br/>
--<br/>
View this message in context: <a href="http://osgeo-org.1560.x6.nabble.com/gdal-dev-Change-of-DECIMAL-PRECISION-in-AAIGrid-tp5075524p5075585.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/gdal-dev-Change-of-DECIMAL-PRECISION-in-AAIGrid-tp5075524p5075585.html</a><br/>
Sent from the GDAL - Dev mailing list archive at Nabble.com.<br/>
_______________________________________________<br/>
gdal-dev mailing list<br/>
gdal-dev@lists.osgeo.org<br/>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></div>
</div>
</div>
</div></div></body></html>