在R中组合点shapefile和道路shapefile
在R中，要组合点shapefile和道路shapefile，可以使用sf包来处理空间数据。sf包提供了一套功能强大的工具，用于读取、处理和分析空间数据。
首先，我们需要加载sf包，并读取点shapefile和道路shapefile。假设点shapefile的文件名为points.shp，道路shapefile的文件名为roads.shp，可以使用以下代码读取这两个shapefile：
library(sf)
# 读取点shapefile
points <- st_read("points.shp")
# 读取道路shapefile
roads <- st_read("roads.shp")接下来，我们可以使用st_join函数将点shapefile和道路shapefile进行组合。st_join函数可以根据空间关系将两个shapefile中的要素进行匹配。假设点shapefile中有一个属性字段ID，道路shapefile中有一个属性字段Road_ID，我们可以使用以下代码进行组合：
# 根据空间关系将点shapefile和道路shapefile进行组合
combined <- st_join(points, roads, join = st_intersects)上述代码将返回一个新的sf对象combined，其中包含了组合后的点和道路要素。
在云计算领域中，可以将上述操作应用于各种场景，例如地理信息系统（GIS）、交通规划、位置分析等。腾讯云提供了一系列与地理信息相关的产品和服务，例如地理位置服务（https://cloud.tencent.com/product/lbs）、地理信息系统（https://cloud.tencent.com/product/gis）等，可以根据具体需求选择适合的产品。
请注意，以上答案仅供参考，具体的实现方法和腾讯云产品选择应根据实际情况进行调整。
相关·内容
我想知道如何组合两个shapefile文件，一个是点shapefile，另一个是R软件中的roads shapefile。在附件中，我留下了一张我在Arcgis软件中结合了这两个shapefile的图像。 为了进行测试，这两个shapefile都可以从以下网站下载：https://github.com/JovaniSouza/JovaniSouza5/blob/master/Example.zip 非常感谢！ ?
浏览 10提问于2020-06-30得票数 1
回答已采纳
我在R中学习spatstat，作为第一步，我正在转换我的shapefile以适合spatstat算法。当我使用as.linnet转换成线性网络(下面代码中的最后一行)时，我会得到一个警告network not connected
roads<-readShapeSpatial("OnlyHigwayroads.shp")
spatstat.roads<-as.psp(roads)
final_roads<-as.linnet(spatstat.roads)
Warning message:
Network is not connected
我不明白这是什么意
浏览 0提问于2016-10-28得票数 0
回答已采纳
1回答
我有推文存储在MySQL数据库中，我已经下载了墨尔本公路网的形状文件。现在，我只想绘制那些位于shapefile中的tweet，即在道路上。
我该怎么做呢？
浏览 0提问于2016-10-08得票数 0
我正在从事一个大学项目，它把数据科学和地理信息系统结合起来。我们需要找到一个开源的解决方案，能够从一个庞大的GPS坐标数据集中获得更多的信息。显然，我不能使用任何有每日请求限制的API。
数据
在这里，您可以找到教授提供给我们的数据集的示例：
longitude <- c(10.86361, 10.96062, 10.93032, 10.93103, 10.93212)
latitude <- c(44.53355, 44.63234, 44.63470, 44.63634, 44.64559)
longlat <- data.frame(longitude,
浏览 0提问于2017-12-06得票数 1
回答已采纳
2回答
R中的最短距离
、、
我想知道如何为下面的代码计算两个属性(点)之间的最短距离。有两个shapefile文件，一个是points shapefile，另一个是roads shapefile。 为了进行测试，这两个shapefile都可以从以下网站下载：https://github.com/JovaniSouza/JovaniSouza5/blob/master/Example.zip library(sf)
roads <- st_read('Roads/Roads.shp')
pts <- st_read('Points/Points.shp') %>%
浏览 67提问于2020-07-01得票数 1
回答已采纳
1回答
我尝试用(Lat，Lon)和每个位置的值来绘制数据。我需要做的是在shapefile上绘制数据，用不同的颜色表示每个值，并平滑这些值，这样地图中就不会出现间隙。
数据结构如下：(经度，纬度，值)，你能帮我吗？
我知道如何在不同颜色的shapefile上绘制数据，但问题是如何平滑这些值。
问题的数据可在以下链接中找到：
其中：
1- (经度、纬度、值)数据在csv文件"cluster__groups“中
2- shapefile文件位于"shapefile_question“文件夹中
我需要在区域中显示它们。输出应类似于下图：
谢谢你的帮助。我试图解决并寻找解决方案，
浏览 7提问于2017-11-20得票数 2
2回答
我有5个加州边界县的shapefile文件，我正在尝试将它们合并成一个大的shapefile文件，并保持这些县的原始边界。然后我想要绘制一些点，比如，这样地图只需要是一个简单的县边界地图。有没有办法在R中做到这一点，或者我必须使用像ArcGIS这样的东西？
有大量关于在shapefile中合并的文档，比如，其中两个shapefile具有相同的坐标。但是没有关于只部分共享坐标的shapefile的文档。
浏览 2提问于2010-06-03得票数 2
回答已采纳
我正在尝试寻找一组位置坐标和一组线(道路或河流)之间的正交距离。这组点的形式是纬度/经度对，而线条是shapefile (.shp)。使用maptools或PBSmapping在地图上绘制它们都不是问题。但我的基本问题是找出一个人从一个地方到一条路或一条河所需的最小距离。有没有办法在R中做到这一点？
浏览 1提问于2012-03-20得票数 13
回答已采纳
使用的图书馆：
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
回答已采纳
下面是我的代码，用来检查多边形中是否存在一个点，但是结果是假的，你能帮我吗？
# Importing Libraries
import shapefile
from shapely.geometry import shape, Point
# Opening the Vector Map
shp_path = "ADMINISTRASIDESA.shp"
# read your shapefile
r = shapefile.Reader(shp_path)
# get the shapes
shapes = r.shapes()
# build a shapel
浏览 5提问于2021-06-10得票数 0
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
我已经用Java语言编写了一个内核密度估计器，它以ESRI shapefile的形式接受输入，并输出估计表面的GeoTIFF图像。为了测试这个模块，我需要一个示例shapefile，不管出于什么原因，我被告知要从R中包含的样本数据中检索一个。问题是，没有一个样本数据是shapefile……
因此，我尝试使用shapefiles包的convert.to.shapefile(4)函数将R中的spatstat包中包含的bei数据集转换为shapefile。不幸的是，事实证明这比我想象的要难。有没有人有这样做的经验？如果你能在这里帮我一把，我将不胜感激。
谢谢你，莱恩
参考文献：，
浏览 0提问于2011-04-27得票数 6
回答已采纳
2回答
如何修复这个线型文件？
、、、
我正在从自动生成的点shapefile创建一条线。当我第一次在ArcGIS中创建这条线时，我得到了一条这样的线，因为点没有按顺序排列： ? 之后，我根据它的位置对点进行了排序，得到了一条如下所示的线： ? 但是无法创建这样的行： ? 请给我在ArcGIS或R编程中修复这个问题的任何解决方案。如果您需要shapefile，我可以为您提供。
浏览 30提问于2020-06-27得票数 0
我希望写一个R程序，读取一个数据帧的lat/long点和一个由13个多边形组成的shapefile，然后识别每个lat/long点位于哪个多边形内。是否有一个R包可以让我完成这种空间连接？我一直在使用rgdal来读取多边形shapefile，但是我没有看到在我的dataframe中使用lat/long列进行空间连接的选项。谢谢!
浏览 4提问于2016-04-25得票数 2
1回答
首先，我要说的是，我已经彻底地查看了处理readOGR错误的其他线程，并为我自己的利益多次尝试实现它们。
我遇到的问题类似于前面提到的问题:也就是说，我不能通过readOGR导入一个shapefile。
我的工作目录是：
("~/Desktop/R_Thesis_Project")
请不要让我意识到这不是完整的工作目录。我也输入了'dns‘，这会产生同样的结果。
shapefile (包括.shp、.dbf、.xml、.shx和.prj )位于：
("~/Desktop/R_Thesis_Project/Shapefile")
请注意，该文件夹没有
浏览 2提问于2016-04-20得票数 4
回答已采纳
我对R相对来说是个新手，所以如果我很难表达我想要做的事情，我很抱歉。我有一个长格式的“空间”面板数据集和一个shapefile。长表单表是一个data.frame，它包括一个日期列(已使用‘as.date’转换为日期)和一个ID列，该ID列与用于识别不同多边形的shapefile中的ID列相同(因此，我的长表单数据集没有长的lat值，只是对应于shapefile中的多边形特性的ID字段)。我想从这两个对象( shapefile和long form dataset)中构造ST类的时空对象。为了做到这一点，我尝试使用stcontruct()和STFDF()，但绝对没有运气。stcontruct()
浏览 3提问于2014-04-22得票数 0
1回答
我正在尝试转换一个具有较小值的多边形shapefile。列propEmp中的值范围从0.000002到0.119419。这是我的尝试：
# Load shapefile
emp_planejado <- shapefile("./planejado/7_wgs/emp_planejado.shp")
# Load raster model
r_bioma<- raster("./_GRID/grid_caatinga_disol_64bit.tif")
# List values from tipologia field
tipologia<
浏览 2提问于2015-09-18得票数 0
在Fiona1.5.0上(我很困惑为什么不同的文件(如.dbf和.gdb)没有打印我的“！”(这是我希望在文件不是.shp的任何时候)退出前的警告。
import fiona
import sys
def process_file(self, in_file, repair_file):
with fiona.open(in_file, 'r', encoding='utf-8') as input:
# check that the file type is a shapefile
if
浏览 0提问于2015-04-23得票数 2
回答已采纳
如何合并具有多个多边形的shapefile中的多边形特征？ rbind和union只是合并shapefile功能的行，它们实际上并不合并多边形本身。 以下示例中的预期结果： 如何获取具有重复ID的以下shapefile _ 2在sptemp中合并为单个多边形？ 以下埃塞俄比亚GADM级别2的示例具有shapefile ID的前两行 _ 2列重复(value=1)。我希望sptemp具有79个功能，将前两行与ID重复的行组合在一起 _ 2. sptemp图 1、 将调出当前sptemp的位置 1、 和sptemp2 2、 没有重复的边界，也就是说多边形也被合并了。 示例代码： 下载、解压缩并加
浏览 40提问于2018-03-19得票数 10
回答已采纳
1回答
在R中绘制地图点
、、、
我已经用maptools成功地将NYC PUMA区域的shapefile加载到R中，我想要在它上面绘制55个点，我在另一个文件中有，如下所示：
X Y pumace10 events_2008 events_2009
-73.9092456917 40.8916125162 3701 2 0
-73.8617096298 40.8899373255 3702 0 0
-73.8010284966 40.8460832277 3703 1 1
但是，这些点将不会打印。
首先，我这样做是为了绘制shapefile：
plot(nycs)
浏览 13提问于2017-02-28得票数 1
我需要导入自定义道路到谷歌地图，(即:新的道路，可能还看不到谷歌地图)，我希望他们被包括在方向服务API中。我首先考虑的方法是使用谷歌地图引擎导入包含我所需道路数据的shapefile，并在此基础上创建地图。然后，我将对该自定义地图中的数据使用方向服务。我非常担心这是否会起作用，据我所知，目前还没有Google Maps Engine的试用版可以用来测试它。如果它不能工作，我真的不想买它。
有没有人知道这种方式是否可行，或者谷歌是否提供了类似的服务，可以满足我的需求？
如果我想这样做，我假设我将需要实现我自己的directions服务，但我想看看是否已经有一个服务可以替代它。
谢谢!
编辑:试
浏览 6提问于2012-08-29得票数 0
2回答
有没有允许交互式显示GIS shapefile的软件包？我正在寻找创建一个简单的图形用户界面，显示简单的shapefile(海岸线等)，但似乎找不到从哪里开始。我最初的任务是用R语言来做这件事，但我想尝试使用python和Qt。
浏览 3提问于2012-08-09得票数 4
回答已采纳
我目前正在探索三个shapefile，每个都有点数据，并且都被限制在同一个窗口中。我也有他们在ppp格式，我已经用来创建内核密度贴图。
plot(density.ppp(smktppp, 0.5, edge=T), main="Supermarket Density")
plot(density.ppp(tptppp, 0.5, edge=T), main="Transport Density")
plot(density.ppp(farmppp, 0.5, edge=T), main="Urban Farm Density")
我想叠加这些图，
浏览 0提问于2017-02-26得票数 1
我目前正在阅读Joel Lawhead的一本Python地理空间分析教程，在尝试编辑shapefile时遇到错误。
我正在使用的shapefile是一个。我在python3的Jupyter Notebook上运行我的代码。
这是我的代码。我只是简单地读入一个shapefile作为Reader对象r，并创建一个与r具有相同shapetype的新写入器对象w。然后，我尝试将记录从r附加到w。
import shapefile
r = shapefile.Reader("NYC_MUSEUMS_GEO")
w = shapefile.Writer(r.shapeType)
w.fie
浏览 24提问于2017-12-23得票数 2
我有一个用R写的代码，用来检测村庄之间的水量(或shapefile多边形中的点)。我正在创建一个超过2379个村庄(村庄有数字名称)的无向图，所以这将带来28亿个组合。这张图只是一个数据框架。我有每个村庄的坐标：
VillageName1, VillageName2
Village1, Village2
Village1, Village3
Village1, Village4
...
Village2379, Village2378
对于每种组合，我都有一个函数，它将以100米的间隔创建从Village1到Village2的点。因此，如果Village1和Village2相距1公里，将创建8
浏览 0提问于2017-09-26得票数 0
因此，我有我的栅格文件
r <- raster('ras')
和一个shapefile
abys <- readShapeSpatial('abys')
我用以下方法计算了shapefile定义的平均值：
r.vals<- extract(r,abys)
r.mean <- lapply(r.vals,FUN=mean)
然而，当我使用几个shapefile时，当我返回输出时，我得到了多个结果，例如：
[[1]]
[1] 9321
[[2]]
[1] 6616
[[3]]
[1] 8348
它应该只返回一个，这是我通常得到的。这是因为
浏览 0提问于2012-10-20得票数 1
回答已采纳
该问题存在于将.csv文件中呈现的列作为新字段附加到现有shapefile中。因此，我已经使用Python和模块pyshp和csv，首先，复制原始shapefile的内容(几何图形和记录)，然后，在这个副本中创建新字段，并在相应的.csv行中迭代，以便在其中插入：
import os, sys
import shapefile, csv
from os.path import basename
filename_full = sys.argv[1]
output_full = sys.argv[2]
name, file_extension = os.path.splitext(file
浏览 0提问于2017-09-12得票数 0
我正在尝试获取我创建的一个县内的分区的shapefile，并将其与县评估区给我的另一个shapefile (地块数据)对齐。当我试图让它们排成一行时，我的街道shapefile就不会与其他所有东西对齐。它们都在相同的坐标系上，我不想为细分重新创建shapefile。有什么想法吗？
浏览 1提问于2008-10-24得票数 0
回答已采纳
我在readShapePoly包中使用maptools读取了一个shapefile，但是不能用readOGR读取相同的文件。我希望有人能帮助我阅读readOGR的shapefile。
我从这里下载了文件orcounty.shp：
我还下载了相关文件：orcounty.shx、orcounty.sbx、orcounty.sbn和orcounty.dbf，并将所有五个文件放在文件夹中：c:/users/mark w miller/gis_in_R/shapefile_example/
下面的代码读取shapefile并显示一些属性：
library(maptools)
setwd('c:
浏览 4提问于2014-03-09得票数 14
回答已采纳
在执行下面的代码行之后，尝试读取shapefile (.shp)并接收到ImportError。在执行下面的代码之前，我甚至还安装了GDAL、geopandas、fiona和shapely。
shapefile = gpd.read_file(r'C:\Users\....\Data Analysis\shapefiles\louisville.shp')
收到的错误如下，
---------------------------------------------------------------------------
ImportError
浏览 6提问于2021-10-11得票数 11
我正在从事一个项目，其中涉及GPS坐标从近海位置。我想为我的每个点测量离海岸的距离。我已经在QGIS中创建了一个有问题的海岸线的shapefile，并使用st_read()函数(本例中名为"biminishore“)成功地将其导入到R中。使用以下代码，我可以在ggplot2中绘制我的shapefile： bplot = ggplot() +
geom_sf(data = biminishore, size = 0.1, color = "black", fill = "green1") +
ggtitle("Bimini, The
浏览 48提问于2021-04-10得票数 0
1回答
我在栅格化shapefile以在0.5*0.5网格上生成点时遇到了问题。shapefile表示全球珊瑚礁的风险级别(低-0，中-100，高-1000，极高-1500)到综合威胁的分类。
我从另一个运行良好的示例中提取了代码，但是当我对我的数据尝试它时，我从plot函数中得不到任何东西。有关shapefile和我的代码的链接，请参见以下内容：
# Read shapefile into R
library(rgdal)
library(raster)
int.threat.2030 <- readOGR(dsn = "Global_Threats/Integrated_
浏览 7提问于2016-01-30得票数 3
回答已采纳
2回答
简化R中的二元栅格化
、、、、
我有一些非常小的国家级多边形和点形状文件，我想在R中进行栅格化，最后的产品应该是一个全局二进制光栅(指示网格单元中心是否被多边形/点覆盖在单元格内)。我的方法是遍历shapefile并对每个shapefile执行以下操作：
# load shapefile
shp = sf::read_sf(shapefile_path)
# create a global raster template with resolution 0.0083
ext = extent(-180.0042, 180.0042, -65.00417, 75.00417)
浏览 2提问于2020-11-05得票数 0
回答已采纳
1回答
我无法在R中计算具有多部分多边形的导入shapefile的面积(一个特征包含两个单独的多边形)。我注意到ArcMap给了我一个与raster::area不同的shapefile面积值。为了找出哪个程序给了我正确的面积，我将shapefile分成了几个部分，并重新计算了两个单独多边形的面积：
library(raster)
> single_part <- shapefile("../Desktop/test/test_sp.shp")
> area(single_part)
[1] 575924.0 433409.8
> sum(area(single
浏览 2提问于2018-09-18得票数 0
我习惯于在ArcMap中工作，正在尝试将分析转换为R。
我有4个shapefile，代表4个评估参数。每个shapefile都有一个名为Parameter_Score的唯一列。每个shapefile也有一个列，其中包含与该像素相关联的Parcel_ID。每个shapefile具有不同的行数，因为每个参数对某些地块的影响不同，因此列出了不同的地块。
我想创建一个shapefile，其中包含Parcel_ID信息以及每个参数的唯一分数。我的问题是shapefile都有不同数量的元素。我想根据Parcel_ID将它们合并在一起，并让值0填充未包含在某个参数中的地块的分数。
例如，我想使用这两个sha
浏览 4提问于2017-03-18得票数 0
我试图使用来自zonal_stats包的函数rasterstats从.shp文件中的每个形状的.tif文件中获取光栅统计信息。我设法在QGIS中做到了这一点，没有任何问题，但我必须对200多个文件进行同样的操作，这将需要很长时间，所以我正在尝试Python方法。文件和复制代码都在我的中。
我的剧本是：
import rasterio
import geopandas as gpd
import numpy as np
from rasterio.plot import show
from rasterstats import zonal_stats
from rasterio.transfor
浏览 20提问于2022-09-13得票数 1
1回答
具有搜索功能的交互式R图
、、、、
尝试在R中构建一个非常基本的交互式地图。一般来说，我想做以下工作： 根据属性(如名称)将ArcGIS -> Plot中的shapefile导出为R -> Search shapefile中的交互式地图 我一直在使用tmap： shapefile = st_read("C:/Projects/_SANDBOX/R/Shapefiles/XXX.shp")
map = tm_shape(shapefile)+tm_fill(col="blue",alpha=0.3) +tm_borders(col="black")
tmap_mode(
浏览 17提问于2021-09-30得票数 1
我试图使用管道工R包将一个函数包装到REST中。作为输入，函数接受一个shapefile，并在转换后以.shp格式和GeoJSON格式返回shapefile。现在，在装饰师的帮助下，我需要将其增强为web服务。
R文件：
#* Spatial Polygon Object
#* @param pathtoshpfile:character Path to Shapefile with Name
#* @param design:character one or two
#* @post /sprayermap
function(pathtoshpfile, design = c("
浏览 2提问于2020-09-08得票数 1
2回答
我不明白graphfrom是如何工作的，它是不是像数组中的“from”位置一样检查值？或者'from‘是对象的属性名？我在我的终端中尝试了这两种方法，当属性名为数字时，使用from作为属性名。使用它作为数组索引会给我一个“from is undefined”的错误。还有graphfrom.push(to)。它是如何工作的？
function buildGraph(edges) {
let graph = Object.create(null);
function addEdge(from, to) {
if (graph[from] == null) {
浏览 12提问于2018-07-27得票数 0
回答已采纳
我已经成功地使用GIS扩展将某个城市的shapefile加载到Netlogo中，如下所示。
使用以下代码
extensions [ gis ]
globals [ countries-dataset min-map max-map]
patches-own [ mapa ]
to setup
reset-ticks
clear-turtles
clear-patches
clear-drawing
clear-all-plots
clear-output
; Note that setting the coordinate syste
浏览 2提问于2016-03-30得票数 0
2回答
我是R的新手，但不是ArcView的新手。我正在绘制一些双模式数据，并希望将该绘制转换为shapefile。具体地说，如果可能的话，我想转换顶点和边，这样我就可以在ArcView中显示相同的绘图以及属性。
我已经安装了"shapefiles“包，并且我看到了convert.to.shapefile命令，但是帮助中并没有讨论如何将XY坐标分配给顶点。
谢谢,
时间
浏览 0提问于2011-05-18得票数 6
我想要生成1000个随机点在一个特定的ZIP代码表区格式文件使用pyshp。我的代码是：
import shapefile
zctashape = shapefile.Reader('C:/mypath/tl_2019_us_zcta510.shp')
shapefile_len = len(zctashape.shapes())
#identify the index for zcta 84049 and designate number of points.
zcta_to_use = '84049'
pointcount = 1000
for i
浏览 2提问于2020-03-04得票数 0
回答已采纳
1回答
直线与多边形相交坐标
、、、、
我正在使用Python、Shapely和Fiona。考虑到有两个shapefile可用，一个线形shapefile和一个多边形shapefile。
如何获得由交点(用Q标记表示)及其各自坐标组成的最终结果shapefile？
浏览 6提问于2014-12-31得票数 1
我有一个shapefile (SpatialLinesDataFrame)，包含古龙水的所有街道，可以从下载。我将这个@数据与来自外部源的数据合并。我如何绘制这些街道(如果可能的话，在谷歌地图上使用ggmaps)，使每条街道都有不同的颜色(或厚度)，取决于其各自的价值？
到目前为止，我已经这样做了：
shapefile <- readOGR(shapfile, "Strasse", stringsAsFactors=FALSE,
encoding="latin-9")
shp <- spTransform(shapefile, CRS("+
浏览 0提问于2014-07-15得票数 6
回答已采纳
1回答
我有一个点形状文件，表示每个村庄的名称、住户和总人口。我已经成功地将它上传到了Netlogo上，但我很难找到一个解决方案，将每个村庄的家庭和种群与不同的海龟群体联系起来。我希望能够利用每个村庄的人口价值来培育海龟。下面的代码是我目前使用的代码，但希望了解如何关联shapefile中的属性值来创建我的模型。
要显示shapefile，我将使用以下内容：
to display-villages
gis:create-turtles-from-points villages houses [
set size 4
set color red
set shape "
浏览 10提问于2022-02-24得票数 0
我是草丛中的初学者，但我想在两个多边形之间找到最便宜的路径。更准确地说，我想从位于一个多边形(多边形A)边缘的任何一点到位于另一个多边形(多边形B)边缘的任何一点，得到最小的代价。
到目前为止，我使用了函数CostDistance和CostPath of ArcGIS，使用了一个成本栅格，其中每个单元格都有一个成本值，一个用于第一个多边形的shapefile，以及一个用于第二个多边形的shapefile。我也想用草做同样的事。我认为函数r.cost允许这样做。但是我不知道如何在参数中指定草中的两个多边形？
您有一个例子说明如何使用r.cost在R中的两个多边形与包spgrass6
非常感谢你的
浏览 2提问于2014-09-01得票数 1
回答已采纳
我有一个R中的shapefile和一个土地覆盖类的tiff文件。如何在shapefile中找到土地覆盖类别的百分比？例如，在该shapefile中，研究区域中有多少土地？
#example data for cover
rat <- levels(raster2 )[[1]]
rat[["landcover"]] <- c("land","ocean/lake", "rivers")
levels(raster2 ) <- rat
#and a raster
raster1 <- raster(matr
浏览 0提问于2020-08-13得票数 0
1回答
我正在使用多边形shapefile从栅格中提取平均值。栅格和矢量的CRS是相同的，但提取函数也会返回NA值。
library(rgdal)
library(raster)
options(stringsAsFactors = FALSE)
Shapefile <- readOGR("Fishnet_geolocated.shp",layer="Fishnet_geolocated")
plot(Shapefile)
CHM_Napi<-raster("Napi_CRS.tif")
crs(Shapefile)
crs(CHM_Napi
浏览 16提问于2020-05-13得票数 0
回答已采纳
我有两个形状文件。第一种是区域的多边形形状文件，每个区域都有一个唯一的ID，第二种是在shapefile1区域发生事件的点型文件。
我想给点文件的每个事件(行)分配多边形shapefile的相应区域ID。
对如何用R做这件事有什么想法吗？
提前谢谢。
浏览 2提问于2016-06-17得票数 0
回答已采纳
1回答
我想使用ggplot向sf地图添加点。请看下面的代码： ggplot(data = shapefile) +
geom_sf()+
geom_point(x = -70.67,y =-33.45, size = 10, colour = "red") 这段代码适用于我的一个shapefile文件，但不适用于另一个shapefile文件，我不知道为什么。下面是我的第一个shapefile的代码输出： ? 下面是第二个shapefile的代码输出： ? 为什么第二个调用不能识别坐标的潜在原因？我看到的两个图之间唯一的区别是，在第一个图中，经度和纬度是用数字
浏览 46提问于2020-11-07得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
