2017-04-11 10 views
0

나는 잠시 동안 가지고 있었던 문제로 나를 도울 수 있기를 바랍니다. MaxEnt에 바이어스 파일을 만들어야하는데,이 튜토리얼은 https://scottrinnan.wordpress.com/2015/08/31/how-to-construct-a-bias-file-with-r-for-use-in-maxent-modeling/이며이 파일을 내 자신의 상황으로 바꿨다. 그러나, 지금 붙어있어 ...R의 커널 밀도 함수가 x와 y의 해상도가 같지 않음

kde2d 함수를 사용하여 2 차원 커널 밀도 추정을 작성한 다음이를 래스터로 변환해야합니다. 그러나 생성 된 래스터는 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 

biasraster의 원래 해상도는 모두 x와 y에 대한 0.00833333이지만, dens.ras의 해상도는 0.0104052, 0.00833333 (X, Y) (변경되었습니다 :

내가 무슨 짓을 y 해상도가 올바른 것입니다.)

질문에서 볼 수 있듯이, 나는 코딩에 관해서는 완전히 멍청하다. 나는 지금 일주일 동안해야할 일을 알아 내려고 노력했지만 아무런 대답도 찾을 수 없으므로 여기 누군가가 나를 도울 수 있기를 바랍니다.

+0

예제 파일이 없으면 인덱싱이 방향 중 하나의 범위를 줄여 다른 해상도로 이어질 것이라고 추측 할 수 있습니다. 아마도'MASS :: kde2d'에서'h ='인수를 지정하려고 시도 할 것입니다. – nya

답변

0

나는 동일한 문제가있었습니다. 난 단지 사용하여 해결할 수 :

<이

dens.ras - 재 샘플 (dens.ras는, climdat는 방법 = "선형") #make 바이어스 climdat

의 동일한 해상도의 파일 만이 additionaly, 나는 또한 필요 dens.ras를 자르십시오.

densmod < 수확물 (dens.ras, 정도 (climdat))이 추가 단계로

내 바이어스 파일은 잘 작동.

+0

네, 이것이 완벽하게 작동했습니다! 정말 고맙습니다! – Nadja