2018吉林省数学建模竞赛A题__自己做的相关结果和代码

2018吉林省数学建模竞赛A题__自己做的相关结果和代码 用到的知识:广义霍夫变换 霍夫变换 最小外接矩形

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\n‘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_script.m
     文件        4517  2016-07-22 15:12  GeneralHough.py
     目录           0  2018-04-30 14:37  gh\
     文件        1204  2011-03-27 18:54  gh\1.jpg
     文件       10728  2011-03-26 19:10  gh\2.jpg
     文件       20331  2011-03-27 19:15  gh\3.jpg
     文件        1631  2011-03-27 19:16  gh\4.jpg
     文件        1051  2011-03-27 19:26  gh\5.jpg
     文件        1850  2011-03-27 19:28  gh\6.jpg
     文件       15243  2011-03-20 17:23  gh\7.jpg
     文件        2309  2011-03-20 17:24  gh\8.jpg
     文件        1106  2018-04-30 14:45  gh\ght.m
     文件         837  2018-04-30 14:47  gh\main.m
     文件        1056  2011-03-27 21:14  gh\main1.m
     文件       39831  2011-03-27 16:27  gh\noise.jpg
     文件        6731  2011-03-27 16:29  gh\original.jpg
     文件        1059  2011-03-26 17:36  gh\rtable.m
     文件        2479  2011-03-26 17:48  gh\template1.jpg
     文件        2335  2011-03-27 15:13  gh\template2.jpg
     文件        2265  2011-03-27 15:17  gh\template3.jpg
     目录           0  2018-05-01 13:40  ght\
............此处省略229个文件信息