Missing support for projection used in geotiff

[TAlonglong]

In line https://github.com/metno/diana/blob/master/src/miRaster/satgeotiff.cc#L368 and following the projection information from a geotiff file is read.

I have a file here ( this might disappear, but all *.tif files in this directory should have the same behaviour:
/lustre/storeB/project/metproduction/products/satdata_polar/production/antenna-viirs-mitiff/natural_color_sun_20220630_090501.tif

And Diana reports:

2022-06-30 09:41:16,776 ERROR metno.GeoTiff.head_diana: Projection CT 9815 not yet supported
2022-06-30 09:41:16,776 ERROR diana.GEOtiff.readGEOtiffHeader: GEOTIFF_head_diana returned false:/lustre/storeB/project/metproduction/products/satdata_polar/production/antenna-viirs-mitiff/natural_color_sun_20220630_090501.tif
2022-06-30 09:41:16,776 ERROR metno.GeoTiff.head_diana: Projection CT 9815 not yet supported
2022-06-30 09:41:16,776 ERROR diana.GEOtiff: GEOTIFF_read_diana returned false:/lustre/storeB/project/metproduction/products/satdata_polar/production/antenna-viirs-mitiff/natural_color_sun_20220630_090501.tif
2022-06-30 09:41:16,776 ERROR diana.SatManager.setData: Failed readSatFile

What I want is DIANA to read the proj4 information out of this geotiff file.

gdalsrsinfo /lustre/storeB/project/metproduction/products/satdata_polar/production/antenna-viirs-mitiff/natural_color_sun_20220630_090501.tif

PROJ.4 : '+proj=omerc +lat_0=0 +lonc=65.01172638348839 +alpha=-8.6347579930128 +k=1 +x_0=0 +y_0=0 +gamma=0 +a=6370997 +b=6370997 +units=m +no_defs '

OGC WKT :
PROJCS["unknown",
    GEOGCS["unknown",
        DATUM["Unknown_based_on_Normal_Sphere_r_6370997_ellipsoid",
            SPHEROID["Normal Sphere (r=6370997)",6370997,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],
    PARAMETER["latitude_of_center",0],
    PARAMETER["longitude_of_center",65.01172638348839],
    PARAMETER["azimuth",-8.6347579930128],
    PARAMETER["rectified_grid_angle",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]

So the geotiff contains enough info to produce a proj4 string.

Any suggestion how solve this?

[alexander-buerger-met-no]

The header file geo_ctrans.inc from libgeotiff-dev has this projection with name CT_HotineObliqueMercatorAzimuthCenter. I need to find some documentation for the parameters.

[alexander-buerger-met-no]

Closed by e0d6b29 and 80edc69.