/data/coastal.tif should be 6.9GB I get weird output from gdal_merge??
gdal_merge.py -o /data/coastal.tif -n -9999 -co BIGTIFF=YES -co COMPRESS=LZW -co BLOCKXSIZE=128 -co BLOCKYSIZE=128 -co NUM_THREADS=ALL_CPUS --config CPL_VSIL_USE_TEMP_FILE_FOR_RANDOM_WRITE YES --config GDAL_SWATH_SIZE 1000000 /vsis3/summer/1/coastal.tif /vsis3/summer/2/coastal.tif & disown -h [1] 507311 ubuntu@ip-172-31-31-16:~$ 0Traceback (most recent call last): File "/usr/bin/gdal_merge.py", line 11, in <module> sys.exit(main(sys.argv)) File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 530, in main fi.copy_into(t_fh, band, band, nodata, verbose) File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 284, in copy_into return raster_copy(s_fh, sw_xoff, sw_yoff, sw_xsize, sw_ysize, s_band, File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 56, in raster_copy return raster_copy_with_nodata( File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 96, in raster_copy_with_nodata data_src = s_band.ReadAsArray(s_xoff, s_yoff, s_xsize, s_ysize, File "/usr/lib/python3/dist-packages/osgeo/gdal.py", line 3653, in ReadAsArray return gdal_array.BandReadAsArray(self, xoff, yoff, File "/usr/lib/python3/dist-packages/osgeo/gdal_array.py", line 422, in BandReadAsArray buf_obj = numpy.empty([buf_ysize, buf_xsize], dtype=typecode) numpy.core._exceptions._ArrayMemoryError: Unable to allocate 377. GiB for an array with shape (225000, 450000) and data type float32 -- Regards, Clive Swan --
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev