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
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
