2017-11-13 21 views
1

R 및 QGIS에서 래스터 파일을 가져 오는 방법에 대한 기본적인 질문이 무엇인지 궁금합니다.R과 QGIS 사이의 래스터 값 차이

단일 대역 래스터가 있습니다. I 래스터 패키지의 "래스터"기능을 이용하여, R로 가져올 때, 화소 값의 범위를 얻을 : I 래스터 벽돌이 층을 적층하는 경우

class  : RasterLayer 
dimensions : 10980, 10980, 120560400 (nrow, ncol, ncell) 
resolution : 10, 10 (x, y) 
extent  : 6e+05, 709800, 5590200, 5700000 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : /data/MTDA/CGS_S2_RADIOMETRY/2017/10/15/S2B_20171015T104525Z_31UFS_TOC_V100/S2B_20171015T104525Z_31UFS_TOC-B02_10M_V100.tif 
names  : S2B_20171015T104525Z_31UFS_TOC.B02_10M_V100 
values  : -32768, 32767 (min, max) 

, I이 최소 - 최대 값을 얻을 (느린 실제 : I가 QGIS 동일한 래스터를로드하는 경우

class  : RasterLayer 
band  : 2 (of 11 bands) 
dimensions : 10980, 10980, 120560400 (nrow, ncol, ncell) 
resolution : 10, 10 (x, y) 
extent  : 6e+05, 709800, 5590200, 5700000 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : /tmp/Rtmp882dZS/raster/r_tmp_2017-11-10_172819_11532_86514.grd 

names  : S2B_20171015T104525Z_31UFS_TOC.B02_10M_V100 
values  : -1129, 9994 (min, max) 

그러나 최소값은 228이고 최대 값은 907 (I는 옵션 "넓이 풀"로이 값을 산출하고 "정확도) ".

그렇다면이 차이는 어디에서 비롯된 것입니까? 정확히 무엇이 R과 QGIS는 다음 작업을 수행 중입니다 ...

답변

-1

결국 나는 차이점을 발견했습니다! R에 질문하면 실제 최소값/최대 값을 얻습니다. 대신 QGIS는 누적 개수로 최소/최대 값을 계산합니다. "최소/최대 값로드"(래스터 속성 창에서)를 "최소/최대"로 설정하면 R 값과 동일한 값을 얻었습니다.

1

첫 번째 개체의 경우 파일이 제공하지 않거나 올바르게 입력되지 않아서 최소 및 최대 값을 알 수 없습니다. RasterLayerr하면

r <- setMinMax(r) 

은 그들이 무엇인지 확인하려면 할 수 있습니다. 두 번째 레이어와 동일하지 않으면 무언가를 섞어 보았을 것입니다. 이러한 것들을 명확히 한 후에 QGIS와 비교하는 것이 유용 할 수 있습니다. 이를 위해서는 예제 파일을 제공해야 할 것입니다.