Mercurial > repos > ecology > map_shp
annotate Map_shp.R @ 0:ce9979488ebb draft default tip
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
author | ecology |
---|---|
date | Sat, 25 Nov 2023 15:18:33 +0000 |
parents | |
children |
rev | line source |
---|---|
0
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
1 library(ggplot2) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
2 library(sf) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
3 library(dplyr) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
4 library(RColorBrewer) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
5 library(ggspatial) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
6 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
7 args = commandArgs(trailingOnly=TRUE) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
8 if (length(args)==0) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
9 { |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
10 stop("This tool needs at least one argument") |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
11 }else{ |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
12 dataMap <- args[1] |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
13 dataEvo <- args[2] |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
14 title <- args[3] |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
15 legend <- args[4] |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
16 coord <- args[5] |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
17 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
18 } |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
19 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
20 title <- gsub("\\\\n", "\n", title) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
21 legend <-gsub("\\\\n", "\n", legend) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
22 #read data |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
23 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
24 data_map = st_read(dataMap) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
25 data_evo = read.delim(dataEvo,header=TRUE,sep="\t") |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
26 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
27 #bring together data |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
28 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
29 data_fin = bind_cols(data_map,data_evo[2]) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
30 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
31 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
32 # define the data intervals |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
33 intervals <- cut(data_fin$Evolution_rate, breaks = c(-Inf, 0, 9, 20, Inf), labels = c("Moins de 0", "0 à 10", "10 à 20", "Plus de 20")) |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
34 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
35 # Make the map with ggplot2 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
36 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
37 if (coord == "true"){ |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
38 ggplot(data_fin) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
39 geom_sf(aes(fill = intervals)) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
40 scale_fill_manual(values = c('#D9F0D3',"#A6DBA0","#5AAE61","#1B7837")) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
41 labs(title = title, fill = legend) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
42 annotation_scale() |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
43 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
44 #outuput |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
45 ggsave("map.pdf", device = "pdf") |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
46 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
47 }else{ |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
48 ggplot(data_fin) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
49 geom_sf(aes(fill = intervals)) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
50 scale_fill_manual(values = c('#D9F0D3',"#A6DBA0","#5AAE61","#1B7837")) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
51 labs(title = title, fill = legend) + |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
52 theme_void()+ |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
53 annotation_scale() |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
54 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
55 #outuput |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
56 ggsave("map.pdf", device = "pdf") |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
57 } |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
58 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
59 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
60 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
61 |
ce9979488ebb
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Geom_mean_workflow commit 3f11e193fd9ba5bf0c706cd5d65d6398166776cb
ecology
parents:
diff
changeset
|
62 |