基于 gma 的 GIS 空间参考系统探索

最新推荐文章于 2024-07-27 13:31:57 发布
最新推荐文章于 2024-07-27 13:31:57 发布
阅读量344 收藏 1
点赞数 1
文章标签： java 数据库 javascript
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42155937/article/details/129904375
版权

0 了解空间参考

为了能够正确地描述要素的位置和形状，需要引入一个用于定义位置的框架，也就是空间参考。空间参考用于存储各类要素坐标属性的坐标系系统，借此描述一个要素的真实位置。主要包括两大部分：坐标系统 和 坐标域。

坐标系统
坐标系统是一个二维或者三维的参照系，用于定位坐标点，通过坐标系统可以确定要素在地球上的真实位置。常用的坐标系主要是以下两种：大地坐标系和投影坐标系。

坐标域
坐标域是一个要素中坐标的允许取值范围。

1 了解 gma （1.1.5）中的空间参考功能

1.1 gma 空间参考功能在库中的位置

gma 中空间参考功能位于 gma.algorithm.core.spat，用以初始化各种坐标系/投影，以方便检索和使用参数和查询相关信息，这是基础的空间参考；基于此构建的扩展（增强）空间参考功能整合在地图工具 gma.map.rcs 中 SpatRef 类中 。gma 库（1.1.5）及其空间参考功能分布如下：

1.2 SpatRef 类支持初始化的坐标系格式

格式示例 （以 WGS 84 为例）
EPSG4326 (或 'EPSG:4326' )
MICoordSys'Earth Projection 1, 104'
PCI['LONG/LAT D000', 'DEGREE', (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)]
PROJJSON'{\n "$schema": "https://proj.org/schemas/v0.4/projjson.schema.json",\n "type": "GeographicCRS",\n "name": "WGS 84",\n "datum": {\n "type": "GeodeticReferenceFrame",\n "name": "World Geodetic System 1984",\n "ellipsoid": {\n "name": "WGS 84",\n "semi_major_axis": 6378137,\n "inverse_flattening": 298.257223563\n }\n },\n "coordinate_system": {\n "subtype": "ellipsoidal",\n "axis": [\n {\n "name": "Geodetic latitude",\n "abbreviation": "Lat",\n "direction": "north",\n "unit": "degree"\n },\n {\n "name": "Geodetic longitude",\n "abbreviation": "Lon",\n "direction": "east",\n "unit": "degree"\n }\n ]\n },\n "id": {\n "authority": "EPSG",\n "code": 4326\n }\n}'
PrettyWkt'GEOGCS["WGS 84",\n DATUM["WGS_1984",\n SPHEROID["WGS 84",6378137,298.257223563,\n AUTHORITY["EPSG","7030"]],\n AUTHORITY["EPSG","6326"]],\n PRIMEM["Greenwich",0,\n AUTHORITY["EPSG","8901"]],\n UNIT["degree",0.0174532925199433,\n AUTHORITY["EPSG","9122"]],\n AXIS["Latitude",NORTH],\n AXIS["Longitude",EAST],\n AUTHORITY["EPSG","4326"]]'
Proj4'+proj=longlat +datum=WGS84 +no_defs'
USGS[0, 0, (0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), 12]
Wkt'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]'
XML'<gml:GeographicCRS gml:id="ogrcrs1">\n <gml:srsName>WGS 84</gml:srsName>\n <gml:srsID>\n <gml:name codeSpace="urn:ogc:def:crs:EPSG::">4326</gml:name>\n </gml:srsID>\n <gml:usesEllipsoidalCS>\n <gml:EllipsoidalCS gml:id="ogrcrs2">\n <gml:csName>ellipsoidal</gml:csName>\n <gml:csID>\n <gml:name codeSpace="urn:ogc:def:cs:EPSG::">6402</gml:name>\n </gml:csID>\n <gml:usesAxis>\n <gml:CoordinateSystemAxis gml:id="ogrcrs3" gml:uom="urn:ogc:def:uom:EPSG::9102">\n <gml:name>Geodetic latitude</gml:name>\n <gml:axisID>\n <gml:name codeSpace="urn:ogc:def:axis:EPSG::">9901</gml:name>\n </gml:axisID>\n <gml:axisAbbrev>Lat</gml:axisAbbrev>\n <gml:axisDirection>north</gml:axisDirection>\n </gml:CoordinateSystemAxis>\n </gml:usesAxis>\n <gml:usesAxis>\n <gml:CoordinateSystemAxis gml:id="ogrcrs4" gml:uom="urn:ogc:def:uom:EPSG::9102">\n <gml:name>Geodetic longitude</gml:name>\n <gml:axisID>\n <gml:name codeSpace="urn:ogc:def:axis:EPSG::">9902</gml:name>\n </gml:axisID>\n <gml:axisAbbrev>Lon</gml:axisAbbrev>\n <gml:axisDirection>east</gml:axisDirection>\n </gml:CoordinateSystemAxis>\n </gml:usesAxis>\n </gml:EllipsoidalCS>\n </gml:usesEllipsoidalCS>\n <gml:usesGeodeticDatum>\n <gml:GeodeticDatum gml:id="ogrcrs5">\n <gml:datumName>WGS_1984</gml:datumName>\n <gml:datumID>\n <gml:name codeSpace="urn:ogc:def:datum:EPSG::">6326</gml:name>\n </gml:datumID>\n <gml:usesPrimeMeridian>\n <gml:PrimeMeridian gml:id="ogrcrs6">\n <gml:meridianName>Greenwich</gml:meridianName>\n <gml:meridianID>\n <gml:name codeSpace="urn:ogc:def:meridian:EPSG::">8901</gml:name>\n </gml:meridianID>\n <gml:greenwichLongitude>\n <gml:angle uom="urn:ogc:def:uom:EPSG::9102">0</gml:angle>\n </gml:greenwichLongitude>\n </gml:PrimeMeridian>\n </gml:usesPrimeMeridian>\n <gml:usesEllipsoid>\n <gml:Ellipsoid gml:id="ogrcrs7">\n <gml:ellipsoidName>WGS 84</gml:ellipsoidName>\n <gml:ellipsoidID>\n <gml:name codeSpace="urn:ogc:def:ellipsoid:EPSG::">7030</gml:name>\n </gml:ellipsoidID>\n <gml:semiMajorAxis uom="urn:ogc:def:uom:EPSG::9001">6378137</gml:semiMajorAxis>\n <gml:secondDefiningParameter>\n <gml:inverseFlattening uom="urn:ogc:def:uom:EPSG::9201">298.257223563</gml:inverseFlattening>\n </gml:secondDefiningParameter>\n </gml:Ellipsoid>\n </gml:usesEllipsoid>\n </gml:GeodeticDatum>\n </gml:usesGeodeticDatum>\n</gml:GeographicCRS>\n'
部分主要坐标系名称WGS84

2 探索空间参考（坐标系/投影）系统

from gma.map import rcs
SR = rcs.SpatRef('WGS84')
类型名称描述使用结果
属性AeraOfUse坐标系使用的经纬度范围。无法识别范围的坐标系则根据中央经线自动生成最大范围 （gma 对全球范围的坐标系进行了修改，其最东端与最西端永远相差0.02° ）SR.AeraOfUse[-179.99, -89.99, 179.99, 89.99]
属性AngularUnits角度单位值SR.AngularUnits0.017453292519943295
属性AngularUnitsName角度单位名称SR.AngularUnitsName‘degree’
属性EPSG坐标系 EPSG 编码。没有 EPSG 编码的坐标系则值为 0SR.EPSG4326
方法Export导出坐标系。支持 MICoordSys，PCI，PROJJSON，PrettyWkt，Proj4，USGS，Wkt，XML 格式SR.Export(Format = 'Proj4')‘+proj=longlat +datum=WGS84 +no_defs’
方法GetProjParms获取坐标系的投影参数SR.GetProjParms(){}
方法GetUseFeature获取坐标系的使用区域 FeatureSR.GetUseFeature()<gma.algorithm.core.dataio.Feature at 0x…>
方法GetUseFeatureInWGS84获取坐标系的使用区域 Feature（WGS84）SR.GetUseFeatureInWGS84()<gma.algorithm.core.dataio.Feature at 0x…>
属性InvFlattening椭球体反扁率SR.InvFlattening298.257223563
属性LinearUnits线形单位值SR.LinearUnits1.0
属性LinearUnitsName线形单位名称SR.LinearUnitsName‘Meter’
属性Name坐标系名称SR.Name‘WGS 84’
属性PRIMEM坐标系的本初子午线SR.PRIMEM{‘Greenwich’: 0.0}
属性ProjWKT 格式的坐标系SR.Proj ‘GEOGCS[“WGS 84”,DATUM[“WGS_1984”,SPHEROID[“WGS 84”,6378137,298.257223563,AUTHORITY[“EPSG”,“7030”]],AUTHORITY[“EPSG”,“6326”]],PRIMEM[“Greenwich”,0,AUTHORITY[“EPSG”,“8901”]],UNIT[“degree”,0.0174532925199433,AUTHORITY[“EPSG”,“9122”]],AXIS[“Latitude”,NORTH],AXIS[“Longitude”,EAST],AUTHORITY[“EPSG”,“4326”]]’
属性SemiMajor椭球体长半轴SR.SemiMajor6378137.0
属性SemiMinor椭球体短半轴SR.SemiMinor6356752.314245179
属性Type坐标系类型SR.Type‘Geographic’

3 在地图上查看坐标系的适用范围

结合 gma 绘图工具和 GetUseFeatureInWGS84 方法，可以很容易的查看某个坐标系使用区域在地图上的位置。本文以 :

  • EPSG:4326（ WGS84）
  • EPSG:4490（China Geodetic Coordinate System 2000）
  • EPSG:2163（NAD27 / US National Atlas Equal Area）
  • EPSG:32650（WGS 84 / UTM zone 50N）

四个坐标系为例进行展示。

from gma.map import plot, rcs
WorldDS = plot.GetWorldDataSource()
WorldLayer = WorldDS.GetLayer(0)

3.1 全球性坐标系（WGS84）

MapF = plot.MapFrame()
MapL1 = MapF.AddLayer(WorldLayer, FaceColor = 'none', EdgeColor = 'black', LineWidth = 0.1, Zorder = 1)

# 绘制坐标系范围
Spat = rcs.SpatRef('WGS84')
UseFeature = Spat.GetUseFeatureInWGS84()
MapF.AddFeature(UseFeature)

3.2 地方坐标系（EPSG:4490 / China Geodetic Coordinate System 2000）

MapF = plot.MapFrame()
MapL1 = MapF.AddLayer(WorldLayer, FaceColor = 'none', EdgeColor = 'black', LineWidth = 0.1, Zorder = 1)

# 绘制坐标系范围
Spat = rcs.SpatRef(4490)
UseFeature = Spat.GetUseFeatureInWGS84()
MapF.AddFeature(UseFeature)

3.3 地方坐标系（EPSG:2163 / NAD27 / US National Atlas Equal Area）

MapF = plot.MapFrame()
MapL1 = MapF.AddLayer(WorldLayer, FaceColor = 'none', EdgeColor = 'black', LineWidth = 0.1, Zorder = 1)

# 绘制坐标系范围
Spat = rcs.SpatRef(2163)
UseFeature = Spat.GetUseFeatureInWGS84()
MapF.AddFeature(UseFeature)

3.4 地方坐标系（EPSG:32650 / WGS 84 / UTM zone 50N）

MapF = plot.MapFrame()
MapL1 = MapF.AddLayer(WorldLayer, FaceColor = 'none', EdgeColor = 'black', LineWidth = 0.1, Zorder = 1)

# 绘制坐标系范围
Spat = rcs.SpatRef(32650)
UseFeature = Spat.GetUseFeatureInWGS84()
MapF.AddFeature(UseFeature)

3.5 换一个视角（底图：EPSG:6931 / WGS 84 / NSIDC EASE-Grid 2.0 North）

MapF = plot.MapFrame(BaseMapProj = 6931)
MapL1 = MapF.AddLayer(WorldLayer, FaceColor = 'none', EdgeColor = 'black', LineWidth = 0.1, Zorder = 1)

# 绘制坐标系范围
Spat = rcs.SpatRef(32650)
UseFeature = Spat.GetUseFeatureInWGS84()
MapF.AddFeature(UseFeature)

洛的地理研学
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
基于通用滞回模型的GMA加工系统滑模变结构控制
06-23
作为精密加工的核心驱动源—超磁致伸缩致动器(GMA)具有推力大、应变大、响应快等优点;然而超磁致伸缩材料的滞回特性阻碍了超磁致伸缩致动器的广泛应用。基于物理或数学方法建立的超磁致伸缩滞回模型都比较复杂,具有...
英特尔gma950显卡驱动 官方最新版
07-03
英特尔gma950是一款实用性非常广的显卡，此驱动程序可以有效的解决显卡不兼容和显卡不能正常安装的问题，让电脑能够正常的进行使用，欢迎有需要的朋友下载使用！官方介绍全新英特尔#174;图形媒体加速器950（英特尔#...
gma 2.0.9 (2024.05.12) | GmaGIS alpha-2 更新日志
weixin_42155937的博客
05-22 467
请使用几何自检测方法：io.ReadVector().GeomSelfTest() 替代。链接：https://pan.baidu.com/s/1P0nmZUPMJaPEmYgixoL2QQ?链接：https://pan.baidu.com/s/1HCalK2iTF3BEIGyHUK7b_g?：修复一个针对自定义坐标系出现的WKID获取失败而导致坐标系无法预览的bug（但不影响使用）。：新增坐标系查询、创建、保存、读取等功能。：地图制图、气候气象、遥感指数、栅格代数。：内容工具、预览工具。
gma 2.0.2 (2023.10.15) 更新日志
weixin_42155937的博客
10-15 244
添加 AutoSimplify 自动简化图层（默认自动简化，如果不需要可以关闭）（与绘制栅格数据集ForcedZoom强制缩放参数类似），以大幅提高超大矢量数据的绘图速率（小概率会自动简化失败，可按需处理）。将 DataFrame 转换为矢量图层Layer，其中WKT列将被识别为几何体（若存在），其他列则为属性表字段。计算所有要素的几何属性。包括长度（当前坐标系）、面积（当前坐标系）、折点数量、几何体数量、几何中心或重心。将 矢量图层（属性表、几何体和几何属性等）转换为DataFrame。
基于 gma 绘制古代洛阳 5 大都城遗址空间分布地图
weixin_42155937的博客
03-17 1436
gma 是一个基于 Python 的地理、气象数据快速处理和数据分析函数包（Geographic and Meteorological Analysis，gma）。
gma 2 教程（二）数据操作：8.矢量操作
weixin_42155937的博客
09-17 262
gma继承了OGR全部的矢量数据读取支持，并按照矢量数据特征进行了归类。gft模块的VectorSummary类，完整统计了gma内所有栅格格式。矢量要素（Feature）是gma矢量数据的最小组成单元，其由一个几何体及其对应的字段构成。数据资源（DataSource）由多个矢量图层组成，一般的矢量文件仅含有一个矢量图层。矢量图层（Layer）是gma矢量数据的基础，其由一个或多个矢量要素构成。gma矢量数据操作类主要包括。矢量格式控制类方法/属性统计。（查看并修属性表编码方式）图 29 矢量数据组成图。
gma 地理空间绘图：(1) 绘制简单的世界地图-3.设置地图框
weixin_42155937的博客
03-06 886
gma 地理空间绘图：(1) 绘制简单的世界地图-3.设置地图框
gma 2 教程（三）坐标参考系统：1.坐标系和坐标参考系统模块简介
weixin_42155937的博客
10-10 262
坐标参考系统是地理空间数据表示和位置定位的基础，它是一种用于描述和测量地球表面位置的标准化框架。其定义了坐标系统、基准面和坐标单位等要素，以确保地球上不同地方的位置可以一致、准确地表示和比较。本章以gma坐标参考系统模块（crs）为基础，详细介绍gma内部坐标参考系统的定义、使用过程，为后续空间绘图、空间数据处理做支撑。
基于gma库计算1990-2022年吐哈盆地栅格SPEI
WangLei_1999的博客
12-11 1483
基于gma库计算1990-2022年吐哈盆地栅格SPEI
基于 gma 的栅格运算思路与应用
weixin_42155937的博客
10-06 438
提供栅格读取、处理等操作过程，并提供标准化的栅格处理思路，方便自主进行栅格数据运算等操作。返回 array，方便结合 numpy、scipy、sklearn、pytorch 等其他库进行分析运算。为方便运算后结果生成 GTiff 时添加属性，这里记录一下原始数据属性。（以 GEBCO 2020 年海洋与陆地中国范围内地形 GTiff 格式数据为例）请联系：Luo_Suppe（微信号）。
gma.rar_gma_gma算法
09-23
**gma.rar_gma_gma算法** 遗传算法（Genetic Algorithm，简称GA）是一种借鉴生物进化原理的全局优化方法，由John H. Holland在20世纪60年代提出。GMA，即“Genetic Matrix Algorithm”，是遗传算法的一个具体实现或...
Gma.QrCodeNet
08-08
《Gma.QrCodeNet：二维码编码库的深度解析与应用》 Gma.QrCodeNet，这是一个专注于二维码编码的开源库，适用于.NET框架。它提供了高效、灵活且易于使用的API，使得开发者能够轻松地在自己的应用程序中集成二维码...
GMA302KU Datasheet
02-28
**GMA302KU 数据表详解** GMA302KU 是一款高性能的G-Sensor（重力感应传感器），广泛应用于移动设备、...参考提供的“GMA302KU-中文简洁版.pdf”和“GMA302KU Datasheet V1.0.pdf”文件，可以获取更详尽的技术信息。
golang 接口
m0_70982551的博客
07-24 965
接口定义了一组方法，这些方法没有具体的实现。接口类型的变量可以存储任何实现了这些方法的类型的值。// 更多方法...在Go语言中，接口值是指存储了实现某个接口的具体类型值的变量。接口值由两部分组成：1.动态类型：这是接口值当前存储的实际类型。2.动态值：这是实际存储的值，该值必须实现了接口中定义的所有方法。
Swagger使用Map接受参数时，页面如何显示具体参数及说
a1058926697的博客
07-26 518
后端使用Map接受参数，要求在swagger页面上显示具体的参数名称、类型及说明。当Map接受参数数量少时，可以使用Swagger自带的注解。自定义注解 @ApiGlobalModel。编写处理注解对应的插件。
springboot之自动装配
weixin_50153914的博客
07-23 658
Spring Boot 通过一系列注解和条件配置实现了自动装配机制，使得开发者无需手动配置大量的 XML 文件或 Java 配置类。自动配置机制利用文件中的自动配置类，并结合条件注解和组件扫描，实现了灵活且强大的自动装配功能。这样，开发者可以专注于业务逻辑的实现，而无需处理繁琐的配置细节。
Java的@DateTimeFormat注解与@JsonFormat注解的使用对比
最新发布
訾博（ZiBo）的博客
07-27 595
在Spring和Jackson框架中，日期和时间格式化是一个常见需求。注解主要用于Spring的表单绑定，而注解则用于Jackson的JSON序列化和反序列化。了解这两个注解的使用场景和方法，可以帮助开发者更高效地处理日期和时间。和是处理日期和时间格式化的两个重要注解。主要用于Spring MVC的请求参数绑定，而主要用于Jackson的JSON序列化和反序列化。了解它们的使用场景和功能，可以帮助开发者更高效地处理日期和时间格式化需求。通过本文的介绍，希望读者能够更清晰地理解和。
GMA302KU：16g三轴数字加速度计技术规格
GMA302KU还配备了一个中断引脚，可以根据需要设置极性和为推挽输出，便于系统根据特定阈值或事件进行实时响应。这一功能对于需要实时监控和快速响应的应用非常有用，例如运动检测或者跌落检测。 GMA302KU是一款集高...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

洛的地理研学

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值