search

LittoSim / LittoSim_model

Participatory simulation model of coastal flooding management
http://littosim.github.io/LittoSim_model/
7 stars 3 forks source link

génération du fichier de rugosité à partir du corine landcover #62

Open gnoubi opened 3 years ago

gnoubi commented 3 years ago

tout est dans le titre mais ca permettrait de finaliser l'automatisation

laatabix commented 3 years ago

Ça se fait déjà dans le script de preprocessing : le script R utilise le fichier "clc_manning.csv" pour remplacer les valeurs du fichier land_cover.shp avec les bons coefficients de manning :

        ```
        # the rugosity grid with values of land_cover
        rugo <- rasterize(land_cover, ras, field='ID');
        ids <- as.vector(as.numeric(as.character(land_cover$ID)))
        vals <- as.vector(as.numeric(as.character(land_cover$cover_type)))
        for(i in(1:length(ids))){
          rugo[rugo == ids[i]] <- vals[i]
        }
        # replacing land cover values with manning coefficients
        clc_manning_file <- xmlValue(xmlElementsByTagName(el = rugosity, 'manning')[[1]]);
        mannings <- read.csv(clc_manning_file, header=TRUE, sep=";")
        for(xman in(1:nrow(mannings))){
            rugo[rugo == as.integer(mannings[xman,]$clc)] <- as.double(mannings[xman,]$manning);
        }
        writeRaster(rugo, paste(output_folder,"rugosity",sep="/"), format = "ascii", overwrite=TRUE);
        print(paste('files rugosity.asc and start.asc have been created in',output_folder));