一次栅格化shapefile的多个字段
栅格化（Rasterization）是将矢量数据转换为栅格数据的过程。在地理信息系统（GIS）中，栅格化常用于将矢量地图数据转换为栅格图像，以便进行空间分析和可视化。
栅格化shapefile的多个字段是指将shapefile文件中的多个属性字段转换为栅格数据的过程。在栅格化过程中，每个属性字段的值将被映射到栅格图像的像素值上，从而实现对属性数据的空间分析和可视化。
栅格化shapefile的多个字段的优势在于可以将多个属性字段的信息整合到一个栅格图像中，方便进行综合分析。例如，可以将地形高度、土地利用类型和人口密度等属性字段栅格化到同一个图像中，以便分析它们之间的关系和空间分布。
栅格化shapefile的多个字段的应用场景包括但不限于以下几个方面：
- 地理空间分析：通过将多个属性字段栅格化到同一个图像中，可以方便地进行地理空间分析，如地形分析、土地利用变化分析等。
- 环境监测：将多个环境指标（如温度、湿度、气压等）栅格化到同一个图像中，可以实现对环境变化的实时监测和分析。
- 城市规划：将人口密度、交通流量、建筑物高度等属性字段栅格化到同一个图像中，可以为城市规划提供决策支持。
- 农业管理：将土壤类型、作物类型、降雨量等属性字段栅格化到同一个图像中，可以实现对农田的精细管理和决策支持。
腾讯云提供了一系列与栅格化相关的产品和服务，包括：
- 腾讯云地理信息服务（Tencent Cloud GIS）：提供了丰富的地理信息处理和分析功能，包括栅格化、空间分析、地图可视化等。详情请参考：腾讯云地理信息服务
- 腾讯云人工智能图像处理（Tencent Cloud AI Image）：提供了图像处理和分析的能力，可以用于栅格化图像的预处理和后处理。详情请参考：腾讯云人工智能图像处理
- 腾讯云大数据分析（Tencent Cloud Big Data）：提供了强大的大数据处理和分析能力，可以用于处理大规模的栅格数据。详情请参考：腾讯云大数据分析
以上是关于栅格化shapefile的多个字段的概念、分类、优势、应用场景以及腾讯云相关产品的介绍。希望对您有所帮助！
相关·内容
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
回答已采纳
我试图使用来自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
2回答
我想在我的shapefile中添加一个z字段。Z值是排水沟中心的高程。我不知道我该怎么做。质心的x和y值在表中。我有这个地区的dem。 ?
浏览 14提问于2019-04-18得票数 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
回答已采纳
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
回答已采纳
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
我从我的同事那里得到了一个地图图形文件，我还得到了一些观测数据，他们的坐标记录在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转换为一个光栅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
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回答
我正在尝试将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
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
我有一个光栅文件，它由三个独立的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
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
回答已采纳
我有一个区域的shapefile，我想要获得该区域不同材料的光谱剖面，例如水，道路，或植被。看起来像这样的东西？ ? 做这件事最好的方法是什么？ 我已经在r中加载了shapefile，但不确定如何从Landsat堆栈中提取值。 ```{r} landsat<-堆栈(B1、B2、B3、B4、B5、B6、B7) roi <- readOGR(“roi”) 形状 提前谢谢你
浏览 7提问于2020-08-30得票数 0
回答已采纳
我有一个我的研究区域的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
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
我正在尝试提取数字化树冠中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
我有一个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
我想要栅格化一个多边形(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
我已经使用包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
4回答
我想知道是否可以在arcmap中使用带有点特征的shapefile作为web地图？
我见过一些例子，但大多数都是将shapefile投影到底层地图上，例如Google地图或ESRI地图。我只是想使用shapefile作为支持缩放的web地图。这个是可能的吗？如果是这样，我该如何实现这一点呢？
到目前为止，我能想到的唯一方法是将shapefile作为服务保存到服务器的arcGIS中……
任何指导和帮助都会非常有帮助。
提前谢谢。
浏览 1提问于2014-05-05得票数 1
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
1回答
我在R中运行物种分布模型，并希望为非洲大陆创建可变栅格，而不是岛屿。我只能找到非洲及其岛屿的shapefiles，而不仅仅是大陆。
1)哪里可以只下载大陆的shapefile？
2)如果没有shapefile，我想手动从我的栅格中删除孤岛。有办法做到这一点吗f.e.是否将某些坐标之间的部分栅格设置为NA？
浏览 16提问于2020-05-27得票数 0
回答已采纳
在android或java中读写shapefile
我如何在Android或java中读写shapefiles？
.shp、.shx等
我找了很久，但什么也没有
我希望读取文件和编辑信息，并写入同一文件或新文件。
这一点很重要，我可以写到文件中。我能读文件，但这不是我的问题
谢谢
浏览 3提问于2017-02-27得票数 2
回答已采纳
我使用modelform_factory生成带有额外表单的模型。对于带有实例的窗体，类型字段是禁用的，对于额外的表单，则启用该字段。
On save()表单不进行验证，因为禁用字段的POST中没有数据。即使是自定义的清洁方法也不能工作().。我想跳过禁用字段的验证，或者保留此字段的实例数据。
models.py
class Attribute(models.Model):
shapefile = models.ForeignKey(Shapefile)
name = models.CharField(max_length=255)
type = models.Integ
浏览 6提问于2013-09-25得票数 1
回答已采纳
1回答
我正在构建一个无头脚本来预处理哨兵图像。有关该会话的信息，设置如下：
import sys
import os
import qgis
from qgis.core import *
#from PyQt4.QtGui import *
QgsApplication.setPrefixPath("/usr", False)
app = QgsApplication([], False)
app.initQgis()
sys.path.append('/usr/share/qgis/python/plugins')
from processing.core.Pr
浏览 1提问于2018-06-06得票数 2
回答已采纳
2回答
我想创建我自己的地图的自定义版本，但是我不知道在哪里可以获得正确的shapefile来创建地图。
多亏了这篇文章，我理解了flex4中的路径，但是在哪里可以获得美联储褐皮书的实际地图却让我摸不着头脑。请帮帮我！
谢谢
浏览 0提问于2011-03-21得票数 0
1回答
我有一个栅格，其中每个属性字段表示火灾风险模型中的不同输入。我正在试着测试我的模型预测火场位置的能力。我用多边形表示两年的火圈(没有用来训练原始模型的年份)。
我想在基于着火多边形的栅格中创建一个二项式存在/缺失字段，其中存在表示由着火多边形覆盖的单元，而不存在表示在研究期间未燃烧的区域。我想包括完全和部分被火多边形覆盖的单元作为存在。在ArcMap或R中有没有好的方法可以做到这一点呢？
我考虑过栅格化我的多边形，但这看起来很笨拙。如果可能的话，我希望避免丢失分辨率。
浏览 5提问于2020-11-14得票数 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
我正在为C#使用电子地理信息系统库，并从一个ShapeFile中获取一个日期字段，该字段被格式化为yyyy-MM-dd，尽管问题是某些ShapeFiles工作正常，而有些则以yyyyMMdd格式返回。只使用一行代码导入此字段：
tempVals.Add(sf.GetAttributeFieldValues(i)[i2].Trim());
ShapeFile的属性表肯定包括这些破折号。
浏览 3提问于2016-03-30得票数 0
回答已采纳
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
我有一个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
我正在使用GDAL将.tif转换为.shp文件。
gdal_polygonize.py input.tif -f “ESRI Shapefile” out.shp
但是，这不起作用，并导致：
Traceback (most recent call last):
File "/usr/local/bin/gdal_polygonize.py", line 166, in <module>
drv = ogr.GetDriverByName(format)
File "/Library/Frameworks/GDAL.framework/Vers
浏览 3提问于2016-05-12得票数 0
现在，我正试图在软件中实现某种深度缓冲区，当我给它写信时，我遇到了一个很大的问题。拥有一个互斥绝对是过分的。所以我创建了一个等于线程数量的互斥变量。我锁定了一个基于当前像素(pixel_index % mutexes_number)的互斥锁，它工作得更好，但仍然非常慢。我想知道在一个真正的GPU里是怎么做到的？有一个聪明的算法或硬件处理吗？
浏览 0提问于2016-01-24得票数 11
回答已采纳
我想要实现一个天气层，它使用GeoServer在地图中显示从-40到40度的温度。我有一个SQL Server表，它看起来类似于：
create table temperature (id int identity primary key, geom geometry, temperature float)
Geom包含一个适用于给定温度测量的点。我在GeoServer中将表添加为图层，但我仍然无法将它呈现出来。现在，我得到的最好的是一个大的红色正方形:D，它似乎呈现为一个纯点，而不是光栅。
如何为光栅的良好显示创建正确的绑定等？
是在SLD级别还是在其他地方完成的？
浏览 1提问于2022-03-14得票数 0
回答已采纳
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
我正在尝试使用KMeans质心来标记/聚集像素来进行土地覆盖分析。我希望只使用sklearn和matplotlib来完成这个任务。目前，我的代码如下所示：
kmeans.fit(band_5)
centroids = kmeans.cluster_centers_
plt.scatter(centroids[:, 0], centroids[:, 1])
band_5的形状是(713,1163)，但是从散点图上我可以看出质心坐标的值远远超过那个形状。
据我所知，KMeans提供的质心需要转换成正确的坐标，然后转换为shapefile，然后在监督的过程中使用它来标记/分组像素。
如何将这些质心转
浏览 9提问于2016-11-02得票数 1
回答已采纳
我已经编写了一个QGIS插件，它显示用户指定的shapefile，然后当用户选择该选项时，插件会自动向shapefile添加一个新字段，用数据填充新字段，然后根据新字段中的数据显示shapefile的副本。
shapefile的新副本将在属性表中显示新字段，但原始版本的属性表不会刷新，因此不会显示新字段。
有没有一种方法可以“刷新”属性表的链接，使其显示更改，而不必从QGIS中删除shapefile，然后重新添加？
谢谢
浏览 43提问于2015-06-05得票数 2
4回答
还有其他条件吗？
、、、、
我使用lambda语句来设置foreach循环：
foreach (Shape shape in shapefile.Where(x=>x.GetMetadata(IgnoreField) != IgnoreValue))
IgnoreField和IgnoreValue是可选的参数。
如果这些字段是空白的(没有使用)，我如何才能改变我的预测来解释这一点呢？还有别的陈述之类的吗？
浏览 5提问于2013-07-18得票数 2
回答已采纳
我想用分水岭作为边界来剪辑光栅数据。我用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
我有数据想要创建一个GIS类型的应用程序，它将具有添加和删除不同类型的层的典型功能。什么是最好的架构方法？
这些数据包括东和诺丁岛的房产位置。我也有GML和Shapefile的军械调查数据。
我知道这是一个非常广泛的问题，但在我看来，主题范围似乎也很广，我不知道该往哪个方向走。
我正在考虑使用Sql2008space和Bing Silverlight控件来可视化maps.To这样做，我是否需要将东方和北部地区转换成GWS84地理数据类型？但是，如果我将shapefile转换为GML，并使用GeomFromGML将所有GML文件导入到sql中，它们将是几何数据类型。这两种类型不是不相容吗？
此外
浏览 3提问于2011-06-18得票数 1
数据集：
我正在处理0.1 x 0.1度规则网格上的空气污染数据的栅格数据集。我还有一个shapefile，其中包含不同土地用途的分类数据。
我正在尝试建立土地利用和空气污染之间的关系，为此，我需要计算R中每个栅格单元中不同的土地利用面积。
有人知道我是怎么做到的吗？
我已经在上面的链接中包含了shapefile和我正在使用的netcdf。
浏览 3提问于2021-05-08得票数 0
我想手动向admin字段集添加一个字段。
函数calcGeometryField返回我想要添加到管理字段集中的字段。
但是，我在管理界面中添加了两次字段。为什么？
class FeatureAdmin(admin.OSMGeoAdmin):
fieldsets = [
('Info', {'fields':['shapefile', 'id_relat', 'attribute_value']})
]
list_display = ['shapefile',
浏览 2提问于2014-03-17得票数 3
回答已采纳
我在网上搜索了三天后发布了这个问题，但没有成功。希望能在这里得到答案。请不要删除这篇文章，因为我在这里也没有找到答案。谢谢。
我有两个档案：
光栅图像文件(即气温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
2回答
如何对面积变量的矢量进行栅格化？
、、、、
我多次使用了同样的光栅化过程，效果相当好：
raster <- rasterize(vect(shapefile.shp), base_grid, "my_variable")
其中光栅是光栅化的形状文件，shapefile.shp是原始矢量，base_grid是光栅骨架，"my“是要考虑的变量。对于与多边形面积无关的变量，这种方法是令人满意的，因为它使用平均计算来重新排列数据(例如:人口、生产、产量、温度、降水量)。
然而，现在我试图从矢量转换为栅格多边形，这些多边形具有可变的收获面积，这不是严格意义上的多边形的区域，而是可以被认为是与整个多边形面积成正比的区域
浏览 17提问于2022-04-11得票数 0
回答已采纳
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
3回答
融合R中面要素的边界
、、、、
我一直在寻找一个现有的R函数，用于聚合同一层中共享公共边界的多边形特征(即，生成类似于ArcGIS中的‘融合边界’工具的输出)。
我使用gdal_polygonizeR ()从光栅文件创建了一个多边形层。一些多边形形状由单个栅格单元分隔，因此在shapefile中存储为不同的特征。我想将这样的多边形特征组合成一个多边形特征，并创建一个新的shapefile (减少多边形元素的总数)，理想情况下具有用于融合的阈值距离。
有没有人知道在R中有这样做的现有方法？
更新:我认为解决方案可能涉及aggregate，然后是disaggregate。我目前正在探索这一点，特别注意确保带孔的多边形特征与父多边
浏览 1提问于2018-10-21得票数 1
1回答
绘制大型sf数据
、、、、
我需要工作10分钟(1/12度)的全球网格所有的土地面积。网格是在R中使用软件包(sf)生成的。电网将限于世界陆地地区。进一步的下游分析需要网格ID。生成网格的代码如下：
library(tidyverse)
library(sf)
sf_use_s2(F)
birds <- st_read('BOTW_breeding_valid_union.gpkg') ## This is just an examplary shapefile I use to set a bbox.
bbox <- st_bbox(birds)
bbox[1] <- -180
bb
浏览 8提问于2022-08-23得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
