问题标签 [pyproj]

我正在尝试在 Ubuntu 16.04 上运行 anaconda python 2.7.15 的远程服务器上安装 pyproj。运行后（在https://anaconda.org/conda-forge/pyprojconda install -c conda-forge pyproj列出的第一个选项），包安装成功（pyproj 2.0.2）。但是，在运行我的程序时，我得到以下回溯：

在我的 python 脚本中，我使用 导入from pyproj import Proj, transform，相关的代码行是：

显然，pyproj 无法找到data目录。搜索目录 ~/anaconda2/lib/python2.7/site-packages/pyproj，我确认没有data目录。我有另一台服务器，在 anaconda3 上安装了 pyproj，它确实包含一个data目录（带有一个epsg文件）。我将此目录复制到 anaconda2 服务器，并按照回溯中的说明使用环境变量声明路径：

我仍然得到相同的错误回溯。

我data在 github repo 上找不到 pyproj 的目录（https://github.com/pyproj4/pyproj/）。

如何解决此错误？

我正在尝试使用 PyProj 构建一个 CRS 实例，以将地理形状与图像边界相匹配。使用以下代码

我收到错误消息：

除了Proj4 wiki 上显然不再存在的文档页面之外，似乎没有很多关于什么使投影有效或无效的文档。谁能帮我理解这里出了什么问题？代码（不是我写的）似乎在过去有效。

我正在尝试自己解决我的问题，但我不能，我正在尝试以您可以想象的每种格式运行此代码，并且在 ArcGIS pro 软件中它是相同的我在任何其他中都找不到此错误消息问题。从类似的问题来看，似乎某些数据文件可能丢失了？

我得到这个错误

我也在 arcgis pro 中尝试过，我得到了同样的结果

我正在使用带有 LCC 投影的 matplotlib-basemap。

使用 .proj4string 后，我在地图数据中发现 x_0, y_0 具有默认值

我想将 x_0, y_0 设置为 0。

但我找不到关于假东和假北的变量

如何将 x_0,y_0 值设置为 0？

https://matplotlib.org/basemap/api/basemap_api.html - 没有关于 x_0、y_0 的参数数据

我的代码在 .py 时运行良好，但是当我使用 pyinstaller 创建 .exe 文件时，它显示 ModuleNotFoundError。

我在制作.exe时写了这个

我怎么解决这个问题？

我读了

pyinstaller 和 pyproj 的问题

我将尝试在 Pyinstaller 的“hooks”文件夹中制作并定位 hook-pyproj.py。但是 hook-pyproj.py 已经有了相同的代码。

之后，我安装了“底图”并尝试从中使用 pyproj。但是它向我显示了同样的错误。

这是我执行 .exe 文件时显示的错误。（我用*****覆盖了一些信息）

这是我的代码。

我一开始试过这个。

但我得到同样的“ModuleNotFoundError: No module named 'pyproj._datadir'”错误。

我正在尝试使用 Pyproj 将坐标转换为经度和纬度，但是我收到一个我无法解决的错误：

b'初始化列表中没有参数'

知道如何解决这个问题吗？

非常感谢您的帮助！

I am trying to plot a map of Spain using Geopandas and Matplotlib. I am using the GeoJSON file you see below in which the geometry is in geographical coordinates (EPSG 4326)

I would like to plot the map in Mercator projection but when I try to convert the units using

I receive an error saying:

I have followed the instructions in http://geopandas.org/projections.html. There it says that if the geopandas dataframe does not count with the information of the CRS, you must define it. But this is not the problem here.

This is the code I have so far

If I run the following command:

I receive the following response (which means that the geopandas DataFrame have the information about CRS)

Then

Returns an error

What am I missing? Thanks in advance.

我想以米为单位计算地球表面两点之间的距离

我已经尝试过底图和cartopy，但两者都导致不同的数字。

底图：

卡托普：

其中 coord0 和 coord1 都是大小为 2 的 numpy 数组，具有坐标的纬度和经度。

卡托普输出：990.6094719605074

底图输出：1072.3456344712​​142

在计算多边形的面积时，与 google earth engine 或 geojson.io 相比，使用 shapely/pyproj 的结果大不相同。我认为这与投影有关，在研究它时，我发现了诸如（如何使用 Python 获得 GeoJSON 多边形的面积）之类的示例，这些示例建议使用 Alders 等面积投影，但这些似乎不是一个选项在 PyProj 的 2.1.3 版本中。我也想知道谷歌地球引擎和geojson.io是否有一种自动方式来选择区域适当的区域投影给定多边形的区域，但我不确定他们将如何选择这些。

例如：

45573.993884405005返回m^2的面积。

Google Earth Engine 回归23944.14737277293，Geojson.io 回归23997.77

我有一个包含十进制度的纬度/经度坐标的数据框。

我的目标是在 1 平方公里的矩形网格上汇总数据。为此，我根据将纬度、经度转换为距赤道的距离（以千米为单位）并四舍五入到最近的千米中描述的方法将坐标转换为千米

该方法包括计算从参考点到点 (lat=0, lon) 和 (lat, lon=0) 的距离。

但它不起作用，因为它似乎取决于参考点。

通过将我的参考点设为 (lon_ref=mean(lon), lat_ref=mean(lat))，我最终聚集在彼此相距 120 公里的相同瓦片点中。

这是我正在使用的代码：

这适用于一些纬度/经度，但不适用于其他人，

例子：

使用上面的代码计算距离和圆：

两点将聚合在一起 但是，两点之间的距离是 85 公里！

我怎么解决这个问题？鉴于我的数据框中有 1000 万个纬度/经度，还有其他有效的方法来进行聚合吗？