Kuinka lataan Syken WCS-rajapinnalta GeoTIFF-muotoisen rasterin?
Kuinka lataan Syken WCS-rajapinnalta GeoTIFF-muotoisen rasterin?
SYKEn satelliitihavaintoja on saatavilla GeoTIFF-datarastereina WCS-latausrajapinnan kautta. Rajapinta on käytettävissä sekä ohjelmallisesti että esimerkiksi internet-selaimella. Luettelo Syken avoimista rajapinnoista.
WCS-rajapinnan käyttöohjeet
Capabilities-dokumentti
Ennen WCS-pyynnön muodostamista on varmistettava, missä koordinaatistossa haluttu rajapinta on julkaistu. Satelliittidatat ovat joko WGS84 (EPSG:4326) tai ETRSTM35-FIN (EPSG:3067) -koordinaatistossa. Tämä selviää rajapinnan Capabilities-dokumentista, jonka saa seuraavalla pyynnöllä:
https://geoserver2.ymparisto.fi/geoserver/eo/wcs?request=GetCapabilities
Palautteesta voi etsiä haluamansa datasetin (tason) nimen perusteella kuvauksen, jossa koordinaattijärjestelmä on kerrottu elementissä ows:BoundingBox.
Datapyyntö
WFS- ja WCS-palvelujen sovelluskäyttöä varten suosittelemme ottamaan käyttöön sovelluskohtaisen yksilöivän tunnisteen, jonka saa pyytämällä Syken gistuesta gistuki@syke.fi. Pyytämällä tunnisteen hyväksyt, että yhteystietosi tallennetaan Syken tietosuojakäytäntöjen mukaisesti. Tietosuojailmoitus Muista lisätä tunnistusavaimesi datapyynnön loppuun …&sykeuserid=“oma avain” kuten ensimmäisessä esimerkissä on havainnollistettu.
Rajapintakutsuun on määritettävä
- Rajapinnan palvelimen osoite: https://geoserver2.ymparisto.fi/geoserver/ows
- Palvelun tyyppi: service=WCS
- Palvelun versionumero: version=2.0.1
- Pyynnön tyyppi: request=GetCoverage
- Alue:
- SUBSET=Lat(y1,y2)&SUBSET=Long(x1,x2) , esim. Lat(61.01847,61.33005)&SUBSET=Long(20.77862,21.50516)
- Lat & Lon, kun kyseessä on WGS84-koordinaattijärjestelmästä (EPSG:4326) oleva aineisto.
- Lat: leveyspiirin koordinaatit, ensin eteläisempi, sitten pohjoisempi
- Long: pituuspiirin koordinaatit, ensin läntisempi ja sitten itäisempi
- Nämä voi poimia esim. Tarkkan karttanäkymästä
- SUBSET=N(y1,y2)&SUBSET=E(x1,x2), esim. SUBSET=N(6704362.2,6710330.2)&SUBSET=E(496609.7,514108.3),
- N & E kun kyseessä on ETRS89/TM35FIN -koordinaattijärjestelmässä (EPSG:3067) oleva aineisto.
- N: pohjoiskoordinaatit, E: itäkoordinaatit
- Näitä koordinaatteja voi poimia esim. kansalaisen karttapaikasta.
- SUBSET=Lat(y1,y2)&SUBSET=Long(x1,x2) , esim. Lat(61.01847,61.33005)&SUBSET=Long(20.77862,21.50516)
- Aineiston tasonimi: coverageId, esim. coverageId=eo:EO_SST, jossa eo viittaa aineiston sijaintiin palvelimella ja EO_SST on aineiston tasonimi
- Koordinaattijärjestelmä, jossa aineistoa halutaan ladata: subsettingcrs, esim. subsettingcrs=http://www.opengis.net/def/crs/EPSG/0/4326 (kun halutaan ladata wgs84:ssa)
- Aineiston lataustyyppi: format=image/tiff, kun halutaan ladata aineisto geotiffeinä
- Aika: SUBSET=time(), esim. SUBSET=time(“2021-04-20T00:00:00Z”), jos kyseessä on aineisto, jossa on esim. päivittäistä dataa
- sykeuserid sykeuserid=“oma avain”
Esimerkkejä
WGS84-koordinaattijärjestelmässä oleva päivän 20.4.2021 satelliittihavainnoista määritetty veden pintalämpötila (EO_SST) määritetyltä alueelta:
WGS84-koordinaattijärjestelmässä olevan päivän 26.4.2021 satelliittihavainnoista määritetty lumenpeittoala (EO_FSC) määritetyltä alueelta:
WGS84-koordinaattijärjestelmässä olevan päivän 16.7.2020 satelliittihavainnoista määritetty pintalevä (EO_HR_WQ_S2_ALGAE) märitetyltä alueelta:
ETRS89/TM35FIN-koordinaattijärjestelmässä olevan päivän 1.6.2021 satelliittihavainnoista määritetty sameusaineisto (EO_HR_WQ_S2_TURB) määritetyltä alueelta:
Osa Tarkkan aineistoista sijaitsee Ilmatieteenlaitoksen rajapinnoilla. Näille rajapintakutsuun määritetään
Rajapinnan palvelimen osoite: https://data.nsdc.fmi.fi/geoserver/ows
Esimerkiksi:
WGS84-koordinaattijärjestelmässä olevan päivän 31.8.2019-satelliittihavainnoista määritetty kasvillisuusindeksi (NDVI) määritetyltä alueelta:
ETRS89/TM35FIN–koordinaattijärjestelmässä olevan päivän 31.8.2019-satelliittihavainnoista määritetty kasvillisuusindeksi (NDVI) määritetyltä alueelta: