Peter的代码(感谢btw!)在您的图像是灰度的情况下工作。对于彩色图像,您需要做的是将所有字节读取为一个1D数组,因为我们不能使用fread读取3D矩阵中的数据(至少据我所知...)。之后,我们对其进行整形,使其成为3D矩阵。重塑矩阵形状后,请注意这是转置结果。因此,我们将不得不独立地对每个通道进行转置,或者您可以巧妙地组合使用flipdim和imrotate为我们进行3D转置。我会旋转图像,使其顺时针旋转90度,但这会使图像在各列上反射。然后,我将flipdim的第二个参数设置为2因为我想在各列之间进行镜像反射以获取原始图像。
因此,您实际上需要做的是:
row=576; col=768;
fin=fopen('m-001-1.raw','r');
I=fread(fin, col*row*3,'uint8=>uint8'); %// Read in as a single byte stream
I = reshape(I, [col row 3]); %// Reshape so that it's a 3D matrix - Note that this is column major
Ifinal = flipdim(imrotate(I, -90),2); % // The clever transpose
imshow(Ifinal);
fclose(fin); %// Close the file
我得到的是这张图片:
另外,您当然可以使用Peter的代码,但是您将使用它,以便一次将图像重建为一个彩色平面。换句话说,您可以执行以下操作:
row=576; col=768;
fin=fopen('m-001-1.raw','r');
I1=fread(fin, [col row],'uint8=>uint8'); %// Red channel
I2=fread(fin, [col row],'uint8=>uint8'); %// Green channel
I3=fread(fin, [col row],'uint8=>uint8'); %// Blue channel
I1 = I1.'; I2 = I2.'; I3 = I3.'; %// Transpose each channel separately
Ifinal = cat(3, I1, I2, I3); %// Create 3D matrix
imshow(Ifinal);
fclose(fin);
您将获得与上面完全相同的图像。