2013-11-14 4 views
1

나는R에서 비 구조화 메쉬의 데이터를 시각화하는 방법은 무엇입니까?

id lat lon value 

어쩌면 형태의 unstructered 그리드 적용되는 숫자 해양 모델에서 ASCII 데이터를 가지고 s.th. 그와 같은 I를 시각화 할 지금

dsize <- 100 
data=data.frame(id=seq(1:dsize), 
      lat=sort(runif(dsize)*10+2570740), 
      lon=runif(dsize)*10+6926070, 
      val=sin(2*pi*(1:dsize)/dsize) + runif(1,-1,1)*10) 

는 (X = 위도, 경도 = Y, Z = 브로)

  1. 3D 예컨대 함께 초래 RGL 라이브러리
  2. 2D, 등고선 플롯
  3. 위도 나 경도를 따라
  4. 섹션
내가 x와 y의 값을 증가 달성하지 않기 때문에 내가 갇혀

... 나는 그것이 경우에 도움이 될 것입니다 생각

+0

검색 –

답변

2

예를 들어 val 변수를 보간하면 akima 패키지의 interp() 함수. 그게 네가 쫓아 온거야? `akima` 관련된 답변 SO

enter image description here

require(rgl) 
require(akima) 

dsize <- 100 
data=data.frame(id=seq(1:dsize), 
      lat=sort(runif(dsize)*10+2570740), 
      lon=runif(dsize)*10+6926070, 
      val=sin(2*pi*(1:dsize/dsize + runif(1,-1,1)*10))) 

#use akima to interpolate a *val* surface 
surf<-interp(data$lat,data$lon,data$val) 

persp3d(x=surf$x, y=surf$y, z=surf$z, 
    theta=-50, phi=-25, expand=0.75, 
    alpha=0.3, col="green", ticktype="detailed", 
    xlab="lat", ylab="long", zlab="level",axes=TRUE) 

rgl.viewpoint(scale=c(1,1,2)) 
관련 문제