Je voudrais utiliser la fonction image pour représenter une zone spécifique...
Voici un code reproductible:
Code : Tout sélectionner
#Génération de ma grille de points
r <- 1
scale = 0.025
x <- seq(-r, r, scale)
nrow.m <- ncol.m <- length(x)
Xg <- expand.grid(X1 = 1:nrow.m, X2=1:ncol.m)
Xg$X1g <- (Xg$X1-1)*scale-r
Xg$X2g <- (Xg$X2-1)*scale-r
Xg$sphere <- sqrt(Xg[,3]^2 + Xg[,4]^2)
Xg <- Xg[Xg$sphere<r,] #grille de points
Xg$D <- as.numeric(runif(nrow(Xg))<0.1) # ajout d'une colonne avec des 0 et des 1
m <- matrix(nrow=nrow.m,ncol=ncol.m)
m[as.matrix(Xg[1:2])] <- Xg$D
image(x,x,m,asp=1,col=c("gray","white")) # représentation selon 0 et 1
Au lieu d'ajouter une condition comme pour la colonne D, j'ai appliqué plusieurs traitements sur cette matrice Xg qui ont permis d'éliminer un certain nombre de points (Xg').
Je voudrais utiliser la fonction image pour représenter la matrice Xg en gris et en blanc les points restants (Xg').
Quelqu'un sait-il comment faire?
Je vous remercie par avance,
Diane