gma 1.1.6 | 2.0.0a1 (2023.05.14) 更新日志

 从现在开始，gma 开始同步更新 gma 2 的前期预览版，并与 gma 1 的正式版同步发布。

gma 1.1.6 更新日志

新增

 1、添加 3组 新的指北针样式。

修复

 1、【map】高分辨栅格数据集内存溢出问题。
 　添加高分辨率栅格数据集且进行强制数据缩放时 偶现内存溢出错误，此问题是由于掩膜过程设计的不太合理导致。

优化

 1、TVDI。TVDI 所有计算结果的存储类型更换为 具名元组（以前为元组），方便抽取对应数据。
 2、坐标参考系统。不再依附于 map 模块，新增 gma.crs 模块，用于统一 gma 库整体的坐标参考系统。

待开放

 1、已含有超过10种其他的遥感植被指数。

简单示例

from gma.map import plot, inres

MapF = plot.MapFrame(BaseMapProj = 3035, Extent = None)

# 1.添加图层（可以重复添加多个图层）
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = None, EdgeColor = 'gray', LineWidth = 0.1)

# 2.添加经纬网
Grid = MapF.AddGridLines()

# 3.绘制四种指北针
Compass = MapF.AddCompass(LOC = (0.1, 0.9), Color = 'black', Style = 'GMACompass1')
Compass = MapF.AddCompass(LOC = (0.2, 0.9), Color = 'blue',  Style = 'GMACompass2')
Compass = MapF.AddCompass(LOC = (0.1, 0.8), Color = 'green', Style = 'GMACompass3')
Compass = MapF.AddCompass(LOC = (0.2, 0.8), Color = 'red',  Style = 'GMACompass4')

# 4.设置地图框（包括刻度和经纬度的标注）
Frame = MapF.SetFrame()

gma 2.0.0a1 更新日志

整体架构调整

重构

 1、底层算法逻辑彻底重构，类继承关系更加明确。
 2、功能函数重新归类，部分函数进行更名或移除。

重要函数更新

 1、Open 栅格文件。
 　现在，打开的栅格文件可以直接进行重投影、重采样、另存为栅格文件或转为矢量图层等操作。
 2、Open 矢量文件。
 　现在，打开的矢量图层可以直接进行裁剪、擦除、重投影，另存为矢量文件或转为栅格数据集等操作。
 3、更多投影支持。
 　现在，gma.crs.ProjMethod 添加了超过 30 种投影方法以供使用。

其他更新

 1、控制台更清爽。此版本开始 gma 会直接抑制 GDAL/OGR 的内部警告。
 2、稳定版 gma 更新。gma 1.1.6
的新功能也会包含在此版本。

简单示例

from gma import crs

############### 构建自定义投影的主要思路
# 1.定义一个 Bonne 投影方法，中央子午线为东经 112°
ProjMethod = crs.ProjMethod.Bonne(CentralMeridian = 112)

# 2.创建一个地理坐标系
## 2.1 创建一个自定义椭球体。椭球体名称：'New Ellips'，不使用默认椭球体；自定义椭球体长半轴 6378137.0m，反扁率 298.257223563
Ellipsoid = crs.Ellipsoid(Name = 'New Ellips', UseDefault = False, SemiMajor = 6378137.0, InvFlattening = 298.257223563)
## 2.2 生成一个地理坐标系。地理坐标系名称：'GMA GCS'，角度单位为 'Degree'，中央经线为 'Greenwich'，椭球体为 2.1 创建的自定义椭球体
GCS = crs.GeogCS(Name = 'GMA GCS', AngularUnit = 'Degree', PRIMEM = 'Greenwich', Ellipsoid = Ellipsoid)

# 3.创建投影坐标系。投影坐标系名称：'GMA PCS'，线性单位为 'Meter'，投影方法为 1 定义的投影方法，地理坐标系为 2.2 创建的地理坐标系
Proj = crs.ProjCS(Name = 'GMA PCS', LinearUnit = 'Meter', ProjMethod = ProjMethod, GCS = GCS)
print(Proj.Export())

PROJCS[“GMA PCS”,GEOGCS[“GMA GCS”,DATUM[“New Ellips”,SPHEROID[“New Ellips”,6378137,298.257223563]],PRIMEM[“Greenwich”,0],UNIT[“Degree”,0.0174532925199433]],PROJECTION[“Bonne”],PARAMETER[“central_meridian”,112],PARAMETER[“standard_parallel_1”,60],PARAMETER[“false_easting”,0],PARAMETER[“false_northing”,0],UNIT[“Meter”,1],AXIS[“Easting”,EAST],AXIS[“Northing”,NORTH]]

from gma.map import plot, inres

MapF = plot.MapFrame(BaseMapProj = Proj, Extent = None)

# 1.添加图层（可以重复添加多个图层）
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = None, EdgeColor = 'gray', LineWidth = 0.1)

# 2.添加经纬网
Grid = MapF.AddGridLines()

# 3.设置地图框（包括刻度和经纬度的标注）
Frame = MapF.SetFrame(ShowLeft = False, ShowBottom = False, ShowRight = False, ShowTop = False)