Difference between revisions of "Höjdkalibering vid beståndsregisterimport"
| Line 2: | Line 2: | ||
Vid import av beståndsregister simulerar programmet träddata baserat på diameterfördelningsfunktioner. För varje träd som skapas beräknas en trädhöjd med hjälp av höjdfunktioner. Dessa funktioner (<ref name = "ClaessonEtAl2001">{{:Bibliografi:Söderberg1992}}</ref>) baseras på provträdsdata från Riksskogstaxeringen och täcker alla regioner i Sverige. För specifika bestånd kan de verkliga trädhöjderna naturligtvis avvika från de beräknade. Därför finns möjlighet att låta systemet kalibrera de beräknade höjderna så att medelhöjden för de skapade träden motsvarar den faktiska medelhöjden. Nyttan med höjdfunktionen är då att den fördelar trädhöjderna på träden så att inte alla träd får samma höjd. | Vid import av beståndsregister simulerar programmet träddata baserat på diameterfördelningsfunktioner. För varje träd som skapas beräknas en trädhöjd med hjälp av höjdfunktioner. Dessa funktioner (<ref name = "ClaessonEtAl2001">{{:Bibliografi:Söderberg1992}}</ref>) baseras på provträdsdata från Riksskogstaxeringen och täcker alla regioner i Sverige. För specifika bestånd kan de verkliga trädhöjderna naturligtvis avvika från de beräknade. Därför finns möjlighet att låta systemet kalibrera de beräknade höjderna så att medelhöjden för de skapade träden motsvarar den faktiska medelhöjden. Nyttan med höjdfunktionen är då att den fördelar trädhöjderna på träden så att inte alla träd får samma höjd. | ||
| − | Kalibreringen görs i två steg, först trädslagsvis om trädslagsvisa medelhöjder finns i importfilen, och sedan för beståndet i sin helhet om en medelhöjd för beståndet har angivits. | + | Kalibreringen görs i två steg, först trädslagsvis om trädslagsvisa medelhöjder finns i importfilen, och sedan för beståndet i sin helhet om en medelhöjd för beståndet har angivits. Kalibreringen görs separat för olika trädslag och skikt (huvudskikt och överståndare). |
==Kalibreringskvoten räknas i princip ut på följande sätt (läs om viktning längre ner).== | ==Kalibreringskvoten räknas i princip ut på följande sätt (läs om viktning längre ner).== | ||
| − | ====Steg 1 - Trädslagsvis för de trädslag som medelhöjd har angivits i registerfilen==== | + | ====Steg 1 - Trädslagsvis för de trädslag som medelhöjd har angivits för i registerfilen==== |
Kalibreringskvoten för trädslag s beräknas som | Kalibreringskvoten för trädslag s beräknas som | ||
<math> | <math> | ||
| − | {q(s)} = \frac{H_{ | + | {q(s,k)} = \frac{H_{register}(s,k)} {H_{funktion}(s,k)} |
</math> | </math> | ||
där | där | ||
| − | <br>''H''<sub> | + | <br>''H''<sub>register</sub> (''s'') = Medelhöjd for trädslag ''s'' (och skikt ''k'') enligt register. |
| − | <br>''H''<sub> | + | <br>''H''<sub>funktion</sub> (''s'') = Grundytevägt medelvärde för trädens beräknade trädhöjder för träd av trädslag ''s'' (coh skikt ''k''). |
| − | Beräknade höjder för trädslag s justeras (multipliceras) sedan med faktorn q(s). T ex om beräknad medelhöjd är 10 m och höjden enligt registret är 12 m, blir kvoten 12/10 = 1.2, vilket betyder att trädens beräknade höjder justeras upp med 20 %. | + | Beräknade höjder för trädslag s justeras (multipliceras) sedan med faktorn q(s,k). T ex om beräknad medelhöjd är 10 m och höjden enligt registret är 12 m, blir kvoten 12/10 = 1.2, vilket betyder att trädens beräknade höjder justeras upp med 20 %. |
====Steg 2 - Kalibrering på beståndsnivå.==== | ====Steg 2 - Kalibrering på beståndsnivå.==== | ||
<br>I detta steg beräknas motsvarande kvot som i steg 1, men för alla träd oberoende av trädslag. Det är de uppdaterade trädhöjderna från justeringen i steg 1 som används vid beräkningen av H<sub>calc</sub>. | <br>I detta steg beräknas motsvarande kvot som i steg 1, men för alla träd oberoende av trädslag. Det är de uppdaterade trädhöjderna från justeringen i steg 1 som används vid beräkningen av H<sub>calc</sub>. | ||
| + | <math> | ||
| + | {q(s)} = \frac{H_{register}(s)} {H_{funktion}(s)} | ||
| + | </math> | ||
| + | |||
| + | där | ||
| + | <br>''H''<sub>register</sub>(''k'') = Medelhöjd för skikt ''k'' i beståndet enligt register. | ||
| + | <br>''H''<sub>funktion</sub> (''k'') = Grundytevägt medelvärde för trädens beräknade trädhöjder för alla skapade träd som tillhör skikt ''k''. | ||
Revision as of 13:11, 28 August 2014
Vid import av beståndsregister simulerar programmet träddata baserat på diameterfördelningsfunktioner. För varje träd som skapas beräknas en trädhöjd med hjälp av höjdfunktioner. Dessa funktioner ([1]) baseras på provträdsdata från Riksskogstaxeringen och täcker alla regioner i Sverige. För specifika bestånd kan de verkliga trädhöjderna naturligtvis avvika från de beräknade. Därför finns möjlighet att låta systemet kalibrera de beräknade höjderna så att medelhöjden för de skapade träden motsvarar den faktiska medelhöjden. Nyttan med höjdfunktionen är då att den fördelar trädhöjderna på träden så att inte alla träd får samma höjd.
Kalibreringen görs i två steg, först trädslagsvis om trädslagsvisa medelhöjder finns i importfilen, och sedan för beståndet i sin helhet om en medelhöjd för beståndet har angivits. Kalibreringen görs separat för olika trädslag och skikt (huvudskikt och överståndare).
Kalibreringskvoten räknas i princip ut på följande sätt (läs om viktning längre ner).
Steg 1 - Trädslagsvis för de trädslag som medelhöjd har angivits för i registerfilen
Kalibreringskvoten för trädslag s beräknas som
där
Hregister (s) = Medelhöjd for trädslag s (och skikt k) enligt register.
Hfunktion (s) = Grundytevägt medelvärde för trädens beräknade trädhöjder för träd av trädslag s (coh skikt k).
Beräknade höjder för trädslag s justeras (multipliceras) sedan med faktorn q(s,k). T ex om beräknad medelhöjd är 10 m och höjden enligt registret är 12 m, blir kvoten 12/10 = 1.2, vilket betyder att trädens beräknade höjder justeras upp med 20 %.
Steg 2 - Kalibrering på beståndsnivå.
I detta steg beräknas motsvarande kvot som i steg 1, men för alla träd oberoende av trädslag. Det är de uppdaterade trädhöjderna från justeringen i steg 1 som används vid beräkningen av Hcalc.
där
Hregister(k) = Medelhöjd för skikt k i beståndet enligt register.
Hfunktion (k) = Grundytevägt medelvärde för trädens beräknade trädhöjder för alla skapade träd som tillhör skikt k.
Då funktionerna är tämligen robusta, medan registerdata inte sällan innehåller data av dålig kvalitet, finns det en inbyggd begräsning för hur mycket trädhöjderna kalibreras. Förvalda värden är max 30 % upp eller ned. Det betyder att kvoten som beräknas sätts till som lägst 0.7 och som högst 1.3. T ex, om en medelhöjd enligt registret är 14 m, men enligt höjdfunktionerna är 10 m, blir kvoten 14/10 = 1.4. Eftersom denna är större än 1.3, sätts kvoten till 1.3. De tillåtna avvikelserna kan ändras av användaren i kontrolltabellen ProductionModel.
- ↑ Söderberg, U. 1992. Funktioner för skogsindelning - Höjd, formhöjd och barktjocklek för enskilda träd (Functions for forest management - Height, form height and bark thickness of individual trees). Institutionen for skogstaxering, SLU, Umeå. ISNN 0348-0496 (in Swedish with English summary).