close all; clear all; clc; imgray= rgb2gray(imread('1.png')); %% edge bw=edge(imgray,'sobel'); % figure,imshow(bw); I=double(imgray); [x y]=gradi……