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圖像配準的方法-醫(yī)學圖像配準目前,圖像配準技術在醫(yī)學領域應用得相當廣泛。醫(yī)學圖像配準具有很重要的臨床應用價值。對各種使用相同或者不同的成像手段所獲得的醫(yī)學圖像進行配準不僅可以應用于醫(yī)療診斷,也可以應用于手術計劃的制定、治療計劃的制定、病理情況的跟蹤和治療效果的評價等各個方面。下面給出一個具體的圖像配準的例子來形象說明配準的應用。【2】醫(yī)學圖像配準。(1 讀取圖像數(shù)據(jù)。因為源圖像以矩陣形式存在一個二進制的文件里,用fread可將其

4、讀取到變量矩陣中。將讀取文件編制成一個子函數(shù)(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 選取匹配點。根據(jù)預定的配準方法,選定足夠的匹配點對。運行下列代碼:cpselect(input,base; %please select 15 points for test.注意:這里允許選擇15個點進行測試。操作很簡單,只需注意選點要均勻布開,以增加其代表性。選定完畢,再選擇FileSave Points

7、 to Workspace命令將數(shù)據(jù)保存到工作區(qū)中。Workspace立刻多出兩個N×2的數(shù)組(其中N為選定的匹配點對數(shù),分別為input_points和base_points。(3 利用十字相關法調(diào)整選定了的匹配點。如果對上面一個步驟的選點工作不滿意,可以進行這一步。如果感覺尚可,則跳過這一步,以減小運算量。運行代碼:input_points_corr = cpcorr(input_points,base_points,input,base;%優(yōu)化選擇點的作用input_points_corr為優(yōu)化后在輸入圖片的對應匹配點。(4 計算變換公式的參數(shù)。利用cp2tform選定變換類型

8、(配準方法,計算變換參數(shù)。以下只需選定一種即可。% (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