MATLAB中对比两幅图像各种直方图的程序-

1、MATLAB中对比两幅图像各种直方图的程序文章来源:中国科学院计算技术研究所作者:刘伟民- 该文章讲述了MATLAB中对比两幅图像各种直方图的程序. MATLAB中对比两幅图像各种直方图的程序,程序如下。figure(1;hold on;aa = 'e:d1.jpg'bb = 'e:d4.jpg'imhist(rgb2gray(imread(aa;hold off;figure(2;hold on;imhist(rgb2gray(imread(bb;hold off;imga = imread(aa;imgb = imread(bb;figure(3;hold

2、on;title('imga1'imhist(imga(:,:,1;hold off;figure(4;hold on;title('imga2'imhist(imga(:,:,2;hold off;figure(5;hold on;title('imga3'imhist(imga(:,:,3;hold off;figure(6;hold on;title('imgb1'imhist(imgb(:,:,1;hold off;figure(7;hold on;title('imgb2'imhist(imgb(:,:,

3、2;hold off;figure(8;hold on;title('imgb3'imhist(imgb(:,:,3;hold off;3 matlab图像配准的方法-医学图像配准目前,图像配准技术在医学领域应用得相当广泛。医学图像配准具有很重要的临床应用价值。对各种使用相同或者不同的成像手段所获得的医学图像进行配准不仅可以应用于医疗诊断,也可以应用于手术计划的制定、治疗计划的制定、病理情况的跟踪和治疗效果的评价等各个方面。下面给出一个具体的图像配准的例子来形象说明配准的应用。【2】医学图像配准。(1 读取图像数据。因为源图像以矩阵形式存在一个二进制的文件里,用fread可将其

4、读取到变量矩阵中。将读取文件编制成一个子函数(RTIread.m。源代码如下:function imMatrix=RTIread(FILENAME,SIZE%RTIread Read the image matrix from binary "Registration Test Image" file.% imMatrix=RTIread(FILENAME,SIZE opens the file FILENAME, and reads the% number of elements specified by SIZE.% FILENAME is a string conta

5、ining the name of the file to be opened.% Valid entries for SIZE are:% N read N elements into a column vector.% inf read to the end of the file.% M,N read elements to fill an M-by-N matrix, in column order.% N can be inf, but M can't.% It returns the image matrix.fid=fopen(FILENAME,'r'im

6、Matrix=fread(fid,SIZE,'uint8=>uint8'fclose(fid;运行以下代码读取图像矩阵:base=RTIread('sag1.tif',256,256;input=RTIread(sag2.tif,256,256;(2 选取匹配点。根据预定的配准方法,选定足够的匹配点对。运行下列代码:cpselect(input,base; %please select 15 points for test.注意:这里允许选择15个点进行测试。操作很简单,只需注意选点要均匀布开,以增加其代表性。选定完毕,再选择FileSave Points

7、 to Workspace命令将数据保存到工作区中。Workspace立刻多出两个N×2的数组(其中N为选定的匹配点对数,分别为input_points和base_points。(3 利用十字相关法调整选定了的匹配点。如果对上面一个步骤的选点工作不满意,可以进行这一步。如果感觉尚可,则跳过这一步,以减小运算量。运行代码:input_points_corr = cpcorr(input_points,base_points,input,base;%优化选择点的作用input_points_corr为优化后在输入图片的对应匹配点。(4 计算变换公式的参数。利用cp2tform选定变换类型

8、(配准方法,计算变换参数。以下只需选定一种即可。% (1 not Fine-tune pointsTlinear = cp2tform(input_points,base_points,'linear conformal'Taffine = cp2tform(input_points,base_points,'affine'Tprojective = cp2tform(input_points,base_points,'projective'Tpolynomial2 = cp2tform(input_points,base_points,'polynomial',2;Tpolynomial3 = c