clc;
clear;
RGB=imread（‘E:\新建文件夾\1.jpg‘‘jpg‘）;??????????????????%讀入彩色圖片
%figure（1）imshow（RGB）title（‘彩色圖‘）;?????????%顯示彩色圖片
I=rgb2gray（RGB）;??????????????????????????%彩色轉化成灰度圖
%figure（2）imshow（I）title（‘灰度圖‘）;???%顯示灰度圖片
%邊界曲線提取
ihd=im2double（I）;
[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘ihd）;
ixc=wdencmp（‘gbl‘ihd‘sym4‘2thrsorhkeepapp）;
%figureimshow（ixc）title（‘消噪后圖像?‘）;
k2=medfilt2（ixc[3?3]）;
%figureimshow（k2）title（‘中值濾波‘）;
isuo=imresize（k21‘bicubic‘）;
%sobert、robert、prewitt、Laplacian、canny算子檢測圖像邊緣
esobel=edge（isuo‘sobel‘）;
erob=edge（isuo‘roberts‘）;
eprew=edge（isuo‘prewitt‘）;
elog=edge（isuo‘log‘）;
ecanny=edge（isuo‘canny‘）;
L=flipud（ecanny）;