如何在MATLAB中读取RAW图像？

Question

我想在MATLAB中打开并阅读.raw图像。我的文件可以是downloaded here。我尝试了以下三个代码片段，但都没有给出预期的结果。

代码段＃1

    row=576;  col=768;
    fin=fopen('m-001-1.raw','r');
    I=fread(fin,row*col,'uint8=>uint8'); 
    Z=reshape(I,row,col);
    Z=Z';
    k=imshow(Z);

它显示了这张图片：

代码段＃2

    f=fopen('m-001-1.raw');
    a=fread(f);
    input_img = reshape(a,768, 576, 3);
    input_img = imrotate(input_img, -90);
    imwrite(input_img, 'm-001-1.jpg'); 

以.jpg格式保存空白（仅白色）图像。

代码段＃3

    id = fopen('m-001-1.raw', 'r');
    x = fread(id, [576,768], 'short');

当我使用imshow(x)时，此图片显示：

如何正确阅读此图片？

Answer 1

您的行/列大小相反。由于MATLAB数组是列主要的，并且光栅图像通常存储为行主要，因此您需要将图像作为[col row]矩阵读取，然后转置它。

row=576;  col=768;
fin=fopen('m-001-1.raw','r');
I=fread(fin, [col row],'uint8=>uint8'); 
Z=I';
k=imshow(Z)

图像复制正在发生，因为每行短768-576 = 192像素，因此每行都会逐渐关闭该行。在4行之后，你已经弥补了差异（4 * 192 = 768），所以你有4个图像复制。

Answer 2

彼得的代码（谢谢顺便说一句！）假设你的图像是灰度的。对于彩色图像，您需要做的是将所有字节作为单个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

我得到的是这张图片：

---

或者，您当然可以使用彼得的代码，但您可以使用它，以便您可以一次重建一个颜色平面的图像。换句话说，你可以这样做：

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);

---

您将获得与上面完全相同的图像。