固体推进剂广泛应用于火箭、导弹等发动机作为动力源，是一种含能的复合材料。由于火箭燃气发生器内压强较低，固体推进剂通常在贫氧状态下燃烧，燃速偏低，不易满足使用性能要求[1]。为了提高燃速，通常采用嵌银丝的方法。金属丝有很好的导热性，在药柱中嵌入金属丝可以有效将热量传递给药柱各部分[2]。为了增大燃面，缩短燃面爬升坡段，在药柱端面上沿银丝方向钻盲孔。盲孔采用复合加工系统加工，需要配合机器视觉系统对银丝位置进行识别，并使系统满足加工精度要求。

机器视觉系统是指通过图像传感器对图像进行摄取，并转换成数字信号，根据像素分布、亮度和颜色等信息，完成对目标特征的提取，实现检测、定位、监控等功能[3]。CMOS(互补金属氧化物半导体)和CCD(电荷耦合器件)作为半导体是目前最常见的两种图像传感器，两者各有优缺点。CCD相比于CMOS具有灵敏度高，成像品质好等优点，因此工业相机多数应用CCD传感器[4]。随着视觉技术的发展，机器视觉系统对图像的处理更加智能化，在工业应用中可以进行有效的缺陷识别、测量、环境建模等工作，具有效率高、精度高等特点，对工业制造领域中智能化生产起到了非常重要的作用。

本研究目的是利用 CCD相机捕捉图像，通过图层优化、边缘检测、灰度阈值处理、感兴趣区域提取、面积识别、递归贪心算法多种方式进行药柱银丝位置的识别与精确定位。本文研究的细银丝位置识别技术对其他类似图像中筛选特征点具有一定的参考价值。国内对于药柱金属丝的检测未有详细说明，文献[5-6]主要集中于药柱燃烧的速率和杀伤力，国外由于技术封锁对于药柱检测报道很少。

在 130 例患者中，CYP2C19*2、CYP2C19*3 基因型和等位基因频率分布见表1。两个位点均符合Hardy-Weinberg遗传平衡吻合度计算 （P>0.05），说明样本具有群体代表性。

1 药柱复合加工系统

复合加工系统是以通用部件为基础，可以同时拥有几种组合方式，配以按工件特定外形和加工工艺设计的专用部件和夹具组成的半自动或自动专用机床，生产效率比单工位加工系统高几倍至几十倍。由于通用部件已经标准化和系列化，可根据需要灵活配置，能缩短设计和制造周期。多工位加工系统兼有低成本和高效率的优点，在大批、大量生产中得到广泛应用，并可用以组成自动生产线,目前技术难点在于如何保证多工位加工系统的加工精度，既能实现加工自动化，又保证加工的精度，需要精确的识别定位。

1.1 复合加工系统的功能及要求

药柱复合加工系统具有整主形端面、银丝位置自动识别处理与精确定位、盲孔加工、孔内药屑清理、孔深及直径自动测量、刀具温度检测等功能。复合加工系统能根据相机检测的银丝位置，自动找准，对中心，钻孔定位精度：不大于Ф0.35 mm；孔径中心轴线与工件中心轴线平行度：≤0.2 mm。钻孔的定位精度主要由银丝位置的识别精度保证，钻头轴线与药柱银丝轴线平行度由系统设计与安装调试保证。机械臂、复合加工系统及打孔刀具刚度是保障此项精度的重要因素。

1.2 复合加工系统的组成

设平面任意位置椭圆方程为

图1 复合加工系统示意图

银丝位置的识别与定位由高精度彩色CCD工业相机完成，安装在整合加工模块的安装板的固定位置。CCD相机具有图像捕捉快速，图像清晰，相机位置摆放灵活，检测方法精度易控制等特点，支持快速连续曝光，动态图像无拖影，能够在任何时刻完成图像捕捉任务。彩色CCD相机有利于计算机从色度上筛选银丝所在区域，同时，自适应的处理图像亮度以及自动增益和自动曝光效果，在原图的基础上优化了对用户显示的效果。同时操作人员对加工准确与否的辨识度较高，通过改写程序可适应多种类型的工件。

方案规定在0～600m井段用PTB-64钻具钻进，用PTB—480钻进。如果说利用PTB—64钻进没有大的问题、直径508mm管柱下放和注浆跟设计没有误差的话；而采用PTB—480钻进时，则发现顶角在600m深处增大 1°15'， 在 880m 处达 4°15’， 即每 100m 为 1°30’；这需要在钻井钻进工艺方面做大的改变，采取了停止使用 PTB消除下降梯度不大于每 100m为 1°30’2总曲度的方案。这一任务通过钻进直径295mm先导孔完成。成套工具中采用了偏心距17mm的偏心接头，接头位于钻头上方。而后，井眼直径扩至480mm。

工业变焦镜头采用全新超低失真、高分辨率设计，具有清晰度高，倍率可变，光学无畸变，高对比度等特点。采用焦阑法，测量误差少，能配合百万像素以上1/3、1/2、2/3工业相机使用，用于高精确检测、高精度尺寸测量、目标位置定位等机器视觉系统。

2 银丝位置多重算法识别技术

因此得方程组：

图2 药柱图片

由于药柱端面上银丝极细，肉眼不易观察，为准确找到银丝位置，使用CCD相机获取图像，通过数字图像处理系统，得到当前状态下银丝的位置，需要对药柱断面的表面形貌进行图像处理。

视觉系统采用六重识别算法，包括图层优化、边缘检测、灰度阈值处理、感兴趣区域提取、面积识别、递归贪心算法，利用此方法可以在银丝位置快速识别的前提下保证检测结果准确无误。

2.1 图层优化处理

基于canny算法的边缘检测实现过程如下：

Gray=R+B+C

(1)

三基色图层优化处理及灰度处理后如图3所示。

图3 视觉图层优化处理示意图

2.2 边缘检测

由于药柱断面的表面形貌的图像处理受表面粗糙度、光照环境等因素影响，故选择合适的边缘检测算子对图像进行边缘特征的提取，并且尽可能减少由噪声引起的伪边缘的干扰。本系统采用算法相对简单，应用广泛的Canny算子进行边缘检测[7]。

由于自然光是混合光源，通过自然打光下的图像提取三基色图层(RGB)，根据银丝在不同图层中的明显程度，对灰度图像进行加权求和，得到最有效的银丝识别图像。加权求和即将RGB三个分量分别作为该点的灰度值，然后相加得到一个灰度值：

1) 利用高斯滤波对图像进行平滑降噪处理

采用高斯函数对药柱断面图像进行处理，函数形式如下。

(2)

σ为高斯曲线标准差，影响图像中检测的边缘平滑程度。

若固守男女授受不亲之礼，则不应施以援手救嫂子；此时施以援手，则是权的体现。而孟子懂得通权达变，用权变思维处理嫂溺问题。在孔子看来，当处在道德两难的特殊境遇时，道德主体可以暂时牺牲乃至违背较低层次的道德规范，以保证更高层次的道德规范的实现，而这种行为仍具有真正的道德价值，这就是“权”。在处理两难问题时，就需要灵活变通，不能固守于经。权变是人的道德选择，也是道德实践的前提。

石蜡组织切片常规脱蜡和水化，加3%过氧化氢溶液室温孵育10 min以阻断内源性过氧化物酶；高温高压抗原修复；PBS缓冲液冲洗3×3 min；滴加一抗(鼠抗人单隆抗体CerbB-2)，4 ℃冰箱过夜；室温下孵育20 min，PBS缓冲液冲洗3×3 min，滴加二抗(生物素化山羊抗小鼠IgG)，室温孵育20 min；PBS缓冲液冲洗3×3 min，滴加新配制的DAB显色液，显微镜下观察，阳性为棕黄色颗粒状、定位于细胞膜；自来水冲洗，苏木精复染，脱水透明，封片、观察结果。

图像I(i, j)的边缘包括梯度幅值M(i, j)和梯度方向θ(i, j)，可采用2×2邻域一阶偏导的有限差分进行近似计算。本文取x、y方向一阶近似偏导模板如下：

(3)

3.内容不具体，中心不明确。内容不具体，这是很多学生的一个共同问题。写作时无从下手，不知道写什么，以至于文章空洞苍白，泛泛而谈，不知所言。文章中心具体表现为作者的感情倾向和主题思想，不论是直露的还是含蓄的，都应明确突出、贯穿始终。但是在文章中心的表现上，很多学生都有不同程度的问题。

(4)

(5)

3) 对梯度幅值进行非极大值抑制

非极大值抑制(NMS)是指保留局部极大值，将非极大值的像素点的灰度值设置为0，这样可以有效去除大多数非边缘点，达到细化边缘的目的。

非极大值抑制的实现是利用梯度方向，与相邻像素梯度赋值进行比较，即以梯度值为M(i, j)的点(i, j)为中心，与3×3邻域内沿梯度方向θ(i, j)的元素比较，如果M(i, j)比沿梯度方向相邻元素大，则将其设为候选边缘点；若相反则将其灰度值设为零，这样可以得出候选边缘图像C(i, j)。

（5）缸筒表面外观质量 激光淬火后，缸筒淬火区无明显的氧化脱碳现象，表面粗糙度值较低，经磁粉无损检测后表面无微小裂纹现象。并经专业的检测手段检测缸筒淬火区没有变形缺陷产生。

4) 检测和连接边缘

31年后的2015年，昆明新空港长水国际机场近旁的小哨乡牧场，云南省草地动物科学研究院院长黄必志和他的同事来到草场的牛舍内，一头头体格壮硕的牛围了过来，“这就是‘云岭牛’的种牛，它们1/4的血缘来自十九世纪八十年代的莫累灰牛。”黄必志说。

Canny算法采用双阈值法，即选取高阈值T1和低阈值T2分别与非极大值抑制后图像进行比对。若像素点梯度幅值大于高阈值，则将该点确定为边缘点，若梯度幅值低于高阈值，则将该点确定为非边缘点；若梯度幅值在高阈值和低阈值之间，但是在该点邻域像素中有已被确认的边缘点，则将该点也确定为边缘点，否则为非边缘点。进行双阈值法检测和连接的边缘三维效果图及边缘示意图如图4。

梯度幅值和梯度方向计算公式如下：

图4 边缘检测处理示意图

2.3 灰度图像阈值处理

若要从药柱断面中识别银丝位置，需要对药柱断面的表面形貌进行分析。为了提取银丝位置的特征点，需要对药柱断面图像进行处理。通过将灰度图进行二值化处理，即通过给定阈值，以该阈值为门限将灰度图转换成黑白图像[9]。

(6)

式中，T为给定阈值，由于是批量生产，可凭经验选取。

2) 计算梯度的幅值和方向

阈值T取10后的灰度阈值处理图如图5所示。由于银丝所在位置与其他部分会显示出不同的灰度值，可以通过筛选类似银丝所表现的像素点确定银丝位置。

图5 灰度阈值处理示意图

2.4 ROI感兴趣区域提取

通过提取图像中边缘样本坐标点，采用最小二乘椭圆拟合的方法，求得中心点位置坐标。

复合加工系统包括药柱夹具、直角坐标驱动单元、整形加工模块。其中整形加工模块整合了CCD相机、打孔钻头、铣刀、毛刷清根、高压气针吹屑、内窥镜探测、磁滞尺测深、红外温度传感器，可实现系统要求的多种功能，满足自动化要求。 复合加工系统结构如图1所示。

Ax2+Bxy+Cy2+Dx+Ey+F=0

(7)

Pi(xi,yi)(i=1,2,…,n)为椭圆上的测量点，依据最小二乘法，所拟合的目标函数为

(8)

欲求中心点坐标即：

(9)

在药柱定长切割结束后，利用工业相机完成药柱断面拍照，完成端面上圆周阵列的6个银丝位置的精确识别，获得银丝相对于钻头的位置。未处理药柱图片如图2所示。

(10)

结合约束条件可以求得方程系数A,B,C,D,E,F的值。结合椭圆的几何知识，可以计算出椭圆的5个参数：位置参数(θ,x0,y0)以及形状参数(a,b,)，即可求得中心坐标。

当先导阀口未开启时，阻尼孔R1与R2中无油液流动，主阀芯上下两腔压力相等，主阀芯在复位弹簧的作用下保持关闭。当系统压力增大至使先导阀口开启时，压力油经R1和先导阀口流回油箱。因阻尼孔R1、R2的作用，主阀芯上腔的压力p2小于主阀芯下腔的压力p0，由压差(p0-p2)形成的液压力克服主阀弹簧力使主阀芯向上移动，主阀口开启，实现溢流作用。

由于银丝位置分布是圆周阵列的形式，故感兴趣区域可以设置为以药柱边缘计算得到的中心点作为中点的环形区域，其外内圆直径为R±p，R为药柱上银丝可能所在直径，但由于药柱生产时银丝位置会有一定误差，故留出一定裕度p，以便使感兴趣区域将6个银丝位置全部包含进去。

其他区域通过漫水填充法进行填充，将灰度都设置为0，这样可以有效生成掩码区域，通过只处理掩码区域，可以得到感兴趣区域。

4.HIV职业暴露后预防性用药原则：（1）治疗用药方案：首选推荐方案为TDF/FTC+RAL或其他INSTIs；根据当地资源，如果INSTIs不可及，可以使用PIs如 LPV/r和DRV/r；对合并肾脏功能下降者，可以使用AZT/3TC。（2）开始治疗用药的时间及疗程：在发生HIV暴露后尽可能在最短的时间内（尽可能在2 h内）进行预防性用药，最好不超过24 h，但即使超过24 h，也建议实施预防性用药。用药方案的疗程为连续服用28 d。

我国金融利率市场遵循“先农村后城镇，先贷款后存款”的规律，在2015年正式放开城镇居民的存款利率，但对贷款利率仍然实行严格的管控。目前，农村信贷资金紧张，农村金融改革乏力，严重影响了乡村振兴战略的实施。要逐步放宽农村信贷利率浮动区间，满足农业建设对信贷资金的需求。同时要建立农村金融市场风险投资机制，完善农村金融财政补贴政策，加强农业担保基金管理，设立农村共管基金，从而增强农村金融市场的资金周转量[6]。

当年8月的一天中午，廉小花带着小云天午睡。小云天醒来后，扒开廉小花的衣服吸吮她的乳头。廉小花被惊醒后推开他，小云天瘪着小嘴哭。想到他自出生5个月后就没见过妈妈，廉小花有些不忍，只好又把乳头放到他嘴里。母性就这样在她心中开始发芽。

通过ROI的提取，能有效提高计算机处理速度，并且避免干扰产生的误差。药柱断面图像中感兴趣区域提取后如图6所示。

图6 感兴趣区域提取示意图

2.5 特征区域面积识别

特征区域面积识别是通过比对可疑银丝位置特征区域的面积，筛选面积部分过小或过大的区域，排除绝对不可能是银丝的点。

由于图像是由一个一个像素组成的，因此在对图像逐行采用区域标记算法将各特征区域进行标记，区域标记法[11]是把连续区域附上相同标记，常见的有四连通域和八连通域标记算法。本文采用八连通域算法[12]，在对特征连通区域进行标记的同时对各特征区域中像素个数进行统计。然后对各区域像素个数进行比对，去除不可能是银丝位置的特征区域。设银丝位置特征点像素个数范围是[n1,n2]，像素个数不在该范围内的标记区域的灰度值设为0。干扰特征消除后的效果图如图7。

图7 干扰特征消除后示意图

2.6 结合递归和贪心算法确定银丝坐标

药柱端面图像最终的处理包括图像检索及坐标值计算。图像检索的方式采用以区域块为单位进行图像扫描，区域块是指利用递推算法合并四邻域(上下左右)相同灰度值像素从而得到若干一定尺寸的区域块[13]。检索出包含所有疑似银丝位置的区域块，区域块中心位置坐标为(xi,yi)。

假设图像I中还存在n个疑似银丝位置点，从上到下的区域块依次标记为i=1,2,…,n，利用贪心算法，即通过解局部最优解来达到整体最优解的方法，比对两两区域块中心坐标之间的距离，在图像I中寻找6个两两之间近似等距的圆周阵列的银丝位置点(阵列中相邻银丝位置之间距离设为l±a)。如符合银丝相对位置关系，将此两点标记为可能是银丝的点，如果某一点被重复标记3次以上，则认为它就是银丝所在位置，反之则将该区域块的灰度值全设置为0，最后提取最优结果，并通过坐标变换得到世界坐标系下银丝的绝对位置，以便进行盲孔的精确定位加工。

2015年下半年，煤炭行业形势更加严峻，煤炭企业经济运行面临的压力进一步增大。在这种条件下，对风选项目的推进工作更加迫切，进一步优化项目方案并缩减项目总投资成为风选项目实施关键所在。

通过递归贪心算法得到的最终银丝位置示意图如图8。

图8 最终处理后的银丝位置示意图

3 结论

本文采用CCD相机，利用图层优化处理、边缘检测、灰度阈值处理、感兴趣区域提取、面积识别、递归贪心算法六重算法成功对整形后的圆柱型药柱断面进行银丝位置识别与定位，有效去除了非银丝位置的特征区域。本文研究的药柱银丝位置识别的方法识别速度较快，可以满足批量生产的要求。

参考文献：

[1] 于宁,杨涛,刘巍,等.嵌银丝端燃药柱燃气发生器理论性能分析[J].弹箭与制导学报,2009,29(1):180-183.

[2] 张有为,王晓宏,杨举贤.嵌金属丝端燃药柱燃烧过程的数值研究[J].推进技术,2007,28(1):4-8.

[3] 卞正岗.机器视觉技术的发展[J].中国仪器仪表,2015(6)：71-75.

[4] 房加强,何建萍,王付鑫,等.视觉传感技术的应用研究与发展趋势[J].电焊机,2013,43(4):46-50.

[5] 王元有,王新华.含金属丝端燃药柱的燃烧机理[J].兵工学报,1980,1(2):42-53.

[6] 李凤生.嵌涂层金属丝药柱的稳态燃烧模型[J].南京理工大学学报,1984(3):65-74.

[7] 李牧,闫继红,李戈,等.自适应Canny算子边缘检测技术[J].哈尔滨工程大学学报,2007,28(9):1002-1007.

[8] JOHN C.A computationai approach to edge detection[J].IEEE Transaiction Pattern Anaiysis， 1986(8)：679 - 698.

[9] 李平,张廷安,汪秉宏,等.灰度阈值对图像分形特征参数提取的分析[J].东北大学学报(自然科学版),2006,27(1):57-60.

[10] 许浩均,黄茜,梁兆敏,等.基于MCSA的感兴趣区域定位方法[J].现代电子技术,2015,38(17):103-106.

[11] 曹长虎,李亚非.一种二值图像连通区域标记快速算法[J].科学技术与工程,2010,33(10):8168-8171,8180.

[12] 李仪芳,刘景琳.基于连通域算法的区域测量[J].科学技术与工程,2008,9(8):2492-2494.

[13] 卢敏,伍铁军.一种基于区域特征的图像检索方法及实现[J].图学学报,2004,25(1):42-46.