数据切片是 Python 中非常实用的方法,Numpy、Pandas 等第三方库的切片操作为数据处理提供了不少便利。如果能对栅格/矢量数据进行切片,那会使地理数据处理也变得方便和快捷。
基于此,自 gma 2.0.6.10 开始,gma 针对 打开的栅格数据 设计了切片操作方法。
gma 网站:https://gma.luosgeo.com
PyPi 项目:https://pypi.org/project/gma/
示例数据:
- 链接:https://pan.baidu.com/s/1zr8p0MZUDD3nWYxe2A9jzA?pwd=gj2w
- 提取码:gj2w
栅格数据切片
from gma import io
DataSet = io.ReadRaster("NE2_50M_SR_W.tif")
一些基本方法
另存栅格数据集
DataSet.SaveAs('TestDS.tif', Format='GTiff')
栅格信息
预览
切片方法说明
DataSet[波段序号/波段序号切片, 行切片, 列切片] --> DataSet
第一个参数:整数型列表(波段序号,从 0 开始)或切片器(slice),操作波段
第二个参数:切片器(slice),操作行(Y)
第三个参数:切片器(slice),操作列(X)
可以依次配置三个参数,被忽略的参数默认为全部