问题标签 [maptools]

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 投票
2 回答
3472 浏览

r - 如何将夏威夷和阿拉斯加添加到 R 中的空间多边形?

如何将夏威夷和阿拉斯加添加到以下代码(取自 Josh O'Brien 的回答:Latitude Longitude Coordinates to State Code in R)?

latlong2state(ak)and代码返回 NA ,latlong2state(hi)但如果代码修改正确,阿拉斯加和夏威夷将作为结果返回。

任何帮助表示赞赏!

0 投票
1 回答
1181 浏览

r - 通过子集得到一个空的 SpatialPolygonsDataFrame?

我希望通过属性对 SpatialPolygonsDataFrame 进行子集化,但我想允许它返回一个空的 SpatialPolygonsDataFrame。

如果我们要将 SpatialPolygonsDataFrame 类型的对象像 data.frames 一样处理,如这里所讨论的,我们应该能够获取并使用空对象。

我很感兴趣,因为我想将它合并到一个函数中,该函数可能会尝试通过一个本质上不会选择任何特征的属性来进行子集化。

有没有办法使这项工作?如果没有,有没有办法初始化一个空的 SpatialPolygonsDataFrame 对象?我想在这样的对象上执行的未来操作涉及在现有地图上过度绘制,所以我希望无论如何都能生成图像,即使是空白的。

0 投票
1 回答
1224 浏览

r - 单个 spatialPolygonsDataFrame 行的子集多边形

我想拆分spatialPolygonsDataFrame由多个多边形几何组成的单行,但我不确定实现这一目标的最明智的方法。例如,如果我们想从usa下面移除阿拉斯加,我们可以过滤掉质心超过 50 度纬度的多边形。有谁知道一种相对无痛的方法来做到这一点?

0 投票
1 回答
1257 浏览

r - 将纬度/经度或 URM 坐标添加到 R 中的 shapefile

我是 R 中的 GIS 新手,我正在尝试将 lat/long 或 UTM 坐标添加到 shapefile。我从这里下载了芝加哥市的边界 (City_Boundaries.shp):http ://www.cityofchicago.org/city/en/depts/doit/supp_info/gis_data.html

我加载了 maptools 和 rgeos 库:

我将数据带入 R 并尝试为 16T 区添加 UTM 代码:

但是,结果数据对我来说没有意义。当我查看“city”中的“coords”槽时,坐标值是,例如,X=1092925 和 Y=1944820。我使用了一个外部工具(http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html)来查找 GPS 坐标。结果是:lat=17.511,long=-81.42。这是牙买加和洪都拉斯海岸之间的某个地方。在我看来,shapefile 坐标存在于它们自己的宇宙中;正如 shapefile 的名称所暗示的那样,坐标提供了对城市形状的准确描述,但这些坐标不会自动映射到地球上。

我的最终目标是确定是否在芝加哥发生了一些地理标记事件。我很乐意将按纬度/经度列出的这些点转换为 UTM,如下所示:

如果这更有意义,我也愿意尝试以原始格式处理事件数据。最终,我需要帮助将芝加哥城市边界变成一个多边形,我可以对照我的事件数据进行检查。像这样的东西:

如何将我的 shapefile 和事件数据放入一个通用(且准确)的参考框架中,以便我可以检查点是否在城市中?您可以提供的任何帮助将不胜感激。

0 投票
1 回答
329 浏览

r - 栅格未裁剪为 shapefile

我正在尝试从 a 中合并一些状态shapefile,并生成raster我可以使用的 a downstream。我已经合并了状态,但是当我创建一个空栅格以使用裁剪功能进行栅格化时似乎失败了。我是新手GIS featuresR非常感谢您的帮助。 Shapefile is from http://www.arcgis.com/home/item.html?id=f7f805eb65eb4ab787a0a3e1116ca7e5

当我然后调用'Corners.crop'我收到的范围时:

我对让这个工作所缺少的东西感到困惑。我也希望有一个 1Km 的分辨率,并且很好奇在空栅格上或在我栅格化形状之后更改分辨率是否会更好。

0 投票
1 回答
883 浏览

r - R将等高线转换为高程图

我希望能够从 R 中的等高线创建一个高程图。我对使用形状文件非常陌生

目前我已经从这里下载了数据, 这些数据提供了适用于整个英国的 .shp 文件。

它还提供了等高线,总结了英国的拓扑结构。

对于高程图,我想要一个data.framedata.table多个均匀间隔的点(彼此相距 100m)来产生一个数据输出,给出一个 x、y 和 z 值。其中 x 和 y 代表纬度和经度(或东向和北向),z 代表高度(以米为单位)。

我认为可能有一些工具会自动为您执行插值,但不确定它如何处理地理空间数据。

这是我的基本开始...

0 投票
0 回答
62 浏览

r - 导航 SpatialPolygonsDataFrame

这可能是一个简单的问题,但是,我想查看地图上的每个输出区域,从这里下载了寄宿生。选择 2011 OA 并获得完整的分辨率,领域的范围。

我已经对其进行了压缩...并使用以下命令读取了 shp 文件

如何绘制或查看每个输出区域的坐标?

0 投票
0 回答
285 浏览

r - 形状文件的 R 交集

我有一个英国道路网络的形状文件。我还有一个连接到车辆的 GPS 跟踪设备的经纬度。我可以为每次旅行将 GPS 经纬度转换为 SpatialLine。然而,GPS 跟踪设备中存在一些错误的数据点。每个连续的经纬对彼此相距约 2 分钟。在地图上绘制经纬度时,很容易在视觉上注意到“错误”点。有没有办法通过评估创建的空间线与道路网络形状文件的交集来以编程方式执行此操作?即,鉴于我们知道这些点应该大致遵循道路的形状,是否有任何 R 包中的函数能够评估这一点?

我知道gInterserct,但我不只是在寻找真假......结果......我试图找出一系列点中的哪些点是“错误的”......

0 投票
1 回答
583 浏览

r - R空间映射:使用类区间函数创建区间

我是 R 新手,目前正在尝试创建等值线图。我的问题是我的休息没有意义。我想做的是创造相等的休息时间:

59% 以下 60-69% 70-79% 80-89% 90% 以上

然而我得到的是:
低于 60% 60-70% 70-80% 80-90% 超过 90%

这是否意味着 values=80 在这两个类别中都有特色?

我创建这些中断的代码是:

0 投票
1 回答
325 浏览

r - 从 maptools 执行 checkPolygonsHoles 后,多边形中的孔消失

注意:在 Edzer Pebesma 的建议下,这个问题已被交叉发布到 R-sig-geo,在这里,它得到了一些很好的回应。


我遇到了以下意外结果checkPolygonsHoles

这里一个恼人的副作用是原始对象wrld_simpl也被改变了。这个结果在我看来就像一个错误,还是我错过了什么?

PS:之前shape_with_hole编辑的对象checkPolygonsHoles,继续表现奇怪: