如何使用MATLAB将图像分割成多个部分

Question

我获得了一个大小为[17935 10968]像素的卫星图像，我想平均切割图像并在各个部分上处理所需的算法（例如：我需要将上面的像素范围切成4个等分）。

如何在不丢失中间像素的情况下分割图像？我的要求类似于（1到5600和5601到末尾像素）。

有人知道如何在MATLAB中分割这么大的图像吗？

- ravi raja

---

那么问题是处理大数据出了问题还是切割图像的方式有误？ - Ray

@Ray 我都有了.. 有什么想法吗？ - ravi raja

Answer 1

方法一

如果您拥有图像处理工具箱，则这是首选和最有效的方法。它利用了非常实用的blockproc函数，该函数专门设计用于以块为单位处理大型图像。例如，当您的图像无法均匀分成相同大小的块时，它会进行填充，并将块处理的结果连接成一个结果矩阵。

最好查看官方文档，但以下是在您的情况下它可能会是什么样子：

vSize = [17935 10968];
imBig = rand([vSize 3]);
nParts = [2 2]; %means divide into 4 parts, 2 horizontal, 2 vertical

blockproc(imBig, ceil(vSize ./ nParts), @yourAlgorithm);


function res = yourAlgorithm(blockStruct)
   %do your processing of the block here and 
   %optionally return a result in 'res'

   %for example, just return the RGB vector of the first pixel
   res = blockStruct.data(1,1,:);
end

方法二

如果您没有图像处理工具箱，可以使用mat2cell函数。首先确定所需的块大小，然后获取包含不同块的单元数组。对于这样大的图像，速度和内存可能成为问题。此代码来自此 Matlab Central答案。

vSize = [17935 10968];
imBig = rand([vSize 3]);
nParts = [2 2]; %means divide into 4 parts, 2 horizontal, 2 vertical

%figure out the size of "regular" block and the last block

vRegBlockSize = ceil(vSize ./ nParts);
vLastBlockSize = vSize - vRegBlockSize .* (nParts - 1);

%put the sizes into a vector
vSplitR = [vRegBlockSize(1)*ones(1,nParts(1)-1), vLastBlockSize(1)];
vSplitC = [vRegBlockSize(2)*ones(1,nParts(2)-1), vLastBlockSize(2)];

%split the image
C = mat2cell(imBig, vSplitR, vSplitC, 3);

%access RGB pixel (x,y) in top left {1,1} block
p = C{1,1}(x, y, :);

- stav

Answer 2

upperLeft = theImage(1:5600, 1:5600, :);
upperRight = theImage(1:5600, 5601:end, :);
lowerLeft = theImage(5601:end, 1:5600, :);
lowerLeft = theImage(5601:end, 1:5601:end, :);

- user3128961

Answer 3

你可以使用 reshape 将图像分成四个矩阵：

A=reshape(Img, 17935 , 10968/4,[]);

然后处理A(:,:,1)，以此类推...

- bla

---

谢谢您的回答，非常有帮助；但是是否有任何方法可以获得RGB合并图像的输出，因为我正在为一个输入RGB图像（17935x10968x3）获取12张图片。 - ravi raja

Answer 4

使用以下代码将图像分割为4个不同的图像：

A=reshape(Img, 17935 , 10968/4, 3, []);

那么A(:,:,:,1)就是第一张图片。

- Lokesh A. R.

Answer 5

假设A是你的17935x10968x3矩阵，我认为你可以这样做：

B = reshape(A, 17935, 10968 / 4, 4, 3);

这样，最后一个维度仍然代表RGB。唯一的区别是它变成了4-D数组。

- Ray

---

我需要按照X和Y维度将其分割成4个相等的正方形部分。谢谢帮忙。 - ravi raja

@ravi 那么即使“17935”无法使它们相等，也要有四个矩形吗？ - Ray

我的意思是，如果图像像A4纸一样大，我需要在中心剪开变成两半，然后再将其中的一半再次切成两半。这样我就可以得到4个部分。 - ravi raja