自动地图包：负值和变差拟合

我前两天开始使用自动映射封装自动化变差函数建模，并发挥它周围的（这样做手工是繁琐和不来太多的优势......虽然我觉得奇怪的是，automap似乎真的很喜欢Matérn函数来模拟变差函数）。我想到的是数据集中存在负值，而这些数据集中没有显示这些数值。自动地图包：负值和变差拟合

这里我只是粘在一起的人物：NR。 1是使用autoFitVariogram和autoKrige后的输出。它似乎没有负面价值。由于我不喜欢图形输出，我使用gstat执行普通克里格法，它使用autoFitVariogram（来自automap）的变差函数（编号2），它清楚地显示负降水区域。我将从autoKrige生成的地图导出到QGIS并进行了双重检查，并且还有许多负值（第3号;黑色区域）。

所以我的问题：

1）由autoKrige产生的克里格预测图显然产生了如下的方案，例如“如果var.pred < 0，则0”的地图。我可以看到这对于制作这样的地图非常有用，但是这对于实际数据有什么意义？我想用我的数据集进行进一步的分析，负降雨量值显得有些不现实。如果将这个“如果var.pred < 0然后0”方案也应用于数据，如果这会使其物理上合理吗？

2.）您可确定与广义最小二乘（GLS）作为自动地图包自身参数附加变差函数参数。但是，我想使用gstat提供的REML适配器（fit.variogram.reml）。我尝试用我的REML估计值替代variogram$var_model（这也产生了我的修改后的变差函数），但是当使用autoKrige执行克里金函数时，将再次使用原始变差函数模型。有没有一种方法可以通过包装提供的其他方式来适应模型？（在autoKrige手册中有一个参数'data_variogram'，但是这似乎做了其他事情）。

谢谢！

来源

2016-01-21 GeoEki

我敢肯定自动映射不限制预测至零以上。这可能是您正在使用的变差函数的乘积。另外，你现在不能使用fit.variogram.reml。 –

我不认为是这样的话;你似乎通过查看图表来得出结论，但请验证所生成的数字。图1在我看来并不包含负值。
autoKrige使用fit.variogram，而不能使用fit.variogram.reml;这是如何编码的，但随时可以修改它。

来源

2016-01-22 11:17:47

1.图1和图3来自相同的输出（来自'autoKrige'输出的第1号，第3号输出为ascii文件并输入到GIS中的相同地图，其中我用黑色给负值着色，而图2（不是基于'autoKrige'，而是'autofitVariogram'，所以相同的数据）暴露了负值的区域，如果图1也包含负值，那么配色方案或签名确实是误导性的。 2与0（'ifelse（prediction_spdf $ var1.pred <0,0，prediction_spdf $ var1.pred）'），它看起来与图1相同。 – GeoEki

2.这样一个可怜的人。不幸的是，我不能写我自己的包（但），但保罗可能会把这作为一个可能的未来实现:-) – GeoEki

自动地图包：负值和变差拟合

回答

相关问题