问题标签 [epsg]

如何使用 R 从 PROJ4 字符串接收 EPSG 代码？

我有一堆栅格数据集（随机分布在世界各地，北半球和南半球），坐标参考系统 UTM WGS84。现在我想使用 R 找出每个数据集的特定 EPSG 代码。你有什么想法吗？

也许可以以动态方式创建 EPSG 代码，因为 EPSG 代码似乎总是北半球的 <326> 和特定区域编号的组合，例如 <32>（对于德国），这会导致EPSG 代码 <326> + <32> = <32632>。

不幸的是，我发现您似乎从南半球的 <327> 开始并添加特定的区域编号，例如 <32>（用于尼日利亚），结果为 <327> + <32> =<32732>。

现在我不确定这个理论是否正确（需要验证，包括来源），以及是否有更简单的方法可以从 PROJ4 字符串接收 EPSG 代码。

我发现的有用链接：

网址：https ://spatialreference.org/ref/epsg/?page=9&search=wgs+84+utm

网址：https ://epsg.io/32632

网址：https ://de.wikipedia.org/wiki/UTM-Koordinatensystem#/media/Datei:Utm-zones.jpg

编辑：

非常感谢，Explorer

I want to reproject my local-geoJson file in epsg:28992 to the OSM EPSG:4326 in OpenLayers. I have the feeling I am close to the solution but I don't know what the next step is. I've tried and looked for multiple examples here at SO but I've the feeling I'm missing a certain line of code somewhere.

For now the code below shows me my local gjson file at null-island. How do I tell Open Layer to reproject it to the Netherlands?

thanks in advance.

我必须在 Python 中编写一个函数，该函数采用坐标列表和坐标参考系统（作为 EPSG 代码）。如果坐标在坐标系中有效，则返回 True，否则返回 False。我怎样才能做到这一点？

我的想法是获取最小和最大纬度和经度，创建一个粘合多边形并检查点是否落在里面。问题是我不知道如何从 EPSG 代码中获取最小和最大纬度和经度。我唯一的想法是将它们手动写入函数，这是毫无意义的。

这是正确的方法还是我想多了而有更简单的方法？

到目前为止，我有：

但这毫无意义，因为函数应该采用任何坐标和任何 EPSG

我正在尝试将数据帧中的 sf 点（CRS：4326，基于度数）转换为本地度量投影（EPSG：23867，DGN95 / UTM 区域 47N）。

所以我有一个数据框“航路点”，其中包含我的所有属性以及 lon / lat 作为列。我在数据框中创建了一个 sf 点几何，其中“lat”和“lon”作为输入（waypoints_sf）。（另外我将坐标保存为 Northing 和 Easting 在两个新列中）

现在，当我尝试对这个数据帧进行 st_transform 时，它不会引发错误或警告，但无论我尝试什么，sf 点都是空几何c(NaN, NaN) 。

那么我做错了什么，也许我真的很想念简单的东西，我真的不知道为什么它不起作用？如果它不与 st_transform 一起使用可能使用 N 和 E raw，那么简单的解决方法是什么？

现在对我有用的解决方法（“航点”是简单的df）：

我需要将 .tiff 文件从 EPSG:4326 - WGS 84 - Geographic 重新投影到 EPSG: 32638 (UTM Zone 38N)。当我运行下面的代码时，输​​出文件是一个重新投影的 .tiff 大小的原始文件的一半。

我正在尝试将 shapefile 绘制为 Folium 地图的点，但我正在拉扯头发以使投影正确。数据在 EPSG:28992 系统，RD 新坐标系中。这是一个投影坐标系，因此它以米为单位给出坐标。要在folium 地图上绘制它，我需要将其转换为 EPSG:3857，如果我是正确的，它适用于度数。不幸的是，在用函数转换它之后，我完全是胡说八道to_crs()。

首先我加载数​​据：

print(WTG.head())显示了这一点：

哪个是对的。它还具有正确的 CRS，如以下输出所示print(WTG.crs)：

现在要将其转换为其他 CRS，我使用以下简单代码：

但是现在检查结果print(WTG.head())对于几何列来说完全是一派胡言。

几何点应该以度为单位。换句话说，坐标应该在 52 和 5 左右。print(WTG.crs)虽然的输出是正确的：

你们中有人知道如何解决这个问题吗？还是我错过了什么？提前致谢！

我收到此错误：FutureWarning: '+init=:' syntax is deprecated。':' 是首选的初始化方法。进行更改时，请注意轴顺序更改：https : //pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6 return _prepare_from_string(" ".join(pjargs ))

我从一个嵌入了谷歌地图的网站上抓取了美国城市的数据。我的数据没有crs，所以我做了研究，发现谷歌地图的epsg是3857。所以我把它设置为3857

最后，我想创建一张标有城市的美国地图。我知道该怎么做。我拥有的地图位于 epsg:4269 中。这就是为什么我用以下方法更改数据的crs：

然后，当我使用以下命令创建 excel 文件时：

points_within = gp.sjoin(gdf, US, op='within') writer = pd.ExcelWriter('locationsUS3.xlsx', engine='xlsxwriter') points_within.to_excel(writer, sheet_name='Sheet1', index=False)
writer.save()

excel文件为空。

不是当我不执行该gdf = gdf.to_crs({'init': 'epsg:4269'})步骤时，而是我收到“用户警告：正在加入的帧的 CRS 不匹配！” 错误

那么这个 FutureWarning 错误描述了什么，我该如何分别解决这个问题，我如何在这个 excel 文件中获取我的数据？

感谢您在这里帮助我！

PS：'US'是来自美国的数据

我想用它L.imageOverlay来显示全球温度，我当前的图像如下图所示：

但是当我使用传单覆盖它时，它显示如下：

也许是因为错误的投影，这是我的代码：

我认为问题是从 Equirectangular 到墨卡托方形投影的转换，但我不知道如何解决这个问题。

这是我的代码：

它会生成这个图像： image，它没有显示背景图并且看起来有点失真。

我的坐标数据最初是使用 RD 坐标标准 (EPSG:28992) 制作的，我使用以下代码将其转换为 EPSG:4326：

通过与一些在线转换器进行比较，我对经纬度输出进行了完整性检查，输出指向正确的位置。

我尝试遵循此解决方案：https://gis.stackexchange.com/questions/348339/using-crs-epsg3857-but-misalignment-between-stamen-background-and-coordinates-o以防我的转换缺少“towgs84 “部分，但图像看起来仍然相同，但颜色略有不同。

我得到了带有 ["x","y"] 列的数据框。

“x”和“y”是 epsg:6991 以色列网格的坐标

例如我喜欢转换x=[212513], y=[617119]坐标，

至： epsg:4326 wgs 84

x=[35°7'2.734"] , y=[31°38'56.771"]

如何在两种类型的坐标之间进行转换？