最近在做CT配准效果的评估,通过比较未配准前两幅图像的差异图和配准后的差异图,很明显的呈现出配准效果。这里的差异图由灰度图转伪彩色图,然后两图融合实现。
方法:1 把CT图像转成BMP,这个自行搜索相关资料完成。
2 把BMP图像转成伪彩色图,MATLAB实现:
clear all;close all;clc;
I=imread('peizhun2.bmp');
imshow(I);
I=double(I);
[m,n]=size(I);
L=256;
for i=1:m
for j=1:n
if I(i,j)<L/4
R(i,j)=0;
G(i,j)=4*I(i,j);
B(i,j)=L;
else if I(i,j)<=L/2
R(i,j)=0;
G(i,j)=L;
B(i,j)=-4*I(i,j)+2*L;
else if I(i,j)<=3*L/4
R(i,j)=4*I(i,j)-2*L;
G(i,j)=L;
B(i,j)=0;
else
R(i,j)=L;
G(i,j)=-4*I(i,j)+4*L;
B(i,j)=0;
end
end
end
end
end
for i=1:m
for j=1:n
G2C(i,j,1)=R(i,j);
G2C(i,j,2)=G(i,j);
G2C(i,j,3)=B(i,j)