我试图得到两个形状文件的交汇点(属于某些大都市区范围内的人口普查区)。我能够成功地获得相交特性,但是当我试图将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<-
好吧,所以我在前一个问题的基础上。我有一组原始数据,我遵循一个特定的协议(因为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
我正在运行一个django应用程序,我正在为我的db使用PostGis扩展。当我发送坐标时,我试图更好地理解引擎盖下发生了什么,特别是因为我使用的是不同的坐标系统,它们会转换成不同的SRID。我的问题有三个:
是django/postgis在DB中创建点或多边形时处理转换。
假设我有这样一个模型(请注意,我正在设置标准SRID=4326):
class MyModel(models.Model):
name = models.CharField(
max_length=120,
)
point = models.PointField(
感谢您抽出时间来阅读。 我不知道为什么这段代码不能工作。总而言之(更改变量名),我有圣路易斯餐馆的位置(坐标、经度和经度),以及便利店的位置(也是lat )。我想看看便利店离餐厅有多远。 #Bring in data
# MO county data
county.sf <- get_acs(state = "MO",
county = c("St. Louis County", "St. Louis City"),
geography = "tract",
我有一组点和,我想从几个大型光栅中提取值,作为围绕这些点的缓冲区。栅格太大,无法保存在内存中(> 1e10细胞)。我在下面说明我目前的方法,但如果有更快的方法,我会感兴趣的。
library(maps)
library(sf)
library(raster)
library(dplyr)
library(parallel)
# sf object with polygones for which we want values
crs <- "+proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +n
我为一份报告而制作的地图出了点问题。我知道问题,但不知道怎么解决。
我的形状文件中的四个多边形是无效的--它们是环自相交。形状文件是纽约大都会区的。当我在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(
我把R中的sf搞得一团糟;我正试图从一个大型的坐标数据框中创建一个特性集合。 我的数据是x和y坐标的data.frame,ID决定了它们的分组(即所有相同的ID都是一个多边形)。我想用这个来创建一个sf集合; # create some dummy data, coordinates for 3 squares (v. similar to provided data)
df <- data.frame(x = c(0,0,1,1,2,2,3,3,4,4,5,5), y = c(1,2,2,1,1,2,2,1,1,2,2,1), ID = rep(1:3, each = 4))
#