Multi Reclass Tool (MRT)
- Door: Hans Roelofsen (WEnR team Biodiversiteit & Beleid)
- Idee en eerste prototype: Henk Meeuwsen (WEnR team Biodiversiteit & Beleid)
- Begonnen: oktober 2020
- Laatste update: december 2022
Inhoud
Inleiding
Terminologie
Het combined raster
Beslisregels uitgelegd
Beslisregels toegepast
Bruto- en netto areaal
Downstream volgens newval
Downstream volgens firstrule
MRT installatie
Combined raster
Beslisregels opstellen
Command Line Interface
MRT output
Meer informatie
Inleiding
WENR ontwikkeld en onderhoud een portfolio aan ecologische modellen, zoals het MNP model. Veel modellen gebruiken GIS kaarten als input. Meestal zijn dit raster-kaarten waarin bepaalde (landgebruiks-)categorieen zijn aangeven. Sommige veelgebruikte rasters, zoals de Top10NL en de SNL Beheertypenkaart, zijn door WENR voorbewerkt tot raster-kaarten via het Basisbestand Natuur & Landschap.
Regelmatig blijkt dat een kaartlaag onvoldoende informatie bevat voor de gevraagde toepassing. Dit kan betrekking hebben op het aantal categorieen in de kaart en de omschrijvingen daarvan, of de vorm en omvang van de ruimtelijke contouren in de kaart. Een mogelijke oplossing hiervoor is twee of meer raster-kaarten te 'stapelen'. Binnen de stapel kaarten worden specifieke combinaties van categorieen geidentificeerd die de gewenste informatie verschaffen.
Een simpel voorbeeld is de combinatie van de Top10NL kaart met de SNL Beheertypenkaart. Eerstgenoemde onderscheidt naald-, loof- en gemengde bossen, laatstgenoemde onderscheidt natuur- en productiebos. In hun combinatie verschijnen zes nieuwe categorieen:
- natuur-loofbos
- gemengd natuurbos
- natuur-naaldbos
- productie-loofbos
- gemengd productiebos
- productie-naaldbos
De MRT is ontwikkeld om dit proces herhaaldelijk en gemakkelijk toe te passen. Het hoofddoel van de MRT is produceren van nieuwe raster-kaarten waar de nieuwe categorieen in staan. De MRT standaardprocedure is:
- identificeer de gewenste bronkaarten
- combineer de bronkaarten met de ArcGIS Pro Combine Tool
- formuleer de nieuwe categorieen en bijbehorende beslisregels
- activeer de MRT en genereer een nieuwe raster-kaart met de nieuwe categorieen.
In dit document wordt de theorie, ontwerp en toepassing van de MRT nader toegelicht.
Terminologie
- raster: geospatial raster bestand zoals gebruikelijk in GIS toepassingen. Tenzij anders vermeld wordt altijd een categorisch raster bedoelt (ie nominaal, discreet) en nooit een continue raster.
- categorie: een van de discrete categorieen in een categorisch raster. Refereert in MRT context wordt vrijwel altijd naar een landgebruikscategorie.
- value: een getal in een raster. Elke value in een raster refereert naar een bepaalde categorie.
- description: de omschrijving van een landgebruikscategorie, zoals zoals bos, water, stad etc.
- legenda: de koppeling tussen values en description voor een categorisch raster.
- categorisatie: alle categorieen binnen een raster
- combined-raster: output van de ArcGIS Pro Combine Tool.
- bron-raster: een van de rasters die is ingevoerd in de ArcGIS Pro Combine tool.
- combined-value: een integer getal in het combined-raster dat symbool staat voor een bepaalde combinatie van categorieen in alle source rasters.
- RAT of VAT: Raster (of Value) Attribute Table. Een tabel waarin de values van een raster worden geassocieerd met aanvullende informatie, zoals een description. In de RAT van een combi-raster staan de unique integer values assigned to each unique combination of input values (ESRI).
- beslisregel: een bepaalde combinatie van categorieen in een of meerdere source-rasters die resulteert in een nieuwe categorie.
- newval: de raster value die een reclass-rule toekent aan categorie behorende bij de reclass-rule
- upstream: het geheel aan reclass-rules die tezamen het combi-raster vertalen naar een output raster met newvals en bijbehorende descriptions.
- downstream: doorvertaling van newval kaart naar een andere categorisatie, waarbij newval categorieen samengeraapt kunnen worden tot een enkele downstream categorie.
-
model: algemene term voor een categorisatie die uit de MRT komt, bijvoorbeeld het
newval model
, of eendownstream-model
. - firstrule: de reclass-rule die als eerste van toepassing is op een pixel in het combined-raster en wiens newval wordt toegekend aan die pixel.
Theorie & ontwerp van de MRT
Het combined-raster
Een combined-raster is een raster waarin de values verwijzen naar specifieke value-combinaties van een n aantal bronrasters. Hieronder wordt het principe van een combined raster geillustreerd aan de hand van twee categorische rasters: Beheertypen en Top10. De beheertypen legenda is:
value | description |
---|---|
10 | N08.02 Open duin |
11 | N16.03 droog bos met productie |
12 | N16.04 vochtig bos met productie |
13 | N17.02 Drooghakhout |
De legenda van de Top10 kaart is:
value | description |
---|---|
100 | grasland |
142 | gemengd bos |
143 | loofbos |
144 | naaldbos |
160 | duin |
161 | duin met riet |
In de figuur hieronder zijn de beheertypen- en Top10 rasters weergegeven (links en midden). Rechts staat het combined_raster. De getallen in het combined_raster verwijzen naar dezelfde combinaties van waarden in de twee bronrasters.
De RAT van het combined_raster is dan:
value | beheertypen | top10 |
---|---|---|
1 | 10 | 160 |
2 | 10 | 100 |
3 | 10 | 143 |
4 | 11 | 142 |
5 | 11 | 143 |
6 | 11 | 144 |
7 | 11 | 100 |
8 | 12 | 161 |
9 | 13 | 144 |
De getallen in het combined_raster verklaren de combinaties van getallen in de twee bronrasters. Dit voorbeeld gebruikt twee bronrasters, maar in de praktijk worden meerdere bronrasters gecombineerd. De ArcGIS Pro Combine tool is gelimiteerd tot 20 bronrasters.
Beslisregels uitgelegd
Beslisregels identificeren gewenste categorieen uit een of meerdere bronrasters en koppelen daaraan een nieuwe categorie. De nieuwe categorie wordt gekenmerkt door een numerieke waarde (value) en een omschrijving (description).
Hieronder worden drie beslisregels als voorbeeld gegeven, gebaseerd op de twee eerder genoemde bronrasters. Het doel is om beheertype N16.03 droog bos met productie
te verfijnen naar loof-, naald- en gemengd bos, gebaseerd op Top10 informatie.
newval | description | remark | volgnr | Beheertypen | Top10 |
---|---|---|---|---|---|
10 | Gemengd droog bos met productie | 1 | 11 | 142 | |
11 | Droog loofbos met productie | 2 | 11 | 143 | |
12 | Droog naalbos met productie | 3 | 11 | 144 |
De beslisregels zijn opgebouwd uit de volgende componenten:
- newval: de waarde (value) van de nieuwe categorie in het resulterende raster
- description: de omschrijving van de nieuwe cateogire
- remark (optioneel): ruimte voor de gebruiker om de beslisregel nader toe te lichten
- volgnr: nummer ter identificatie van de beslisregel
- Top10: de gewenste categorieen uit het Top10 bronraster
- beheertypen: de gewenste categorieen uit het Beheertypen bronraster.
Het is mogelijk dat N16.03 droog bos met productie
niet altijd samenvalt met loof-, naald- of gemengd bos uit de Top10. Beslisregel #4 identificeert die plekken, door te stellen dat de Top10 value niet belangrijk is. Dit wordt aangeduid met een n
.
newval | description | remark | volgnr | Beheertypen | Top10 |
---|---|---|---|---|---|
13 | Overig droog bos met productie | 4 | 11 | n |
Beslisregels worden van boven naar beneden toegepast: zodra een locatie in het raster is toegekend aan een beslisregel wordt die niet meer overschreven door een lagere beslisregel. Dit geldt voor beslisregel #4, die van toepassing is op alle gebieden met N16.03. Beslisregels #1-3 zijn echter specifieker, of met ander woord strenger, omdat ze een aanvullende eis in de Top10 kaart bevatten. Omdat deze eerder worden toegepast, resteert voor beslisregel #4 alleen N16.03 buiten Top10 loof-, naald- en gemengd bos. Hierom is het belangrijk om strengere beslisregels te prioriteren boven lossere regels.
Meerdere beslisregels kunnen leiden tot dezelfde newval, bijvoorbeeld als informatie over dezelfde categorie afkomstig is uit verschillende bronrasters. Beslisregel #5 identificeert duinen uit Top10 en beslisregel #6 identifeert duinen uit de Beheertypenkaart. Beiden resulteren in dezelfde categorie: 'Duinen' (aangeduid met value 14).
newval | description | remark | volgnr | Beheertypen | Top10 |
---|---|---|---|---|---|
14 | Duinen | duinen Top10 | 5 | n | 160,161 |
14 | Duinen | duinen uit beheertypen | 6 | 10 | n |
Als beslisregels #5 en #6 gecombineerd zouden worden in een enkele beslisregel, geeft dat een wezenlijk ander resultaat:
newval | description | remark | volgnr | Beheertypen | Top10 |
---|---|---|---|---|---|
14 | Duinen | 7 | 10 | 160,161 |
Als er meerdere categorieen binnen een bronraster gespecifeerd worden, beschouwt de MRT deze als optioneel ten opzichte van elkaar. 160, 161
voor Top10 in beslisregel #5 wordt geinterpreteerd als 160 OF 161. Categorieen uit verschillende bronrasters zijn 'dwingend' ten opzichte van elkaar. Beslisregel #1 wordt dus geinterpreteer als 142 uit Top10 EN 142 uit beheertypen.
Dezelfde newval
moet altijd gekoppeld zijn aan dezelfde description
. De MRT controleert hierop en geeft feedback waar nodig. Tevens controleert de MRT of de opgegeven values in de de beslisregels, daadwerkelijk voorkomen in het betreffende bronraster. Als dat niet het geval is, wordt de MRT gestopt.
De beslisregels worden aan de MRT aangeboden via een Excel-sheet. Het betreffende tabblad wordt vaak aangeduid als het upstream
gedeelte van de MRT.
Beslisregels toegepast
De MRT vertaald elke beslisregel naar een database-query die wordt toegepast op de RAT van het combined_raster. Voor alle beslisregels ("regel") wordt hieronder aangegeven of het combined_raster voldoet aan de query (1/0).
value | beheertypen | top10 | regel01 | regel02 | regel03 | regel04 | regel05 | regel06 | regel07 |
---|---|---|---|---|---|---|---|---|---|
1 | 10 | 160 | 0 | 0 | 0 | 0 | 1 | 1 | 1 |
2 | 10 | 100 | 0 | 0 | 0 | 0 | 0 | 1 | 0 |
3 | 10 | 143 | 0 | 0 | 0 | 0 | 0 | 1 | 0 |
4 | 11 | 142 | 1 | 0 | 0 | 1 | 0 | 0 | 0 |
5 | 11 | 143 | 0 | 1 | 0 | 1 | 0 | 0 | 0 |
6 | 11 | 144 | 0 | 0 | 1 | 1 | 0 | 0 | 0 |
7 | 11 | 100 | 0 | 0 | 0 | 1 | 0 | 0 | 0 |
8 | 12 | 161 | 0 | 0 | 0 | 0 | 1 | 0 | 0 |
9 | 13 | 144 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Vervolgens voegt de MRT 'beslisregel #999' toe die van toepassing is op alle waarden in het combined_raster. Dit is de terugval optie als geen van de door de gebruiker opgegeven beslisregels van toepassing is. Tenslotte bepaalt de MRT welke belisregel als eerste (dwz de meest linker) van toepassing was. Dit wordt vastgelegd in een nieuwe kolom first_rule
.
value | beheertypen | top10 | regel01 | regel02 | regel03 | regel04 | regel05 | regel06 | regel07 | regel999 | first_rule |
---|---|---|---|---|---|---|---|---|---|---|---|
1 | 10 | 160 | 0 | 0 | 0 | 0 | 1 | 1 | 1 | 1 | regel05 |
2 | 10 | 100 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | regel06 |
3 | 10 | 143 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 1 | regel06 |
4 | 11 | 142 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | regel01 |
5 | 11 | 143 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | regel02 |
6 | 11 | 144 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | regel03 |
7 | 11 | 100 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | regel04 |
8 | 12 | 161 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | regel05 |
9 | 13 | 144 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | regel999 |
De vertaling van het combined_raster naar het output raster verloopt dan volgens onderstaande tabel (waarbij 0
wordt aangehouden als newal
voor beslisregel #999):
combined_value | newval |
---|---|
1 | 14 |
2 | 14 |
3 | 14 |
4 | 10 |
5 | 11 |
6 | 12 |
7 | 13 |
8 | 14 |
9 | 0 |
De kerntaak van de MRT is om het combined-raster te vertalen volgens bovenstaande tabel en het resultaat weg te schrijven naar een nieuw rasterbestand. Het resulterende raster wordt aangeduid als het newval
raster. Daarnaast is de MRT in staat om een nieuw raster te maken waarin niet de newval
, maar de van-toepassing-zijnde beslisregel staat aangegeven. Dit is het firstrule
raster.
Bruto en netto areaal
De RAT van het combined_raster heeft ook een Count
kolom, waarin het aantal pixels met die pixelwaarde staan. Gegeven de oppervlakte van een pixel, kan deze informatie eenvoudig omgezet worden naar een areaal in hectare. De RAT in dit voorbeeld ziet er dan als volgt uit (veronderstelende een oppervlakte van 1 per pixel):
value | beheertypen | top10 | oppervlak |
---|---|---|---|
1 | 10 | 160 | 2 |
2 | 10 | 100 | 2 |
3 | 10 | 143 | 1 |
4 | 11 | 142 | 3 |
5 | 11 | 143 | 1 |
6 | 11 | 144 | 2 |
7 | 11 | 100 | 1 |
8 | 12 | 161 | 1 |
9 | 13 | 144 | 1 |
Met deze informatie kan het bruto en netto areaal berekend worden dat van toepassing is op elke beslisregel.
newval | description | remark | volgnr | Beheertypen | Top10 | bruto | netto | verschil |
---|---|---|---|---|---|---|---|---|
10 | Gemengd droog bos met productie | 1 | 11 | 142 | 3 | 3 | 0 | |
11 | Droog loofbos met productie | 2 | 11 | 143 | 1 | 1 | 0 | |
12 | Droog naalbos met productie | 3 | 11 | 144 | 1 | 1 | 0 | |
13 | Overig droog bos met productie | 4 | 11 | n | 7 | 1 | 6 | |
14 | Duinen | 5 | n | 160,161 | 3 | 3 | 0 | |
14 | Duinen | 6 | 10 | n | 5 | 3 | 2 | |
14 | Duinen | 7 | 10 | 160,161 | 2 | 0 | 2 | |
0 | overig | 999 | - | - | 14 | 1 | 13 |
Het bruto areaal is datgene waar de beslisregel in beginsel op van toepassing is, voordat beslisregels met hogere prioriteit worden toegepast. Het netto areaal is datgene waar de beslisregel uiteindelijk van toepassing is; de firstrule. Per definitie zijn bruto- en netto gelijk voor de beslisregel met hoogste prioriteit.
Netto arealen kunnen gesommeerd worden per nieuwe categorie:
newval | description | netto |
---|---|---|
10 | Gemengd droog bos met productie | 3 |
11 | Droog loofbos met productie | 1 |
12 | Droog naaldbos met productie | 1 |
13 | Overig Droog bos met productie | 1 |
14 | Duinen | 6 |
Downstream: doorvertaling van de nieuwe categorieen
Beslisregels resulteren in nieuwe categorieen. Het kan wenselijk zijn om die vervolgens te vertalen naar een andere categorisatie. Dit proces wordt aangeduid als 'downstream' (ie benedenstrooms, of afgeleid van de nieuwe kaart). Downstream categorisaties kunnen opgesteld worden voor zowel de newvals
van de nieuwe categorieen, of op de individuele beslisregels (de firstrule
). Hieronder worden beide situaties toegelicht.
volgens newval
In onderstaand voorbeeld worden de categorieen van belisregels #1-7 doorvertaald naar twee downstream-modellen: 'algemeen_landgebruik' en 'bostypen'.
newval | description | algemeen_landgebruik | bostype |
---|---|---|---|
10 | Gemengd droog bos met productie | Bos | Gemengd |
11 | Droog Loofbos met productie | Bos | Loof |
12 | Droog naaldbos met productie | Bos | Naald |
13 | Overig droog bos met productie | Bos | Gemengd |
14 | Duinen | Zand |
Downstream modellen hoeven niet elke newval
te adresseren: bostypen
blijft leeg voor newval
14. Dit resulteert in een NoData
.
De downstream regels staan in een apart tabblad van de Excel. Hierin staan doorgaans minder regels in dan in het upstream tabblad: elke newval wordt immers eenmalig benoemt in Downstream. Voor elke kolom in het downstream tabblad kan de MRT een kaart produceren.
In het Downstream tablad worden de categorie descriptions ingevuld in plaats van numerieke values (zoals in het upstream tabblad). De value die wordt toegekend aan elke categorie in een downstream-model (ie Gemengd
en Loof
voor bostypen
) kan gedefinieerd worden in een tabblad met de naam van het model:
Description | Value |
---|---|
Gemengd | 141 |
Loof | 142 |
Naald | 143 |
De gebruiker is vrij om nieuwe values te kiezen voor de downstream categorieen. In dit voorbeeld wordt voor bostypen
de values van Top10 gereconstrueerd. Categorieen in downstream modellen zonder bijbehorend tabblad worden automatisch genummerd vanaf 1, zoals in dit geval voor algemeen_landgebruik
.
Description | Value |
---|---|
Bos | 1 |
Zand | 2 |
De MRT kan voor elk downstream-model een nieuw raster wegschrijven. Dit gebeurt door de combined_value te koppelen aan de values van het downstream model. Voor bijvoorbeeld het bostype
model verloopt dat volgens onderstaande tabel:
value | bostype_value |
---|---|
1 | 0 |
2 | 0 |
3 | 0 |
4 | 141 |
5 | 142 |
6 | 143 |
7 | 141 |
8 | 0 |
9 | 0 |
volgens firstrule
Het kan wenselijk zijn om dezelfde newval anders te behandelen in de downstream, bijvoorbeeld om de beslisregels te groeperen naar een bepaalde gemeenschappelijke deler. Hiertoe kunnen downstream-modellen toegevoegd worden als extra kolommen aan de upstream sheet (deze kolomnamen moeten beginnen met een _
). In onderstaande tabel wordt een downstream model _oorsprong
gedefinieerd.
newval | description | remark | volgnr | Beheertypen | Top10 | _oorsprong |
---|---|---|---|---|---|---|
10 | Gemengd droog bos met productie | 1 | 11 | 142 | top10&SNL | |
11 | Droog loofbos met productie | 2 | 11 | 143 | top10&SNL | |
12 | Droog naalbos met productie | 3 | 11 | 144 | top10&SNL | |
13 | Overig droog bos met productie | 4 | 11 | n | SNL | |
14 | Duinen | 5 | n | 160,161 | top10 | |
14 | Duinen | 6 | 10 | n | SNL | |
14 | Duinen | 7 | 10 | 160,161 | top10&SNL |
Beslisregels #5-7 resulteren volgens bovenstaand schema in andere categorieen voor _oorsprong
, ondanks hun gedeelde newval
. Net als bij downstream volgens newval kunnen de numerieke waardes de categorieen in het _oorsprong
model gedefinieerd worden in een apart tabblad. Als dat niet gebeurt, kent de MRT automatisch nummers toe.
Voor het _oorsprong downstream-model kan de MRT een nieuw raster wegschrijven, gebaseerd op onderstaande tabel (veronderstellende dat de _oorsprong categorieen automatisch genummerd worden vanaf 1):
value | _oorsprong_value |
---|---|
1 | 1 |
2 | 1 |
3 | 1 |
4 | 1 |
5 | 1 |
6 | 1 |
7 | 2 |
8 | 3 |
9 | 0 |
Op dezelfde wijze zoals eerder beschreven, kan voor elke categorie in een downstream model het bruto- en netto areaal bepaald worden.
Toepassing van de MRT
Hieronder volgt instructie om de MRT te installeren, de benodigde bestanden te bereiden en om de MRT te gebruiken.
Installatie
- installeer git
- clone deze repository naar een lokale directory (bijv
c:/apps/mrt
) mbv git:git clone git@git.wur.nl:roelo008/mrt.git
. WUR-gebruikers kunnen een ssh key toevoegen aan hun wur-gitlab account. - installeer miniconda
- mbv miniconda creeer een Python environment met de volgende packages:
- open Anaconda Prompt en activeer de Python environment:
activate <environment_name>
. - navigeer naar de installatie-directory van de MRT:
cd c:/apps/mrt/
en vervolgens naar desrc
subdirectory: - verifieer dat MRT werk:
python MRT.py --test
Optioneel kan de MRT repository gedownload worden vanaf git.wur.nl, in dat geval is installatie van git niet nodig.
Creeer een combined-raster
Om daadwerkelijk met de MRT aan de slag te gaan, is de eerste stap het vezamelen van alle bronrasters. Deze moeten van integer data-type en identiek vormgegeven zijn:
- aantal rijen en kolommen
- pixel groote (ie resolutie)
- oorsprong in geprojecteerd coordinaten-referentiesysteem.
Een veelgebruikt format voor rasters binnen WEnR is:
- 2.5m resolutie
- 130.000 rijen, 112.000 kolommen
- extent in RD-New (epsg:28992):
- left: 0 m
- bottom: 300.000 m
- right: 280.000 m
- top: 625.000 m
Gebruik de ArcGIS Pro Combine Tool om alle bronrasters te combineren in een combined raster. Het output-format moet gdal
compatible zijn; in de praktijk is GeoTiff
het beste format. Pas LZW compressie toe om de bestandsgrootte te beperken. De output ziet er dan uit als:
-
combi_raster_v1.tif
: het raster zelf -
combi_raster_v1.tif.vat.dbf
: de RAT -
combi_raster_v1.tif.ovr
: overviews (optioneel)
Let op dat NoData
in een source-raster doorwerkt in het combi-raster. Daarom heeft het de voorkeur geen NoData
te hebben in de source-rasters zodat het combi-raster volledig gevuld is met geldige data. Source-rasters pixels met ongeldige data kunnen aangeduid worden met een 0
.
Opstellen beslisregels
De beslisregels en eventuele downstream doorvertalingen worden vastgelegd in een Excel document. De beslisregels staan in een apart tabblad met de volgende kolommen:
newval | description | remarks | reclass | volgnummer | <bronraster_1> | <bronraster_2> | <bronraster_n> |
---|
Verifieer dat de bronraster_ kolomnamen corresponderen met de kolomnamen in de combi-raster RAT. (merk op dat dat kolomnamen in een RAT beperkt zijn tot 10 karakters).
Maak een reclass-rule per regel in de excel:
-
newval
: het getal dat de nieuwe landgebruikscategorie vertegenwoordigt -
description
: de omschrijving van de nieuwe categorie -
remark
: ruimte voor opmerkingen -
reclass
staat op1
als de reclass-rule toegepast moet worden, of op0
als de reclass-rule genegeerd moet worden -
volgnummer
: nummer ter identificatie van de beslisregel. Wordt niet gebruikt in de MRT. - in elke source-raster kolom
- de waardes uit dat bron-raster die wel meedoen (hele getallen, komma gescheiden)
- OF welke waardes uit het source-raster niet meedoen (hele getallen met prefix
-
) - OF indifferentie naar dat source-raster, aangeduid met
n
Command Line Interface
De MRT wordt aangestuurd via de Command Line. Er zijn drie vereiste parameters en diverse optionele. De vereiste parameters dienen als eerste en in de juiste volgorde ingevoerd te worden.
-
xls
volledig adres naar de Excel sheet met de reclass rules (bijv.c:\docs\scenario1.xlsx
) -
rast
volledig adres naar het gecombineerde raster (bijv.c:\apps\geodata\combi\combiraster_v1.tif
) -
us
naam van het Excel tabblad waarin de reclass-rules staan
De optionele parameters zijn:
-
--help
laat help zien en sluit programma weer af -
--ds
naam van tabblad met downstream regels -
--test
verifieer of alle input in orde is, doet verder niks -
--report
maak rapport (Excel sheet) met alle informatie over de uitkomsten van de multi-Reclass. Rapporteert oa over bruto- en netto areaal van elke beslisregel. -
--write
schrijf de nieuwe landgebruik categorieën en/of een of meer downstream categorieën weg als een geospatial raster. Meerdere waardes mogelijk. -
--xtnt
ruimtelijke extent van het output raster. Keuze uitBNL
(0, 300.000, 280.000, 625.000),LCEU
(10.000, 300.000, 280.000, 620.000) ofsrc
(gelijk aan de extent vanrast
). DefaultBNL
. Deze optie is momenteel niet goed geimplementeerd, beter niet gebruiken! -
--of
output bestandstype. Keuze uitgeotiff
,flt
ofmap
. Defaultgeotiff
. -
--base_out
directory waar de output aangemaakt wordt. -
--verbose
True
/False
, geef meer feedback tijdens draaien van de tool. DefaultFalse
. -
--nodata
specify which value in the output raster should be tagged asnodata
. By default, none of the output pixel values are tagged asnodata
. -
--mnp
Schrijf een *.csv legenda voor het output raster. Bedoeld voor neergeschaalde beheertypen. Dit is nodig om de kaart te koppelen aan MNP. Hiervoor is de benb_utils module vereist. -
--ovr
create overviews when writing new rasters. Default False
Onderstaand voorbeeld draait de MRT en creëert een rapport en een nieuw Geotiff raster met LCEU extent met daarin de nieuwe landgebruik categorieën (aangeduid als 'newval').
python MRT.py c:/apps/scenario_reclass.xlsx c:/apps/combinatie_raster.tif upstream_sheet --report --write newval --xtnt LCEU --base_out c:/apps/output
Onderstaand voorbeeld creëert 3 output rasters. Merk op hoe de --xtnt
en --nodata
parameters voor elk van de drie output rasters gespecificeerd wordt. Gebruik X
om terug te vallen op de default waarde van de parameter.
python MRT.py c:/apps/scenario_reclass.xlsx c:/apps/combinatie_raster.tif upstream_sheet --write newval algemeen_landgebruik Bostype --xtnt LCEU src BNL --nodata X 0 0
In onderstaand voorbeeld wordt de --mnp
functie aangeroepen. Hiermee wordt een *.csv legenda gemaakt voor het output raster. Deze optie is specifiek bedoeld voor (neergeschaalde) beheertypen, maar werkt in principe ook voor andere downstream modellen.
python mrt.py c:/apps/scenario_reclass.xlsx c:/apps/combinatie_raster.tif upstream_sheet --ds downstream_sheet --write neergeschaalde_beheertypen --mnp neergeschaalde_beheertypen
Potentiele fouten in de Excel worden ondervangen en teruggelegd aan de gebruiker. Als er geen downstream sheet wordt gespecificeerd voor de MRT, kunnen enkel de newval en firstrule modellen weggeschreven als nieuw raster.
MRT Output
Hieronder worden de producten omschreven die de MRT kan produceren. De daadwerkelijke output van de de MRT is afhankelijk van de opdracht van de gebruiker.
Rasters
De MRT kan nieuwe rasterbestanden produceren. In het geval alleen een upstream
sheet is opgegeven, kan de MRT het newval
of het firstrule
raster produceren. Als de beslisregels zijn uitgebreid met downstream modellen volgens de firstrule
, dan kunnen deze downstream modellen ook geproduceerd worden. De rasters die de MRT produceert bestaan uit de volgende bestanden:
- <basename>_<model_name>.tif: het GeoTiff bestand met het daadwerkelijke raster
- <basename>_<model_name>.qml: QGIS layer style file. In dit bestand staan de legenda en kleurstelling van de nieuwe categorieen Als het *.tif bestand in QGIS wordt geopend, wordt dit bestand automatisch herkend en toegepast.
- <basename>_<model_name>.tif.vat.dbf: de raster attribute table (RAT). Dit is een digitale tabel waarin de koppeling tussen de rastervalue en de categorie is vastgelegd. Dit bestand wordt automatisch herkend door ArcGIS Pro en door QGIS mits de juiste plugin is geïnstalleerd
De basename van MRT producten is gelijk aan de de bestandsnaam van de Excel met de beslisregels, gevolgd door de datum & tijdstip waarop dit bestand voor het laatst gewijzigd is. model_name verwijst naar het 'model' dat de MRT heeft geproduceerd. Dit kan de newval
of firstrule
zijn, of een van de downstream modellen.
Rapport
De MRT produceert standaard een Excel rapport als het newval
model wordt weggeschreven, in andere gevallen is het rapport optioneel. Het rapport kan ook gemaakt worden als er geen rasters worden geproduceerd. Het rapport bevat de volgende tabbladen.
metadata
Hierin zijn alle details van de MRT run vastgelegd, zoals de gebruikte bronbestanden en de datum-tijdstip van de run.
ReclassRules
De beslisregels worden hier nogmaals genoemd, zoals ze zijn opgegeven door de gebruiker. Daarnaast zijn extra kolommen toegegoegd:
- mrt_query: de formele database query die wordt toegepast op de combined_raster RAT
- name: naam van de beslisregel
- priority: numerieke indicatie van de prioriteit van de beslisregel, waarbij 0 de hoogste prioriteit heeft
- gross_match_ha: bruto areaal in hectare van het combined_raster dat voldoet aan de beslisregel.
- net_match_ha: netto areaal in hectare van het combiend_raster dat voldoet aan de beslisregel.
- difference: verschil tussen bruto en netto areaal.
Het bruto areaal van een beslisregel slaat op het areaal dat voldoet aan de beslisregel, voordat hogere-prioriteit beslisregels worden toegepast. Het netto-areaal is het areaal dat resteert na toepassing van hogere-prioriteit beslisregels. Het netto areaal van een beslisregel is daarom nooit kleiner dan het bruto areaal.
De bruto, netto en verschil kolommen helpen de gebruiker te verifieren dat de beslisregels naar verwachting werken. Een netto-areaal van 0 hectare kan bijvoorbeeld betekenen dat de beslisregel te laag ge-prioriteerd is, of te 'streng' is.
Reclassed
In dit tabblad staat de RAT van het combined_raster. Voor elke waarde is aangegeven aan welke newval
het is toegekend. Wegens beperkingen van Excel wordt dit tabblad alleen aangemaakt als de RAT niet meer dan 1.048.576 rijen bevat.
NotClassified
In dit tabblad staan de values uit de verschillende bronrasters die niet gebruikt zijn in een beslisregel.
source_raster
Voor elk bronraster wordt een apart tabblad aangemaakt waarin de values genoemd worden, plus areaal daarvan.
mnp_tabel
Optioneel maakt de MRT een *.csv bestand waarin de legenda van het geproduceerde raster staat. De legenda is opgemaakt volgens de vereisten van de land type raster legenda van het MNP model. Deze optie is expliciet bedoeld voor rasters met SNL Beheertypen.
Meer weten?
Neem contact op met Hans Roelofsen.