如何调试使用单个shapefile屏蔽多个栅格的for循环？
调试使用单个shapefile屏蔽多个栅格的for循环的过程如下：
- 首先，确保你已经安装了相应的GIS软件，如ArcGIS、QGIS等，并且熟悉其基本操作。
- 导入shapefile和栅格数据：使用GIS软件打开shapefile和栅格数据，确保它们在同一个坐标系统下，并且能够正确显示。
- 编写代码：选择你熟悉的编程语言，如Python，使用相应的GIS库（如ArcPy、GDAL等）来编写代码。首先，导入所需的库和模块，然后使用循环遍历shapefile中的每个要素。
- 获取要素的几何信息：对于每个要素，使用相应的库函数获取其几何信息，如边界框、几何类型等。
- 根据要素的几何信息，确定要屏蔽的栅格范围：根据要素的几何信息，计算出要屏蔽的栅格的范围。这可以通过栅格的行列号、坐标等方式来确定。
- 屏蔽栅格：使用相应的库函数，将确定的栅格范围设置为无效值或者将其屏蔽。
- 保存结果：根据需要，将屏蔽后的栅格数据保存为新的文件。
- 调试和测试：运行代码，并根据需要进行调试和测试，确保代码能够正确地屏蔽指定的栅格。
需要注意的是，以上步骤中的具体实现方式和库函数可能因使用的GIS软件和编程语言而有所不同。因此，在实际操作中，你需要根据自己的情况进行相应的调整和修改。
推荐的腾讯云相关产品：腾讯云地理信息服务（Tencent Cloud Location Service），该产品提供了丰富的地理信息服务能力，包括地理编码、逆地理编码、路径规划等，可以帮助开发者更方便地处理地理数据。
产品介绍链接地址：https://cloud.tencent.com/product/lbs
相关·内容
1回答
我是r的新手，所以我还在摸索基础知识。我正在尝试获取多个shapefile的列表，并使用它们来裁剪几个栅格堆栈。然后使用裁剪的光栅堆栈，并从它们创建直方图。我知道我必须列出shapefiles的列表，我已经想好了怎么做。我还想出了如何单独裁剪光栅和制作直方图。不幸的是，这个循环把我搞糊涂了。有人有什么建议吗？
这是我到目前为止所知道的：
library('rgdal')
library('raster')
library('landsat')
library('ggplot2')
library('rasterVis
浏览 0提问于2016-01-11得票数 0
我想用分水岭作为边界来剪辑光栅数据。我用Fiona和Rasterio试过了。
这里有我的代码：
import fiona
import rasterio
from rasterio.rio.clip import clip
with fiona.open("oreto_bacino2.shp", "r") as shapefile:
geoms = [feature["geometry"] for feature in shapefile]
with rasterio.open("cn.asc") as src:
浏览 0提问于2018-03-09得票数 4
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
回答已采纳
1回答
我正在尝试使用30米(1角秒)分辨率NED (国家高程数据集)，按州提取5州区域5000英尺以下的区域。我有一个州的shapefile和一个大的(大约4 4GB) NED的光栅。我使用的是velox包，它可以显著提高速度，但我觉得必须有一种更快的方法(这仍然需要很长时间)。
使用velox包，我将栅格转换为velox对象，然后使用velox的提取函数提取像元。因为我知道每个单元格的分辨率，所以我计算每个州5000英尺以下的单元格数量，然后乘以每个单元格的面积。
DEM_velox <- velox(NED)
DEM_cells <- DEM$extract(DEM_velox,sp
浏览 0提问于2020-02-15得票数 0
我从我的同事那里得到了一个地图图形文件，我还得到了一些观测数据，他们的坐标记录在WGS84中。我读取了R中的shapefile：
# read shapefile
> shpmap <- readOGR(dsn = "soil_type", layer = "soil_type_data")
OGR data source with driver: ESRI Shapefile
Source: "soil_type", layer: "soil_type_data"
with 18542 features
It ha
浏览 0提问于2015-10-13得票数 1
我想要栅格化一个多边形(shapefile；1和2的值；1 = outside a protected area，2 = inside a protected area)，到一个1分钟分辨率的栅格。
我使用的是栅格化。
有人能告诉我使用fun= first或fun= last是什么意思吗？每种方法产生的结果略有不同。我到处都找不到任何解释。如果我使用fun=min, max, mean，结果也会略有不同。
浏览 4提问于2016-04-28得票数 1
我有一个区域的shapefile，我想要获得该区域不同材料的光谱剖面，例如水，道路，或植被。看起来像这样的东西？ ? 做这件事最好的方法是什么？ 我已经在r中加载了shapefile，但不确定如何从Landsat堆栈中提取值。 ```{r} landsat<-堆栈(B1、B2、B3、B4、B5、B6、B7) roi <- readOGR(“roi”) 形状 提前谢谢你
浏览 7提问于2020-08-30得票数 0
回答已采纳
1回答
我从worldwildlife.org下载了一个用于世界陆地生态区域的shapefile。可以在此处加载该文件：。
它是一个标准的图形文件，我想用它做两件事。首先:从我的本地目录中获取shapefile并将其裁剪到北美东部的某个范围(ext=范围(-95，-50，24，63))
# Read shapefile using package "maptools"
eco_shp <- readShapeLines("F:/01_2013/Ecoregions/Global/wwf_terr_ecos.shp",
浏览 5提问于2013-02-21得票数 9
回答已采纳
1回答
我正在尝试将csv文件转换为geotiff。我的.csv文件在1列中包含多边形坐标，在另一列中包含值。例如： id,polygon,val
1,[(2.43534,23.353453),(2.13534,22.353453),(2.53534,22.053453),...],12.15
2,[(12.43534,23.353453),(12.13534,22.353453),(12.53534,22.053453),...],16.15
... 我尝试的方法是-我为每个多边形(每行)创建一个shapefile，并尝试将所有shapefile组合到一个栅格中-这是可能的吗？ 另外，请建议其他可
浏览 30提问于2019-03-20得票数 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
我有一个我的研究区域的shapefile，我想创建一个栅格，其中属于该shapefile的每个单元格都有value=1，其他单元格都是value=NA。 我想创建这个新的栅格，与我将在分析中使用的环境栅格具有相同的原点和像元大小。 这是我正在使用的代码，但我总是得到这样的错误："Error in rep(value，length.out = 2)：尝试复制类型为‘S4’的对象“ #Create raster from shapefile
library(raster)
library(rgdal)
x <- readOGR(dsn="D:/Mestrado/Tese
浏览 12提问于2021-02-09得票数 0
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 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
回答已采纳
2回答
我想在我的shapefile中添加一个z字段。Z值是排水沟中心的高程。我不知道我该怎么做。质心的x和y值在表中。我有这个地区的dem。 ?
浏览 14提问于2019-04-18得票数 0
4回答
我想知道是否可以在arcmap中使用带有点特征的shapefile作为web地图？
我见过一些例子，但大多数都是将shapefile投影到底层地图上，例如Google地图或ESRI地图。我只是想使用shapefile作为支持缩放的web地图。这个是可能的吗？如果是这样，我该如何实现这一点呢？
到目前为止，我能想到的唯一方法是将shapefile作为服务保存到服务器的arcGIS中……
任何指导和帮助都会非常有帮助。
提前谢谢。
浏览 1提问于2014-05-05得票数 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
我有一个光栅文件，它由三个独立的2D-numpy数组(数据、经度、经度)和一个shapefile (.shp)组成，它们是从一个NetCDF文件导入的。我想使用shapefile作为掩码来计算掩码中所有数据点的平均值、中位数和标准差。在Python中做这件事最好的地方是什么？ This is how it looks like on the map import matplotlib.pyplot as plt
import geopandas
plume = geopandas.read_file(shapefile)
fig = plt.figure()
ax = fig.add_sub
浏览 13提问于2021-10-12得票数 0
我正在尝试将一个矢量shapefile转换为一个光栅tiff文件。我改编了Michael Diener的Python地理空间分析Cookbook中的源代码。该代码适用于线和多边形shapefile，但仅显示点shapefile的黑屏。 def main(shapefile):
#making the shapefile as an object.
input_shp = ogr.Open(shapefile)
#getting layer information of shapefile.
shp_layer = input_shp.GetLayer()
#pixel_size det
浏览 303提问于2020-01-20得票数 3
1回答
我有一个五大湖的海岸线形状文件和一张安大略湖的地理参考湖岸照片。我想把海岸线的shapefile叠加成一张tiff照片作为一个新的栅格层。具体地说，我希望将shapefile和tiff照片的像素与相同的坐标进行匹配。结果将是一个m x n( tiff照片的一个通道的尺寸)栅格层，其中，如果坐标与tiff文件匹配，则值为1，如果不匹配，则值为0。阅读完这些文件后，我应该怎么做？我可以在Python中使用基本的gdal、rasterio和geopandas。我搜索了一周的现有答案，但没有找到一个。 # read tiff file
ds = gdal.Open(path)
band = ds.G
浏览 52提问于2021-09-12得票数 0
回答已采纳
我有一个.shp文件，我想把它转换成GEOTIFF。我的形状文件由一个大多边形组成，里面有许多多边形。我使用了下面的代码，但是输出的TIF只包含一个大的多边形。 from osgeo import ogr, gdal
import subprocess
InputVector = Shapefile
OutputImage = OutputfileName
gdalformat = 'GTiff'
datatype = gdal.GDT_Byte
burnVal = 1
Shapefile = ogr.Open(InputVector)
Shapefile_layer
浏览 38提问于2019-03-05得票数 2
回答已采纳
我正在尝试提取数字化树冠中DSM(CHM)中的像素值。 首先，我在shapefile和raster中设置了我的工作目录。 TreeCrowns <-shapefile("plot1sag_shape/plot1sag.shp")
CHM <- raster('272280split4.tif') 然后我尝试提取像素值 pixel <- raster::extract(CHM, TreeCrowns, method= 'simple', weights=FALSE, fun=NULL) 但是我得到了一个空的列表，每个多边形的值都
浏览 25提问于2020-08-27得票数 0
1回答
我正在尝试使用下面的代码在多个光栅图上绘制相同的shapefile： ras <- list.files("filepath", pattern = "\\.tif$", full=TRUE)
s <- stack(ras)
plot(s)
plot(world, add=TRUE) 不幸的是，使用下面的代码添加shapefile 'world‘似乎不起作用。 plot(world, add=TRUE) 任何帮助都将不胜感激。 谢谢
浏览 28提问于2020-01-10得票数 0
我已经使用包raster和rgdal从CSV制作了shapefile，但我想将这些shapefile裁剪到特定的区域，这是我可以做到的。
然而，当我编写try to write.csv时，它不起作用。
另外，我在writeOGR之后丢失了X和Y。
datansd = read.csv("..csv", header=T,sep=",")
### tell R what the coordinates are and define the coordinate system ###
coordinates(datansd) = c("X","
浏览 0提问于2015-04-10得票数 0
我有一个rasterstack，我需要使用shapefile从每个栅格中提取值。
我知道使用“普通”栅格堆栈很容易，但我有一个RasterStackTS对象。不幸的是，我必须管理一个RasterStackTS对象，因为我已经从以前的计算中获得了if。
这里有一个可重复的例子，取自Babak Naimi的。
library(rts)
# location of files
path <- system.file("external", package="rts")
# list of raster files:
lst <- list.files(
浏览 2提问于2015-07-17得票数 0
在android或java中读写shapefile
我如何在Android或java中读写shapefiles？
.shp、.shx等
我找了很久，但什么也没有
我希望读取文件和编辑信息，并写入同一文件或新文件。
这一点很重要，我可以写到文件中。我能读文件，但这不是我的问题
谢谢
浏览 3提问于2017-02-27得票数 2
回答已采纳
我写了一个字符串匹配测试的程序，用来测试cpu的性能。
我只是通过<<<1,1>>>调用内核，一个块包含一个线程，执行时间是430ms，然后我使用一个块两个线程<<<1,2>>>调用内核，执行时间是303ms，最后我通过<<<2,1><<调用内核，两个块和一个线程，时间只有430ms的一半(即215ms)。
块中的线和经线有什么不同？是什么让一个包含两个线程的块比两个块慢，每个块一个线程？
浏览 0提问于2013-05-19得票数 3
1回答
我在R中运行物种分布模型，并希望为非洲大陆创建可变栅格，而不是岛屿。我只能找到非洲及其岛屿的shapefiles，而不仅仅是大陆。
1)哪里可以只下载大陆的shapefile？
2)如果没有shapefile，我想手动从我的栅格中删除孤岛。有办法做到这一点吗f.e.是否将某些坐标之间的部分栅格设置为NA？
浏览 16提问于2020-05-27得票数 0
回答已采纳
1回答
Spatstat中的样方
、、、
我现在正在做探索性分析，目标是绘制一个样本地图，确定是否存在完全的空间随机性(视觉上和使用卡方检验)等。然而，我在绘制样本地图时遇到了问题。我之前问过，这导致我重新投影了我的数据。
下面是我的代码：
library(rgdal) #Brings Spatial Data in R
library(spatstat) # Spatial Statistics
library(maptools)
library(raster)
# Load nyc zip code boundary polygon shapefile
s <- readOGR("/Users/my_name/D
浏览 3提问于2020-12-15得票数 0
2回答
在我的分析中，我必须遍历光栅的单元格id (即单元格号)，这是非可选的，并且是在不同的软件设置(BUGScode)中完成的，这一点我不太熟悉。我的光栅必须完全符合形状文件(我的学习区域)，否则我会循环我没有价值或兴趣的单元格号码，然后.好吧，你知道接下来会发生什么;)
可复制示例
### get a shapefile
GER<-getData("GADM", country="DEU",level=1,download=T,path=getwd())
### make a raster
raster<-raster(ext=extent(GER)
浏览 2提问于2021-04-20得票数 0
回答已采纳
我有一个web应用程序，它应该从用户的磁盘中读取shapefile。我使用一个MultipartFile类(https://docs.spring.io/spring/docs/current/javadoc-api/org/springframework/web/multipart/MultipartFile.html)来上传它。据我所知，从MultipartFile恢复文件路径是不可能的。它给我带来了一个问题，因为我使用以下代码来解释shapefile内容(我附加了一个至关重要的构造函数)： private ShapeFile(final File srcFile) {
浏览 17提问于2019-12-27得票数 1
2回答
我正在使用新的velox提取函数来加速shapefile的栅格提取。
默认情况下，旧栅格软件包的提取函数会返回单元格值的列表，例如，当您使用以下格式时：
val.list <- raster::extract(raster, shapefile)
新的velox包需要一个fun=参数，我无论如何也不能让它返回值：
vx.raster <- velox(raster)
vx.vals <- vx.raster$extract(shapefile, fun=??????)
I have tried:
fun=values (returns error Error during
浏览 0提问于2017-04-17得票数 2
我正在试着做和这篇文章一样的事情： Geojson to shapefile convertion using shapefile library 用户回答了他自己的问题，但我不完全理解他的答案。有没有人能解释一下他是如何将多重多边形转化为多边形的？
浏览 28提问于2019-05-07得票数 0
1回答
使用bbox裁剪R中的图像
、、、、
我有一个shapefile的边界框。我想使用这个bbox来裁剪一个图像，我有一个光栅和shapefile。 我尝试过使用gIntersection和crop (如下所示)。不确定是否需要将bbox转换为多边形才能执行此操作。 raster = county raster
plot_parcel_bb = bbox
county_shp = county shape
#4.Use bounding box to clip the appropriate county imagery tile (raster&vector in QGIS)
#Check class
cl
浏览 48提问于2021-06-26得票数 0
下面的代码实现了我的目标，计算一个家庭范围重叠我的治疗区域的百分比。不过，我有190个家范围。对于创建一个for循环或其他方法来自动化所有归属范围的这个过程，有什么建议吗？谢谢
我不知道如何在这种情况下创造一个可复制的例子，但如果有人能指导我完成这一过程，我会的。
library(sf)
library(plyr)
library(amt)
library(raster)
library(mapview)
######################################################
#read in treatment area shapefile
treat
浏览 6提问于2022-01-19得票数 0
回答已采纳
2回答
显示Shapefile
、、、
我有一个想要显示的。我尝试使用matplotlib来显示它，但得到的结果如下：
然而，当我尝试使用在线网站显示时，我得到了这样的结果；
怎样才能得到第二张图片？
下面是我的代码：
import shapefile
import matplotlib.pyplot as plt
print("Initializing Shapefile")
sf = shapefile.Reader("ap_abl")
apShapes = sf.shapes()
points = apShapes[3].points
print("Shapefile Init
浏览 0提问于2015-05-26得票数 3
回答已采纳
1回答
我有一个充满整个屏幕的HTML5画布，如下所示：
在这种状态下，它是完全缩小的，但我允许用户像在google地图中那样放大和平移。每次视口更改时，我都会将视口内可见的所有内容重新绘制到画布上。(因此，当您放大时，需要绘制的内容较少)。
现在缩放和平移有很大的延迟，因为每次重绘都需要30毫秒。在google chrome中分析时，大部分时间来自context.fill() (用于绘制彩色圆圈)和context.drawImage() (用于在圆圈上绘制有光泽的透明覆盖层)。
所以我的渲染效率很低，但由于绘制只在视口边界改变时发生，所以我不能一次只重绘屏幕的一小部分。
有没有人知道另一种我可
浏览 1提问于2014-08-02得票数 0
我想要一张大桌子，然后分成小一点的。我有下面的代码片段，它通过手动将"NAME"替换为ui00000bvbb.lad15nm中的唯一名称来工作
CREATE TABLE "NAME" AS
SELECT parcels_all_shapefile.* AS parcels
FROM ui00000bvbb INNER JOIN parcels_all_shapefile ON ST_Intersects(ui00000bvbb.wkb_geometry, parcels_all_shapefile.wkb_geometry)
WHERE
浏览 1提问于2018-06-25得票数 4
回答已采纳
如何使用开源库以编程方式将shapefile从AGD66转换为GDA94？我不想使用arcgisscripting，因为我没有许可证。
这需要是自动化的。可以使用bash或python脚本。
这比正常的重投影稍微复杂一些，因为在这些坐标参考系统之间使用了不同的椭球，因此需要使用失真网格。
浏览 0提问于2010-06-04得票数 0
回答已采纳
我正在寻找使用Inner的永久原始SQL查询。
我有这些模型：
class EzMap(models.Model):
layers = models.ManyToManyField(Shapefile, verbose_name='Layers to display', null=True, blank=True)
class Shapefile(models.Model):
filename = models.CharField(max_length=255)
class Feature(models.Model):
shapefile = m
浏览 0提问于2013-07-22得票数 2
回答已采纳
我在网上搜索了三天后发布了这个问题，但没有成功。希望能在这里得到答案。请不要删除这篇文章，因为我在这里也没有找到答案。谢谢。
我有两个档案：
光栅图像文件(即气温2020-01-01.tif)
世界各国边界形状文件(即World_Countries_base_map.shp)
目标：--我想在光栅文件的顶部绘制shapefile，然后以Jpeg文件格式保存这个图，以便最终得到这样的结果：
我是Python中的新手，我使用Spyder编写了这段简单的代码：
# Import needed packages
import os
import rasterio
import
浏览 2提问于2020-05-24得票数 1
1回答
我刚开始接触GDAL，只是把脚弄湿了。
我试着把储存在netCDF中的光栅和shapefiles进行比较。shapefiles是netCDF覆盖区域的子部分，数据集使用的投影略有不同。我将shapefile数据集转换为netCDF投影。netCDF文件包含lat、lon、x、y的栅格数组和一维数组。
现在，我的代码使用gdal.RasterizeLayer将shapefile栅格化为tiff，然后使用gdal.ReprojectImage将其重新投影到新的tiff中。我的问题是，我无法确定如何确定第二个tiff的范围-我需要它来选择netCDF数据的子部分。
下面是我代码的相关部分：
#Ext
浏览 6提问于2014-10-22得票数 2
回答已采纳
1回答
我正在使用gdal库加载一个tiff文件并创建一个shapefile。
当我使用QGIS GUI加载shapefile时，没有关于标高的信息。
我想在变换的同时保持高度。
import os
from osgeo import gdal,ogr,osr,gdalnumeric
import numpy as np
# this allows GDAL to throw Python Exceptions
gdal.UseExceptions()
print "reading tif file..."
try:
ds = gdal.Open( "file.t
浏览 2提问于2014-07-30得票数 1
1回答
如何在shapefile对象后面绘制"raster“对象？这两个点都可以自己打印，但点不能在光栅上打印：
require(rgdal)
require(maptools)
require(raster)
myproj = "+proj=utm +zone=12 +north +ellps=WGS84 +units=m"
shp = readShapeSpatial(fn.shp, proj4string = CRS(myproj))
ras = raster(fn.tif)
plot(ras)
plot(shp, bg="transparent",
浏览 0提问于2012-02-15得票数 12
回答已采纳
1回答
摘要：
我正在尝试计算R中大量多边形的面积，我已经阅读了一些关于如何计算的帖子( & )，但我遇到的问题是我的shapefile太大(1.7 my )，无法导入。如果我不能导入这个文件，我就不能计算多边形的面积。
扩展说明：
我实际上是在计算澳大利亚维多利亚的房产面积。多边形表示这些属性。我从为整个维多利亚下载了VicMaps的。然而，考虑到shapefile的大小，我不得不将搜索范围缩小到一个地方政府区域(LGA)，并计算了多边形面积(仅用于测试)。shapefile大小为15.5MB。
library(raster)
x <- shapefile("D:/Downlo
浏览 1提问于2019-07-14得票数 0
数据集：
我正在处理0.1 x 0.1度规则网格上的空气污染数据的栅格数据集。我还有一个shapefile，其中包含不同土地用途的分类数据。
我正在尝试建立土地利用和空气污染之间的关系，为此，我需要计算R中每个栅格单元中不同的土地利用面积。
有人知道我是怎么做到的吗？
我已经在上面的链接中包含了shapefile和我正在使用的netcdf。
浏览 3提问于2021-05-08得票数 0
3回答
融合R中面要素的边界
、、、、
我一直在寻找一个现有的R函数，用于聚合同一层中共享公共边界的多边形特征(即，生成类似于ArcGIS中的‘融合边界’工具的输出)。
我使用gdal_polygonizeR ()从光栅文件创建了一个多边形层。一些多边形形状由单个栅格单元分隔，因此在shapefile中存储为不同的特征。我想将这样的多边形特征组合成一个多边形特征，并创建一个新的shapefile (减少多边形元素的总数)，理想情况下具有用于融合的阈值距离。
有没有人知道在R中有这样做的现有方法？
更新:我认为解决方案可能涉及aggregate，然后是disaggregate。我目前正在探索这一点，特别注意确保带孔的多边形特征与父多边
浏览 1提问于2018-10-21得票数 1
1回答
我可以读一个rastre和exrec多边形:如果我有这个代码来读取光栅和shapefile：
library(raster)
library(geojsonsf)
library(sf)
library(exactextractr)
r <- raster(matrix(rnorm(10*12), nrow=10), xmn = -180, xmx= 180, ymn = -90, ymx= 90)
myurl <- "http://eric.clst.org/assets/wiki/uploads/Stuff/gz_2010_us_050_00_5
浏览 5提问于2022-04-29得票数 0
1回答
我有一个问题要处理图像。我从我的栅格中裁剪了一些多边形，然后提取了带值：
library(raster)
library(sf)
#Read raster
EX1<-raster::stack("ueg.tif")
#Plot raster
plotRGB(EX1, r = 1, g = 2, b = 3)
#Read shapefile
cxs <- sf::st_read("polygons_cx.shp")
#Clip polygon in raster
rr <- mask(EX1, cxs)
#Extract values
浏览 0提问于2021-06-08得票数 0
2回答
我在一个文件夹里有多个光栅。我需要在一个多边形形状文件上提取每个光栅的平均值(有更多的2500个多边形)。
我遇到了两个函数zonal和extract。它说，提取也可以用于点，线和多边形。这是唯一的区别吗？(预期为“是/否”)
我如何从这些多个栅格中提取平均值，并根据它们的文件名为这些提取的平均值指定不同的列名？
编辑：：
我在某个地方找到了一个代码并实现了它。但这需要花费很长时间，一点进展也没有。
grids <- list.files("my_path", pattern = "*.tif$")
#check the number of files
浏览 0提问于2018-04-20得票数 5
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
