Аппроксимация графика по фотографии

imdata=imread('C:\Users\lomt\Desktop\321.jpg');
imagesc(imdata); %загрузили картинку 
grid minor
[x,y]=ginput(60); %функция на выходе даёт координаты 60 точек, которые 
% пользователь сам выбирает мышкой
x=x-x(1); % переносим начало кривой в начало координат
y=y-y(1);
[p]=polyfit(x,y,6); % аппроксимация полиномом 6-ой степени
x_new=0:0.1:100;
f=polyval(p,x_new);
figure(1);
plot(x_new,f,'r'); %строим получившуюся кривую
grid on

imdata=imread('C:\Users\lomt\Desktop\321.jpg');
figure(1)
image(imdata); %загрузили картинку 
grid minor
[x,y]=ginput(60); %функция на выходе даёт координаты 60 точек, которые 
% пользователь сам выбирает мышкой
x=x-x(1,1); % переносим начало кривой в начало координат
y=y-y(1,1);
p=polyfit(x,y,6); % аппроксимация полиномом 6ой степени
x_new=0:0.1:150;
f=polyval(p,x_new);
figure(2);
plot(x_new,-f,'r'); %строим получившуюся кривую
grid on

clc; close; clear;
img = imread('data.jpg');
img1 = rgb2gray(img);
image(img1)