2012-09-06 2 views
0

나는이 코드의 결과가 모양을netcdf 파일에서 이진 파일을 작성하는 방법은 무엇입니까?

library(ncdf) 
download.file("http://gswp/Fixed/SoilDepth.nc", destfile="SoilDepth.nc") 
soil <- open.ncdf("SoilDepth.nc") 
soil$var[[3]] -> var3 
get.var.ncdf(soil, var3) -> SoilDepth 
download.file("http://gswp/Fixed/landmask_gswp.nc", destfile="landmask.nc") 
landmask <- open.ncdf("landmask.nc") 
landmask$var[[3]] -> varland 
get.var.ncdf(landmask, varland) -> land 
land = t(land) 
land[land==1] <- SoilDepth 
land[land==0] <- NA 
land = t(land) 
image(land) 

netCDF의 파일의 바이너리 파일을 작성하고 싶습니다! [토양의지도] [1]

지금은에 그것을 쓰고 싶어 이진 파일 :

결과가 거꾸로 된 이미지입니다.

+0

어디에서 쓸 수 있습니까? 또한 모든 전이와 함께 무엇입니까? 불필요한 코드를 없애기 위해 각 단계별로 신중하게 작업하면 도움이됩니다. [행/열 순서를 바꾸려면 – mdsumner

+1

을 사용하십시오. mdsumner는 아마도'y <- y [nrow (y) : 1,]'와 같은 의미 였기 때문에 행렬 y의 행을 뒤집을 수 있습니다. – MvG

답변

1

ylim=c(1,0)을 지정하여 직접 역 오리지나 날을 요청했습니다. 간단히

image(y) 

에 마지막 라인을 변경하고 이중 전치 자원을 먹고 제외하고 아무것도하지 않는 당신은 괜찮을거야.

상세 내용 : 일본어와 다시 판독 된 데이터의 유일한 차이점은 전자가 (불가) NA 있던 장소에서 후자 NaN (즉 숫자가 아님)을 갖는다는 사실이다. 본 변형 수율 완전히 동일한 데이터 취소이 후

y[is.nan(y)] <- NA 

, y하면 land에서 구별된다.

관련 문제