GIS中的投影坐标系

投影说明

GIS操作过程中,不可避免的涉及到处理数据的问题

而数据中有一个极为重要的东西就是其空间参考(见下图)

在这里插入图片描述

地理坐标系与投影坐标系

这里给出一些重要概念的简单定义

地理坐标:就是用经纬度表示地面点位的球面坐标。

地理坐标系进行地图投影后就变成了投影坐标系

地图投影:是按照一定的数学法则将地球椭球面上的经纬网转换到平面上,使地面点位的地理坐标(φ,λ)与地图上相应的点位的平面直角坐标(x,y)或平面极坐标(δ,ρ)间,建立起一一对应的函数关系,能够实现由地球椭球面向地图平面的科学转变。

简单来说,坐标系数据或地图的属性,而投影坐标系的属性

img

地理坐标系经过投影后变成投影坐标系,立体向平面的转变,投影坐标系因此由地理坐标系投影组成,投影坐标系必然包括有一个地理坐标系

一个坐标系可以有投影也可以没投影,但有投影坐标系就一定有地理坐标系

img

投影坐标系的应用

由上介绍我们不难得出这样的结论

要计算要素的面积长度这样的信息时我们需要做的一个点就是必须应用投影坐标系进行计算

在本文的第一张图里可以看到数据属性的空间参考中只有地理坐标系,这显然印证了该数据未进行投影处理

当然,在GeoScene Pro中,我们可以通过直接右击相应数据属性查看其对应投影信息(如下图)
在这里插入图片描述

分析该图,同样不难发现其投影参数

投影方法:墨卡托投影

中央经线:0°

每单位米数:1

应用过众多GIS软件,大家都清楚GIS软件提供了众多投影坐标系,那是不是每一个投影坐标系我们都能用

国际上通用最常见的地理坐标系WGS84

我国基本比例尺地形图除1:100万采用兰勃特投影(Lambert)外,其他均采用高斯-克吕格投影

为减少投影变形,高斯-克吕格投影分为3度或6度带投影。

下图用的是WGS_1984_UTM_Zone_46N投影坐标系,中央经线93°,相对于上面的图,这张图**“歪了”**,这就是不同投影带的差别,靠近中央经线的地区形变最小,越远形变越大,所有根据经纬度选择对应的投影带也是很重要的哟

img

选用投影坐标系

经过如上的讲解,接下来我们就不得不解决一个问题(如何应用最合适投影坐标系呢)

国际通用的UTM等轴墨卡托投影

计算UTM投影带很简单,计算公式如下:

北半球地区,选择最后字母为“N”的带,南半球“S”

带号=(经度整数位/6)的整数部分+31

例如汶川县东经102°51′~103°44′

103/6≈17+31=48

即:WGS_1984_UTM_Zone_48N

当然,计算显然比较麻烦,这里直接放个表

img

国内特有的CGCS2000

具一些投影坐标系的例子

①CGCS2000_3_Degree_GK_CM_105E

②CGCS2000_3_Degree_GK_Zone_35

③CGCS2000_GK_CM_105E

④CGCS2000_GK_Zone_18

对应的说明

①3°分带下的35带投影(横坐标没有带号)

②3°分带下的35带投影(横坐标有带号)

③6°分带下的18带投影(横坐标没有带号)

④6°分带下的18带投影(横坐标有带号)

一般,在我国标准地形图上,X坐标的前两位代表所在分度带带号,其余的表示X坐标

显然,很难很复杂,那3°和6°之间如何计算呢

3度带带号=(经度+1.5°)/3取整

6度带带号=(经度+3°)/6取整

(🐔该处转换公式此前因 笔者疏忽 出现了一个 little mistake,感谢好友 🌐 于YXS 的指正)

太复杂了,直接上图

img

GIS软件中的投影

GeoScene Pro以及其他GIS软件里面有三个常用的工具与坐标有关的工具,分别是(如下图)
在这里插入图片描述

①定义投影

此工具对于数据集的唯一用途是定义未知坐标系。使用此工具前,必须已获知该数据集的正确坐标系

在这里插入图片描述

注意:这里定义的可以是栅格数据也可以是矢量数据,只要它本身没有坐标,但是你知道它正确的坐标,就可以用这个工具去定义。 什么时候会出现这种情况呢,比如,CAD的dxg等数据加载到ArcGIS里面是没有坐标系的,但是在CAD画图的时候是有坐标系的,这时候就可以用这个工具去定义正确的坐标系,所以当你拿到一个CAD文件的时候就一定要问画图小哥哥(可能是学测绘的)用的是什么坐标系。

②投影

将空间数据从一种坐标系投影到另一种坐标系。矢量数据专用。 注意哦,这里的矢量数据必须要有坐标系。
在这里插入图片描述

③投影栅格

栅格数据集从一种投影变换到另一种投影。同样的,这里的栅格数据必须要有坐标系。
在这里插入图片描述

当数据(无论矢量、栅格,如果是一张已出图的图像,这又是另一种情况,后续具体举例)没有坐标系又知道其正确坐标系时,用“定义投影”工具定义坐标系

当矢量数据需要转换坐标时,用“投影”工具

当栅格数据需要转换坐标时,用“栅格投影”工具

### GIS中的高程投影坐标系 #### 高程投影的概念及其重要性 在地理信息系统(GIS)中,除了水平位置外,还需要考虑垂直方向的位置信息即高度。为了精确表示三维空间内的物体或地形特征,在二维的投影坐标系基础上引入了第三维度——高程。这种带有高程坐标的系统称为高程投影坐标系[^1]。 #### 使用方法 当涉及到具有海拔差异的数据集时,如地形图、DEM(数字高程模型),就需要采用合适的高程参考框架。通常情况下,会先定义一个基准面作为零点参照物,比如平均海平面;接着利用特定算法计算其他地点相对于该基准的高度值。对于具体的软件操作而言: - **ArcGIS**: 用户可以在创建要素类时指定Z轴单位以及对应的垂向坐标参考系统(VCS),也可以后期编辑属性表添加相应的字段存储高程数值。 - **QGIS**: 支持导入含z值得矢量文件(.shp,.geojson等), 并允许设置项目CRS(Project CRS)以确保不同来源之间的兼容性和准确性。 ```python from osgeo import ogr,osr driver = ogr.GetDriverByName('ESRI Shapefile') dataSource = driver.Open('example.shp', 0) layer = dataSource.GetLayer() spatialRef = layer.GetSpatialRef() # 获取当前的空间参考对象 new_srs=osr.SpatialReference() new_srs.ImportFromEPSG(4978)# WGS84 geographic 3D coordinate system EPSG code transform = osr.CoordinateTransformation(spatialRef,new_srs) for feature in layer: geom = feature.GetGeometryRef() transformed_geom=geom.Clone().Transform(transform) ``` 这段Python代码展示了如何使用GDAL库读取shapefile并将其几何结构从原始坐标系转换到WGS84大地测量系统的例子。这里选择了EPSG:4978代表了一个包含椭球体参数在内的完整的三维坐标体系[^2]。 #### 常见问题及解决方案 1. 数据不一致:如果两个数据源使用的VRS不一样,则可能导致叠加分析出现问题。解决办法是在加载前统一所有层所依赖的标准; 2. 显示异常:某些可视化工具可能无法正确渲染带高程的信息。此时应检查是否启用了支持3D视图的功能选项,并确认已安装必要的插件; 3. 计算误差:由于地球并非完美的圆形,因此即使在同一套规则下也可能存在细微差别。可以通过选择更贴近实际地貌特性的局部优化方案来减小偏差影响范围。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值