9 克里金法的开源方法？ 我有一个点数据集，想要使用Krige，最好使用开源软件包。如果可能的话，我也想在过程中选择半变异函数模型以改善估计。 42 open-source-gis interpolation spatial-statistics kriging
1 选择IDW与Kriging插值进行DEM创建？ 我正在尝试使用非常规则地间隔约10m的点数据创建DEM。我要插值的区域是一个教育机构，有许多平坦的停车场和足球场，但仍然有一些相当陡峭的山丘，经常会停在停车场。由于这些已知的高原，我排除了Splining方法。但是，我仍然不确定在使用IDW和Kriging方法之间。尝试了这两种方法后，我看不出太大的区别，但经过一点研究仍未做出决定。 任何人都有智慧的话可以为我解决这个问题？ 23 dem interpolation kriging inverse-distance-weighted
3 投影航空照片时应使用哪种重采样技术？ 我正在对航空照片进行一些时间密集的投影，而且很好奇-哪种最佳重采样技术最适合用于航空照片？在ArcMap中，我的选项是NEAREST，BILINEAR，CUBIC和MAJORITY。 建议将最近邻和多数用于分类数据，而将三次卷积和双线性插值用于连续数据。 我很好奇，是否有用于投影航空照片的常用算法。我刚刚完成了使用最近邻投影的图像，看起来似乎不错，但是航拍照片不是分类数据，因此接下来我将尝试使用Bilinear。 编辑 我不认为航空照片与DEM或降水数据是连续数据一样，但是whuber指出它们是连续的，应该这样处理。再次感谢。 23 raster interpolation resampling aerial-photography
3 使用GDAL在Python中平滑/插值栅格？ 我正在用Python开发，并使用OSGEO的GDAL来处理栅格和shapefile并与之交互。 我想获取一个具有点特征的shapefile并将其插入到表面栅格中。现在，我使用的是'RasterizeLayer'方法，该方法将点要素中的值刻录到栅格中（使用所有nodata值进行设置），但将所有未触及的像素保留为'nodata'值。因此，我剩下的是棋盘式栅格。 使用RasterizeLayer之后，我得到了什么： 我想要的最终产品是： 我相信我要寻找的功能在arcgisscripting导入中称为'Spline_sa（）'。 GDAL是否具有类似的功能，或者有其他方法来获得所需的输出？ 20 python raster gdal interpolation
2 克里格插值的最小样本数 我得到了一些带有标本数量的数据，并要求使用克里金法对其进行插值。 经过一些调查，似乎克里金法结果（在ArcGIS Geostatistical Analyst中使用默认参数执行）不令人满意。内插值与测量值（尤其是顶部值）相差很大，并且表面看起来不可靠。这是图片： 我想主要问题是样本数量不足。 我们应该使用多少点来获得可靠的结果？ 也许克里金法不适用于这种多样化的价值？ 19 arcgis-10.1 interpolation kriging geostatistical-analyst sampling
1 哪种插值技术适合小湖的测深？ 这是一个约13公顷的湖泊，在约10条样线中有81个采样深度点： 以前在ARCgis中有一个50公顷的湖泊和大约100个深度采样点，我使用该TopotoRaster工具获得了不错的输出。但是，QGIS或开源软件似乎不存在这样的精确等效项。 在QGIS中，我使用Raster Interpolation插件中的TIN方法获得了下面的插值栅格。但是，这似乎是一个幸运的结果，当我尝试重复时，尽管没有真正改变任何参数，但我得到了不同的结果（不令人满意）。 这些比城市湖泊更多的城市（构造）水箱，因此，它们的河床轮廓相当规则，因此采样点密度低。一些文章似乎建议使用逆距离加权（IDW）作为最合适的技术，但是我似乎获得了最差的结果（可能是由于低点密度）。 考虑到湖泊的大小，其规则轮廓和点采样密度，可以采用任何启发式方法来得出哪种插值方法合适？（在克里金，IDW，双线性，三次卷积，TIN或样条之间） 还是总是反复试验？ 16 qgis raster interpolation bathymetry lakes
7 寻找插值工具以生成剖面的河流测深法（地形） 按照标题，我正在寻找一种可以从横截面产生河流测深的工具。尽管可能知道这种工具的人可能不需要任何解释，为什么这种插值有点复杂，但是对于那些可能知道并将其用于其他目的的人，这里有一些解释。 为了接收带有河床测深法的河谷DEM，我需要在要创建测深法的原始DEM上用栅栏围住/标记。基本上，在湿河床的整个区域，TLS都没有机会测量底面（测深法）。河流的形状各不相同，有时更曲折，有时更平直，并且底部的深度及其对齐方式变化很大，这使得插值更加复杂，因为它必须遵循河流的中心线并进行插值垂直和水平。 尽管我确实遇到了一些问题，但Venkatesh Marwade博士已经准备了对ArcGIS 9.2的此类扩展，并且只是想知道是否有人会知道可以产生类似内容的软件或代码。 （编辑：请参阅本教程中的图形。） 15 interpolation bathymetry
3 如何在PostGIS中内插GPS位置 我每五秒钟就有一张PostGIS GPS位置表： 2011-01-01 00:00:05, POINT(x1,y1) 2011-01-01 00:00:10, POINT(x2,y2) 2011-01-01 00:00:15, POINT(x3,y3) ... 我正在寻找一个查询，该查询将每秒返回值（时间戳和点）。可以假定点是通过直线连接的。 我专门在寻找一种在数据库内部执行此操作的方法，而不是通过编写一些外部脚本来实现。 13 postgis sql interpolation
3 租金价格建模-使用哪种插值方法？ 我有约140万个家庭的全国数据集。那里我有关于租金，面积（房间数和平方米）以及每个家庭的一些其他特征的信息。 我想使用此数据来创建整个国家的租金价格表，并使用此信息作为对拥有或没有租金信息的剩余约150万个家庭的价值估算的代理。 这里有几个问题： 这样的方法是否完全适合这种问题？ 哪种插值方法最适合在此使用？ 另外，是否可以考虑有关家庭规模的信息？ 我使用的是ArcInfo许可的ArcGIS 9.3。 13 arcgis-desktop arcgis-9.3 interpolation modelling
2 PostGIS插值/三角剖分选项 我正在寻找一个开源的工作流程来自动化DEM的构建。我们有一系列属于LIDAR数据集的站点。我们想要为每个站点创建站点特定的DEM，并且正在寻求使过程自动化。 到目前为止，我们已经实现了自动化： 数据加载到PostGIS中（激光雷达和站点点） 创建站点边界（st_buffer和st_envelope的组合） 现在，我们正在寻找在每个站点上插值数据子集并将其导出到Surfer7网格文件的选项。 当前，我们正在使用QuantumGIS查询PostGIS数据库中的数据，导出到csv文件，然后在Surfer7中手动导入和网格化数据。希望我们也可以使它自动化。 因此，对于此用例，我们希望将激光雷达数据三角剖分为DEM。我们可能想到的其他潜在情况涉及不同的插值方法-因此，如果有用于反距离和克里金法的选项-我们很感兴趣！ 对我们来说，这是一个非常学习的练习-我们正在逐步进行中！ 12 postgis dem interpolation lidar
4 用Python在栅格上对点数据进行双线性插值？ 我有一个栅格，想要对它进行点插值。这是我的位置： from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array = source.GetRasterBand(1).ReadAsArray() # Close raster source = None # Compute mid-point grid spacings ax = array([gt[0] + ix*gt[1] + gt[1]/2.0 for ix in range(nx)]) ay = … 12 python algorithm interpolation
1 您如何使用GRASS的v.kernel？ 我对如何使用GRASS的v.kernel感到困惑。 我的向量层约为250万点。我想使用v.kernel制作一个热图以显示浓度，因为我有可变实例，它们的点重叠，有时重叠很大。 我已经在GRASS中获得了此矢量层，并且显示得很好。 我已经尝试根据我在这里和其他论坛上看到的内容使用GRASS的v.kernel命令，除了输出仅是一个粉红色正方形的栅格外，我无法执行任何其他操作。 这是我正在使用的命令： v.kernel --verbose input=master_grass7 output=master_grass7a_heatmap stddeviation=.0001 我已经将标准差从1000000更改为.000001，并且没有任何效果。 我反复阅读了v.kernel文档，但并不太了解它的含义。至少，这些说明是基于深奥的概念，没有实际意义。我也检查了源代码，但我也不是很了解。是的，我可以阅读C。问题在于它取决于GRASS GIS中其他地方定义的许多内容。 我也做了很多Google搜索，但是找不到完整的指南。我所得到的只是v.kernel文档/手册页的零散副本，或者显然是让它正常工作的人。 我还检查了内核密度估计（KDE）的概念，即使这样，我也看不到如何使用v.kernel命令。该命令似乎是对KDE的特定解释；它的开关似乎与通用KDE概念不太吻合。 回到这里的主要问题：不熟悉GRASS产品开发的人如何使用v.kernel命令？有没有简单的语言翻译？ 12 grass interpolation density
3 R或ArcGIS中的时空插值？ 我正在尝试使用ArcGIS 9.3中的反加权距离工具从多个点计算平均降雨量。 我的问题是：每个点都有自己的时间序列，因此插值过程应该能够进行所有年份（可以说是迭代的一种）。 以下是样本属性表： ID X Y Name Rain1990 Rain1991 Rain1992 Rain1993 .... Rain2010 1 xx1 yy1 AA 1210 1189 1863 1269 ...... 2 xx2 yy2 BB 1492 1502 2187 1923 ...... ...... 有人可以告诉我该怎么做吗？ 编辑1：我最终使用C ++代码完成了此操作，这需要ArcGIS遮罩网格，数据文件和所有点的位置。 编辑2：我最近使用R来执行此插值任务。您可以使用hydroTSM，gstat或spacetime包。以下是几个示例链接： http://spatial-analyst.net/wiki/index.php?title=Spatial_interpolation_exercises_%28NL%29 http://www.geostat-course.org/Topic_Bivand_2012 编辑3：在下面为以后的读者添加了一个工作示例 12 arcgis-desktop r interpolation spatial-statistics time-series
4 如何正确插入温度？ 我试图对年平均温度进行插值以产生一个“逼真的”表面。在QGIS中，我使用了Raster-Interpolation-Interpolation。TIN和IDW两种方法均未提供“真实的”表面（例如，与地图集中的良好地图相比）。 IDW（因子3）： TIN（还显示插值点）： 有什么提示如何获得“更好，更现实”的插值？ 11 qgis interpolation climate