| 查看: 1056 | 回复: 0 | |||
zhongliuhua1新虫 (初入文坛)
|
[交流]
最小二乘法拟合直线请教
|
|
W5是用canny算子提取后的边界，用最小二乘法对提取得边界进行直线拟合，然后再求两直线夹角，但是程序运行错误。 最小二乘法代码如下： j=0; [xa,ya]=size(BW5); for i=1:xa q=find(BW5(i,  ==1);if~isempty(q);%判断能否找到边界点 j=j+1 [q1,q2]=size(q);%q1，q2分别是q矩阵横坐标和纵坐标的个数 e0(j, =[i,q(1),q(q2)];%i是边界点的横坐标；q（1）是第一条边界曲线纵坐标；q（q2）是第2条边界曲线纵坐标；e0是由i，q1，q（q2）组成的矩阵end end p1=polyfit(e0(:,1),e0(:,2),1);%p1是左边界曲线的拟合直线系数 p2=polyfit(e0(:,1),e0(:,3),1);%p2是右边界曲线的拟合直线系数 x=1:xa;%横向扫描点 y1=polyval(p1,x);%左边界拟合所得直线的值范围 y2=polyval(p2,x);%右边界拟合所得直线的值范围 plot(x,y1,x,y2);%画出拟合直线 a=p1(1);b=p2(1); theta=atan(abs((a-b)/(1+a*b))); alfa=theta*(180/pi); |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复