Skip to content

Instantly share code, notes, and snippets.

@bertt
Last active June 12, 2023 15:47
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save bertt/1ba1eaff55c99f0deef36c1c245c81d8 to your computer and use it in GitHub Desktop.
Save bertt/1ba1eaff55c99f0deef36c1c245c81d8 to your computer and use it in GitHub Desktop.
gdal_warp_issue
// following works without gdalwarp:
$ wget https://ns_hwh.fundaments.nl/hwh-ahn/AHN3/DTM_5m/M5_65HZ2.zip
$ unzip M5_65HZ2.zip
$ gdalbuildvrt ahn.vrt M5_65HZ2.TIF
$ mkdir tiles
$ docker run -it -v D:/aaa/testdownload:/data tumgis/ctb-quantized-mesh ctb-tile -f Mesh -C -N -o tiles ahn.vrt
// following does not work with gdalwarp:
$ wget https://ns_hwh.fundaments.nl/hwh-ahn/AHN3/DTM_5m/M5_65HZ2.zip
$ unzip M5_65HZ2.zip
$ gdalwarp -s_srs EPSG:7415 -t_srs EPSG:4326+4979 M5_65HZ2.TIF M5_65HZ2_4326.TIF
$ gdalbuildvrt ahn.vrt M5_65HZ2_4326.TIF
$ mkdir tiles
$ docker run -it -v D:/aaa/testdownload:/data tumgis/ctb-quantized-mesh ctb-tile -f Mesh -C -N -o tiles ahn.vrt
Error: The source dataset does not have a spatial reference system assignedError: The source dataset does not have a spatial reference system assigned
But following works (without VRT):
$ docker run -it -v D:/aaa/testdownload:/data tumgis/ctb-quantized-mesh ctb-tile -f Mesh -C -N -o tiles M5_65HZ2_4326.TIF
Fix: set the srs somewhere? how?
Attempts:
- $ gdal_edit -a_srs EPSG:4326+4979 ahn.vrt -> same error
- $ gdal_translate -a_srs EPSG:4326+4979 ahn.vrt ahn1.vrt -> same error
- $ gdalbuildvrt ahn.vrt -a_srs EPSG:4326+4979 M5_65HZ2_4326.TIF -> same error
Fix?
- $ gdalbuildvrt ahn.vrt -a_srs EPSG:4326 M5_65HZ2_4326.TIF
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment