本文主要是介绍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')
这篇关于gma 2 教程(二)数据操作:4. 普通栅格操作的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!