地理信息系统

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


6
什么是变更检测?如何使用开源工具执行此类分析?
从维基百科页面: GIS(地理信息系统)的变更检测是一个过程,用于测量特定区域的属性在两个或多个时间段之间如何发生变化。变化检测通常涉及比较在不同时间拍摄的区域的航空照片或卫星图像。该过程最常与环境监测,自然资源管理或衡量城市发展相关 比较如何完成?用什么工具?我觉得说明不完整。或缺少某些东西。 在哪里或在哪本书中可以找到有关变更检测的更多信息? 我应该使用哪些工具来使用shapefile中的数据执行此类分析?(请仅开放源代码) 编辑 有关变更检测的一些论文(理论和技术) 变更检测技术(D. LU,E。BRONDI,ZIO和E. MORAN,2004,pdf) NDVI时间序列中的趋势变化检测:年际变化和方法学的影响Forkel,M.,Carvalhais,N.,Verbesselt,J.,Mahecha,MD,Neigh,C.,Reichstein,M.(2013)遥感5 (2013)5。-ISSN 2072-4292-第 2113-2144。 全球植被活动趋势的变化Jong,R.de,Verbesselt,J.,Zeileis,A.,Schaepman,ME(2013)遥感5(2013)3。-ISSN 2072-4292-第 1117年-1133年。 根据MODIS时间序列(2000-11),夏季海冰下降,温度升高和西伯利亚北极冻原植被变化之间的关系。Dutrieux,LP,Bartholomeus,HM,Herold,M.,Verbesselt,J.(2012)环境研究快报7 (2012)4。-ISSN 1748-9326-页。12 使用卫星图像时间序列进行近实时干扰检测Verbesselt,JP,Zeileis,A.,Herold,M.(2012)遥感环境123(2012)。-ISSN 0034-4257-第 98-108。全球绿化和褐变趋势变化:短期趋势对长期变化的贡献Jong,R. de,Verbesselt,J.,Schaepman,ME,Bruin,S. de(2012)全球变化生物学18 (2012)2。-ISSN 1354-1013-第 642-655。 物候变化检测,同时考虑卫星图像时间序列中的突变趋势和渐进趋势Verbesselt,J.,Hyndman,R.,Zeileis,A.,Culvenor,D.(2010)遥感环境114(2010)12。-ISSN 0034-4257-第 2970-2980。 检测卫星图像时间序列中的趋势和季节变化Verbesselt,J.,Hyndman,R.,Newnham,G.,Culvenor,D.(2010)遥感环境114(2010)1。-ISSN 0034-4257-第 106-115。 (我会在以后添加更多,就像发现更多著名论文一样)

7
使用QGIS识别字段中的重复属性?
我有一个包含数千个点的shapefile文件。它具有一个应该唯一的ID码字段。数据输入服务员有时会错误地键入ID,从而创建重复项。现在,我正在手动滚动字段以查找重复项。 使用搜索查询生成器还有另一种方法吗?

2
R将SF对象转换回SpatialPolygonsDataFrame
如何将sf对象转换回SpatialPolygonsDataFrame? 我正在运行R 3.4.0和SF版本0.4.1。该SF文件中提到的as_Spatial功能,但这不是提供我的安装(出于某种原因)。
27 r  sf 

2
使用QGIS 3将打印布局设置为纵向?
我无法使打印版式以纵向显示。它以横向打开,并且“布局/页面设置...”允许纵向显示,但布局画布仍以横向显示。 在2.18中,“组合”面板的第一项为“页面大小”,并带有一个横向/纵向控制框。 在3.0.0中,“布局”面板没有此功能,在“项目属性”面板中找不到它。 我可以在某处设置肖像吗?

3
外形不正确
我试图在新计算机上运行Shapely,并尝试从源代码安装它,使用其网页上的安装程序并使用pip install shapely。我尝试从库中进行的每个引用都会from shapely import *引发错误WindowsError: [Error 126] The specified module could not be found。在尝试利用Shapely并弄清楚如何使其正常运行时,是否有人有类似的经验? 引发的Python错误: WindowsError: [Error 126] The specified module could not be found 环境: Windows 7 64位 Python 2.7


6
在Ubuntu 16.04 Xenial上安装QGIS 2.14 Essen
我在Xenial上安装QGIS时遇到问题。 `debian`发布 当我尝试使用debian发行版(即没有ubuntugis依赖项)时,出现以下错误(请注意,没有xenial发行版,因此我wily在设置时使用了sources.list): ~$ sudo apt install qgis python-qgis qgis-plugin-grass Reading package lists... Done Building dependency tree Reading state information... Done Some packages could not be installed. This may mean that you have requested an impossible situation or if you are using the unstable distribution that some required packages have …

3
椭球是数学上的必要吗?
文献通常只是告诉我们,大地水准面太复杂而无法用数学方法来描述,因此我们使用不同的椭球体对其进行近似。 这些椭圆体在数学上是否必要,或者我们也可以定义从大地水准面模型到平面坐标的投影吗?

1
重力/霍夫模型工具
我正在寻找一种使用基于点的图层来模拟重力模型的方法。 我所有的点都分配有一个z值,该值越高,其“影响范围”就越大。这种影响与到中心的距离成反比。 这是一个典型的霍夫模型,每个点都是局部最大值,它们之间的谷表明它们之间影响区域的极限。 我尝试了Arcgis(IDW,成本分配,多项式插值)和QGIS(热图插件)中的几种算法,但没有发现任何对我有帮助的算法。我也找到了这个线程,但是对我不是很有帮助。 或者,如果有一种方法可以通过相应点的z值影响每个像元的大小,那么我也可以通过生成Voronoi图的方式来满足。

4
卸载所有ArcGIS产品的最简单方法是?
Windows计算机可能会一次或多或少地卸载任何和所有ArcGIS产品的最简单,最直接的方法是什么?也许可以编写脚本? 从“控制面板”中的“添加/删除程序”非常简单,但是当安装了多个产品时,这是单击-单击-单击的过程,等待产品1的卸载,这可能需要一些时间,然后单击-单击-单击以获取# 2,依此类推。对于一台机器这没关系,但是当有几台机器时,这确实很痛苦。


3
gdalwarp cutline和shapefile
我正在尝试根据shapefile几何形状裁剪我的Raster文件。我正在使用以下代码 gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif 但这会在shapefile几何图形范围之外产生黑色。我在这里举个例子。在我要剪切的第一个图像中。第二个图像是结果栅格,但我想要第三个栅格,其中几何图形的外部为空。 第一张图片:INPUT Raster 第二张图片:OUTPUT Raster 第三张图片:我想要这样的输出
26 gdal  gdalwarp 

2
计算UTM区域外的面畸变?
我的一位同事正在处理分布在两个UTM区域中的数据。大部分数据位于一个区域中,而一些异常值位于另一个区域中。他想知道如果这些离群值位于主UTM区域中,那么它们的区域失真会是什么。 是否知道已知特征在另一个UTM区域内的距离,可以使用公式来计算面畸变?

5
是否将栅格完全加载到numpy数组中?
我一直在尝试检查我在DEM栅格上的过滤器以进行模式识别,并且它总是导致缺少最后一行和最后一列(如..20)。我已经尝试了PIL库,图像加载。然后用numpy。输出是相同的。 我以为我的循环出了点问题,当检查数组中的值时(只是用ArcCatalog中的Identification选择像素),我意识到像素值未加载到数组中。 因此,只需打开,放入数组并保存数组中的图像即可: a=numpy.array(Image.open(inraster)) #raster is .tif Float32, size 561x253 newIm=Image.new(Im.mode, Im.size) Image.fromarray(a).save(outraster) 导致删除最后的行和列。抱歉,无法发布图片 任何人都可以帮助理解为什么?并建议一些解决方案? 编辑: 因此,在伙计们的帮助下,我成功地将小型栅格加载到numpy数组中,但是当图像较大时,我开始出现错误。我想这是关于numpy数组的限制,因此数组会像这样自动重塑或变形……所以ex: Traceback (most recent call last): File "<pyshell#36>", line 1, in <module> ima=numpy.array(inDs.GetRasterBand(1).ReadAsArray()) File "C:\Python25\lib\site-packages\osgeo\gdal.py", line 835, in ReadAsArray buf_xsize, buf_ysize, buf_obj ) File "C:\Python25\lib\site-packages\osgeo\gdal_array.py", line 140, in BandReadAsArray ar = numpy.reshape(ar, [buf_ysize,buf_xsize]) File …

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.