태그 보관물: raster

raster

* .adf 파일을 R로 읽는 방법? = WGS84 + towgs84 =

.adf 파일을 R로로드하고 싶습니다. 데이터는 다음 페이지에 있습니다.
http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057

인터넷에서 조사 한 결과 다음 코드를 시도했습니다. 문제는 RasterLayer 클래스에서 거기에 있으면 안되는 부정적인 가치를 얻는다는 것입니다. 왜 이런 일이 발생하는지 모르겠으므로 누군가가 나를 도울 수 있기를 바랍니다!?

암호:

library(rgdal)
library(RColorBrewer)
dpath<- path...

x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)

‘r’의 출력은 다음과 같습니다.

r 클래스 : RasterLayer 크기 : 2160, 4320, 9331200 (nrow, ncol, ncell) 해상도 : 0.08333333, 0.08333333 (x, y) 범위 : -180, 180, -90, 90 (xmin, xmax, ymin, ymax) 좌표 심판. : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0 + no_defs 데이터 소스 : 메모리 이름 : band1 값 : -997, 16 (최소, 최대)

값에서 ’16’은 16 개 클래스의 성장 기간 길이를 나타냅니다. 그러나 나는 그 ‘-997’이 어디에서 왔는지 궁금합니다. 아마도 코드에 문제가있을 수 있습니다. 심판?

다음은 ‘xx’의 데이터 요약입니다.

데이터 요약 : Min. 1 분기 중간 평균 3 쿼터 최대 NA의 -997 3 5-9 8 16 7123158

그리고 우리가 xx의 데이터를 더 자세히 보면 :

테이블 (xx $ band1)

-99712 34 5678 9 10 11 12 31711 429643 83011 166674207228270161 240958 183342 118608 98795 88473 73743 56022 13 14 15 16 30104 45521 52216 31832

실제로이 ‘-997’이 있습니다. 나는 NA가 오존 인 것 같아서 데이터 로딩에 문제가 있거나 데이터를 이해하지 못합니까?



답변

당신은 거의 옳습니다 :

해양의 경우 NODATA는 -32768로 설정됩니다. 또한 -997은 해안선에서 제외되지 않은 큰 호수에 설정됩니다.

호수에서는 픽셀 내용 (성장 기간)이 의미가 없으므로 -997을 NODATA로 안전하게 처리 할 수 ​​있습니다.


답변