问题标签 [proj]

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 投票
0 回答
317 浏览

r - 如何让 rgdal/proj 4 为 mac 工作

我在尝试在 R 中使用 adehabitatHR 包时遇到问题。当我加载 rgdal 包时,我收到以下错误消息。

我还没有找到适合计算机编程技能不高的人的答案(例如如何为 rgdal 升级 proj4)。关于这个主题的现有答案不是以外行可以理解的形式编写的——例如,我不知道如何从其他线程中链接的 GitHub 代码中“编译”、“克隆”等。这些答案预设了(在我看来是)大量不适用于我的研究领域的编程知识。如果有人能以直截了当的方式向具有基本计算机知识的人解释它,我将不胜感激。

0 投票
1 回答
1237 浏览

excel - proj.xll 加载项添加到 Excel 后不会加载

我已经搜索并搜索以找到无法正常工作的 .xll excel 加载项的解决方案。

加载项是 proj.xll ( https://github.com/jbuonagurio/proj.xll )

说明只是说;

将 proj.xll 复制到 %USERPROFILE%\AppData\Roaming\Microsoft\AddIns。

转到 Excel 选项、加载项、管理 Excel 加载项。单击浏览,然后选择 proj.xll。

之后我加载我的电子表格并收到此错误;

“proj.xll”的文件格式和扩展名不匹配。该文件可能已损坏或不安全。除非你相信它的来源,否则不要打开它。还是要打开它?

我单击“是”,然后我得到带有“ LËÝbÉÜRichcÉÜPEL ”等乱码文本的单元格。

一些网站表明缺少加载项所需的 .dll 文件,但我不确定这是什么 dll 文件,如何找出我需要的 dll 文件以及如果找到了 dll 文件应该放在哪里.

我试图让这个加载项在多台计算机和安装的不同版本的 excel 上工作。

如果有人能提供帮助,那就太好了。

谢谢,

0 投票
0 回答
143 浏览

java - Proj4j 投影结果非常不准确

我正在尝试使用 EquidistantAzimuthalProjection 根据其相对距离和与已知位置的方位角来获取新位置的纬度/经度。这是我的源代码:

由于我已经将单位设置为米,而 10m 是一个非常小的距离。结果不应离起点的纬度/经度太远(43.694195,-79.262262)。然而 inverseProject 的结果是:

文档明确指出 inverseProject 的输出单位是度,输入单位与 Projection 本身一致(已经设置为米)。我在调用这个库时错过了什么吗?

非常感谢您的经验

0 投票
0 回答
310 浏览

r - R coord_map 与 coord_proj proj4 正射投影

试图重现

https://en.wikipedia.org/wiki/Central_African_Republic#/media/File:Central_African_Republic_-_Location_Map_(2013)_- CAF -_UNOCHA.svg

我的第一个挑战是重现没有伪影的正交投影:

coord_map在球体边缘附近的 shapefile 中返回一些伪影(例如中国。(直线而不是弧线)

这可以通过使用 proj4 来解决coord_proj吗?

有没有办法给地球/海洋一个背景颜色?

0 投票
1 回答
835 浏览

geopandas - geopandas 投影问题

我现在正在处理shp文件,但遇到了投影问题。让我在下面给你我的代码。

将 pandas 作为 pd 导入 geopandas 作为 gpd 从 geopandas 导入 GeoSeries、GeoDataFrame 导入 o​​s

这是我的数据构建

数据构建

数据构建

当我执行时print(base_map.crs), print(Aelly_to_crs.crs),我得到如下结果。

结果

数据构建

上图是执行Aelly.plot(ax=base_map.plot()) 的结果,可以看到两张图不匹配。我该如何解决这个问题?

- - - - - -编辑

我想要的输出如下图。

在此处输入图像描述

这是我的数据链接:http: //blog.naver.com/khm2963/220929301892

下面的图片是下载我的苍蝇的过程

在此处输入图像描述 在此处输入图像描述

0 投票
1 回答
1705 浏览

react-leaflet - 如何在 react-leaflet 中使用 proj4lefleat?

我有一个工作反应传单地图,我需要重新定义以保存自定义投影(ETRS89 / UTM 区域 32N,EPSG:25832),因为传单发送的默认 EPSG:3857 WMS 请求正在生成一个移位地图(所以市政府我正在使用的数据告诉我使用他们的原生投影)。

我可以像这样导入 proj4leaflet 并定义投影:

当我将此添加到我的地图时:

之后我的 WMSTileLayer 和 TileLayer 都不会渲染。调试这个的最好方法是什么?请注意,由于我使用的是开放式政府数据 WMS,因此我无法访问我正在查询的服务器,因此我看不到那里的日志。

0 投票
1 回答
743 浏览

r - 通过源码安装 rgdal

堆栈上已经有很多这样的问题......我发现了很多,但它并没有解决我的问题。

这些是我迄今为止尝试过的事情。访问 CRAN 网站 cran.r-project.org/web/packages/rgdal/index.html

我们看到为了构建这个包,我们需要

  • GDAL >= 1.6.3,来自 trac.osgeo.org/gdal/wiki/DownloadSource 的库和
  • 来自 download.osgeo.org/proj/ 的 PROJ.4 (proj >= 4.4.9); 

此外,我们需要依赖项

  • 取决于:R (≥ 3.3.0)、方法、sp (≥ 1.1-0)

但是这些已经安装好了。

我们已经安装了系统依赖项,因为

但是,当我们尝试通过

configure.args 选项来自 rgdal 包安装中给出的答案

但是我们在安装时仍然收到错误...这是日志的一部分

请注意,此处触发的警告显然没有显示停止器(如https://groups.google.com/forum/#!topic/batchjobs/TvwTbqrWFyU上所说)所以似乎系统依赖项(GDAL 和 PROJ)是两者都满足...包依赖项(方法,sp)也在那里。但是,安装无法定位某些文件。我从 stackoverflow 知道,这样您就可以为 R 的 install.packages() 命令提供某些参数...如上面的代码片段所示。所以它可能是我没有指定足够的路径......或者它的其他东西

Atlast,旧edgenode上gdal和proj的位置由下式给出

  • gdal:/usr/local/lib/libgdal.so.20
  • 项目:/usr/local/lib/libproj.so.0 11:26:32

你们能告诉我从这里去哪里吗?

0 投票
0 回答
126 浏览

r - 问题从 R Leaflet 中的 Javascript 传单地图重新创建坐标转换

底部是 Javascript 代码,它设法在荷兰坐标系 (Rijksdriehoekstelsel) 中添加地图。我试图在 R 的 Leaflet 中重新制作它,但由于某种原因,我只得到一个灰屏......

这些选项中的大多数都适用于我拥有的这个不同的 Tilelayer,但是那个在某些分辨率 lvl 上很古怪,所以我真的需要让这个工作。

proj4def、分辨率、平铺、平铺选项、投影边界、crsClass 和代码都是相同的......唯一不同的是 Project Unproject 部分,但我不知道 R 中的等价物是什么,我也不知道认为这很重要。

0 投票
2 回答
5535 浏览

r - 如何在 ggplot2 中正确绘制投影网格数据?

多年来,我一直使用 usingggplot2绘制气候网格数据。这些通常是投影的 NetCDF 文件。单元格在模型坐标中是方形的,但根据模型使用的投影,在现实世界中可能并非如此。

我常用的方法是首先在合适的规则网格上重新映射数据,然后进行绘图。这引入了对数据的小修改,通常这是可以接受的。

但是,我认为这还不够好:我想直接绘制投影数据,而不需要重新映射,因为ncl如果我没记错的话,其他程序(例如)可以这样做,而无需触及模型输出值。

但是,我遇到了一些问题。我将在下面逐步详细说明可能的解决方案,从最简单到最复杂,以及它们的问题。我们能克服它们吗?

编辑:感谢@lbusett 的回答,我得到了这个包含解决方案的好功能。如果您喜欢它,请点赞@lbusett 的回答

初始设置

我们为每个模型单元创建了两个数据框,一个带有模型坐标,一个带有真实的经纬度交叉点(中心)。

可选:使用较小的域

如果您想更清楚地看到单元格的形状,您可以对数据进行子集化并仅提取少量模型单元格。请注意,您可能需要调整点大小、绘图限制和其他便利设施。您可以像这样子集,然后重做上面的代码部分(减去load()):

如果你想完全理解问题,也许你想尝试大域和小域。代码是相同的,只是点大小和地图限制发生了变化。下面的值适用于大的完整域。好的,现在让我们开始吧!

从瓷砖开始

最明显的解决方案是使用瓷砖。咱们试试吧。

这是结果: 在此处输入图像描述

好的,现在更高级的东西:我们使用真正的 LAT-LON,使用方形瓷砖

在此处输入图像描述

好的,但那些不是真正的模型方块,这是一个黑客。此外,模型框在域的顶部发散,并且都以相同的方式定向。不太好。让我们自己投影正方形,即使我们已经知道这不是正确的做法……也许它看起来不错。

在此处输入图像描述

首先,这需要很多时间。不能接受的。此外,再次:这些不是正确的模型单元。

尝试点,而不是瓷砖

也许我们可以用圆点或方点代替瓦片,也可以投影它们!

在此处输入图像描述

我们可以使用方点...和投影!我们走得更近了,尽管我们知道它仍然不正确。

在此处输入图像描述

不错的结果,但不是完全自动的,并且绘图点不够好。我想要真实的模型细胞,它们的形状被投影改变了!

多边形,也许?

因此,正如您所看到的,我正在寻求一种正确绘制模型框的方法,该模型框以正确的形状和位置投影。当然,模型中为正方形的模型框,一旦投影出来,就不再是规则的形状了。所以也许我可以使用多边形并投影它们?我尝试使用rasterToPolygonsfortify关注这篇文章,但没有这样做。我试过这个:

在此处输入图像描述

好的,让我们尝试用 lat-lons 代替...

在此处输入图像描述

(对不起,我改变了图中的色标)

嗯,甚至不值得尝试投影。也许我应该尝试计算模型单元角的纬度,并为此创建多边形,然后重新投影?

结论

  1. 我想在其原生网格上绘制投影模型数据,但我无法这样做。使用瓦片是不正确的,使用点是骇人听闻的,并且由于未知原因,使用多边形似乎不起作用。
  2. 通过 投影时coord_map(),网格线和轴标签是错误的。这使得投影的 ggplots 无法用于出版物。
0 投票
0 回答
69 浏览

geospatial - 从 GCP 创建仿射或其他类型的投影并转换为 proj4 字符串

我想知道像http://www.georeferencer.com这样的工具是如何制作的,如果可以采用一组 GCP,创建转换矩阵并将其转换为 proj4 字符串,这样它就可以用作投影ol3 以与他们相同的方式在邮件顶部显示静态图像。