GeoServer voltea mis GeoTIFFs

Question

Tengo algunos archivos GeoTIFF que quiero renderizar como una capa de mosaico de imágenes en GeoServer. Los archivos fueron creados con gdal_grid y se ven bien en QGIS:

Y se ven bien en GeoServer cuando se cargan como almacenes GeoTIFF individuales. Aquí está el mayor de los dos renderizados en la vista previa de GeoServer OpenLayers:

Pero si intento utilizar un almacén de mosaico de imágenes (utilizando exactamente los mismos archivos, en el mismo directorio), GeoServer los muestra al revés. Cada archivo es volteado dentro de su propia caja (correcta):

¿Hay alguna manera de conseguir que el Mosaico de Imágenes los represente en la orientación adecuada? Preferiría no tener que regenerar las imágenes, pero si hay un ajuste en el gdal_grid línea de comandos que lo haga funcionar, también estaría bien.

Ambas imágenes están en el sistema de coordenadas "EPSG:4326".

---

gdalinfo salida de una de las imágenes:

$ gdalinfo bath-SCI-3sec.tif
Driver: GTiff/GeoTIFF
Files: bath-SCI-3sec.tif
Size is 1801, 601
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (-118.500416666666993,32.499583333333298)
Pixel Size = (0.000833333333334,0.000833333333333)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-118.5004167,  32.4995833) (118d30' 1.50"W, 32d29'58.50"N)
Lower Left  (-118.5004167,  33.0004167) (118d30' 1.50"W, 33d 0' 1.50"N)
Upper Right (-116.9995833,  32.4995833) (116d59'58.50"W, 32d29'58.50"N)
Lower Right (-116.9995833,  33.0004167) (116d59'58.50"W, 33d 0' 1.50"N)
Center      (-117.7500000,  32.7500000) (117d45' 0.00"W, 32d45' 0.00"N)
Band 1 Block=1801x1 Type=Float32, ColorInterp=Gray

---

Línea de comandos utilizada para hacer la imagen:

$ gdal_grid \
    -ot Float32 \
    -outsize 1801 601 \
    -txe -118.500416666667 -116.999583333333 \
    -tye 32.4995833333333 33.0004166666667 \
    -a nearest \
    -a_srs EPSG:4326 \
    bath-SCI-3sec.vrt \
    bath-SCI-3sec.tif

---

Contenido de bath-SCI-3sec.vrt :

<OGRVRTDataSource>
    <OGRVRTLayer name="bath-SCI-3sec">
        <SrcDataSource relativeToVRT="1">bath-SCI-3sec.csv</SrcDataSource>
        <GeometryType>wkbPoint</GeometryType>
        <GeometryField encoding="PointFromColumns" x="field_1" y="field_2" z="field_3" />
    </OGRVRTLayer>
</OGRVRTDataSource>

---

Aquí está el archivo de propiedades (generado por GeoServer) para el mosaico de imágenes:

#-Automagically created from GeoTools-
#Wed Jun 17 10:00:07 CDT 2020
MosaicCRS=EPSG\:4326
Levels=8.333333333351735E-4,8.333333333335202E-4
Heterogeneous=false
Name=bathymetry
HeterogeneousCRS=false
TypeName=bathymetry
Caching=false
ExpandToRGB=false
LocationAttribute=location
PathType=RELATIVE
SuggestedSPI=it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReaderSpi
CheckAuxiliaryMetadata=false
LevelsNum=1

---

Esto es lo que dice GeoServer sobre el CRS de la capa Mosaico de Imágenes:

Answer 1

El mosaico de imágenes no puede manejar imágenes de entrada con ningún tipo de transformación afín, esto incluye tanto la rotación como el volteo.

La situación podría corregirse cambiando el código del mosaico de imágenes.

Si estás interesado en donar algo de tu tiempo y trabajo, te sugiero que te suscribas a la lista de desarrolladores de GeoTools y pidas indicaciones (el código del mosaico de imágenes vive en GeoTools).

Si por el contrario está interesado en financiar la mejora, consulte la página de apoyo comercial en http://geoserver.org/support/ (aviso de que trabajo para una de esas organizaciones).