[Dutch] 3-uur neerslagaccumulaties van het KNMI

[Paul Meems]

Ik probeer om de 3-uur neerslagaccumulaties van het KNMI (
https://data.knmi.nl/datasets/radar_tar_corr_accum_03h/2.0) te gebruiken.

Al onze geo-data is in WGS84, dus mijn eerste stap is het converteren van
het HDF5-bestand van KNMI naar een andere projectie en dan gelijk ook naar
GeoTiff wat toch wat makkelijker te verwerken is.
Hier gebruiken we de GDAL tools voor, v2.0.
Als ik rechtstreeks GDAL_Warp gebruik krijg ik een foutmelding over
ontbrekende transformatie matrix. Via Google heb ik soortgelijke problemen
gevonden en daar was de oplossing om eerst GDAL_Translate te doen en
dan de *georeferenced
bounds* toe te voegen.

Dat doe ik moet dit commando:

gdal_translate -of GTiff -a_nodata 65535 -a_srs "+proj=stere +lat_0=90
+lon_0=0 +lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0" -a_ullr 0
55.973602 9.0092793 48.895298
"HDF5:\"RAD_NL25_RAC_03H_201612041200.h5\"://image1/image_data"
"test-translate.tif"


-a_ullr moet dan de bounds zijn. Die haal ik uit de meta-data van het
bestand.

De projectiestring komt ook uit de meta-data.


Vervolgens kan ik dan een reprojectie doen met:

gdalwarp -overwrite -s_srs "+proj=stere +lat_0=90 +lon_0=0 +lat_ts=60
+a=6378.14 +b=6356.75 +x_0=0 y_0=0" -t_srs EPSG:3857 -dstnodata 65535
-of GTiff "test-translate.tif" "test-warp.tif"


Het resultaat had dan in NL moeten liggen met coördinaten van 52 07' 25'',
5 15' 46'' (Amersfoort)

Maar ik krijg 89 29' 23'', 173 44' 12''


Volgens mij zijn er drie mogelijke oorzaken:


   1. Ik doe de conversie en reprojectie niet goed
   2. Ik geef niet de juiste bounds op in -a_ullr
   3. Ik gebruik niet de juiste projectiestring

Ik ben vast niet de eerste die dit wil doen, weet iemand de juiste workflow?


Alvast bedankt.


Met vriendelijke groet,

Paul



[image: TopX Geo-ICT] <http://www.topx-geo-ict.nl>

Paul Meems
Senior GIS consultant
06-53989481

Wij bieden ondersteuning voor MapWindow GIS <http://www.mapwindow.org/>

Download the latest MapWinGIS mapping engine.
<https://github.com/MapWindow/MapWinGIS/releases>

Download the latest MapWindow 5 open source desktop application.
<https://github.com/MapWindow/MapWindow5/releases>
------------- volgend deel ------------
Een HTML-bijlage is gescrubt...
URL: <http://lists.osgeo.org/pipermail/dutch/attachments/20161212/38fa77dd/attachment.html>