避免在raster::extract(rst，shp)中出现for循环

文章/答案/技术大牛

发布

1回答

r、performance、raster、sf

我当前的代码如下所示： for (b in c(1:nrow(buildings_shp))){ building <- buildings_shp[b,] buffered<- st_buffer(building, 3) raster_cropped <- crop(raster, extent(buffered)) mean <- extract我想知道是否有方法可以通过以下方式加速这一过程： a)避免<

浏览 16提问于2020-10-29得票数 0

回答已采纳

2回答

从光栅中提取值

r、raster、shapefile、sf

sfc_POINT”的函数‘提取’的继承方法# Clean environment library(sf) shp = st_read("PATH_TO_MY_SHP") #Random samplingRdmsamp = st_s

浏览 2提问于2022-01-22得票数 -1

回答已采纳

2回答

将多个栅格的像素值提取到R中的相同csv

r、loops、csv

我有多个栅格散布在一个地区和一个shp与一些多边形的一些栅格。like this image, there are many others 我希望在不丢失shp属性值的情况下，将每个图像中这些多边形的中值提取到一个相同的csv中。我尝试使用for循环，但返回的csv只有列表中图像的值1。dir<-"C:/Users" listImages <- lapply(list.files(dir, pattern =

浏览 38提问于2021-02-04得票数 0

回答已采纳

1回答

使用由R中的多个多边形组成的shapefile从NetCDF中提取数据？

r、dataframe、mask、shapefile、netcdf

但是，我感兴趣的是根据shapefile的每个多边形提取数据(在本例中为12)，并使用多边形的名称保存数据帧。下面是我的示例代码library(rgdal)library(maptools)Subbasin_SHP=readOGR("my_Shapefile.shp")crs(Subbasi

浏览 2提问于2019-06-25得票数 3

1回答

整理shps inR的列表

r-raster、rgdal、rasterize

我一直试图对在R中导入的一个shp列表进行光栅化。当我调用这个列表时，会出现一个错误，因为它不是识别为dataframe，而是作为一个字符。错误:无法找到用于签名“字符”、“RasterLayer”的函数“栅格化”的继承方法shp_ELC <- readOGR("EconomicLandConsesions/ELC.shp") utm <

浏览 0提问于2018-02-15得票数 1

回答已采纳

1回答

将嵌套的for-循环工作到R代码中以减少冗余

r、for-loop

我的县档案是肯塔基州的，可以从埃斯里的ArcGIS在线数据服务中获得。我有可以工作的代码，但我想对其进行改进。$") 如何将嵌套的for循环放置到现有代码中以处理一些冗

浏览 3提问于2019-10-22得票数 2

回答已采纳

1回答

使用过多内存从{raster}包中提取

r、memory、extract、raster

我一直在使用raster包中的extract函数，使用shapefiles定义的区域从光栅文件中提取数据。但是，我遇到了这个进程现在所需的内存量问题。栅格文件较大(约1.6 The )shp <- mclapply(list.files(pattern="*.shp",full.names=TRUE), readShapePoly,mc.cores=6) ndvi <- raster("NDVI.dat"

浏览 2提问于2013-03-29得票数 3

回答已采纳

1回答

一种提取R中栅格面积的有效方法

r、raster、sp、rgdal

DEM_velox <- velox(NED)有没有更好更快的方法来做这件事呢？我知道raster包中有area函数，rgeos包中有gArea函数，但据我所知，这些函数是用于多边形(Shapefile)的，而不是光栅。这方面的任何帮助都是很好的。

浏览 0提问于2020-02-15得票数 0

1回答

从栅格提取平均值期间返回的NAs

r、extract、raster

library(rgdal)options(stringsAsFactors = FALSE)Napi_extract <- raster::extract(CHM_Napi, # the raster that you wish to extract values from) summary(Napi_extra

浏览 16提问于2020-05-13得票数 0

回答已采纳

3回答

使用R完成简单的图像/模式识别任务？

r、image-processing、plot、ocr

(我的下一个问题将是-当有许多线的图像时，是否有一种方法可以检测这些线在哪里相交/“相互接触”)#riddle 1 (find

浏览 2提问于2010-06-29得票数 4

回答已采纳

1回答

使用shapefile提取多个栅格像素值并使用半径列表对其进行循环

python、r、gis、raster、r-raster

我试图解决的问题是，我想要使用包含多个点位置的shapefile提取多个栅格的像素值，对于每个点，我想提取多个栅格的值“时间序列”，并循环通过要在提取函数中输入的缓冲区半径列表，最终结果是每个缓冲区的每个站点的csv 到目前为止，我用栅格库在R下编写了一个程序，我的代码由两个循环组成，第一个循环进入每个栅格，并提取与shapefile中多个位置对应的像素值，第二个循环包含要循环的半径列表，以提取每个像素半径的平均值(path='J:/

浏览 19提问于2020-07-05得票数 0

回答已采纳

2回答

在多个地块中将shapefile叠加到栅格上

r、plot、gis、raster、shapefile

我尝试使用如下代码将南亚的shapefile覆盖在多个光栅图的顶部：'a‘是一个多层光栅文件。filepath/", pattern = "\\.tif$", full=TRUE)south_asia <- readOGR('/filepath/south_asia.shp') #to import shapefile cropped <- crop(x = s, y = extent(south_asia))

浏览 3提问于2020-01-10得票数 3

1回答

创建R循环以从目录中读取shapefile并对每个shapefile执行分区统计

r、for-loop、gis、lapply

我在一个名为"Counties“的目录中有120个县shapefiles。<- raster("C:/path/to/raster/noaa_normal.tif") files <- list.files(dir, pattern="[.]shp$", full.names=TRUE) all

浏览 11提问于2019-10-01得票数 0

回答已采纳

3回答

计算多边形各部分到最近点的距离。

r、leaflet、gis、r-leaflet

在下面的代码中，我使用来自gCentroid()包的rgeos来计算多边形质心，然后在质心周围绘制一个缓冲区。例如，在多边形单元A中，我显示了两个假设的光栅单元，并表示了我所要计算的直线距离。library(rgdal) url <- "https://www.dropbox.com/s/25n9

浏览 2提问于2018-01-16得票数 4

回答已采纳

2回答

raster:：extract生成一个空列表

r、raster、lidar

我正在尝试提取数字化树冠中DSM(CHM)中的像素值。首先，我在shapefile和raster中设置了我的工作目录。TreeCrowns <-shapefile("plot1sag_shape/plot1sag.shp") CHM <- raster('272280split4.tif') 然后我尝试提取像素值 pixel<- raster::extract(CH

浏览 25提问于2020-08-27得票数 0

2回答

只输出最后结果的循环。

r、tidyverse、geo

(path="~/world/", pattern="*.tif$", full.name=TRUE, all.files=TRUE)lon 在数据帧上创建一个新列，该列中插入文件名我尝试过使用这个for循环，虽然理论上看上去还不错，但我不明白为什么它只输出给funvar最后的结果。我不喜欢它覆盖结果而不是附加结果。如果我用一个类似的循环来处理变异对象和更简单的对象，它会追加它们，所以我不确定问题是什么frast

浏览 4提问于2020-08-02得票数 1

回答已采纳

1回答

将HDF5转换为栅格格式时出现问题

raster、sp、rhdf5

R代码和hf5文件以及转换后的光栅图像显示在link (附件)中。Extract Area of Interest e = extent(shp) plot(r1_mask,axes = TRUE,ext = extent(shp)) # E

浏览 39提问于2020-07-08得票数 1

1回答

如何在R中将多个Shapefiles写入ASCII

loops、ascii

我有附加的代码，用于转换为ASCII格式的单个多边形的工作，但我有27000个多边形，所以我需要能够循环通过他们所有。目前，这个循环获取相同的多边形，并将其写入27084次，每次都会覆盖输出。我需要做什么更改才能让它接受目录中的每个多边形，并为每个多边形编写一个不同的*.ascii文件？library("rgdal")library(&quo

浏览 0提问于2018-07-01得票数 0

2回答

如何从年度光栅层堆栈中提取值并使用时间最近的栅格

r、datetime、stack、extract、raster

6 1 13 f 4143885 2761944 1 2010-11-18 08:59:00#create raster

浏览 8提问于2022-02-14得票数 0

1回答

如何在R中使用for循环将多个栅格裁剪到相同的程度？

r、loops、crop、raster

在这个脚本的未来迭代中，我将拥有超过4个栅格，所以我试图编写一个循环，将目录中的所有栅格裁剪到相同的程度。<- paste0(here::here("data", "s2_rasters", raster_files)) e <- extent(wp_shp)

浏览 25提问于2021-04-17得票数 0

回答已采纳

点击加载更多