Questions tagged «geometry»

(1)抽象几何;(2)将地理数据存储在投影坐标中;(3)矢量特征的通用ESRI术语;“几何”对象的类。

6
将shapefile转换为CSV,包括属性和几何?
我有一个shapefile,其中包含60k多个条目,所有这些条目都是具有相应属性(总计,土地所有者名称,税ID#等)的多边形。我最终需要的是一个具有所有这些属性及其对应几何形状的CSV文件(采用KML兼容的xyz格式,而不是WKT格式)。 我知道我可以在Excel中打开.dbf文件并获取属性。我也知道我可以在QGIS中打开shapefile并将数据复制到Excel,这可以获取属性和WKT几何形状。 是否有一种简单的方法可以将shapefile转换为具有属性和Google Earth友好几何图形的CSV(可在Excel中打开)?

3
SELECT DISTINCT在PostGIS几何列上的精度是多少?
我想知道SELECT DISTINCT操作员在PostGIS几何图形上的精度如何。在我的系统上,以下查询给我一个5的计数,这意味着如果插入的点相差小于1e-5,则认为插入的点相等,并且我不确定这是否是PostGIS的功能,这是我安装的问题或错误。 有人知道这是否是预期的行为吗? CREATE TEMP TABLE test (geom geometry); INSERT INTO test VALUES (St_GeomFromText('POINT (0.1 0.1)')), (St_GeomFromText('POINT (0.001 0.001)')), (St_GeomFromText('POINT (0.0001 0.0001)')), (St_GeomFromText('POINT (0.00001 0.00001)')), (St_GeomFromText('POINT (0.000001 0.000001)')), (St_GeomFromText('POINT (0.0000001 0.0000001)')), (St_GeomFromText('POINT (0.00000001 0.00000001)')), (St_GeomFromText('POINT (0.000000001 0.000000001)')); SELECT COUNT(*) FROM (SELECT DISTINCT geom FROM test) AS test; count ------- 5 …

2
如何计算PostGIS中两条线相交的角度?
我想计算在PostGIS中相交的两条线之间的角度。 PostGIS中角度计算的起点似乎是ST_Azimuth-但这需要点作为输入。我的第一个想法是采用相交线的端点,并对这些端点执行方位角计算。这还不够好,因为大多数线要素不是笔直的,而且我对相交处的角度感兴趣。所以我想出了一个嵌套的操作,它通过以下步骤进行: 识别两个线要素表之间的所有交点。 在相交点周围创建一个非常小的缓冲区 确定线要素与缓冲区外部相交的点(如果有多个点,则取第一个点-我真的只对角度是接近0、90或180度感兴趣) 计算这两点的ST_Azimuth。 完整的SQL有点想在这里发布,但是如果您有兴趣的话,我可以在这里进行介绍。(顺便说一句,有没有比继承WITH语句下的所有字段更好的方法了?) 结果看起来不正确,所以我显然做错了: 编辑我在EPSG:3785中重新计算,结果有些不同,但仍然不正确: 我的问题是此过程中的缺陷在哪里。我是否误解了ST_Azimuth的功能?是否存在CRS问题?还有其他东西吗?也许有很多简单的方法可以做到这一点?
19 postgis  geometry 

2
寻找隧道中心线?
我有一些地图文件由代表隧道的“折线”(每条线只是一个顶点列表)组成,我想尝试找到隧道的“中心线”(下面用红色粗略显示)。 过去,我使用Delaunay三角剖分取得了一些成功,但我想避免使用该方法,因为它(通常)不允许轻松/频繁地修改我的地图数据。 关于如何执行此操作的任何想法? 我正在使用相当原始的C ++。

5
比较ArcPy中的两个几何?
我正在尝试比较两个单独的要素类,以识别它们之间的差异(某种diff函数)。我的基本工作流程: 我使用SearchCursor提取几何 使用修改后的两个要素类的几何形状另存为GeoJSON __geo_interface__(从valveLondon获取 return {'type': 'Polygon', 'coordinates': [[((pt.X, pt.Y) if pt else None) for pt in part] for part in self]})。这是为了避免ESRI与游标一起使用的共享几何对象以及无法进行深拷贝的操作(有关gis.stackexchange的某些讨论在此进行了讨论)。 根据唯一标识符检查两个要素类的几何。例如,将FC1 OID1几何与FC2 OID1几何进行比较。为了获得几何形状作为ESRI对象实例,呼叫arcpy.AsShape()(改性具有孔(见上面第2点)与读取的多边形return cls(Array([map(lambda p: Point(*p) if p is not None else Point(), part) for part in coordinates]))。该比较简单地geom1.equals(geom2)为在所指示的Geometry类。 我期望找到约140个几何形状的变化,但是我的脚本坚持认为有430个变化。我试图检查那些GeoJSON表示形式,它们是相同的,但是Geometry类equals()拒绝这么说。 下面是一个示例: >>> geom1geoJSON {'type': 'Polygon', 'coordinates': [[(-122.8423481559999, 47.060497293000083), (-122.84239755599992, …
18 arcpy  geometry 

1
Postgis-将“ geometry_dump”数据类型转换为“ geometry”类型
我正在使用已在其上运行ST_Union()的属性边界的数据集,因此所有多边形在表的一行中都显示为多多边形。我想将它们分开,以便为每种形状得到单独的行/多边形。我试过在其上运行ST_Dump(),它似乎可以分离出数据集中的多边形,但是我无法再在GIS中查看它(我目前正在使用QGIS)。此列的数据类型称其为'geometry_dump',我认为这就是为什么我无法对其进行可视化的原因。当我从QGIS访问我的Postgres数据库时,此新表将列在没有几何图形的数据集的标题下。 有谁知道我可以做些什么将“ geometry_dump”转换为“ geometry”?

6
使用ArcObjects在Esri Geometry和WKT之间转换?
在ArcGIS ArcObjects API中,似乎无法将ArcObjects几何图形转换为知名文本表示(反之亦然)。我唯一能找到的就是转换为WKB(IWkb接口)。 有没有办法在几何对象和WKT之间执行转换,还是我必须自己实现?我的主要目标是实现.NET且无太多外部依赖。


3
我如何从几何字段获取SRID
从几何字段获取SRID的功能是什么。 我已将一个shapefile插入定义了Geometry字段的SQL Server表中。我知道shapefile的投影为Swreff99_1800(EPSG:3011)。 现在我的问题是如何通过使用几何字段的某些OGC方法获得SRID:3011?

3
在ArcMap中删除具有相同几何形状的重复要素?
是否可以唯一标识要素图层中的每个几何? 本质上,我有一个多边形要素图层,其中许多行由相同的几何图形表示。我希望每个独特的几何图形有1行。因此,我正在寻找一些开箱即用的方法来根据几何形状组合属性。 我来自Postgres世界,可以在那里解决此问题,但我不知道如何(或是否有可能)使用ESRI工具以非编程方式完成此任务。

3
仅将某些列导出到ArcGIS for Desktop中的CSV文件?
我已经使用arcpy编写了python脚本,该脚本将面要素类输出到File Geodatabase中。我添加了将属性导出到单独的CSV文件的功能。我正在使用在这篇文章中找到的代码,效果很好。但是,该代码将导出要素类中的每一列。我只是想导出不具有下列名称字段:OBJECTID,Shape,或Shape_Length。 我的CSV文件生成成功,并且正确不包含OBJECTID或Shape_Length字段。但是,该Shape字段被写入文件。写入该字段的示例值是: <geoprocessing describe geometry object object at 0x28CB90A0> 我添加了一行以打印遍历它们的字段名称,令人惊讶的Shape是,它没有被打印。就像ArcGIS正在隐藏它或给它起其他名字一样。 我的函数的代码如下: def exportToTable(): """ Exports the final outputs to a CSV File. """ # Create path to CSV File (note the varialbe outputPath is declared elsewhere). CSVFile = outputPath+'\\FinalOutput.csv' arcpy.AddMessage("Created CSV File: %s" %CSVFile) # Get all fields in …

2
校正间隙和重叠的多边形(地块)几何
我遇到了我们几乎每天都需要的资源,它非常需要维修。这不是我们自己的产品,供应商通常对解决该问题不感兴趣。 我们使用的多边形shapefile包含多达2000个宗地边界。宗地本身并不是一般意义上的宗地,但是尽管如此,它还是非常有效地描述了数据。“宗地”不得重叠或在宗地之间不产生间隙。但是,数据质量非常低,几乎好像是在禁用捕捉的情况下手工绘制了包裹!见下图: 在某些情况下,它比照片中的差很多,或者间隙/重叠是如此之小,以至于在执行空间过程(例如相交)之前,它们是无法识别的。 如何解决这类问题?我已经附加了shapefile,它代表照片中的8个多边形。这8个多边形同时显示了我的差距和重叠问题。可以自动修复2000多个包裹吗? http://temp-share.com/show/dPf3mpihW 注意:我已经探索了ArcMap中的“修复几何”工具,但这只能纠正重叠的多边形。尽管重叠的几何形状是两个问题中最糟糕的一个,但我想尽可能解决这两个问题!

5
术语:多边形还是区域?
在撰写有关GIS几何的文章时,我不确定是否使用术语“多边形”或“区域”。 特别是,“多边形”一词是否暗示特定的几何类型?例如,一位同事向我建议,“多边形”不适用于甜甜圈功能! 但是即使那是真的,有人在乎吗?在这方面是否有任何强烈的共识? 另外,哪个单词可以更好地翻译成其他语言,例如德语和法语?我正在尝试考虑的不仅仅是英语使用者。

2
无效几何的含义是什么
我已在Postgis数据库中导入了一些数据,并且报告了某些几何形状无效(ST_IsValidReason报告了自相交或环自相交)。 我正在执行的查询似乎不受这些几何的无效方面的影响(我仅使用ST_Distance查询)。 几何无效时会破坏什么? 是否可以“自动”修复这些几何(缓冲区(geom,0)或ST_SimplifyPreserveTopology(geom,0.0001))?

2
如何从python控制台绘制多边形?
我一直在访问和重新访问PyQGIS Cookbook中有关几何处理的页面:http ://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html,但似乎无法弄清楚如何获取从Python控制台绘制的多边形。有人可以帮忙吗?

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.