Questions tagged «interpolation»

我想知道是否有人可以澄清使用ArcGIS IDW空间插值时是否应将点数据投影到等距投影？ 我正在研究来自北美西部的数据集，其跨度约为30度。数据当前为纬度/经度（NAD83）。 在插值过程中计算到样本点的距离时，ArcGIS是“即时投影”还是以某种方式调整了纬度？或者我应该在保留距离的投影中提供所有内容吗？

11 arcgis-desktop  coordinate-system  interpolation 

我使用包gstat，variogram()function 在R中制作了半变异函数。我想看看我的模型残差中是否存在空间自相关（使用glmm，在相距几公里至900公里的站点之间，物种的丰富度是栖息地的函数）。 我的单位是千米，所以我的解释是，距离只有100多公里，直到空间自相关不再是“问题”。我想知道是否有人可以解释为什么金块看起来这么高？这是否意味着即使在相似的位置，仍然存在相对较高的差异？或者，这个波浪形的变异函数是否意味着我应该调整滞后次数和滞后距离，直到获得更典型的形状？ 为了进一步研究，我还使用了variog()geoR包中的函数，并使用breaks=seq(0,100,10)来尝试查看更近的距离（使用相同的点和相同的模型残差）。这表明最接近的点有更多不同，这也没有意义。也许这表明不存在空间自相关，并且我的模型已经对此进行了说明。 我找到了这个极好的资料来源，“没有眼泪的地静止器”，并且在第51页上有一些关于拟合方差图的好的建议。根据这个建议，我的第一个似乎具有正确的范围。因此，这回到了第一个问题-我该如何解释？

10 r  interpolation  spatial-statistics  autocorrelation  variogram 

这是我在此的头一篇博文。 我从土壤深度调查中得出了几点，我需要对整个区域的深度进行插值。最合适的技术是什么？使用该技术时我需要考虑什么？指向教程的指针会很棒！ 谢谢

10 interpolation 

我正在寻找一种网格分类数据的方法。我从航海图和领域表中提取了一些点集合，这些点指定了海洋底部表面的性质。这些数据是分类数据，而不是数字数据，它们不是定期甚至随机抽样的。创建航海图以帮助导航和锚定；它们不是为绘制栖息地而创建的。因此，在岸边附近会有更多的探测声，在这些探测声中，相对浅的深度可能会对航行造成危险，并且船只倾向于停泊。离海岸较远，那里的深度足以进行导航和锚固是不切实际的，因此发声的频率要低得多。 有没有其他人试图从航海图创建网格化的底图？ 我查看了蒂森（Vornoi）多边形，但是沿海岸的测深集中导致沿海岸，远海的大型多边形以及在远处呈扇形延伸的长饼形多边形之间的精细“蜂窝”。使用最近邻居的网格产生的结果几乎相同。 我需要一种方法来限制浅水，近岸点的影响-一种方法来限制那些长的饼形多边形。在更深的水域中，我不认为底部的性质将是近岸底部的延续。我开始沿着两条线思考-都使用深度。一种方法是使用网格像元与相邻点之间的深度差来加权“最近”邻居的选择。另一个是取消选择深度大于某些指定公差的相邻点。或者，也许不是预先指定的公差，我可以对深度范围进行分类，然后将相邻点的选择限制为相同深度范围或分类中的点。 关于如何实现这两个选项中的任何一个想法？ 自从在其他论坛上与同事交谈以来，我一直在研究其他两种方法。第一种方法是使用障碍物（深度为100m的轮廓）来限制近岸数据的影响。这种方法的挑战在于，可以使用障碍的任何ESRI插值例程都旨在处理连续数据而不是不连续数据。在创建Thiessen多边形之前，我可以使用障碍物将点分解为较浅的近岸点和较深的点。但是，由于ArcGIS为矩形区域而不是复杂区域创建了Thiessen多边形，因此我预计边缘效果会很猖ramp。 由一些同事建议的第二种方法是克里金法。最初，我一度不理会克里金法，因为我只考虑了连续数据。克里金法的挑战在于它也不是为分类数据而设计的。现在，我正在研究表面深度和性质的协同克里金法，但是任何类型的克里金法都将涉及使用整数数值表示表面性质。之后，必须将所得的浮点数字代码还原为原始整数编码。不漂亮。 任何人都可以建议其他路线吗？（也许可以使用地形分析。例如，比休止角陡的斜坡不可能是沉积物。我正在寻找更简单的东西，无论如何，我没有足够的空间分辨率的数据。） 问候，

10 raster  interpolation  kriging  voronoi-thiessen 

我正在尝试在沿海地区进行插值，将数条河流中的数据点排入海湾和半岛的相对两侧。 我创建了海岸线图层以将其用作障碍，但是无法弄清楚如何在ArcGIS 10中执行此操作。 我没有在任何克里金工具中看到障碍物。 我试着只是进行克里金法，然后将其裁剪到土地上，但是我得到了不切实际的结果，数据点靠在一起但被土地分开会造成问题。 我在9.3版的地理处理脚本中看到了用于执行此操作的在线文档，但对于10版则没有任何相似之处。

10 arcgis-desktop  arcgis-10.0  interpolation  kriging  barriers 

当我想执行克里金法时，它有时只能工作，这取决于我在数据表中使用的值。作为krige函数的结果，我得到了var1.pred: NA NA NA ...and var1.var: NA NA NA ...（但仅当我在数据表中使用“错误”值时）。 例如： 它的工作原理总是（到目前为止）时，我只用10个值 它的工作原理，当我用50个值，但只有特定的一些 当我使用50个值和“错误”值时，它不起作用 当我使用25个值和前面提到的“错误”值时，它起作用 我不明白为什么它有时起作用而有时却不起作用。奇怪的是，当我添加Zwiesel;49.02999878;13.22999954;2.2到数据表中时，当我使用少于20个值时，它就起作用了；但是当我使用了50个以上的值时，它就不起作用了。 我的错误在哪里？ myWeatherTable.csv： Place;Latitude;Longitude;Temperature Aachen;50.77999878;6.09999990;3 Abbikenhausen;53.52999878;8.00000000;7.9 Adelbach;49.04000092;9.76000023;3.1 Adendorf;51.61999893;11.69999981;1.9 Alberzell;48.45999908;11.34000015;4.6 ... ... 我的代码执行克里金插值 WeatherData <- read.csv(file="myWeatherTable", header = TRUE, sep ";") coordinates(WeatherData) = ~Longitude + Latitude vario <- variogram(log(Temperature) ~1, WeatherData) vario.fit <- fit.variogram(vario, vgm("Sph")) min_lon …

10 r  interpolation  kriging  gstat 

我想创建风图。这里有关于ArcGIS和QGIS的介绍。由于我无权访问Spatial Analyst，因此尝试了QGIS版本。但是，我不知道如何对度数据进行插值。 例如，风向可以在一个点处为359°，在第二点处为1°。此点之间的值可能应为359°，0°或1°。 由于插值不知道它是度数据，因此看起来插值在180°左右。 如何使用QGIS插值学位数据？

9 qgis  interpolation 

在的文档中，sp我看到coords该SpatialPolygons函数的参数仅接受2列。无法创建3维多边形吗？我正在通过使用package gstat通过将时间添加为a的第三个空间维度来进行时空克里金法SpatialPointsDataFrame。我想将其扩展为以多边形为块的克里金法。有什么办法吗？

9 r  interpolation  spatial-statistics  kriging  gstat