dasymetric: Create Dasymetric Maps

## code to prepare `corine_18` dataset goes here

load("data/districts.rda")
temp = tempfile()
temp2 = tempfile()
#download.file("https://daten.gdz.bkg.bund.de/produkte/dlm/clc5_2018/aktuell/clc5_2018.utm32s.shape.zip",temp)
temp = "C:/Users/janni/Downloads/clc5_2018.utm32s.shape.zip"
require(sf)
unzip(zipfile=temp, exdir = temp2)
file_list = list.files(temp2, pattern = "xx.shp$", recursive = TRUE, full.names=TRUE)
shp_list = lapply(file_list,read_sf)
shp_list[[1]] = st_crop(shp_list[[1]],st_bbox(districts))
shp_list[[2]] = st_crop(shp_list[[2]],st_bbox(districts))
shp_list[[3]] = st_crop(shp_list[[3]],st_bbox(districts))
shp_list[[4]] = st_crop(shp_list[[4]],st_bbox(districts))
shp_list[[5]] = st_crop(shp_list[[5]],st_bbox(districts))
corine_18 = rbind(shp_list[[1]], shp_list[[2]], shp_list[[3]], shp_list[[4]], shp_list[[5]])
# corine_18 = st_crop(SHP_file, st_bbox(districts))

sf::st_crs(corine_18)$wkt = "PROJCRS[\"ETRS89 / UTM zone 32N\",\n    BASEGEOGCRS[\"ETRS89\",\n        DATUM[\"European Terrestrial Reference System 1989\",\n            ELLIPSOID[\"GRS 1980\",6378137,298.257222101,\n                LENGTHUNIT[\"metre\",1]]],\n        PRIMEM[\"Greenwich\",0,\n            ANGLEUNIT[\"degree\",0.0174532925199433]],\n        ID[\"EPSG\",4258]],\n    CONVERSION[\"UTM zone 32N\",\n        METHOD[\"Transverse Mercator\",\n            ID[\"EPSG\",9807]],\n        PARAMETER[\"Latitude of natural origin\",0,\n            ANGLEUNIT[\"degree\",0.0174532925199433],\n            ID[\"EPSG\",8801]],\n        PARAMETER[\"Longitude of natural origin\",9,\n            ANGLEUNIT[\"degree\",0.0174532925199433],\n            ID[\"EPSG\",8802]],\n        PARAMETER[\"Scale factor at natural origin\",0.9996,\n            SCALEUNIT[\"unity\",1],\n            ID[\"EPSG\",8805]],\n        PARAMETER[\"False easting\",500000,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8806]],\n        PARAMETER[\"False northing\",0,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8807]]],\n    CS[Cartesian,2],\n        AXIS[\"(E)\",east,\n            ORDER[1],\n            LENGTHUNIT[\"metre\",1]],\n        AXIS[\"(N)\",north,\n            ORDER[2],\n            LENGTHUNIT[\"metre\",1]],\n    USAGE[\n        SCOPE[\"Engineering survey, topographic mapping.\"],\n        AREA[\"Europe between 6°E and 12°E: Austria; Belgium; Denmark - onshore and offshore; Germany - onshore and offshore; Norway including - onshore and offshore; Spain - offshore.\"],\n        BBOX[38.76,6,83.92,12]],\n    ID[\"EPSG\",25832]]"

usethis::use_data(corine_18, overwrite = TRUE)