地理信息系统

有人知道ESRI WKT和OGC WKT投影格式字符串之间的区别的确切清单吗？ 我知道有各种工具可以帮助您从ESRI WKT转换为OGC WKT，包括GDAL实用程序和各种网站服务。但是我的问题不是实用的，我只是想了解这些服务使用的格式/语法差异。先前的Stackexchange问​​题仅讨论了特定示例之间的差异，或者仅涉及可用工具和服务。 即使您只知道一个区别，也可以将其发布。根据我自己的经验，应该只有很少的差异。我知道的区别是： esri定义中的大多数文本项使用下划线，而ogc使用空格。 esri wkt中定义基准的文本与ogc wkt相同，只是以“ D_”开头。 有时，某些预定义的PROJCS，PROJECTION，GEOCCS和DATUM的文本标识符以不同的方式编写（例如，“ NAD83”为“ North_American_1983”）。我猜想知道哪些标识符的拼写不同的唯一方法是拥有列表或查找表，因此请给任何您知道不同的名字命名。 除了ogc每个单词的大写字母大写而esri所有小写的字母之外，各种PARAMETER文本值都相同。但是，我见过没有使用此规则的情况，有人知道在尝试加载它们的软件中标题是否真正重要？ UNIT的类型在ogc中拼写为大写字母，在esri中拼写为小写，例如“ Degree”与“ degree”。在某些情况下，我看到ogc的“米”被拼写为“ meter”和“ m”，而在其他情况下，法语的拼写为“ metre”。有人知道这两种格式或任何其他单位类型的正确约定是什么吗？

9 coordinate-system  wkt  ogc 

有谁知道如何执行这样的overpass-api请求？ “告诉我德国的所有银行”？ 我可以在一个特定的边界框中显示所有“便利设施=银行”。但是德国（或所有大多数国家）在这样的bbox中不太适合。 有什么办法吗？你有例子吗？服务器性能或运行时是次要的，因为我们可以托管自己的服务器或类似的东西:)

9 openstreetmap  filter  overpass-api  search 

该PyQGIS食谱介绍如何设置空间索引，但只说明了它的用法的一半： 创建空间索引—以下代码创建一个空索引 index = QgsSpatialIndex() 向索引添加功能-索引接受QgsFeature对象，并将其添加到内部数据结构中。您可以手动创建对象，也可以使用上一个调用提供程序的nextFeature（）的对象。 index.insertFeature(feat) 一旦空间索引中填充了一些值，就可以进行一些查询 # returns array of feature IDs of five nearest features nearest = index.nearestNeighbor(QgsPoint(25.4, 12.7), 5) 获取属于返回的特征ID的实际特征的最有效步骤是什么？

9 qgis  pyqgis  spatial-index 

如何在不进入现场收集数据的情况下从城市的某个部门获取建筑物的高度值？无论如何，是否可以通过卫星图像，SRTM，ASTER GDEM或任何DEM进行处理？ 值得一提的是，这不是美国城市，也不是具有免费LIDAR数据的城市（位于罗马尼亚）。

9 dem  remote-sensing  3d  image  romania 

我在一个小型计划委员会中工作，是2个GIS人员中的1个，那里有2个ArcGIS许可证。最近有一位志愿者开始与我们合作。他正在使用QGIS并主要编辑shapefile。我可能应该首先考虑这个问题。...ArcGIS和QGIS在一起玩得很好吗？最终，所有内容最终都将在ArcGIS中出现，在QGIS中编辑shapefile或地理数据库，然后在ArcGIS中打开和/或编辑它们是否会遇到任何问题？我不知道是否要在这里提出另一个问题...如果在QGIS中编辑了地图标注，可以在ArcGIS中打开该地图吗？

9 qgis  arcgis-desktop  compatibility 

我有一张等值线的postgis表，其定义如下： CREATE TABLE myisolines ( gid serial NOT NULL, isotime timestamp without timezone, val numeric(10,4), geom geometry(LineString,4326) ); 在外观上，此线串对象如下所示： 我知道数据的空间范围，因此可以添加一个Bbox，这样LineStrings可以是封闭的。 我想myisopolygons从myisolines具有多边形的表格中创建一个等值多边形表格，该表格不会重叠，而是创建一个连续的曲面，并具有一列等值线val最低val的列，由此形成多边形。我知道它可以由自封闭的等值线（孤岛）形成，也可以由用bbox封闭的等值线形成，在这种情况下，val应从特定的等值线中获取。在视觉上，它应该看起来像这样： 我以为可以以某种方式创建拓扑，然后将面转换为多边形，但是我不知道如何正确地进行。如何才能做到这一点？ 另一个选择是在bbox和每个创建的多边形之间递归使用差异函数，但我想这不是正确的方法，而且绝对不快。

9 postgis  postgresql  polygon  convert  linestring 

我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 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 3年前关闭。 我有两层：定义属性的点层（95个对象）和定义与属性相关的土地图的多边形层（211个对象）。每个点都涉及一个或多个多边形。我要生成的是一张表，列出每个点到与该点关联的所有多边形的距离（最近的边或多边形的质心-两者都会做）。在QGIS和ArcGIS中，计算最近距离相对容易，但是这些计算至少在我采用的方法中都忽略了所有遥远的多边形。理想情况下，我希望输出 Point_ID | distance to polygon 1 | distance to polygon 2 | distance to polygon 3 | etc. ArcGIS 10或QGIS 2.2+中的任何指针将最有帮助。

9 qgis  arcgis-desktop  arcgis-10.0  distance 

我正在尝试为丹麦市政当局查找shapefile的质心，然后找到它们之间的行驶时间。我使用R的readShapePoly从功能maptools与组合gCentroid，从功能rgeos，和一切正常。但是，我得到诸如 SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA 我显然不能在Google中使用这来获取旅行时间。我正在寻找一种将这些数字转换为经度-纬度的方法，但不知道如何。 当我readOGR从rgdal库中读取数据时，我得到了相同的坐标，但是它告诉我以下有关我假设是投影的信息（但是坐标是相同的） Slot "proj4string": CRS arguments: +proj=utm +zone=32 +ellps=intl +units=m +no_defs 可重现的示例：我将示例的数据放在这里：https : //github.com/sebastianbarfort/shapefiles 这应该重现该问题： library(maptools) library(rgdal) library(rgeos) map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp") centroid = gCentroid(map) centroid

9 shapefile  r  coordinate-system 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我正在使用OGR / GDAL在python中构建脚本。 我有一组shapefile和一组GeoTiff栅格文件。 如果脚本不与栅格区域相交，我希望脚本忽略它们。 shapefile不是矩形，因此我不能简单地比较layer.GetExtent（）返回的xmin / xmax，ymin / ymax值。我需要代表其整体形状的实际多边形，然后通过某种方式确定该多边形是否与栅格正方形相交。 我当时想我可以以某种方式将shapefile中的所有多边形合并为一个要素，然后读取该要素上的几何图形，然后将该信息与栅格范围进行比较。但是，我不确定具体如何执行此操作。 如何从shapefile中提取边界多边形信息？ 如何确定该多边形是否与给定的正方形相交？

9 python  raster  shapefile  gdal  ogr 

我有两个shapefile，点和折线。我正在尝试用点数据ID（systemID）填充我的折线图层，这些点从线条开始到结束为止。例如。创建两列（从和到）。线路已经分配了流向。我正在使用ArcHydro工具，但Arcpy也很好。 在ArcHydro中，我尝试使用“属性工具”>“为线生成从/到节点”，但是它创建的“从”和“从”值并非来自我的点数据。我的另一种选择是对顶点使用特征，然后将结果连接到我的点数据表中以获得systemID。但是然后我也必须更新折线层。我敢肯定，有一种更简单的方法可以做到这一点。 注意：并非所有线在顶点处都具有点，因此空白值适合它们。

9 arcgis-desktop  arcgis-10.1  arc-hydro 

我有一个带有2个属性列的Excel文件，1个功能信息2个当前功能应符号化的颜色数（RGB）。可以使用唯一的“功能信息”字段（通过属性表连接）将该文件文件与图层连接。因此，我具有每个多边形要素的颜色编号属性作为属性。 我如何使用此信息来符号化每个功能，这些功能的特征色已从excel工作表导入？

9 arcgis-desktop  arcgis-10.1  symbology  excel  rgb 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 8个月前关闭。 我想使用GeoTools创建地图并将其保存到图像（例如JPEG）。我的要求很简单： 创建一个包含两层的世界地图：政治边界和标线。这些图层来自不同的来源和不同的投影。 将地图输出到不同的投影（例如“ EPSG：5070”，“ EPSG：4326”，“ EPSG：54012”，“ EPSG：54009”等） 将输出裁剪到不同的AOI（例如-124.79至-66.9 lon，24.4至49.4 lat）。 我想通过API以编程方式执行此操作。到目前为止，我取得的成就有限。我学会了使用这种方法创建地图并以各种投影形式输出： //Step 1: Create map MapContent map = new MapContent(); map.setTitle("World"); //Step 2: Set projection CoordinateReferenceSystem crs = CRS.decode("EPSG:5070"); //Conic projection over US MapViewport vp = map.getViewport(); vp.setCoordinateReferenceSystem(crs); //Step 3: Add layers to map CoordinateReferenceSystem mapCRS = …

9 geotools 

使用QGIS 2.4 python控制台，我试图仅显示图层功能的子集。我已经设法从图层中获得了两个功能： featuresIterator = layer.getFeatures(QgsFeatureRequest().setFilterExpression ( u'"field"=4')) 并从dataProvider中： featuresIterator = layer.dataProvider().getFeatures(QgsFeatureRequest().setFilterExpression ( u'"field"=4')) （我不确定是否会有所区别）。 我现在想要实现的是仅显示这些功能。使用与GUI中的查询生成器完全相同的方法。我不知道是否必须在渲染级别（将渲染某些功能而不是其他功能）或数据提供者的级别（仅向图层提供仅包含数据的视图的视图）进行管理所需的子集），或介于两者之间的某处。

9 qgis  python 

我正在尝试让QGIS与Python的Anaconda发行版一起工作，但是我不知道如何使它工作。 在一个答案中，提到了QGIS不能与Anaconda一起使用，我想仔细检查一下，以防万一最近几个月有任何更新。 如果可能的话，请您指导我一份资源，说明如何进行此操作？ 我正在使用带有conda 3.7.1和python 2.7.6的Windows 7（64位）。

9 qgis  windows  anaconda