地理信息系统

制图师,地理学家和GIS专业人员的问答

3
GDAL和Python:如何获取具有特定值的所有单元格的坐标?
我有一个Arc / Info Binary Grid(具体地说,是ArcGIS流累积栅格),并且我想识别所有具有特定值(或值范围)的像元。最终,我想要一个表示这些单元的点的shapefile。 我可以使用QGIS打开hdr.adf并获得此结果,工作流程为: QGIS>栅格菜单>栅格计算器(用目标值标记所有点) QGIS>栅格菜单>多边形 QGIS>矢量菜单>几何子菜单>多边形形心 编辑质心以删除不需要的多质心(这些质心= 0) 这种方法“完成了任务”,但是对我来说并不有吸引力,因为它创建了两个我必须删除的文件,然后我必须从形心的shapefile中删除不需要的记录(即那些= 0)。 一个现有的问题可以解决这个问题,但是它是针对ArcGIS / ArcPy量身定制的,我想留在FOSS领域。 是否有人拥有可查询栅格像元值的现有GDAL / Python配方/脚本,并且在找到目标值(或目标范围内的值)时,将记录添加到shapefile中?这样不仅可以避免与UI的交互,而且还可以通过一次创建干净的结果。 我通过与Chris Garrard的一个演示文稿进行对比而对此进行了射击,但是光栅工作不在我的操盘手,并且我不想用我的弱代码使问题变得混乱。 如果有人希望使用确切的数据集,我将它作为.zip放在这里。 [编辑注释]留给后代保留。请参阅与om_henners的评论交流。基本上,x / y(行/列)值被翻转。最初的答案是这样的: (y_index, x_index) = np.nonzero(a == 1000) 倒像这样: (x_index, y_index) = np.nonzero(a == 1000) 当我第一次遇到屏幕快照中说明的问题时,我想知道我是否正确地实现了几何图形,并通过在此行中翻转x / y坐标值进行了实验: point.SetPoint(0, x, y) ..如.. point.SetPoint(0, y, x) 但是那没有用。而且我不认为尝试翻转om_henners的Numpy表达式中的值,是错误地认为在任一行上翻转它们都是等效的。我认为真正的问题分别与x_size和y_size值有关,当行和列索引用于计算单元格的点坐标时,将应用和。30-30 [原始编辑] …
12 python  raster  gdal 

1
在QGIS目录中对图层进行排序?
当我从PostGIS将一堆图层加载到QGIS中时,它们似乎以相反的字母顺序加载。 是否有任何工具或插件可让您以其他顺序快速对QGIS TOC进行排序?例如,我可以想到一个可能需要的几个选项: 按几何对图层进行排序,点在最上面,然后是线,多边形和栅格。 按名称对图层排序(字母,反向字母)。 首先按几何排序,然后按名称排序。这将是#1&#2的组合。点图层排在顶部(与#1相同),但随后按字母顺序对点图层进行排序(a点,b点,c点等)。 加载大量图层时,这将有助于快速建立正确的地图生成顺序。
12 qgis  layers 


5
是否可以在我自己的服务器上托管ArcGIS Online的开源替代方案?
这个问题的答案是社区的努力。编辑现有答案以改善此职位。它目前不接受新的答案或互动。 我知道这可能是太广泛的话题。但是我正在寻找与ArcGIS Online类似的Web-GIS门户,该门户可自行安装(可配置)以托管在我的私有服务器中。主要功能应包括: 用户可以上传空间数据(shp,CSV,KML,JSON等)并与之交互 用户可以对其他人提交的图层进行评分,评论和下载 我认为,QGIS服务器,QGIS网络客户端和/或QGIS云可以发展成为可以在任何服务器上部署的全功能云服务。

1
如何将多边形与点图层上的最近点关联?(或者……如何将邮政区域分配给业务员?)
我有一个邮政编码区域的多边形层,我想用它们根据我们的销售人员所在的位置来创建销售区域。因此,我假设我需要为销售人员点图层上最接近的每个多边形进行求解。我想我可以在多边形或质心的最近部分上执行此操作……我认为这不会有太大的不同。 有人可以指出我这种分析类型的方向吗? 谢谢

4
使用pymssql将点插入SQL Server?
我正在通过python脚本从Web服务(Twitter)收集X和Y值。从长远来看,这将持续几个月,我打算停在600万点左右。 最初输入的坐标是WGS84地理信息,但我需要将其转换为预计的WGS Web Mercator。稍后,该表​​将发布到ArcGIS Server地图服务并进行缓存。 这是一个学习Python的私人项目,没有任何期限,并且想知道仅使用SQL Server 的本机空间类型是否是一个好主意? 我当前未经测试的计划: 使用SSMS创建带有GEOMETRY字段设置(和其他一些属性)的表 在我的python脚本中,利用arcpy或pyproj将WGS84中的纬度/经度转换为WGS84 Web Mercator(或者我可以以某种方式避免这种情况,以及使用SQL可以实现的所有方法吗?) 利用pymssql插入记录,然后将这些点插入表的GEOMETRY字段中。 我的问题是,在WGS84中采用一对经纬度,然后将其插入使用SQL Server空间类型的SQL Server表中,并在其中具有结果点层的一种好,简单且有效的方法是什么? WGS84 Web Mercator,以便可以在ArcGIS Desktop 10.1中渲染/查询它们? 我确实可以访问arcpy / ArcSDE 10.1,但希望以此作为不需要ArcSDE的示例。


2
使用还是不使用关系类?
我正在根据大量数据构建政府数据库。现在,我遵循ESRI地方政府数据模型的标准方法,但经过修改以适合特定政府的需求。 我注意到在地方政府数据模型中,ESRI使用了很多关系类。我过去曾经使用过关系类,但是倾向于避免使用它们。这是因为我通常看不到为什么我需要与要素类分开的表中有关形状的属性的原因。 对于我收集的几乎所有数据,我仍然看不到这样做的原因。对于Parcel Fabric,我可以看到一些有用的原因,但我仍然认为最好还是直接将属性保留在要素类中。 现在要清楚一点,我将使用关系类将要素类引用到要素类,这个问题更具体地说明了为什么我将数据存储在可以轻松存储在要素类中的非空间表中。 无论如何,感谢您的投入!

1
如何自动设置图例标签的格式?在QGIS中
我必须使用分位数模式的“样式”->“研究生”来做10张地图(作为开始)。 我需要格式化不同范围的标签,以使图例看起来易于理解:我希望输入1%-7%而不是1.0000-7.0000,而无需输入。 在arcgis中,可以格式化标签。但是我似乎在QGIS 1.8中找不到任何类似的东西。某些插件或其他QGIS版本是否可以? 几年后更新:新版本似乎已解决了该问题,如下所示。
12 qgis  style  legend 

2
SP1 for ArcGIS 10.1具有64位Python扩展,用于背景地理处理-对其他Python使用的影响
如果安装了用于ArcGIS Desktop 10.1sp1的64位背景地理处理,这将如何影响一般使用的Python配置? Esri为ArcGIS 10.1产品发布了SP1,并且为Desktop组件提供了一个64位Python 2.7.2实例,该实例在安装时将Esri 32位Python 2.7调用替换为用于背景地理处理的64位调用(发送的处理)活动的ArcGIS Desktop GUI会话的一侧)。 安装是可选的,即如果没有安装,则ArcGIS Desktop只能继续使用32位用户存储空间。使用它,在Esri程序文件夹中创建一个\ bin64条目,并在“用于桌面背景地理处理(64位)的HKEY_LOCAL_MACHINE \ SOFTWARE \ ESRI \ ArcGIS 10.1 SP1”处建立了用于64位处理的注册表项,然后进行了背景地理处理在64位内存空间中完成。 Esri的32位和现在的64位Python均已安装完整的2.7.2。 在后台地理处理Esri的注册表期间,路径和配置文件将负责调用64位Python可执行文件。在需要时,ArcGIS Desktop会话中保留使用32位Python。 因此,留给每个用户的是决定是否将已安装的64位Python 2.7.2集成到ArcGIS之外的他们的计算需求中。当32位和64位版本的Python都驻留在系统上时,如何最好地做到这一点? 问了另一种方式-Esri负责在其ArcGIS Desktop中集成64位Python,那么对于ArcGIS之外的其余用例,将有哪些方法可以对其进行配置?

4
使分段的道路在QGIS中显示为一条线?
我正在研究QGIS和英国OS Strategi数据集。如果我使用简单的线型渲染道路,它们看起来会很好: 但是,如果我采用一些更高级的方法,则它看起来是错误的,因为该样式应用于道路的每个路段,而不是整条道路(如我们所知),如以下蓝色高速公路上的白线所示: 如何将这些段合并在一起,否则无法将QGIS视为一条线? 我没有将数据用于导航,所以我不介意这些段是否丢失。
12 qgis  road 

5
Windows上的GEOS和Shapely安装
尝试在Windows 64位计算机上安装Shapely时,无法使GEOS库正常工作。 到目前为止,我已经运行了OSGeo4W安装程序,从中安装了GDAL(我相信该软件包中包含geos库)。之后,我检查了C:\OSGeo4W\bin目录中是否有geos_c.dll ,但是我错过了一些配置步骤,或者库不起作用。 我需要Shapely来工作,因此我也pip install shapely安装了GDAL之后也运行了,显然可以使用(尽管找不到GEOS的C库)。 在我的代码中,我可以导入Shapely,但是当我尝试使用它时,出现一条错误消息,告诉我"geos.dll"找不到。
12 qgis  python  gdal  shapely 


2
使用inkscape抛光qgis地图以进行发布
我在修改在QGIS中创建的地图时遇到问题。我需要在地图和轴标签上添加一些文本。显然,标签无法在qgis中旋转,因此我需要在inkscape或相关程序中执行此操作。 如果我将地图导出为pdf,则无法根据要向其提交文章的期刊的指南调整“刻度”的字体大小(纬度和长值)。 我发现了simplesvg工具,该工具提供了非常好的分辨率地图,但是随后我没有沿着地图边界的坐标(或者该工具有办法吗?) 有没有人可以帮助我?我快要死了..:s 谢谢!爱伦
12 qgis 

6
使用ArcGIS Desktop将线划分为线段?
我有美国各州的shapefile。我提取了两个状态之间的边界,得到了一条线段。例如,佐治亚州和北卡罗来纳州之间的边界。 有没有办法将该边界进一步划分为25英里段的单位? 同样,有一种方法可以将线段分成相等长度的段。例如,取该边界线段,并将其分成等长的3个线段。 我使用ArcGIS 10.1和ArcMap。 如果要分割几行(几千行)?

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.