Je débute dans la cartographie avec R et j'aimerais solliciter votre aide afin de compléter mon script.
Ce script a pour but de créer une carte type couranto (le package ncdf4 va lire dans mes fichiers netcdf).
Je suis dans le système wgs et ce que j'aimerais faire pour compléter mon script est de pouvoir ajouter des points en fonction de leurs coordonnées (latitude, longitude), dans ce même système et les visualiser sur ma carte de courants.
J'aimerais également pouvoir "zoomer" sur une fenêtre moins large, par exemple plus proche du trait de côte.
Voici le lien vers mon fichier en input pour reproduire le code > https://drive.google.com/file/d/15zx3ZQ ... sp=sharing
Ci-dessous, mon code qui produit un plot couranto, j'aimerais votre aide/vos conseils pour pouvoir ajouter des points en fonction de leurs coordonnées (importer un .csv éventuellement?) et pouvoir zoomer sur certaines parties de cette carte.
Les packages sont nombreux et pas tous utiles pour ce que je souhaite faire là, je ne les ai pas supprimés car ce code me sert également au traitement d'autres fichiers netcdf.
Code : Tout sélectionner
#repertoire
setwd("C:/Users/cebar/Desktop/sorties_CMS")
dir()
#packages
library(raster)
library(sp)
library(rgdal)
library(dplyr)
library(ggplot2)
library(devtools)
library(rJava)
library(ncdf4)
library(reshape2)
library(spatstat)
library(maptools)
library(RColorBrewer)
library(chron)
library(lattice)
library(rasterVis)
library(RColorBrewer)
library(ggthemes)
library(viridis)
#systeme de projection
wgs<-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
#lecture du fichier couranto et extraction variables
ncin<-nc_open("nest_1_20170511030000.nc")
lat <- ncvar_get(ncin,"Latitude")
nlat <- dim(lat)
lon <- ncvar_get(ncin,"Longitude")
nlon <- dim(lon)
time <- ncvar_get(ncin,"Time")
tunits <- ncatt_get(ncin,"Time","units")
nt <- dim(time)
zv_array <- ncvar_get(ncin,"zv")
dlname <- ncatt_get(ncin,"zv","long_name")
dunits <- ncatt_get(ncin,"zv","units")
fillvalue <- ncatt_get(ncin,"zv","_FillValue")
dim(zv_array)
zu_array <- ncvar_get(ncin,"zu")
dlnameU <- ncatt_get(ncin,"zu","long_name")
dunitsU <- ncatt_get(ncin,"zu","units")
fillvalue <- ncatt_get(ncin,"zu","_FillValue")
dim(zu_array)
speed_array<-sqrt((zv_array^2)+(zu_array^2))#calcule de la vitesse = sqrt(uv^2+zv^2)
dim(speed_array)
speed_array[speed_array==fillvalue$value] <- NA
length(na.omit(as.vector(speed_array[,,1])))
zv_slice <- speed_array[,,1]
zv_slice<-raster(zv_slice)
uvel<-flip(t(zv_slice), direction='y')
plot(uvel)
#add extent information
extent(uvel)<-c(min(lon), max(lon), min(lat), max(lat))
crs(uvel)<-wgs