CN116804865B - 一种三轴自动化程编特征识别与刀轨生成方法 - Google Patents
一种三轴自动化程编特征识别与刀轨生成方法 Download PDFInfo
- Publication number
- CN116804865B CN116804865B CN202311083327.5A CN202311083327A CN116804865B CN 116804865 B CN116804865 B CN 116804865B CN 202311083327 A CN202311083327 A CN 202311083327A CN 116804865 B CN116804865 B CN 116804865B
- Authority
- CN
- China
- Prior art keywords
- image
- contour
- web
- tool path
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000003754 machining Methods 0.000 claims abstract description 15
- 230000006870 function Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000013135 deep learning Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 8
- 241001270131 Agaricus moelleri Species 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 abstract description 9
- 230000000694 effects Effects 0.000 description 3
- 239000003086 colorant Substances 0.000 description 2
- 238000004040 coloring Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 101150026109 INSR gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35349—Display part, programmed locus and tool path, traject, dynamic locus
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种三轴自动化程编特征识别与刀轨生成方法,属于数控加工领域,包括以下核心步骤:三维建模软件载入零件模型图并设置模型的显示参数,获取零件二维平面图像,制作训练样本、学习训练、得到模型权重参数,基于一阶梯度计算识别结果轮廓,对每一特征图进行三维轮廓扫描得到邻接关系图,特征的精细化识别并结合三维建模软件得到刀轨程序。本发明可自动识别图像中对应的驱动面与支撑面,并基于识别结果生成刀轨,识别结果的准确性较现有的拓扑规则判断方式有着明显提升,模型的适应性与鲁棒性更强,可显著的降低NC程序编制时间提升效率,同时可大幅降低人工成本,为进一步在线提供零件刀轨程序生成服务奠定了基础。
Description
技术领域
本发明涉及零件的数控加工与特征编程技术,具体涉及一种三轴自动化程编特征识别与刀轨生成方法。
背景技术
在数控加工领域,要完成对所设计零部件的加工,数控加工刀轨的生成是一切工作的核心。然而现今在世界范围内还没有真正实现全自动的刀轨生成方法,基于毛坯要得到理想中的零件,需要手动的选择一组加工面即驱动面(一般为腹板)与引导面(一般为侧壁与圆角),然后基于选择的特征生成对应区域的刀轨。因此想要加工完整个零件,在设计的三维图中需将所有特征元素进行手动选择,基于此得到加工刀轨的信息。
对于熟练的数控加工程序设计工艺师来说,手动选择加工特征的主要优势有,可根据加工面的特征选择适宜的刀轨生成策略。然而手工选择加工面的缺陷也较为明显,即对于多特征(面、槽腔等)元素零件需要手动的点选所有特征面效率较低,不能满足自动化、高效率以及低成本的需要。
为了实现零件自动化程编,部分学者完成了前期的探索工作,现有的成果主要基于所设计零件所有面之间的拓扑关系与零件的几何参数生成刀轨。现有的方法主要步骤包含:导入待加工零件三维图、获取所有面与边的空间拓扑信息(拓扑图的结点表示一个面,面与面之间的连接用相邻的线表示)、选择法矢方向确定在三维空间中最底面、基于选择的面为起始逐步生成刀轨信息直至基于设定的规则得到所有面的加工刀轨信息。分析现有的基于最底面为起始刀轨生成方法主要的缺点有只能适用于零件结构较为简单的特征,若零件中最底面不存在或是零件最底面存在层层嵌套则不适用。
为实现自动化程编高效与准确的需要,同时减少在零件程编中人为的介入,准确的特征(驱动面+引导面)识别是实现自动化程编的关键,也是行业急需解决的问题。
发明内容
本发明的目的在于提供一种三轴自动化程编特征识别与刀轨生成方法,以解决现今人工生成零件加工程序效率低、高成本以及基于规则的拓扑关系判断驱动面与引导面导致的高误判问题,实现从零件输入自动判断特征面并得到数控机床加工NC程序,同时为进一步在线程序生成服务奠定了基础。
为了实现上述发明目的,本发明的技术方案如下:
一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,包括如下步骤:
S1、在三维建模软件中载入零件模型图,设置其显示参数以区分前景和背景;
S2、采集不同视角的零件表面二维投影图像;
S3、基于深度学习的分割网络粗识别图像特征区域;
S4、针对粗识别结果得到的图像进行灰度化处理;
S5、采用最大间方差策略实现灰度化后图像的二值化,得到包含前景和背景的二值图像;
S6、通过二值图像的全局一阶梯度计算获取图像特征区域的轮廓图;
S7、基于随机灰度赋值法得到图像特征区域轮廓图中的独立轮廓区域;
S8、基于激光扫描得到每个独立轮廓区域的轮廓线扫描数据,并生成二维邻接关系图;
S9、建立邻接准则判断函数实现快速程编驱动面与引导面特征组的精细识别;
S10、基于精细识别结果自动生成加工刀轨程序。
进一步的,零件表面二维投影图像为自动采集:在零件模型图中选择一特征面,以其对应的法矢方向为整个零件的视角方向,并截取屏幕图像;循环得到不同视角对应的特征面,使得零件在屏幕方向上的投影面积最大。
进一步的,零件表面二维投影图像为手动采集:利用三维建模软件菜单栏的中的工具,依次选择“图像”、“捕获”,在“选择模式”中框选出整个零件区域,点击“选项”,在弹出的“捕获选型”中选择像素,在“更多”中设置图像的分辨率宽和长;设置好捕获参数后,点击“捕获”按钮截取图像,并将截取的图像进行保存。
进一步的,所述的针对粗识别结果得到的图像进行灰度化处理包括:针对深度学习识别结果得到的图像,根据各个通道灰度值以加权求和的方式灰度化,得到灰度化后的图像。
进一步的,所述的采用最大间方差策略实现灰度化后图像的二值化包括:
对于大小为M*N包含sum个灰度级的图像,对应的计算式子为:
; (3)
; (4)
; (5)
灰度图中小于阈值像素点的平均灰度值为:
;(6)
;(7)
; (8)
其中,p i 表示灰度级为i的像素数占总像素的个数;n i 表示灰度级为i的像素数;p background 表示像素点被分配为background的概率累加和;background与foreground分别表示前景与背景;threshold表示前景与背景分割的阈值;表示图像中所有像素点的灰度值方差。
进一步的,所述的基于随机灰度赋值法得到图像特征区域轮廓图中的独立轮廓区域包括:
针对公式计算得到的图像特征区域轮廓图,采用全局像素种子点填充策略得到每一连通区域的语义信息;
每一填充区域采用随机赋值的方式进行处理,得到具有不同灰度的图像特征区域;
提取具有相同灰度信息的像素点并进行分类,从而得到独立轮廓区域;
针对每一个具有随机灰度的独立轮廓区域,分别进行灰度化和二值化处理,再基于Canny算子计算得到独立轮廓区域的轮廓线。
进一步的,所述的基于激光扫描得到每个独立轮廓区域的轮廓线扫描数据包括:
选择轮廓线任意位置为起始位置,腹板面法矢方向为激光点的移动方向,竖直扫描原始三维图,激光点发出的光线垂直于腹板轮廓线指向外侧,发出的光线遇到平面即反射对应的距离为扫描距离L;竖直方向的移动距离为零件三维图像中垂直于腹板方向的最大值H max ,当轮廓线上一个点扫描完成后沿着轮廓线顺时针方向移动激光点,直到将一个独立腹板的轮廓线全部扫描完成为止。
进一步的,所述的由轮廓线扫描数据生成的二维邻接关系图中,X轴表示不同的腹板区域轮廓线段;Y轴表示腹板法矢方向移动的距离;针对每个腹板区域独立的轮廓区域后建立激光扫描点,扫描轮廓线对应的三维图。
进一步的,所述的建立邻接准则判断函数实现快速程编驱动面与引导面特征组的精细识别,对应的判断函数为:
;(13)
; (14)
其中,sum表示腹板封闭轮廓线中存在的直线段总数量;i表示直线段的标号;flag li 表示对应线段的状态标识;、/>、/>、/>分别表示腹板轮廓线l在第i段与y di 在二维坐标系中对应区域S中图像的R、G、B灰度值;di表示底角区域;up表示沿腹板法矢方向高于底角的区域。
综上所述,本发明具有以下优点:
1、本发明解决了现今零件刀轨程序需基于人工选择驱动面加支撑面方式生成刀轨,可自动的识别图像中对应的驱动面与其对应的支撑面并基于识别结果生成刀轨;
2、由于采用了深度学习加图像处理方法,识别结果的准确性较现有的拓扑规则判断方式有着明显的提升,同时对不同模型的适应性与鲁棒性更强。采用提出的方式生成刀轨可显著的降低NC程序编制时间并提升效率,同时可大幅降低人工成本,为进一步在线提供零件刀轨程序生成服务奠定了基础;
3、本发明可自动的识别图像中对应的驱动面与其对应的支撑面,并基于识别结果生成刀轨,识别结果的准确性较现有的拓扑规则判断方式有着明显提升,模型的适应性与鲁棒性更强,可显著的降低NC程序编制时间提升效率,同时可大幅降低人工成本,为进一步在线提供零件刀轨程序生成服务奠定了基础。
附图说明
图1为本发明所述方法包含的核心步骤;
图2为自动截屏获得的零件表面图像;
图3为腹板区域训练样本图像;
图4为基于深度学习网络识别的腹板区域结果图像;
图5为基于式子计算得到的识别结果轮廓图;
图6为识别结果不同腹板区域上色处理效果;
图7为所有腹板信息的图像中提取出每一个单独的腹板;
图8为激光点扫描方式与轮廓线之间的关系;
图9为基于轮廓线扫描数据生成的二维邻接关系图。
具体实施方式
下面结合实施例对本发明作进一步地详细说明,但本发明的实施方式不限于此。
实施例1
本实施例提供了一种三轴自动化程编特征识别与刀轨生成方法,可自动识别图像中对应的驱动面与其对应的支撑面,并基于识别结果生成刀轨,识别结果的准确性较现有的拓扑规则判断方式有着明显提升。
本方法包含的核心步骤如图1所示:
步骤一、使用三维建模软件打开设计好的零件图,在软件中设置显示参数,以得到利于分析的图像。
在三维建模软件中设置背景灰度对应的RGB值,所有面的显示灰度对应的RGB数值,所有线条对应的RGB数值,设置线型和对应的线宽。具体设置方式为:设背景灰度对应的RGB值为x1、x2、x3,所有面的显示灰度对应的RGB数值为y1、y2、y3,所有线条对应的RGB数值为z1、z2、z3,线型为line,对应的线宽为w(mm),其中x1、x2、x3∈[0,255], y1、y2、y3∈[0,255]且至多两个数与x1、x2、x3相等,z1、z2、z3∈[0,255] 且至多两个数与x1、x2、x3相等,z1、z2、z3 至多两个数与y1、y2、y3相等。例如,设置背景灰度对应的RGB值为255、255、255,所有面的显示灰度对应的RGB数值为178、178、178,所有线条对应的RGB数值为0、0、0,线型为1,对应的线宽为1:0.13mm。
步骤二、获取零件表面二维投影图像
本步骤包含手动模式与自动模式。手动模式采集方式为:利用三维建模软件中的工具,依次选择“图像”——“捕获”,在“选择模式”中框选出整个零件区域,点击“选项”,在弹出的“捕获选型”中选择像素,在“更多”中设置图像的分辨率宽和长。设置好捕获参数后,点击“捕获”按钮截取图像,并将截取的图像进行保存,对应的格式为.bmp。
自动模式采集方式为:在零件图中选择一特征面,以其对应的法矢方向为整个零件的视角方向,并截取屏幕图像。循环得到不同视角对应的特征面,使得零件在屏幕方向上的投影面积最大,采用三维建模软件VB宏代码实现自动截屏,获得的零件表面图像如图2所示。
步骤三、基于深度学习的分割网络对图像特征区域(即腹板区域)进行识别。
本步骤中,训练样本的制作方式为:采用灰度填充法对腹板区域进行标注,如对应的灰度值RGB为α1、α2、α3;然后基于OpenCV库遍历图像中所有像素点的灰度值,若对应的R、G、B灰度值为y1、y2、y3,则其灰度值设置为y1、0、0,这样设置是为了对训练腹板特征进行训练,该灰度值可不唯一。而不满足此灰度值关系(即y1、y2、y3)的像素点灰度值统一设置为0、0、0,并按原始输入图像的尺寸大小对图像进行保存,使处理后的图像大小保持不变,对应的图像格式为.png。
本步骤中,图像特征区域的识别网络模型采用VGG-16提取特征,使用交叉熵损失函数对识别结果进行评价,类别包含了前景(即腹板)与背景(即腹板以外的其他区域)两大类,训练时设置的模型超参数epoch设为num(>200)。采用图像增强(旋转、缩放)的策略获取足够数量的训练样本。腹板区域训练样本图像如图3所示。
训练完成后得到记录有epoch、准确率、损失量等参数数值的.log日志文件,基于日志文件中loss最小对应的epoch保存的.pth权重文件为测试时初始化网络的参数。识别结果中腹板区域标注的灰度RGB为h1、h2、h3,除腹板之外的位置对应的RGB灰度值为k1、k2、k3(对应的灰度与h1、h2、h3皆不相同)。将识别结果与输入原图进行映射以保证识别结果图像的尺度保持一致。基于深度学习网络识别的腹板区域结果图像如图4所示,该结果为粗判断结果。
步骤四、针对识别结果得到的图像,根据识别结果图像各个通道灰度值以加权求和的方式灰度化,得到灰度化后的图像。
本步骤基于Visual Studio 与OpenCV实现,灰度化式子为:
h(x,y)=H R (x,y)*a 1+H G (x,y)*a 2+H B (x,y)*a 3 (1)
a 1+a 2+a 3=1.2(2)
其中,H R (x,y)、H G (x,y)、H B (x,y)分别表示与原始输入图像等尺度的预测图像中(x,y)位置像素点对应的红色、绿色、蓝色通道像素数值;a 1、a 2、a 3分别表示灰度化处理时红色、绿色、蓝色通道数值所占的权重;h(x,y)灰度化后对应的单通道像素灰度值。
步骤五、采用最大间方差策略实现灰度化后图像的二值化。
对于大小为M*N包含sum个灰度级的图像,对应的计算式子为:
; (3)
; (4)
; (5)
灰度图中小于阈值像素点的平均灰度值为:
;(6)
;(7)
; (8)
其中,p i 表示灰度级为i的像素数占总像素的个数;n i 表示灰度级为i的像素数;p background 表示像素点被分配为background的概率累加和;background与foreground分别表示前景与背景;threshold表示前景与背景分割的阈值;表示图像中所有像素点的灰度值方差。
步骤五、基于上述步骤可获得包含前景(即腹板)和背景(即腹板之外区域)的二值图像。其中前景为1,背景为0,然后针对所得的二值图像计算其全局一阶梯度,式子为:
(9)
梯度大小计算式子:
(10)
梯度方向计算式子:
(11)
(12)
其中,grad表示梯度;表示偏导数;/>与/>分别表示在X轴与Y轴方向上的偏导数。
得到二值图像中所有点处的导数后,设定阈值grad threshold ,将梯度值小于该阈值的像素点梯度值置为零,根据八邻域最近邻原则,将所有满足关系的像素点进行连接,得到若干独立的封闭区域,封闭区域即为识别的腹板区域。
基于公式计算得到的识别结果轮廓图如图5所示。图5为若干独立的封闭轮廓线组成图像,采用全局像素种子点填充的策略得到每一连通区域的语义信息(即灰度信息)。
具体实施时,每一个填充的区域采用随机赋值的方式,结果轮廓图以左上角为开始点,对起点像素进行随机赋RGB数值,值不超过255,且为正整数,并对与起点连通的所有点赋以相同的数值;从左至右,从上至下依次对每一个点进行扫描并随机赋RGB数值,若对应的值为黑色且不为白色则直接跳过进入下一个点,反之则采用随机的方式进行赋值,且对应的数值与之前出现过的赋值不一致;直到处理完所有的点为止。识别结果中不同腹板区域上色处理效果如图6所示。
步骤六、基于上述步骤得到具有不同颜色的腹板区域,采用循环遍历的方式得到具有相同灰度信息的像素点,并基于像素点的灰度值对其进行分类。
本步骤基于OpenCV获取图像像素,使用C++中的容器进行保存,实现的方式为,char ch;start = clock();int num = 1;for (int i = 0; i<img.size[0]; i++){for(int j = 0; j<img.size[1]; j++){if (img.at<Vec3b>(i, j)[0] == 255&&img.at<Vec3b>(i, j)[1] == 255&&img.at<Vec3b>(i, j)[2] ==255)continue;if (img.at<Vec3b>(i, j)[0] == 0&&img.at<Vec3b>(i, j)[1] == 0&&img.at<Vec3b>(i, j)[2] ==0)continue;Vec3b vc3 = img.at<Vec3b>(i, j);int vec_int = vc3[0] * 256 +vc3[1]* 257 + vc3[2]* 259; //flagint length, first = 0;for (int k = 0; k<all_location_info.size(); k++){//if (all_location_info[k].vec_int == vec_int) if(all_location_info[k].cha[vec_int] == 1){all_location_info[k].length += 1;Location loc;loc.row = i;loc.col =j;all_location_info[k].locationArray.push_back(loc);first = 1;break;}}if (first== 0){VecHead head;head.length = 1; //首个录入长度为1,head.vec = vc3;//图像位置信息head.cha[vec_int] = 1;//int te =vc3[1];std::cout<<"输出第"<<num<<"个不同灰度"<<vc3<<endl;num++;head.vec_int =vec_int;Location loc;loc.row = i;loc.col = j;head.locationArray.push_back(loc);all_location_info.push_back(head);}}}。
基于代码得到所有像素结果,提取满足灰度一致关系的像素点坐标及其灰度值,即可从所有腹板信息的图像中提取出每一个单独的腹板。
步骤七、得到具有随机灰度的腹板区域后,对每一个独立的灰度区域,进行灰度化处理。本实施例中,R、G、B通道的权重分别为a 1=0.3、a 2=0.5、a 3=0.2,然后基于OpenCV函数threshold实现二值化threshold(img, img, 0, 255, THRESH_BINARY + THRESH_OTSU)。二值化后基于Canny算子计算得到图像的轮廓线,对应的像素宽度为1,处理后的效果如图7所示。
步骤八、得到每个腹板区域独立的轮廓区域后建立激光扫描点。
操作时,选择轮廓线任意位置为起始位置O(x,y),腹板面法矢方向为激光点的移动方向,竖直扫描原始三维图,移动的步长可以为0.1mm,激光点发出的光线垂直于腹板轮廓线指向外侧,发出的光线遇到平面即反射对应的距离为扫描距离L。竖直方向(由轮廓线上的点构成的线垂直于腹板的方向)的移动距离为零件三维图像中垂直于腹板方向的最大值H max ,当轮廓线上一个点扫描完成后沿着轮廓线顺时针方向移动激光点,对应的步距可以为0.2mm,直到将一个独立腹板的轮廓线全部扫描完成为止。激光点扫描方式与轮廓线之间的关系如图8所示。
步骤九、基于扫描得到的轮廓线不同位置不同高度的扫描数据,建立由轮廓线扫描数据生成的二维邻接关系图。
部分轮廓线二维邻接关系的表示如图9所示。其中,X轴表示不同的腹板区域轮廓线段;Y轴表示腹板法矢方向移动的距离;在图像中,相同轮廓线段对应的x值在不同移动距离y上具有不同的扫描距离(扫描距离为激光点至其前方对应特征面距离,移动距离为激光点沿着轮廓线移动的长度),用不同颜色的区域表示。二维坐标图像中,不同区域的颜色表示方式,以腹板区域轮廓线的垂线扫描距离开始,若扫描区域对应的槽为非开口槽,激光点线经过底角区域对应的扫描距离从0增加到R值(R为刀具的半径,一般为3mm),该区域采用统一的灰度表示,继续沿腹板的法矢方向移动激光点对应的扫描距离保持为定值R,对应的区域采用统一的灰度表示。沿腹板区域轮廓线的垂直方向开始,若为开口槽则对应的扫描距离为零表示。
步骤十、采用迭代循环处理的方式,将所有腹板区域进行激光点的三维扫描,并得到各个区域的二维邻接关系图。
步骤十一、建立邻接准则判断函数实现快速程编驱动(腹板面)加引导(侧壁面与圆角面)特征组的精细识别,对应的判断函数为:
;(13)
; (14)
其中,sum表示腹板封闭轮廓线中存在的直线段总数量;i表示直线段的标号;flag li 表示对应线段的状态标识;、/>、/>、/>分别表示腹板轮廓线l在第i段与y di 在二维坐标系中对应区域S中图像的R、G、B灰度值;di表示底角区域;up表示沿腹板法矢方向高于底角的区域。
步骤十二、基于上述判断函数,可以得到基于深度学习识别的所有腹板区域对应轮廓线的邻接关系,若一独立腹板轮廓线中存在一段区域其两侧各存在满足判断准则的面,即对应的面的flag为1其可与该面组成驱动加引导的组合关系。由此组合关系将识别结果返回映射到三维图像中,基于三维建模软件即可自动的得到NC(数字控制加工)刀轨程序。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。
Claims (6)
1.一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,包括如下步骤:
S1、在三维建模软件中载入零件模型图,设置其显示参数以区分前景和背景;
S2、采集不同视角的零件表面二维投影图像;
S3、基于深度学习的分割网络粗识别图像特征区域;
S4、针对粗识别结果得到的图像进行灰度化处理;
S5、采用最大间方差策略实现灰度化后图像的二值化,得到包含前景和背景的二值图像;
S6、通过二值图像的全局一阶梯度计算获取图像特征区域的轮廓图;
S7、基于随机灰度赋值法得到图像特征区域轮廓图中的独立轮廓区域;
S8、基于激光扫描得到每个独立轮廓区域的轮廓线扫描数据,并生成二维邻接关系图:选择轮廓线任意位置为起始位置,腹板面法矢方向为激光点的移动方向,竖直扫描原始三维图,激光点发出的光线垂直于腹板轮廓线指向外侧,发出的光线遇到平面即反射对应的距离为扫描距离L;竖直方向的移动距离为零件三维图像中垂直于腹板方向的最大值H max ,当轮廓线上一个点扫描完成后沿着轮廓线顺时针方向移动激光点,直到将一个独立腹板的轮廓线全部扫描完成为止;二维邻接关系图中,X轴表示不同的腹板区域轮廓线段;Y轴表示腹板法矢方向移动的距离;针对每个腹板区域独立的轮廓区域后建立激光扫描点,扫描轮廓线对应的三维图;
S9、建立邻接准则判断函数实现快速程编驱动面与引导面特征组的精细识别;
对应的判断函数为:
;
;
其中,sum表示腹板封闭轮廓线中存在的直线段总数量;i表示直线段的标号;flag li 表示对应线段的状态标识;、/>、/>、/>分别表示腹板轮廓线l在第i段与y di 在二维坐标系中对应区域S中图像的R、G、B灰度值;di表示底角区域;up表示沿腹板法矢方向高于底角的区域;
S10、基于精细识别结果自动生成加工刀轨程序:基于上述判断函数,可以得到基于深度学习识别的所有腹板区域对应轮廓线的邻接关系,若一独立腹板轮廓线中存在一段区域其两侧各存在满足判断准则的面,即对应的面的flag为1其可与该面组成驱动加引导的组合关系;由此组合关系将识别结果返回映射到三维图像中,基于三维建模软件即可自动的得到NC刀轨程序。
2.根据权利要求1所述的一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,零件表面二维投影图像为自动采集:在零件模型图中选择一特征面,以其对应的法矢方向为整个零件的视角方向,并截取屏幕图像;循环得到不同视角对应的特征面,使得零件在屏幕方向上的投影面积最大。
3.根据权利要求1所述的一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,零件表面二维投影图像为手动采集:利用三维建模软件中的工具,依次选择“图像”—“捕获”,在“选择模式”中框选出整个零件区域,点击“选项”,在弹出的“捕获选型”中选择像素,在“更多”中设置图像的分辨率宽和长;设置好捕获参数后,点击“捕获”按钮截取图像,并将截取的图像进行保存。
4.根据权利要求1所述的一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,所述的针对粗识别结果得到的图像进行灰度化处理包括:针对深度学习识别结果得到的图像,根据各个通道灰度值以加权求和的方式灰度化,得到灰度化后的图像。
5.根据权利要求1所述的一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,所述的采用最大间方差策略实现灰度化后图像的二值化包括:
对于大小为M*N包含sum个灰度级的图像,对应的计算式子为:
;
;
;灰度图中小于阈值像素点的平均灰度值为:
;
;
;
其中,p i 表示灰度级为i的像素数占总像素的个数;n i 表示灰度级为i的像素数;p background 表示像素点被分配为background的概率累加和;background与foreground分别表示前景与背景;threshold表示前景与背景分割的阈值;表示图像中所有像素点的灰度值方差。
6.根据权利要求1所述的一种三轴自动化程编特征识别与刀轨生成方法,其特征在于,所述的基于随机灰度赋值法得到图像特征区域轮廓图中的独立轮廓区域包括:
针对公式计算得到的图像特征区域轮廓图,采用全局像素种子点填充策略得到每一连通区域的语义信息;
每一填充区域采用随机赋值的方式进行处理,得到具有不同灰度的图像特征区域;
提取具有相同灰度信息的像素点并进行分类,从而得到独立轮廓区域;
针对每一个具有随机灰度的独立轮廓区域,分别进行灰度化和二值化处理,再基于Canny算子计算得到独立轮廓区域的轮廓线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311083327.5A CN116804865B (zh) | 2023-08-28 | 2023-08-28 | 一种三轴自动化程编特征识别与刀轨生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311083327.5A CN116804865B (zh) | 2023-08-28 | 2023-08-28 | 一种三轴自动化程编特征识别与刀轨生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116804865A CN116804865A (zh) | 2023-09-26 |
CN116804865B true CN116804865B (zh) | 2023-12-08 |
Family
ID=88079789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311083327.5A Active CN116804865B (zh) | 2023-08-28 | 2023-08-28 | 一种三轴自动化程编特征识别与刀轨生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116804865B (zh) |
Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0373854A2 (en) * | 1988-12-12 | 1990-06-20 | General Electric Company | Detecting internal structures contained within the interior region of a solid object |
WO1994024640A1 (en) * | 1993-04-13 | 1994-10-27 | Daniel Karron | System and method for surface rendering of internal structures within the interior of a solid object |
US6078642A (en) * | 1998-02-11 | 2000-06-20 | Analogice Corporation | Apparatus and method for density discrimination of objects in computed tomography data using multiple density ranges |
EP1062555A1 (en) * | 1998-02-11 | 2000-12-27 | Analogic Corporation | Computed tomography apparatus and method for classifying objects |
CN1541364A (zh) * | 2001-06-15 | 2004-10-27 | ��ѧ�뼼������˾@UNM | 具有选择运算符的数字电路 |
CN104360634A (zh) * | 2014-11-12 | 2015-02-18 | 南京航空航天大学 | 基于特征的蒙皮镜像铣削数控程序快速生成方法 |
TW201536390A (zh) * | 2014-03-27 | 2015-10-01 | Univ Nat Taiwan | 幾何結構解析方法、其系統及其電腦程式產品 |
EP3196863A1 (en) * | 2014-08-01 | 2017-07-26 | Shenzhen Cimc-tianda Airport Support Ltd. | System and method for aircraft docking guidance and aircraft type identification |
CN107767382A (zh) * | 2017-09-26 | 2018-03-06 | 武汉市国土资源和规划信息中心 | 静态三维地图建筑物轮廓线的自动提取方法及*** |
CN107845134A (zh) * | 2017-11-10 | 2018-03-27 | 浙江大学 | 一种基于彩色深度相机的单个物体的三维重建方法 |
CN109241628A (zh) * | 2018-09-08 | 2019-01-18 | 西北工业大学 | 基于图谱理论和聚类的三维cad模型分割方法 |
CN111310558A (zh) * | 2019-12-28 | 2020-06-19 | 北京工业大学 | 一种基于深度学习和图像处理法的路面病害智能提取方法 |
CN111815503A (zh) * | 2020-06-03 | 2020-10-23 | 上海交通大学 | 一种基于线激光扫描的螺柱特征参数测量方法 |
CN111815776A (zh) * | 2020-02-04 | 2020-10-23 | 山东水利技师学院 | 综合机载和车载三维激光点云及街景影像的三维建筑物精细几何重建方法 |
CN112258491A (zh) * | 2020-10-30 | 2021-01-22 | 国网河南省电力公司电力科学研究院 | 一种电力设备缺陷三维模拟方法 |
CN112308895A (zh) * | 2019-07-31 | 2021-02-02 | 天津工业大学 | 一种构建真实感牙列模型的方法 |
US11010943B1 (en) * | 2020-12-18 | 2021-05-18 | Ivan Bajic | Method and system for digital coloring or segmenting of multi-color graphics |
CN114882272A (zh) * | 2022-04-22 | 2022-08-09 | 成都飞机工业(集团)有限责任公司 | 一种航空制造件全角度投影图像面片属性融合分析方法 |
CN115204609A (zh) * | 2022-06-20 | 2022-10-18 | 成都飞机工业(集团)有限责任公司 | 数控加工工艺程编质量评价方法、装置、设备及介质 |
CN115345822A (zh) * | 2022-06-08 | 2022-11-15 | 南京航空航天大学 | 一种面向航空复杂零件的面结构光自动化三维检测方法 |
CN115439840A (zh) * | 2022-08-22 | 2022-12-06 | 成都飞机工业(集团)有限责任公司 | 一种航空件槽区识别方法、装置、设备及介质 |
CN115482522A (zh) * | 2022-09-21 | 2022-12-16 | 成都飞机工业(集团)有限责任公司 | 飞机结构件转角特征识别方法、装置、设备及存储介质 |
CN116229017A (zh) * | 2022-11-09 | 2023-06-06 | 南京大学 | 一种基于逆向投影的管件三维模型重建方法 |
CN116276938A (zh) * | 2023-04-11 | 2023-06-23 | 湖南大学 | 基于多零位视觉引导的机械臂定位误差补偿方法及装置 |
CN116542968A (zh) * | 2023-06-29 | 2023-08-04 | 中国铁路设计集团有限公司 | 一种基于模板匹配的钢筋智能计数方法 |
CN116596824A (zh) * | 2023-07-17 | 2023-08-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 多相机三维面模型融合方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070292033A1 (en) * | 2006-06-19 | 2007-12-20 | Chao-Wang Hsiung | Passive and interactive real-time image recognition software method |
US8306302B2 (en) * | 2008-09-29 | 2012-11-06 | Carestream Health, Inc. | Noise suppression in diagnostic images |
US8346007B2 (en) * | 2009-12-22 | 2013-01-01 | Carestream Health, Inc. | Noise suppression in cone beam CT projection data |
US9538336B2 (en) * | 2012-12-28 | 2017-01-03 | Trimble Inc. | Performing data collection based on internal raw observables using a mobile data collection platform |
CN106796246B (zh) * | 2014-02-24 | 2021-11-26 | 布鲁克纳米公司 | 自动扫描探针显微镜***中的精密探针部署 |
US9971335B2 (en) * | 2014-04-11 | 2018-05-15 | Thomas Marshall Tucker | Hybrid dynamic tree data structure and accessibility mapping for computer numerical controlled machining path planning |
JP6697501B2 (ja) * | 2018-03-26 | 2020-05-20 | ファナック株式会社 | 工作システム |
US11748942B2 (en) * | 2020-08-13 | 2023-09-05 | Siemens Mobility Pty Ltd | System and method for automatically generating trajectories for laser applications |
-
2023
- 2023-08-28 CN CN202311083327.5A patent/CN116804865B/zh active Active
Patent Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0373854A2 (en) * | 1988-12-12 | 1990-06-20 | General Electric Company | Detecting internal structures contained within the interior region of a solid object |
WO1994024640A1 (en) * | 1993-04-13 | 1994-10-27 | Daniel Karron | System and method for surface rendering of internal structures within the interior of a solid object |
US6078642A (en) * | 1998-02-11 | 2000-06-20 | Analogice Corporation | Apparatus and method for density discrimination of objects in computed tomography data using multiple density ranges |
EP1062555A1 (en) * | 1998-02-11 | 2000-12-27 | Analogic Corporation | Computed tomography apparatus and method for classifying objects |
CN1541364A (zh) * | 2001-06-15 | 2004-10-27 | ��ѧ�뼼������˾@UNM | 具有选择运算符的数字电路 |
TW201536390A (zh) * | 2014-03-27 | 2015-10-01 | Univ Nat Taiwan | 幾何結構解析方法、其系統及其電腦程式產品 |
EP3196863A1 (en) * | 2014-08-01 | 2017-07-26 | Shenzhen Cimc-tianda Airport Support Ltd. | System and method for aircraft docking guidance and aircraft type identification |
CN104360634A (zh) * | 2014-11-12 | 2015-02-18 | 南京航空航天大学 | 基于特征的蒙皮镜像铣削数控程序快速生成方法 |
CN107767382A (zh) * | 2017-09-26 | 2018-03-06 | 武汉市国土资源和规划信息中心 | 静态三维地图建筑物轮廓线的自动提取方法及*** |
CN107845134A (zh) * | 2017-11-10 | 2018-03-27 | 浙江大学 | 一种基于彩色深度相机的单个物体的三维重建方法 |
CN109241628A (zh) * | 2018-09-08 | 2019-01-18 | 西北工业大学 | 基于图谱理论和聚类的三维cad模型分割方法 |
CN112308895A (zh) * | 2019-07-31 | 2021-02-02 | 天津工业大学 | 一种构建真实感牙列模型的方法 |
CN111310558A (zh) * | 2019-12-28 | 2020-06-19 | 北京工业大学 | 一种基于深度学习和图像处理法的路面病害智能提取方法 |
CN111815776A (zh) * | 2020-02-04 | 2020-10-23 | 山东水利技师学院 | 综合机载和车载三维激光点云及街景影像的三维建筑物精细几何重建方法 |
CN111815503A (zh) * | 2020-06-03 | 2020-10-23 | 上海交通大学 | 一种基于线激光扫描的螺柱特征参数测量方法 |
CN112258491A (zh) * | 2020-10-30 | 2021-01-22 | 国网河南省电力公司电力科学研究院 | 一种电力设备缺陷三维模拟方法 |
US11010943B1 (en) * | 2020-12-18 | 2021-05-18 | Ivan Bajic | Method and system for digital coloring or segmenting of multi-color graphics |
CN114882272A (zh) * | 2022-04-22 | 2022-08-09 | 成都飞机工业(集团)有限责任公司 | 一种航空制造件全角度投影图像面片属性融合分析方法 |
CN115345822A (zh) * | 2022-06-08 | 2022-11-15 | 南京航空航天大学 | 一种面向航空复杂零件的面结构光自动化三维检测方法 |
CN115204609A (zh) * | 2022-06-20 | 2022-10-18 | 成都飞机工业(集团)有限责任公司 | 数控加工工艺程编质量评价方法、装置、设备及介质 |
CN115439840A (zh) * | 2022-08-22 | 2022-12-06 | 成都飞机工业(集团)有限责任公司 | 一种航空件槽区识别方法、装置、设备及介质 |
CN115482522A (zh) * | 2022-09-21 | 2022-12-16 | 成都飞机工业(集团)有限责任公司 | 飞机结构件转角特征识别方法、装置、设备及存储介质 |
CN116229017A (zh) * | 2022-11-09 | 2023-06-06 | 南京大学 | 一种基于逆向投影的管件三维模型重建方法 |
CN116276938A (zh) * | 2023-04-11 | 2023-06-23 | 湖南大学 | 基于多零位视觉引导的机械臂定位误差补偿方法及装置 |
CN116542968A (zh) * | 2023-06-29 | 2023-08-04 | 中国铁路设计集团有限公司 | 一种基于模板匹配的钢筋智能计数方法 |
CN116596824A (zh) * | 2023-07-17 | 2023-08-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 多相机三维面模型融合方法 |
Non-Patent Citations (9)
Title |
---|
Collision cross section prediction of deprotonated phenolics in a travelling-wave ion mobility spectrometer using molecular descriptors and chemometrics;Gonzales, Gerard Bryan;ANALYTICA CHIMICA ACTA;68-76 * |
Encounter trajectories for deep space mission ASTER to the triple near Earth asteroid 2001-SN263. The laser altimeter (ALR) point of view;Hussmann, H;ADVANCES IN SPACE RESEARCH;648-661 * |
Trajectory Generation and Tracking Using the AR.Drone 2.0 Quadcopter UAV;P.Vílez;2015 12th Latin American Robotics Symposium and 2015 3rd Brazilian Symposium on Robotics (LARS-SBR);73-78 * |
Trajectory planning for vascular navigation from 3D angiography images and vessel centerline data;A. Azizi;2017 International Conference on Manipulation, Automation and Robotics at Small Scales (MARSS);1-6 * |
基于加工特征的舵机壳体孔系特征数控加工编程技术;尚佳策;组合机床与自动化加工术;93-97+101 * |
基于特征的快速编程技术在飞机整体肋类零件中的应用;惠昕;制造技术与机床;41-45 * |
基于特征的数控编程技术在航空企业的推广应用;***;李迎光;刘旭;汤立民;;航空制造技术(第06期);54-59 * |
基于特征识别的整体壁板快速编程推理算法;乐毅;航空制造技术;80-83+92 * |
*** ; 李迎光 ; 刘旭 ; 汤立民 ; .基于特征的数控编程技术在航空企业的推广应用.航空制造技术.2016,(第06期),54-59. * |
Also Published As
Publication number | Publication date |
---|---|
CN116804865A (zh) | 2023-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109615611B (zh) | 一种基于巡检影像的绝缘子自爆缺陷检测方法 | |
CN112785643A (zh) | 一种基于机器人平台的室内墙角二维语义地图构建方法 | |
US20040155877A1 (en) | Image processing apparatus | |
EP1091320A2 (en) | Processing multiple digital images | |
CN106651752A (zh) | 三维点云数据配准方法及拼接方法 | |
CN114140679A (zh) | 一种缺陷融合方法、装置、识别***及存储介质 | |
CN110060338B (zh) | 一种基于bim模型的预制构件点云识别方法 | |
JP2020106467A (ja) | 欠陥検査装置、欠陥検査方法、及びそのプログラム | |
EP1081648B1 (en) | Method for processing a digital image | |
TWI773063B (zh) | 圖像產生系統及圖像產生方法 | |
RU2541917C2 (ru) | Устройство и способ для идентификации автора произведения искусства | |
CN107230203A (zh) | 基于人眼视觉注意机制的铸件缺陷识别方法 | |
CN109540925B (zh) | 基于差影法与局部方差测量算子的复杂瓷砖表面缺陷检测方法 | |
CN115239951B (zh) | 一种基于点云数据处理的墙面分割识别方法及*** | |
CN109978940A (zh) | 一种sab安全气囊尺寸视觉测量方法 | |
CN114882272A (zh) | 一种航空制造件全角度投影图像面片属性融合分析方法 | |
CN115953550A (zh) | 针对线结构光扫描的点云离群点剔除***及方法 | |
CN113506276A (zh) | 一种测量结构位移的标志物及方法 | |
CN116824273B (zh) | 一种航空制造件任意视角二维投影图像面片属性判断方法 | |
CN116804865B (zh) | 一种三轴自动化程编特征识别与刀轨生成方法 | |
CN112581424B (zh) | 一种光学元件表面与亚表面缺陷的分类提取方法 | |
CN116645351A (zh) | 一种复杂场景的缺陷在线检测方法及*** | |
US20200074685A1 (en) | System and method for representing and displaying color accuracy in pattern matching by a vision system | |
CN113096099B (zh) | 基于颜色通道合并的透水沥青混合料连通空隙识别方法 | |
CN116125489A (zh) | 一种室内物体三维检测方法、计算机设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |