[Dutch] 3-uur neerslagaccumulaties van het KNMI

[Siebe Bosch]

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 
> <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
>
> **
>
> TopX Geo-ICT <http://www.topx-geo-ict.nl>
>
> 	
>
> Paul MeemsSenior GIS consultant06-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 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/b141b5f0/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/b141b5f0/attachment-0001.gif>