有没有一种灵活的方法可以通过从shapefile读取的多边形或用一组顶点创建的多边形来裁剪栅格？
是的，可以通过使用GDAL库来实现从shapefile读取的多边形或一组顶点创建的多边形来裁剪栅格。GDAL（Geospatial Data Abstraction Library）是一个开源的地理空间数据处理库，它提供了许多用于处理栅格和矢量数据的功能。
在GDAL中，可以使用以下步骤来裁剪栅格：
- 打开栅格数据集：使用GDAL的Open函数打开需要裁剪的栅格数据集。
- 创建裁剪区域：根据从shapefile读取的多边形或一组顶点创建裁剪区域。可以使用GDAL的Geometry类来表示多边形，并使用其方法来添加顶点或从shapefile读取多边形。
- 设置裁剪参数：使用GDAL的TranslateOptions类来设置裁剪参数，包括裁剪区域、输出栅格数据集的格式等。
- 执行裁剪操作：使用GDAL的Translate函数执行裁剪操作，将裁剪后的栅格数据保存到指定的输出路径。
以下是一个示例代码，演示如何使用GDAL来裁剪栅格：
from osgeo import gdal, ogr
# 打开栅格数据集
dataset = gdal.Open('input.tif')
# 创建裁剪区域
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(10, 10)
ring.AddPoint(20, 10)
ring.AddPoint(20, 20)
ring.AddPoint(10, 20)
ring.AddPoint(10, 10)
polygon = ogr.Geometry(ogr.wkbPolygon)
polygon.AddGeometry(ring)
# 设置裁剪参数
options = gdal.TranslateOptions(cutlineDSName='polygon.shp', cutlineLayer='polygon', cutlineWhere="ID=1", format='GTiff')
# 执行裁剪操作
gdal.Translate('output.tif', dataset, options=options)在上面的示例中，我们创建了一个矩形裁剪区域，并将其保存为shapefile文件。然后，我们使用GDAL的Translate函数将输入栅格数据集裁剪为指定的区域，并将裁剪后的数据保存为output.tif文件。
这里推荐腾讯云的云服务器（https://cloud.tencent.com/product/cvm）和对象存储（https://cloud.tencent.com/product/cos）作为云计算相关产品，可以提供稳定可靠的计算和存储资源，以支持裁剪栅格的需求。
相关·内容
我正尝试用python中的多边形来裁剪光栅(ascii或tif文件)。多边形可以从shapefile中读取，也可以由一组顶点创建，如POLYGON((0 0, 1 0, 1 1 , 0 1, 0 0))。 我知道如何使用参数cutlineDSName来裁剪栅格，该参数表示一个shapefile名称。 但是，如果我想使用由一组点定义的简单多边形来剪切栅格呢？
浏览 15提问于2019-05-11得票数 1
我试着用一个小的多边形来裁剪一个大的多多边形形状文件。它使用st_intersection工作，但是这需要很长时间，所以我尝试将多多边形转换为栅格，然后用较小的多边形进行裁剪。
## packages - sorry if I've missed any!
library(raster)
library(rgdal)
library(fasterize)
library(sf)
## load files
shp1 <- st_read("pathtoshp", crs = 27700) # a large multipolygon shapefile to c
浏览 7提问于2022-10-26得票数 0
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
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
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
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
2回答
我有一个shapefile (以后将被称为源文件)，我需要用一个多多边形的shapefile来剪辑它，这样我就可以为每个多边形剪裁一个shapefile。我尝试了geopandas，虽然我能够通过从多个多边形文件中分别选择多边形来裁剪源文件，但是当我试图在多边形上循环以自动化裁剪过程时，我会得到以下错误：
错误: TypeError：“掩码”应该是GeoDataFrame，GeoSeries或(多)多边形，得到<class 'tuple'>
代码：
import geopandas as gpd
source = ('source-shapefile.sh
浏览 9提问于2022-04-17得票数 2
回答已采纳
1回答
在一个关于OpenGL的教程中，提到了这两个词，不幸的是没有一个很好的解释。
为了找出答案，我读了三次这篇文章，但恐怕我是个太普通的人.我的理解是：
视频框是矩形的，摄像机不能拍摄完整的场景，所以相机的矩形视图之外的所有东西都是“剪切”或“剪裁”。
另一方面，“剔除”意味着在进一步的计算中，由于前景中的其他多边形，这些多边形被排除在无法看到的范围内。
我想这就是链接文章试图在摘要项目1和2中描述的。对吗？
我仍然不明白的是以下几点：
4.与裁剪相比，裁剪建立了最高的多边形。在电子游戏中，剪裁方法比裁剪技术使用得更多。
所以..。我有选择，是使用剔除还是剪裁？那我上面的理解是错的.？
5.
浏览 0提问于2019-10-04得票数 5
回答已采纳
我有一个.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
如果我们有一个栅格，比方说一个国家的整数高程数据，和一个多边形形状文件，比方说在那个国家离散了300个河流流域，每个流域都有一个唯一的名称，我们如何才能最容易地为它们得到像这样的输出呢？
basinID, gridcellelev
a, 320
a, 321
a, 320
b, 17
b, 18
b, 19
最繁琐的方法似乎是过滤/转换单个shapefile为300个shapefile，将栅格裁剪300次为300个uniqueID栅格，读回它们，为每个流域生成单独的表格，然后将它们全部rbind()在一起。
另一方面，理想的方式似乎是跳过文件生成，而不是保存xy数据，只使用一个栅格和一个sha
浏览 86提问于2019-07-19得票数 0
回答已采纳
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回答
泰森多边形问题
、、、
我正在尝试运行一个thiessen多边形。但每次我运行它都会得到第二张图片。第一个是每个点的降水量栅格。如何让泰森多边形看起来像德克萨斯州的形状，而不是矩形。
谢谢!
浏览 1提问于2014-03-04得票数 0
我想得到一个Voronoi多边形图，在50X50域内的一组点。所有的voronoi脊都应该裁剪在边界框的边缘，而voronoi顶点应该位于边框的边缘或边框内。
浏览 2提问于2021-06-25得票数 0
3回答
融合R中面要素的边界
、、、、
我一直在寻找一个现有的R函数，用于聚合同一层中共享公共边界的多边形特征(即，生成类似于ArcGIS中的‘融合边界’工具的输出)。
我使用gdal_polygonizeR ()从光栅文件创建了一个多边形层。一些多边形形状由单个栅格单元分隔，因此在shapefile中存储为不同的特征。我想将这样的多边形特征组合成一个多边形特征，并创建一个新的shapefile (减少多边形元素的总数)，理想情况下具有用于融合的阈值距离。
有没有人知道在R中有这样做的现有方法？
更新:我认为解决方案可能涉及aggregate，然后是disaggregate。我目前正在探索这一点，特别注意确保带孔的多边形特征与父多边
浏览 1提问于2018-10-21得票数 1
我有两个多边形shapefile，我想一个接一个地裁剪。我在谷歌上搜索，但我只能找到由边界框裁剪或由多边形裁剪的点，而这不是我需要的。我还在其他编程语言中找到了一些东西，除了R ()。你能帮我一下吗？
谢谢蒂亚戈
浏览 0提问于2015-08-27得票数 1
3回答
我有一些正在尝试优化的OpenGL绘图代码。在决定是否将渲染数据发送到OpenGL之前，它目前正在测试所有图形对象的可见性客户端。(这比听起来容易。它绘制的是2D场景，因此裁剪非常简单:只需针对视口矩形的当前坐标进行测试。)
我突然想到，通过将整个场景传递给OpenGL并让图形处理器负责裁剪，可以极大地简化整个模型。但有时总数可能非常非常复杂，涉及多达100,000个精灵，其中大多数永远不会被渲染，因为它们不在镜头外，为了简单起见，我不希望最终杀死帧率。
我使用的是OpenGL 2.0，我有一个非常简单的顶点着色器和一个复杂得多的片段着色器。是否可以保证，如果顶点着色器运行并为多边形的所有顶点
浏览 2提问于2011-11-01得票数 2
回答已采纳
1回答
裁剪是在齐次剪辑空间中进行的，因此在将顶点乘以模型视图矩阵之后。但是，当我们在剪辑空间中进行裁剪时，我们插值的新顶点的位置在剪辑空间中，那么如何将这些新的顶点坐标转换回视图空间呢？
我可以进行逆变换，但代价似乎很高。那么，在图形管道中，裁剪到底是如何工作的呢？我的意思是，我们怎样才能得到新生成的顶点在视图空间中的位置？
例如，使用OpenGL，我将在顶点着色器中执行以下操作：
layout (location = 0) in vec3 pos; // vertex position in model space
out vec3 posViewSpace;
// M V P are the
浏览 0提问于2020-08-08得票数 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
2回答
我在UTM32中有一个UTM32，在UTM32中有一个矩形的坐标。(这种预测可能并不总是这样，但预测总是一样的)
我只需要使用矩形裁剪图像。
矩形由：(xmin，xmax，ymin，ymax)
699934.584491, 700160.946739, 6168703.00544, 6169364.0093
我知道如何从点生成多边形，如何从多边形生成shapefile，我知道如何使用这些点创建一个蒙面的numpy数组。但是，我不知道如何使用多边形、shapefile或掩码来实际裁剪图像。
我已经看过了描述：
然而，我并不真正理解它，它似乎过于复杂。(就像我不知道什么直方图拉伸应该在那里做，除了
浏览 1提问于2016-07-07得票数 4
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
回答已采纳
2回答
部分凸壳
、
我有一个顶点P1，P2，P3，.....P11的多边形。我的顶点坐标数据类型是双倍的。
我在P1和P7之间也有一条线。
我想要在P1和P7之间创建一个部分凸包，并在P7之后保留原来的多边形顶点。
因此，最终多边形如下；
到目前为止，我将整个多边形转换为凸包，删除凸包中的顶点，并添加壳顶点。它可以很好地处理小多边形，但当顶点数增加时，这样的管理方式并不容易。
我试图寻找适用于部分凸壳的c#算法，但是除了一些研究之外，我什么也找不到。
有什么想法吗？
浏览 2提问于2018-07-22得票数 0
回答已采纳
6回答
我有一个详细的2D多边形(代表一个地理区域)，它是由一组非常大的顶点定义的。我正在寻找一种算法，将简化和平滑多边形，(减少顶点的数目)的约束，区域的结果多边形必须包含详细多边形的所有顶点。
对于上下文，下面是一个复杂多边形边缘的示例：
我的研究：
我发现Ramer Peucker算法可以减少顶点的数量，但是得到的多边形不会包含所有原始多边形的顶点。参见本文
我考虑扩大多边形(我相信这也被称为外向多边形偏移)。我发现了以下问题：和。但我不认为这会大大减少我的多边形的细节。
谢谢你给我的任何建议！
浏览 4提问于2011-02-18得票数 37
2回答
缩放一个多边形以触摸另一个多边形
、、、、
我期待缩放一个凹多边形(特别是，应用一个相对于形状的质心位置的缩放仿射变换到两个轴)，这样它就会与另一个凹多边形相交/接触。每个多边形都由一组坐标定义。
下图显示了一种迭代方法:逐步缩放多边形，直到两个多边形的最近点之间的距离等于零(我正在使用库的DistanceOp.nearestPoints() )。
有没有一种非迭代的方法来做到这一点？一种立即生成所需缩放因子的方法，而无需迭代缩放和检查？
浏览 6提问于2021-03-23得票数 4
我正在尝试将一个矢量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回答
我正在尝试使用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
进口货包：
from rasterio.mask import mask
import geopandas as gpd
打开一个shapefile：
gdf = gpd.read_file(shpfilepath+clipshape)
打开一个栅格文件：
img = rasterio.open(f'{rstfilepath}raw_immutable/SuperView/{SV_filename}{ext}')
然后执行行动：
for poly_gon in gdf.geometry:
out_image, out_transform = mask(img, poly
浏览 7提问于2022-04-19得票数 -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
1回答
这感觉像是已经有一个函数可以轻松完成了，但我找不到一个。
我最终要做的是:我有一个带有3个显示边界框的特征的shapefile；我想使用其中一个特征在另一个shapefile或多边形中选择其中的所有特征。我使用ogr2ogr，使用-clipsrc和-clipsrcwhere标志从我的边界框shapefile中选择我的边界框。这很有效，但我不能包含或排除掉边界的多边形，我也不想裁剪它们。因此，我认为我应该使用-spat标志，只需导入边界框多边形范围。
我知道我可以通过以下命令获得多边形的范围：
polygon = 'mouth'
inDriver = ogr.GetDriverB
浏览 30提问于2017-07-08得票数 11
回答已采纳
1回答
我在Ubuntu17.10上运行Python2.7，通过apt安装osgeov2.2.1。
我的代码加载osgeo，1)尝试创建带有1层和1字段的shapefile，2)创建多边形(4点几何)，然后3)将多边形添加到shapefile的层。在第三个函数使用layer.CreateFeature(特性)之前，一切都没有问题地运行：
driver = ogr.GetDriverByName('ESRI Shapefile')
ds = driver.Open(shapefile_name,-1)
srs = osr.SpatialReference()
srs.ImportFrom
浏览 0提问于2018-03-08得票数 1
当我在R中使用栅格库中的函数rasterize时，控制台中出现以下信息：
"Found 6085 region(s) and 6409 polygon(s)".
shapefile中有6085个多边形，所以我想知道这两个术语是什么意思？
浏览 4提问于2013-04-24得票数 0
3回答
顶点上的快速隐曲面确定
、、、
我有一个非常昂贵的顶点着色器，因此我想尽量减少多边形渲染的数量。
我可以在cpu上进行合理的裁剪，但它不会下降到每个多边形或每个顶点级别，这是我想要的。我的目标是表现。
有什么快速的技术可以做到这一点吗？
编辑:我知道顶点着色器可以改变顶点的位置，但是这并不是顶点着色器昂贵的原因。因此，我被允许计算顶点的位置，在某些管道中，然后在可能的情况下计算其余的每个顶点。
浏览 7提问于2014-08-16得票数 0
回答已采纳
2回答
如何对面积变量的矢量进行栅格化？
、、、、
我多次使用了同样的光栅化过程，效果相当好：
raster <- rasterize(vect(shapefile.shp), base_grid, "my_variable")
其中光栅是光栅化的形状文件，shapefile.shp是原始矢量，base_grid是光栅骨架，"my“是要考虑的变量。对于与多边形面积无关的变量，这种方法是令人满意的，因为它使用平均计算来重新排列数据(例如:人口、生产、产量、温度、降水量)。
然而，现在我试图从矢量转换为栅格多边形，这些多边形具有可变的收获面积，这不是严格意义上的多边形的区域，而是可以被认为是与整个多边形面积成正比的区域
浏览 17提问于2022-04-11得票数 0
回答已采纳
2回答
查找多个多边形的公共轮廓
、、、
我正在尝试寻找一种在多个多边形之间绘制公共轮廓的算法。我的意思是如图所示：
我们有两个矩形(在我的例子中，它们不是矩形，而是大多数角度为直角的多边形)，我正在寻找图像第二部分上的共同轮廓，如红色路径。在我看来，最大的问题是找到新的点，我在图像的第二部分标记为黄色。多边形永远不会与自身相交或接触。我将一个多边形存储为逆时针顺序的点。
我正在寻找一些线索，来源，甚至关键字，我应该谷歌，这可以让我的任务更容易…
编辑:它类似于凸包，但看的是边，而不是顶点，黄点可能在边的连续性上，因为我看着它。
EDIT2:好的，我需要在多边形周围画一个给定大小的边界，但是如果两个多边形比边界大小更近，那么
浏览 2提问于2011-07-10得票数 6
回答已采纳
1回答
多边形图像
所有的多边形都很简单，它们上面没有洞。电路板多边形(P0到P7)
红色多边形(R0到R6)
绿色多边形(G0 G1 P2 G3)
黄色多边形(Y0到Y3)
我想得到新的四个多边形，标记为1到4，多边形1的坐标是(J7 J10 R5 R4)。当我使用多边形裁剪算法时，我可以很容易地得到结果，板的差异(红并绿并黄)。但是当我有超过10,000个多边形时，我需要很长时间才能得到结果。我的多边形是简单的，我的结果多边形也是简单的，结果多边形中也没有洞。
你知道，我可以很容易地用眼睛找到图像中的四个多边形，但如何使用算法来找到它们。谢谢。
浏览 2提问于2013-07-14得票数 0
我的问题很简单。
在着色器代码中是否有HLSL着色器指令或任何方式来访问GPU在索引缓冲区数组中的当前位置？我使用的是XNA Game Studio4.0，它通过GraphicsDevice.DrawIndexedPrimitives()方法将VertexBuffers和IndexBuffers传递给每个网格的图形管道。我知道GPU会在内部遍历IndexBuffer，但我希望可以确定GPU当前在IndexBuffer中的哪个条目：
为了表达我在psuedocode中的意思: IndexBufferi，我想知道GPU当前正在处理的i的值。
这样做的原因是我想知道当前渲染的是哪个多边形，而仅仅知道
浏览 3提问于2013-02-22得票数 2
我一直在使用R中的栅格软件包中的intersect()函数将一个空间多边形数据框(HUC-4流域)剪裁到另一个空间多边形数据框(由科罗拉多州、爱达荷州、蒙大拿州、犹他州和怀俄明州组成的区域)的范围。 我希望保留与要裁剪到的空间数据框重叠的空间多边形的整个范围。使用intersect()裁剪HUC-4分水岭，以便它们不会超出要裁剪到的状态的范围。 我使用的分水岭数据可以从ftp://rockyftp.cr.usgs.gov/vdelivery/Datasets/Staged/Hydrography/WBD/National/GDB/ (WBD_National_GDB.zip)下载。 包括科罗
浏览 13提问于2019-07-02得票数 0
回答已采纳
我正试图为我的计算机图形学课程弄清楚这一点，但我发现它很难理解。
我相信从n个角度的凸多边形裁剪得到的最小角度数是0，例如，如果多边形在裁剪矩形之外，但我不确定这是否有意义。
至于最大角度数，我相信它是2n，就好像每个角尖都在矩形之外。然而，我也不确定这一点。
浏览 0提问于2016-06-26得票数 3
回答已采纳
2回答
假设我们有一组多边形，可以改变摄像机的视角，并且可以在3D环境中转换摄像机。从某些角度来看，其中一些多边形是完全，被一个或几个其他多边形所遮挡。对于每个绘制的帧，我们知道每个多边形的精确坐标，并可以按“对摄像机的距离增加”或“减少到摄像机的距离”的顺序迭代它们。
现在我的问题是：
什么是有效的方法来确定一个多边形是否被其他人完全遮挡，这样我们就可以在绘图过程中跳过它来提高性能。
浏览 0提问于2011-10-02得票数 4
回答已采纳
1回答
我已经成功地使用了DotSpatial.Contains函数来测试我的每个点(180万点)是否位于我的shapefile中。然而，算法非常慢，因为我正在测试大量的点和非常复杂的多边形。下面是一个例子：
图像中的德国边界是我的shapefile，用于选择点(简化，但仍然有14.000个顶点)，红色矩形是我的180万点所在的区域。
为了寻找空间纬度/经度坐标的更快的多边形点测试，我遇到了Raycasting算法：。
我将代码翻译为VB.Net，它运行时没有出错，但是它没有找到任何点/shapefile的交集。我知道lat/long坐标的困难，但在德国地区，lat/long坐标与标准笛卡尔坐标系相匹
浏览 3提问于2015-03-26得票数 1
回答已采纳
1回答
如何将矩形多边形拟合成不规则多边形？或者更好；如何将最大矩形多边形拟合成不规则多边形？
背景：我有一堆重叠的光栅文件，它们代表了星载传感器系统的不同数据点。为了进一步的分析，我想得到最大的矩形多边形(所有四个角与90度角)，它符合所有这些数据集的足迹/轮廓，以掩蔽和裁剪我的数据集。换句话说，如何将最大尺寸的矩形多边形插入绿色轮廓(见下图)？
Idea：首先获取每个数据集的轮廓或足迹，然后溶解所有的脚印以获得轮廓，最终将最大范围的矩形多边形拟合到不规则的多边形中。
浏览 4提问于2020-03-01得票数 1
我在一个文件夹里有多个光栅。我需要在一个像素上提取每个光栅的平均值(我做点和注意多边形)。
它说，提取也可以用于点，线和多边形。但是对于函数来说(有趣=平均)。我想知道通过这种方法，我是否必须乐在其中=“刻薄”？因为我的形状文件里有点而不是多边形..。？点在我的栅格的每个像素上。
你能帮帮我吗?谢谢
这就是我想要做的代码：
test <- list.files("my_path", pattern = "\\.tif$", full=TRUE)
STACK <- stack(ff)
SHAPE <- shapefile("my_pat
浏览 3提问于2021-05-17得票数 0
回答已采纳
我在GeoServer中有imageMoasic栅格层。我直接将我的API中的geotiffs处理成GeoServer实例数据，然后从那里创建栅格层。我以WMS的形式访问栅格层。我的问题是:我的geotiffs覆盖了可以是5-6个州的大多边形，但在客户端，我想要裁剪栅格来显示特定客户端的特定州。例如，如果我的芝加哥客户登录到我的应用程序，我不希望他们看到大多边形的其他部分，而只看到芝加哥。这种情况下有哪些选项？
浏览 18提问于2018-04-27得票数 0
我正在使用Mapbox和Leaflet来绘制，编辑和删除多边形等。我还使用上传shapefiles并将它们添加到地图中。
对于某个应用程序，我需要将层转换为GeoJSON对象，然后使用以下代码将其转换为JSON字符串：
var fcgeojson = layer.toGeoJSON();
var fcgeojsonstring = JSON.stringify(fcgeojson);
我在中有一个可以工作的版本。而且对于大多数shapefile来说，它都运行得很好。然而，在某些情况下，比如我使用时，其中一个多边形的JSON字符串实际上是不完整的。您可以通过单击jsfiddle中的"Te
浏览 7提问于2017-01-27得票数 1
回答已采纳
我是Qgis的新手。我试图导出一个多边形的shapefile，这样我就可以读取它的顶点(例如，用python包，比如fiona或shapely)作为lat和lon坐标(按度表示)。我现在得到的是以米为单位的坐标。我使用的EPSG:4326投影坐标参考系统，虽然我可以改变。任何小费都欢迎。
浏览 0提问于2018-11-27得票数 2
回答已采纳
3回答
如果在三维坐标空间中有5个顶点，如何确定这些顶点的有序性。即顺时针或逆时针方向。
如果我再详细说明一下，
我有一个3D模型，它由一组多边形组成。每个多边形都是顶点的集合，我想要计算多边形曲面的范数。为了计算范数，我必须考虑按逆时针方向排列的顶点。我的问题是给出一组顶点，我如何确定它是顺时针的还是逆时针的？
这是为导航网格生成，我想删除多边形，不能走的代理。为此，我的方法是计算曲面范数(多边形的垂直向量)，并根据二维平面的角度去除多边形。为了计算范数，我应该知道排列顺序点的顺序。因此，对于多边形中给定的一组点，如何确定点的排列顺序。
例如。
polygon1包括Vertex1 = -21.84
浏览 10提问于2010-01-01得票数 1
1回答
部分离屏纹理-性能差异
、、、、
在图片上有两种情况的例子，其中一个纹理多边形正在呈现。全部由Opengl ES2完成。
( A)多边形部分远离视口B)多边形完全在其内部
我的问题是：
是情况‘A’将消耗更少的系统/gpu资源*，因为纹理部分脱离屏幕，或者它的执行是否与我在视图中呈现它一样，为什么？？
*“资源”--意思是速度，而不是记忆。
我知道opengl在呈现纹理之前先计算顶点，如果顶点脱离屏幕/视图端口，它会跳过任何进一步的计算，但是对于部分脱离屏幕的纹理对象，情况是否相同呢？会省略纹理中不可见的部分吗？
浏览 2提问于2014-04-20得票数 1
回答已采纳
1回答
我需要在一个非常规则的多边形上裁剪一条线。多边形将始终是一个矩形，方向沿欧几里德轴，并可能有角“切断”。看看这张图片，就明白我的意思了。
拐角可以被切割，以便在其中一条边上形成一个点，但不能再多了。在某些情况下(大多数情况下，实际上)不会切割角落，我已经实现了Cohen-Sutherland算法来处理这些情况，但我不知道如何在切割角落时做到这一点。我想有一些通用的算法可以用来对多边形裁剪直线，但我对此了解不多。而且，对于我的情况来说，这可能是过度杀伤力。我处理的多边形总是简单的、凹面的和单调的。有没有什么好的(最好是简单的)算法呢？
浏览 1提问于2013-01-26得票数 0
回答已采纳
我需要在gstat库中创建一个用于空间插值的遮罩网格。具体地说，我有不同的随机分布的采样点，我需要创建包围这些点的最小凸多边形。然后，我必须创建一个空间网格，它应该被计算的外壳裁剪，以便将插值限制在这个多边形的范围内。如果有人能给我解释一下详细的程序并提供一些例子，我将不胜感激。提前谢谢你。
浏览 0提问于2012-10-18得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
