--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/crea_carte_G.R	Wed Sep 11 09:18:37 2024 +0000
@@ -0,0 +1,38 @@
+#Author : Seguineau Pauline
+
+
+#Create a map from cluster
+
+library(sf)
+library(tmap)
+library(dplyr)
+
+args = commandArgs(trailingOnly=TRUE) 
+if (length(args)==0)
+{
+    stop("This tool needs at least one argument")
+}else{
+    data <- args[1]
+}
+
+clus <- read.table(data, sep="\t" ,header=TRUE, na.strings = "na")
+
+#tmap method
+
+sf_data <- st_as_sf(clus, coords = c("long", "lat"), crs =4326)
+
+grouped_data <- sf_data %>%
+  group_by(cluster) %>%
+  summarize()
+
+map <- tm_shape(grouped_data) + 
+  tm_dots(col = "cluster", palette = "Accent", n=256, size = 0.1, title = "ecoregions")+
+  tm_scale_bar(position = c("right","top"))+
+  tm_compass(position = c("right","top"))+
+  tm_layout(frame = FALSE,legend.position = c("left","bottom"))+
+  tm_xlab("Longitude")+
+  tm_ylab("Latitude")+
+  tm_grid(alpha = 0.2)
+
+#Save the map 
+tmap_save(map, "ecoregions.png")