OGC Maps

[guillemc23]

Hi!

I'm testing loading a OGC Maps response for a STAC Item and a COG. They both work when querying them in PNG format, but, are they supposed to be georreferenced when querying for them as Tiff? I checked the OGC Maps spec and, to my understanding, they should.

The result I'm getting shows no georreferentiation, even though /tiles shows them in the right place.

[vincentsarago]

@guillemc23 can you share an example, I cannot reproduce on my side

# Without BBOX
gdalinfo /vsicurl/https://titiler.xyz/stac/map\?url\=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400\&assets\=visual\&f\=tif
Driver: GTiff/GeoTIFF
Warning 1: HTTP response code on https://titiler.xyz/stac/map?url=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400&assets=visual&f=tif.vat.dbf: 422
Files: /vsicurl/https://titiler.xyz/stac/map?url=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400&assets=visual&f=tif
Size is 1024, 943
Coordinate System is:
GEOGCRS["WGS 84",
    ENSEMBLE["World Geodetic System 1984 ensemble",
        MEMBER["World Geodetic System 1984 (Transit)"],
        MEMBER["World Geodetic System 1984 (G730)"],
        MEMBER["World Geodetic System 1984 (G873)"],
        MEMBER["World Geodetic System 1984 (G1150)"],
        MEMBER["World Geodetic System 1984 (G1674)"],
        MEMBER["World Geodetic System 1984 (G1762)"],
        MEMBER["World Geodetic System 1984 (G2139)"],
        MEMBER["World Geodetic System 1984 (G2296)"],
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]],
        ENSEMBLEACCURACY[2.0]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    USAGE[
        SCOPE["Horizontal component of 3D system."],
        AREA["World."],
        BBOX[-90,-180,90,180]],
    ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (47.549420289846744,-23.757848863896069)
Pixel Size = (0.000051792116221,-0.000051756981555)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  47.5494203, -23.7578489) ( 47d32'57.91"E, 23d45'28.26"S)
Lower Left  (  47.5494203, -23.8066557) ( 47d32'57.91"E, 23d48'23.96"S)
Upper Right (  47.6024554, -23.7578489) ( 47d36' 8.84"E, 23d45'28.26"S)
Lower Right (  47.6024554, -23.8066557) ( 47d36' 8.84"E, 23d48'23.96"S)
Center      (  47.5759379, -23.7822523) ( 47d34'33.38"E, 23d46'56.11"S)
Band 1 Block=1024x2 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA 
Band 2 Block=1024x2 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA 
Band 3 Block=1024x2 Type=Byte, ColorInterp=Blue
  Mask Flags: PER_DATASET ALPHA 
Band 4 Block=1024x2 Type=Byte, ColorInterp=Alpha

# With BBOX
 gdalinfo /vsicurl/https://titiler.xyz/stac/map\?url\=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400\&assets\=visual\&f\=tif\&bbox\=47.5,-23.8,47.6,-23.75\&crs\=EPSG:3857
Driver: GTiff/GeoTIFF
Warning 1: HTTP response code on https://titiler.xyz/stac/map?url=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400&assets=visual&f=tif&bbox=47.5,-23.8,47.6,-23.75&crs=EPSG:3857.vat.dbf: 500
Files: /vsicurl/https://titiler.xyz/stac/map?url=https://stac.eoapi.dev/collections/MAXAR_cyclone_emnati22/items/38_302000112320_105001002A31C400&assets=visual&f=tif&bbox=47.5,-23.8,47.6,-23.75&crs=EPSG:3857
Size is 1024, 560
Coordinate System is:
PROJCRS["WGS 84 / Pseudo-Mercator",
    BASEGEOGCRS["WGS 84",
        ENSEMBLE["World Geodetic System 1984 ensemble",
            MEMBER["World Geodetic System 1984 (Transit)"],
            MEMBER["World Geodetic System 1984 (G730)"],
            MEMBER["World Geodetic System 1984 (G873)"],
            MEMBER["World Geodetic System 1984 (G1150)"],
            MEMBER["World Geodetic System 1984 (G1674)"],
            MEMBER["World Geodetic System 1984 (G1762)"],
            MEMBER["World Geodetic System 1984 (G2139)"],
            MEMBER["World Geodetic System 1984 (G2296)"],
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]],
            ENSEMBLEACCURACY[2.0]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4326]],
    CONVERSION["Popular Visualisation Pseudo-Mercator",
        METHOD["Popular Visualisation Pseudo Mercator",
            ID["EPSG",1024]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["Web mapping and visualisation."],
        AREA["World between 85.06°S and 85.06°N."],
        BBOX[-85.06,-180,85.06,180]],
    ID["EPSG",3857]]
Data axis to CRS axis mapping: 1,2
Origin = (5287675.812680494971573,-2722973.976836008951068)
Pixel Size = (10.871044022780552,-10.860947851911520)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 5287675.813,-2722973.977) ( 47d30' 0.00"E, 23d45' 0.00"S)
Lower Left  ( 5287675.813,-2729056.108) ( 47d30' 0.00"E, 23d48' 0.00"S)
Upper Right ( 5298807.762,-2722973.977) ( 47d36' 0.00"E, 23d45' 0.00"S)
Lower Right ( 5298807.762,-2729056.108) ( 47d36' 0.00"E, 23d48' 0.00"S)
Center      ( 5293241.787,-2726015.042) ( 47d33' 0.00"E, 23d46'30.01"S)
Band 1 Block=1024x2 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA 
Band 2 Block=1024x2 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA 
Band 3 Block=1024x2 Type=Byte, ColorInterp=Blue
  Mask Flags: PER_DATASET ALPHA 
Band 4 Block=1024x2 Type=Byte, ColorInterp=Alpha

[vincentsarago]

oh but it doesn't with with tiff because of https://github.com/developmentseed/titiler/blob/main/src/titiler/core/titiler/core/utils.py#L111

I'll going to push a fix for this

[guillemc23]

I can confirm it works with tif with a single f 😀

The fix looks great, I'll try to test it soon