clc;
clear;

pic?=?imread（‘../../pic/Fig0340（a）（dipxe_text）.tif‘）;
figure（1）;
imshow（pic[0255]）;
pic?=?double（pic）;
[mn]=size（pic）;

h?=?zeros（55）;
sum?=?0;
%ρ?=?3?%高斯濾波
for?i?=?1:?5
????for?j?=?1?:?5
????????h（ij）?=?floor（exp（-（（i-3）^2?+?（j-3）^2）/（2*3*3））*1024）;
????????sum?=?sum?+?h（ij）;
????end
end


pic_exp?=?zeros（m+4n+4）;
pic_exp（3:end-23:end-2）?=?pic;
pic_exp（1:）?=?pic_exp（5:）;
pic_exp（2:）?=?pic_exp（4:）;
pic_exp（:1）?=?pic