问题标签 [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 投票
1 回答
3658 浏览

r - 读取带有字符编码的 shapefile

有没有办法读取具有特定字符编码的形状文件?我正在尝试阅读加拿大 shapefile,其中某些名称中包含特殊(法语)字符。我可以手动转换它们,但如果在某个地方有一个我很盲目的设置,我不希望这样做。

0 投票
1 回答
256 浏览

r - 在 R 中使用 for 循环绘制小倍数

我正在尝试绘制一个地图小型多重网格,该网格显示自 1900 年以来与佛罗里达州相交的飓风/热带风暴。我使用一些空间查询对该项目的所有大西洋风暴的数据库进行了子集化。

我现在正在佛罗里达州、一些毗邻的州、佛罗里达州的一些主要城市,当然还有奥基乔比湖的多边形顶部绘制我有限数量的飓风轨迹的线shapefile 。这是简单的代码:

我遇到的三个大问题:

1)循环给了我每场风暴的地图。我希望代码在网格中为每年(即使在没有风暴的年份)和当年的所有风暴生成佛罗里达/南部地图,最好带有标签。

2)我想在所有风暴中设置风速的颜色,而不仅仅是循环中每个特定行中的那些。这样,即使是一年中唯一的风暴,强风暴(如 1992 年的安德鲁)也会显得更暗。也许我可以通过重新编码一个可以相应设置样式的分类(H1、H2 等)变量来解决这个问题。

3) 假设我能找出第一条,我很难在每条风暴路径上渲染标签。maptools 文档不是很好。

无论如何,这是目前的输出(标题是 shapefile 中两个字段的串联):

在此处输入图像描述

真正的问题是第 1 号。提前感谢您的帮助。

0 投票
1 回答
1369 浏览

r - R:在 shapefile 中查找/匹配 NUTS 代码

我在哪里可以在 R 的 shapefile 中找到 EU NUTS 代码?

NUTS = 统计地域单位命名法

每个 NUTS 地区都有一个 NUTS 代码:

德国的示例图片.

您还可以在此 Excel 文件中找到详细的 NUTS 代码。

当我下载 NUTS shapefile 时,我找不到这个 NUTS 代码。

我从以下位置下载了 shapefile NUTS_2013_20M_SH.zip

ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts#nuts13

下载后,我使用以下代码将 shapefile 上传到 R 中:

笔记:

  • dsn =是保存 shapefile 的文件夹的路径。
  • layer =是不带文件扩展名的文件名(例如,.shp)

如果我使用summary(NUTS_shape),我只找到一个“NUTS_BN_ID”但没有 NUTS 代码。此外,如果我使用str(NUTS_shape[4461,])查看示例项目 4461 的列表结构,似乎没有 NUTS 代码。

有谁知道 NUTS shapefile 是否包含 NUTS 代码,或者 NUTS_BN_ID 是否以某种方式链接到此代码?

0 投票
3 回答
2132 浏览

r - 经纬度集合的边界多边形

我有一个表,其中包含一个名为 queryResult 的城市中某些位置的所有纬度和经度,我执行以下操作:

1 - 获取城市的栅格地图[例如布莱克浦]

这将返回地图上的以下点

在此处输入图像描述

现在我想绘制所有这些经纬度的外边界[城市边界线],类似于以下预期结果

布莱克浦市多边形

更新 1:提供输入数据并应用建议的ahull解决方案:

我应用了@spacedman 和@cuttlefish44 建议的 ahull 解决方案,得到的结果与预期的多边形大不相同:

将 ahull 应用于数据

您可以从以下链接下载包含所有纬度和经度的 .csv 文件:Blackpool Lat,Lon

Google 建议的区域边界如下所示:

在此处输入图像描述

0 投票
2 回答
3052 浏览

r - 在 R 中使用带有“maps”包的 xlim/ylim 投影会产生更大的地图

我正在使用 R 和 maps/mapproject/maptools 包来绘制一些地图,并注意到一种对我来说似乎违反直觉并且实际上限制了我能做的行为。

在不指定任何投影的情况下绘制欧洲地图(限制取自ETRS89 / ETRS-LCC,因此没有冰岛,也剪裁在东部):

未指定投影的欧洲地图

结果与预期的一样,正在使用限制并且生成的地图遵循它们。

根据帮助,地图默认使用的投影是:

这对我的需要来说不是一个很好的投影,我将使用 LCC 投影,其平行线如上面的 spatialreference.org 链接所示:

指定投影的欧洲地图

结果出乎意料,因为它包括更大的区域(向北走很远,包括俄罗斯),基本上使地图无法使用。

奇怪的是,在使用网格时,会清楚地考虑原始限制:

带有 LCC 投影和网格的欧洲地图

我想要的(以及我认为是上述代码的结果)是一个矩形裁剪,其中包括我指定的限制(由于使用的投影进行了调整,因此预期的面积比上面的矩形更大) . 此外,每当使用 map() 进行投影时,地图和边界都会出现空白区域。

问题是:在使用 map/mapproj/maptools 时有没有办法得到这个结果?我试图人为地更改 xlim/ylim 没有好的结果,因为它似乎在很大的间隔内工作(即更改它们不会产生效果,直到突然一半的欧洲随着下一次递减而消失)。

0 投票
1 回答
5924 浏览

r - 图例范围和颜色分布的 spplot 问题

我的情节和图例的正确颜色范围有问题。

这是我使用的代码:

这是包含数据的文件:https ://www.sendspace.com/file/hjtatp (压缩后通常重 57 mb)

这里开始地图(但地图有次要优先级,可以跳过)

这是没有任何比例修改的样子: 在此处输入图像描述

所以一切都是蓝色的。显然有很大的距离,从最小值到最大值。我想修复标度,例如最后一个值将“高于 x”。我试着这样做:

在此处输入图像描述

所以现在这看起来好多了。我是这样做的:

所以我手动添加了 at 值(但不是准确的比例)。一切看起来好多了,但是..

如您所见,右侧的比例不是均匀分布的。我看不到任何蓝紫色,只有橙色和黄色。

地图上的一些地点也是亮黄色(德国地区),因为这里的值最高,但遗憾的是比例尺上没有这种颜色。

可能我没有正确地做到这一点。我不知道如何将比例设置为看起来不错。我想有这样的规模:

在此处输入图像描述

我通过添加:

但同样,这只是假秤,它不会起作用。

第二个快速问题:如何在散点数据之上添加国家等高线?因为现在轮廓被隐藏在彩色数据下:c

0 投票
1 回答
1318 浏览

r - Leaflet addPolylines 没有绘制数据对象的所有值(或线)

我正在使用在生成的 web 地图addPolylines上覆盖特定值的轮廓,但是,并没有绘制它提供的所有数据。根据第一张图像,至少有两个值为125的区域是明显的,但传单生成的地图仅显示一个区域。数据可以从这里下载,生成 SpatialLinesDataFrame 对象和生成地图的可重现脚本如下:leafletaddPolylines
addPolylinesleaflet

降水模式

c.linesSP具有三条线的坐标,但addPolylines仅绘制一个区域的等高线


具有一条等高线的 WebMap

具有两条等高线的 Web 地图

目前尚不清楚为什么addPolylines不绘制它最初提供的所有坐标。非常感谢建议。

0 投票
1 回答
2313 浏览

r - 在 R 中绘制地图点

我已经使用 maptools 成功地将 NYC PUMA 区域的 shapefile 加载到 R 中,我想在它上面绘制 55 个点,我在另一个文件中如下所示:

但是,这些点不会绘制。

首先,我这样做是为了绘制 shapefile:

它绘制 shapefile这是形状文件

然后我尝试在顶部绘制点,但无论我执行以下哪一项,它总是失败:

或者

或者

(最后一个将数据绘制在一个只是 XY 散点图而不是覆盖在 shapefile 上的新图上)

任何想法如何做到这一点?

编辑,添加文件(希望修改为工作文件!):

形状文件信息:https ://www.sendspace.com/file/wbqrpb 点文件:https ://www.sendspace.com/file/9yrrbu

0 投票
0 回答
59 浏览

r - 从 R 中的形状文件中查找沿海和国际边界

我想从形状文件中自动创建两个变量:1. 一个区域是否有国际边界的虚拟指标和 2. 一个区域是否有沿海边界的虚拟指标。

例如,对于几内亚来说,变量 1 将是下面带有红点的区域,变量 2 将是带有蓝点的区域(我是通过肉眼看到的)。

在此处输入图像描述

@data插槽中似乎没有这些类型特征的任何信息:

也许他们在其他地方(我对形状文件几乎没有经验)?是否有某个函数至少可以让我创建一个变量来指示一个区域是否没有外部边界(即上面地图中所有没有点的区域)?

0 投票
1 回答
365 浏览

r - 将空间多边形与外部数据合并并在属性上选择多边形

我正在尝试在 R 中绘制我的研究区域。到目前为止,我已经下载了人口普查 OA 边界数据,并将属性数据合并到我希望查看的 OA 中。
因此,让我们从总共 1,000 个 OA 中说。我只想查看其中的 500 个 OA,我有外部数据集。Land 是我的空间多边形数据框,其中包含 1,000 个 OA。它有一列是 OA 名称“地理代码”。但是,它有它的坐标,因为你运行的是:

它产生所有坐标。

Houseprice 是一个包含房价和带有地理编码的列的数据框。

远代码:

这表明我只有 500 个条目(这是我想要的)。

这最终绘制了整个 1,000 个 OA,而不仅仅是我想要的 500 个。

对于仅要绘制数据的 500 个 OA,我必须做什么?感谢帮助