search

swisstopo / topo-satromo

Erdbeobachtungs-Satellitendaten fürs Trockenheitsmonitoring (SATROMO)
BSD 3-Clause "New" or "Revised" License
13 stars 3 forks source link

_metadata.json : add GDAL command to generate RGB display #74

Open davidoesch opened 5 months ago

davidoesch commented 5 months ago

Is your feature request related to a problem? Please describe. There is a need to add a field with a GDAL command to each _metadata.json file as for the gdal enhancement function, so BGDI can use this info to generate VRT see https://jira.swisstopo.ch/browse/BGDIDIC-2772

Describe the solution you'd like A field should be added to each _metadata.json file for the _10m band containing the GDAL command: gdal_translate /vsicurl/https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif gdal_translate_byte.vrt -a_nodata 9999 -scale -exponent 0.5 -of VRT. having -a_nodata 9999 and https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif changed according to the file

Describe alternatives you've considered Currently FSDI is generating the VRT files for each file and stores it internally

Additional context Once this implemented, KOGIS will use the provided command and as fallback their json

davidoesch commented 5 months ago

maybe even try generate the VRT conent and store it in the json: more complex band math could be applied and the real TRUE COLOR result achieved:

Start with `

PROJCS["CH1903+ / LV95",GEOGCS["CH1903+",DATUM["CH1903+",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6150"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4150"]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",46.9524055555556],PARAMETER["longitude_of_center",7.43958333333333],PARAMETER["azimuth",90],PARAMETER["rectified_grid_angle",90],PARAMETER["scale_factor",1],PARAMETER["false_easting",2600000],PARAMETER["false_northing",1200000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","2056"]] 2.4744300000000000e+06, 1.0000000000000000e+01, 0.0000000000000000e+00, 1.3032400000000000e+06, 0.0000000000000000e+00, -1.0000000000000000e+01 Area DEFLATE PIXEL 9999 Gray /vsicurl/https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif 1 9999 0.5 0 12964 0 255 9999 /vsicurl/https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif 2 9999 0.5 0 13200 0 255 9999 /vsicurl/https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif 3 9999 0.5 0 13556 0 255 9999 /vsicurl/https://data.geo.admin.ch/ch.swisstopo.swisseo_s2-sr_v100/2024-03-19t102639/ch.swisstopo.swisseo_s2-sr_v100_mosaic_2024-03-19t102639_bands-10m.tif 4 9999 0.5 0 12625 0 255

`

davidoesch commented 2 months ago

Feedback from planet/sinergsie "EO Browser is powered by a proprietary technology Sentinel Hub. " -> no info how they operate their system