地理信息系统

问题： 我正在尝试遍历所有图层，以查找名为的图层"CADAnnotation"。 如果该层存在，则从mxd中删除该层 笔记： 从独立脚本运行（即不在Arcmap中） "CADAnnotation" 数据类型是CAD注释要素类 "CADAnnotation" 不是在地理数据库中，而是从AutoCAD .dwg创建的 "CADAnnotation" 在名为 "ACAD" 如果"ACAD"可以删除组图层，那也可以删除它"CADAnnotation"。 到目前为止的代码： for item in mxds: print (item) mxd = arcpy.mapping.MapDocument(item) df=arcpy.mapping.ListDataFrames(mxd,"Project Area")[0] for lyr in arcpy.mapping.ListLayers(mxd, "*",df): if lyr.name == "CADAnnotation": print(lyr.dataSource) arcpy.Delete_management("CADAnnotation") print("Layer Deleted") else: pass 代码注释： 我可以找到图层没问题 该行arcpy.Delete_management("CADAnnotation")不起作用会引发错误。 题： 如何删除"CADAnnotation"和/或"ACAD"分组图层？

9 arcpy  arcmap  layers 

我严格要求使用OpenLayers检索GeoServer存储中所有图层的名称。我在这里浏览了这篇文章，但是没有回答我的问题。我还介绍了GeoServer的REST API，如果我使用的是另一种语言，例如C＃，Java或PHP，这将是理想的选择。但是我想完全用JavaScript做到这一点。 更新： 从WFS GetCapabilities文档中可以看到，我可以通过使用过滤器发送WFS GetCapabilities请求来检索所有层的名称。但是我不确定是否可以在GetCapabilities请求中包含过滤器。我也不知道这样的过滤器会是什么样子。 对于如何实现此目标的任何建议将不胜感激。任何替代方法也将被理解。 更新： 我相信我在这里找到了解决方案；“可以将WFS GetCapabilities请求过滤为仅返回与特定名称空间相对应的图层。为此，请将以下代码添加到您的请求中：namespace = <namespace>”。我希望这对其他人有帮助。

9 geoserver  openlayers  layers  ows 

我在ArcGIS Desktop中有一个GeoTIFF栅格，可以从以下两个链接下载该栅格： http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif 如果将它们加载到ArcGIS Desktop中，将会看到较大的统一区域。使用“信息”工具单击这些区域将显示像素值为“无数据”。但是，当您选择在图层的“符号”首选项中将“无数据”值显示为透明（或任何颜色）时，“无数据”值将不遵循该规则。 有趣的是，图像主要区域之外的“无数据”值确实遵循该规则（将颜色设置为红色即可看到该规则），但图像内部的值则不然。 我真的很难尝试解决此问题-根本无法解决该问题。 有人有什么想法吗？

9 arcgis-desktop  raster  arcmap  arcgis-10.0  spatial-analyst 

我在想是否有办法将所有层的特定项目另存为从属文件，以便与协作者轻松共享。 例如，创建一个包含与项目相同结构的压缩或非压缩文件夹： Root_folder/ my_project.qgs \____________myfolder \___Layer 1 \___Layer 2 \___Layer 3 \____________myfolder2 \___Layer 4 \___Layer 5 \___Layer 6 \____________myfolder3 \________myfolder4 \___Layer 6 \___Layer 7 \___Layer 8 如果不存在，就当前的API进行开发是否会很棘手？

9 qgis  qgis-plugins 

我正在尝试为以下距离捕获3种场景： 换乘地铁站，具有2个或更多的相邻站。即，所讨论的站点连接2条或更多的主要路线，并具有2条或更多的相邻站点。 终点站地铁站，只有1个相邻站。这是行尾的车站。 嵌入式地铁站，其中有2个相邻站，两种方法之一。 我正在尝试计算一个可能称为“相邻站点之间的平均距离”的值 该arcpy.GenerateNearTable_analysis()可处理两种选择：距离最接近的特征，并且距离所有功能之间。 有没有人有解决这些情况的聪明方法？注意，在属性表的“ StationType”字段下，每个站都被指定为“交换”，“端子”或“内联”。 添加： 这是一些基于@whuber在评论中建议的伪代码。我还没有时间弄清楚这个问题，因此，如果有人想刺中它，您将获得复选标记！;） 我看了一下NetworkX库，它似乎可以按我的意愿工作。 给定图： A —― B ―― C ―― D | E 以及节点和链接： Nodes = ["A", "B", "C", "D", "E"] Links = [("A", "B"), ("B", "C"), ("C", "D"), ("B", "E")] def myFunction(node): identify the links that node belongs to count the …

9 python  arcpy  network 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我是当地一家小型历史博物馆的实习生。我们的预算为$ 0。项目概况： 我有一堆旧纸质地图。当我说老时，我的意思是1800s以上。它们已被扫描到TIFF中。我对它们进行了地理引用，并导出为GEOTIFF。我已经数字化了我们关注的地图的组成部分（铁路，道路等）。 目标是构建一个Web应用程序（我将ESRI的ArcGIS Online用于带有CSV数据的其他Web地图项目），该应用程序使用户可以选择（过去几年）的图层来查看我们县交通的发展情况。除了图层选择器之外，图层不透明度还需要更改一个选项。 我已经研究了MapBox API，但缺少不透明度滑块。QGIS Web客户端看起来很有希望，并且与ArcGIS Server客户端类似，但是我想知道对于我们的最终用户来说它是否太复杂了。我没有太多的时间来解决这个问题，也没有太多的编程经验。

9 qgis  web-mapping  software-recommendations  open-source-gis  tilemill 

给定一条在地球表面的线，我如何绘制一条垂直于它的线？ 抱歉，这是一个非常简单的问题。我认为这将是一项简单的任务，但事实证明这是违反直觉的。 我从下图中的蓝线开始（请参阅链接-我尚无法上传数据）。我通过计算蓝线（m）的坡度找到一条垂直线，然后绘制另一条坡度为-1 / m的线（绿）。当我在Matlab中绘制线时（使用“图”和“轴相等”），它们看起来像预期的那样垂直。 http://imgur.com/7qMkx 但是，当我将这些线导出到Google Earth（使用KML工具箱）时，它们看起来不再是垂直的（请参见下面的链接；较短的线是上图中的蓝线）。 http://imgur.com/ncJQ7 我知道奇怪的事情发生在曲面上，但是我认为这些线至少应该局部垂直。我怀疑这与Google Earth中的投影有关-特别是网格单元的边长看起来大致相似，而纵向边长为1度，而纬度边长为0.5学位。 因此，总而言之： 寻找垂直线在曲面上有效的方法吗？（即，绘制渐变为-1 / m的线） 在Google Earth图片中，垂直线看起来是否像预期的那样，还是有些奇怪？ 更新： 为了提供更多背景信息：我正在查看从飞机上获取的雷达数据。彩色区域是“斯瓦特”，已记录了观察结果。我在上面的说明中开始的蓝线与条幅平行：这是飞机的飞行线（飞机沿大致西南方向移动）。雷达朝左与飞行线正交的方向看。我试图画一条垂直于飞行线的线。这应该是雷达所看的方向，并且应该整齐地割开条幅。如您所见，事实并非如此。

9 algorithm  spherical-geometry  matlab  matlab-mapping-toolbox 

对于那些对GIS编程方面不熟悉的人来说，我觉得许多出色的工具和资源是不可能的。我对基于GUI的向导感兴趣，这些向导用于使用gdal和ogr之类的库处理地理处理。有小费吗？

9 gdal  ogr  gui 

在一家咨询公司的工作向我介绍了许多热门话题技术和关键词。其中之一是BIM，即建筑信息模型。在肯定其实用性的建筑师，结构和土木工程师以及精通数据的软件开发人员和建筑师中，这绝对是一个受欢迎的话题。关于开发GIS工作流以支持BIM项目可交付成果的讨论很明确。 GIS中BIM的情况如何，我们目前处于什么地位？BIM将如何改善GIS工作流程，可交付成果的质量等。从什么时候开始，GIS分析就开始关注“建筑物内部”。传统上讲，让AutoCAD /微站技术人员负责内部设计，而GIS用于那些外部应用程序（分区，土地使用管理）以及宏。

9 bim 

我刚刚开始使用QGIS。我的目标是根据人口创建乔治亚州人口普查区域的热图。我拥有所有数据-乔治亚州每个人口普查区的人口。 通过导入普查区shapefile并遵循此教程，我已经能够成功创建热图：http : //qgis.spatialthoughts.com/2012/03/using-tabular-data-in-qgis.html 但是，我坚持如何在其下添加一个外观漂亮的地图图层，以显示常见的地理特征-公路，城市，水道等。我试图从cloudmade中导入OSM数据，但细节太多了。该地图应该是相当高的级别，但是那里的shapefile显示了每条道路。我正在尝试替换MS MapPoint（无2010年人口普查区）。我尝试获得的功能将根据我缩放的距离来提供不同级别的细节。通常只显示主要道路/城市等。我不在乎这是自动还是手动。 有人能指出我正确的方向吗？任何帮助将非常感激。

9 qgis  openstreetmap 

我希望下面选择（突出显示）与网格相交的“项目”层，然后从选定的“网格”中提取一条记录。我不需要保存该层。 注意：在Python控制台窗口的ArcMap中运行脚本。 import arcpy # Make a layer from the feature class arcpy.MakeFeatureLayer_management("project","project_lyr") # Then add a selection to the layer based on location to features in another feature class arcpy.SelectLayerByLocation_management("project_lyr","INTERSECT","Grid") 该代码运行并且没有给出错误，但是它没有选择（突出显示）任何功能。知道为什么吗？运行也真的很慢。 基本上我正在尝试这样做： 这几乎很快就可以运行，但是我想使用python来做同样的事情。我只需要选定要素表中的值，就不需要该层的副本。 解决方案的一部分可以是使用在存储器的工作区，其是 将数据写入内存工作区通常比写入其他格式快得多... 只是不确定如何实现内存中的工作区

9 python  arcmap  select-by-attribute 

我想出了以下表达式，可将5.1234转换为5°7`24.24“ toint（$ x）|| '°'|| toint（（（（$ x）-toint（$ x））* 60）||'`'|| substr（（tostring（（（（（（$ x））-toint（$ x））* 60）-toint（（（（$ x）-toint（$ x））* 60））* 60），1,5）|| “” 问题是：对于5.1234之类的某些点，它可以工作。但是对于其他人却不是。我怀疑问题在于整数转换会舍入小数位数而不是将其截断。 还有其他选择吗？谢谢。

9 qgis  latitude-longitude  calculator 

我已经从加拿大统计局下载了一个GML文件，并希望将其转换为GeoJSON以在Web的D3可视化中使用它的一部分。但是，我找不到能为我做到这一点的转换器（我尝试了http://converter.mygeodata.eu/，但是返回的GeoJSON不包含任何实际的几何形状）。 将GML转换为GeoJSON的最佳方法是什么？

9 geojson  visualisation  gml  d3 

我的背景：经验丰富的程序员，业余历史学家，GIS初学者（所以我正在使用QGIS）。我确实学会了使用QGIS地理参考工具。 输入：假设我有一个类似Wikipedia的图像：File：VoyagesOfRabbanBarSauma.jpg。假定它是一个“裸”的JPG文件，没有任何地理元数据，并且作者无法提供任何此类信息。输入图像已获得GNU FDL的许可，我很高兴在此发布我的衍生作品。 目标：将Wikipedia地图中显示的路径导入QGIS，以便我可以更改投影，叠加高程，为自然地理添加标签等。 警告：一种解决方法可能只是将地图中位置的纬度/经度输入到新的GIS项目中。我不想走这条路，因为我想记录该特定作者对地理信息的估计。 我尝试过的操作：我已将QGIS Georeferencer应用于栅格JPG文件，并在各大洲的边缘有很多点。大陆内部的点对齐一般。 我想知道的是：这张特定地图的扭曲在欧亚大陆的北部边缘非常明显，而且看起来像是标准的投影。在我看来，应该有一种算法方法来识别生成该图的投影和投影参数，并恢复图像像素与经纬度之间的准确非插值转换？ 我对此有两个问题： 在未知时识别Shapefile的坐标系？ 旧图纸上的未知坐标系 这些人似乎在谈论涉及GDAL或QGIS的蛮力方法，而对于GIS初学者来说，似乎有些劳动强度（基本上是地理配准，并且使其适合已知的投影），但我愿意尝试一下，看看是否可以使它们自动化。 这两个GIS.SE上的说明都比我的经验所允许的要高一些，关于如何完全使用这种蛮力方法的任何提示将不胜感激。例如，上面的问题2显示了如何使用gdaltransform将坐标列表转换为经/纬度。输入坐标是xy像素位置吗？

9 qgis  raster  coordinate-system  georeferencing 

shapefile的最大字段长度为254（https://en.wikipedia.org/wiki/Shapefile#Limitations），因此，我遇到输入属性的问题，该属性的大小大于254。 有更长的字段长度有哪些替代方案？

9 qgis  shapefile  limitations