Raster projection retrieval is confusing?

Open ansonl opened this issue 6 months ago • 2 comments

TouchTerrain's projection retrieved from the raster seems wrong starting sometime on or after the change at https://github.com/ChHarding/TouchTerrain_for_CAGEO/commit/d9cb1ce299c8525d8aa3778c104b8cac7adf32ad

I haven't looked closely at the code changes for the projection retrieval yet but it might be a change in the get_GDAL_projection_and_datum() function that was added.

The previous commit at https://github.com/ChHarding/TouchTerrain_for_CAGEO/commit/a779fb188038a271dd9ff435a3443c032f5e501d has the correct results.

Expected output (from the commit https://github.com/ChHarding/TouchTerrain_for_CAGEO/commit/a779fb188038a271dd9ff435a3443c032f5e501d before the linked commit above): projection: Asia_Lambert_Conformal_Conic

Actual output: DEM projection & datum: Lambert_Conformal_Conic_2SP WGS_1984

I attached a test file and config test.zip

Jun 30 '25 07:06 ansonl