clear?allclose?allclc
f1=imread（‘bone.tif‘）;
f2=imfilter（im2double（f1）-fspecial（‘laplacian‘）‘replicate‘）;
f3=im2double（f1）+f2;
w1=fspecial（‘sobel‘）;
f41=imfilter（f1-w1‘replicate‘）;
f42=imfilter（f1-w1‘‘replicate‘）;
f4=abs（f41）+abs（f42）;
f5=imfilter（f4fspecial（‘average‘[5?5]）‘replicate‘）;
f6=f3.*im2double（f5）;
f7=im2double（f1）+f6;
f8=f7.^0.5;

figure（1）imshow（f1[]）;%fig3.46a
figure（2）imshow（f2[]）;%fig3.46b
figure（3）imshow（f3[]）;%fig3.46c
figure（4）imshow（f4[]）;%fig3.46d
figure（5）imshow（f5[]）;%fig3.46e
figure（6）imshow（f6[]）;%fig3.46f
figure（7）imshow（f7[]）;%fig3.46e
figure（8）imshow（f8[]）;%fig.46g

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-03?22:27??lab_左國義\
?????文件??????403645??2003-06-28?22:59??lab_左國義\bone.tif
?????文件?????????742??2012-11-03?16:00??lab_左國義\complex_spatialfilter.asv
?????文件?????????645??2012-11-03?22:26??lab_左國義\complex_spatialfilter.m