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

[Even Rouault]

Terra,

you didn't specify how you got GDAL and PROJ, and which versions of 
those you use.

Make sure to solve the "proj_create_from_database: Cannot find proj.db" 
error. Check if the PROJ_LIB environment is set and if it points to the 
directory where proj.db is placed

The "free(): invalid pointer" error could also suggest that you have 
several PROJ versions that are linked at runtime by GDAL.

Even

Le 09/11/2022 à 17:54, Terra Frost a écrit :
> 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?
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20221109/3cee6ea3/attachment-0001.htm>