地理信息系统

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

2
克里格插值的最小样本数
我得到了一些带有标本数量的数据,并要求使用克里金法对其进行插值。 经过一些调查,似乎克里金法结果(在ArcGIS Geostatistical Analyst中使用默认参数执行)不令人满意。内插值与测量值(尤其是顶部值)相差很大,并且表面看起来不可靠。这是图片: 我想主要问题是样本数量不足。 我们应该使用多少点来获得可靠的结果? 也许克里金法不适用于这种多样化的价值?





6
如何从代理后面访问QGIS插件存储库?
我一直在尝试通过使用此论坛和其他人来帮助自己QGIS。我的大多数问题都可以通过下载第三方存储库中可用的插件来回答,但是由于某些原因,我无法与这些插件建立连接。任何人都知道为什么会这样吗?

2
使用QGIS裁剪栅格图层?
如何在QGIS中裁剪栅格图层? 我想从图层的一部分制作轮廓,包括整个图层将创建比我需要的数据更多的数据,并且减慢了渲染速度。 我发现可以通过使用栅格计算器并指定新的图层范围来做到这一点,但是我如何知道要为行和列输入哪些值? 我也尝试在Photoshop中裁剪该geotiff,但似乎丢失了坐标数据。 有什么方法可以从栅格图层的一部分创建轮廓而不必先裁剪吗?

2
如何计算PostGIS中两条线相交的角度?
我想计算在PostGIS中相交的两条线之间的角度。 PostGIS中角度计算的起点似乎是ST_Azimuth-但这需要点作为输入。我的第一个想法是采用相交线的端点,并对这些端点执行方位角计算。这还不够好,因为大多数线要素不是笔直的,而且我对相交处的角度感兴趣。所以我想出了一个嵌套的操作,它通过以下步骤进行: 识别两个线要素表之间的所有交点。 在相交点周围创建一个非常小的缓冲区 确定线要素与缓冲区外部相交的点(如果有多个点,则取第一个点-我真的只对角度是接近0、90或180度感兴趣) 计算这两点的ST_Azimuth。 完整的SQL有点想在这里发布,但是如果您有兴趣的话,我可以在这里进行介绍。(顺便说一句,有没有比继承WITH语句下的所有字段更好的方法了?) 结果看起来不正确,所以我显然做错了: 编辑我在EPSG:3785中重新计算,结果有些不同,但仍然不正确: 我的问题是此过程中的缺陷在哪里。我是否误解了ST_Azimuth的功能?是否存在CRS问题?还有其他东西吗?也许有很多简单的方法可以做到这一点?
19 postgis  geometry 


5
如何创建PostGIS模板?
我正在按照教程进行安装。但是,当它说 启动pgAdmin并基于PostGIS模板创建一个新数据库 我在哪里可以得到这个Postgis模板?我知道我已经使用自制软件安装了PostGIS。但是我找不到模板。我现在正在运行macOS。
19 postgis  homebrew 



1
莫兰I的强大替代品
Moran的I(一种衡量空间自相关的指标)并不是一个特别可靠的统计数据(它可能对空间数据属性的偏斜分布敏感)。 有哪些更健壮的技术来测量空间自相关?我对像R这样的脚本语言中易于使用/可实现的解决方案特别感兴趣。如果解决方案适用于独特的情况/数据分布,请在答案中指定。 编辑:我正在用一些示例扩展问题(以回应对原始问题的评论/答案) 有人建议,置换技术(使用蒙特卡洛程序生成Moran's I采样分布)提供了一种可靠的解决方案。我的理解是,这种测试消除了对Moran's I分布进行任何假设的需要(假设测试统计量可以受数据集的空间结构影响),但是,我看不到置换技术如何正常地校正分布式属性数据。我提供两个示例:一个示例说明了偏斜数据对局部Moran I的统计影响,另一个示例对整体Moran I的影响-即使在置换测试下也是如此。 我将使用Zhang 等。的(2008)分析为第一个例子。在他们的论文中,他们使用置换测试(9999个模拟)显示了属性数据分布对局部Moran I的影响。我已经使用原始数据(左图)和对相同数据的对数转换(右图)在GeoDa中重现了作者针对铅(Pb)浓度(在5%置信水平)的热点结果。还显示了原始Pb和对数转换后的Pb浓度的箱线图。在这里,转换数据时,重要热点的数量几乎翻了一番。此示例表明,即使使用蒙特卡洛技术,本地统计信息对属性数据的分布也很敏感! 第二个示例(模拟数据)演示了偏斜数据可能会对全局Moran I产生的影响,即使使用置换测试也是如此。R中的示例如下: library(spdep) library(maptools) NC <- readShapePoly(system.file("etc/shapes/sids.shp", package="spdep")[1],ID="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) rn <- sapply(slot(NC, "polygons"), function(x) slot(x, "ID")) NB <- read.gal(system.file("etc/weights/ncCR85.gal", package="spdep")[1], region.id=rn) n <- length(NB) set.seed(4956) x.norm <- rnorm(n) rho <- 0.3 # autoregressive parameter W <- nb2listw(NB) # …

3
如何计算给定大地纬度下的地球半径?
(我看到Wikipedia上有一个方程式可以完全满足我的要求,但没有参考文献。我无法确认该方程式的有效性!) 我已经了解了地心纬度和大地纬度之间的区别。 假定已知的半大半圆a和半小半圆b半径。如何在给定的大地纬度下计算半径? 我需要某种专家确认(派生,派生链接,专家确认,解释等)。

6
如何在OpenLayers / Leaflet中创建动画群集标记?
由于Openlayers似乎无法比拟的平滑过渡,我想使用Leaflet。我特别需要能够将标记与自定义样式的标记组聚类的功能。 可以在Redfin上找到证明我正在努力实现的最佳示例。 请注意,将鼠标悬停在群集上可提供流畅的动画,并通过单击标记将群集分割成动画效果。当达到单个点不与给定公差重叠的缩放级别时,实际标记将以动画方式拆分到其位置。 我真的很想将类似的效果整合到其中,并希望就最佳方法或您可能知道的其他示例提供一些指导。 注意: Redfin似乎正在使用Flash,我想远离它。在理想的情况下,我想尽可能通过Javascript实现此功能,但认为可能需要通过HTML5 / canvas完成。 也许使用p。-Raphael.js或d3.js 有任何想法吗? 顺便说一句-找到了这个例子,但它似乎特别适用于Google Maps。

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.