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

gdalsrsinfo flanders_bel72.tif

PROJ.4 : '+proj=lcc +lat_1=51.16666723333333 +lat_2=49.8333339
ₒ

+lat_0=90 +lon_0=4.367486666666666 +x_0=150000.013

ₒ +y_0=5400088.438 +ellps=intl

ₒ

+towgs84=-106.869,52.2978,-103.724,0.3366,-0.457,1.8422,-1.2747

ₒ

+units=m +no_defs'

OGC WKT :
PROJCS["Belge 1972 / Belgian Lambert 72",
GEOGCS["Belge 1972",
DATUM["Reseau_National_Belge_1972",
SPHEROID["International 1924",6378388,297.0000000000014,
AUTHORITY["EPSG","7022"]],

TOWGS84[-106.869,52.2978,-103.724,0.3366,-0.457,1.8422,-1.2747],
AUTHORITY["EPSG","6313"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4313"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",51.16666723333333],
PARAMETER["standard_parallel_2",49.8333339],
PARAMETER["latitude_of_origin",90],
PARAMETER["central_meridian",4.367486666666666],
PARAMETER["false_easting",150000.013],
PARAMETER["false_northing",5400088.438],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","31370"]]

The geotransform has been updated to:

• adfGeotransform[0]: 19081.2

• adfGeotransform[1]: 24.9417

• adfGeotransform[2]: 1.92022

• adfGeotransform[3]: 231729

• adfGeotransform[4]: 1.92829

• adfGeotransform[5]: − 24 . 9058

From these coefficients, we can calculate the rotation angle α of the affine

4 ⓦ ). This can also be verified visually by open-
ing both the original and transformed images in an image viewer (e.g., openev or
QGIS).

transformation (approximately −