从堆栈中提取图像组

Question

我有一个图像堆栈，对于每个图像，我需要一个定义，由相邻的10 (5之前和5之后)定义的子堆栈。图像堆栈被定义为IM(m,:,:) (下面的详细信息)。

m = 0;
for j = 1:projections
    for k = 1:number_images
        u = sprintf('images_%03i_%05i.fits',j-1, k-1);
        IM(m,:,:) = fitsread(u);
    end
end

我想做的是(代码不起作用)

for k = 6:number_images
    rolling_interval_IM = [IM(k-5,:,:):IM(k+5,:,:)];

有更好的方法吗？

Answer 1

number_images = 60;
offset = 5;

for n=offset+1:1:number_images-offset
    rolling_interval_IM = squeeze(IM(n-offset:1:n+offset,:,:));        
end

您没有说明是否要存储每个滚动间隔，因此这只会覆盖前面的内容。

我认为把时间维度作为你的第一个维度有点不寻常。这通常是我的最后一次经历。