标签[sf] - 堆栈内存溢出

将多边形 class sf 转换为 class owin in r - transform a polygon class sf to class owin in r

我的多边形具有以下结构：我正在尝试使用spatstat.geom package 将t1转换为 class owin object： t1中的坐标是纬度/经度：我如何将t1投影到 UTM 才能获得我的owin object？！ ...

tmap 中按日期列的颜色点 - Color points in tmap by date column

我正在尝试使用tmap tm_dots通过时间变量为一组地理点着色，但它一直将时间变量视为分类变量。有没有办法将POSIXct列用作连续变量？演示该问题的一些示例代码是：library(sf) library(tmap) df = data.frame(lon=1:20,lat=1:20,ts= ...

更新 sf linestring 的（第一个和最后一个）坐标 object - update (first and last) coordinates of sf linestring object

我有一个带有 LINESTRING 几何（~ 100000 行）的大型 sf 对象。在这个 object 中，我想更改每个线串的第一个和最后一个坐标/点。我知道我可以通过将线串转换为点，更改每个线串 ID 的第一行和最后一行/点的几何形状，然后转换回线串来做到这一点。但这似乎有很多额外的（不 ...

过滤并排除 stat_density2d_filled 中的级别 - Filter and exclude a level in stat_density2d_filled

我有类似 map 美国某些地点的代码和数据：我有这样的 output：我希望排除一些级别，例如 (0,0002,0004] 和 (0,0004,0006)。此外，我想手动为级别着色。我有此代码scale_fill_manual(values = c("khaki1", "khaki2", "k ...

使用 SF 和 ggplot 到 map USA state map 有问题 - Having problem to map USA state map with SF and ggplot

我有一些点（纬度和经度）是这样的：要map这个，我有这样一段代码：代码1 但是，结果没有正确显示分数我可以用这段代码 plot 它但我正在寻找像上面的代码一样的 map 有一条曲线 map 代码 2 你能修复代码1吗？ ...

`devtools::check()` 期间失败的示例 - Examples failing during `devtools::check()`

由于这不是我第一次在开发一些小的 R package 时遇到dplyr::filter()与存储在/data/下的 sf object 的问题，我想现在可能是寻求帮助的时候了。上次执行devtools::check()时测试出现问题，现在示例执行出现问题。在这两种情况下，测试和示例在手动执行时都 ...

如何将 SpatRaster 的比例尺单位从 km 更改为 mi？ - How can I change the units of the scale bar from km to mi for a SpatRaster?

我正在尝试使用 SpatRaster 底图（CRS 26913 - NAD83，UTM ZONE 13N）在 ggplot 中创建一个 map。当我尝试使用 function annotation_scale 添加比例尺时，scale_bar 表示的距离以千米为单位显示正确。但是，当我尝试将 ...

不存在简单的要素几何：返回表的 data.frame 或 tbl_df 问题 - No simple feature geometries present: returning a data.frame or tbl_df problem with a table

我的大学有一项任务，我需要从表格中转换坐标。因此，我有一张坐标类型错误的表格，代码可以正常工作，并且这些坐标已成功转换。一位朋友给我发了一张坐标正确的新表，但当我尝试对其进行转换时，R 工作室控制台给我这条消息：不存在简单的特征几何：返回 data.frame 或 tbl_df 我无法使用 st ...

R 的坐标对 (UTM) 之间的距离，单位为公里 - Distance between pairs of coordinastes (UTM) in km with R

我有一个带有 2 对 UTM (32N) 坐标的 dataframe，我需要计算它们之间从起点到终点的公里差。我正在尝试使用 sf 库，使用“by_element”，但我收到一条错误消息“st_distance(data, by_element = TRUE): .missing_y non è ...

将 stat_density2d_filled 中的级别数字减少到更低的数字 - Decreasing level's numbers in stat_density2d_filled to lower numbers

我有很多点到map（密度）。我想将 map 中的级别数修改为较低的数字，因此我不会在 map 中使用多种颜色，而是使用较低的数字：要设置边框和 plot 我使用此代码的数据：正如您在 plot 中看到的那样，我有 11 个级别，我如何降低级别以在我的地图（3 或 4）中具有较低的 colors ...

从 R 中的列表中提取 sf 对象的名称 - Extract the names of sf objects from a list in R

我正在尝试提取列表中对象的名称，它们都是 sf。所以我不能使用名字。 NULL 我需要向量中的 object 个名称来替换列名称。我几乎尝试了一切都无济于事。 object_names <- c() ...

如何在 r 中的光栅 map 之上正确组合多边形 map？ - How to properly combine a polygon map on top of a raster map in r?

我想在使用库 ggmap 获得的光栅 map 之上组合美国一些州的形状文件ggmap 。这是我的代码这就是我得到的：还收到以下警告，我不确定这是什么意思：如您所见，US state shapefile 和栅格 map 之间存在不匹配。shapefile map us.state.long ...

如何从 geom_sf 生成的星号 map 可视化中删除随机对角线？ - How to remove random diagonal lines from star map vizualizations produced by geom_sf?

语境我正在开发一个名为 starBliss 的自定义星星 map package 。虽然我认为我已经成功地从我的问题中弄清楚了我需要什么，但看起来有些边缘情况开始出现问题。有关更多背景信息，请在此处查看 GitHub 问题问题我目前正在研究一个 function，它采用星座线 shapefil ...

使用 R 的 sf 查找两个不同数据集中的几何图形之间的交集 package - Find intersection between geometries in two different datasets using R's sf package

使用从该网站下载的 shapefile，我试图获取几何图形的（区域）交集，其中每个几何图形都在不同的文件中。但是，当我尝试该操作时出现错误，我一直无法解决。下面的代码有助于说明我遇到的问题：上面的交集操作不起作用；它抛出一个错误说我该如何解决？ ...

将网格剪辑到 shapefile 但保持单元格完整 - Clip a grid to a shapefile but keep cells whole

我有两个 shapefile：一个六角网格，一个边界 Shapefile：我正在尝试将 Grid 剪辑到 Shapefile 但无论我如何尝试st_intersect() 、 st_overlap() 、 st_contains()等，结果始终相同：我想要的是剪辑网格并保持剩余的单元格完整，如下 ...

如何在 map 上重新定位 position (ggplot/sf/geometry) - How to relocate the position on the map (ggplot/sf/geometry)

我用R画了一个中美map。我想 map 跨越这两个国家的“整体思维”。你能帮我解决以下问题吗？搬迁这两个国家？两国之间的空间太大了。我想缩小两国之间的距离，放大两国的面积。 “阿拉斯加”的大小更小？我认为阿拉斯加的面积相对于美国的面积来说太大了。这是我的代码。这是 plot。在 ...

无法使用 sf() 和 rnaturalearth() 提取 R 中的湖泊几何数据 - Unable to pull lake geometry data in R with sf() and rnaturalearth()

我在尝试使用rnaturalearth package 提取大湖数据时收到一条错误消息。我已经使用相同的代码一段时间了，直到今天才收到错误消息。我不确定该站点是否已关闭，但这种情况已经持续了几个星期。我在他们的 github 上发帖，但没有收到回复。其他任何人都可以复制此错误或知道如何解决它或 ...

优化代码以将大向量拆分为较小的文件？ - Optimise code to split large vector into smaller files?

在我下面的玩具/示例代码中，我制作了一个覆盖世界的网格，然后使用此网格将一个大型复杂的全球数据集拆分为每个网格单元格的一个文件。在我的实际工作中，这是一个瓶颈，需要很长时间。我会很感激一些关于优化它的想法和想法。我通过利用并行处理取得了一些成功，但我也认为它可以做得“更聪明”。 ...

当我尝试使用 sf package 中的汇总 function 来溶解多边形的边界时，R 返回错误。如何解决？ - R returns an error when I try to dissolve polygons' boundaries using the summarize function in the sf package. How to solve it?

我在sf package 中使用 R 的summarize function 来消除一些多边形的边界，但我遇到了问题，因为它在某些情况下会返回错误。举一个具体的例子，当我试图在第 80 届美国国会的形状文件中溶解加利福尼亚州 state 的地区边界时，我得到以下错误：使用的代码如下：它在其他 ...

支持 terra 和 sf 中的 vsi7z（直接读取压缩的 7z 文件）？ - support for vsi7z (direct reading of compressed 7z files) in terra and sf?

目前的terra和sf版本不支持直接读取7zip压缩文件，因为GDAL还不支持。但似乎该功能将包含在定于 2023 年 5 月 1 日发布的GDAL 3.7.0 版中（ https://github.com/OSGeo/gdal/milestone/42 ）。知道terra和sf目前正在使用 ...