matlab读取红外raw文件并聚类增强

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