Je débute avec R en essayant d'adapter certains de mes scripts Matlab, et je suis en train de me casser les dents sur une représentation graphique.
Dans le but de représenter des paramètres environnementaux en 3D (latitude, longitude et profondeur), je suis en train d'essayer d'ajouter un shapefile de trait de côte pour une meilleure visualisation.
Pour l'instant, j'utilise plot 3D pour le mapping (en codage 0 pour le continent et NA pour l'océan), mais mon problème est que cette approche me fait apparaître des lignes bizarres qui gâche le visuel. Je ne sais pas comment les enlever, ou si une autre fonction serait plus adaptée pour ce que je souhaite faire. J'ajoute ici la figure en l'état et les lignes de code qui correspondent.
Merci d'avance des conseils et aides!
Code : Tout sélectionner
rm(list=ls())
library(maptools)
library(plot3D)
library(rgdal)
library(rgl)
library(raster)
SPDF <- readOGR(dsn = ".", layer = "ne_10m_land")
## Example RasterLayer
r <- raster(nrow=1e4, ncol=1e4, crs=proj4string(SPDF))
r[] <- 1:length(r)
## crop and mask
r2 <- crop(r, extent(SPDF))
r3 <- mask(r2, SPDF)
## values to 1
r3@data@values[!is.na(r3@data@values)] <- 0
r3@data@values
## Set adjust=FALSE to fix it
decorate3d()
myPal <- colorRampPalette(brewer.pal(1, 'Greys'))
plot3D(r3, col=myPal,adjust=TRUE, zfac=1,rev=TRUE)