Deadwood-ai / deadwood-api

Main FastAPI application for the deadwood backend
GNU General Public License v3.0
0 stars 0 forks source link

some cogs are not visible. #46

Open JesJehle opened 2 months ago

JesJehle commented 2 months ago

ids: 2146, 2036, 2158, 2175

example id: 2183: https://data.deadtrees.earth/cogs/v1/7b1b130e-e6c1-47a6-94f0-a5e21d03467a_uavforsat_DDH005_ortho/7b1b130e-e6c1-47a6-94f0-a5e21d03467a_uavforsat_DDH005_ortho_cog_jpeg_ts_web-optimized_q75.tif

info of geotif:

Driver: GTiff/GeoTIFF
Files: 7b1b130e-e6c1-47a6-94f0-a5e21d03467a_uavforsat_DDH005_ortho.tif
Size is 35935, 25287
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)"],
        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 = (13.844565149956994,51.104451666057670)
Pixel Size = (0.000000305057217,-0.000000305057217)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=DEFLATE
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  13.8445651,  51.1044517) ( 13d50'40.43"E, 51d 6'16.03"N)
Lower Left  (  13.8445651,  51.0967377) ( 13d50'40.43"E, 51d 5'48.26"N)
Upper Right (  13.8555274,  51.1044517) ( 13d51'19.90"E, 51d 6'16.03"N)
Lower Right (  13.8555274,  51.0967377) ( 13d51'19.90"E, 51d 5'48.26"N)
Center      (  13.8500463,  51.1005947) ( 13d51' 0.17"E, 51d 6' 2.14"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
  NoData Value=255
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
  NoData Value=255
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
  NoData Value=255

The cog:

Driver: GTiff/GeoTIFF
Files: 7b1b130e-e6c1-47a6-94f0-a5e21d03467a_uavforsat_DDH005_ortho_cog_jpeg_ts_web-optimized_q75.tif
Size is 33024, 36864
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)"],
            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 = (1541161.582799855619669,6639798.008277200162411)
Pixel Size = (0.037322767717371,-0.037322767717371)
Metadata:
  AREA_OR_POINT=Area
  OVR_RESAMPLING_ALG=NEAREST
  TILING_SCHEME_NAME=WebMercatorQuad
  TILING_SCHEME_ZOOM_LEVEL=22
Image Structure Metadata:
  COMPRESSION=YCbCr JPEG
  INTERLEAVE=PIXEL
  JPEGTABLESMODE=1
  JPEG_QUALITY=75
  LAYOUT=COG
  SOURCE_COLOR_SPACE=YCbCr
Corner Coordinates:
Upper Left  ( 1541161.583, 6639798.008) ( 13d50'40.16"E, 51d 6'16.17"N)
Lower Left  ( 1541161.583, 6638422.142) ( 13d50'40.16"E, 51d 5'48.23"N)
Upper Right ( 1542394.130, 6639798.008) ( 13d51'20.02"E, 51d 6'16.17"N)
Lower Right ( 1542394.130, 6638422.142) ( 13d51'20.02"E, 51d 5'48.23"N)
Center      ( 1541777.856, 6639110.075) ( 13d51' 0.09"E, 51d 6' 2.20"N)
Band 1 Block=512x512 Type=Byte, ColorInterp=Red
  Overviews: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288
  Mask Flags: PER_DATASET
  Overviews of mask band: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288
Band 2 Block=512x512 Type=Byte, ColorInterp=Green
  Overviews: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288
  Mask Flags: PER_DATASET
  Overviews of mask band: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288
Band 3 Block=512x512 Type=Byte, ColorInterp=Blue
  Overviews: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288
  Mask Flags: PER_DATASET
  Overviews of mask band: 16512x18432, 8256x9216, 4128x4608, 2064x2304, 1032x1152, 516x576, 258x288

with no data:0: image

with no nodata: image

Seems that the overview is corrupted, zooming to a higher zoom level, its visible again: image

JesJehle commented 2 months ago

I canot figure this out for now. meaning some datasets canot be displayed until i found a solution. But