Questions tagged «geometry»

有哪些工具可以从较旧的扫描地图中得出用于现代GIS系统的几何图形？ 我指的是生成实际数据的库，而不是提供用于手动编辑或跟踪功能的界面。 如果没有这样的工具，是否有来自相关领域的学术论文或研究（或无耻的骇客！）可能会阐明该问题？

10 geometry  algorithm  software-recommendations  vectorization 

我想从代表苏格兰城镇和村庄的一系列观点中找到一个“难以接近的观点”。使用FOSS GIS工具（QGIS，SAGA，GRASS，Postgres / PostGIS ...）。 这是离人口中心最远的地方。这将是不包含城镇的最大圆圈的中心，但边缘将具有城镇。 有一种有效的栅格方法-GDAL 接近滤镜。 有基于VECTOR的方法吗？不需要栅格化的一项？ 这是我对栅格的尝试... 从一个点层开始，以米为单位的合适投影（我正在使用OSGB，27700） 添加一个名为“ one”的字段，并将其设置为常数值为1的整数。 向量到栅格，使用-a并使用GEO而不是像素来编辑命令。确保像素为正方形！ 在邻近栅格（在本例中为1 km宽）上创建轮廓 这是它的样子。我用手画了一些候选点；-) 我以为可以解决重复的缓冲区（使用Multi ring缓冲区插件），但是我怀疑我没有足够的内存来使它起作用。

9 qgis  geometry  proximity 

我得到了一个具有多边形形状的数据集，其中的几何由QGIS插件“ Geometry Updater ” 更新。现在我有一个问题，在PostGIS中，我得到了无效的几何图形，并显示了错误消息“孔位于外壳外部”。但是，当在QGIS中观察几何图形时，实际上有两种形状，其中一种是主要几何形状，另一种是（可能是）孔。 我以为原因是几何类型设置为ST_Polygon。但是即使将其更改为ST_MultiPolygon，错误仍然存​​在。有人知道如何解决此问题？也可以在QGIS中手动进行。 编辑： 这是图片中形状的原始WKT（将其更改为之前ST_MultiPolygon）： POLYGON（（681328.211640639 5312556.84866055,681333.876018639 5312566.85338355,681344.126811639 5312584.17146655,681360.13777864 5312605.81168455,681368.571629639 5312620.16794654,681375.15817664 5312633.84801255,681380.74133464 5312647.74367554,681384.01769064 5312658.41482554,681392.01169764 5312694.63471254,681393.07199764 5312701.83700354,681389.16681464 5312705.85459254,681379.79126564 5312704.26135554,681393.93193764 5312715.39845154,681401.81727164 5312723.43382354,681398.09962564 5312710.95653254,681395.19739764 5312693.68448654 ，681387.20535464 5312657.40285554,681383.87846764 5312644.38168654,681378.69327164 5312632.35257854,681367.00831364 5312610.47778954,681346.401239639 5312581.83289054,681336.365640639 5312564.27437955,681340.354365374374365679365468163768563752225764 5312565.30077054,681395.29993964 5312563.90629854,681410.98915364 5312545.67508754,681413.79585564 5312536.98783954,681415.77221864 5312530.93185454,681416.991529641 5312523.86320254,681406.38747164 5312515.12414654,681403.96790864 5312507.57002154,681395.94924164 5312514.30082354,681394.08765664 5312516.71419154,681390.91027864 5312518.71529654,681387.07111564 5312520.63375054,681381.90838164 5312522.38690354,681376.13571264 5312525.04785255,681366.362498639 5312530.79640254,681349.071127639 5312540.57236955,681332.287335639 5312550.05537054 ，681321.992430639 …

9 qgis  postgis  geometry 

给定形状的两个不同分区（出于论证的考虑，一个国家的两个不同行政区划），我如何才能找到一个新的分区，使这两个分区都适合，从而允许（并优化）某些错误？ 例如，忽略该错误，我想要一个执行此操作的算法： 也许有助于用设定的术语表达这一点。使用以下编号： 我可以将上面的分区表示为： A = {{1}，{2}，{3,4,7,8}，{5}，{6}，{9,10,13,14}，{11}，{12}，{15} ，{16}} B = {{1,2,5,6}，{3}，{4}，{7}，{8}，{9}，{10}，{13}，{14}，{11,15} ，{12,16}} 点B = {{1,2,5,6}，{3,4,7,8}，{9,10,13,14}，{11,15}，{12,16}} 并且产生A点B的算法似乎很简单（类似，如果两个元素在A（B）中的分区中合并在一起，则将它们在B（A）中的分区合并-重复直到A和B相等。 但是，现在想象一下，这两个分区之间的某些行略有不同，因此不可能找到完美的答案，而是我希望最佳答案受制于最小化某些错误准则。 举一个新的例子： 在左列，我们有两个没有公共线的分区（除了外部边框本身）。上述唯一可能的解决方案是琐碎的解决方案，即右栏。但是，如果我们允许“模糊”解决方案，那么中间一栏可能是允许的，比如说竞争总面积的5％（即分配给每个粗化分区中的不同子区域）。因此，我们可以将中间列描述为表示“误差小于5％的最小粗略公共分区”。 那么，实际答案是顶部行，中间列还是中间行，中间列中的分区，或者介于两者之间的分区，就不那么重要了。

9 geometry  algorithm 

我SpatialLines在R：中创建了两个对象。 这些对象是通过以下方式创建的： library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) 现在，我想以某种方式得出结论，这是旋转和翻转的同一条线，并且它们之间的差等于0（即形状相等）。 为此，可以使用maptools包装并旋转第1行，例如： spl180 <- maptools::elide(spl1, rotate=180) 然后必须使用rgeos包装将每条旋转的线与2号线进行比较，例如： hdist <- rgeos::gDistance(spl180, spl2, byid=FALSE, hausdorff=TRUE) 但是，这是计算SpatialLines对象匹配的昂贵方法，尤其是在对象数约为1000的情况下。 有什么聪明的方法可以完成这项工作吗？ PS此外，上述方法不能保证所有可能的旋转和翻转。 P.S2。如果＃1线相对于＃2缩小，则＃1和＃2之间的差必须仍等于0。 更新：

9 r  geometry  similarity 

我有一个问题，可以使用ArcGIS Spatial Analyst工具箱中的“ 区域几何”工具来解决。但是我没有Spatial Analyst的许可证，因此我正在寻找替代方案。可能使用QGIS。 如何找到适合不规则多边形的圆的最大半径？ 请注意，多边形可以是凸包或凹包（如下所示），并且解决方案必须同时解决这两个问题。 我尝试了约瑟夫的解决方案，但不幸的是结果不是我想要的。 首先，我确实有非常不规则的多边形，例如： 如果我遵循约瑟夫的描述，结果将如下所示： 这肯定是该解决方案后的结果，但这不是我的问题的答案。 对我来说重要的是要回答一个问题，即圆的最大半径有多大，以便无论圆的中心在哪里，圆都仍完全位于多边形内部。 例如，多边形的北部有更多的空间，因此可以放置一个比多边形的南部大得多的圆。但是，这个圆圈有多大？

9 qgis  geometry  spatial-analysis 

作为背景，此问题源于试图回答有关如何从“表格到椭圆”工具输出的线生成重叠多边形的问题？ 使用ArcMap 10.1 GUI，使用椭圆构造工具将真实的曲线多边形数字化为文件地理数据库要素类非常容易，但是... 从包含中心点，长轴，短轴和方位角的表中读取行（arcpy.da.SearchCursor）时，是否可以编写真正的曲线椭圆多边形？ 我曾希望arcpy.da.InsertCursor可以使用一个令牌来执行此操作，但是SHAPE@似乎受到Geometry对象支持的功能的限制，并且似乎不包含真实的曲线。

9 arcpy  arcgis-10.1  geometry  ellipse  true-curves 

我在墨卡托投影（google maps）上有一条线（Ax，Ay-Bx，By）和最接近该线的随机点（Cx，Cy），我想知道最近的点（图像上的透明蓝色）指向那条线（图像中为蓝色） 编辑：澄清这是在墨卡托投影（球面投影）

9 google-maps  coordinate-system  geometry  php  spherical-geometry 

我使用ENVI EX进行了基于对象的图像分析，并将每个分类特征导出为shapefile。我将这些shapefile转移到ArcMap 10中，并尝试将它们与从GeoBase获得的相同区域的其他shapefile相交。 当我尝试运行相交工具时，收到一条错误消息，指出从ENVI EX提取的shapefile的拓扑无效。我试图打开拓扑工具来修复这些shapefile的拓扑，但是按钮显示为灰色，我似乎无法单击它们。拓扑工具不适用于我拥有的任何shapefile。 有什么办法吗？如何修复拓扑并运行相交？

9 arcgis-desktop  geometry  intersection  topology  envi 

从OSM数据提取道路几何形状的最佳方法是什么？目的是显示由于天气原因封闭的高速公路路段。气象数据由主要城市的位置（沿道路测得的KM）和受影响路段的长度组成。主要问题是基于这些距离来修路。之后渲染似乎很简单。

9 geoprocessing  openstreetmap  geometry 

我有一系列代表建筑物轮廓的经度和纬度坐标 例如 -0.5485381346101759,53.2285150736142 -0.5482220594232723,53.22842450827133 -0.5482298619861881,53.22841205254449 ...（中间点未列出）... -0.5483123769301657,53.22882101914848 如何计算中点？我发现了一些教程，显示了如果您有三个坐标（例如，http : //mathforum.org/library/drmath/view/68373.html），该如何做，但是在许多情况下，我有三个以上的坐标。 谢谢

9 geometry  spherical-geometry  coordinates  mathematics 

由于无效的几何形状，我经常花费很多时间来处理怪异的bug。症状差异很大，包括： 数据库变慢 错误\空查询结果 日志文件超载了错误和警告 与几何算子不一致的结果（例如，交集） 我想对无效几何的类型进行分类，以帮助这里的每个人清理其数据。 请以您最喜欢的无效几何类型回答（每个答案一种）。欢迎使用屏幕截图，说明和解决方案-但请尽量简短。

9 geometry  invalid-data 

我指的是本文在 2.6.1节中，关于两个TIN的加法和减法： 可以精确确定两个TIN的相加，并将其存储到新的TIN中，因为分段线性函数的相加又会产生分段线性函数。通过对T1和进行叠加T2来完成加法，为此有几种算法。之后，我们得到一个细分，其中所有面都有3、4、5、6边。现在，我们必须填写叠加层顶点的高度信息。 尽管我可以理解段落中的每个单词，但我不知道如何在实践中执行上述步骤来获得两个TIN的填充/填充。 更具体地说，我想知道如何执行两个TIN的叠加。本文末尾提供了参考资料，但由于我不在大学图书馆内，所以我无法访问它们。因此，非常感谢任何易于访问的在线参考（或代码示例）！

9 geometry  tin 

从维基百科文章： 在几何中，仿射变换或仿射图或两个向量空间（严格来说是两个仿射空间）之间的亲和力（来自拉丁语，仿射，“与……相连”）由线性变换和平移组成。 谁能给出在GIS中何时以及如何使用它的真实示例？

9 coordinate-system  geometry  affine-transformation 

我想复制一个面要素类，并将所有面在x和y方向上偏移约10英尺。我问上周是否有任何方法可以做到这一点，并且得知我很可能需要使用arcpy创建自己的python脚本。我使用arcpy制作了自己的脚本，但无法正常工作： import arcpy from arcpy import env import os env.overwriteOutput = True # Get arguments: # Input polygon feature class # Output polygon feature class # inputFeatureClass = arcpy.GetParameterAsText(0) outputFeatureClass = arcpy.GetParameterAsText(1) xShift = arcpy.GetParameterAsText(2) yShift = arcpy.GetParameterAsText(3) shapeName = arcpy.Describe(inputFeatureClass).shapeFieldName # Create the output feature class with the same …

9 arcpy  arcgis-10.0  geometry