gma 2 教程（二）数据操作：4. 普通栅格操作

普通栅格定义

 　如下图所示，gma将普通栅格定义为1个（记录数据的）数组+1组（表征栅格属性信息的）元数据组成的数据集（DataSet）。

打开栅格文件

 　gma打开各种普通栅格文件（GTiff、ENVI、PNG、JPEG等），返回一个栅格数据集（DataSet）类。

from gma import io
DataSet = io.Open("HYP_50M_SR_W.tif")
print(DataSet)

<gma.algos.dataio.rads.DataSet object at 0x…>

获取栅属性

 　打开的栅格数据集（DataSet）包含的属性及示例结果如下表：

重投影

 　（一）方法简介

 　（二）示例

DataSet6931 = DataSet.Reproject('EPSG:6931')
print(DataSet6931)

<gma.algos.dataio.rads.DataSet at 0x…>

重采样

 　（一）方法简介

 　（二）示例

DataSet5 = DataSet.Resample(5)
print(DataSet5)

<gma.algos.dataio.rads.DataSet at 0x…>

读取数据到NumPy数组

 　（一）方法简介

 　（二）示例

Data = DataSet.ToArray()
print(Data)

[[[132 128 128 … 132 132 132]
[119 118 120 … 118 120 118]
[120 121 123 … 119 121 121]
…] …]

转为矢量图层

 　（一）方法简介

 　（二）示例（矢量图层过大，建议用其他数据尝试）

Layer = DataSet.ToLayer()
print(Layer)

<gma.algos.dataio.veds.Layer at 0x…>

另存为栅格文件

 　（一）方法简介

 　（二）示例（另存为ENVI文件）

DataSet.SaveAs('HYP_50M_SR_W.dat', Format='ENVI')