有没有办法添加丢失的CRS？我已经得到了一些旧的Esri shapefile，其中很大一部分有CRS: NA，我希望能够批量设置丢失的CRS。我可以在导入QGIS时手动设置它，但似乎不能使用R设置它。我尝试了st_transform()，但如果一开始没有设置CRS，它似乎不起作用。
浏览 18提问于2020-04-21得票数 0
回答已采纳
背景plot(us_shps, main = "Original")问题[41] "40" "41" "42" "43" "
浏览 1提问于2017-09-05得票数 4
回答已采纳
25.7298909 8281.33216 46.6298909 8298.9321 V1是x，V2是y，我想知道是否有办法将其转换为sf以下是前两个多边形的可复制数据集： > dput(head(newdata, 38))
structure(list(.id = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
浏览 21提问于2021-08-16得票数 1
回答已采纳
我目前正在处理R中的点的sf数据集，我想将该数据集转换为SpatialPointsDataFrame，以便进行一些下游分析。这看起来很简单，所以我使用了as_Spatial()函数，但它抛出了一个我以前从未见过的错误： Error in sp::CRS(SRS_string = from$wkt) : unused argument(SRS_string = from$wkt) 我尝试转换的对象对于我添加到问题中来说
浏览 24提问于2020-10-06得票数 1
回答已采纳
1回答
给定一组具有全球分布的坐标，我想在世界地图上绘制一个多边形(最好是一个矩形)，然后提取出这些图中的所有坐标。多边形的位置将根据点密度手动放置，但它们的大小必须保持不变。在这里，我提供了一个在欧洲随机分布的点的最小可重现性示例。我还添加了一个定向图像，希望它有助于理解所需的输出。首先，我想在地图上添加多边形，然后提取该区域内的所有点。谢谢你提前提供任何可能的帮助。lat = coord$lat)
coordinate
浏览 6提问于2022-01-29得票数 1
回答已采纳
2回答
我想在美国或美国各地区建立网格(在x坐标和y坐标的数据帧的意义上)，抛出在美国边界以外的原始网格中的任何点。我有一些代码似乎可以工作，但是当我将网格增量缩小到1km (1e3)左右时，它会非常慢。也许有一种方法可以在没有lapply或循环的情况下构建我需要的简单特性集合，或者这可以用MULTIPOINT来完成，而不是使用POINT的简单特性集合。library(sf)
crs.us.atlas = 2163 # h
浏览 2提问于2021-01-12得票数 1
回答已采纳
所以我有两个SF对象，一个多边形和一个多边形类，代表一个城市的居民区和子居民区。两者都来自ArcGIS，并且覆盖了完全相同的区域。然而，当我绘制这两个对象时，子邻域在我的单叶图中稍微对齐了一点。原来的预测是：Simple feature collection with 92 features and 19 fields
geometry type:5.38763888888889 +k=0.999
浏览 1提问于2018-02-13得票数 1
这就是我如何将它转换为spatialPolygons (包sp)的方法。7550903, 7550903, 7550503, 7550503, 7550903, 7550803, 7550803, 7550403, 7550403, 7550803))
如何直接从my.d
浏览 3提问于2020-04-14得票数 8
回答已采纳
2回答
我有一个栅格，我需要根据sf对象中的多边形属性进行编辑。这两个对象具有相同的范围。填充了9的栅格 r <- raster( nrow=10, ncol=10, xmn=0, xmx=10, ymn=0, ymx=10 )classobject# Set same
浏览 20提问于2020-01-11得票数 2
回答已采纳
我最近开始使用机载LiDAR数据，该数据具有一个由水平(投影)和垂直分量组成的复合坐标参考系统。下面显示了一个示例，它使用代码从WKT描述创建复合CRS对象。我从LiDAR点云派生出各种光栅层，我只想将复合CRS的水平组件分配给每一个(示例中的EPSG:7856)。是否有人知道现有的包函数将可靠地提取水平PROJCRS组件，即允许各种新旧CRS定义？更新2021-11-01:调整了WKT字符串的原始示例，以提供在R<
浏览 2提问于2021-10-28得票数 0
回答已采纳
我在一个目录中有1000多个shape文件，我只想选择其中10个我已经知道其名称的文件，如下所示：
15TVN44102267_Polygons.shp, 15TVN44102275_Polygons.shp我将感谢来自社区的任何帮助。
浏览 1提问于2019-06-15得票数 1
2回答
我有一个关于纽约市黄色出租车服务区的shapefile : taxi_zones.shp。下面是我使用的R代码： library(sf) "x" = c(-73.97817,-74.00668,0,500),
pnts_sf <
浏览 31提问于2020-10-09得票数 1
回答已采纳
我有一组纬度/经度点，试图转换成R中的一个简单特征(sf)，我的目标是用从urbnmapr库检索到的多边形在美国地图上绘制这些位置。
用我们的地理参考绘制，如代码所示，结果显示所有的点。当使用geom_sf()绘制这些点时，它们最终会到达南达科他州。尽管我认为st_as_sf()函数的正确使用是正确的，但纬度/经度点似乎没有被转换成正确的坐标参考系统。turbi
浏览 4提问于2020-08-29得票数 0
回答已采纳
2回答
library(randomcoloR)
my.shp$ID<- 1:nrow(my.shp) main = list(label = "TEST"))
我的</em
浏览 0提问于2019-01-20得票数 3
回答已采纳
3回答
有没有一种简单的方法可以将一个简单的特征(sf对象)的边界转换到另一个CRS？ 函数st_bbox()的结果属于bbox类。无法使用st_transform()将其转换为另一个CRS。我使用的是基于EPSG: 28992的计算边界框： sf::st_bbox(xmin = 187470, xmax =194587,
ymin = 409753, ymax
浏览 13提问于2019-03-08得票数 4
回答已采纳
我试图弄清楚为什么我不能让两个sf对象连接在一起。邻域赏金对象。为了获得bu_code变量，我尝试使用与邻域代码sf对象相同的crs创建一个sf对象。= st_crs(neigbourhood_boundries_sf)) Simple feature collection with 80840 feature
浏览 1提问于2021-10-06得票数 0
回答已采纳
1回答
我有一个名为Lon和Lat的dataframe a.fil。name1 43.37390 132.9703 name2 43.35311 132.7493我创建对象点point2.sfg <- st_point(unname(unlist(a.fil[2, 2:3])))i <- 1
浏览 1提问于2020-06-29得票数 1
1回答
我有一个光栅和一个形状文件：library(sf)mtq <- st_read(system.file("gpkg/mtq.gpkg", package="cartography"), quiet = TRUE)
我想将光栅r与shapefile mtq相交，并将相应的像素作为NA (
浏览 1提问于2022-03-02得票数 0
回答已采纳
我有两个带有sf几何的数据框，其中一个有两个坐标，另一个有三个坐标，即一个坐标基于c(1,2)，另一个基于c(3,4,0)。我想绑定这两个表，但失败了，因为R“期望三个坐标”。我想知道在sf中将几何图形从两个坐标转换为三个坐标的代码是什么？换句话说，我希望将第一个坐标从c(1，2)更改为c(1，2，0)。 非常感谢您提前给予的帮助:)
浏览 18提问于2021-08-10得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
