地理信息系统

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

3
集成块GIS的替代方案?[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗?更新问题,使其仅通过编辑此帖子来关注一个问题。 3年前关闭。 我父亲想寻找歧管的替代品的问题是: “我是一名独立的环境顾问,他是GIS的大量用户,负责空间数据处理,地理参考自数字化地图,使用LIDAR进行地形建模以及分析GPS现场数据。我拥有一个工作站8核Intel i7 3.4 GHz,16 Gb内存,NVIDIA GeForce GTX 570运行64位Windows 7旗舰版以及用于系统和缓存文件的专用HD。多年来,由于直观的GUI,广泛的功能,可编程性,本机64位和较低版本,Manifold GIS一直是我的首选软件价格。但是,目前的版本(集成块8)现在是4岁,现在是许多任务非常缓慢。见这里用于最近的矢量处理比较,其中Manifold的输出速度比ArcGIS慢40倍。尽管Manifold具有令人印象深刻的曲面渲染和分析功能套件,但它处理LIDAR数据的速度非常慢,几乎无法用于除非常小的数据集以外的任何数据。歧管启用了CUDA,可用于有限的地形处理功能集(请参见手册>“索引”>“ N”>“ NVIDIA”),而不是实际渲染。下一个版本Manifold 9花费大量时间的一个原因可能是因为正在为CUDA编码大量软件,每个人都承认这很困难。已经承诺了Manifold的忠实用户群将大大加快和改进Manifold 9,但是在2年内没有出现此新版本的迹象。我和许多歧管用户一样,对软件的任何开发都开始感到失望。为了保护我的业务并保持领先地位,我被迫寻找其他GIS解决方案。 歧管的替代方案必须具有以下特征: 以64位本机运行,因此可寻址大量RAM(即<4Gb,通常为16Gb及以上)。 能够使用GPU进行图像/表面密集处理,例如Nvidia图形卡上的CUDA。 地形建模和LIDAR处理能力,能够“叠加”矢量和栅格地图组件,并“飞越”视频演示。 是一个稳定的版本。 带有“表面扩展件”的歧管8专业版价格为590美元(370英镑)。到目前为止,我所看到的所有选项都无法胜任。这些包括: ArcGIS –我刚刚被报价为每年1495英镑+ 425英镑以进行升级,另外还有2495英镑的3D Analyst扩展。 QGIS –似乎仅在稳定的1.8.0 Windows版本中支持32位,有故障,并且缺乏内部数字地形映射功能。 R的空间包–对于我想做的事情,此选项似乎太笨拙和静态。 期待面临类似困境的人们的回应–这些选择中最好的解决方案是什么?或者我应该完全追求其他东西?”


1
如何在R中创建一个有吸引力的Choropleth贴图?
我一直在尝试使用R中的Choropleth映射技术,获得了ggplot2在功能强大的数据分析包中创建漂亮图形的好处。 当涉及到映射时,我无法产生与我通常从QGIS获得的结果相当的结果。使用Sheffield中的流数据(此处是完全可复制的数据和.qgs文件),QGIS轻松生成了以下内容: 使用R(使用此处描述的代码)可以产生的最好的效果是这两张地图,第一张使用GISTools,第二张使用ggplot2: 对我来说,如果只能解决行错误的问题(可能是fortify()命令有问题,或者不使用此处介绍的readOGR()读取shapefile),则ggplot2选项似乎更具吸引力。 所以问题是两方面的:ggplot2选项是R中最佳的Choropleth映射解决方案吗?如果是的话,我该如何解决白线错误的问题? 找出我所做的可复制代码在这里。 编辑-由于使choropleth()输出选项更具吸引力: 可以通过以更高的分辨率导出并删除图例来减少此选项的丑陋性(其他R图例可用)。仍然无法改变线的粗细,但它越来越好了!红线代表流向雇用超过5000名员工的区域。
18 r  thematic-map 

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 

4
如何在QGIS中使用物理屏障/边界进行核密度估计?
有谁知道一个程序,该程序将允许内核在计算时将边界考虑在内,而不是简单地掩盖不可能的区域? 到目前为止,我已经发现:Arc10的GME插件(Hawthorne Beyer)-每次指定边界shapefile时都会出错。我尝试了很多不同的shapefile类型和边界复杂性等。当我不指定边界时,它可以正常工作。 R中的AdehabitatHR软件包(Calange,2011年)-效果很好,但您指定的边界必须非常简单-线段的长度是内核带宽的3倍,而不是太曲折。对于我的数据,这是一个极大的简化。 所以我想知道是否有其他软件可以做到这一点,例如GRASS或QGIS。 谢谢
18 qgis  grass 

4
从ArcGIS REST服务下载数据集
在此信息页面上查看esri.com上托管的数据集:http ://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ImageCat_NLT/MapServer/layers 试图弄清楚如何访问原始数据(原始经纬度坐标或SHP)。 看起来数据可用,我只是想不出如何从此接口中提取数据。

1
我们可以在Postgis2.0中的栅格表上创建视图吗?
是否可以在包含栅格数据的表上创建PostgreSQL视图? 我已经成功创建了视图,但是输入raster_columns不正确(缺少所有栅格信息)。另外,我无法使用QGIS打开SQL视图(很可能是因为raster_columns不正确)。 这与gemetry_columns和视图上的问题的情况相同,但是对于栅格列,我找不到在创建视图期间强制typemod的正确语法。 我所能找到的就是: 光栅(ARRAY ['8BUI','8BUI','8BUI'],4326,0.5,-0.5,...); 但是在“ ...”之后出现的是一个谜。 是否有人设法在raster_columns中创建带有正确条目的视图? 一个简短的例子将不胜感激。

6
从QGIS中的加权点生成密度图的最佳方法是什么?
我真的很想重复ArcGIS Kernel Density工具的功能-我该怎么做? 我想将答案限制在易于安装的免费和开源软件上,即QGIS可以在所有平台上轻松地通过GRASS进行安装,因此QGIS + GRASS可以,但是SAGA不能(不幸的是,因为它具有是理想的工具)。 我正在尝试制作QGIS保护区内的野生动物密度图。以下示例中的密度图是在ArcGIS中根据野生动植物的点观测生成的,原始大小的组大小(矢量文件中的字段)用于加权内核密度(SpatialAnalyst)工具中每个网格单元的密度,选定的搜索半径和网格像元大小: 在先前关于密度的问题中,建议使用GRASS v.kernel工具来模仿ArcGIS Kernel Density工具,但是v.kernel不能完成相同的工作。看完手册并(成功地)生成密度图后,似乎v.kernel函数仅适用于点密度,并且没有机会为每个点提供变量(例如原始计数)以加权每个点。 更新 *.surf.*GRASS中似乎有各种工具可以帮助创建密度曲面-并且这些工具确实接受加权列或z值,或者在栅格上执行。@underdark建议v.surf.rst-“ zcolumn”将是我的权重(计数)变量-但我不知道如何要求该工具创建特定的网格大小或使用特定的半径。 有关如何使用v.surf.rst或其他任何方法的建议? 样本数据 x,y,count 431250,8707500,0 418750,8707500,5 413750,8707500,3 411250,8707500,1 408750,8685000,0 411250,8685000,0 416250,8685000,0 416250,8682500,6 411250,8682500,3 418750,8680000,0 433750,8677500,3 421250,8677500,0 423750,8675000,1 431250,8672500,0 428750,8672500,2 426250,8672500,2 423750,8670000,0
18 qgis  grass  density 

2
如何将带有自定义数据的KML文件导入Postgres / PostGIS数据库?
按照惯例,我并不是说KML文件的结构不同,但是在本<description></description>节中,我具有带有此类HTML标签的数据: <?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Document> <Placemark> <name>New York City</name> <description><![CDATA[ <table> <tr><td colspan='2'>Attributes</td></tr> <tr><td>CITY_ID</td><td>150335</td></tr> <tr><td>POPULATION</td><td>8244910</td></tr> <tr><td>AREA</td><td>468.48</td></tr> </table>]]> </description> <Point> <coordinates>40.712964,-74.003886,0</coordinates> </Point> </Placemark> </Document> </kml> 我想轻松地将每个地标的所有这些数据导入数据库中具有相应列的表中,该怎么做?
18 postgis  kml  ogr2ogr 


3
如何在R中泛化SpatialPolygons和守恒拓扑?
我想对SpatialPolygonsDataFrame R中的A进行泛化,以使绘图更加有效。我所产生的多边形会产生大量的pdf图,这些图加载缓慢且难以集成到文档中。 我尝试使用程序包中的gSimplify函数rgeos,但不幸的是,它不能保存多边形的拓扑。 有没有我可以使用的替代方法?


10
使用ArcGIS Desktop重新排序文件地理数据库中的字段(永久)?
使用ArcGIS 10.0,是否可以对文件地理数据库中的字段重新排序,并且该顺序可以持久保存到数据库中? 我们考虑了导出为XML,并编写了一个解析器来更改字段顺序,然后根据XML重新创建一个新数据库。这似乎是一个费力的过程。 另一个想法是将现有字段定义导入Excel,对其重新排序,导出为CSV,然后编写Python脚本以在新的地理数据库中创建这些字段。在我们从头开始编写脚本之前,是否存在这样的脚本? 不幸的是,ArcGIS Diagrammer似乎不支持字段重新排序。 是否可以使用任何脚本或技巧来简化此过程?



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.