读文件
clc;clear all;
id = fopen('file.raw','r');
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);
循环读文件,判断数字是几位数字
for i = 0:730
temp = num2str(i);
if(length(temp) == 1)
str1 = ['000',temp];
elseif(length(temp) == 2)
str1 = ['00',temp];
elseif(length(temp) == 3)
str1 = ['0',temp];
end
root = [obj_root,'\',str1,'.raw'];
fileID = fopen(root,'r');
A = fread(fileID,'uint16');
fclose(fileID);
rows = Slice;%根据mhd文件的图像大小设定
clos = Detector_num;
imgsize = size(A,1);%读出的图像的size为n*1大小
single_image = reshape(A,[rows,clos]);
% % single_image = reshape(A,[clos,rows]); wrong
% % imtool(single_image);
obj_I_120kv(i+1,:,:) = single_image;
end