利用Matlab提取圖片曲線
更新時間：2021年08月20日 15:34:23 作者：陸嵩
一張圖片，怎樣才可以提取里面曲線的數(shù)據(jù)，從而協(xié)助我們完成其他需要的數(shù)據(jù)呢？本文通過示例來進(jìn)行說明，文中大量的代碼以及圖片都可以幫助小伙伴們了解
利用 MATLAB 提取圖片曲線
給你一張圖片，如何提取里面曲線的數(shù)據(jù)，從而利用這些數(shù)據(jù)進(jìn)行圖像重繪、加工處理、測距、擬合得到函數(shù)表達(dá)式等操作呢？
行文動機(jī)
前段時間，有個朋友問了我一個問題，大概意思就是要給圖像的流線測距離，在我的印象里面，MATLAB 是似乎沒有這種直接的功能的。
那么換個角度來理解一下這個問題，如果給你一張圖像，如何提取里面點的數(shù)據(jù)？其實，有了曲線的數(shù)據(jù)，后面想干嘛就干嘛了。
一直沒空弄這個，今天偷閑，安排?。。?/p>
圖像的讀入與裁剪
以下面的圖像作為例子。
我們先導(dǎo)入圖像，進(jìn)行簡單的裁剪。為什么要裁剪呢？其實不裁剪也沒關(guān)系，因為我后面是基于像素點的顏色來提取的曲線。如果你想提取的曲線不能通過顏色區(qū)分，那么，最好通過裁剪，把你不想要的部分盡可能地剪掉。
%% 讀入圖片，展示，有必要的話可以適當(dāng)做一些裁剪
A = imread('a.jpg');%讀取到一張圖片
imshow(A);
A = imcrop(A);%使用鼠標(biāo)裁剪一波
imshow(A);
[low_num,col_num,~] = size(A);
顏色拾取
觀察圖像發(fā)現(xiàn)，我們要提取的曲線是藍(lán)色的，所以我希望通過顏色把它區(qū)分出來。那么我們就要知道這個曲線的 RGB 值。我希望通過鼠標(biāo)點選的方式獲取到顏色值。
這里我偷個懶，直接采用了 slandarer 開源的顏色提取工具。這個模塊不是我寫的，特此聲明，請尊重原創(chuàng)。
%% 顏色提取 getcolor(); color = color_list_temp(1,:);
function getcolor
global control;
global ima;
global GUI;
global x_limit;
global y_limit;
global color;
global a;
global color_number;
global color_list;
global page;
global total_page;
global color_list_temp;
color=[];
page=1;
total_page=2;
color_list=[0 0 0];
color_list(1,:)=[];
color_number=1;
rgb_type=1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GUI.fig=figure('units','pixels',...
'position',[350 100 800 500],...
'Numbertitle','off',...
'menubar','none',...
'resize','off',...
'name','getcolor',...
'color',[0.95 0.95 0.95]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
uh1=uimenu('label','設(shè)置');
uimenu(uh1,'label','RGB類型設(shè)置','callback',@RGBset)
function RGBset(~,~)
GUI.rgbfig=figure('units','pixels',...
'position',[360 370 180 200],...
'Numbertitle','off',...
'menubar','none',...
'name','RGBset',...
'resize','off');
GUI.axes=axes('Units','pixels',...
'parent',GUI.rgbfig,...
'PlotBoxAspectRatio',[1 1 1],...
'Color',[0.95 0.95 0.95],...
'Box','on', ...
'XLim',[0 500],...
'YLim',[0 500], ...
'XColor',[0.95 0.95 0.95],...
'YColor',[0.95 0.95 0.95],...
'YDir','reverse', ...
'xtick',[],'ytick',[]);
GUI.checkbox1=uicontrol('parent',GUI.rgbfig,...
'style','checkbox',...
'string','范圍:0－1',...
'position',[45 150 400 30],...
'fontsize',10,...
'value',rgb_type,...
'callback',@ifon1);
GUI.checkbox255=uicontrol('parent',GUI.rgbfig,...
'style','checkbox',...
'string','范圍:0－255',...
'position',[45 120 400 30],...
'fontsize',10,...
'value',~rgb_type,...
'callback',@ifon255);
GUI.makesurebutton=uicontrol('parent',GUI.rgbfig,...
'style','pushbutton',...
'string','確定設(shè)置',...
'position',[45 70 100 25],...
'fontsize',10,...
'callback',@settype);
function ifon1(~,~)
if(get(GUI.checkbox1,'value')==1)
set(GUI.checkbox1,'value',1);
set(GUI.checkbox255,'value',0);
else
set(GUI.checkbox1,'value',1);
end
end
function ifon255(~,~)
if(get(GUI.checkbox255,'value')==1)
set(GUI.checkbox255,'value',1);
set(GUI.checkbox1,'value',0);
else
set(GUI.checkbox255,'value',1);
end
end
function settype(~,~)
rgb_type=get(GUI.checkbox1,'value');
if ~isempty(color)
set(GUI.text2,'string',['[',num2str((color/255).*rgb_type+color.*(~rgb_type)),']']);
end
show_color(page);
close(GUI.rgbfig)
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
uh2=uimenu('label','保存');
uimenu(uh2,'label','儲存為mat','callback',@saveas_mat)
uimenu(uh2,'label','儲存為txt','callback',@saveas_txt)
uimenu(uh2,'label','儲存為excel','callback',@saveas_exl)
uimenupic=uimenu(uh2,'label','儲存為對照圖');
uimenu(uimenupic,'label','儲存全部頁碼','callback',@saveas_pic_all);
uimenu(uimenupic,'label','儲存當(dāng)前頁碼','callback',@saveas_pic_now);
function saveas_mat(~,~)
try
[filename, pathname] = uiputfile({'*.mat','mat'});
color_list_temp=(color_list/255).*rgb_type+color_list.*(~rgb_type);
save([pathname,filename],'color_list_temp');
catch
end
end
function saveas_txt(~,~)
try
[filename, pathname] = uiputfile({'*.txt','記事本'});
color_list_temp=(color_list/255).*rgb_type+color_list.*(~rgb_type);
[m,n]=size(color_list_temp);
fid=fopen([ pathname,filename],'w');
for ii=1:m
for jj=1:n
if jj==n
fprintf(fid,'%d\r\n',color_list_temp(ii,jj));
else
fprintf(fid,'%d\r\t',color_list_temp(ii,jj));
end
end
end
fclose(fid);
catch
end
end
function saveas_exl(~,~)
[filename, pathname] = uiputfile({'*.xlsx','記事本'});
color_list_temp=(color_list/255).*rgb_type+color_list.*(~rgb_type);
xlswrite([ pathname,filename],color_list_temp)
end
function saveas_pic_all(~,~)
page_with_color=total_page-1;
px=50;
gap_px=10;
pic=ones(9*px,page_with_color*px+(page_with_color-1)*gap_px,3);
for p=1:page_with_color
for ii=(p-1)*9+1:p*9
for kk=1:3
if ii<=length(color_list)
pic((ii-(p-1)*9-1)*px+1:(ii-(p-1)*9)*px,(p-1)*(px+gap_px)+1:(p-1)*(px+gap_px)+px,kk)=color_list(ii,kk)/255;
end
end
end
end
[filename, pathname] = uiputfile({'*.jpg;*.png','All Image Files';...
'*.jpg','JPG';'*.png','PNG' });
imwrite(pic,[pathname,filename]);
end
function saveas_pic_now(~,~)
try
[m,~]=size(color_list);
m=m-(page-1)*9;
m(m>9)=9;
px=50;
pic=ones(9*px,1*px,3);
if m>0
for ii=(page-1)*9+1:(page-1)*9+m
for kk=1:3
pic((ii-1)*px+1:ii*px,1:px,kk)=color_list(ii,kk)/255;
end
end
else
end
[filename, pathname] = uiputfile({'*.jpg;*.png','All Image Files';...
'*.jpg','JPG';'*.png','PNG' });
imwrite(pic,[pathname,filename]);
catch
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%uh3=uimenu('label','導(dǎo)入');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GUI.text=uicontrol('parent',GUI.fig,...
'style','text',...
'string','色彩識別',...
'horizontalalign','center',...
'position',[50 440 400 30],...
'backgroundcolor',[0.85 0.89 0.85],...
'foregroundcolor','k',...
'fontsize',15);
GUI.text1=uicontrol('parent',GUI.fig,...
'style','text',...
'string','',...
'horizontalalign','center',...
'position',[460 330 100 100],...
'backgroundcolor',[1 1 1],...
'foregroundcolor','k',...
'fontsize',10);
GUI.text2=uicontrol('parent',GUI.fig,...
'style','text',...
'string','',...
'horizontalalign','center',...
'position',[350 440 210 30],...
'backgroundcolor',[1 1 1],...
'foregroundcolor','k',...
'fontsize',10);
GUI.savecolorbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','儲存顏色',...
'position',[460 290 100 30],...
'backgroundcolor',[0.85 0.89 0.85],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@save_color);
GUI.deletedatabutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','清空數(shù)據(jù)',...
'position',[460 230 100 30],...
'backgroundcolor',[0.8 0.9 0.9],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@clear_data);
GUI.deletepicbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','刪除圖片',...
'position',[460 180 100 30],...
'backgroundcolor',[0.8 0.9 0.9],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@delete_pic);
GUI.getcapbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','屏幕截圖',...
'position',[460 130 100 30],...
'backgroundcolor',[0.8 0.9 0.9],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@get_capture);
GUI.getpicbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','讀取圖片',...
'position',[460 80 100 30],...
'backgroundcolor',[0.8 0.9 0.9],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@getImage);
GUI.getcolorbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'tag','recc',...
'string','獲取顏色',...
'position',[460 30 100 30],...
'backgroundcolor',[0.8 0.9 0.9],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@get_color);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:9
GUI.text=uicontrol('parent',GUI.fig,...
'tag',num2str(i),...
'style','text',...
'string','',...
'horizontalalign','left',...
'position',[600 440-40*(i-1) 30 30],...
'backgroundcolor',[1 1 1],...
'foregroundcolor','k',...
'fontsize',10);
end
for i=1:9
GUI.text=uicontrol('parent',GUI.fig,...
'tag',[num2str(i),'t'],...
'style','text',...
'string','',...
'horizontalalign','center',...
'position',[640 440-40*(i-1) 150 30],...
'backgroundcolor',[1 1 1],...
'foregroundcolor','k',...
'fontsize',8);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GUI.inputbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','清除最后一個顏色',...
'position',[600 80 190 30],...
'backgroundcolor',[0.85 0.89 0.85],...
'foregroundcolor','k',...
'fontsize',15,...
'callback',@delete_last);
%GUI.inputbutton=uicontrol('parent',GUI.fig,...
%'style','pushbutton',...
%'string','導(dǎo)出數(shù)據(jù)',...
%'position',[600 30 190 30],...
%'backgroundcolor',[0.85 0.89 0.85],...
%'foregroundcolor','k',...
%'fontsize',15,...
%'callback',@output_data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GUI.lpbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','<上一頁',...
'position',[600 30 70 30],...
'backgroundcolor',[0.85 0.85 0.85],...
'foregroundcolor','k',...
'fontsize',12,...
'callback',@lastpage);
GUI.npbutton=uicontrol('parent',GUI.fig,...
'style','pushbutton',...
'string','下一頁>',...
'position',[720 30 70 30],...
'backgroundcolor',[0.85 0.85 0.85],...
'foregroundcolor','k',...
'fontsize',12,...
'callback',@nextpage);
GUI.page=uicontrol('parent',GUI.fig,...
'style','text',...
'string',[num2str(page),'/',num2str(total_page)],...
'horizontalalign','center',...
'position',[670 30 50 27],...
'backgroundcolor',[0.95 0.95 0.95],...
'foregroundcolor','k',...
'fontsize',12);
function lastpage(~,~)
page=page-1;
page(page<1)=1;
set(GUI.page,'string',[num2str(page),'/',num2str(total_page)]);
show_color(page);
end
function nextpage(~,~)
page=page+1;
page(page>total_page)=total_page;
set(GUI.page,'string',[num2str(page),'/',num2str(total_page)]);
show_color(page);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GUI.axes=axes('Units','pixels',...
'PlotBoxAspectRatio',[1 1 1],...
'Position',[50 30 400 400],...
'Color',[0.98 0.98 0.98],...
'Box','on', ...
'XLim',[0 500],...
'YLim',[0 500], ...
'XColor','w','YColor','w',...
'YDir','reverse', ...
'Tag','picbagaxes',...
'xtick',[],'ytick',[]);
hold on;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function output_data(~,~)
%if color_number>1
%disp(color_list)
%end
%end
function show_color(cur_page)
len_list=size(color_list,1);
for ii=(cur_page-1)*9+1:(cur_page-1)*9+9
if(ii<=len_list)
set(findobj('tag',num2str(ii-(cur_page-1)*9)),'backgroundcolor',color_list(ii,:)/255)
set(findobj('tag',[num2str(ii-(cur_page-1)*9),'t']),'string',['[',num2str((color_list(ii,:)/255).*rgb_type+color_list(ii,:).*(~rgb_type)),']'])
else
set(findobj('tag',num2str(ii-(cur_page-1)*9)),'backgroundcolor',[1 1 1])
set(findobj('tag',[num2str(ii-(cur_page-1)*9),'t']),'string','')
end
end
end
function save_color(~,~)
if ~isempty(get(GUI.text2,'string'))
%set(findobj('tag',num2str(color_number)),'backgroundcolor',color/255)
%set(findobj('tag',[num2str(color_number),'t']),'string',['[',num2str((color/255).*rgb_type+color.*(~rgb_type)),']'])
color_list(color_number,:)=color;%(color/255).*rgb_type+color.*(~rgb_type);
color_number=color_number+1;
if color_number-1>(total_page-1)*9
page=ceil(color_number/9);
total_page=total_page+1;
set(GUI.page,'string',[num2str(page),'/',num2str(total_page)]);
end
show_color(page)
end
end
function delete_last(~,~)
if color_number>=1
%set(findobj('tag',num2str(color_number-1)),'backgroundcolor',[1 1 1])
%set(findobj('tag',[num2str(color_number-1),'t']),'string','')
color_list(end,:)=[];
color_number=color_number-1;
if color_number-2<=(total_page-2)*9
page=ceil((color_number-1)/9);
total_page=total_page-1;
set(GUI.page,'string',[num2str(page),'/',num2str(total_page)]);
end
show_color(page)
end
end
function delete_pic(~,~)
control=0;
set(findobj('Tag','picbagaxes'),...
'XLim',[0 500],...
'YLim',[0 500],...
'Position',[50 30 400 400],...
'Color',[0.98 0.98 0.98]);
delete(a);
end
function clear_data(~,~)
control=0;
set(GUI.text1,'backgroundcolor',[1 1 1]);
set(GUI.text2,'string','');
set(findobj('Tag','picbagaxes'),...
'XLim',[0 500],...
'YLim',[0 500],...
'Position',[50 30 400 400],...
'Color',[0.98 0.98 0.98]);
set(findobj('tag','recc'),'string','獲取顏色');
delete(a);
end
function get_color(~,~)
if control==0
set(GUI.text1,'backgroundcolor',[1 1 1]);
set(GUI.text2,'string','');
end
control=1;
set(gcf,'WindowButtonMotionFcn',@whilemovefcn)
set(gcf,'WindowButtonDownFcn',@whileclickfcn)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function whilemovefcn(~,~)
xy=get(gca,'CurrentPoint');
x=xy(1,2);y=xy(1,1);
if x<=x_limit&&y<=y_limit&&x>=0&&y>=0
x(x>x_limit)=x_limit;
y(y>y_limit)=y_limit;
x(x<1)=1;
y(y<1)=1;
x=round(x);
y=round(y);
if control==1
color=double([ima(x,y,1),ima(x,y,2),ima(x,y,3)]);
set(GUI.text1,'backgroundcolor',color/255);
set(GUI.text2,'string',['[',num2str((color/255).*rgb_type+color.*(~rgb_type)),']']);
end
else
if control==1
set(GUI.text1,'backgroundcolor',[1 1 1]);
set(GUI.text2,'string','');
end
end
end
function whileclickfcn(~,~)
xy=get(gca,'CurrentPoint');
x=xy(1,2);y=xy(1,1);
if x<=x_limit&&y<=y_limit&&x>=0&&y>=0
control=0;
set(GUI.text1,'backgroundcolor',color/255);
set(GUI.text2,'string',['[',num2str((color/255).*rgb_type+color.*(~rgb_type)),']']);
set(findobj('tag','recc'),'string','繼續(xù)取色');
%disp(color/255)
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function get_capture(~,~)
screensize=get(0,'screensize');
screensize=1.5*screensize;
robot=java.awt.Robot();
rectangle=java.awt.Rectangle();
rectangle.x=0;
rectangle.y=0;
rectangle.width=screensize(3);
rectangle.height=screensize(4);
image=robot.createScreenCapture(rectangle);
data=image.getData();
temp=zeros(screensize(3)*screensize(4)*3,1);
temp=data.getPixels(0,0,screensize(3),screensize(4),temp);
temp=uint8(temp);
R=temp(1:3:end);
G=temp(2:3:end);
B=temp(3:3:end);
R=reshape(R,[screensize(3),screensize(4)]);
G=reshape(G,[screensize(3),screensize(4)]);
B=reshape(B,[screensize(3),screensize(4)]);
R=R';
G=G';
B=B';
x_limit=screensize(4);
y_limit=screensize(3);
leng=max([x_limit,y_limit]);
set(findobj('Tag','picbagaxes'),...
'XLim',[0 leng],...
'YLim',[0 leng]);
ima=cat(3,R,G,B);
delete(a);
a=imshow(ima);
end
function getImage(~,~)
warning off;
try
[filename, pathname] = uigetfile({'*.jpg;*.tif;*.png;*.gif','All Image Files';...
'*.*','All Files' });
ima = imread([ pathname,filename]);
[x,y,~]=size(ima);
x_limit=x;y_limit=y;
leng=max([x_limit,y_limit]);
set(findobj('Tag','picbagaxes'),...
'XLim',[0 leng],...
'YLim',[0 leng]);
delete(a);
a=imshow(ima);
catch
end
end
end
顏色轉(zhuǎn)換與色差計算
用 RGB 比較顏色之間的相似度時，存在很大的問題，不建議直接使用，因為往往一個通道的一點改變，會導(dǎo)致最后融合在一起的顏色發(fā)生巨大變化，而如果三個通道的同時改變，卻只會使最后的明暗發(fā)生變化，色調(diào)并不會產(chǎn)生巨大變化。而這也是H系列色彩空間普遍存在的問題。
所以，經(jīng)過思考，我決定 RGB 的值轉(zhuǎn)為 HSV 顏色空間，再進(jìn)行色差計算。
色差計算，直接利用兩個點的顏色在椎體上的歐式距離，即通過計算如下坐標(biāo)點的歐式距離。我們通過編寫 color_dist 函數(shù)計算。
%計算兩個 HSV 顏色之間的距離 function d = color_dist(color_std,color) [x0,y0,z0] = getPos(color_std(1),color_std(2),color_std(3)); [x,y,z] = getPos(color(1),color(2),color(3)); d = sqrt((x-x0).^2+(y-y0).^2+(z-z0)^2); end function [x,y,z] = getPos(H,S,V) r = 1; h = sqrt(3); x = r*V*S*cos(H); y = r*V*S*sin(H); z = h*(1-V); end
%% 顏色轉(zhuǎn)為 HSV，再進(jìn)行色差計算
A2 = rgb2hsv(A);
color2 = rgb2hsv(color);
D = ones(low_num,col_num)*2;
for i=1:low_num
for j=1:col_num
D(i,j) = color_dist(color2,A2(i,j,:));
end
end
mesh(D);
分離曲線
通過調(diào)節(jié)閾值參數(shù)，可以把我們想要的坐標(biāo)軸過濾掉。閾值越小，過濾效果越明顯。
%% 根據(jù)色彩，把想要的曲線分離出來
threshold = 0.5;%可以調(diào)整閾值使分離效果變好
I = (D<threshold);
for i=1:3
RGB= A(:,:,i);
RGB(~I)=255;
A3(:,:,i) = RGB;
end
imshow(A3)
二值化，提取數(shù)據(jù)
二值化圖像是以矩陣形式存儲的。我們根據(jù)色素點在矩陣中的位置，利用行列指標(biāo)和坐標(biāo)軸的標(biāo)準(zhǔn)化關(guān)系，提取數(shù)據(jù)，重建坐標(biāo)系。
%% 二值化，提取數(shù)據(jù) leftUp = [0,1];%標(biāo)注截取的圖片的左上角和右下角 rightDown = [1,0]; B = rgb2gray(A3); imshow(B); [pos_row,pow_col] = find(B~=255); Ps = [pos_row-1,pow_col-1]; Ps(:,1) = Ps(:,1)./(low_num-1); Ps(:,2) = Ps(:,2)./(col_num-1); x = (rightDown(1)-leftUp(1)).*Ps(:,2)+leftUp(1); y = leftUp(2)-abs(rightDown(2)-leftUp(2)).*Ps(:,1); scatter(x,y,0.38); X = [x,y];
看著八九不離十了，但是注意這里的坐標(biāo)點是以散點的形式畫出來的。是無序且沒有區(qū)分度的，我們甚至無法使用plot。
數(shù)據(jù)點分類與排序
肉眼可見，這幾條曲線是分隔開的。我們?nèi)绾伟堰@些數(shù)據(jù)按曲線分開且其上的點按順序排好呢？我的做法是，使用鼠標(biāo)選中你在意的曲線的一端，利用距離延拓法，還原整條曲線。
%% 數(shù)據(jù)分類與排序 N = size(X,1); gfrom = ginput(1); [~,minI] = min(sum((X - repmat(gfrom,N,1)).^2,2)); from = X(minI,:); X(minI,:) = []; X = [from;X]; tol = 0.01; X1 = findcurvepath(X,tol); plot(X1(:,1),X1(:,2))
其中用到了一個自己寫函數(shù) findcurvepath，它可以從一個數(shù)據(jù)點出發(fā)，把所有的數(shù)據(jù)點按距離遠(yuǎn)近，像串珍珠一樣串在一塊。
function ps1 = findcurvepath(ps0,tol)
%這個函數(shù)將點就近連接起來
ps1(1,1:2) = ps0(1,1:2);%將第一個取出來
ps0 = ps0(2:end,:);%ps0重置為剩下的部分
p = ps1(end,:);%p是取出來的最后一個
while ~isempty(ps0)
i = 1;
while 1
distances = p2psdistance(p,ps0);
inds = find(distances == min(distances));%%%%%%%修復(fù)一個bug
if length(inds) > 1
if size(ps1,1)<i+1
break;
end
p = ps1(end-i,:);%p是取出來的最后一個
i = i+1;
else
break;
end
end
distances = p2psdistance(p,ps0);
if min(distances)>tol
break;
end
inds = find(distances == min(distances));%%%%%%%修復(fù)一個bug
ind = inds(1);
ps1(end+1,1:2) = ps0(ind,1:2);
ps0(ind,:) = [];
p = ps1(end,:);
end
end
function distances = p2psdistance(p,ps)
distances = sqrt((p(1)-ps(:,1)).^2+(p(2)-ps(:,2)).^2);
end
后話
我不喜歡做 GUI，因為亂七八糟的界面用起來確實令人心煩。另外一方面，不管是 GUI 或者說進(jìn)一步的 exe 桌面程序，都依賴于 MATLAB 的環(huán)境。也就是說，你要運(yùn)行代碼，就必須先裝 MATLAB 的環(huán)境。試想，一個安裝過和簡單用過 MATLAB 的人，都是有能力可以通過代碼修改一些參數(shù)的，那么你做 GUI 不是畫蛇添足么。
到此這篇關(guān)于利用Matlab提取圖片曲線的文章就介紹到這了,更多相關(guān)Matlab提取圖片曲線內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
django-rest-swagger對API接口注釋的方法
今天小編就為大家分享一篇django-rest-swagger對API接口注釋的方法，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2019-08-08
這篇文章主要介紹了基于python 取余問題(%)詳解，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2020-06-06
python實現(xiàn)將中文日期轉(zhuǎn)換為數(shù)字日期
這篇文章主要介紹了python實現(xiàn)將中文日期轉(zhuǎn)換為數(shù)字日期，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2020-07-07
大家好，本篇文章主要講的是一款開源的Python一鍵搶票神器，感興趣的同學(xué)趕快來看一看吧，對你有幫助的話記得收藏一下2022-02-02
Python遠(yuǎn)程桌面協(xié)議RDPY安裝使用介紹
這篇文章主要介紹了Python遠(yuǎn)程桌面協(xié)議RDPY安裝使用介紹,本文講解了RDPY的安裝、RDPY的簡單使用兩部份內(nèi)容,需要的朋友可以參考下2015-04-04
詳解python opencv、scikit-image和PIL圖像處理庫比較
這篇文章主要介紹了詳解python opencv、scikit-image和PIL圖像處理庫比較，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2019-12-12
這篇文章主要介紹了基于CUDA out of memory的一種神奇解決方式，具有很好的參考價值，希望對大家有所幫助。如有錯誤或未考慮完全的地方，望不吝賜教2023-02-02
這篇文章主要介紹了4款Python 類型檢查工具的相關(guān)資料，幫助是及早檢查，提前發(fā)現(xiàn)類型的錯誤，增強(qiáng)代碼的一致性與可維護(hù)性。（還有防止脫發(fā)，喵），感興趣的朋友可以了解下2020-10-10
Python使用asyncio實現(xiàn)異步操作的示例
本文主要介紹了Python使用asyncio實現(xiàn)異步操作的示例,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2025-01-01
Python?Pandas學(xué)習(xí)之series的二元運(yùn)算詳解
二元運(yùn)算是指由兩個元素形成第三個元素的一種規(guī)則，例如數(shù)的加法及乘法;更一般地,由兩個集合形成第三個集合的產(chǎn)生方法或構(gòu)成規(guī)則稱為二次運(yùn)算。本文將詳細(xì)講講Pandas中series的二元運(yùn)算，感興趣的可以了解一下2022-09-09