问题标签 [spatial-data-frame]

I would like to find an efficient way to combine certain character + numeric column values in a list of SpatialPolygonsDataFrame objects. Here is reproducible data:

What I would like to do is add another column that is a combination of the FIPS, REGION, and SUBREGION columns, separated by an underscore ('_'). I know how to add+name a new column to each SPDF object in the list as done in the loop below, but I don't know how to get the desired column row entry:

The row entries for the new unique.id column would be in this format: FIPS_REGION_SUBREGION. For example, for the ATG polygon feature in spdfl[[1]], I would like the 'unique.id' column to have an entry like this:

Please advise on how to do this for all features in the SPDF list.

我有一个带有经纬度坐标的data.frame：

和一个国家（西班牙）的边界，

我想df只选择.border

我怎样才能做到这一点？

注意：在我的可重现示例中df，第一个入口点在内部，第二个显然在外部，第三个在外部但靠近海岸。

我在 R 中加载了一个 shapefilereadOGR并命名为data.map

通过使用，我可以很好地看到形状

问题是多边形的 ID（绘图的 ID）与@dataID 不同。我的意思是，当我跑步时：

然后运行

图中绘制的多边形不是数据行所指的多边形。

我目前正在尝试将北卡罗来纳州的所有邮政编码映射到传单地图上。最终，我想用另一个数据集的值对它们进行颜色编码。

虽然我已经设法使用这些值将标记放在地图上，但以彩色方式使用邮政编码这样做确实会为我的项目增加价值。

我在 R 中下载了 USAboundaries 库，所以我认为我有足够的数据在传单中创建多边形，但数据不是 SpatialPointsDataFrame 格式。

USAboundaries 数据集的示例输出如下所示：

在我看来，这应该足以转换为 SpatialPointsDataFrame，但我不知道如何在 R 中做到这一点！任何提示将非常感谢！

我想从空间中删除具有某些值的行，data.frame因为标准 R似乎没有通过类对象grep()的所有槽传播。sp

bd@data[- grep("xcluded", bd@data$Notes),]

它给了我一个错误：

错误：试图从不是 S4 对象的对象（类“data.frame”）中获取插槽“数据”

我读到sp.na.omit可以删除NA但不知道如何删除具有特定值的行（例如，删除在 bd@data$Notes 中具有“排除”的行）。有什么建议么？

抱歉，我无法通过 R 创建 shapefile 示例，但我会尝试提供有关我的 sp 数据框的更多信息：

> class(bd)

1 “空间线数据帧”

attr(,"包")

1个“sp”

bd 属性

类：空间线数据帧

特点：8855

范围：172.6811、174.5966、-36.36374、-34.42634（xmin、xmax、ymin、ymax）

crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

变量：134

名称：OBJECTID_1, Name_in_us, Unique_ID, RAMM_Road_, iSequence, SH, OBJECTID, ROADID, RoadLength, RoadNameAn, Displaceme, road_id, AvgWidth, pave, CJNEX_urba, ...

最小值：1, 014-0000, 1188496, 0, 0, No, 0, 0, 0, 012-0132 (805), 0-1008m, 0, 0, 混凝土, 农村, ...

最大值：5299, ZEALANDIA ST, 2044000001, 2065, 480, Yes, 999, 1683, 45826.7818765, ZIDICH ROAD, 9925-9966m, 3520, 18.2, Unsealed, Urban, ...

df(bd@lines) 内的折线

$行[[1000]]

“线条”类的对象

插槽“线”：

[ 1 ]

“线”类的对象

插槽“坐标”：

[1,] 174.3629 -35.77290

[2,] 174.3627 -35.77281

[3,] 174.3624 -35.77276

bd@data[1000,]

OBJECTID_1 Name_in_us Unique_ID RAMM_Road_ iSequence SH OBJECTID ROADID RoadLength RoadNameAn Displaceme road_id AvgWidth ...

1000 1000 域 RD 1.16e+08 116 2 否 89 116 94.75686 域 RD (116) 0-95m ...

解决方案：1.作为下面的答案/评论 2.bd[- grep("xcluded", bd$Notes),]改用

我想在 R中将空间数据框投影到 EPSG 25833，但 QGIS 似乎不知道（为了重现性，我使用在他/她对这个问题的回答中创建的代码 jazzurro稍作改动）

现在，当我将 shapefile 加载到 QGIS 中时，它不知道投影。

有任何想法吗？

我正在尝试使用 R 程序将 SpatialDataFrame (shapefile) 的两列合并为一列，但是在两列中都有空格，当它们与名称加 NA 一起时，但是我希望 NA 不会出现在我的新专栏。我使用了粘贴功能。像这样的东西：

obs。我的原始数据不是这个并且有更多的列

我用这个：

1.

事实上，这将是期望的结果：

我有一个数据表，其中包含大约 1700 列和 100 行数据，带有唯一标识符。它是调查数据，一个组织的每个员工都回答了相同的 9 个问题，但它被编译成每个组织的一行数据。python/pandas 中是否有一种方法可以垂直整合这些数据，而不是它已经在 x 轴上的细长格式？我目前正在剪切和粘贴。

我在 R (no2df) 中有一个数据框，它看起来像这样：

由于上面的矩阵表明每四个点（具有相同的 id）形成一个多边形（或者在我的情况下是一个像素或单元格），单元格的边界显示在“latb”和“lonb”列以及中心单元格在“lat”和“lon”列下重复四次，因为它是每个单元格的相同中心。最后，最重要的是，我在“no2tc”列下描绘了变量 NO2，它遵循与 lat 和 lon 相同的推理。

问题如下：

我应该怎么做（请给出步骤）将此数据框转换或操作为空间数据框，以便最终我能够导出到包含具有变量“no2tc”作为单元格属性的单元格的 ESRI shapefile，所以我可以在 QGIS 中查看吗？

我设法从 HDF 文件中提取 MODIS 土地覆盖数据并将其放入栅格中。

我想把它放到一个数据框中，但是从原来的正弦曲线中重新投影+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs，我认为。

对于与我正在分析的其他数据集兼容的普通椭圆体/WGS84。

这是我尝试过并且似乎有效的方法：

但是，当我随后将我的土地覆盖数据放入这个新投影中的数据框时，所有土地覆盖值都变为 NA。

这就是数据框在正弦投影中的样子（4 是土地覆盖分类）

通过我的重新投影，它看起来像这样：

关于我做错了什么或如何让土地覆盖坐标正确重新投影的任何建议？

干杯！！！！

我还读到有一个 NASA MODIS 重投影工具，但它不再存在/可用。有人知道吗？