R中的内核密度函数给出不相等的x和y分辨率

Question

希望您能为我解决已有一段时间的问题。 我需要为MaxEnt制作一个偏差文件，为此我使用了本教程： https ://scottrinnan.wordpress.com/2015/08/31/how-to-construct-a-bias-file-with-r- 在maxent模型中使用/并将其更改为我自己的情况。 但是，我现在被困住了...

我需要使用kde2d函数创建2d内核密度估计，然后将其转换为栅格。 但是，创建的栅格对于x和y具有不同的分辨率。 这是一个问题，因为我必须在不会接受不相等的x和y分辨率的MaxEnt中使用它。

这是我所做的：

biasraster <- raster("file.tif") #load raster with all the occurrences  
presences <-which(values(biasraster)==1)  
pres.locs<- coordinates(biasraster)[presences,]  
dens <-kde2d(pres.locs[,1],pres.locs[,2],n=c(nrow(biasraster),ncol(biasraster))) #2d kernel density function on the biasraster  
dens.ras<-raster(dens) #create raster from kde2d function 

对于x和y，biasaster的原始分辨率均为0.00833333，但dens.ras的分辨率已更改为0.0104052，即0.00833333（x，y）（因此y分辨率是正确的）。

从这个问题可能可以看出，当涉及到编码时（在r中），我是个菜鸟。 大约一周以来，我一直在努力寻找解决方法，但找不到任何可行的答案，因此希望这里有人可以帮助我。

Answer 1

我有同样的问题。 我只能使用以下方法解决它：

dens.ras <-resample（dens.ras，climdat，method =“ bilinear”）＃以相同的分辨率设置climdat的偏差文件

但另外，我还需要剪切dens.ras。

densmod <-crop（dens.ras，程度（climdat））

通过这些额外的步骤，我的偏差文件效果很好。