2013-01-21 6 views
3

자동 매핑을 시도하고 있습니다. 나는 불확실성이 isolines에 의해 표시되는 불확실성 맵을 원합니다. 나는이 isolines과 그들과 연관된 값을 export 한 다음 웹 맵 서비스를 통해 표시하려고합니다. 가능한가? 코드를 따라 당신이 당신의 질문에 대한 직접적인 대답없는 등고선WMS를 통해 윤곽선을 내보내려면

require(gstat) 
require(lattice) 
data(meuse) 
data(meuse.grid) 
coordinates(meuse) = ~x+y 
gridded(meuse.grid) = ~x+y 
m <- vgm(.59, "Sph", 874, .04) 
# ordinary kriging: 
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m) 
l=as.data.frame(x) 
contourplot(var1.var ~ x+y,data=l,cuts=10) 
+2

윤곽선을 폴리곤으로 변환해야하는 이유는 무엇입니까? 왜 당신의 웹지도 서비스를 통해 라인으로 그들을 제공합니까? – Spacedman

+0

@ 스페이스 맨. 나는 제목을 편집했다. –

+0

@Spacedman 안돼요? –

답변

1

를 생성 할 수 있습니다,하지만 어쩌면 그것은 올바른 방향을 가리킬 수 있습니다 : 당신이 당신의 contourplot

에서 개체를 생성 할 수 있습니다

CP<- contourplot(var1.var ~ x+y,data=l,cuts=10) 

그리고 것은 그 구조

fix(CP) 

지방 검사 보는 시도 등각에 대한 타 여기에 저장됩니다

x <- CP$panel.args.common[1]$x 
y <- CP$panel.args.common[2]$y 
z <- CP$panel.args.common[3]$z 

그러나 나는 그것을 추출하는 방법을 지금 알아낼 수 없습니다. 아마도 약간의 도움이 될 수 있습니다.

+1

친애하는 Andrey, 답변 주셔서 감사합니다.하지만이 값은 등고선 값이 아닙니다. 이것들은 객체 x를 출력합니다. –