[Dutch] 3-uur neerslagaccumulaties van het KNMI

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

[Dutch] 3-uur neerslagaccumulaties van het KNMI

Paul Meems (Top-X)
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




_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch
Reply | Threaded
Open this post in threaded view
|

Re: [Dutch] 3-uur neerslagaccumulaties van het KNMI

SiebeBosch

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

_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch
--

_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch
Reply | Threaded
Open this post in threaded view
|

Re: [Dutch] 3-uur neerslagaccumulaties van het KNMI

Paul Meems (Top-X)
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:

Bedankt voor het meedenken.


Met vriendelijke groet,

Paul




Op 13 december 2016 15:50 schreef Siebe Bosch <[hidden email]>:

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

_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch
--

_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch


_______________________________________________
Dutch mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/dutch