Topografinen kosteusindeksi (Beven & Kirkby, 1979) on laskettu koko Suomelle 16 metrisessä hilassa, joka on yhtenevä Luonnonvarakeskuksen monilähteisen kansallisen metsäinventaarion kartta-aineiston kanssa (Luke MS-NFI, 2015). Laskennassa hyödynnettiin resoluutioltaan 2 m 2 m pintamallia (kesäkuun 2016 saatavuuden mukaan) ja 10 m 10 m pintamallia muilta osin (Maanmittauslaitos 2017). Pintamalli ositettiin 500 metrin puskurilla SYKEn kolmannen jakovaiheen valuma-alueiden rajojen mukaan (SYKE, 2010) ja aggregoitiin keskiarvolla 16 metrin hilaruutuihin. Pintamallin esikäsittelyt ja kosteusindeksin laskenta toteutettiin valuma-alueittain python-koodilla Whitebox GIS-ohjelmistossa (Lindsay, 2014). Teiden ja uomien risteyskohdat muokattiin pintamalliin ensin, ja sen jälkeen uomaelementit poltettiin pintamalliin 1m maksimiarvolla niin että huomioitiin mahdolliset erimielisyydet maastotietokannan uomien sijainnin ja pintamallin välillä vaimennuskertoimella. Piantamlli esikäsiteltiin edelleen kaiverrusmenetelmällä (’fast breaching’), joka täyttämisen sijaan käsittelee mahdolliset veden virtauksen estävät ja pintamallin epätarkkuudesta johtuvat väärät kuoppakohdat kaivertamalla sopivan virtausreitin pintamalliin (Lindsay, 2016). Esikäsittelyn jälkeen virtaussuunnat ja virtaaman kerääntymisrasterit laskettiin ’D-infinity’ -menetelmällä (Tarboton, 1997), ja kaltevuus määritettiin alkuperäisestä muokkaamattomasta pintamallista. Topografinen kosteusindeksi laskettiin ja aineiston siirtämisen sujuvoittamiseksi TWI-indeksin arvot kerrottiin 1000:lla jotta arvot voitiin muuttaa kokonaisluvuiksi indeksin tarkkuuden vähentymättä.
Hila: 16*16m
Koordinaattijärjestelmä: ETRS89 / ETRS-TM35FIN (EPSG:3067)