clear;
clc;
f1=fopen('5-15.raw','r');
data = fread(f1, 'uint16');
fclose(f1);
len = length(data);
min0 = min(data);
max0 = max(data);
interal = (max0 - min0)/10;
for i = 1:len
if data(i) < min0 + interal
data(i) = min0;
elseif data(i) < min0 + 2*interal
data(i) = min0 + 2*interal;
elseif data(i) < min0 + 3*interal
data(i) = min0 + 3*interal;
elseif data(i) < min0 + 4*interal
data(i) = min0 + 4*interal;
elseif data(i) < min0 + 5*interal
data(i) = min0 + 5*interal;
elseif data(i) < min0 + 6*interal
data(i) = min0 + 6*interal;
elseif data(i) < min0 + 7*interal
data(i) = min0 + 7*interal;
elseif data(i) < min0 + 8*interal
data(i) = min0 + 8*interal;
elseif data(i) < min0 + 9*interal
data(i) = min0 + 9*interal;
else
data(i) = max0;
end
end
k = len/(336*256);
image = reshape(data, 336, 256, k);
image = imrotate(image,270);
% imshow(image);
% 镜像变换
[Height,Width] = size(image); %获取原图像的高度和宽度
T1 = affine2d([-1 0 0;0 1 0;Width 0 1]); %构造空间变换结构T1，这里为水平镜像变换矩阵
image = imwarp(image,T1); %对原图像I进行水平镜像变换
figure,imshow(image/255,[]);
%热力图显示
cmap = colormap(jet(256));
rgb = ind2rgb(image,cmap);
