1

我正在构建一个使用 PolyLineZ (ESRI Shapefile) 数据并重写外围 Z 值的 Windows 窗体应用程序。最小和最大 Z 值由用户通过界面定义

我们以以下为例,假设最小值为 0,最大值为 10:

XY      Z
1,1     0
1,3     1
1,5     7
1,7     11*
1,10    10

需要对带有 11 的值进行插值,因为它不属于用户定义的范围。这显然是一个非常简化的例子。某些折线可能会丢失更多值。

我做了什么:

我研究过线性插值。看了一些 youtube 视频,我很难理解它。

我需要的:

来自任何语言的代码示例或线性/双线性/三线性插值背后理论的“英语”解释,以便我可以将其实现到我的程序中。我的数学技能不是最好的,所以我很难理解维基百科对它的定义。

我还假设线性插值是我需要研究的,

编辑:目前正在实施以下内容,如果我错了,请阻止我

我正在使用我认为是毕达哥拉斯理论类型的方法。我还没有让它捕获异常(即,确保左边的点实际上是离开的,确保列表没有超出范围等),这可能会在以后出现

internal static double calculateDistance(XYPoints a, XYPoints b)
{
    double xd = b.X - a.X;
    double yd = b.Y - a.Y;
    return Math.Sqrt(xd * xd + yd * yd);
}


for (var i = 0; i < polylinez.ZPoints.Count;i++)
{
    if (polylinez.ZPoints[i] > maxValue || (polylinez.ZPoints[i] < minValue))
    {
        //polylinez.ZPoints[i] = (((1 - polylinez.XYpoints[i].X) * polylinez.ZPoints[i - 1]) + (polylinez.XYpoints[i].X * polylinez.ZPoints[i + 1]));
        double prevdistance = calculateDistance(polylinez.XYpoints[i - 1], polylinez.XYpoints[i]);
        double nextdistance = calculateDistance(polylinez.XYpoints[i], polylinez.XYpoints[i + 1]);
        double fraction = prevdistance / nextdistance;
        double diffsBetweensZ = polylinez.ZPoints[i + 1] - polylinez.ZPoints[i - 1];
        Console.WriteLine(polylinez.ZPoints[i - 1] + (diffsBetweensZ * fraction));
    }
}

return polylinez;

它返回 9.12 作为上述示例表的答案。这听起来对我来说是正确的。我用互联网上的样本数据检查了我的距离计算器方法,它似乎在起作用。

4

1 回答 1

0

第一步,创建一个计算距离的例程:

internal static double calculateDistance(XYPoints a, XYPoints b)
{
    double xd = b.X - a.X;
    double yd = b.Y - a.Y;
    return Math.Sqrt(xd * xd + yd * yd);
} 

我将变量名称更改为更合乎逻辑的名称(我的变量名称不同)

//get distance frpm previous point to point in question
double prevdistance = calculateDistance(prevXYpoint, currentXYPoint);
//get distance frpm point in question to the next point
double nextdistance = calculateDistance(currentXYPoint, nextXYPoint);
//generate a ratio
double fraction = prevdistance / (nextdistance + prevdistance);
//find out the difference between the two known points
double diffsBetweensZ = nextZpointValue - prevZpointValue;
//interpolate!
double newZvalue = (prevZpointValue + (diffsBetweensZ * fraction));

我在几组数据上检查了这一点,这是我能找到的最准确的东西......让我大吃一惊的是,我在任何地方都找不到任何现有的代码来做到这一点。

于 2014-02-26T14:24:50.880 回答