R中具有两个几何的sf上的几何运算
在R中，sf包是一个用于空间数据处理的强大工具。它提供了许多函数和方法，用于在sf对象上进行几何运算。
几何运算是指对空间数据进行各种操作，如计算距离、相交、合并、裁剪等。sf包中的几何运算函数可以应用于sf对象，这些对象可以表示点、线、面等几何要素。
以下是一些常用的sf包中的几何运算函数：
- st_distance：计算两个几何对象之间的距离。
- 优势：可以方便地计算不同几何对象之间的距离，如点到线的距离、线到面的距离等。
- 应用场景：在空间分析中，常常需要计算不同要素之间的距离，如寻找最近邻、计算空间自相关等。
- 推荐的腾讯云相关产品：腾讯云地理信息服务（GIS）。
- 产品介绍链接地址：https://cloud.tencent.com/product/gis
- st_intersection：计算两个几何对象的交集。
- 优势：可以方便地计算不同几何对象之间的交集，如线与面的交集、面与面的交集等。
- 应用场景：在空间分析中，常常需要计算不同要素之间的交集，如寻找共享边界、计算面积重叠等。
- 推荐的腾讯云相关产品：腾讯云地理信息服务（GIS）。
- 产品介绍链接地址：https://cloud.tencent.com/product/gis
- st_union：计算两个几何对象的并集。
- 优势：可以方便地计算不同几何对象之间的并集，如线与面的并集、面与面的并集等。
- 应用场景：在空间分析中，常常需要计算不同要素之间的并集，如合并多个区域、计算总面积等。
- 推荐的腾讯云相关产品：腾讯云地理信息服务（GIS）。
- 产品介绍链接地址：https://cloud.tencent.com/product/gis
- st_difference：计算两个几何对象的差集。
- 优势：可以方便地计算不同几何对象之间的差集，如面与面的差集、线与面的差集等。
- 应用场景：在空间分析中，常常需要计算不同要素之间的差集，如提取不同区域、删除重叠部分等。
- 推荐的腾讯云相关产品：腾讯云地理信息服务（GIS）。
- 产品介绍链接地址：https://cloud.tencent.com/product/gis
需要注意的是，以上只是sf包中几何运算函数的一小部分示例，sf包提供了更多功能丰富的几何运算函数，可以根据具体需求选择合适的函数进行使用。同时，腾讯云地理信息服务（GIS）是一个强大的云计算产品，提供了丰富的地理信息处理和分析功能，可以与R中的sf包结合使用，实现更复杂的空间数据处理任务。
相关·内容
我在R中有一个相对简单的任务:我有两个多边形层，一个包含39个焦点多边形的“补丁”层和一个包含一个焦点地貌类型的多多边形的“陆地覆盖”层。我需要剪辑这两层，以便我只保留的范围，重点覆盖类型内的焦点多边形。有时，这在sf::st_intersection中很好，有时使用sf::st_difference和一个“负”的陆地覆盖层(包含所有非震源覆盖类型的范围)很好，有时这两种方法都不起作用。起初，我认为这些不同的行为取决于由此产生的地形肤色，但情况似乎并非如此。
我得到的错误是形式上的错误
CPL_geos_op2中的错误(op，st_geometry(x)，st_geometry(y))：评估
浏览 1提问于2021-03-11得票数 2
我正在建立一个州的形状文件，在那里阿拉斯加和夏威夷被描绘成位于德克萨斯州南部的某个地方，以便于绘制一张说明性地图。通过使用maptools包和来自的一些代码，我已经能够使用TIGER中的形状来完成这个任务。
然而，我现在遇到了麻烦，因为我想把城市添加到我的地图上。将我的形状制作成sp对象，然后使用maptools::elide以前在阿拉斯加很好，但是带比例的elide在集合点上的工作方式不一样，所以我的城市最终在错误的地方结束了：
library(maptools)
library(sf)
library(tmap)
library(tidyverse)
library(tidygeocode
浏览 4提问于2021-06-10得票数 0
回答已采纳
可能是一个非常基本的问题，但我在简单特性R包的文档中什么也没有发现。
我正在寻找本机sf函数来动态提取sf对象的所有列，而不需要几何图形。就像SP@data与sp对象一样。
以下函数完成了这项工作，但我更喜欢使用本机函数：
st_data <- function(SF) { SF[, colnames(SF) != attr(SF, "sf_column"), drop = TRUE]}
一个典型的用法是当我想按属性合并两个sf数据集时(merge不工作于两个sf对象)：merge(SF1, st_data(SF2))。
在这种情况下，使用st_geometry(SF2)
浏览 2提问于2017-09-06得票数 4
回答已采纳
我觉得自己缺乏直觉，不知道sf对象将如何与tidyverse工具，特别是dplyr交互。下面是一个例子。我们希望取两个点的缓冲器，然后减去较小的缓冲区，在每个点周围留下一个环/环。我提供了两种方法。
第一个中断几何列，使用map2获取环，然后将其锁定回原始的sf数据。这种方法是有意义的；我知道所有这些部件是如何工作的，但并不觉得“整洁”或整洁，因为我们创建了中间对象，我们必须手动替换crs，而且感觉似乎很容易将几何图形与属性不匹配。
第二种方法是我想出的办法，就是把它保存在一个管道里。在这里，我希望group_by和mutate能够工作--某种程度上来说--管理至少没有错误地运行，并且似乎创建
浏览 0提问于2018-04-10得票数 3
回答已采纳
我目前正试图在R3.3.3 (Windows)下实现sf::st_make_valid()函数的工作。
我的问题是，sf包不是针对liblwgeom库构建的。为了在windows下安装liblwgeom，我尝试安装osgeo4w，根据它的安装选项列出liblwgeom库。尽管如此，在通过osgeo4w安装程序安装liblwgeom之后，sf仍然不识别liblwgeom库。
关于有效功能的sf文档:只有当包与liblwgeom链接时才可用，二进制CRAN发行版目前不是这种情况；有关如何安装liblwgeom的说明，请参阅包源代码存储库。
有什么建议吗？
install.packages("
浏览 4提问于2017-09-16得票数 2
回答已采纳
我试图得到两个形状文件的交汇点(属于某些大都市区范围内的人口普查区)。我能够成功地获得相交特性，但是当我试图将sf_intersect的输出转换为SpatialPolygonsDataframe时，我会得到以下错误：
“as_Spatial( from )中的错误:不支持将功能类型sfc_GEOMETRY转换为sp”
这是我的密码：
library(sf)
library(dplyr)
library(tigris)
library(sp)
#download shapefiles corresponding to metro areas
metro_shapefiles<-
浏览 0提问于2018-02-14得票数 2
2回答
缺失美学几何
、、
我希望使用以下数据集创建一个地图：
具有以下形状的文件：
这是我的R码：
library("tidyverse")
library("dplyr")
library("sf")
library("ggplot2")
library("tmap")
library("tmaptools")
library("RColorBrewer")
data = readr::read_csv("C:/Users/amanm/Desktop/annualincomedata
浏览 11提问于2022-08-11得票数 -1
回答已采纳
1回答
我从PostGIS-DB读取了一个sf对象，如下所示：
parking <- st_read(dsn = con, query='SELECT * FROM "osm"."planet_osm_point"')
parking具有以下类：
> class(parking)
[1] "sf" "data.frame"¨
我想将data.table功能添加到属性表中。
当我这样做的时候
parking <- as.data.table(parking)
sf对象丢失了，我只有一个dat
浏览 0提问于2021-04-14得票数 1
我必须在R中形成2个简单的特性集合(一个用于地区，另一个用于区域)，我需要将它们组合成一个具有不同属性名称的R中的单个GeoJSON FeatureCollection (这样我可以在MapBox地图中过滤数据层，参见)。例如，一层可能给出每个地区的人口，而另一层可能给出每个区域的人口。这似乎与：有关
到目前为止，这就是我所拥有的
library(sf)
library(dplyr)
library(sfheaders)
library(geojsonio)
path_to_data <- "Ghana_261"
sf_use_s2(FALSE)
districts
浏览 6提问于2022-11-30得票数 1
1回答
我有一个来自的1弧秒SRTM光栅层。
我试着把那个栅格变成一组多边形，每100米高的波段。然后我想要溶解或结合相同属性的多边形在一起。
我的问题是，大多数多边形在这个过程中得到的是一个无效的几何。
我不知道问题出在哪里，所以我将描述我所做的所有步骤：
library(tidyverse)
library(sf)
library(stars)
# First I open the raster, remove NAs and set the min and max values
r <- raster::raster(my_raster)
raster::NAvalue(r) <-
浏览 2提问于2020-04-02得票数 0
回答已采纳
2回答
我有一个包含两个几何列的sf数据集。它看起来是这样的：
> trip_geo
dstid sourceid dest_geom source_geom
1 1 1 MULTIPOLYGON (((-2.607985 5... MULTIPOLYGON (((-2.607985 5...
2 1 2 MULTIPOLYGON (((-2.607985 5... MULTIPOLYGON (((-2.57022 51...
3 1 3 M
浏览 20提问于2020-03-08得票数 0
回答已采纳
以下代码： library(rnaturalearth,rnaturalearthdata)
library("sf")
#load world map and crop the rectangle I want
blank2 <- ne_countries(scale = "large",type = 'countries', returnclass = "sf")
blank1 <- st_crop(blank2 , xmin = -22, xmax = 44, ymin = 30, ymax = 66)
#
浏览 16提问于2021-07-11得票数 0
回答已采纳
2回答
我已经为这件事困惑了好几天了。我有一个栅格，我需要根据sf对象中的多边形属性进行编辑。这两个对象具有相同的范围。 填充了9的栅格 r <- raster( nrow=10, ncol=10, xmn=0, xmx=10, ymn=0, ymx=10 )
values(r) <- 9
> r
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 1, 1 (x, y)
extent : 0, 10, 0, 10 (xmin, xmax, ymin,
浏览 20提问于2020-01-11得票数 2
回答已采纳
1回答
我使用pyshp库检索形状的坐标。
sf = shapefile.Reader(r"{}".format(boundary_file))
shapes = sf.shapes()
fields = sf.fields
records = sf.records()
for record in records:
if record['NAME'] in cities_list:
city = record['NAME']
s = sf.shape(record_id)
geom = s.__ge
浏览 1提问于2021-03-08得票数 0
我在R中使用了一些空间数据，并且想知道是使用依赖于旧的空间格式(package sp)的包/函数还是使用新的包sf。我做了这个测试，基于找到的代码。
这样做的目的是“确定空间点数据集中每个点的最大距离在xx米范围内的所有点”。
library(tictoc)
# define a buffer distance and a toy data
maxdist <- 500
df<-data.frame(x = runif(10000, 0, 100000),y = runif(10000, 0, 100000),id = 1:10000)
# doing the analysis
浏览 1提问于2018-03-13得票数 7
1回答
你好啊！我正在尝试使用R和一系列纬度/经度点在世界地图上绘制多艘船的路线。
一切都很好，直到船穿过国际日期变更线(-180/180)时，路径跳过地图。
我尝试在R中应用st_wrap_dateline()函数，但它看起来只适用于具有两个点的数据集，即开始和结束坐标。
我的R代码在下面，还有一个绘图截图--任何帮助都是非常感谢的！
# Download background Blue Marble Globe image
download.file("https://www.researchvessels.org/images/nasa_base_v2.png", &#
浏览 18提问于2019-11-23得票数 4
回答已采纳
1回答
我想使用R中的sf包将我的数据组织为一个简单的特征对象，它有几个几何列，例如多边形、它的质心和它的缓冲区，所有这些都共享相同的属性。 我很难将这样的对象保存到磁盘上。例如，如果我尝试使用st_write()将其另存为shapefile，则会得到以下结果 Error in clean_columns(as.data.frame(obj), factorsAsCharacter) :
list columns are only allowed with raw vector contents
Warning:
1: In abbreviate_shapefile_names(obj) :
浏览 50提问于2021-09-09得票数 0
回答已采纳
我使用R中的sf包和arcpullr包从ArcGIS REST服务中提取数据，并将其作为sf对象处理。我遇到了MULTIPOLYGON的一个问题，其中sf只缓冲MULTIPOLYGON的一部分(即它缓冲一个多边形，但只缓冲另一个多边形的小条)。在缓冲找到的示例时，我无法复制这个问题。
这是一个MRE (对不起，如果您没有它，就必须安装arcpullr )。
# install.packages("arcpullr")
library(arcpullr)
library(sf) # redundant since loaded with arcpullr, but here fo
浏览 7提问于2021-12-07得票数 2
回答已采纳
library(sf)
DT <- data.frame(
lat1=c(-26.3, -27.1, -26.8),
long1=c(27.3, 28.5, 28.2),
lat2=c(-26.5, -26.9, -27.2),
long2=c(28.3, 27.5, 28.1)
)
DT$p1 <- st_geometry(st_as_sf(DT,coords = c("lat1","long1")))
DT$p2 <- st_geometry(st_as_sf(DT,coords = c("lat2",&
浏览 0提问于2018-09-28得票数 5
回答已采纳
我正在处理一个空间问题，将sf包与dplyr和purrr结合使用。
我更喜欢在mutate调用中执行空间操作，如下所示：
simple_feature %>%
mutate(geometry_area = map_dbl(geometry, ~ as.double(st_area(.x))))
我喜欢这种方法，它允许我使用%>%和mutate运行一系列空间操作。
我不喜欢这种方法似乎大大增加了sf函数的运行时间(有时令人望而却步)，我希望听到关于如何克服这种速度损失的建议。
这里是一个reprex，它详细地说明了速度损失问题。
请注意：这不是一个很小的例子，需要从ESRI下
浏览 0提问于2018-01-31得票数 6
回答已采纳
1回答
好吧，所以我在前一个问题的基础上。我有一组原始数据，我遵循一个特定的协议(因为IRB的原因我需要遵循这个协议)来随机抵消这些数据。
每个点都有自己的随机偏移量。城市点为0~2 5km，农村为0~5km，随机选取1%的偏移量为0~10 5km。
这是我的代码，到目前为止，我使用的是投影坐标参考系统
displaceData <- function(sfDataset) {
dataset.sf <- sfDataset
# 2. Generate a random direction by generating angle between 0 and 360, a
浏览 4提问于2022-04-03得票数 0
回答已采纳
https://www.r-spatial.org/r/2017/06/22/spatial-index.html上解释说“R树是在第一个参数(x)上构造的，用于匹配二进制函数的第二个参数(y)上的所有几何图形”。但是，最终没有解释可以指导用户选择sf对象以放入第一个参数的标准是什么。这个选择有没有一些经验法则？
浏览 15提问于2021-02-10得票数 2
我在试着检查两个多边形在R中是否相交，在绘图时，它们显然没有。在检查交集时，rgeos::gIntersects()当前返回FALSE，而sf::intersects()返回TRUE。我想这是因为多边形是(1)大的和(2)紧密的，所以当它们在平面上看起来不是相交，而是在一个球面上它们看起来是相交的？
理想情况下，我可以将我的工作流全部保存在sf中--但我想知道是否有一种方法可以在这里使用sf::intersects() (或其他sf函数)返回FALSE？
下面是一个例子：
library(sf)
library(rgeos)
library(leaflet)
library(leaflet.e
浏览 6提问于2022-08-31得票数 5
回答已采纳
1回答
我为一份报告而制作的地图出了点问题。我知道问题，但不知道怎么解决。
我的形状文件中的四个多边形是无效的--它们是环自相交。形状文件是纽约大都会区的。当我在QGIS中打开它时，这四个有问题的多边形是有效的。只有当我使用ggplot或tmap在R中打开它们时，我才有问题。
这是我的密码：
# remove all objects in R and clear R
remove(list=ls())
# Load packages
library(tidyverse) # tidy verse package for data maninpulation
library(sp)
library(
浏览 1提问于2020-05-06得票数 0
1回答
我正在寻找一种快速的方法在我的多边形在R.我想把我的多边形根据它的中心移动到一个新的中心。
library(sf)
library(ggplot2)
theta <- (0:6) * pi / 3
hexagon <- data.frame(
x = sin(theta),
y = cos(theta))
poly <- hexagon %>%
st_as_sf(coords = c("x", "y")) %>%
mutate(geometry = st_combine(geometry))%>%
s
浏览 6提问于2022-09-21得票数 0
回答已采纳
1回答
所以我试图绘制地球上的坐标，并跟踪每个国家有多少坐标。我绘制了地图并进行了精细的坐标，但是当我试图使用交点来计算每个国家(多边形)内有多少坐标时，就会产生误差。我尝试过使用st_make_valid函数来修复地球形状文件，但是它破坏了几何图形。我刚开始使用R，所以任何帮助都是非常感谢的。
我使用了以下代码来绘制地球形状文件和顶部的坐标：
library(tidyverse)
library(sf)
library(rmapshaper)
library(rnaturalearth)
library(rnaturalearthdata)
library(sp)
library(raster)
浏览 16提问于2022-11-22得票数 1
回答已采纳
我正在尝试将下面代码的输出(这是一个的list列表)转换为一个数据格式，其中的列包含每个观察的简单功能集合 -one多边形。
以下是可重复例子的产出和预期产出：
# Reproducible example:
library(tidyverse)
library(sf)
library(magrittr)
# define radius for circle
radius <- 40
r <- units::set_units(radius, units::as_units("nmile"), mode = "standard"
浏览 1提问于2018-02-27得票数 7
回答已采纳
1回答
通过非sf列内连接两个sf对象
、、、、
我试图连接两个sf数据帧使用内连接或左连接。这些数据文件中都有几何图形列。我不断地发现错误：
check_join(x，y)中的错误:y应该是data.frame；对于空间连接，使用st_joinFALSE
可复制的例子如下：
df1 <- data.frame(
var = c("a", "b", "c"),
lon1 = c(20,35,45),
lat1 = c(50,10,15)
) %>% st_as_sf(coords = c("lon1", "lat1"), dim =
浏览 0提问于2018-02-28得票数 5
回答已采纳
使用的图书馆：
library(sp)
library(sf)
library(ggplot2)
library(ggmap)
创建一个名为"coordinate.data“的数据文件，其中以经度和纬度作为列名，以及气象站位置作为行名。
longitude <- c(-73.964482,-73.953678,-73.893522,-73.815856,-74.148499)
latitude <- c(40.767544,40.631762,40.872481,40.734335,40.604014)
coordinate.data <- data.frame(
浏览 3提问于2020-07-02得票数 2
回答已采纳
假设我有下面的数据和代码。代码返回point数据，但我想要一个polygon。
如何执行spatial join，使其返回具有point和polygon属性的polygon？(基本上，数据将根据多边形中的点进行匹配/连接)
代码+数据
library(sf)
library(tidyverse)
# Sample poly
poly = st_read(system.file("shape/nc.shp", package="sf")) # included with sf package
# Sample points
pts = data.frame(
浏览 4提问于2022-05-13得票数 0
回答已采纳
我有一份数据文件，里面有龙/拉特的信息。目的是在rad距离特定点st_p的范围内找到所有的点。实际上，我已经有了R中的代码，但我需要在python上也这样做。
这里我要做的是，我将数据格式转换成sf对象，定义一个缓冲区，即我与缓冲区的交集。
这是R码。我只是不知道在Python中使用什么库才能做到这一点。
within_radius <- function(df, st_p, rad) {
# Transform to an sf object and cahnge from lon/lat to utm
sf_df <- st_transform(st_as_sf(
浏览 1提问于2020-09-25得票数 0
回答已采纳
1回答
网格设计中的随机地理空间采样
、、、、
我试图创建一个空间调查设计，首先在一个研究地点设置一个网格，然后在每个网格单元中随机抽样。下面是一个虚构的例子，它在蒙大拿州创建了一个六边形网格。问题是st_sample在整个形状上提供随机点，而不是在网格中。我在st_sample中尝试了另一个st_sample，它们不起作用。我怀疑st_sample不认为我的网格是多个多边形？
有什么想法，或者使用这些包/函数或其他在R？
States <- raster::getData("GADM", country = "United States", level = 1)
Montana <- S
浏览 15提问于2022-07-23得票数 2
3回答
我想知道是否有一种简单的方法来编写包含坐标的点sf对象(sf R包)的CSV。
你可以使用st_write(input, "output.csv")，它会写一个没有坐标的CSV。我写一个带有坐标的文件的技巧是：
coords <- st_coordinates(input)
input_dat <- input %>% st_set_geometry(., NULL)
input_dat <- cbind(input_dat, coords)
但似乎必须有一种更简单的方法。
按照要求，下面是上面代码的设置：
input <- data.frame(
浏览 4提问于2017-08-05得票数 6
1回答
使用R包sf，我试图确定某些点是否出现在一个shapefile的范围内(在本例中，Hawai‘s，EEZ)。有关的shapefile可以找到。不幸的是，这个区域的边界跨度为+/-180度，我认为这就是困扰我的原因。(我在sf网站上读到了一些关于球面几何学的新版本，但是我还没有能够安装那个版本。我认为我正在处理的多边形是足够“平坦”，以避免任何这些问题)。部分问题似乎是我的shapefile包含多个几何图形，这些几何图形被日期线分解了，但我不知道如何组合它们。
如何使用sf判断某些点是否在shapefile中某个对象的边界内(恰好跨越日期线)？
我尝试过st_shift_longitude的各种
浏览 4提问于2021-03-25得票数 1
回答已采纳
我正在尝试建立一张德国联邦州的地图，带有特殊的边界(投票区)：
install.packages("OpenStreetMap")
install.packages("sf")
install.packages("osmdata")
install.packages("tmap")
library(OpenStreetMap)
library(sf)
library(osmdata)
library(tmap)
## I use this because the other overpass server didnt wor
浏览 7提问于2021-05-14得票数 0
我希望将sf对象转换为dataframe，并将其还原到原始状态。但是，当我进行st_as_text(st_sfc(stands_sel$geometry))转换时，很难再次检索它。在我的例子中：
library(sf)
# get AOI in shapefile
download.file(
"https://github.com/Leprechault/trash/raw/main/sel_stands_CMPC.zip",
zip_path <- tempfile(fileext = ".zip")
)
unzip(zip_path, ex
浏览 2提问于2021-10-12得票数 1
回答已采纳
我试图使用R包tmap绘制一个SpatialLinesDataFrame。然而，我一直收到这样的错误消息：
CPL_geos_is_empty中的错误(st_geometry(X))：计算错误: IllegalArgumentException:点阵必须包含0或>1个元素。
我第一次尝试用tmap绘制"sp“对象。然后，我将sp对象转换为"sf“对象，仍然得到相同的错误消息。我做了一点谷歌搜索，我认为这可能与SpatialLinesDataFrame中的线路没有完全连接的事实有关。
所以我给出了一组连在一起的线，它用tmap绘制得很好。
我把这里使用的形状文件放在
浏览 1提问于2019-08-28得票数 3
2回答
R中的Graphml到Sf对象
、、、、
我正在尝试将一个graphml文件的节点转换为R中的sf对象，直到我创建了sfc对象。一旦我添加了属性，R就不再能够对其进行绘图或执行任何操作。 library(igraph)
library(sf)
#read the graph
gi <- read.graph(
"https://jospueyo.github.io/root/girona.graphml",
format = "graphml")
#Extract coordinates and create id vector
x <- as.numeric(V(gi)$x)
浏览 21提问于2020-07-24得票数 0
我正在使用tidycensus软件包来提取人口普查数据和几何图形。我希望能够计算人口密度，并使结果与我在上看到的结果相匹配。我注意到从tidycenus返回的地理变量与我使用sf包sf::st_area()函数计算的结果不同。
library(tidyverse)
library(tidycensus)
census_api_key("my_api_key")
library(sf)
options(tigris_use_cache = TRUE)
pop_texas <-
get_acs(geography = 'state',
var
浏览 11提问于2022-06-18得票数 0
回答已采纳
下面是我在C++中使用SFML库的示例程序。我想要创建一个自定义函数'draw_func‘，它在提供的坐标处绘制一些东西(如本例中的矩形)。我将返回变量的类型设置为sfml对象矩形(这是我返回的和绘制的)，但是屏幕是黑色的。
#include <iostream>
#include <math.h>
#include "SFML/OpenGL.hpp"
#include <SFML/Graphics.hpp>
sf::RectangleShape draw_func(int x, int y)
{
sf::Rectangle
浏览 0提问于2018-03-26得票数 1
回答已采纳
使用osmdata绘制从osm下载的道路数据时，在geom_sf中使用较大的size值时，生成的打印中会有间隙(请参见下图)。
这里是一个使用伦敦西南部的一段道路的可重复使用的例子。如何在绘图时消除线条上的白色间隙？
library(tidyverse)
library(sf)
library(osmdata)
# define bounding box for osm data
my_bbox <-
matrix(c(-0.2605616, -0.2605616,
-0.2004485, -0.2004485,
-0.2605616,
浏览 1提问于2018-08-13得票数 3
回答已采纳
我有两个sf对象，我希望将它们并排绘制为ppp。根据的问题，绘图一开始运行良好，绘图也很好，但随后我单击了“清除所有绘图”，生成了两个并排的绘图。在此之后，坐标不再显示在任一轴上，其中一个绘图显示字符，如下所示。我甚至重新启动了R，但结果仍然相同。下面是代码，示例数据可以从下载。我不知道清除绘图窗口后出了什么问题。如何解决此问题？
sf = st_read("~", "sf")
# Regrouping data
sf = sf::st_as_sf(sf, coords = c("Long", "Lat")) %>%
浏览 21提问于2021-03-04得票数 0
1回答
我有一个包含多个序列(T1、T2、T3 )的sf多边形数据帧，它们都在相同的尺度上:它们是不同时间点的观测值。我可以用来描述T1
ggplot(map)+geom_sf(aes(fill=T1))
我想要做的是将所有这三个(T1，T2和T3)作为面(单独的地图)绘制在同一个图形上。我确信有办法做到这一点，但我找不到。有人能告诉我是怎么做到的吗？谢谢!
添加:关于此问题的两个附加说明。首先，上面描述的数据结构可以使用spplot绘制，其中T是spplot的zcol参数的参数。因此，在这方面，我的问题相当于询问如何将spplot结构转换为geom_sf可用的结构。
其次，假设我使用sf读入一个sh
浏览 0提问于2021-02-06得票数 0
1回答
将json对象转换为空间线数据帧
、、、、
我从EnergyData下载了一个伊拉克电网的geojson文件。您可以自己下载文件，并访问网页。
我试图使用geojsonio::geojson_read将文件读入R中，但是它抛出了错误。另一位用户帮助指出，能源数据网站提供的文件格式不正确，不符合geojson格式，因此文件名具有误导性。所以他们推荐我使用jsonlite::read_json。虽然该函数确实将文件读取到R中，但它不是空间对象。这是一个多层次的名单。我也尝试过sf::st_read，但是对象是一个凌乱的数据框架。我需要我的最后一个对象是一个空间线数据帧，这样我就可以将它与另一个空间行对象合并并裁剪，并写出生成的shapefi
浏览 9提问于2021-12-29得票数 1
回答已采纳
当我运行exact_extract(rasterlayer, sf_object)时，我得到了以下错误： in CPP_exact_extract(x, weights, wkb, default_value, default_weight, : Unsupported geometry type. 我不确定发生了什么，因为我的其他提取几乎相同的sf对象的工作。我想通了，所以在这里记录我的答案。
浏览 17提问于2021-09-11得票数 1
回答已采纳
我目前正在尝试添加标签到美国县的地图，显示县的名称。我使用geom_sf可视化了地图，一切都很好。我指定了一些其他的东西，比如填充、美观等等，它会像我想要的那样绘制出来。但是如果我尝试通过geom_sf_label添加标签，我总是得到以下错误：stat_sf_coordinates requires the following missing aesthetics: geometry。
以下是工作的代码：
df %>%
group_by(county_fips)%>%
count()%>%
full_join(geo_data, by = c("count
浏览 22提问于2019-11-28得票数 2
回答已采纳
我是C++的新手。我正在用SFML和Box2D开发一个中断克隆，在编译时我得到了这个错误。错误的详细信息：
c:\program files (x86)\visual studio express 2013\vc\include\xutility(2420): error C2280: 'Tile &Tile::operator =(const Tile &)' : attempting to reference a deleted function
c:\users\harry\documents\visual studio 2013\projects\cpp
浏览 7提问于2014-07-24得票数 2
回答已采纳
我有五个经度和纬度，形成这样的形状。
df <- c(order=1:5,
lon=c(119.4,119.4,119.4,119.5,119.5),
lat=c(-5.192,-5.192,-5.187,-5.187,-5.191))
如何使用这样的sf包轻松地将它们转换为sf多边形数据帧？
## Simple feature collection with 1 feature and 0 fields
## geometry type: POLYGON
## dimension: XY
## bbox: xmin: 1
浏览 0提问于2018-01-22得票数 17
回答已采纳
2回答
sf对象的df <- df[complete.cases(df), ]等价物是什么？
我知道sf对象的行为很像带有“粘性”geometry列的dataframes，但是运行上面的代码会返回以下错误：
df <- df[complete.cases(df), ]
Error in complete.cases(df) : invalid 'type' (list) of argument
这可能是因为sf对象将其geometry变量存储为列表。非常感谢您的意见！
浏览 13提问于2018-09-05得票数 1
回答已采纳
1回答
R中的全球地图
、、、
我试图使用R中的地图包创建一个全球地图，我试图使用来自wrld_simpl的数据。我期待的是一个地图，其中的人口("Pop2005“变量)是绘制的。我开发的代码如下：
data(wrld_simpl)
world<-wrld_simpl
world_sf = st_as_sf(world)
world_sf_proj = st_transform(world_sf, crs = 3785)
world_cartogram <- cartogram_cont(world_sf_proj, "POP2005")
plot(world_cartogram)
尽
浏览 8提问于2022-02-02得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
