第４期 微处理机 ＮＯ．４ ２０１７年８月 ＭＩＣＲＯＰＲ０ＣＥＳＳＯＲＳ Ａｇｕ．，２０１７ 基于ＣＴ技术的图像边缘检测及ＯｐｅｎＣＶ实现术 刘洲洲 ，闫 伟 （１．西安航空学院电子工程学院，西安陕西７１００７７， ２．延安大学附属医院东关放射科，延安陕西７１６０００） 摘要：计算机断层技术（ＣＴ）的出现，是人类医学史上的一大突破，标志着现代医学正式跨入数 字化时代。ＣＴ图像的三维重建是医学影像学中的一个重要领域，同时口腔ＣＴ又是使用率最高的领 域之一。作为图像分割的一种方法，利用开源计算机视觉库（ＯｐｅｎＣＶ）具有多平台支持，多语言接口， 强移植性，快处理速度等优点，使用不同的边缘提取算子提取ＣＴ图像的轮廓并进行比较，通过ＣＴ机 扫描而得的图像进行双边滤波器平滑滤波，边缘保留，Ｃａｎｎｙ算子边缘提取等图像预处理。在算法实 现上提出了方便调试的轨迹条，可以很好地对各种参数进行可视的实时调节，便于找出最优方案。 关键词：ＣＴ技术；数字图像处理；边缘检测；开源计算机视觉库；平滑滤波；边缘提取 ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００２－２２７９．２０１７．０４．０２１ 中国分类号：ＴＮ３９３ 文献标志码：Ａ 文章编号：１００２—２２７９一（２０１７）０４—００８７—０４ Ｉｍａｇｅ Ｅｄｇｅ Ｄｅｔｅｃｔｉｏｎ ａｎｄ ＯｐｅｎＣＶ Ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ， Ｂａｓｅｄ ｏｎ ＣＴ Ｔｅｃｈｎｏｌｏｇｙ Ｌｉｕ Ｚｈｏｕｚｈｏｕ ，Ｙａｎ Ｗｅｉ ．Ｓｃｈｏｏｌ ｏｆ Ｅｌｅｃｔｒｏｎｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｘｉ’帆Ａ ｅｒｏｎａｕｔｉｃａｌ Ｕｎｉｖｅｒｓｉｔｙ，Ｘｉ’ａｉｒ ７１００７７，Ｃｈｉｎｏ， ２．Ｙｏｎ’觎Ｕｎｉｖｅｒｓｉｔｙ Ａｆｉｆｌｉａｔｅｄ Ｈｏｓｐｉｔａｌ Ｒａｄｉｏｌｏｇｙ，Ｙｏｎ’帆７１６０００，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅ ａｄｖｅｎｔ ｏｆ ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ ｔｅｃｈｎｉｑｕｅ（ＣＴ）ｉｓ ａ ｂｒｅａｋｔｈｒｏｕｇｈ ｏｆ ｈｕｍａｎ ｍｅｄｉｃａｌ ｈｉｓｔｏｒｙ，ｍａｒｋｉｎｇ ｔｈｅ ｍｏｄｅｌＴｌ ｍｅｄｉｃｉｎｅ ｏｆｆｉｃｉａｌｌｙ ｅｎｔｅｒｅｄ ｔｈｅ ｄｉｇｉｔａｌ ａｇｅ．ＣＴ ｉｍａｇｅ ３Ｄ ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ ｉｓ ａｎ ｉｍｐｏｒｔａｎｔ ｆｉｅｌｄ ｏｆ Ｍｅｄｉｃａｌ Ｉｍａｇｅｏｌｏｇｙ，ａｎｄ ｏｒａｌ ＣＴ ｉｓ ｏｎｅ ｏｆ ｔｈｅ ｆｉｅｌｄｓ ｈａｖｉｎｇ ｈｉｇｈｅｓｔ ｕｔｉｌｉｚａｔｉｏｎ ｒａｔｅ．Ａｓ ａ ｍｅｔｈｏｄ ｏｆ ｉｍａｇｅ ｓｅｇｍｅｎｔａｔｉｏｎ，ｕｓｉｎｇ ＯｐｅｎＣＶ ｈａｓ ｓｏｍｅ ａｄｖａｎｔａｇｅｓ ｓｕｃｈ ａｓ ｍｕｌｔｉｐｌｅ ｐｌａｔｆｏｒｍ ｓｕｐｐｏｒｔｉｎｇ， ｍｕｈｉｌｉｎｇｕａｌ ｉｎｔｅｒｆａｃｅ，ｓｔｒｏｎｇ ｐｏｒｔａｂｉｌｉｔｙ，ａｎｄ ｑｕｉｃｋ ｐｒｏｃｅｓｓｉｎｇ．Ｕｓｉｎｇ ｄｉｆｆｅｒｅｎｔ ｅｄｇｅ ｄｅｔｅｃｔｉｏｎ ｏｐｅｒａｔｏｒｓ ｔｏ ｅｘｔｒａｃｔ ｔｈｅ ｏｕｔｌｉｎｅ ｏｆ ＣＴ ｉｍａｇｅ ａｎｄ ｃｏｍｐａｒｅ，ｔｈｅ ｉｍａｇｅ ｐｒｅｐｒｏｃｅｓｓｉｎｇ，ｓｕｃｈ ａｓ ｔｏｗ－ｓｉｄｅ ｆｉｌｔｅｒ ｓｍｏｏｔｈｉｎｇ ｉｆｌｔｅｒｉｎｇ，ｅｄｇｅｓ ｐｒｅｓｅｒｖ！ｎｇ，ａｎｄ Ｃａｎｎｙ ｏｐｅｒａｔｏｒ ｅｄｇｅ ｅｘｔｒｃａｔｉｏｎ，ｉｓ ｅｘｅｃｕｔｅｄ ｔｏ ｔｈｅ ｉｍａｇｅｓ ｆｒｏｍ ＣＴ ｓｃａｎｎｉｎｇ． Ａ ｔｒａｃｋ ｂａｒ ｅａｓｙ ｔｏ ｔｅｓｔ ｉｓ ｐｕｔ ｆｌｏｒｗａｎｄ ｆｏｒ ａｌｇｏｒｉｔｈｍ ｉｍｐｌｅｍｅｎｔａｔｉｏｎ，ｗｈｉｃｈ ｃａｎ ｍａｋｅ ｖｉｓｉｂｌｅ ｔｒｕｅ ｔｉｍｅ ｃｏｎｔｒｏｌ ｔｏ ａ ｖａｒｉｅｔｙ ｏｆ ｐａｒａｍｅｔｅｒｓ ｔｏ ｃｏｎｖｅｎｉｅｎｔｌｙ ｆｉｎｄ ｏｐｔｉｍａｌ ｓｃｈｅｍｅ． Ｋｅｙ ｗｏｒｄｓ：ＣＴ ｔｅｃｈｎｏｌｏｇｙ；Ｄｉｇｉｔａｌ ｉｍａｇｅ ｐｒｏｃｅｓｓｉｎｇ；Ｅｄｇｅ ｄｅｔｅｃｔｉｏｎ；ＯｐｅｎＣＶ；Ｓｍｏｏｔｈｉｎｇ ｆｉｌｔｅｒ； Ｅｄｇｅ ｅｘｔｒａｃｔｉｏｎ １ 引 言 的地方。边缘是图像的重要信息，是图像分割的一个 重要特征。对人来说，图像的边缘部分是首先映人脑 图像分割作为图像分析、理解的基础，在诸多领 海的，可以区分物体与物体之间的界限，边缘是一个 域具有广泛的应用，如基于内容的图像检索、机器视 区域的结束，又是另一个区域的开端，利用边缘的特 觉、指纹识别以及生物医学图像处理方面的病变检 性就可以对图像进行有效的分￣Ｊ［３１。 测、军事图像处理方面的地形匹配和目标锁定［１＿２ｌ。 边缘检测是基于边界识别的一类图像分割方法。边 ２边缘检测 缘检测的基本思路是通过寻找图像中不同部分的边 计算机处理得出的边缘并不一定就是图像的真 界，从而达到突出边缘。边缘是指像素灰度有阶跃性 正边缘，只是在图像信息方面引起计算机注意的部 变化的区域，在图像中最直观的边界就是颜色突变 分，具体参数要结合具体情况进行适当的整改凹。由 基金项目：国家自然科学基金资助项目（６１６０１３６５） 作者简介：刘洲洲（１９８１～），男，山西省运城市人，博士生，研究方向：智能优化算法。 收稿日期：２０１６—０９—２６ ・８８． 微处理机 于实际物体的信息是三维的，而得到的图像数据是 着边缘方向的像素变化比较小，较为平缓，垂直于边 二维，从三维降到二维的投影必然会造成信息的缺 缘方向的像素变化较大，边缘上的这种信息变化可 失，再加上环境和噪声对图像采集带来的影响，使得 以用微分算子检测出来，通常使用一阶或者两阶导 很多边缘不能被有效的识别，识别出的边缘也不一 数来进行边缘检测。边缘大致可分成两种：一种是阶 定为真边缘。方向和幅度是边缘的两个重要属性，沿 跃状的如图１、３、５；另一种是屋顶状的如图２、４、６。 图１ 阶跃函数原函数 图２屋顶函数原函数 图３ 阶跃函数一阶导数 图４屋顶函数一阶导数 图５阶跃函数二阶导数 图６屋顶函数二次导数 边缘检测首先要对图像进行滤波，因为边缘检 数构成的。Ｊｏｈｎ Ｃａｎｎｙ在１９８６年提出Ｃａｎｎｙ算子，它 测的算法主要是对图像信息进行一阶导数或者二阶 与使用局域窗口梯度的Ｓｏｂｅｌ和Ｌａｐｌａｃｅ算子不同， 导数，但是导数对噪声很敏感，如果图像中有噪点， 对噪声的敏感度不强，是先对图像进行平滑然后求导 那边缘检测会放大这种噪声效果，通常显示为发光 数的，Ｇａｕｓｓｉａｎ函数与图像的卷积可以达到抗噪的目 亮点。因此必须采用滤波器来模糊噪声带来的影响。 的，导数又可以对边缘进行检测。所以Ｃａｎｎｙ算子在 然后是对图像进行增强，增强算法可以将邻域内比 实际进行图像处理的时候效果比较好，是最常用的边 较明显的灰度变化较为明显地表现出来。经过增强 缘检测手段之一。由于Ｇａｕｓｓｉａｎ函数是一个对称函 的图像在邻域中可能会有一些点的幅度值比较大， 数，图像的边缘在边缘方向是对称的，在垂直于边缘 但是在大部分隋况下，这些幅度值较大的点并不是 的方向上是反对称的，所以Ｃａｎｎｙ算子对边缘两侧的 所需要的边缘，在编程实现中往往通过阈值法来进 变化特别敏感，而对平行于边缘方向的变化不敏感。 行选择性取舍。 设二维Ｇａｕｓｓｉａｎ函数： １ ３ Ｃａｎｎｙ算子 ＝ ｅｘｐ（一 在实际边缘提取的应用中，对噪声的抑制和对 最优阶跃边缘检测算子是基于卷积ＶＧ＋，（ ， ）， 边缘的提取是很难同时完全满足的嘲，例如高精度零 边缘强度为ＩＶＧ ｆ（ｘ，圳，边缘方向为 件和ＣＴ等对细节边缘有着极高的要求的实例在普 ＶＧ通边缘检测步骤下较难实现。因此选择适当的边缘 ｐ ＊丽ｆ（ｘ ，ｙ） （２） 检测算子、步骤和参数可以使得到的结果接近所需 基于Ｇａｕｓｓｉａｎ函数的可分性，可以把 ６的两 要的边缘提取目的。 个滤波卷积模板分解成两个一维行列滤波器： 有一个线性的算子可以在边缘检测和抗噪声中 ＯＧ（ｘ，ｙ）＝ｋｘｅｘｐ（～ｘ＿ Ｌ） ｅ （一取得一个平衡，能很好地在提取边缘的同时保留一些 －￣－ｒ）＝ｏｘ ２ ｄ ２＾（ ） （ ）（３） 细节，这个算子就是Ｇａｕｓｓｉａｎ函数的一阶导数。而 ａＧ（ｘ ，ｙ）Ｃｍａｎｙ边缘检测算子就是由Ｇａｕｓｓｉａｎ函数的一阶导 ＝ ｅ （一番）ｅ ～ ： ） ）（４） ４ 止 儿 一寺）． ＿）＝．Ｊｋ． ｘ 砉 ） ｃ州一 ¨＝ ｅｘ 熹） 刮 ｈ，（ １＝ｘｈ （． ），ｈ．（１，）＝ｖｈ （１，） 、ｆ ８ （５） ／ｒ＾、… ｒｆ＿ｌＧ（一 ）址边缘。Ｉ 数； （ ）是带宽乃 的低 通滤波器脉冲响 ； ／Ｊ：（ ｕｓｓ 噪， ｉ门均力 差，ＳＮＲ他越人，提取的边缘质 越高 （６） 定化精度凇ｆＪ！ｌＪ： Ｌ＝ 可把』ｉＩｊ个一维行列滤波器 ／’（ ． ，）进仃岱币ｊ！： ：．ｌｒ Ｇ （＿ｌＹ）／／（ ）Ｕｘｌ ｊ （９） ），Ｅ ： （－Ｌ （１ ，（ ¨ （７） 。 ｉｌ ｒｌｘ）ｄｘ 令 ．Ｇ （— ）为Ｇ（一 ）ｆｆ 阶导数，Ｄ，（ Ｙ）为ｈ（－Ｙ）的 ４（ｉ 厢．．．．，ｆ』ｌｆ，／）＝ａｒｃｔａｎ ， 一阶导数；Ｌ为川术 价定位精确度的 ．Ｌ ｆ 地 单边缘响应准则： 确保对于边缘像素只仃 个单独的响心，钟：Ｊ 的脉７Ｉ｛Ｊ响 导数的零交义， ａ（ｉ，ｉ） 爪边缘强度， （，， ）址垂ｎ边缘的力‘ 向一 大，定他的精确度越岛 平 ＣａＩＩＩ１Ｖ提…他Ｊ｝ｊ找川ｊ Ｇ 算子与Ｉ纠像，， （ ． ）的 卷积在边缘力‘向㈨，Ｊ最大值的方法来定位边缘这 种疗法可以以任意一点的Ａ（ｉ ，）值是否为“（ ）邻城 均距离 满 ： 内的最大值来划断陔点是 为边缘 、ｆ，一个 一 满 足 个条什ｌｊ、ｆ．可’ｆＪ芰认为足冈像的一个边缘＾ ４边缘提取的评价 Ｃａｔ ｖ给…ｒ 个评价边缘伶测性能ｆ，ｆ，Ｊ指抓 Ｊ： 一外源汁ｌ ：机 觉 （（）ｆ １ Ｓ（１Ｉ… ＣｏｌｌｌｌｌｌｌＩ¨．－ｉ— 足低火 慨 昝． 幢多的包含虹实边缘，减少将 ； 址 ＳｉＯｉｑ １５１）ｒⅢ’ｙ，（） １（： ）ｌ 址 特尔公刊 ２０００ ｆｌｑｆｉ！ Ｉ的一个软什 多语 包 ｒｉ下多针埘 像处 （） ｎ（：ｖ具 多 ‘ 的 支¨ ， 库 数，＿Ｊ『以Ａ接 ・Ⅵ川足边缘的部分０ ｌ为边缘的错误结 ； 址Ｉ‘￣’，Ｊ Ｉ＂１ Ｊ确度， 得…的检测结 ｊ、 ｉ发都存 Ｆ确的实际边 单像素宽．刈边缘 一个唯一的响 ， 模 边抖 ｛－ 找ｆ１１最优边界针埘以Ｌ 个指怀，ＩＪＩｌＩ川ＣⅢ…、，捉… 了Ｈ ｒ一阶微分泌 芝器的边缘愉 ＩＪ！ｉＪ￣：ｔｑｆｉ＂准则：ｆ 【 】匕准ＩＪ！ｌＪ： Ｉ，姒ｆ｛； Ｉ＇ｔ－：．怏处理速度等优点 所他川１ 的是ＣＰ 乃ｌｎｔｅｌ¨ ｔ ｉ５—２４５０Ｍ ｎ０ ６４化｝】（：ｆＪ１． ＩＩ）Ｅ为、’ｉｓｕａｌ‘ｌｕｄｉｔ｝２（）ｌ ５，（ ＩｌＰｌｌ ’３．０ 边缘滤波、 －５ｉ＇￣９，１ｔｑ像如 ７所爪． ８ ｌ ７ 边｝ｆｇ 攻』， 的 Ｉ冬ｌｆ ｉ冬Ｉ ８ ｌ，ａｌｌｌⅢ Ｉ ：ｒ边 挚｝ １ ＝ ．ｃ熊塞 娑 阈■Ｉ一 磐ＩＩ照：０ 【 鲤蛹一 参蠹值：’ 参戢谴：囊—■ ．９０． ２０１ ７ 示了Ｉ ａｐ］ａｃ・ｅ钎：于边缘提取 像， ９表玎 Ｓｒ ｌ ｒ Ｊ 子边缘愉测 像 比较 ９（ａ）、罔９（Ｉ１）、网９（ｆ、）参数为０的时候 特点，参数为Ｉ时边缘比较清晰．但是细节处存在部 分粘连，都不适合进行重建 箅子进行边缘提取 １０表示使用Ｃａｎｎｙ ＳｏＭｅ算子细 较 清楚，但是存在冈像边缘不清等 参数值：：７０……～ 使川Ｃａｎｔｏ算子对 像进行边缘检测对于细节 ｒ要’ 霸嚼 函圜■■■■■■ 参数值：：４４ 图１０（ａ）ｃａｎｎｙ算子‘街罔像 罔１０（ｈ）Ｃａｎｎｙ算子保留颌骨 Ｌｉｎ Ｈｕｉ，Ｚｈａｏ Ｃｈａｎｇｓｈｅｎｇ，Ｓｈｕ Ｎｉｎｇ．Ｄｅｔｅｃｔｉｏｎ ａｎｄ Ｅｖａｌｕ— 的体现比较其他方法更高．参数的可调性也更优秀 图１０（ａ）显示在参数为７０时罔像可以只保留牙齿 图像，图１０（ｂ）则可以在参数４４时保留颌骨部分的 ａｔｉｏｎ ｏｆ Ｃａｎｎｙ Ｏｐｅｒａｔｏｒ Ｂａｓｅｄ ｏｎ［Ｊ１．Ｊｏｕｍａｌ ｏｆ Ｈｅｉｈｍｇｊｉａｎｇ Ａｃａｄｅｍｙ ｏｆ Ｅｎｇｉｎｅｅｒｉｎｇ，２００３，（０２）：６－１ ６． 图像＝Ｃａｎｎｙ算子可以在实际应用中根据需要调节 参数达到控制边缘部分的功能。 ４】余洪山，王耀南．一种改进型Ｃａｎｎｙ边缘检测算法ｌＪ１．计 算机工程与应川，２００４．２０：２７—２９． Ｙｕ Ｈｏｎｇｓｈａｎ．Ｗａｎｇ Ｙａｏｎａｒｔ．Ａｎ Ｉｍｐｒｏｖｅｄ Ｃ￣ｍｎｙ Ｅｄｇｅ Ｄｅ— ６结束语 ｔｅｔ‘ｌｉｏｎ Ａｌｇｏｒｉｌｈｍ…．Ｃｏｍｐｕｔｅｒ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ａｐｐｌｉｃａ— ｔｉｏｎｓ．２００４．（２０）：２７—２９． 对数字图像处理及其ＯｐｅｎＣＶ实现进行了学 习，在了解ＣＴ成像技术的原理进行深入了解之后 ５１曹璐．戴青云．潘晴．一种改进型Ｃａｎｎｙ联合阈值分割的 图像背景去除在外观专利闭像检索中的应用ｆＪ１．计算机 针对热门的口腔ＣＴ罔像进行系统的整理优化，将 ＣＴ机扫捕而得的图像进行Ｃａｎｎｙ算子边缘提取等 图像预处理。在算法实现上提 了方便调试的轨迹 条，可以很好地对各种参数进行可视的实时调节．便 于找出最优方案：为口齿排列网像的重建工作做ｒ 较为完善的准备Ｔ作一 参考文献： ｒ程与应．【｝ｊ，２０１２，（０１）：２０８—２１ １ Ｃａｏ Ｌｕ．Ｄａｉ Ｑ．１１ｇｙｔｌＩ１．Ｐａｎ Ｑｉｎｇ． Ｆｈｅ Ｂａ＜。ｋｇｒｏｕｎ＜Ｉ Ｉｍａｇｅ ｏｆ ａ Ｍｏｄｉｉｆｅｄ Ｃａｎｎｙ Ｃｏｍｂｉｎｅｄ ｗｉｔｈ Ｔｈｒｅｓｈｏｌｄ Ｓｅｇｍｅｎｉｆａｉｆｉｏｎ ｉｎ ｔｈｅ Ｒｅｍｏｖａｌ ｏｆ Ｄｅｓｉｇｎ Ｐａｉｆｅｎｌ ｈｎａｇｅ Ｒｅｔｒｉｅ￣ｒａｌ Ｕｓｉｎｇ【Ｊ１． Ｃｏｍｐｕｔｅｒ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ａｐｐｌｉｃａｔｉｏｎｓ，２０ ｌ ２，（Ｏ ｌ １：２０８—２ ｌ １ ６１刘荣科，张晓林，廖小涛，聂振钢．基于可逆提升小波变换 平１ｌ上下文预测的ＳＡＲ冈像无损压缩算法ＩＪ１．航空学报． ２００２．ｆ０６）：５３４－５３７． Ｌｉｕ Ｒｏｎｇｋｅ．Ｚｈａｎｇ Ｘｉａｏｌｉｎ，Ｌｉａｏ Ｘｉａｏｔａｏ．Ｎｉｅ Ｚｈｅｎｇａｎｇ． ｆｌ１姜璜．基于 绘制的 像三维可视化算法研究【Ｄ１．成都： 电子科技大学．２００７． Ｊｉａｎｇ ＹＩＩＩ１．Ｓｕｒｆａ＜・Ｐ Ｒｅｎｄｅｒｉｎｇ／ｈｅ Ｉｍａｇｅ Ｂａｓｅｄ ｏｎ Ｌｏｓｓｌｅｓｓ Ｃｏｍｐｒｅｓｓｉｏｎ Ａｌｇｏｒｉｔｈｍ ｏｆ ＳＡＲ ｌｎｍｇｅ Ｂａｓｅｄ ｏｆｆ Ｗａｖｅｌｅｔ Ｔｒａｎｓｉｆｔｒｍ ｆｉｎ（１ Ｒｅｖｅｒｓｉｈｌｅ Ｃｏｎｔｅｘｌ Ｐｒｅｄｉｔ、ｔｉｏｎ…． Ａｉｒｅｒａｆｌ ２００２．（０６）：５３４—５３７． Ｔｈ，‘ＰＰ—ｄｉｍｅｎｓｉ．ｎａ１、’ｉｓｕａｌｉｚａｌｉｏｎ ＡＩｇｏｒｉｆｉｈｍ ＩＤＩ．Ｃｈｅｎｇｄｕ： 【Ｊ 、’Ｐｒｓｊｆｖ ｏｆ Ｅｌｅｃｔｒｏｎｉ（’Ｓｔ‘ｉｅｍ ｒ ａｎｔｉ ｒｅ（、ｈｎｏＩｏｇｙ．２００７ ７ｌ袁勇．一种改进的中值滤波方法 成都理１ 大学学报 【２ｌ 群．ＣＴ断层 像三维重建巾捕值技术研究 应Ｊ｛ｌＩＤ］． 苏州大学．２００７． ￣＇ａｎｇ Ｑｎｎ．Ｓｔｕｄｙ ａｎ（Ｉ ＡＩ）ｐｌｉ， ａｉｆｉｏｎ ｎｆ Ｉｎｌｅｔ ｐｎｌａｔｉｏｎ ｒＰＰｈ— ｎｉｑｕｅ ｉｎ ３Ｄ Ｒｅｅｏｎｓｆｉｆｌｌｅｉｆｉｏｎ ｏｆ ＣＴ Ｆａｕｈ Ｉｍａｇｅｓ【Ｄ１．Ｓｕｚｈｏｕ： Ｓｕｚｈｏｕ［ｉｎｉｖｅｒｓｉｆｉｙ，２００７ （ｆｌ然科学版），２０１３＇Ｙ．４０；Ｎｏ．１７７０２：１２５—１２９． Ｙｕａｎ Ｙｏｎｇ．Ａｎ ｉｍｐｒｏ￣ｒｅｄ Ｍｅｄｉａｎ Ｆｉｈｅｒｉｎｇ Ｍｅｔｈｏｄ【ＪＪ．Ｊｏｕｒ— ｎａｌ ｏｆ Ｃｈｅｎｇｄｕ Ｕｎｉｖｅｒｓｉｆｉｖ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ（Ｎ ＡＴＩ ＪＲＡＬ ＳＣＩ— ＥＮ（：Ｅ ＥＤＩ一１ＩＯＮ）．２０１３．ｖ．４０；Ｎｏ．Ｉ７７０２：１２５－１２９． ８】陈学奇．罔像测量中的算法研究【Ｄ】．长沙：中南大学． ２０ｌ３． 林卉．赵长胜．舒宁．基于Ｃａｎｎｙ算子的边缘检测及评价 …．黑龙ｆ１ Ｉ 程学院学报．２００３．（０２）：６—１ ６． Ｃｈｅｒｔ Ｘｕｅｑｉ．Ａｌｇｏｒｉｌｈｍ Ｒｅｓｅａｒｅｈ ｉｎ Ｉｍａｇｅ Ｍｅａｓｕｒｅｍｅｎｉｆ【Ｄ】 ．ＣｈＨｎ ｓｈＨ・《 ｅｎｌｆｒ８ｌ Ｓｎｌｌｌｈ【７ｎｉｖｅｒｓｉ！ｖ ２０１ ３
