gma 2 教程(二)数据操作:4. 普通栅格操作

普通栅格定义

  如下图所示,gma将普通栅格定义为1个(记录数据的)数组+1组(表征栅格属性信息的)元数据组成的数据集(DataSet)。
gma 2 教程(二)数据操作:4. 普通栅格操作_第1张图片

打开栅格文件

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

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

获取栅属性

  打开的栅格数据集(DataSet)包含的属性及示例结果如下表:
gma 2 教程(二)数据操作:4. 普通栅格操作_第2张图片

重投影

  (一)方法简介
gma 2 教程(二)数据操作:4. 普通栅格操作_第3张图片
  (二)示例

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

gma 2 教程(二)数据操作:4. 普通栅格操作_第4张图片

重采样

  (一)方法简介
gma 2 教程(二)数据操作:4. 普通栅格操作_第5张图片
  (二)示例

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


gma 2 教程(二)数据操作:4. 普通栅格操作_第6张图片

读取数据到NumPy数组

  (一)方法简介
gma 2 教程(二)数据操作:4. 普通栅格操作_第7张图片
  (二)示例

Data = DataSet.ToArray()
print(Data)

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

转为矢量图层

  (一)方法简介
gma 2 教程(二)数据操作:4. 普通栅格操作_第8张图片
  (二)示例(矢量图层过大,建议用其他数据尝试)

Layer = DataSet.ToLayer()
print(Layer)

另存为栅格文件

  (一)方法简介
gma 2 教程(二)数据操作:4. 普通栅格操作_第9张图片
  (二)示例(另存为ENVI文件)

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

你可能感兴趣的:(gma,教程)