%function [ output_args ] = TestDraw( input_args )
clear all;close all;clc;
load('rabbit294X13.mat.');x1=120;x2=380;y1=110;y2=370;
Img=rabbit_Av_diastole;
I = Img;%gaussianBlur(Img,2);
I8=uint8(I);
Id=double(I);
I_norm=Id/(max(Id(:)));
imagesc(I_norm); colormap(gray);axis ([x1 x2 y1 y2]);axis('square', 'on');set(gcf,'outerposition', get(0,'screensize'));hold on;
button=1;ii=1;
while button~=3
[a,b,button]=ginput(1);
hold on;plot(a,b,'r');
if ii~=1
x=[x a];
y=[y b];
plot(x(ii-1:ii),y(ii-1:ii),'r');
else
x=a;y=b;
end
ii=ii+1;
end
plot([x(1) x(ii-1)],[y(1) y(ii-1)],'r');
[x,y] = snakeinterp(x,y,4,1);
TrueEdgeR294X13_Intra=x+i*y;
savefile = 'TrueEdgeR294X13_Intra';
save(savefile,'TrueEdgeR294X13_Intra');
clear all;close all;
load('rabbit294X13.mat.');x1=120;x2=380;y1=110;y2=370;
Img=rabbit_Av_diastole;
I = Img;%gaussianBlur(Img,2);
I8=uint8(I);
Id=double(I);
I_norm=Id/(max(Id(:)));
imagesc(I_norm); colormap(gray);axis ([x1 x2 y1 y2]);axis('square', 'on');set(gcf,'outerposition', get(0,'screensize'));hold on;
button=1;ii=1;
while button~=3
[a,b,button]=ginput(1);
hold on;plot(a,b,'r');
if ii~=1
x=[x a];
y=[y b];
plot(x(ii-1:ii),y(ii-1:ii),'r');
else
x=a;y=b;
end
ii=ii+1;
end
plot([x(1) x(ii-1)],[y(1) y(ii-1)],'r');
[x,y] = snakeinterp(x,y,4,1);
TrueEdgeR294X13_Extra=x+i*y;
savefile = 'TrueEdgeR294X13_Extra';
save(savefile,'TrueEdgeR294X13_Extra');
% clear all;close all;