Nieuwe regioindeling Veen

[BartRijken]

Deze shapefiles bevatten de grenzen per landschap, voor de export van de indicatoren. We gaan de door Claude voorgestelde/gemaakte grenzen dus niet gebruiken voorlopig.

Waarom deze grenzen? Omdat de landschapsteams zelf ook berekeningen hebben gedaan. Die hebben ze geaggregeerd op basis van deze grenzen. Wij willen hierop aansluiten.

Merk op dat deze grenzen soms overlappen. Dat is voor dit doel niet erg, zo lang we maar niet over de regio's aggregeren. Des te belangrijker om alle indicatoren ook op NL schaal te exporteren, vanuit hun oorspronkelijke grid domain. Elke indicator wordt dus zowel op landschap niveau als voor NL geexporteerd (en waar mogelijk/aangegeven als kaart).

[BartRijken]

• Veen: https://deltares.data.surf.nl/index.php/f/31732645
• Rivieren: https://nl2120.sharepoint.com/:u:/r/sites/Werkstromen/Gedeelde%20documenten/Werkstroom%20overstijgende%20producten/Toekomstbeelden/Kernteam/Synthese%20RuimteScanner/Begrenzing%20riviergebieds.zip?d=wbf9b456a4abd43ec9f1cc6dc6d91435d&csf=1&web=1&e=NqgzZW
• Zand: https://nl2120.sharepoint.com/:u:/r/sites/Werkstromen/Gedeelde%20documenten/Werkstroom%20overstijgende%20producten/Toekomstbeelden/Kernteam/Synthese%20RuimteScanner/?d=we86ac03f48104ca38ee1d3d65b83660f&csf=1&web=1&e=TXdGuw)
• Kust:https://nl2120.sharepoint.com/:u:/r/sites/Werkstromen/Gedeelde%20documenten/Werkstroom%20overstijgende%20producten/Toekomstbeelden/Kernteam/Synthese%20RuimteScanner/begrenzing_zand.zip?d=we86ac03f48104ca38ee1d3d65b83660f&csf=1&web=1&e=TXdGuw

Ad Veen: Hiermee kun je dus de aggregaties doen voor de verschillende toekomstbeelden. Het bestand bestaat uit verschillende polygonen, omdat het de peilvakken beslaat. Ik had ook een dissolved versie gemaakt, maar dan komen er kleine gaatjes in, dus denk dat dat binnen de Ruimtescanner sneller op te lossen is.

Ad Rivieren: Voor rivieren kan je de shape gebruiken uit je mail; en bij overlap met veen mag veen bovenop.

[jipclaassens]

waarom nu een anderen dan #705?

[jipclaassens]

Ingebouwd en doorgerekend. De vier aangeleverde begrenzingen zijn nu de indeling waarop de indicatoren worden gerapporteerd; de indeling die ik eerder voor #705 maakte blijft in de configuratie staan maar wordt daar niet meer voor gebruikt.

Hoe de overlap is opgelost

De opmerking dat de grenzen elkaar overlappen is de kern van het ontwerp geworden. Gemeten op het raster: kust en zand 113.646 ha, kust en veen 68.805 ha, rivieren en veen 34.464 ha, veen en zand 33.331 ha. Kust en rivieren raken elkaar niet.

Daardoor kan er geen gedeelde kolom zijn die per cel zegt in welke regio hij ligt, want een cel kan bij twee landschappen horen. Elk landschap staat daarom als een eigen eenheid van een regel met alleen een masker, en de indicatorentabel wordt per landschap uitgedraaid met dezelfde machinerie als voor Nederland. Optellen over de landschappen is daarmee per constructie onmogelijk, wat aansluit bij de waarschuwing in dit issue.

Kiezen gaat met ModelParameters/Advanced/IndicatorRegio_ref, die vier opties erbij heeft gekregen. Die parameter is nu ook met de omgevingsvariabele IndicatorRegio te zetten, net als ExportZichtjaar, zodat de vier runs te scripten zijn. Elke run schrijft twee tabellen: die van het gekozen landschap en de landelijke. Vier runs geven dus de vier landschappen plus Nederland, en dat kost ongeveer een kwartier per landschap.

De bestanden

Drie ervan staan op de share met datum in de naam. Zand niet: die levering is geometrisch identiek aan Zand_Claimregios_20260708.shp die er al stond, tot op de hectare en zelfs dezelfde bestandsgrootte, dus de configuratie wijst daarnaar in plaats van een tweede kopie te maken die uit elkaar kan gaan lopen.

Voor veen zaten er twee bestanden in de levering die een factor twee schelen. Ik houd de peilvakversie aan van 354.173 ha, want daarop aggregeert het veenteam zelf. Het andere bestand is de deelgebiedenshape van 732.079 ha, die geometrisch gelijk is aan wat de waterbergingsregio in het model al leest.

De 2.459 gaten in die peilvakversie beslaan samen 30.238 ha, 8,5 procent. Ze zijn niet allemaal gedicht: er zitten echte uitsparingen tussen tot 3.072 ha. Een district-analyse op het complement dicht alleen aaneengesloten stukken onder een drempel, standaard 10 ha, wat 2.826 ha oplevert. Dat haalt de slivers tussen aangrenzende peilvakken weg en laat steden en plassen open. De drempel is een parameter.

Controle

De rastering is naast een meting op de shapes zelf gelegd. Rivieren 351.657 tegen 351.655 ha, zand 1.668.976 tegen 1.668.981, veen ongedicht 354.144 tegen 354.173. Kust komt op 1.811.344 tegen 1.836.437, en dat verschil van 25.094 ha is exact het deel dat ten oosten van x 280.000 buiten het modelraster valt, dus op zee en over de grens.

Twee landschapstabellen zijn doorgerekend voor WLO_hoog_BAU in Y2040, veen en kust, de twee lastigste. Veen beslaat 10,2 procent van Nederland en de algemene indicatoren liggen daar dicht bij: verstedelijking 15,2 procent, nieuwe natuur 11,9, overstromingsschade 10,4. De veenspecifieke indicatoren steken er ver bovenuit, met 76,1 procent van de SOMERS-emissies en 91,7 procent van de veenreserve. Bij kust valt op dat de buitendijkse overstromingsschade er met 33,9 procent hoog uitkomt, wat past bij een begrenzing waarvan de helft op zee ligt.

De overlap is daarmee ook meteen zichtbaar: de SOMERS-emissies zijn voor 76,1 procent aan veen toegerekend en voor 28,2 procent aan kust, samen 104,3 procent. Dat is geen fout maar precies wat overlappende grenzen doen, en het is de reden dat deze tabellen niet bij elkaar opgeteld mogen worden.

Twee dingen om te weten

Van Nederland valt 469.871 ha buiten alle vier de landschappen. Dat is geen restklasse maar gebied waarover geen enkel team iets zegt. Die hectares staan in de landelijke tabel en in geen enkele landschapstabel. Vier landschapstabellen naast elkaar geven dus in twee richtingen niet heel Nederland: er ontbreekt gebied en er is dubbeltelling.

Kust is een begrenzing van het kustsysteem en niet van de kustzone op land. Van de 1.811.344 ha binnen het modelraster ligt 901.219 ha binnen de provinciegrenzen; de rest is zee. Voor de landgebonden indicatoren maakt dat niets uit, want daar staat op zee niets, maar voor een oppervlaktevraag wel.

De rivieren- en zandtabel zijn nog niet uitgedraaid. Dat is dezelfde machinerie met een ander masker, en beide maskers zijn hierboven al tegen de vector getoetst, dus dat is een run en geen bouwwerk. Ze komen mee in de eerstvolgende opleveringsronde.

De documentatie staat op de wiki onder Landschapsgebieden.

[BartRijken]

@jipclaassens hierbij de juiste begrenzing van veen:: https://deltares.data.surf.nl/index.php/f/31758650

Net ontvangen van team Veen.

Graag de vorige vervangen.