search

threedi / hhnk-threedi-plugin

QGIS plugin gebruikt voor analyses van 3Di modellen bij HHNK.
1 stars 0 forks source link

Netcdf essentials #53

Open wvangerwen opened 1 year ago

wvangerwen commented 1 year ago

Resultaten klaar zetten voor meerdere toepassingen. Legger, rasters.

Python class

wvanesse commented 1 year ago

Voorstel naamgeving en indeling: Twee methoden, één voor 1D berekening en één voor 2D bereningen.

Voor 1D berekening (0d1d_test / HT)

Laagnaam: resultaat_polder_0d1d_rev0_1D_punt.shp (ja shape omdat het ook in esri kan) Bevat de kolommen:

  1. id: id van het rekenpunt
  2. type: nuttig? in result zit 1d/2d, maar het is allemaal 1d...
  3. H_ini: initiele waterstand op punt
  4. H_P0: waterstand voor de bui begint
  5. H_Pend: waterstand aan het einde van de bui
  6. H_Pend_1h: waterstand 1 uur voor het einde van de bui
  7. H_end: waterstand aan einde simulatie
  8. dH_ini_P0: (H_P0 - H_ini) verschil in waterstand tussen initieel en start neerslag (check op uitzakken/lekken peil)
  9. dH_P0_end: (H_Pend - H_P0) verschil in waterstand tussen start neerslag en einde neerslag (check op maximale opstuwing)
  10. dH_Pend: (H_Pend - H_Pend_1h) verschil waterstand tussen einde regen en einde regen min 1 uur (check op stationaire situatie tijdens berekening)
  11. dH_end: (H_end - H_P0) verschil waterstand einde simulatie en voor start bui (check herstel peil)
  12. geometry; geometrie van het punt

Laagnaam: resultaat_polder_0d1d_rev0_1D_lijn.shp Bevat de kolommen:

  1. id: id van het rekenpunt
  2. type: type is hier watergang of een type kunstwerk (, maar ook 1d2d of 2d2d.)
  3. Q_norm: afvoer gemiddeld tussen einde bui en 1 uur voor einde bui (zo is het nu toch?) m3/s
  4. U_norm: stroomsnelheid op dezelfde moment/manier als Q_norm m3/s
  5. nodID_up: bovenstrooms rekenpunt id
  6. nodID_dow: benedenstrooms rekenpunt id
  7. H_ini_up: bovenstrooms initiele waterstand
  8. H_ini_dow: benedenstrooms initiele waterstand
  9. H_Pend_up: maximale waterstand bovenstrooms
  10. H_Pend_do: maximale waterstand benedenstrooms
  11. dH_m: absoluut verschil in waterstanden
  12. length: lengte segment
  13. S_cm_km: verhang in cm per kilometer
  14. flow_dir: richting van stroming tov tekenrichting

Ik twijfel wat te doen met de bui eigenschappen. Ze zijn belangrijk voor checks, maar het is gewoon heel vaak dezelfde waarde. Dus in aparte tabel of als kolommen met telkens zelfde getal? Laagnaam: bui_polder_0d1d_rev0_1D.dbf

  1. modelname
  2. modelslug
  3. sim_dur: duur van de berekening
  4. T_P0: moment waarop neerslag begon
  5. T_Pend: moment waarop neerslag stopt
  6. Pmax: maximale neerslagintensiteit omgerekend naar mm/uur (als het goed is altijd 100)
  7. Ptot: totaal neerslagvolume (in mm?) Nb. kunnen we hier ook de 'pumplines' aan toevoegen? Die ontbreken, hoort blijkbaar niet bij '1D_all'?

Alternatief is alleen de lijnen als shape te exporteren zodat het naar de leggertool en geoweb kan.

Voor 2D berekeningen

Deze hoeft in principe niet naar esri dus zouden we gepackage van kunnen maken. Is we mooi als de 2D vlakken en 1D punten in één bestand kunnen stoppen. Wat denk jij? Geopackage naam: resultaat_polder_rev01D2D.gpg type berekening kan zijn 1D2D test of TXXXX blok/piek gxg Laagnaam: polder_rev0_1D_punt Bevat de kollommen:

  1. id
  2. type
  3. H0: initiele waterstand
  4. Hmax: maximale waterstand
  5. V0: volume aan begin berekening
  6. Vmax: maximale volume
  7. Su_max: maximale natte oppervlak? Wil jij dit aanvullen met waarden die je nodig hebt voor de correctie en interpolatie van de waterdiepterasters?

Laagnaam: polder_rev0_lijn.shp Bevat de kolommen:

  1. id: id van het rekenpunt
  2. type: type is hier watergang, een type kunstwerk, 1d2d of 2d2d.
  3. Q_max: maximale afvoer absoluut m3/s
  4. U_max: maximale stroomsnelheid absoluut m3/s
  5. nodID_up: bovenstrooms rekenpunt id
  6. nodID_dow: benedenstrooms rekenpunt id

Laagnaam: bui_polder_rev0.dbf Zelfde als hierboven

Deze oplossing wilde ik jullie niet onthouden -_-

@wvangerwen ben benieuwd naar jouw commentaar

wvangerwen commented 1 year ago

Ik kom op zo iets uit. Zou in eerste instantie alles in gpkg zetten met leesbare kolomnamen. Als het dan echt nodig blijkt om het over te zetten naar shapefile dan zou ik nog eigenlijk liever csv zien. Dingen in 1 bestand hebben is een stuk overzichtelijker.

resultaat_polder_rev01D2D.gpkg

wvanesse commented 1 year ago

Ja klopt. Een bestand heeft voordelen. Het wordt wel per som een apart bestand neem ik aan hè?

Groet, Wouter

Op 2 mei 2023 om 08:54 heeft wvangerwen @.***> het volgende geschreven:



Ik kom op zo iets uit. Zou in eerste instantie alles in gpkg zetten met leesbare kolomnamen. Als het dan echt nodig blijkt om het over te zetten naar shapefile dan zou ik nog eigenlijk liever csv zien. Dingen in 1 bestand hebben is een stuk overzichtelijker.

resultaat_polder_rev01D2D.gpkg

— Reply to this email directly, view it on GitHubhttps://github.com/threedi/hhnk-threedi-plugin/issues/53#issuecomment-1530974678, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ALVEVP7G5LHKD2E7PXJJAIDXECVTFANCNFSM6AAAAAAXCJNLGU. You are receiving this because you were assigned.Message ID: @.***>

[HHNK]

Deze e-mail geldt alleen als formeel besluit als dat specifiek benoemd is in de mail of in de bijlage daarbij. Heeft u een formeel besluit nodig of twijfelt u over de rechtsgeldigheid van deze mail, neem dan telefonisch contact met ons op of kijk op onze website

wvanesse commented 7 hours ago

Plan NetCDF essentails

Extent existing python class (hhnk_threedi_tools.core.result_rasters.netcdf_to_gridgpkg import NetcdfToGPKG) to include:

Containing these layers with listed colums/data below.

points_1d

nodes1d_all

  1. id: id van het rekenpunt
  2. type: 1d connected/isolated/embedded
  3. H_ini: initiele waterstand op punt
  4. H_P0: waterstand voor de bui begint
  5. H_Pend: waterstand aan het einde van de bui
  6. H_Pend_1h: waterstand 1 uur voor het einde van de bui
  7. H_end: waterstand aan einde simulatie
  8. dH_ini_P0: (H_P0 - H_ini) verschil in waterstand tussen initieel en start neerslag (check op uitzakken/lekken peil)
  9. dH_P0_end: (H_Pend - H_P0) verschil in waterstand tussen start neerslag en einde neerslag (check op maximale opstuwing)
  10. dH_Pend: (H_Pend - H_Pend_1h) verschil waterstand tussen einde regen en einde regen min 1 uur (check op stationaire situatie tijdens berekening)
  11. dH_end: (H_end - H_P0) verschil waterstand einde simulatie en voor start bui (check herstel peil)
  12. geometry; geometrie van het punt

grid_2d

nodes2d_all -> 2d grid polygon

  1. id
  2. type
  3. H0: initiele waterstand
  4. Hmax: maximale waterstand
  5. V0: volume aan begin berekening
  6. Vmax: maximale volume
  7. Su_max: maximale natte oppervlak? Wil jij dit aanvullen met waarden die je nodig hebt voor de correctie en interpolatie van de waterdiepterasters?
  8. HTXXX: Maximum waterlevel every hour
  9. geometry (polygon)

lines_1d

From these sources:

Containing these properties:

  1. id: id van het rekenpunt
  2. type: nuttig? in result zit 1d/2d, maar het is allemaal 1d...
  3. H_ini: initiele waterstand op punt
  4. H_P0: waterstand voor de bui begint
  5. H_Pend: waterstand aan het einde van de bui
  6. H_Pend_1h: waterstand 1 uur voor het einde van de bui
  7. H_end: waterstand aan einde simulatie
  8. dH_ini_P0: (H_P0 - H_ini) verschil in waterstand tussen initieel en start neerslag (check op uitzakken/lekken peil)
  9. dH_P0_end: (H_Pend - H_P0) verschil in waterstand tussen start neerslag en einde neerslag (check op maximale opstuwing)
  10. dH_Pend: (H_Pend - H_Pend_1h) verschil waterstand tussen einde regen en einde regen min 1 uur (check op stationaire situatie tijdens berekening)
  11. dH_end: (H_end - H_P0) verschil waterstand einde simulatie en voor start bui (check herstel peil)
  12. geometry; geometrie van het punt

lines_2d

Types 1d2d, 2d2d

Containing these properties:

  1. id: id van het rekenpunt
  2. type: type is hier watergang, een type kunstwerk, 1d2d of 2d2d.
  3. Q_max_pos: maximale afvoer m3/s
  4. Q_max_neg
  5. U_max_pos: maximale stroomsnelheid m3/s
  6. U_max_neg
  7. nodID_up: bovenstrooms rekenpunt id
  8. nodID_dow: benedenstrooms rekenpunt id
  9. geometry; geometrie

settings

Some summary of important settings.

rain

Rainfall timeseries if applicable

breaches

List of active breaches with these properties

  1. id of breach
  2. breach_start
  3. Max_breach_depth
  4. Time_vertical_growth
  5. Max breach width
  6. Q_max_pos: maximale afvoer m3/s
  7. Q_max_neg
  8. U_max_pos: maximale stroomsnelheid m3/s
  9. U_max_neg
  10. nodID_up: bovenstrooms rekenpunt id
  11. nodID_dow: benedenstrooms rekenpunt id
  12. geometry; geometrie
  13. QTXXX: discharge every hour
  14. UTXXX: flow velocity every hour
  15. BWTXXX: breach width every hour