我有一个image_data
单元格数组,其中包含400张图像:
image_data = {img1, img2, ... img400}
每个图像(例如img1
)是一个128 x 128
像素值数组:
img1 = [0.2, 0.5, .....]
我想获取存储在image_data
中的所有图像的所有第i个像素,但是不使用for循环。
例如:
image_data = {img1, img2, img3}
img1 = [0.1, 0.2, 0.3]
img2 = [0.4, 0.5, 0.6]
img3 = [0.7, 0.8, 0.9]
我想快速计算出以下方差:
var([img1(1),img2(1),img3(1)]) = var([0.1,0.4,0.7])
var([img1(2),img2(2),img3(2)]) = var([0.2,0.5,0.8])
var([img1(3),img2(3),img3(3)]) = var([0.3,0.6,0.9])
请尝试避免使用循环,因为它很慢。
答案 0 :(得分:2)
第一步是通过沿第3维连接图像将笨拙的单元格数组转换为3维数组:
image_array = cat(3, image_data{:});
现在,我们要做的就是使用下标将每个图像中 ith 像素的向量转换为向量,并将其传递给var
:
pixel_var = var(image_array(1,1,:)); % var for pixel 1
或者一次完成整个图像并索引结果:
image_var = var(image_array, 0, 3); % take var along 3rd dimension
答案 1 :(得分:0)
只需添加到@beaker的解决方案中即可。
img1 = [0.1, 0.2, 0.3];
img2 = [0.4, 0.5, 0.6];
img3 = [0.7, 0.8, 0.9];
image_data = {img1, img2, img3};
image_array = cat(1, image_data{:}) % concatinate the cell array to 2-d array
pix_array = var(image_array,0,1) % apply var() to all columns
我明白了:
image_array =
0.1000 0.2000 0.3000
0.4000 0.5000 0.6000
0.7000 0.8000 0.9000
pix_array =
0.0900 0.0900 0.0900