Image (Re-)projections and Merging - Open Source Geospatial Tools: Applications in Earth Observation - page 124

•

adfGeotransform[4]: 0

• −

25
Because the transform is performed on the input image itself, we create a copy of
the original image. Then, try to update the file with an error threshold of 1pixel:

adfGeotransform[5]:

cp flanders.tif flanders_bel72
gdalmove.py -t_srs EPSG:31370 -et 1 flanders_bel72.tif

___Corner___ ________Original________ _______Adjusted_________
ₒ ______ Err (geo) ______ _Err (pix)_
Upper Left 3799000.000 3168000.000 18997.758 231789.028

ₒ

0.000 0.000 0.00 0.00

Lower Left 3799000.000 3075000.000 26139.940 139142.128

ₒ

0.000 0.000 0.00 0.00

Upper Right 4032000.000 3168000.000 251445.027 249758.435

ₒ

-0.000 0.000 0.00 0.00

Lower Right 4032000.000 3075000.000 258550.648 157052.136

ₒ

36.561 59.399 -1.64 2.26

Center 3915500.000 3121500.000 138806.212 194453.992

ₒ -13.729 -3.711 0.56 -0.11
Maximum check point error is 2.79065 pixels which exceeds the
ₒ error threshold so the file has not been updated.

The reported maximum point error is more than 1pixel (2 . 79065). To update the

file, we set the error threshold to 3pixels:

gdalmove.py -t_srs EPSG:31370 -et 3 flanders_bel72.tif

___Corner___ ________Original________
ₒ

_______Adjusted_________ ______ Err (geo) ______ _Err

ₒ (pix)_
Upper Left 3799000.000 3168000.000 18997.758 231789.028

ₒ

0.000 0.000 0.00 0.00

Lower Left 3799000.000 3075000.000 26139.940 139142.128

ₒ 0.000 0.000 0.00 0.00

Upper Right 4032000.000 3168000.000 251445.027 249758.435

ₒ

-0.000 0.000 0.00 0.00

Lower Right 4032000.000 3075000.000 258550.648 157052.136

ₒ 36.561 59.399 -1.64 2.26

Center 3915500.000 3121500.000 138806.212 194453.992

ₒ

-13.729 -3.711 0.56 -0.11

Updating file...
Done.

This time the file has actually been updated. We can check with gdalsrsinfo
that the coordinate system has indeed been changed.