djBdZddlZddlmZdZdZdZdZ d d Z d Z y) u`Ohuet kääreet GDAL-kutsuille. Kaikki GDAL-riippuvuus on tässä moduulissa, jotta terrarium.py ja tilemath.py pysyvät testattavina ilman GDAL-asennusta. KORKEUSJÄRJESTELMÄ (tärkeä): MML:n korkeusmalli on N2000-korkeuksia. EPSG:3067 on 2D-koordinaatisto, joten gdal.Warp muuntaa vain tasokoordinaatit (3067 -> 3857) eikä PROJ yritä pystymuunnosta — pikseliarvot (korkeudet) kulkevat läpi koskemattomina. Tämä on varmistettu testillä tests/test_pipeline_e2e.py::test_heights_survive_warp. Jos lähtöaineistoon joskus ilmestyy 3D-CRS (compound CRS), tämä oletus pitää tarkistaa uudelleen. N)tilemathgcv ddlm}m}|j ||fS#t$r}t d|d}~wwxYw)u=Palauttaa (gdal, osr) tai kaatuu selkeällä asennusohjeella.r)gdalosrzGDAL:n Python-sidokset puuttuvat. Asenna esim.: pip install gdal (viralliset wheelit, GDAL >= 3.9) tai conda-forge: conda install -c conda-forge gdalN)osgeorr ImportError SystemExit UseExceptions)rres )/home/claude/oh9ab/pipeline/gdal_utils.py import_gdalrsL#  9  B   s 8 38c t\}}|jt||Dcgc] }t|c}}|td|z|j |Scc}w)usgdalbuildvrt: yhdistä karttalehdet virtuaalirasteriksi (ei mosaikoi fyysisesti — VRT on vain viittauslista).ugdal.BuildVRT epäonnistui: %s)rBuildVRTstr RuntimeError FlushCache)inputsout_pathr_pvrts r build_vrtr's]mGD! --H '?1A'? @C {;hFGGNN O (@sA" c t\}}|jt|}|j}|j|j }}|j |j}|j }|jd|j|j|j|j|j||} g} d} t| dzD]6} | | z } | | |zdf| |zt|fd| |zft|| |zfgz } 8| Dcgc]3\}}|d||dzz||dzz|d||dzz||d zzf5}}}|Dcgc]\}}| j||d d}}}|Dcgc]}|d }}|Dcgc]}|d }}t|t|t!|t!|fScc}}wcc}}wcc}wcc}w) uLähderasterin kattavuus WGS84 lon/lat -laatikkona (minlon, minlat, maxlon, maxlat). Muunnetaan rasterin reunojen pisteistö (ei vain kulmat), koska projektiomuunnoksessa laatikon ääriarvot voivat osua reunan keskelle. irrN)rOpenrGetGeoTransform RasterXSize RasterYSizeSpatialReference GetProjectionImportFromEPSGSetAxisMappingStrategyOAMS_TRADITIONAL_GIS_ORDERCoordinateTransformationrangefloatTransformPointminmax)src_pathrrdsgtwhsrcdsttredgenitcrptsxyllrlonslatss r source_bounds_lonlatrD2s ID# 3x= !B   B >>2>>qA  r//1 2C   Cts==>s==> % %c3 /B D A 1q5\S E !a%AuQx03A,q1q5@QRRS q! qEA1I BqE )2a51r!u9+NNOO Jc*tj}||z |z}||z |z}|jdj||||}|j t j }|tdzk} | jryt j|| <|S)uLue yhden ruudun 256x256-ikkuna warpatusta VRT:stä. Palauttaa float64-taulukon, jossa nodata on NaN, tai None jos ruutu on kokonaan nodataa (jätetään kirjoittamatta). rg?N) r TILE_SIZE GetRasterBand ReadAsArrayastypenpfloat64rTallnan) warped_dstile_xtile_yrWrXtscolrowarrmasks r read_tile_windowrtzs   B G r !C G r !C  ! !! $ 0 0c2r BC **RZZ C &3, D xxzCI Jrb)bilinear) __doc__numpyrhrFrrTrrrDrartrbr rys9   6F",#Lrb