R: Création d'une covariable sur un système de coordonnées projetées

Question

J'essaye de créer une covariable en utilisant le paquet Spatstat, j'ai un fichier shp en tant que polygone afin d'assigner une valeur à la covariable si les points sont en dehors du polygone. Mon problème est qu'il y a un problème sur la projection des coordonnées parce que quand j'ai fini le processus, j'ai un chiffre étrange, quelqu'un peut-il m'aider s'il vous plaît. Mon code:

#Creating covariate Z: 

W<- owin(xrange=c(767768.3,768773.6),yrange=c(517335.8,518044.8)) 
xp<-c(768773.6) #cambia! 
yp<-c(518044.8) 
X<-ppp(x=xp,y=yp,W) 
par(mfrow=c(2,2)) 
Z<-density(X,500)/(7e-06) 

#reading and projecting the shp: 

Prodes1<-readOGR(dsn="forestsseparate.shp", layer="forestsseparate", dropNULLGeometries=TRUE) 
projection(Prodes1)<-CRS("+proj=utm +zone=18 +ellps=intl +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs") 

#Extracting the shp polygon coordinates: 

xcoords<-rev(Prodes1@polygons[[1]]@Polygons[[1]]@coords[,1]) 
ycoords<-rev(Prodes1@polygons[[1]]@Polygons[[1]]@coords[,2]) 
pruebalu<-owin(poly=list(x=xcoords,y=ycoords)) 

mx<-(Z$xrange[2]-Z$xrange[1])/Z$xstep 
my<-(Z$yrange[2]-Z$yrange[1])/Z$ystep 

xcoorde<-matrix(NA,nrow=1,ncol=mx) 
ycoorde<-matrix(NA,nrow=1,ncol=my) 

#Finding each pixel center 

for(i in 1:mx){ 
    xcoorde[i]<-(Z$xrange[1]+i*Z$xstep) -Z$xstep/2 
} 

for(i in 1:my){ 
    ycoorde[i]<-(Z$yrange[1]+i*Z$ystep) -Z$ystep/2 
} 

#Changing the pixel value if is outside the polygon: 

for(i in 1:length(Z$xcol)){ 
    for(j in 1: length(Z$yrow)){ 

     if(inside.owin(xcoorde[i],ycoorde[j], pruebalu)==FALSE){ 
      whichPixel<-nearest.raster.point(xcoorde[i], ycoorde[j], Z) 
      Z$v[whichPixel$col,whichPixel$row]<-0.4 
     } 
    } 
} 

Voici mes résultats:

Answer 1

Désolé, je fait une erreur avec l'ordre des éléments Z de $ v, voici la correction:

for(i in 1:length(Z$xcol)){ 
    for(j in 1: length(Z$yrow)){ 

     if(inside.owin(xcoorde[i],ycoorde[j], pruebalu)==FALSE){ 
      whichPixel<-nearest.raster.point(xcoorde[i], ycoorde[j], Z) 
      Z$v[whichPixel$row,whichPixel$col]<-0.4 
     } 
    } 
} 

Maintenant c'est parfait!