저는 정기적으로 Python에서 GDAL을 사용하여 자체 GeoTIFF 래스터를 만듭니다.
from osgeo import gdal
from numpy import random
data = random.uniform(0, 10, (300, 200))
driver = gdal.GetDriverByName('GTiff')
ds = driver.Create('MyRaster.tif', 200, 300)
band = ds.GetRasterBand(1)
band.WriteArray(data)
ds = band = None # save, close
그러나 ArcCatalog / ArcGIS로 결과를 볼 때는 통계가 없으므로 검정 또는 회색으로 나타납니다. 래스터를 마우스 오른쪽 버튼으로 클릭하고 ArcCatalog에서 "Calculate Statistics ..."를 선택하거나 (다른 방법으로 수행 할 수있는 방법이 있음) 명령 프롬프트에서 gdalinfo 를 사용 하여이 문제 를 해결 합니다.
gdalinfo -stats MyRaster.tif
MyRaster.tif.aux.xml
ArcGIS에서 래스터를 적절하게 스케일링하기 위해 사용하는를 생성 합니다. PAM (Persistent Auxiliary Metadata) 파일에는 통계, 특히 최소값과 최대 값이 포함됩니다.
<PAMDataset>
<PAMRasterBand band="1">
<Metadata>
<MDI key="STATISTICS_MINIMUM">0</MDI>
<MDI key="STATISTICS_MAXIMUM">10</MDI>
<MDI key="STATISTICS_MEAN">5.0189833333333</MDI>
<MDI key="STATISTICS_STDDEV">2.9131294111984</MDI>
</Metadata>
</PAMRasterBand>
</PAMDataset>
내 질문 : GDAL이 통계 파일을 만들 수있는 기본 제공 방법이 gdalinfo -stats
있습니까 (명령 사용 이외 )? 아니면 내가 직접 써야합니까?
band.GetStatistics(0,1)
실제로 통계를 계산하고 하나의 파일에 GeoTIFF 메타 데이터에 추가합니다. 다른 파일이 필요하지 않습니다. 그러나 Esri 제품으로 테스트하면 ArcGIS 9.3 이전 버전이 아닌 ArcGIS 10.0 이상에서만 작동합니다.