2012-05-26 1 views
4

이것은 후속 조치 this question입니다. 나는 그 질문에 주어진 답을 사용하여 데이터의 변이를 그릴 수 있었다. 거리 축에 enter image description heregstat 배리 그램에서 거리 단위는 무엇입니까?

당신이 볼 수 있듯이, 번호 : imhttps://gist.github.com/2780792에 주어진 dput 데이터를 포함하는 변수가 내가 사용하고있는 코드는 (이 나에게 다음과 같은 플롯을 제공

library(gstat) 
point_data <- as(im, 'SpatialPointsDataFrame') 
v <- variogram(band1 ~ 1, data = point_data) 
plot(v) 

입니다 0에서 약 150 그러나 데이터의 단위는 도의 범위 :

> head(coordinates(point_data)) 
      x  y 
[1,] -1.849353 52.06165 
[2,] -1.759728 52.06165 
[3,] -1.401227 52.06165 
[4,] -1.311602 52.06165 
[5,] -1.221977 52.06165 
[6,] -1.132352 52.06165 

,상의 거리 측정의 단위는 무엇 점을 감안 X 축? 나는 그들이 그들이 코디네이트와 같은 부대에있을 것을 기대할 것이다. 그들이 미터가되는 것은 이해가되지 않습니다. 나는 그들이 km 일 수 있다고 생각합니다.

gstat 매뉴얼에서 이것에 관해서 아무 것도 찾을 수 없기 때문에 기본 데이터 단위를 사용한다고 가정 했었지만 이것이 사실이 아닌 것 같습니다.

아이디어가 있으십니까?

+0

나는 잘 모르지만 km 단위의 큰 원형 거리 일 수 있습니다 (처음 두 관측치 사이의 거리는 예를 들어 ~ 10km입니다). –

+0

거리를 수동으로 확인하려면이 블로그를 통해 해당 계산을 실행하는 기능을 설정하는 데 도움이 될 수 있습니다. http://www.r-bloggers.com/great-circle-distance-calculations-in-r/ 추신 : 여기 있습니다 GIS 관련 질문에 대한 http://gis.stackexchange.com입니다. SO보다 훨씬 작은 커뮤니티이지만, 거기에 R 사용자가 분명히 있습니다! –

답변

9

메일 링리스트에서 유사한 질문을하는 사람들이 있었고 도움말 파일의 세부 사항을 다시 확인했습니다. 데이터가 비 프로젝트입니다 (proj = + longlat). the documentation으로 메일 목록에 나와있는 정보는 다음과 같습니다 : km 단위로 great circle distances을 기반으로해야합니다. 베리오 그램 투영 (논리) 매개 변수를 체크 아웃 :

예상 데이터의 경우, 유클리드 거리, 투영 대원 거리 (km)에 대해 계산된다.

관련 문제