clear all
clc
close all
image=double(imread('fig4a.tif'));
gimage_5=func_bhpfp(image,5,2);
gimage_15=func_bhpfp(image,15,2);
gimage_30=func_bhpfp(image,30,2);
gimage_80=func_bhpfp(image,80,2);
figure
subplot(231);
imagesc(image);
colormap(gray);
title('Original')
subplot(232);
imagesc(gimage_5);
colormap(gray);
title('5')
subplot(233);
imagesc(gimage_15);
colormap(gray);
title('15')
subplot(234);
imagesc(gimage_30);
colormap(gray);
title('30')
subplot(235);
imagesc(gimage_80);
colormap(gray);
title('80')