本质上它被称为二维多元插值。 我有一组位置 (x,y) 和相应的数据,例如:z=颜色值。 对于某些 (x,y),z 列缺少数据。 我想使用最近的两个已知(x1,y1,z1)和(x2,y2,z2)来插值(x3,y3,z3=?)。
理论上这个维基/双线性插值是我喜欢得到的。
JMP有办法做到这一点吗?
另外,我需要首先按另外两列 A、B 对 X、Y 数据进行分组。例子。 COl A、B 代表文件名和帧号。 X,Y 是像素在帧中的位置,Z 是像素颜色。
我想对每个文件中每个帧中缺失像素的 Z 值进行插值。
样本数据集
| 文件 | 框架 | X | 是 | 颜色 |
| 蚂蚁人 | 1 | 0 | 0 | 255 |
| 蚂蚁人 | 1 | 5 | 5 | 195 |
| 蚂蚁人 | 1 | 3 | 2 | ?? |
| 蚂蚁人 | 2 | 0 | 0 | 255 |
| 蝙蝠侠 | 3 | 0 | 0 | 255 |
| 蝙蝠侠 | 3 | 3 | 2 | ?? |
| 蝙蝠侠 | 3 | 5 | 5 | 100 |
如果我理解正确,双线性插值将需要四个点(而不是两个)。 分组列的处理会很混乱,因此对于一组,您可以执行以下操作:
Names Default To Here( 1 );
// Example table: Four complete observations, and one missing 'Colour' to be
// obtained by bilinear interpolation
dt = New Table( "Antman_1",
Add Rows( 5 ),
New Column( "X",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Values( [0, 5, 3, 5, 0] )
),
New Column( "Y",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Values( [0, 5, 2, 0, 5] ),
Set Display Width( 92 )
),
New Column( "Colour",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Values( [255, 195, ., 100, 150] )
)
);
// Using the notation at: https://en.wikipedia.org/wiki/Bilinear_interpolation
x1 = dt[1, "X"];
x2 = dt[2, "X"];
y1 = dt[1, "Y"];
y2 = dt[2, "Y"];
q11 = dt[1, "Colour"];
q21 = dt[4, "Colour"];
q12 = dt[5, "Colour"];
q22 = dt[2, "Colour"];
x = dt[3, "X"];
y = dt[3, "Y"];
// Interpolate in the X direction
r1 = Interpolate(x, x1, q11, x2, q21);
r2 = Interpolate(x, x1, q12, x2, q22);
// Interpolate in the Y direction
p = Interpolate(y, y1, r1, y2, r2);
Wait(3);
dt[3, "Colour"] = p;
这不太可能是“最佳”方法,但为了做得更好,需要一些示例数据(在 JMP 表中)。 您还需要决定如何处理边缘情况。