Onderwerp: Kustfijn-model

FMJ schreef : Zomaar een gokje, kan het raster gebaseerd zijn op de vlieglijnen met laseraltimetrie metingen van RWS.

Die screenshots van jou op z'n tijd laten dat niet echt zien anders. Dan zie ik korrelige plaatjes met heel grote rasterboeien en kleinere vectorboeien uit het RWS-bestand daar dwars overheen. Niet echt elegant. Al helemaal niet als je je bedenkt dat die rasterkaarten door de Hydrografische Dienst gewoon gerenderd worden uit een database, eigenlijk screenshots zijn van een vectorkaart.

Los daarvan wil ik kunnen kiezen of ik voor mijn routering vertrouw op de rekenmogelijkheden van WinGPS of toch liever op de internationaal bewezen rekenmogelijkheden van Expedition, Adrena of QT-VLM... In dat licht bezien zijn de OES-ENC kaarten een grote vooruitgang. De digitale gegevens van de Hydrografische Dienst zijn daardoor eindelijk niet meer alleen beschikbaar voor klanten van Stentec.

FMJ schreef : Uit mijn plaatje is het niet goed te zien, maar die lijnen variëren niet alleen in afstand, maar ze lopen ook mee met de kromming van de kust. Bij Den Helder lopen ze weer bijna O-W en N-Z.
Er zijn wel libraries die interpolatie van punt-data naar een raster kunnen doen, dus daar ga ik nog een keer induiken.


Voor dcsmv6_zunov4_zuno_kf_hirlam moet je Maps2D hebben

Baasklusje schreef : De belangrijkste reden voor IMO om over te gaan naan vector kaarten was tweeledig.
Minder geheugen.(Wel degelijk van belang als je een groter gebied dan alleen NL aan kaarten hebt)
Eenvoudig updaten met minimaal data gebuik. (Ook op zee via een sateliet je wereldwijde kaarten set iedere week updaten)

En het voordeel dat rasterkaarten gemaakt worden door ervaren kaartenmakers, die kunnen bepalen wat waar moet komen en hoe groot staat als een huis. Daar kan (nog) geen computerprogramma tegenop. (Er wordt momenteel de laatste hand gelegd aan een nieuwe versie voor de ENC codering(S57), waarin de kaartenmakers ook meer invloed hebben op het uiteindelijke kaartbeeld.)

Menno schreef : Dat zou gaaf zijn!

Wanneer ik op http://noos.matroos.rws.nl/maps2d/start/ kies voor dcsm-v6 kf hirlam, kan ik rechts kiezen op welk raster ik de stroomvectoren wil zien. RWS zelf heeft die interpolatie dus al in de vingers..

edit: als ik hier kijk: http://noos.matroos.rws.nl/direct/get_map2series.php?source=&unit=&x=&y=&tstart=&tstop=&tinc=&coordsys= dan zou ik dmv de -loc flag gevolgd door heel veel coordinaten de boel op m'n eigen gewenste rooster op kunnen vragen..

edit2: als ik bij maps2d links onderaan kies voor Data ipv Image, dan kan ik zelf m'n grid en grid-spacing opgeven... Daar kan ik dan voor gewenst gebied en gewenst tijdvak een file opvragen met data op het door mij gewenste rooster!

Was er inmiddels ook achter dat die flag verwijst naar een systeem, niet naar een kaart.

Als ik het goed begrijp zit de eerder genoemde 10-minuten data overigens in hmcn_kustfijn.. bovengenoemde gaat per uur. Die hmcn_kustfijn lijkt echter niet meer beschikbaar.

RWS kan volgens mij de interpolatie al voor je doen:

Of zie ik het verkeerd?

Dat denk ik ook, zie mijn 2e edit

NetCDF lijkt vrij eenvoudig converteerbaar naar GRIB(2) mbv CDO.

Hydrodynamische modellen werken dikwijls met een onregelmatig rekenrooster. Op die manier kan een hoge resolutie geboden worden waar nodig ( bvb zeegaten) zonder over het hele modeldomein een hoge resolutie te moeten gebruiken. Zo kan de rekentijd ingekort worden...

FMJ schreef :
Dat gaat waarschijnlijk niet zonder een interpolatieslag als de oorspronkelijke data op een onregelmatig rooster zit.

RWS biedt zelf al een geïnterpoleerde versie aan, zie de post van Christiaan hierboven

Excuus - overheen gelezen. In dat geval lukt het waarschijnlijk wel om van netcdf naar grib te converteren.

Overigens is NetCDF een veel prettiger formaat om mee te werken dan grib(2), maar helaas zijn niet alle weerdata viewers hier op ingericht. Mocht je een viewer zoeken voor netCDF bestanden dan is bijvoorbeeld IDV een optie.

Het grid in die geïnterpoleerde versie lijkt een interpolatie van het curvilinear grid te zijn, en dus niet omgezet naar rectangular lonlat.
Het is me ondertussen wel gelukt met CDO

Voorbeeld:


XyGrib:


OpenCPN:

Wat goed Menno!!

Fantastisch!

Móóóí!!

Super gaaf Menno!

Mooi ook om te zien hoe we met elkaar hier in korte tijd mooie dingen voor elkaar krijgen!


We kunnen de aanpak natuurlijk ook om gaan draaien: XyGrib en OpenCPN NetCDF leren lezen.. Misschien kan sommige moderne routeringssoftware dat al wel?!

Ik geloof dat er hier weer mooie dingen ontstaan

De interpolatie gaat nog niet helemaal goed. Bij de eerdere plaatjes was dat al zichtbaar ten noorden van Vlieland, maar als ik een groter gebied probeer gaat er nog meer mis:

Ja wat jammer zeg! Maar CDO kan best veel, dus wie weet!

Neem aan dat je al met de -setmisstonn en de diverse -setremap* flags gespeeld hebt?

Ik weet niet precies wat CDO precies doet en kan dezer dagen (al een tijd niet meer gebruikt), maar ik vrees dat je voor betere interpolatie zelf wat zult moeten schrijven. Met name nabij de kust is dat vrij lastig goed te krijgen. (niet over eiland-punten heen middelen bijvoorbeeld, en niet extrapoleren over land). Om het goed te doen heb je een landmasker op het nieuwe rooster nodig.

Als je een rooster kiest met een resolutie dicht bij het oorspronkelijke grid, dan kun je met een simpele "nearest" wel een heel eind komen. Bijvoorbeeld Matlab of octave, en dan interpn() gebruiken (regelmatig -> regelmatig) of griddata (onregelmatig -> regelmatig).

Zit hier nog beweging in?

In hoeverre zijn de strömungsdateien van de grib1 en 2 files van BSH hiermee te vergelijken?
Minder fijn misschien maar toch ook met op en afwaai dacht ik?