1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
|
f = figure('name','Selection,'numbertitle','off');
[filename, path] = uigetfile( ...
{'*.mat','MAT-files (*.mat)'; ...
'Pick a file');
fullpathToFile = fullfile(path,filename);
if filename~=0
S=load(fullpathToFile);
current_data=S.db;
xaxis=S.x;
yaxis=S.distance;
subplot(3, 1, 1);
imagesc(xaxis,yaxis,current_data);
% figure
% imshow(current_data)
end
[x,y] = ginput(2)
hln=line('XData',[x(1) x(2) x(2) x(1) x(1)],'YData',[y(1) y(1) y(2) y(2) y(1)],'marker','*','markeredgecolor','b','MarkerSize',10);
xmin=round(min(x));
ymin=round(min(y));
xmax=round(max(x));
ymax=round(max(y));
distancex=(xmin:xmax);
distancey=(ymin:ymax);
new_data=current_data(ymin:ymax,xmin:xmax);
figure
RI=imref2d(size(new_data));
RI.XWorldLimits = [xmin xmax];
RI.YWorldLimits = [ymin ymax];
img2=imagesc(RI.XWorldLimits,RI.YWorldLimits,new_data); |
Partager