[gdal-dev] gdal2tiles: RuntimeError: Cannot find coordinate operations from `ENGCRS["WGS_1984_Web_Mercator_Auxiliary_Sphere"

[Terra Frost]

I have two TIFs that I'm able to run gdal2tiles on without issue, however,
when I do gdal_merge I am unable to do gdal2tiles on the result. When I did
gdal_merge on the two TIFs I got this warning:

Warning 1: PROJ: proj_create_from_database: Cannot find proj.db
Warning 1: The definition of projected CRS EPSG:3857 got from GeoTIFF keys
is not the same as the one from the EPSG registry, which may cause issues
during reprojection operations. Set GTIFF_SRS_SOURCE configuration option
to EPSG to use official parameters (overriding the ones from GeoTIFF keys),
or to GEOKEYS to use custom values from GeoTIFF keys and drop the EPSG code.
Warning 1: PROJ: proj_create_from_database: Cannot find proj.db
Warning 1: The definition of projected CRS EPSG:3857 got from GeoTIFF keys
is not the same as the one from the EPSG registry, which may cause issues
during reprojection operations. Set GTIFF_SRS_SOURCE configuration option
to EPSG to use official parameters (overriding the ones from GeoTIFF keys),
or to GEOKEYS to use custom values from GeoTIFF keys and drop the EPSG code.
0...10...20...30...40...50...60...70...80...90...100 - done.


Despite the warning the two TIFs merged without issue, however, when I try
to run gdal2tiles on the merged TIF I get this error:

Traceback (most recent call last):
  File "/usr/bin/gdal2tiles.py", line 16, in <module>
    sys.exit(main(sys.argv))
  File "/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.py", line
3314, in main
    single_threaded_tiling(input_file, output_folder, options)
  File "/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.py", line
3229, in single_threaded_tiling
    conf, tile_details = worker_tile_details(input_file, output_folder,
options)
  File "/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.py", line
3145, in worker_tile_details
    gdal2tiles.open_input()
  File "/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.py", line
1797, in open_input
    self.warped_input_dataset = reproject_dataset(
  File "/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.py", line
987, in reproject_dataset
    to_dataset = gdal.AutoCreateWarpedVRT(from_dataset,
  File "/usr/lib/python3/dist-packages/osgeo/gdal.py", line 3976, in
AutoCreateWarpedVRT
    return _gdal.AutoCreateWarpedVRT(*args)
RuntimeError: Cannot find coordinate operations from
`ENGCRS["WGS_1984_Web_Mercator_Auxiliary_Sphere",EDATUM[""],CS[Cartesian,2],AXIS["easting",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["northing",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
to `EPSG:3857'
free(): invalid pointer
Aborted


When I do gdalinfo on the merged TIF I get this:

Corner Coordinates:
Upper Left  (-10880978.065, 3538629.874)
Lower Left  (-10880978.065, 3537702.774)
Upper Right (-10880302.765, 3538629.874)
Lower Right (-10880302.765, 3537702.774)
Center      (-10880640.415, 3538166.324)

When I do gdalinfo on the two TIFs that I merged I get this:

Upper Left  (-10880811.355, 3538629.874)
Lower Left  (-10880811.355, 3538052.574)
Upper Right (-10880302.755, 3538629.874)
Lower Right (-10880302.755, 3538052.574)
Center      (-10880557.055, 3538341.224)


...and...

Upper Left  (-10880978.065, 3538273.090)
Lower Left  (-10880978.065, 3537702.790)
Upper Right (-10880474.265, 3538273.090)
Lower Right (-10880474.265, 3537702.790)
Center      (-10880726.165, 3537987.940)


Any ideas?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20221109/f827e616/attachment.htm>