2018吉林省数学建模竞赛A题__自己做的相关结果和代码 用到的知识:广义霍夫变换 霍夫变换 最小外接矩形
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 198919 2018-04-28 14:52 A题--吉林省第九届大学生数学建模竞赛.pdf
文件 134895 2018-04-29 21:52 A题信息.docx
文件 231102 2018-04-28 14:57 B题--吉林省第九届大学生数学建模竞赛.pdf
文件 14538 2018-04-29 09:45 B题信息.docx
文件 1126134 2018-05-01 13:47 DATA1.bmp
文件 2370 2018-04-19 16:21 DATA1.mat
文件 1126134 2018-05-01 13:49 DATA2.bmp
文件 2580 2018-04-19 16:21 DATA2.mat
文件 3925 2018-04-19 16:21 DATA3.mat
文件 1879 2016-04-26 15:15 findcircle.m
文件 383 2018-04-29 14:04 findcircle_script.m
文件 4517 2016-07-22 15:12 GeneralHough.py
目录 0 2018-04-30 14:37 gh
文件 1204 2011-03-27 18:54 gh1.jpg
文件 10728 2011-03-26 19:10 gh2.jpg
文件 20331 2011-03-27 19:15 gh3.jpg
文件 1631 2011-03-27 19:16 gh4.jpg
文件 1051 2011-03-27 19:26 gh5.jpg
文件 1850 2011-03-27 19:28 gh6.jpg
文件 15243 2011-03-20 17:23 gh7.jpg
文件 2309 2011-03-20 17:24 gh8.jpg
文件 1106 2018-04-30 14:45 ghght.m
文件 837 2018-04-30 14:47 ghmain.m
文件 1056 2011-03-27 21:14 ghmain1.m
文件 39831 2011-03-27 16:27 gh
oise.jpg
文件 6731 2011-03-27 16:29 ghoriginal.jpg
文件 1059 2011-03-26 17:36 gh
table.m
文件 2479 2011-03-26 17:48 gh emplate1.jpg
文件 2335 2011-03-27 15:13 gh emplate2.jpg
文件 2265 2011-03-27 15:17 gh emplate3.jpg
目录 0 2018-05-01 13:40 ght
............此处省略229个文件信息
function [para] = findcircle(BWstepRstepAngleminRmaxRp)
[mn] = size(BW);% BW:二值图像;
cntR = round((maxR-minR)/stepR)+1;% stepR:检测的圆半径步长, minR:最小圆半径, maxR:最大圆半径
cntAngle = round(2*pi/stepAngle);% stepAngle:角度步长,单位为弧度
hough_space = zeros(mncntR);% hough_space:参数空间,h(abr)表示圆心在(ab)半径为r的圆上的点数
[rowscols] = find(BW);
cntPoints = size(rows1);
% Hough变换将图像空间(xy)对应到参数空间(abr)
% a = x-r*cos(angle) b = y-r*sin(angle)
for i=1:cntPoints
for r=1:cntR
for k=1:cntAngle
a = round(rows(i)-(minR+(r-1)*stepR)*cos(k*stepAngle));
b = round(cols(i)-(minR+(r-1)*stepR)*sin(k*stepAngle));
if(a>0 && a<=m && b>0 && b<=n)
hough_space(abr) = hough_space(abr)+1;
end
end
end
end
% 寻找满足阈值的圆的参数
max_para = max(max(max(hough_space)));
index = find(hough_space>=max_para*p); % p:以p*hough_space的最大值为阈值,p取0,1之间的数
length = size(index1);
hough_circle=zeros(mn);
for i=1:cntPoints
for k=1:length
par3 = floor(index(k)/(m*n))+1;
par2 = floor((index(k)-(par3-1)*(m*n))/m)+1;
par1 = index(k)-(par3-1)*(m*n)-(par2-1)*m;
if((rows(i)-par1)^2+(cols(i)-par2)^2<(minR+(par3-1)*stepR)^2+5 && (rows(i)-par1)^2+(cols(i)-par2)^2>(minR+(par3-1)*stepR)^2-5)
hough_circle(rows(i)cols(i)) = 1;% hough_circl:二值图像,检测到的圆
end
end
end
for k=1:length
par3 = floor(index(k)/(m*n))+1;
par2 = floor((index(k)-(par3-1)*(m*n))/m)+1; % 圆心y坐标
par1 = index(k)-(par3-1)*(m*n)-(par2-1)*m; % 圆心x坐标
par3 = minR+(par3-1)*stepR; % 圆的半径
% fprintf(1‘Center %d %d radius %d
‘par1par2par3);
para(k:) = [par1par2par3]; % para:检测到的圆的圆心、半径
end
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 198919 2018-04-28 14:52 A题--吉林省第九届大学生数学建模竞赛.pdf
文件 134895 2018-04-29 21:52 A题信息.docx
文件 231102 2018-04-28 14:57 B题--吉林省第九届大学生数学建模竞赛.pdf
文件 14538 2018-04-29 09:45 B题信息.docx
文件 1126134 2018-05-01 13:47 DATA1.bmp
文件 2370 2018-04-19 16:21 DATA1.mat
文件 1126134 2018-05-01 13:49 DATA2.bmp
文件 2580 2018-04-19 16:21 DATA2.mat
文件 3925 2018-04-19 16:21 DATA3.mat
文件 1879 2016-04-26 15:15 findcircle.m
文件 383 2018-04-29 14:04 findcircle_sc
文件 4517 2016-07-22 15:12 GeneralHough.py
目录 0 2018-04-30 14:37 gh
文件 1204 2011-03-27 18:54 gh1.jpg
文件 10728 2011-03-26 19:10 gh2.jpg
文件 20331 2011-03-27 19:15 gh3.jpg
文件 1631 2011-03-27 19:16 gh4.jpg
文件 1051 2011-03-27 19:26 gh5.jpg
文件 1850 2011-03-27 19:28 gh6.jpg
文件 15243 2011-03-20 17:23 gh7.jpg
文件 2309 2011-03-20 17:24 gh8.jpg
文件 1106 2018-04-30 14:45 ghght.m
文件 837 2018-04-30 14:47 ghmain.m
文件 1056 2011-03-27 21:14 ghmain1.m
文件 39831 2011-03-27 16:27 gh
oise.jpg
文件 6731 2011-03-27 16:29 ghoriginal.jpg
文件 1059 2011-03-26 17:36 gh
table.m
文件 2479 2011-03-26 17:48 gh emplate1.jpg
文件 2335 2011-03-27 15:13 gh emplate2.jpg
文件 2265 2011-03-27 15:17 gh emplate3.jpg
目录 0 2018-05-01 13:40 ght
............此处省略229个文件信息
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)