地理信息系统

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

1
从ArcGIS for Desktop的Python窗口导入GDAL？
我正在尝试使用同时使用arcpy和gdal模块的ArcGIS 10.1中的Python窗口运行一些Python代码。但是，当我尝试导入gdal模块时，出现错误： ImportError: No module named osgeo 显然，它找不到该模块，因此我将主要的python site-packages目录添加到了sys.path列表中： sys.path.append(r"C:\Python27\lib\site-packages") 但是，当我尝试运行import osgeo或from osgeo import gdal出现以下错误时： Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <module> _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swig_import_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: …
3
如何处理具有256个以上变量的表？
我正在处理人口普查数据，并下载了几个CSV文件，每个文件都有600列/变量。我想将它们全部存储在可查询的数据库中，但是到目前为止，我尝试过的所有操作（MS Access，Arc地理数据库表）都将表截断为256列。对于非DBA人士可以访问的大型表，是否有解决方案？
10 database  census  csv 
2
如何在Windows中通过Apache代理Geoserver
我在localhost：80上使用apache Web服务器，从localhost：8888上的jetty servlet提供了geoserver服务。我想设置代理以重定向请求，以便将对localhost / geoserver的请求路由到localhost：8888 / geoserver 我正在遵循http://gistutor.com/geoserver/21-intermediate-geoserver-tutorials/38-configuring-geoserver-proxy-for-public-and-remote-data-access.html中的教程，因此我陷入了它说的地方 copy the text below and right click to paste the ProxyPass into your httpd.conf file right above the </VirtualHost> tag. ProxyRequests Off ProxyPreserveHost On <Proxy *> Order deny,allow Allow from all </Proxy> ProxyPass /geoserver http://localhost:8080/geoserver ProxyPassReverse /geoserver http://localhost:8080/geoserver 我找不到在httpd.conf文件中的上方文本插入位置。httpd.conf文件中的虚拟主机指令指向httpd-vhosts.conf文件，如下所示 # # Use …
1
在使用pyqgis遍历矢量层时，如何检查是否选择了要素？
在使用以下代码（从pyqgis cookbook的示例总结）遍历向量层的同时，有没有办法检查是否选择了要素？ provider = vlayer.dataProvider() feat = QgsFeature() allAttrs = provider.attributeIndexes() provider.select(allAttrs) while provider.nextFeature(feat): geom = feat.geometry() attrs = feat.attributeMap() for (k,attr) in attrs.iteritems(): print "%d: %s" % (k, attr.toString()) 另外，我可以使用创建一个选定功能的列表vlayer.selectedFeatures()，但我希望有一种方法可以直接检查每个功能。
10 qgis  python  pyqgis  select 
3
从PostgreSQL数据库中仅选择空间表？
我的数据库同时包含空间表和非空间表，但是我只想从查询中检索空间表。 任何仅选择空间表的建议。“ the_geom”是空间表中的几何列。 否则，可以从其列名中选择表。 我尝试使用此代码select relname from pg_stat_user_tables WHERE schemaname='public'；但是从中我们得到所有的表名。
2
如何在GDAL / OGR Python绑定中捕获PostgreSQL / PostGIS错误？
查询失败后是否可以在OGR中捕获PostgreSQL输出错误消息？ 例如 conn.ExecuteSQL('SELECT * FROM non_existing_table;') 在Postgres中，这将返回错误消息： "Error: relation 'non_existing_table' does not exist." 在ogr中，我们获得了标准： "Error: current transaction is aborted, commands ignored until end of transaction block" 如果我们能够捕获/解析Postgres错误消息，那将是一个额外的好处。考虑到我对SWIG / CPL的了解，我怀疑这有点长远，但是值得一问，看看是否有人有想法。理想情况下，这将是在Python中进行，但是我也会接受C / C ++注释！ 谢谢。
1
如何可视化具有不确定性的方位角数据？
我正在尝试制作一个图，显示在每个点具有不同范围不确定性的方位角数据。这张来自1991年论文的老式人物捕捉了我旨在的“领结情节”构想： 关于如何制作类似数字的任何建议？我是GIS领域的新手，但是我确实可以通过我的大学使用ArcGIS。我的Arc经验仅限于制作地质图，因此我不必做任何过于奇特的事情。 我在Arc和QGIS中浏览了符号系统选项，但没有看到我认为可以完成此工作的任何设置。注意，这不仅仅是旋转方位角的蝴蝶结形状的问题；每个“领结”的角度范围需要不同。 我将我的Python技能评为“高级中级”，而将我的R技能评为“低级中级”，因此如果有必要，我不反对与matplotlib和mpl_toolkits.basemap或类似的库一起破解某些东西。但是我认为在走那条路之前，我会在这里征求建议，以防万一有一个我还没有听说过的GIS领域的简单解决方案。
3
按百分比增长PostGIS几何
我正在寻找一种在PostGIS几何周围添加缓冲区的方法，但是缓冲区的大小应取决于几何的大小。也就是说，我希望每个几何图形都扩大5％。 我的想法是我正在寻找相交的几何，但是我想考虑的每个几何可能有多达5％的误差。 有人知道最好的方法吗？ 该数据库的行数为一百万行，因此我希望它的运行速度相当快。
4
使用python从栅格获取纬度/经度？
我想知道是否有人在不使用ArcGIS的情况下从栅格中获取高程数据，而是以python list或dict？ 我以元组列表的形式获取XY数据： xy34 =[perp_obj[j].CalcPnts(float(i.dist), orientation) for j in range (len(perp_obj))] 我想遍历列表或将其传递给函数或类方法，以获取xy对的相应高度。 我对该主题进行了一些研究，gdal API听起来很有前途。谁能建议我如何处理事情，陷阱和示例代码？ 不能选择GDAL，因为我无法在正在使用的计算机上编辑系统路径变量！ 有人知道其他方法吗？
3
在给定遥测数据的情况下，可以使用哪种开放源代码（命令行/ API）方法对航拍图像进行正射校正和地理参考？
如何使用开源软件拍摄图像并自动对其进行矫正？与每个图像关联的信息是： 相机规格 相机的纬度，经度，海拔 相机方向：面向无人飞行器，给出俯仰角和偏航角 我想输出一个geotiff文件，并希望最终将图像缝合在一起以创建地图。
2
从头开始创建多光谱图像
我想从Cero制作多光谱图像以对其进行一些测试。真的很简单，例如5条完全均匀的带，上面有盐和胡椒粉的噪声，或者中心是不同值的正方形。显然，这只是一堆矩阵，一个多维数组，可以很直接地生成。我想使用python和gdal来实现这一点，但是gdal相当封闭，我根本不了解它。理想情况下，我想创建一个geotiff文件。有人可以帮我吗？一些指针或gdal教程，这是非常温柔的吗？谢谢你们。
10 python  raster  gdal  imagery 
2
如何不自己做数学就将地磁坐标转换为地理坐标？
我有一些地磁数据。也就是说，它们是经度和纬度，但是是指加拿大的磁性北极而不是地理北极。您可能会猜到，我想将它们与其他坐标系中的其他数据进行匹配。 （数据起源于极光椭圆的计算；即，它们起源于地磁，即使在原则上也没有可以参考的地理资源。） 数据将存储在PostGIS中。 这样做的简单方法是找到地磁坐标系的EPSG代码，用该代码标记数据，然后将所有内容神奇地重新投影到我当时需要的任何位置。但是，我在spatialreference.org上找不到包含地磁或磁场的任何EPSG代码。各种Google搜索也都没有结果。 对于我需要的精度水平，任何经修改以具有不同北极的纬度/经度坐标系（例如WGS84）就足够了。因此，如果有一种简单的方法来定义这样的坐标系并将其加载到PostGIS中，那也可以。 澄清度 例如，假设我只是想从UTM 15N转换为WGS84，而不是从地磁坐标转换为地理。我可以将UTM对象未经修改地保存到PostGIS中，并用EPSG代码26915对其进行标记。然后，当我稍后进行查询时，我可以对PostGIS说“请在EPSG代码4326中返回结果”，然后PostGIS将自动转换为WGS84。 同样，我想用一些坐标系标记地磁数据，并让PostGIS在幕后进行重新投影，而无需reproject()手动调用函数。
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.