问题标签 [pyproj]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
1016 浏览

python - pyproj中的datum和ellps有什么区别?

我正在使用从 ECEF 转换为 lla

我不确定ellpsdatum关键字有何不同。文档指向我这里: http: //proj4.org/parameters.html#parameter-list但仍不清楚它们有何不同?

0 投票
0 回答
362 浏览

python - georasters pip 在 Windows 上安装失败

遵循https://github.com/ozak/georasters中给出的步骤,但对于 pip 安装过程的两个步骤具有相同的回溯。“收集 pyproj ..”时因“msvc.py”出现问题而爆炸。

有人可以帮忙吗?谢谢

0 投票
1 回答
750 浏览

python - 从 CA 状态平面坐标到 Lat、Long 的 PyProj 转换不如预期

我正在尝试使用 pyproj 将 Calfornia State Plane Zone 5 XY 坐标转换为纬度、经度。

为了验证,我知道 CA Zone 5 XY 坐标 (6559361.78613, 1834842.95456) 是地址 13422 Ankerton St, Bassett, CA。纬度、经度为 (34.0342403°, -118.0076074°)。

但是,输出为(34.65142393815357°, -65.96879755500356°),位于大西洋某处,远离加利福尼亚 5 区。

0 投票
2 回答
142 浏览

python - 增加numpy 2D数组的列

我有这个玩具示例,它使用 numpy 和一个名为 PyProj 的外部包。Lat 和 lon 是二维数组,包含一些特定领域信息的坐标。我想要做的是从我任意选择的中心点计算球体上的距离。lat_0 和 lon_0 的形状是

但是 API 调用 inv 不喜欢这样。我收到运行时错误 -

它需要一系列形状

所以我希望 lat_0 和 lon_0 的形状与 lon 和 lat 的形状相同,所有常量值都是中心纬度和经度。增加 lon_0 和 lat_0 的列并用中心值填充它以使其与 lon 和 lat 的形状相同而不使用 for 循环的最有效方法是什么?

0 投票
0 回答
611 浏览

python - 最小二乘法给出错误的高度,LLA 到 ECEF 到 Python 中的 LLA

目标是从已知 GPS 位置的地面站找到移动单元的位置,向移动单元报告模糊或有噪声的范围。为此,我希望脚本将从谷歌地图获取的纬度经度高度 (LLA) 转换为最小二乘算法可以使用的 XYZ 坐标,然后将 XYZ 猜测/结果返回为谷歌地图上准确的 LLA 坐标。当我运行脚本时,纬度和经度似乎相当准确,但高度错误得令人难以置信,我无法弄清楚原因。我知道 ECEF 到 LLA 的转换是基于 alt 的,所以如果我能找出它为什么如此不准确,我也可以提高纬度/经度的准确性。

有谁知道为什么高度如此错误以及我该如何解决?

0 投票
1 回答
744 浏览

python - 大圆到点的距离

我有三点:A、B 和 C(WGS84 中的纬度、经度坐标)

我想计算 C 和 A、B 之间的最小路径之间的距离。(这必须适用于任何纬度、经度,所以我无法进行预测)

pyproj 或geographiclib.geodesic 可以为我做到这一点吗?

编辑:不重复,因为我想在 WGS84 中工作,而不是在球体上工作。正如我所说,这必须适用于任何经纬度。我知道这个问题没有简单的解决方案,这就是为什么我问是否有一种简单的方法可以使用 pyproj 或 geodesic 之类的库来做到这一点。

0 投票
1 回答
2487 浏览

python-3.x - 在pyspark中将纬度和经度转换为UTM坐标

我有数据框包含每个点的经度和纬度坐标。我想将每个点的地理坐标转换为 UTM 坐标。

我尝试使用 utm 模块(https://pypi.org/project/utm/

但我收到此错误:

更新:

创建应用 utm 或 pyproj 函数的 udf 后

结果是:

0 投票
1 回答
12214 浏览

python - 如何使用 Python 获取 GeoJSON 多边形的面积

我有 GeoJSON

http://geojson.io显示为

在此处输入图像描述

我想用 Python 计算它的面积(87106.33m^2)。我怎么做?

我试过的

它给出了1.1516745933889345e-05-233827.03300877335第一个没有任何意义是意料之中的,但是我该如何解决第二个呢?(我不知道如何设置pyproj.Projinit 参数)

我想epsg:4326它是 WGS84 (来源)是有道理的,但对于epsg:3857我不确定。

更好的结果

以下更接近:

它给出了 87254.7m^2 - 与 geojson.io 所说的仍然有 148m^2 不同。为什么会这样?

0 投票
2 回答
1736 浏览

python - 计算从地面看到的卫星的天顶角和仰角

我有卫星和地面观察者的纬度、经度和高度。我正在尝试计算从地面看到的卫星的卫星天顶角和卫星方位角。

我目前正在尝试使用astropy. 该文档描述了如何计算太阳的天顶角和方位角。

我试图将卫星直接放在观察者上方,这样天顶角应该是 0°(或仰角 90°)。

但是,显然.transform_to(AltAz(...))并不意味着我认为的意思,因为答案不是我想要的:

我预计仰角接近 90°,而不是 63°。

pyorbital中也有一些例程,但它似乎不包含我正在寻找的内容。

我相信它也应该可以使用pyephem,除了这似乎期望我从目录中获得卫星信息。我已经有了纬度,经度,高度。

我也看过,pyproj但我认为这仅适用于大地水准面的物体。

如何计算从地面观察到的卫星的高度和方位角?

0 投票
1 回答
1080 浏览

python-3.x - 使用 Cartopy 和 Pyproj 重新投影 GOES16 卫星数据

我正在尝试使用 cartopy 或 pyproj 重新投影 GOES16 Full Disk 图像。我想把它变成一个不同的投影。对于此示例,我尝试将其重新投影到墨卡托。然而,当我运行代码时,我得到了一个完整的全球数据图像,而不是墨卡托投影,也没有任何 cartopy 特征。我觉得我错过了一个简单的步骤,但无法弄清楚它是什么。下面是一个可重现的示例 - 我使用的是 Python 3.5。