In the output the pixel sizes are positive Pixel Size = (10.000000000000000,10.000000000000000).
However, this creates some problems with various gdal tools, e.g., creatinge VRTs or Mosaics:
gdalbuildvrt: Warning 1: gdalbuildvrt does not support positive NS resolution. Skipping <filename>
gdal_merge.py:
Traceback (most recent call last):
File "/usr/bin/gdal_merge.py", line 14, in <module>
sys.exit(main(sys.argv))
^^^^^^^^^^^^^^
File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 650, in main
return gdal_merge(argv)
^^^^^^^^^^^^^^^^
File "/usr/lib/python3/dist-packages/osgeo_utils/gdal_merge.py", line 572, in gdal_merge
t_fh = driver.Create(out_file, xsize, ysize, bands, band_type, create_options)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/usr/lib/python3/dist-packages/osgeo/gdal.py", line 3742, in Create
return _gdal.Driver_Create(self, *args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
RuntimeError: Attempt to create 10000x-17999 dataset is illegal,sizes must be larger than zero.
Full GDAL-Info on an output file
Driver: GTiff/GeoTIFF
Files: masked_curve_id69785.tiff
Size is 437, 386
Coordinate System is:
PROJCRS["ETRS89-extended / LAEA Europe",
BASEGEOGCRS["ETRS89",
ENSEMBLE["European Terrestrial Reference System 1989 ensemble",
MEMBER["European Terrestrial Reference Frame 1989"],
MEMBER["European Terrestrial Reference Frame 1990"],
MEMBER["European Terrestrial Reference Frame 1991"],
MEMBER["European Terrestrial Reference Frame 1992"],
MEMBER["European Terrestrial Reference Frame 1993"],
MEMBER["European Terrestrial Reference Frame 1994"],
MEMBER["European Terrestrial Reference Frame 1996"],
MEMBER["European Terrestrial Reference Frame 1997"],
MEMBER["European Terrestrial Reference Frame 2000"],
MEMBER["European Terrestrial Reference Frame 2005"],
MEMBER["European Terrestrial Reference Frame 2014"],
MEMBER["European Terrestrial Reference Frame 2020"],
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]],
ENSEMBLEACCURACY[0.1]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4258]],
CONVERSION["Europe Equal Area 2001",
METHOD["Lambert Azimuthal Equal Area",
ID["EPSG",9820]],
PARAMETER["Latitude of natural origin",52,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",10,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["False easting",4321000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",3210000,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["northing (Y)",north,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["easting (X)",east,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["Statistical analysis."],
AREA["Europe - European Union (EU) countries and candidates. Europe - onshore and offshore: Albania; Andorra; Austria; Belgium; Bosnia and Herzegovina; Bulgaria; Croatia; Cyprus; Czechia; Denmark; Estonia; Faroe Islands; Finland; France; Germany; Gibraltar; Greece; Hungary; Iceland; Ireland; Italy; Kosovo; Latvia; Liechtenstein; Lithuania; Luxembourg; Malta; Monaco; Montenegro; Netherlands; North Macedonia; Norway including Svalbard and Jan Mayen; Poland; Portugal including Madeira and Azores; Romania; San Marino; Serbia; Slovakia; Slovenia; Spain including Canary Islands; Sweden; Switzerland; Türkiye (Turkey); United Kingdom (UK) including Channel Islands and Isle of Man; Vatican City State."],
BBOX[24.6,-35.58,84.73,44.83]],
ID["EPSG",3035]]
Data axis to CRS axis mapping: 2,1
Origin = (4549180.000000000000000,2746020.000000000000000)
Pixel Size = (10.000000000000000,10.000000000000000)
Metadata:
value_type=continuous
AREA_OR_POINT=Area
Image Structure Metadata:
LAYOUT=COG
COMPRESSION=LZW
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 4549180.000, 2746020.000) ( 13d 2'39.28"E, 47d47' 7.87"N)
Lower Left ( 4549180.000, 2749880.000) ( 13d 2'46.70"E, 47d49'12.80"N)
Upper Right ( 4553550.000, 2746020.000) ( 13d 6' 8.95"E, 47d47' 1.98"N)
Lower Right ( 4553550.000, 2749880.000) ( 13d 6'16.51"E, 47d49' 6.91"N)
Center ( 4551365.000, 2747950.000) ( 13d 4'27.86"E, 47d48' 7.40"N)
Band 1 Block=512x512 Type=Float64, ColorInterp=Gray
Description = masked_curve
This is not the case in the input images, where the origins are different and the pixel sizes are negative. Example: Full GDAL-Info on a SIAM-Processed file (notice the negative pixel size):
Description
In the output the pixel sizes are positive
Pixel Size = (10.000000000000000,10.000000000000000)
.However, this creates some problems with various gdal tools, e.g., creatinge VRTs or Mosaics:
gdalbuildvrt:
Warning 1: gdalbuildvrt does not support positive NS resolution. Skipping <filename>
gdal_merge.py:Full GDAL-Info on an output file
This is not the case in the input images, where the origins are different and the pixel sizes are negative. Example: Full GDAL-Info on a SIAM-Processed file (notice the negative pixel size):
Additional context
This is not really a but since it is not wrong per se, but the input/output should be consistent in my opinion.
References https://github.com/Sen2Cube-at/sen2cube.at/issues/848