自适应全变分图像去噪Matlab源代码

% Decompose textures in image with spectral TV, gray-scale images % Script by Guy Gilboa (Jan 2015). % Based on: [1] G. Gilboa, "A total variation spectral framework for scale and texture analysis." SIAM Journal on Imaging Sciences 7.4 (2014): 1937-1961. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% Important params to change according to image / application %% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Max_time = 1; % Maximal scale to be processed (in evolution time) %Max_time = 40; % large scale Num_of_bands = 50; % Number of bands phi(x,t); %% dt = Max_time/Num_of_bands; f = double(imread('fruits.bmp')); %f = f(50:100,50:100); % debug f = double(f); f = f(21:185,71:240); % fruits, melon f = f/255; % pixels are in the range [0,1] figure(1); imshow(f); title('f') % Compute Phi bands and residual f_r %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [S,T,Phi,f_r] = specTV_evolve(f, Max_time, dt); % evolve image % Define Filters hp_i = 10; bp_low_i=13; bp_high_i = 30; % incices of low and high bandpass cutoff H1 = zeros(size(T)); H1(1:lp_i)=1; % high pass, melon H2 = zeros(size(T)); H2(bp_low_i:bp_high_i)=1; % band pass, melon % Reconstruct filtered image given Phi, H and f_r %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f_H1 = specTV_filter( Phi, H1, f_r, dt ); f_H2 = specTV_filter( Phi, H2, f_r, dt ); % Plot S(t) and filter response figure(3); plot(T,S,'b',T(1:hp_i),S(1:hp_i),'g',T(bp_low_i:bp_high_i),S(bp_low_i:bp_high_i),'r'); xlabel('t'); legend('Spectrum S(t)','High Pass Filter H1(t)', 'Band Pass Filter H2(t)') % show a few Phi(t) instances %phi_show = [4 20 30 50]; figure(4); imshow(0.5+f_H1*2); title('High Pass image f_{H1} x 2') figure(5); imshow(0.5+f_H2); title('Band Pass image f_{H2}') figure(6); imshow(f_r); title('Residual f_r (Low-pass)') Ind = 1:10; % Index of phi bands RbyC = [2 5]; % Row by Column boxes fignum =10; % Num of figure Rescale = 1; % Possible image size rescaling (1 - no rescale, >1 larger images) contrast = 5; I = specTV_show_phi( Phi*contrast,dt, Ind, RbyC, fignum, Rescale); title(['A few instances of \phi(x;t), from bands ' num2str(Ind(1)) ' to ' num2str(Ind(end))])