基于CUDA的高速并行高斯滤波算法.docVIP

基于 ＣＵＤＡ 的高速并行高斯滤波算法卢文龙 王建军 刘晓军 基于 ＣＵＤＡ 的高速并行高斯滤波算法 卢文龙 王建军 刘晓军 （华 中 科技大学机械科学与工程学院 ，湖 北 武 汉 ４３００７４） 摘要 为加快表面三维形貌分析中高斯滤波算法的执行速度，提出了一种基于计算统一设备构架（ＣＵＤＡ）的 高斯滤波算法来实现高速并行处理．分析高斯滤波算法原理和 ＣＵＤＡ 并 行 计 算 体 系，将 ＣＵＤＡ 并 行 计 算 技 术引入到表面分析领域．针对高斯滤波 数 据间依赖性弱和 ＣＵＤＡ 采用单指令多线程 （ＳＩＭＴ）执 行 模 型 的 特 点，总结出适合于 ＣＵＤＡ 的并行高斯滤波算法流程．实验证 明：该 方 法 与 ＣＰＵ 串行处理方法相比，其 加 速 比 达到４０ 倍以上，可以有效提高数据处理能力． 关键词 高斯滤波；计算统一设备构架；图像处理器；表面形貌；并行处理 中图分类号 文献标志码 文章编号 １６７１－４５１２（２０１１）０５－００１０－０４ ＴＨ１３２．１ Ａ ＣＵＤＡ－ｂａｓｅｄａｌｇｏｒｉｔｈｍｆｏｒｈｉｇｈ－ｓｐｅｅｄｐａｒａｌｌｅｌＧａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇ ＬｕＷｅｎｌｏｎｇ ＷａｎｇＪｉａｎｊｕｎ ＬｉｕＸｉａｏｊｕｎ （ＳｃｈｏｏｌｏｆＭｅｃｈａｎｉｃａｌＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇ，ＨｕａｚｈｏｎｇＵｎｉｖｅｒｓｉｔｙ ｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，Ｗｕｈａｎ４３００７４，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ ＩｎｏｒｄｅｒｔｏｓｐｅｅｄｕｐｔｈｅＧａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇｉｎｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｓｕｒｆａｃｅｔｅｘｔｕｒｅａｎａｌｙｓｉｓ，ａｎ ｅｆｆｉｃｉｅｎｔｍｅｔｈｏｄｂａｓｅｄｏｎｃｏｍｐｕｔｅｕｎｉｆｉｅｄｄｅｖｉｃｅａｒｃｈｉｔｅｃｔｕｒｅ （ＣＵＤＡ）ｉｍｐｌｅｍｅｎｔｅｄｏｎｇｒａｐｈｉｃｐｒｏ－ ｃｅｓｓｉｎｇｕｎｉｔ（ＧＰＵ）ｗａｓｄｅｓｉｇｎｅｄ．ＣＵＤＡｐａｒａｌｌｅｌｃｏｍｐｕｔｉｎｇｔｅｃｈｎｏｌｏｇｙｗａｓｉｎｔｒｏｄｕｃｅｄｉｎｔｏｔｈｅｓｕｒ－ ｆａｃｅｔｅｘｔｕｒｅａｎａｌｙｓｉｓａｒｅａｂｙａｎａｌｙｚｉｎｇｐｒｉｎｃｉｐｌｅｏｆＧａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇａｎｄＣＵＤＡｃｏｍｐｕｔｉｎｇａｒｃｈｉｔｅｃ－ ｔｕｒｅ．ＰａｒａｌｌｅｌＧａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇａｌｇｏｒｉｔｈｍｕｓｅｄｉｎＣＵＤＡ ｗａｓｇｉｖｅｎｆｏｒｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｌｏｗｄｅ－ ｐｅｎｄｅｎｃｙｏｆ Ｇａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇａｎｄ ＳＩＭＴ （ｓｉｎｇｌｅｉｎｓｔｒｕｃｔｉｏｎ ｍｕｌｔｉｐｌｅｔｈｒｅａｄ）ｅｘｅｃｕｔｉｏｎ ｍｏｄｅｌｏｆ ＣＵＤＡ．ＥｘｐｅｒｉｍｅｎｔｓｐｒｏｖｅｔｈａｔｔｈｅｃａｌｃｕｌａｔｉｏｎｓｐｅｅｄｂａｓｅｄｏｎＣＵＤＡｉｓ４０ｔｉｍｅｓｆａｓｔｅｒｔｈａｎｔｈａｔｏｆ ｔｒａｄｉｔｉｏｎａｌｓｅｑｕｅｎｃｅａｌｇｏｒｉｔｈｍｂａｓｅｄｏｎＣＰＵ，ｃａｎｅｆｆｅｃｔｉｖｅｌｙｉｍｐｒｏｖｅｄａｔａｐｒｏｃｅｓｓｉｎｇｃａｐａｂｉｌｉｔｙ． Ｇａｕｓｓｉａｎｆｉｌｔｅｒｉｎｇ；ｃｏｍｐｕｔｅｕｎｉｆｉｅｄｄｅｖｉｃｅａｒｃｈｉｔｅｃｔｕｒｅ；ｇｒａｐｈｉｃｐｒｏｃｅｓｓｉｎｇｕｎｉｔ；ｓｕｒｆａｃｅ ｔｅｘｔｕｒｅ；ｐａｒａｌｌｅｌｐｒｏｃｅｓｓｉｎｇ Ｋｅｙｗｏｒｄｓ 高斯滤波作为一种零相移的滤波方法，广泛 应用于三维表面形貌分析中，并被国际标准化组 织ＩＳＯ／ＴＣ２１３确认为国际标准．传统高斯滤波算 法通过 ＣＰＵ 来 实 现，在 时 域内测量数据与高斯 权函数进行卷积或两者在频域内相乘再反变换到 时域从而得到评定基准［１－２］．当数据量较大时这种 方法效率非常低下，难以满足实际需求．近年来基 于图像处理 器 （ＧＰＵ）的 大 规 模 通 用 并 行 计 算 为 三维表面形貌的高效分析计算提供了可能，且其 成本和功耗只有传统基于 ＣＰＵ 计 算 的 １／１０ 和 １／２０［３］，例如 ＮＶＩＤＩＡ 公司提出的计算统一设备 构架（ＣＵＤＡ），它 是 一 种 将 ＧＰＵ 作 为 数 据 并 行 计算设备的 软 硬 件 体 系，其 特 点 是 将 ＣＰＵ 作 为 主机，ＧＰＵ 作为协处理器［４］．在这个模型中 ＣＰＵ 和 ＧＰＵ 协同工作，ＣＰＵ 负责 逻辑性强的事务和 串行计算任 务，而 ＧＰＵ 主 要 处 理 高 度 线 程 化 的 并行计算任务． 本文以 ＧＰＵ 为计 算平台提出了 一 种基于 ＣＵＤＡ 的 高 速 高 斯 滤 波 算 法．通过实验将基于 ＧＰＵ 并行处理方法与基于 ＣＰＵ 的串行方法进行 收稿日期 ２０１０－１０－１４． 对比，