umep workflow

import%20marimo%0A%0A__generated_with%20%3D%20%220.17.6%22%0Aapp%20%3D%20marimo.App()%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Example%3A%20Complete%20UMEP%20workflow%20using%20pymdurs%20and%20solweig%0A%0A%20%20%20%20This%20example%20demonstrates%20how%20to%3A%0A%20%20%20%201.%20Collect%20urban%20data%20using%20pymdurs%20(DEM%2C%20buildi%20%20%20%20ngs%2C%20vegetation)%0A%20%20%20%202.%20Process%20data%20for%20UMEP%20analysis%20(DSM%2C%20CDSM)%0A%20%20%20%203.%20Calculate%20Sky%20View%20Factor%20(SVF)%20using%20solweig%0A%20%20%20%204.%20Optionally%20run%20SOLWEIG%20for%20thermal%20comfort%20analysis%0A%0A%20%20%20%20Inspired%20by%3A%20https%3A%2F%2Fgithub.com%2FUMEP-dev%2Fumep-rust%2Fblob%2Fmain%2Fdemos%2Fathens-demo.py%0A%0A%20%20%20%20Required%20dependencies%20(install%20separately)%3A%0A%20%20%20%20%60%60%60bash%0A%20%20%20%20pip%20install%20geopandas%20rasterio%20pyproj%0A%20%20%20%20%60%60%60%0A%20%20%20%20%23%23%23%20umep%20(optional)%20for%20additional%20SOLWEIG%20features%0A%20%20%20%20%60%60%60bash%0A%20%20%20%20pip%20install%20%22solweig%20%40%20git%2Bhttps%3A%2F%2Fgithub.com%2FUMEP-dev%2Fsolweig.git%40main%22%0A%20%20%20%20%60%60%60%0A%20%20%20%20Note%3A%20On%20Apple%20Silicon%20(ARM64)%2C%20solweig%20may%20require%20the%20x86_64%20target%3A%0A%20%20%20%20%60%60%60bash%0A%20%20%20%20rustup%20target%20add%20x86_64-apple-darwin%0A%20%20%20%20%60%60%60%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20os%0A%20%20%20%20from%20pathlib%20import%20Path%0A%0A%20%20%20%20import%20geopandas%20as%20gpd%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20rasterio%0A%20%20%20%20from%20osgeo%20import%20gdal%2C%20gdalconst%0A%20%20%20%20from%20rasterio.features%20import%20rasterize%2C%20shapes%0A%20%20%20%20from%20rasterio.transform%20import%20from_bounds%0A%20%20%20%20from%20shapely.geometry%20import%20shape%0A%0A%20%20%20%20import%20pymdurs%0A%0A%20%20%20%20%23%20Try%20to%20import%20solweig%20for%20SVF%20calculation%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20from%20umepr%20import%20solweig_runner_rust%2C%20svf%0A%0A%20%20%20%20%20%20%20%20HAS_UMEPR%20%3D%20True%0A%20%20%20%20except%20ImportError%3A%0A%20%20%20%20%20%20%20%20HAS_UMEPR%20%3D%20False%0A%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20umepr%20package%20not%20available.%20SVF%20calculation%20will%20be%20skipped.%22)%0A%20%20%20%20%20%20%20%20print(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20%20Install%20with%3A%20pip%20install%20'solweig%20%40%20git%2Bhttps%3A%2F%2Fgithub.com%2FUMEP-dev%2Fsolweig.git%40main'%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20print(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20%20On%20Apple%20Silicon%2C%20you%20may%20need%3A%20rustup%20target%20add%20x86_64-apple-darwin%22%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%23%20Try%20to%20import%20umep%20for%20additional%20utilities%20(optional)%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20from%20umep%20import%20wall_heightaspect_algorithm%20%20%23%20noqa%3A%20F401%0A%20%20%20%20%20%20%20%20from%20umep.functions.SOLWEIGpython%20import%20solweig_runner_core%0A%0A%20%20%20%20%20%20%20%20HAS_UMEP%20%3D%20True%0A%20%20%20%20except%20ImportError%3A%0A%20%20%20%20%20%20%20%20HAS_UMEP%20%3D%20False%0A%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20umep%20package%20not%20available.%20Some%20features%20will%20be%20disabled.%22)%0A%0A%20%20%20%20print(%22HAS_UMEPR%22%2C%20HAS_UMEPR)%0A%20%20%20%20print(%22HAS_UMEP%22%2C%20HAS_UMEP)%0A%0A%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%23%20COSIA%20Color%20to%20Class%20Mapping%0A%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20TABLE_COLOR_COSIA%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22B%C3%A2timent%22%3A%20%22%23ce7079%22%2C%0A%20%20%20%20%20%20%20%20%22Zone%20imperm%C3%A9able%22%3A%20%22%23a6aab7%22%2C%0A%20%20%20%20%20%20%20%20%22Zone%20perm%C3%A9able%22%3A%20%22%23987752%22%2C%0A%20%20%20%20%20%20%20%20%22Piscine%22%3A%20%22%2362d0ff%22%2C%0A%20%20%20%20%20%20%20%20%22Serre%22%3A%20%22%23b9e2d4%22%2C%0A%20%20%20%20%20%20%20%20%22Sol%20nu%22%3A%20%22%23bbb096%22%2C%0A%20%20%20%20%20%20%20%20%22Surface%20eau%22%3A%20%22%233375a1%22%2C%0A%20%20%20%20%20%20%20%20%22Neige%22%3A%20%22%23e9effe%22%2C%0A%20%20%20%20%20%20%20%20%22Conif%C3%A8re%22%3A%20%22%23216e2e%22%2C%0A%20%20%20%20%20%20%20%20%22Feuillu%22%3A%20%22%234c9129%22%2C%0A%20%20%20%20%20%20%20%20%22Coupe%22%3A%20%22%23e48e4d%22%2C%0A%20%20%20%20%20%20%20%20%22Broussaille%22%3A%20%22%23b5c335%22%2C%0A%20%20%20%20%20%20%20%20%22Pelouse%22%3A%20%22%238cd76a%22%2C%0A%20%20%20%20%20%20%20%20%22Culture%22%3A%20%22%23decf55%22%2C%0A%20%20%20%20%20%20%20%20%22Terre%20labour%C3%A9e%22%3A%20%22%23d0a349%22%2C%0A%20%20%20%20%20%20%20%20%22Vigne%22%3A%20%22%23b08290%22%2C%0A%20%20%20%20%20%20%20%20%22Autre%22%3A%20%22%23222222%22%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20%23%20COSIA%20to%20UMEP%20classification%20mapping%0A%20%20%20%20COSIA_TO_UMEP%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22B%C3%A2timent%22%3A%202%2C%20%20%23%20Building%0A%20%20%20%20%20%20%20%20%22Zone%20imperm%C3%A9able%22%3A%201%2C%20%20%23%20Paved%0A%20%20%20%20%20%20%20%20%22Zone%20perm%C3%A9able%22%3A%206%2C%20%20%23%20Bare%20Soil%0A%20%20%20%20%20%20%20%20%22Piscine%22%3A%207%2C%20%20%23%20Water%0A%20%20%20%20%20%20%20%20%22Serre%22%3A%201%2C%20%20%23%20Paved%0A%20%20%20%20%20%20%20%20%22Sol%20nu%22%3A%206%2C%20%20%23%20Bare%20Soil%0A%20%20%20%20%20%20%20%20%22Surface%20eau%22%3A%207%2C%20%20%23%20Water%0A%20%20%20%20%20%20%20%20%22Neige%22%3A%207%2C%20%20%23%20Water%0A%20%20%20%20%20%20%20%20%22Conif%C3%A8re%22%3A%206%2C%20%20%23%20Bare%20Soil%20(trees%20not%20in%20UMEP%2C%20mapped%20to%20soil)%0A%20%20%20%20%20%20%20%20%22Feuillu%22%3A%206%2C%20%20%23%20Bare%20Soil%20(trees%20not%20in%20UMEP%2C%20mapped%20to%20soil)%0A%20%20%20%20%20%20%20%20%22Coupe%22%3A%205%2C%20%20%23%20Grass%0A%20%20%20%20%20%20%20%20%22Broussaille%22%3A%205%2C%20%20%23%20Grass%0A%20%20%20%20%20%20%20%20%22Pelouse%22%3A%205%2C%20%20%23%20Grass%0A%20%20%20%20%20%20%20%20%22Culture%22%3A%205%2C%20%20%23%20Grass%0A%20%20%20%20%20%20%20%20%22Terre%20labour%C3%A9e%22%3A%206%2C%20%20%23%20Bare%20Soil%0A%20%20%20%20%20%20%20%20%22Vigne%22%3A%205%2C%20%20%23%20Grass%0A%20%20%20%20%20%20%20%20%22Autre%22%3A%201%2C%20%20%23%20Paved%0A%20%20%20%20%7D%0A%0A%20%20%20%20%23%20UMEP%20labels%0A%20%20%20%20UMEP_LABELS%20%3D%20%7B%0A%20%20%20%20%20%20%20%201%3A%20%22Paved%22%2C%0A%20%20%20%20%20%20%20%202%3A%20%22Building%22%2C%0A%20%20%20%20%20%20%20%203%3A%20%22Evergreen%20Trees%22%2C%0A%20%20%20%20%20%20%20%204%3A%20%22Deciduous%20Trees%22%2C%0A%20%20%20%20%20%20%20%205%3A%20%22Grass%22%2C%0A%20%20%20%20%20%20%20%206%3A%20%22Bare%20Soil%22%2C%0A%20%20%20%20%20%20%20%207%3A%20%22Water%22%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20def%20hex_to_rgb(hex_color)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Convert%20hex%20color%20to%20RGB%20tuple.%22%22%22%0A%20%20%20%20%20%20%20%20hex_color%20%3D%20hex_color.lstrip(%22%23%22)%0A%20%20%20%20%20%20%20%20return%20tuple(int(hex_color%5Bi%20%3A%20i%20%2B%202%5D%2C%2016)%20for%20i%20in%20(0%2C%202%2C%204))%0A%0A%20%20%20%20def%20geodataframe_to_tif_with_metadata(%0A%20%20%20%20%20%20%20%20gdf%3A%20gpd.GeoDataFrame%2C%0A%20%20%20%20%20%20%20%20output_tif%3A%20str%2C%0A%20%20%20%20%20%20%20%20column%3A%20str%20%3D%20%22type%22%2C%0A%20%20%20%20%20%20%20%20resolution%3A%20float%20%3D%201.0%2C%0A%20%20%20%20)%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Convert%20a%20GeoDataFrame%20to%20TIF%20with%20classification%20metadata.%0A%0A%20%20%20%20%20%20%20%20Args%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20gdf%3A%20GeoDataFrame%20with%20geometries%20and%20classification%20column%0A%20%20%20%20%20%20%20%20%20%20%20%20output_tif%3A%20Output%20GeoTIFF%20path%0A%20%20%20%20%20%20%20%20%20%20%20%20column%3A%20Column%20name%20containing%20classification%20values%0A%20%20%20%20%20%20%20%20%20%20%20%20resolution%3A%20Pixel%20resolution%20in%20meters%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20print(%22%5Cn%F0%9F%93%8A%20Conversion%20du%20GeoDataFrame%20en%20TIF...%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Colonne%3A%20%7Bcolumn%7D%2C%20R%C3%A9solution%3A%20%7Bresolution%7D%20m%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Validate%20GeoDataFrame%0A%20%20%20%20%20%20%20%20if%20len(gdf)%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20ValueError(%22Le%20GeoDataFrame%20est%20vide%2C%20impossible%20de%20cr%C3%A9er%20un%20raster%22)%0A%0A%20%20%20%20%20%20%20%20bounds%20%3D%20gdf.total_bounds%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Bounds%3A%20%7Bbounds%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Calculate%20dimensions%0A%20%20%20%20%20%20%20%20width%20%3D%20int((bounds%5B2%5D%20-%20bounds%5B0%5D)%20%2F%20resolution)%0A%20%20%20%20%20%20%20%20height%20%3D%20int((bounds%5B3%5D%20-%20bounds%5B1%5D)%20%2F%20resolution)%0A%0A%20%20%20%20%20%20%20%20%23%20Validate%20dimensions%0A%20%20%20%20%20%20%20%20if%20width%20%3C%3D%200%20or%20height%20%3C%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20ValueError(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22Dimensions%20invalides%3A%20width%3D%7Bwidth%7D%2C%20height%3D%7Bheight%7D.%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22Bounds%3A%20%7Bbounds%7D%2C%20R%C3%A9solution%3A%20%7Bresolution%7Dm.%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22V%C3%A9rifiez%20que%20la%20r%C3%A9solution%20n'est%20pas%20trop%20grande%20par%20rapport%20%C3%A0%20l'%C3%A9tendue.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Dimensions%20calcul%C3%A9es%3A%20%7Bwidth%7Dx%7Bheight%7D%20pixels%22)%0A%20%20%20%20%20%20%20%20transform%20%3D%20from_bounds(%0A%20%20%20%20%20%20%20%20%20%20%20%20bounds%5B0%5D%2C%20bounds%5B1%5D%2C%20bounds%5B2%5D%2C%20bounds%5B3%5D%2C%20width%2C%20height%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20Rasterize%0A%20%20%20%20%20%20%20%20shapes_iter%20%3D%20((geom%2C%20value)%20for%20geom%2C%20value%20in%20zip(gdf.geometry%2C%20gdf%5Bcolumn%5D))%0A%20%20%20%20%20%20%20%20raster%20%3D%20rasterize(%0A%20%20%20%20%20%20%20%20%20%20%20%20shapes%3Dshapes_iter%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20out_shape%3D(height%2C%20width)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20transform%3Dtransform%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fill%3D0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20dtype%3Dnp.uint8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20all_touched%3DFalse%2C%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20Statistics%0A%20%20%20%20%20%20%20%20print(%22%5Cn%3D%3D%3D%20Statistiques%20du%20raster%20%3D%3D%3D%22)%0A%20%20%20%20%20%20%20%20print(f%22Dimensions%3A%20%7Bwidth%7Dx%7Bheight%7D%20pixels%22)%0A%20%20%20%20%20%20%20%20print(f%22R%C3%A9solution%3A%20%7Bresolution%7D%20m%2Fpixel%22)%0A%20%20%20%20%20%20%20%20print(f%22Superficie%20totale%3A%20%7B(width%20*%20height%20*%20resolution**2)%20%2F%2010000%3A.2f%7D%20ha%22)%0A%0A%20%20%20%20%20%20%20%20print(%22%5CnR%C3%A9partition%20des%20types%3A%22)%0A%20%20%20%20%20%20%20%20for%20val%20in%20sorted(np.unique(raster))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20val%20!%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20count%20%3D%20np.sum(raster%20%3D%3D%20val)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20area_m2%20%3D%20count%20*%20(resolution**2)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20area_ha%20%3D%20area_m2%20%2F%2010000%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20percentage%20%3D%20(count%20%2F%20raster.size)%20*%20100%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%20%3D%20UMEP_LABELS.get(val%2C%20%22Unknown%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20%7Bval%7D%20-%20%7Blabel%3A20s%7D%3A%20%7Barea_ha%3A8.2f%7D%20ha%20(%7Bpercentage%3A5.2f%7D%25)%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Save%20with%20metadata%0A%20%20%20%20%20%20%20%20with%20rasterio.open(%0A%20%20%20%20%20%20%20%20%20%20%20%20output_tif%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22w%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20driver%3D%22GTiff%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20height%3Dheight%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20width%3Dwidth%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20count%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20dtype%3Draster.dtype%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20crs%3Dgdf.crs%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20transform%3Dtransform%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20compress%3D%22lzw%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nodata%3D0%2C%0A%20%20%20%20%20%20%20%20)%20as%20dst%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20dst.write(raster%2C%201)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Add%20metadata%20tags%0A%20%20%20%20%20%20%20%20%20%20%20%20dst.update_tags(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20description%3D%22COSIA%20Land%20Cover%20Classification%20(UMEP%20format)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20resolution%3Df%22%7Bresolution%7Dm%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20classes%3Dstr(UMEP_LABELS)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20print(f%22%5Cn%E2%9C%85%20Fichier%20sauvegard%C3%A9%3A%20%7Boutput_tif%7D%22)%0A%0A%20%20%20%20%20%20%20%20return%20raster%0A%0A%20%20%20%20def%20vectorize_cosia_raster(cosia_tiff_path%3A%20str)%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Vectorize%20COSIA%20raster%20by%20matching%20RGB%20colors%20to%20landcover%20classes.%0A%0A%20%20%20%20%20%20%20%20Args%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20cosia_tiff_path%3A%20Path%20to%20COSIA%20GeoTIFF%20file%0A%0A%20%20%20%20%20%20%20%20Returns%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20GeoDataFrame%20with%20classified%20polygons%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20print(%22%5Cn%F0%9F%94%8D%20Vectorisation%20du%20raster%20COSIA...%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Fichier%3A%20%7Bcosia_tiff_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Create%20RGB%20to%20class%20mapping%0A%20%20%20%20%20%20%20%20rgb_to_class%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20hex_to_rgb(color)%3A%20classe%20for%20classe%2C%20color%20in%20TABLE_COLOR_COSIA.items()%0A%20%20%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%20%20with%20rasterio.open(cosia_tiff_path)%20as%20src%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Read%203%20bands%20(RGB)%0A%20%20%20%20%20%20%20%20%20%20%20%20image%20%3D%20src.read()%20%20%23%20Shape%3A%20(3%2C%20height%2C%20width)%0A%20%20%20%20%20%20%20%20%20%20%20%20transform%20%3D%20src.transform%0A%20%20%20%20%20%20%20%20%20%20%20%20crs%20%3D%20src.crs%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Combine%20RGB%20into%20single%20integer%20per%20pixel%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20R%20*%20256%5E2%20%2B%20G%20*%20256%20%2B%20B%0A%20%20%20%20%20%20%20%20%20%20%20%20combined%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(image%5B0%5D.astype(np.uint32)%20%3C%3C%2016)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(image%5B1%5D.astype(np.uint32)%20%3C%3C%208)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20image%5B2%5D.astype(np.uint32)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Vectorize%0A%20%20%20%20%20%20%20%20%20%20%20%20results%20%3D%20shapes(combined%2C%20transform%3Dtransform)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20geoms%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20rgb_values%20%3D%20%5B%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20geom%2C%20value%20in%20results%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20Decode%20integer%20to%20RGB%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value_int%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20int(value)%20if%20isinstance(value%2C%20(float%2C%20np.floating))%20else%20value%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%20%3D%20(value_int%20%3E%3E%2016)%20%26%20255%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20g%20%3D%20(value_int%20%3E%3E%208)%20%26%20255%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20b%20%3D%20value_int%20%26%20255%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20geoms.append(shape(geom))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20rgb_values.append((r%2C%20g%2C%20b))%0A%0A%20%20%20%20%20%20%20%20gdf%20%3D%20gpd.GeoDataFrame(%7B%22rgb%22%3A%20rgb_values%2C%20%22geometry%22%3A%20geoms%7D%2C%20crs%3Dcrs)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20%7Blen(gdf)%7D%20polygones%20cr%C3%A9%C3%A9s%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Match%20colors%20to%20COSIA%20classes%0A%20%20%20%20%20%20%20%20def%20match_color(rgb)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22Find%20closest%20matching%20COSIA%20class%20by%20RGB%20color.%22%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20min_dist%20%3D%20float(%22inf%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20best%20%3D%20%22Autre%22%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20target_rgb%2C%20classe%20in%20rgb_to_class.items()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dist%20%3D%20sum((a%20-%20b)%20**%202%20for%20a%2C%20b%20in%20zip(rgb%2C%20target_rgb))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20dist%20%3C%20min_dist%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20min_dist%20%3D%20dist%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20best%20%3D%20classe%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20best%0A%0A%20%20%20%20%20%20%20%20gdf%5B%22classe%22%5D%20%3D%20gdf%5B%22rgb%22%5D.apply(match_color)%0A%20%20%20%20%20%20%20%20gdf%5B%22couleur%22%5D%20%3D%20gdf%5B%22classe%22%5D.map(TABLE_COLOR_COSIA)%0A%20%20%20%20%20%20%20%20gdf%5B%22type%22%5D%20%3D%20gdf%5B%22classe%22%5D.map(COSIA_TO_UMEP)%0A%0A%20%20%20%20%20%20%20%20%23%20Drop%20RGB%20column%0A%20%20%20%20%20%20%20%20gdf%20%3D%20gdf.drop(columns%3D%5B%22rgb%22%5D)%0A%0A%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20Vectorisation%20termin%C3%A9e%3A%20%7Blen(gdf)%7D%20polygones%20classifi%C3%A9s%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Classes%20trouv%C3%A9es%3A%20%7Bgdf%5B'classe'%5D.value_counts().to_dict()%7D%22)%0A%0A%20%20%20%20%20%20%20%20return%20gdf%0A%0A%20%20%20%20def%20run_cosia(output_path%3A%20Path)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Main%20workflow%3A%20Download%20COSIA%2C%20vectorize%2C%20and%20convert%20to%20UMEP%20format.%22%22%22%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22%F0%9F%8C%8D%20COSIA%20Workflow%3A%20Download%2C%20Vectorize%2C%20and%20Convert%20to%20UMEP%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%23%20Bounding%20box%20(La%20Rochelle%20area%2C%20France)%0A%20%20%20%20%20%20%20%20%23%20Format%3A%20min_x%2C%20min_y%2C%20max_x%2C%20max_y%20(WGS84%2C%20EPSG%3A4326)%0A%20%20%20%20%20%20%20%20bbox_wgs84%20%3D%20(-1.152704%2C%2046.181627%2C%20-1.139893%2C%2046.18699)%0A%20%20%20%20%20%20%20%20working_crs%20%3D%202154%20%20%23%20Lambert%2093%0A%0A%20%20%20%20%20%20%20%20print(%22%5Cn%F0%9F%93%A6%20Configuration%3A%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Bounding%20box%3A%20%7Bbbox_wgs84%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20CRS%3A%20EPSG%3A%7Bworking_crs%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20Output%20folder%3A%20%7Boutput_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%201%3A%20Download%20COSIA%20from%20IGN%20API%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%201%3A%20Downloading%20COSIA%20from%20IGN%20API...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20cosia%20%3D%20pymdurs.geometric.Cosia(output_path%3Dstr(output_path))%0A%20%20%20%20%20%20%20%20cosia.set_bbox(*bbox_wgs84)%0A%20%20%20%20%20%20%20%20cosia.set_crs(working_crs)%0A%0A%20%20%20%20%20%20%20%20print(%22%E2%8F%B3%20Downloading%20COSIA%20from%20IGN%20API...%22)%0A%20%20%20%20%20%20%20%20cosia%20%3D%20cosia.run_ign()%0A%0A%20%20%20%20%20%20%20%20cosia_tiff_path%20%3D%20cosia.get_path_save_tiff()%0A%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20COSIA%20downloaded%3A%20%7Bcosia_tiff_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20if%20os.path.exists(cosia_tiff_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20size%20%3D%20os.path.getsize(cosia_tiff_path)%20%2F%20(1024%20*%201024)%20%20%23%20MB%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%8A%20File%20size%3A%20%7Bsize%3A.2f%7D%20MB%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%202%3A%20Vectorize%20COSIA%20raster%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%202%3A%20Vectorizing%20COSIA%20raster...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20gdf%20%3D%20vectorize_cosia_raster(cosia_tiff_path)%0A%0A%20%20%20%20%20%20%20%20%23%20Save%20vectorized%20shapefile%0A%20%20%20%20%20%20%20%20landcover_shp%20%3D%20output_path%20%2F%20%22cosia_landcover.shp%22%0A%20%20%20%20%20%20%20%20gdf.to_file(landcover_shp%2C%20driver%3D%22ESRI%20Shapefile%22)%0A%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20Shapefile%20saved%3A%20%7Blandcover_shp%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%203%3A%20Convert%20to%20UMEP%20format%20and%20rasterize%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%203%3A%20Converting%20to%20UMEP%20format%20and%20rasterizing...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%23%20Convert%20to%20working%20CRS%0A%20%20%20%20%20%20%20%20gdf%20%3D%20gdf.to_crs(working_crs)%0A%0A%20%20%20%20%20%20%20%20%23%20Filter%20valid%20geometries%0A%20%20%20%20%20%20%20%20gdf_valid%20%3D%20gdf%5Bgdf.geometry.notna()%5D.copy()%0A%20%20%20%20%20%20%20%20if%20len(gdf_valid)%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Aucune%20g%C3%A9om%C3%A9trie%20valide%2C%20arr%C3%AAt%20du%20traitement%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20return%0A%0A%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%8A%20%7Blen(gdf_valid)%7D%20g%C3%A9om%C3%A9tries%20valides%20sur%20%7Blen(gdf)%7D%20totales%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Rasterize%20to%20UMEP%20format%0A%20%20%20%20%20%20%20%20landcover_tif%20%3D%20output_path%20%2F%20%22landcover.tif%22%0A%20%20%20%20%20%20%20%20raster%20%3D%20geodataframe_to_tif_with_metadata(%0A%20%20%20%20%20%20%20%20%20%20%20%20gdf%3Dgdf_valid%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20output_tif%3Dstr(landcover_tif)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20column%3D%22type%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20resolution%3D1.0%2C%20%20%23%201%20meter%20resolution%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Summary%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22%E2%9C%85%20COSIA%20workflow%20complete!%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22%F0%9F%93%81%20Output%20files%3A%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20-%20COSIA%20raster%3A%20%7Bcosia_tiff_path%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20-%20Landcover%20shapefile%3A%20%7Blandcover_shp%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%20%20%20-%20UMEP%20landcover%20raster%3A%20%7Blandcover_tif%7D%22)%0A%0A%20%20%20%20%20%20%20%20return%20cosia%2C%20gdf%2C%20raster%0A%0A%20%20%20%20def%20fill_nodata_nan(tif_path%2C%20nodata_fallback%3D-9999)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Replace%20NaN%20with%20nodata%20in%20a%20GeoTIFF%20so%20SOLWEIG%20(int%20conversion)%20does%20not%20fail.%22%22%22%0A%20%20%20%20%20%20%20%20path%20%3D%20Path(tif_path)%0A%20%20%20%20%20%20%20%20if%20not%20path.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%0A%20%20%20%20%20%20%20%20with%20rasterio.open(path%2C%20%22r%2B%22)%20as%20src%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20nd%20%3D%20src.nodata%20if%20src.nodata%20is%20not%20None%20else%20nodata_fallback%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20i%20in%20range(1%2C%20src.count%20%2B%201)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20band%20%3D%20src.read(i)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20np.issubdtype(band.dtype%2C%20np.floating)%20and%20np.any(np.isnan(band))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20band%20%3D%20np.where(np.isnan(band)%2C%20nd%2C%20band)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20src.write(band%2C%20i)%0A%20%20%20%20%20%20%20%20%20%20%20%20src.nodata%20%3D%20nd%0A%0A%20%20%20%20def%20run_umep(output_path)%3A%0A%20%20%20%20%20%20%20%20print(%22%F0%9F%8C%86%20Starting%20UMEP%20workflow%20with%20pymdurs%20and%20umepr...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%23%20Configuration%0A%20%20%20%20%20%20%20%20output_folder%20%3D%20output_path%0A%20%20%20%20%20%20%20%20output_path%20%3D%20Path(output_folder).absolute()%0A%20%20%20%20%20%20%20%20output_path.mkdir(parents%3DTrue%2C%20exist_ok%3DTrue)%0A%20%20%20%20%20%20%20%20output_folder_str%20%3D%20str(output_path)%0A%0A%20%20%20%20%20%20%20%20%23%20Bounding%20box%20(La%20Rochelle%20area%2C%20France)%0A%20%20%20%20%20%20%20%20%23%20Format%3A%20min_x%2C%20min_y%2C%20max_x%2C%20max_y%20(WGS84%2C%20EPSG%3A4326)%0A%20%20%20%20%20%20%20%20bbox_wgs84%20%3D%20(-1.152704%2C%2046.181627%2C%20-1.139893%2C%2046.18699)%0A%0A%20%20%20%20%20%20%20%20%23%20Working%20CRS%20(Lambert%2093%20-%20EPSG%3A2154)%0A%20%20%20%20%20%20%20%20working_crs%20%3D%202154%0A%0A%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%A6%20Bounding%20box%3A%20%7Bbbox_wgs84%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%F0%9F%97%BA%EF%B8%8F%20%20Working%20CRS%3A%20EPSG%3A%7Bworking_crs%7D%22)%0A%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%81%20Output%20folder%3A%20%7Boutput_folder_str%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%201%3A%20Collect%20DEM%20from%20IGN%20API%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%201%3A%20Collecting%20DEM%20from%20IGN%20API...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20dem%20%3D%20pymdurs.geometric.Dem(output_path%3Doutput_folder_str)%0A%20%20%20%20%20%20%20%20dem.set_bbox(*bbox_wgs84)%0A%20%20%20%20%20%20%20%20dem.set_crs(working_crs)%0A%20%20%20%20%20%20%20%20dem%20%3D%20dem.run()%0A%0A%20%20%20%20%20%20%20%20dem_source%20%3D%20Path(output_folder_str)%20%2F%20%22DEM.tif%22%0A%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20DEM%20collected%20and%20saved%20to%3A%20%7Bdem_source%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%202%3A%20Load%20LiDAR%20data%20from%20IGN%20WFS%20service%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20dsm_source%20%3D%20Path(output_folder_str)%20%2F%20%22DSM.tif%22%0A%20%20%20%20%20%20%20%20if%20not%20dsm_source.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22Step%202%3A%20Loading%20LiDAR%20data%20from%20IGN%20WFS%20service...%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Create%20Lidar%20instance%0A%20%20%20%20%20%20%20%20%20%20%20%20lidar%20%3D%20pymdurs.geometric.Lidar(output_path%3Doutput_folder_str)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Set%20bounding%20box%20(same%20as%20DEM)%0A%20%20%20%20%20%20%20%20%20%20%20%20lidar.set_bbox(*bbox_wgs84)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Set%20CRS%20(same%20as%20DEM)%0A%20%20%20%20%20%20%20%20%20%20%20%20lidar.set_crs(working_crs)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%F0%9F%93%A6%20Bounding%20box%20set%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20geo%20%3D%20lidar.geo_core%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%F0%9F%97%BA%EF%B8%8F%20%20CRS%3A%20%7Bgeo.epsg%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Generate%20CDSM%20from%20vegetation%20and%20water%20classes%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Classification%3A%202%20%3D%20Ground%2C%203%20%3D%20Low%20Vegetation%2C%204%20%3D%20Medium%20Vegetation%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%205%20%3D%20High%20Vegetation%2C%209%20%3D%20Water%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%F0%9F%8C%B3%20Generating%20CDSM%20from%20vegetation%20and%20water%20classes...%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20classification_list%20%3D%20%5B3%2C%204%2C%205%2C%209%5D%20%20%23%20Vegetation%20and%20water%20classes%0A%20%20%20%20%20%20%20%20%20%20%20%20lidar.run(file_name%3D%22CDSM.tif%22%2C%20classification_list%3Dclassification_list)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9C%85%20CDSM%20generated%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Generate%20DSM%20from%20ground%20and%20buildings%20classes%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%F0%9F%8F%A2%20Generating%20DSM%20from%20ground%20and%20buildings%20classes...%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20classification_list%20%3D%20%5B2%2C%206%5D%20%20%23%20Ground%20and%20buildings%20classes%0A%20%20%20%20%20%20%20%20%20%20%20%20dsm_output_path%20%3D%20lidar.run(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20file_name%3D%22DSM.tif%22%2C%20classification_list%3Dclassification_list%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9C%85%20LiDAR%20processing%20complete!%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%81%20DSM%20GeoTIFF%20saved%20to%3A%20%7Bdsm_output_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Check%20if%20file%20exists%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20os.path.exists(dsm_output_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20size%20%3D%20os.path.getsize(dsm_output_path)%20%2F%20(1024%20*%201024)%20%20%23%20MB%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%8A%20DSM%20GeoTIFF%20file%20size%3A%20%7Bsize%3A.2f%7D%20MB%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22%F0%9F%93%8A%20File%20contains%203%20bands%3A%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22%20%20%20-%20Band%201%3A%20DSM%20(Digital%20Surface%20Model)%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22%20%20%20-%20Band%202%3A%20DTM%20(Digital%20Terrain%20Model)%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22%20%20%20-%20Band%203%3A%20CHM%20(Canopy%20Height%20Model)%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%203%3A%20Warp%20and%20clip%20rasters%20using%20mask%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%203%3A%20Warping%20and%20clipping%20rasters%20with%20mask...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20mask_shp_path%20%3D%20Path(output_folder_str)%20%2F%20%22mask.shp%22%0A%20%20%20%20%20%20%20%20if%20mask_shp_path.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20warp_options%20%3D%20gdal.WarpOptions(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20multithread%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20format%3D%22GTiff%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20xRes%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20yRes%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20outputType%3Dgdalconst.GDT_Float32%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dstNodata%3DNone%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dstSRS%3D%22EPSG%3A2154%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cropToCutline%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cutlineDSName%3Dstr(mask_shp_path)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cutlineLayer%3D%22mask%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Clip%20DEM%0A%20%20%20%20%20%20%20%20%20%20%20%20dem_clip_path%20%3D%20Path(output_folder_str)%20%2F%20%22DEM_clip.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20dem_source.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gdal.Warp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20destNameOrDestDS%3Dstr(dem_clip_path)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20srcDSOrSrcDSTab%3Dstr(dem_source)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20options%3Dwarp_options%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20DEM%20clipped%20to%3A%20%7Bdem_clip_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Clip%20DSM%0A%20%20%20%20%20%20%20%20%20%20%20%20dsm_clip_path%20%3D%20Path(output_folder_str)%20%2F%20%22DSM_clip.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20dsm_source.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gdal.Warp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20destNameOrDestDS%3Dstr(dsm_clip_path)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20srcDSOrSrcDSTab%3Dstr(dsm_source)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20options%3Dwarp_options%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20DSM%20clipped%20to%3A%20%7Bdsm_clip_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Clip%20CDSM%0A%20%20%20%20%20%20%20%20%20%20%20%20cdsm_clip_path%20%3D%20Path(output_folder_str)%20%2F%20%22CDSM_clip.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20cdsm_source%20%3D%20Path(output_folder_str)%20%2F%20%22CDSM.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20cdsm_source.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gdal.Warp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20destNameOrDestDS%3Dstr(cdsm_clip_path)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20srcDSOrSrcDSTab%3Dstr(cdsm_source)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20options%3Dwarp_options%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20CDSM%20clipped%20to%3A%20%7Bcdsm_clip_path%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Clip%20Landcover%0A%20%20%20%20%20%20%20%20%20%20%20%20landcover_clip_path%20%3D%20Path(output_folder_str)%20%2F%20%22landcover_clip.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20landcover_source%20%3D%20Path(output_folder_str)%20%2F%20%22landcover.tif%22%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20landcover_source.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gdal.Warp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20destNameOrDestDS%3Dstr(landcover_clip_path)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20srcDSOrSrcDSTab%3Dstr(landcover_source)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20options%3Dwarp_options%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20Landcover%20clipped%20to%3A%20%7Blandcover_clip_path%7D%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Mask%20shapefile%20not%20found%2C%20skipping%20clipping%22)%0A%0A%20%20%20%20%20%20%20%20%23%20Set%20paths%20for%20later%20steps%0A%20%20%20%20%20%20%20%20dsm_path%20%3D%20Path(output_folder_str)%20%2F%20%22DSM_clip.tif%22%0A%20%20%20%20%20%20%20%20cdsm_path%20%3D%20Path(output_folder_str)%20%2F%20%22CDSM_clip.tif%22%0A%20%20%20%20%20%20%20%20dem_tiff_path%20%3D%20Path(output_folder_str)%20%2F%20%22DEM_clip.tif%22%0A%20%20%20%20%20%20%20%20lc_path%20%3D%20Path(output_folder_str)%20%2F%20%22landcover_clip.tif%22%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%204%3A%20Calculate%20Sky%20View%20Factor%20(SVF)%20using%20umepr%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22Step%204%3A%20Calculating%20Sky%20View%20Factor%20(SVF)%20using%20umepr...%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20if%20not%20HAS_UMEPR%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Skipping%20SVF%20calculation%20-%20umepr%20not%20available%22)%0A%20%20%20%20%20%20%20%20elif%20dsm_path%20and%20os.path.exists(dsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Get%20bounding%20box%20from%20the%20actual%20raster%20bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20This%20ensures%20the%20bbox%20matches%20exactly%20with%20the%20raster%20extent%0A%20%20%20%20%20%20%20%20%20%20%20%20with%20rasterio.open(dsm_path)%20as%20dsm_src%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20Get%20the%20actual%20bounds%20of%20the%20raster%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bounds%20%3D%20dsm_src.bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20total_extents%20%3D%20%5Bbounds.left%2C%20bounds.bottom%2C%20bounds.right%2C%20bounds.top%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%90%20Using%20raster%20bounds%20for%20SVF%3A%20%7Btotal_extents%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20svf_output_dir%20%3D%20str(output_path%20%2F%20%22svf%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20os.makedirs(svf_output_dir%2C%20exist_ok%3DTrue)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20svf.generate_svf(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dsm_path%3Ddsm_path%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bbox%3Dtotal_extents%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20out_dir%3Dsvf_output_dir%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cdsm_path%3Dcdsm_path%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20cdsm_path%20and%20os.path.exists(cdsm_path)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20trunk_ratio_perc%3D25%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20trans_veg_perc%3D3%2C%20%20%23%203%25%20transmissivity%20for%20vegetation%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20use_tiled_loading%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20SVF%20calculation%20complete!%20Output%20in%3A%20%7Bsvf_output_dir%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20except%20Exception%20as%20e%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9A%A0%EF%B8%8F%20%20SVF%20calculation%20failed%3A%20%7Be%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20%20Error%20details%3A%20%7Btype(e).__name__%7D%3A%20%7Be%7D%22)%0A%20%20%20%20%20%20%20%20elif%20dsm_path%20and%20os.path.exists(dsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Skipping%20SVF%20calculation%20-%20umepr%20not%20available%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Skipping%20SVF%20calculation%20-%20DSM%20not%20available%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%205%3A%20Generate%20wall%20heights%20for%20SOLWEIG%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20if%20HAS_UMEP%20and%20dsm_path%20and%20os.path.exists(dsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22Step%205%3A%20Generating%20wall%20heights%20for%20SOLWEIG...%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Get%20bounding%20box%20from%20the%20actual%20raster%20bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20with%20rasterio.open(dsm_path)%20as%20dsm_src%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bounds%20%3D%20dsm_src.bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20total_extents%20%3D%20%5Bbounds.left%2C%20bounds.bottom%2C%20bounds.right%2C%20bounds.top%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20walls_output_dir%20%3D%20str(output_path%20%2F%20%22walls%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20os.makedirs(walls_output_dir%2C%20exist_ok%3DTrue)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20wall_heightaspect_algorithm.generate_wall_hts(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dsm_path%3Ddsm_path%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bbox%3Dtotal_extents%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20out_dir%3Dwalls_output_dir%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9C%85%20Wall%20heights%20generated!%20Output%20in%3A%20%7Bwalls_output_dir%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20except%20Exception%20as%20e%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%E2%9A%A0%EF%B8%8F%20%20Wall%20height%20generation%20failed%3A%20%7Be%7D%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Step%206%3A%20Run%20SOLWEIG%20for%20thermal%20comfort%20analysis%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20if%20HAS_UMEPR%20and%20dsm_path%20and%20os.path.exists(dsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22Step%206%3A%20Running%20SOLWEIG%20for%20thermal%20comfort%20analysis...%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Replace%20NaN%20with%20nodata%20in%20rasters%20to%20avoid%20%22cannot%20convert%20float%20NaN%20to%20integer%22%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20rast%20in%20(dsm_path%2C%20cdsm_path%2C%20dem_tiff_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20rast%20and%20os.path.exists(rast)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20fill_nodata_nan(rast)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20SRR%20%3D%20solweig_runner_rust.SolweigRunRust(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22configsolweig.ini%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22parametersforsolweig.json%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20use_tiled_loading%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20tile_size%3D1024%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20SRR.run()%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9C%85%20SOLWEIG%20run%20complete!%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22%E2%9A%A0%EF%B8%8F%20%20Skipping%20SOLWEIG%20-%20missing%20requirements%20or%20DSM%20not%20available%22)%0A%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20%23%20Summary%0A%20%20%20%20%20%20%20%20%23%20%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%3D%0A%20%20%20%20%20%20%20%20print(%22%5Cn%22%20%2B%20%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(%22%E2%9C%85%20UMEP%20workflow%20complete!%22)%0A%20%20%20%20%20%20%20%20print(%22%3D%22%20*%2060)%0A%20%20%20%20%20%20%20%20print(f%22%F0%9F%93%81%20All%20outputs%20saved%20to%3A%20%7Boutput_folder_str%7D%22)%0A%20%20%20%20%20%20%20%20print(%22%5CnGenerated%20files%3A%22)%0A%20%20%20%20%20%20%20%20if%20dsm_path%20and%20os.path.exists(dsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20-%20DSM%3A%20%7Bdsm_path%7D%22)%0A%20%20%20%20%20%20%20%20if%20cdsm_path%20and%20os.path.exists(cdsm_path)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20-%20CDSM%3A%20%7Bcdsm_path%7D%22)%0A%20%20%20%20%20%20%20%20if%20dem_tiff_path.exists()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20-%20DEM%3A%20%7Bdem_tiff_path%7D%22)%0A%20%20%20%20%20%20%20%20if%20os.path.exists(str(output_path%20%2F%20%22svf%22))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20-%20SVF%3A%20%7Boutput_path%20%2F%20'svf'%7D%22)%0A%20%20%20%20%20%20%20%20if%20HAS_UMEP%20and%20os.path.exists(str(output_path%20%2F%20%22walls%22))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20-%20Wall%20heights%3A%20%7Boutput_path%20%2F%20'walls'%7D%22)%0A%0A%20%20%20%20def%20_main_()%3A%0A%20%20%20%20%20%20%20%20output_folder%20%3D%20%22.%2Foutput%2Fumep_workflow%22%0A%20%20%20%20%20%20%20%20output_path%20%3D%20Path(output_folder).absolute()%0A%20%20%20%20%20%20%20%20output_path.mkdir(parents%3DTrue%2C%20exist_ok%3DTrue)%0A%20%20%20%20%20%20%20%20cosia%2C%20gdf%2C%20raster%20%3D%20run_cosia(output_path%3Doutput_path)%0A%20%20%20%20%20%20%20%20run_umep(output_path%3Doutput_folder)%0A%0A%20%20%20%20_main_()%0A%20%20%20%20return%20(rasterio%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Visualisation%20Tmrt%20(Mean%20Radiant%20Temperature)%0A%0A%20%20%20%20Affichage%20du%20raster%20Tmrt_average.tif%20g%C3%A9n%C3%A9r%C3%A9%20par%20SOLWEIG.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(rasterio)%3A%0A%20%20%20%20from%20matplotlib%20import%20pyplot%20as%20plt%0A%20%20%20%20from%20rasterio.plot%20import%20show%0A%0A%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(15%2C%2015))%0A%20%20%20%20src%20%3D%20rasterio.open(%22output%2FTmrt_average.tif%22)%0A%20%20%20%20plt.imshow(src.read(1)%2C%20cmap%3D%22pink%22)%0A%20%20%20%20show(src%2C%20ax%3Dax%2C%20cmap%3D%22plasma%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A