【正文】 [28] 覃征，[J].西安::7376[29] 張潔，蔣寧，[J].電子科技大學學報，: 5759附錄對醫(yī)學PET圖像和MRI圖像進行小波變換融合的實現(xiàn)程序如下：clc。clear all。I1=imread(39。 TC 39。)。 %讀圖像I2=imread(39。 MR 39。)。 M1=double(I1)。M2=double(I2)。zt=3。 %maximum deposition levelap=1。mp=3。i=fuse_sih(M1, M2, zt, ap, mp)。figure(4)subplot(1,3,1)imshow(i,[0,255])。title(39。 Wavelet fused image39。 ) 。colormap(gray)。axis image off。subplot(1,3,2)imshow(M1,[0,255])。title(39。Original image139。 ) 。colormap(gray)。axis image off。subplot(1,3,3)imshow(M2,[0,255])。title(39。Original image239。 ) 。colormap(gray)。axis image off。 基于小波的CT和MRI的功能圖像與解剖圖像的融合程序：clc。clear all。I1=imread(39。39。)。 %讀圖像I2=imread(39。39。)。 M1=double(I1)。M2=double(I2)。zt=3。 %maximum deposition levelap=1。mp=3。i=fuse_sih(M1, M2, zt, ap, mp)。figure(4)subplot(1,3,1)imshow(i,[0,255])。title(39。 Wavelet fused image39。 ) 。colormap(gray)。axis image off。subplot(1,3,2)imshow(M1,[0,255])。title(39。Original image139。 ) 。colormap(gray)。axis image off。subplot(1,3,3)imshow(M2,[0,255])。title(39。Original image239。 ) 。colormap(gray)。axis image off。%函數(shù)fuse_sih程序：function Y = fuse_sih(M1, M2, zt, ap, mp)%Y = fuse_sih(M1, M2, zt, ap, mp) image fusion with SIDWT, Wavelet is Haar%% M1 input image A% M2 input image B% zt maximum deposition level% ap coefficient selection highpass (see ) % mp coefficient selection base image (see ) %% Y fused image % (Oliver Rockinger )% check inputs [z1 s1] = size(M1)。[z2 s2] = size(M2)。if (z1 ~= z2) | (s1 ~= s2) error(39。Input images are not of same size39。)。end。% cells for selected imagesE = cell(3,zt)。% loop over deposition depth analysisfor i1 = 1:zt % calculate and store actual image size [z s] = size(M1)。 zl(i1) = z。 sl(i1) = s。 % define actual filters (inserting zeros between coefficients) h1 = [zeros(1,floor(2^(i12))), , zeros(1,floor(2^(i11)1)), , zeros(1,max([floor(2^(i12)),1]))]。 g1 = [zeros(1,floor(2^(i12))), , zeros(1,floor(2^(i11)1)), , zeros(1,max([floor(2^(i12)),1]))]。 fh = floor(length(h1)/2)。 % image A Z1 = conv2(es(M1, fh, 1), g1, 39。valid39。)。 A1 = conv2(es(Z1, fh, 2), g139。,39。valid39。)。 A2 = conv2(es(Z1, fh, 2), h139。,39。valid39。)。 Z1 = conv2(es(M1, fh, 1), h1, 39。valid39。)。 A3 = conv2(es(Z1, fh, 2), g139。,39。valid39。)。 A4 = conv2(es(Z1, fh, 2), h139。,39。valid39。)。 % image B Z1 = conv2(es(M2, fh, 1), g1, 39。valid39。)。 B1 = conv2(es(Z1, fh, 2), g139。,39。valid39。)。 B2 = conv2(es(Z1, fh, 2), h139。,39。valid39。)。 Z1 = conv2(es(M2, fh, 1), h1, 39。valid39。)。 B3 = conv2(es(Z1, fh, 2), g139。,39。valid39。)。 B4 = conv2(es(Z1, fh, 2), h139。,39。valid39。)。 % select coefficients and store them E(1,i1) = {selc1(A1, B1, ap)}。 E(2,i1) = {selc1(A2, B2, ap)}。 E(3,i1) = {selc1(A3, B3, ap)}。 % copy input image for next deposition stage M1 = A4。 M2 = B4。 end。% select base coefficients of last depostion stageA4 = selb(A4,B4,mp)。% loop over deposition depth synthesisfor i1 = zt:1:1 % define actual filters (inserting zeros between coefficients) h2 = fliplr([zeros(1,floor(2^(i12))), , zeros(1,floor(2^(i11)1)), , zeros(1,max([floor(2^(i12)),1]))])。 g2 = fliplr([zeros(1,floor(2^(i12))), , zeros(1,floor(2^(i11)1)), , zeros(1,max([floor(2^(i12)),1]))])。 fh = floor(length(h2)/2)。 % filter (rows) A4 = conv2(es(A4, fh, 2), h239。, 39。valid39。)。 A3 = conv2(es(E{3,i1}, fh, 2), g239。, 39。valid39。)。 A2 = conv2(es(E{2,i1}, fh, 2), h239。, 39。valid39。)。 A1 = conv2(es(E{1,i1}, fh, 2), g239。, 39。valid39。)。 % filter (columns) A4 = conv2(es(A4+A3, fh, 1), h2, 39。valid39。)。 A2 = conv2(es(A2+A1, fh, 1), g2, 39。valid39。)。 % add images A4 = A4 + A2。end。% copy imageY = A4。%函數(shù)es程序：function Y = es(X, n, wo)%Y = ES2(X, n) symmetric extension of a matrix on selected borders% % X input matrix% n number of rows/columns to extend% wo where to extend% wo == 1: left and right% wo == 2: up and bottom% wo == 3: all borders%% Y extended matrix% (Oliver Rockinger )[z s] = size(X)。if wo == 1 Y = zeros(z, s+2*n)。 Y(:,n:1:1) = X(:,2:1:n+1)。 Y(:,n+1:1:n+s) = X(:,:)。 Y(:,n+s+1:1:s+2*n) = X(:,s1:1:sn)。 end。if wo == 2 Y = zeros(z+2*n, s)。 Y(n:1:1,:) = X(2:1:n+1,:)。 Y(n+1:n+z,:) = X(:,:)。 Y(n+z+1:1:z+2*n,:) = X(z1:1:zn,:)。 end。if wo == 3 Y = zeros(z+2*n, s+2*n)。 Y(n+1:n+z,n:1:1) = X(:,2:1:n+1)。 Y(n+1:n+z,n+1:1:n+s) = X。 Y(n+1:n+z,n+s+1:1:s+2*n) = X(:,s1:1:sn)。 Y(n:1:1,n+1:s+n) = X(2:1:n+1,:)。 Y(n+1:n+z,n+1:s+n) = X。 Y(n+z+1:1:z+2*n,n+1:s+n) = X(s1:1:sn,:)。 end。25