dj dZddlZdZdejdejfdZdejdejfdZdejd ejdejfd Zdejfd Zy) uLTerrarium-korkeuskoodaus (Mapzen/AWS elevation-tiles -yhteensopiva). Purkukaava (sama kuin selainprototyypissä, kuuluvuus.html): h = (R*256 + G + B/256) - 32768 [m] eli 24-bittinen kokonaisluku v = R*65536 + G*256 + B tulkitaan arvoksi v/256 - 32768. Askel on 1/256 m ~ 3,9 mm, mikä riittää 2 m korkeusmallille reilusti. Terrarium-formaatissa ei ole nodata-arvoa: nodata/meri koodataan korkeudeksi 0 m (N2000-nolla on likimain merenpinta, joten tämä on fysikaalisesti järkevä täyte). Metsäruuduille (Luken MVMI) käytetään omaa kanavapakkausta, ks. encode_forest(). Ng@height_mreturnctj|tj}tjtj||d}tj |t zdz}tj|ddjtj}tj|jdztj}|dz dz|d <|d z dz|d <|dz|d <|S) uKorkeus metreinä (float, NaN = nodata) -> (H, W, 3) uint8 RGB. Kvantisointi pyöristää lähimpään 1/256 m askeleeseen, joten dekoodausvirhe on enintään 1/512 m. dtypep@ri.r..) npasarrayfloat64whereisfiniterint_OFFSETclipastypeuint32emptyshapeuint8)rhvrgbs (/home/claude/oh9ab/pipeline/terrarium.pyencoder%s 82::.A QC(A W%&A 1l#**2995A ((177T> 2C7d"CK6T/CKd(CK Jr#c|djtj}|djtj}|djtj}|dz|z|dz ztz S)u3(H, W, 3) uint8 RGB -> korkeus metreinä (float64).rrrr )rrrr)r#rgbs r$decoder+'sf F 2::&A F 2::&A F 2::&A u9q=1u9 $w ..r& cover_pctctj|tj}tj|tj}tjtj||d}tjtj||d}tj |j dztj}tjtj|ddjtj|d<tjtj|ddjtj|d<|S) uPuustoruudun kanavapakkaus: R = puuston keskipituus metreinä (0..255, kokonaisluku) G = latvuspeittävyys prosentteina (0..100) B = varattu tulevalle käytölle (nyt 0) NaN (nodata, esim. vesistö tai ei-metsämaa) -> 0 molemmissa kanavissa. Huom: Luken MVMI ilmoittaa keskipituuden desimetreinä — muunnos metreiksi tehdään ennen tätä funktiota (build_tiles.py: --height-scale). rrr rr rdr) rrrrrzerosrr rrr)rr,r!cr#s r$ encode_forestr1/s 82::.A 9BJJ/A QC(A QC(A ((177T> 2C''"''!*a-44RXX>CK''"''!*a-44RXX>CK Jr&c|djtj|djtjfS)uD(H, W, 3) uint8 -> (keskipituus_m, peittävyys_pct) float64-pareina.rr)rrr)r#s r$ decode_forestr3Ds3 v;  bjj )3v;+=+=bjj+I IIr&) __doc__numpyrrndarrayr%r+r1r3r&r$r8s  RZZBJJ"/ /rzz/BJJ2::"***JrzzJr&