update land use based allocation

geocaruso · May 30, 2023 · a55b5ce · a55b5ce
1 parent ea24fa1
commit a55b5ce
Show file tree

Hide file tree

Showing 6 changed files with 22 additions and 0 deletions.
diff --git a/.DS_Store b/.DS_Store
diff --git a/R/map.lux.grid.pop.R b/R/map.lux.grid.pop.R
@@ -65,12 +65,16 @@ merged2many$Pop2021EU_sh_of_LAU2<-merged_agg$Pop2021EU/merged_agg$SumPop2021EU
 
 Grid1km_LAU2_Pop2021EU<-merged2many[match(Grid1km_LAU2_popEU$CELLCODE,merged2many$CELLCODE),]
 
+
+#Mapping correspondance shares
 pfi_sh<-ggplot.themap(Grid1km_LAU2_Pop2021EU,
                      "Pop2021EU_sh_of_LAU2",n=7, style="fisher")
 pfi_sh
 
+
 #save
 sf::st_write(Grid1km_LAU2_Pop2021EU,"data/Grid1km_LAU2_Pop2021EU.gpkg", delete_dsn=TRUE)
 write.csv(sf::st_drop_geometry(Grid1km_LAU2_Pop2021EU),"data/Grid1km_LAU2_Pop2021EU.csv")
 
+#print
 
diff --git a/R/map.lux.grid.urban.R b/R/map.lux.grid.urban.R
@@ -0,0 +1,17 @@
+#map.lux.grid.urban.R
+
+#Add urban surfaces to grid and compute share of each cell in each
+# previously assigned communes (usign total surface) based on urban surface
+#Urban surface is built from Urban Atlas soil sealing percentage categories
+
+#urban atlas 2012 downloaded from https://land.copernicus.eu/local/urban-atlas/urban-atlas-2012?tab=download
+#urban atlas 2018 downloaded from https://land.copernicus.eu/local/urban-atlas/urban-atlas-2018?tab=download
+#And unzipped in data folder
+
+UA12<-sf::st_read("data/EXT/LU001L1_LUXEMBOURG_UA2012_revised_v021.gpkg")
+UA18<-sf::st_read("data/EXT/LU001L1_LUXEMBOURG_UA2018_v013.gpkg")
+artifcode <-c("11100", "11210", "11220", "11230", "11240")
+artif12<- UA12[UA12$code_2012 %in% artifcode,]
+artif18<-UA18[UA18$code_2018 %in% artifcode,]
+sf::st_crs(artif12)<- sf::st_crs(grid_estat_lu)
+sf::st_crs(artif18)<- sf::st_crs(grid_estat_lu)
diff --git a/data/.DS_Store b/data/.DS_Store
diff --git a/data/.gitignore b/data/.gitignore
@@ -0,0 +1 @@
+EXT/
diff --git a/data/EXT/.DS_Store b/data/EXT/.DS_Store