| 查看: 1125 | 回复: 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人回复
-
闲聊
已经有6人回复
-
江苏省自然基金 什么时候出结果
已经有7人回复
-
浙江师范大学是怎么坑我的
已经有10人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有8人回复
-
科研求助
已经有4人回复