0

我想在 scilab 中制作一个不重叠的 3D 矩阵,并使用 NxNxN 大小的立方体循环遍历它,将立方体中体素的平均值分配给中心体素。

转到第 2 页,看看我到底需要什么。 https://drive.google.com/file/d/0B5wCiQEnPJYZdGxCcmRBc0NHNGc/edit?usp=sharing

非常感谢你。

PS不用担心平均音量无效的问题。矩阵可以只是 100x100x100 或类似的。

4

1 回答 1

1

您描述的问题与您发布的论文不同。我按照您的描述进行了操作,但将其拆分为多行,以便您轻松适应您的要求。

这与有人在您的其他问题中发布的 MatLab 答案非常相似。我不知道其他答案中使用的任何与 SciLabmat2cell类似的东西。cellfun

clear; clc

K = 100
N = 5

mid = floor(N/2)

volume = rand(K, K, K)
cubeCount = floor( K / N )

for x=0:cubeCount-1
    for y=0:cubeCount-1
        for z=0:cubeCount-1

            // Get a cube of NxNxN size  
            cube = volume((1:N)+N*x, (1:N)+N*y, (1:N)+N*z);

            //Calculate the average value of the voxels in the cube
            avg = sum( cube ) / (N * N * N);

            // Assign it to the center voxel
            volume( N*x+mid+1, N*y+mid+1, N*z+mid+1 ) = avg
        end
     end
 end

disp( volume )
于 2013-10-28T21:27:35.063 回答