《地理坐标（经纬度坐标）和屏幕坐标（xy坐标）间的转换》的读后笔记

VIP文章恒宝乐园

已于 2023-07-21 14:51:16 修改

阅读量2.8k

点赞数

分类专栏： GIS 文章标签：几何学线性代数 gis

于 2022-03-12 08:13:44 首次发布

本文链接：https://blog.csdn.net/hengbao4/article/details/123437848

版权

今天读了一下《地理坐标（经纬度坐标）和屏幕坐标（xy坐标）间的转换》，觉得内容很好，特摘录和附上自己的解释。

背景

在我们的屏幕上，有一张地图，这张地图经过缩放、平移、旋转，最终地理坐标和屏幕坐标的关系大致如下图所示：

这种关系要怎么描述呢？我们可以假设地图是一张纸，而屏幕是一堵墙。只要我们有两个图钉，我们就能把纸定在墙上。我们把这两个点称为锚点。锚点在屏幕坐标系上的坐标是(x1,y1)和(x2,y2)，对应在地理坐标系上的坐标是(lon1,lat1)和(lon2,lat2)。

那现在的问题就变成了，已知两个锚点的坐标，

（1）地理坐标转屏幕坐标：已知任意一点的地理坐标(lon,lat)，求它在屏幕上的坐标(x,y)

（2）屏幕坐标转地理坐标：已知任意一点的屏幕坐标(x,y)，求它的经纬度坐标(lon,lat)

lon 表示 longitude (纬度)
lat 表示 latitude (经度)
x,y 表示屏幕上的只是一个任意位置的坐标点

在一个小范围内（例如是方圆几公里内），我们可以假设地面是平的，而不是弯的。如果经纬度都用弧度表示，那么1纬度方向对应的长度是：

$\cdot lat$

注意：此处的 lat 指的不是我们常用的纬度度数，而是纬度弧度。半径乘以弧度才能得到弧长！
$\frac{度数}{360}=\frac{弧度}{2π}$
故弧度 = 度数 $\cdot \frac{ \pi }{180}$

其中 $R$ 是地球半径。而相同经度间的距离会随着纬度的增加而减少，

在这里插入图片描述

在 $l a t$ 这一纬度下，1经度方向对应的长度是：

$\cdot \cos{(lat)} \cdot lon$

上式中的 $\cdot \cos{(lat )}$ 就是上图中的某一纬度下的切面的半径，然后在半径乘以弧度。

那么，(lon,lat)这个坐标平面化后的坐标就是： $\cdot lon \cdot \cos{(lat)},R \cdot lat)$

故，此坐标可以作为我们在某一个小的平面范围内的类似于x、y的坐标

由于我们是假定在小范围内纬度是不变的，故此时若我们有两个经度不同的点，求经度长度的差时，应该遵循：

$dlen_{lon} = R \cdot lon_2 \cdot cos (lat) - R \cdot lon_1 \cdot cos (lat) = R(lon_2 - lon_1 )\cdot cos (lat)$

其实上式中的求纬度长度差时用到的 $l a t$ ，使用 $lat_1$ 或者是 $lat_2$ 都对，因为我们在计算经度长度时视作一致。

如果求一个纬度长度差与经度长度差比值的话：

$\frac{dlen_{lat}}{dlen_{lon}} = \frac{R \cdot lat_2-R \cdot lat_1}{R(lon_2 - lon_1 )\cdot cos (lat)}$

好，我们上下同时除以 $R$ ，可得

$\frac{lat_2- lat_1}{(lon_2 - lon_1 )\cdot cos (lat)}$

又因为前面提到过，我们计算经度长度时视纬度为不变的，故我们随意地，把 $lat_2$ 代给 $l a t$ .

$\frac{lat_2- lat_1}{(lon_2 - lon_1 )\cdot cos (lat_2)}$

这里需要注意的是：

关注

专栏目录