[Dutch] 3-uur neerslagaccumulaties van het KNMI

[Paul Meems]

Bedankt Siebe en ook Marc Coevoet voor jullie reacties.

Ongeduldig als ik ben had ik de vraag ook op de GDAL-Dev mailinglijst
gepost en daar heb ik de juiste gegevens gekregen, inclusief een mooie
uitleg:
http://osgeo-org.1560.x6.nabble.com/gdal-dev-Convert-precipitation-accumulations-from-HDF5-to-GeoTiff-td5299661.html#a5299706

Bedankt voor het meedenken.


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>



Op 13 december 2016 15:50 schreef Siebe Bosch <siebe op hydroconsult.nl>:

> Paul,
>
> Ik gebruik voor conversie van KNMI naar RD de volgende argumenten in
> GDALWARP (niet gdal_translate):
>
> D:\gdal64\bin\gdalwarp -of netCDF -dstnodata -999 -te 9000 29900 296000
> 651000 -tr 1000 1000 -s_srs "+proj=stere +x_0=0 +y_0=0 +lat_0=90 +lon_0=0
> +lat_ts=60 +a=6378.388 +b=6356.906" -t_srs "EPSG:28992"
> Ik zie wat kleine verschillen met jouw projectiestring.
>
> mvg,
> Siebe Bosch
>
>
>
> On 12-12-16 12:17, Paul Meems wrote:
>
> 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>
>
> _______________________________________________
> Dutch mailing listDutch op lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/dutch
>
> --
>
> _______________________________________________
> Dutch mailing list
> Dutch op lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/dutch
>
------------- volgend deel ------------
Een HTML-bijlage is gescrubt...
URL: <http://lists.osgeo.org/pipermail/dutch/attachments/20161213/37cba19b/attachment-0001.html>
------------- volgend deel ------------
Een niet-tekst bijlage is gescrubt...
Naam: signature_60pct.gif
Type: image/gif
Grootte: 30165 bytes
Omschrijving: niet beschikbaar
URL : <http://lists.osgeo.org/pipermail/dutch/attachments/20161213/37cba19b/attachment-0001.gif>