地理信息系统

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

3
如何将OpenLayers连接到PostGIS数据?
我已经在OpenLayers中阅读了如何在OpenLayers中加载KML,GeoJSON等矢量文件。但是如何连接到PostGIS数据。如果不是,为什么我应该将数据上传到PostGIS,而不是直接使用WebServices放置数据。我有什么想念的吗? 我想加载此矢量数据,以便可以通过单击要素来查看属性。我需要PostGIS,因为将来我想启用查询。我认为只有维护像PostGIS这样的数据库才有可能。但是如何使用OpenLayers? 我已经阅读了《 OpenLayers入门指南》,但是没有提到任何地方。请帮忙!!

1
从栅格属性表中的另一个属性重新分配栅格像元值?
我有一个大型的(大小为48 US),分辨率为100m,具有栅格属性表(RAT)。数据是由其他人创建的。RAT中的属性之一具有我关心的数据值,而实际的栅格像元值毫无意义。我一直在尝试找出如何将我关心的属性的值重新分配给实际的栅格像元值,但还无法弄清楚如何做到这一点。我要这样做的原因是,我有一个相当复杂的自定义Python工具,其他人开发了该工具来执行我需要的计算,但是不幸的是它没有使用RAT属性的能力。它仅使用栅格像元分配的值,正如我在这种情况下所说的那样,这是毫无意义的。 我尝试了一种运行栅格转点工具的round回方式,以将目标字段作为我关心的属性来创建点要素类,以为可以从所得的点要素类创建新栅格。但是,栅格指向工具失败,表明我没有足够的临时内存(我正在运行具有16 GB RAM和1 TB数据驱动器的Alienware计算机)。 我有什么选择吗?我使用的是ArcGIS 10,但如果我可以做我在另一个软件包中容易描述的事情,我将愿意尝试另一种方式(我不嫁给ArcGIS)。

1
如何重新投影shapefile?
我想重新投影一个QGIS项目,但是没有这样的选择。当寻找答案时ogr2ogr,通过重新投影每个矢量层(在此示例中为shapefile),我看到了该任务的建议。但是我没有得到预期的结果: C:\temp>ogrinfo -al -so misc.shp INFO: Open of `misc.shp' using driver `ESRI Shapefile' successful. Layer name: misc Geometry: Line String Feature Count: 10 Extent: (21.267388, 42.015857) - (21.270225, 42.017470) Layer SRS WKT: PROJCS["MGI_Balkans_zone_7_deprecated", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",21], PARAMETER["scale_factor",0.9999], PARAMETER["false_easting",7500000], PARAMETER["false_northing",0], UNIT["Meter",1]] id: Integer (10.0) C:\temp>ogr2ogr -t_srs …

2
使用QGIS将不规则多边形分成相等的区域?
我是测量师,并使用QGIS维护土建项目上的土壤测试记录,并不断获得所需的测试批面积或批号。这些形状通常不是不规则的,并且令人怀疑。 是否有插件/解决方案将包含不规则多边形的层切成相等的区域或具有指定区域的手数,结果是根据提供的面积获得的手数最多?

4
如何从[LINE | MULTILINE] STRING提取POINT
我想在其中创建一个vectorLayer,QGIS以显示形成DB LINESTRING或MULTILINESTRING存储在PostGISDB中的所有点。 我认为我需要将LINESTRING(或MULTILINESTRING)的所有点转换为POINT。 第一个问题: 有PostGIS功能吗? 第二个问题: 如果要从结果POINT表中排除重复的点,该怎么办?

2
计算QGIS中相交多边形的面积
我有两个shapefile集合: 1-邮政编码,行政和政治区 2-各种土地利用类别(例如公园和地段)以及点(例如树) 我想按邮递区号,行政区和政治区以及以相同方式划分的点数对这些土地使用类别的面积进行统计。例如,SE15中有4.1公顷的公园,SE22中有879棵树,依此类推。 我正在使用QGIS。我已经处理了原始数据,以准备好所有shapefile,但是在此网站和其他网站上进行了无济于事的搜索之后,我看不到这样做的方法。

2
在Python脚本工具参数中使用MultiValue?
我的脚本与一起选择的单个要素类可以正常工作clip_fc = arcpy.GetParameter(0)。 但是,当我在脚本属性中将clip_fc设置为MultiValue,然后选择多个要素类时,出现错误“参数无效”。 谁能告诉我我要去哪里错了? import arcpy import os.path from arcpy import env arcpy.env.workspace = r"C:\GIS Home\project_1\clipshps\Combined_mf.gdb" arcpy.env.overwriteOutput = True clip_fc = arcpy.GetParameter(0) clip_name = str(clip_fc) output_dir = r"C:\GIS Home\project_1\clipshps\outputs_clip" for fds in arcpy.ListDatasets(): for fc in arcpy.ListFeatureClasses('','',fds): out = os.path.join(output_dir, fc + "_" + clip_name[40:] + "_clip.shp") arcpy.Clip_analysis(fc, clip_fc, …

2
您如何看待PostGIS中的部分3d功能?
我们具有来自调查数据的功能,其中包含部分3d信息。 最常见的示例是表示道路的2D LineString,其中包含在被调查的某些点的高程信息。其他示例包括屋顶形状-MultiLineString,其中某些关键点具有从建筑平面图中分配的标高,但不是全部。 使用PostGIS,您建议使用哪种数据模型来存储此类信息,以使其尽可能可访问,而不会丢失或生成内插信息?
10 postgis  data  3d  model 

2
将地理配准数据添加到TIFF图像中
我想.tiff从.png文件创建地理参考文件。因此,我.tiff使用以下命令行将.png转换为普通文件: convert image.png image.tiff 在这里,一切正常,我有一个有效的.tiff文件,显示正常。然后,我想在其中添加一些地理配准数据。为此,我找到了一个名为的工具tiff2geotiff,我正在这样使用它: tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff 这给了我黑色图像,该图像在QGis中导入,在适当的地方显示良好,除了在各处都是黑色(证明已正确添加了地理参考数据)。 因此,当我这样做时: gdalinfo georeferenced.tiff 我得到: Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength Driver: GTiff/GeoTIFF Files: araujo_tmp.tiff Size is 1585, 979 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS …
10 qgis  gdal  php  geotiff-tiff 

5
在“数据驱动页面”环境中显示动态表?
我正在创建一个地图系列,每张地图有许多多边形。多边形有注释,但是有些多边形很小,没有注释(注释会完全隐藏多边形)。因此,我有一个属性表,其中存储了那些特定多边形的标签。 如何在表格上显示该表格,每张表格仅显示属性值,而不显示所有2000+个值? 我正在寻找一个简单的动态解决方案。

3
是否可以通过Python命令设置项目CRS?
是否可以通过Python命令在QGIS中设置项目的项目CRS? 在我的情况下,用户应打开QGIS,启动我的书面插件,此插件将加载具有epsg代码31467(高斯Krueger区域3)的空间数据。不幸的是,QGIS的默认CRS是WGS84,所以我遇到了单位和比例的麻烦。 我的第一种方法是使用以下代码: my_crs = core.QgsCoordinateReferenceSystem(31467, core.QgsCoordinateReferenceSystem.EpsgCrsId) self.iface.mapCanvas().mapRenderer().setDestinationCrs(my_crs) 这将CRS更改(在已经可视化数据之后)到高斯Krueger区域3。至少右下角的espg代码已更改。单位仍然是度,刻度是完全错误的。高斯·克鲁格的1:5000是WGS 84中的1:528822376(在我的情况下是QGIS)。在属性中手动设置项目CRS之后,项目crs实际上是GK3,并且比例和单位正确。 所以我的问题是,我可以用不同于完成项目的其他方式设置项目CRS吗?我想避免手动方式。 我通过添加以下代码解决了该问题: qgis.utils.iface.mapCanvas().setMapUnits(0) qgis.utils.iface.mapCanvas().refresh() 这会将单位设置为米,并且刻度会自动更新。0代表米,1代表英尺,2代表度,3代表未知数。
10 qgis  python 

2
通过PyQGIS添加工具栏?
通过教程,我学习了如何通过python将工具按钮添加到plugins-toolbar。现在,我想知道如何通过python添加带有工具栏按钮的完整工具栏。 有人可以提供示例代码吗?
10 pyqgis  toolbar 

3
使用Python获取QGIS中一项选定功能的数据吗?
我对python和QGIS有问题。是否可以使用python获取QGIS中一项选定功能的数据,如果可以,怎么办? 在我的示例中,我有一个图层,其中显示了一些用于城市规划的地块。如果现在选择一个宗地,我想获取该宗地的数据。此数据应通过python脚本读取,然后添加到自动创建的PDF文件中。最后一步很容易,但是我没有从所选功能中获取数据。 可以从活动层获取数据,但是也可以从该层的一个选定要素中获取数据吗? 希望您能够帮助我。问候
10 qgis  python  pyqgis 

4
将重叠爆炸到新的非重叠多边形?
给定多个以多种方式重叠的多边形,我想从这些要素中迭代导出所有不与其他多边形重叠的多边形。 该产品将具有许多不重叠的功能,将它们加在一起便构成了原始功能。 然后可以将乘积用作区域统计的输入,这比在每个多边形上迭代区域统计要快得多。 我一直试图在ArcPy中对此进行编码,但没有成功。 是否已经存在执行此操作的代码?


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.