0

这是在 C 中。

我有两个二维数组 ArrayA 和 ArrayB，它们对相同的空间进行采样。B 对与 ArrayA 不同的属性进行采样的频率低于 ArrayA，因此它比 A 小。

只是想尝试定义一些变量： ArrayA：SizeAX 按 SizeAY，由 indexA 索引，用于位置 posAX，posAY ArrayB：SizeBX，按 SizeAY，由 indexB 索引，用于位置 posBX，posBY

ArrayA 和 ArrayB 是指向数组开头的指针，其中先存储 X 的行，然后递增 Y，然后存储下一行 X（Y=1）

所以我需要从给定的 indexA 设置 indexB，使其成为最近邻样本，以与 indexA 的值相关联。

这是我所在的位置（请更正任何错误！请注意，我从索引 0 开始）：如果 ArrayA 是 9x9 而 ArrayB 是 3x3： (posX,posY) posA 0,0; indexA = 0 posB 0,0; 指数B = 0

POSA 8,0; indexA = 8（第一行结束）posB 2,0；指数B = 2

位置 0,1; indexA = 9 posB 0,0; indexB = 0（仍然更接近底部点）

位置 0,3; indexA = 27 posB 0,1；指数B = 3

POSA 8,8; indexA = 80（最后一点）posB 2,2；指数B = 8

到目前为止我有： indexA = posAX + (posAY * SizeAX)

我尝试过的（当然失败了）： indexB = (int) (indexA * (SizeBX * SizeBY / (SizeAX * SizeAY)) + 0.5) // 似乎只适用于第一行和最后一个值.. 但是这显然是行不通的——但我很好奇它是如何将两者映射在一起的，但我会在修复它后研究它。

我无法访问 posAY 或 posAX，只能访问 indexA，但我应该能够使用 mod 和余数将其分解，对吧？还是有更有效更快的方法？一个

我也试过这个：

indexB = (posAY * SizeBY / SizeAY) * SizeBY + (posAX * SizeBX / SizeAX)

我认为问题是我需要将 X 和 Y 索引分开，然后再使用 SizeBX 和 SizeBY ？

一个额外的警告是 ArrayA 和 ArrayB 来自更大的数据集，它们都对更大的空间进行采样。由于矩形是任意的，因此 ArrayA 或 ArrayB 可能具有最接近矩形边界的点，从而导致其他问题，即最近邻居真正抓取的方式。我也不确定如何解决这个问题。

4

2 回答 2

0

这些是您选择的一些令人反感的名称，但至少它们已被定义。

我认为您想在 [0...1]x[0...1] 真实坐标空间中通过 (x,y) ；也就是说，A 的右下角应该从 B 的右下角获取值，同样对于左上角的中间等。这意味着您应该将数组中点的外边缘视为 0 宽度point 在 [0...1]x[0...1] 框的边缘采样值；即，如果您有一个 3x3 数组，则在 (0.5,0.5) 处有一个点，其余点位于边缘。

我假设您的 2d B 数组中包含真实值，因此插值是有意义的；因为数组的大小不同

这是从a开始的计划

indexA -> (posAX,posAY) -> (x,y) -> (fracBX,fracBY) ->(通过最近邻插值) 来自 ArrayB 的值

重要提示：(fracBX,fracBY) 是框 [0...SizeBX-1]x[0...SizeBY-1] 中的实值坐标。

让我们一步一步来。假设我理解你，值在内存中的顺序是左-> 右、上-> 下（英文阅读），就像标准 C 数组一样。然后：

unsigned posAX=indexA%SizeAX;
unsigned posAY=indexA/SizeAX;

现在，让我们映射到 (x,y)：

double x=posAX/(SizeAX-1.0); // we get double division when we subtract by 1.0
double y=posAY/(SizeAY-1.0);

现在，到 (fracBX,fracBY)，其中 0<=fracBX<=SizeBX 和 0<=fracBY<=SizeBY：

double fracBX=x*(SizeBX-1);
double fracBY=y*(SizeBY-1);

现在，在 B 数组中（最多 4 个）最近的积分点之间进行插值：

unsigned intBX=(unsigned)fracBX;
double aBX=fracBX-intBX;
unsigned intBY=(unsigned)fracBY;
double aBY=fracBY-intBY;
double *bv=ArrayB+(intBX*sizeBY)+intBY;
#define INTERP(alpha,v1,v2) ((1-alpha)*v1+alpha*v2)
#define INTERPI(alpha,i1,i2) (alpha>0 ? INTERP(alpha,bv[i1],bv[i2] : bv[i1])
double v0=INTERPI(aBX,0,1);
double value=fracBY>0 ? INTERP(aBY,v0,INTERPI(aBX,sizeBY,sizeBY+1)) : v0;

价值就是你的答案。需要检查小数位置 aBX 和 aBY 是否为 0，以防止访问数组末尾的值（即使这些值乘以 0 会被忽略，这也可能导致段错误）。或者，您可以通过分配比您需要的多 1 行/列来简化事情。

bv[0] 是 ArrayB[intBX][intBY]，bv[1] 是右一，bv[sizeBY] 是下一，bv[sizeBY+1] 是右下一。(aBX,aBY) 是 [0...1]x[0...1] 中的另一个点，但这次以 ArrayB 中的四个相邻点为界。

于 2009-12-18T21:06:42.590 回答
0

你的意思是你有两个不同间距的网格覆盖相同的区域，并在其中一个中给出一个网格点的索引，你想在第二个中找到最接近的一个。像这样：

int posBX = (int)floorf((float(posAX) / float(sizeAX - 1)) * float(sizeBX - 1) + 0.5f);
int posBY = (int)floorf((float(posAY) / float(sizeAY - 1)) * float(sizeBY - 1) + 0.5f);
int indexB = posBX + posBY * sizeBX;

获取posAX和posAY来自indexA：

posAX = indexA % sizeAX;
posAY = indexA / sizeAX;

做双线性插值：

float bx = (float(posAX) / float(sizeAX - 1)) * float(sizeBX - 1);
float by = (float(posAY) / float(sizeAY - 1)) * float(sizeBY - 1);
int x = min(int(floorf(bx)), sizeBX - 2); //x + 1 must be <= sizeBX - 1
int y = min(int(floorf(by)), sizeBY - 2); //y + 1 must be <= sizeBY - 1
float s = bx - float(x);
float t = by - float(y);
float v[4];
v[0] = arrayB[x + y * sizeBX];
v[1] = arrayB[x + 1 + y * sizeBX];
v[2] = arrayB[x + (y + 1) * sizeBX];
v[3] = arrayB[x + 1 + (y + 1) * sizeBX];
float result = (v[0] * (1.0f - s) + v[1] * s) * (1.0f - t) +
               (v[2] * (1.0f - s) + v[3] * s) * t;
于 2009-12-18T20:56:26.150 回答