3d工具算法c语言描述,C语言中的3D Sobel运算符算法

我目前正在努力用C制作3D Sobel边缘检测器(我还很陌生)。它不能完全按预期工作(突出显示实体3D对象中的非边缘)，我希望有人可以看到我出了问题的地方。 (对不起，这篇文章中的空格不多)

首先，im是已复制到tm的输入图像，每边的边界均为1像素。

我遍历图像：

for (z = im.zlo; z <= im.zhi; z++) {

for (y = im.ylo; y <= im.yhi; y++) {

for (x = im.xlo; x <= im.xhi; x++) {

我制作了一个数组，该数组将容纳x，y和z方向的更改，并循环通过3x3x3立方体：

int dxdydz[3] = {0, 0, 0};

for (a = -1; a < 2; a++) {

for (b = -1; b < 2; b++) {

for (c = -1; c < 2; c++) {

现在是这里的肉，有点棘手。我正在对我的Sobel运算符进行加权，以便如果您想象内核的一个2D表面，则该表面将是{{1,2,1}，{2,4,2}，{1,2,1}}。换句话说，核心像素的权重与其到中心像素的4连通度有关。

为此，我将e定义为3-(| a | + | b | + | c |)，以使其为0、1或2。内核将按3加权。每个像素处的^ e。

内核像素的符号将仅由a，b或c的符号确定。

int e = 3 - (abs(a) + abs(b) + abs(c));

现在，我通过将a，b和c打包成一个数组并从0-1-2循环来遍历。例如，当a为0时，我们不想向x添加任何值，因此我们使用if语句(深度为8级！)来排除该值。

int abc[3] = {a, b, c};

for (i = 0; i < 3; i++) {

if (abc[i] != 0) {

要添加的值应仅为该像素处的图像值乘以该像素处的内核值。 abc[i]仅为-1或1，而(int)pow(3, e)是接近中心的权重。

dxdydz[i] += abc[i]*(int)pow(3, e)*tm.u[z+a][y+b][x+c];

}

}

}

}

}

最后获取x，y和z中平方变化之和的平方根。

int mag2 = 0;

for (i = 0; i < 3; i++) {

mag2 += (int)pow(dxdydz[i], 2);

}

im.u[z][y][x] = (int)sqrt(mag2);

}

}

}

当然，我可以遍历图像并将3x3x3立方体乘以3D内核：

int kx[3][3][3] = {{{-1,-2,-1},{0,0,0},{1,2,1}},

{{-2,-4,-2},{0,0,0},{2,4,2}},

{{-1,-2,-1},{0,0,0},{1,2,1}}};

int ky[3][3][3] = {{{-1,-2,-1},{-2,-4,-2},{-1,-2,-1}},

{{0,0,0},{0,0,0},{0,0,0}},

{{1,2,1},{2,4,2},{1,2,1}}};

int kz[3][3][3] = {{{-1,0,1},{-2,0,2},{-1,0,1}},

{{-2,0,2},{-4,0,4},{-2,0,2}},

{{-1,0,1},{-1,0,1},{-1,0,1}}};

但是我认为循环方法要性感得多。