matlab读取raw文件

clc;clear all;
id = fopen('file.raw','r','b');
A = fread(id,'uint16');
imgsize = size(A);%读出的图像的size为n*1大小
rows = 512;%根据mhd文件的图像大小设定
clos =512;
nums = imgsize(1)/rows/clos; %计算得到图像张数
I = reshape(A,[rows,clos,nums]);
single_image = reshape(I(:,:,2),[rows,clos]);%显示第3张
B=uint16(single_image);%转化为unit8显示
B=B';%求图像
imtool(B);