-
摘要: 肺部CT图像病变区域检测是肺病辅助诊断技术的重要研究内容,其通过自动分析CT图像并输出病变区域的位置和尺寸等信息,帮助放射科医生做出决策,有利于肺病的早期发现与治疗.本文回顾了肺部CT图像中病变区域自动检测方法所取得的进步,并引入一个通用框架表示和描述现有方法,对2012年以来肺部病变区域辅助检测算法进行了系统性分析和性能汇总.最后讨论了目前存在的问题和有待克服的困难,探讨了未来可能的发展方向.Abstract: Automatic detection of lesion regions in lung CT images is an important research topic in computer aided diagnosis of lung diseases. The system can automatically analyze CT images, output the locations and sizes of lesion regions to help radiologists make decisions, and promote early detection and therapy of lung diseases. In this paper we review the achieved progress of automatic detection methods of lesion regions in lung CT image, and introduce a generic structure for expressing and describing existing detection methods. Furthermore, we provide a systematic analysis and comprehensive performance summary of the latest detection algorithms from 2012. Finally, we point out the challenges ahead, and discuss the future direction of computer aided detection of lung lesions.
-
Key words:
- Lung CT /
- lung nodule /
- lung vessel /
- lymph node /
- computer aided detection
-
在控制与决策等领域中, 往往会遇到许多优化问题, 有时我们不仅需要找到全局最优解, 还需要找到局部最优解进行辅助决策, 这类优化问题称为多峰搜索 (Multi-peaksearching) 或多峰优化 (Multimodal optimization) 问题.随着问题的复杂性不断提升, 传统的数学规划理论很难解决.
大规模可并行的群智能算法的出现为解决多峰优化问题带来了新的思路.其中粒子群算法 (Particleswarm optimization, PSO) 是一种基于群体协作的随机搜索算法, 但由于其对局部最优解搜索能力不好并且需要较多的先验知识, 许多学者对其进行了改进[1-4].同样, 差分进化算法 (Differential evolution, DE) 作为进化算法的一个重要分支, 收敛速度比PSO快, 但由于贪婪选择策略, 在求解具有大量局部优点的问题时, 种群极易陷入局部最优.有学者使用临近衰减策略[5]和双目标[6]等策略来解决这个问题.这些算法能够在局部找到接近最优的极值点, 但无法给出所有的极值点.也有学者尝试使用蚁群算法[7]和在极值点附近循环扰动的Nelder-Mead单纯形法[8]来找出所有极值, 但采用的测试函数极值点均较少, 而且搜索时间较长, 不大适合用于搜索有极多极值的函数.
量子算法以其完备的理论和强大的并行能力[9], 被很多学者用于群智能算法中.文献[10]在2010年提出了模拟量子谐振子的概念, 文献[11]中给出了算法实现, 即多尺度量子谐振子算法 (Multi-scale quantum harmonic oscillatoralgorithm, MQHOA).该算法数学原理简单物理含义明确且无需先验知识.完备的采样能力和精确的搜索聚焦能力使MQHOA能较为准确地在多维空间实现对最优中把MQHOA与模拟退火算法 (Simula-ted annealing, SA)、量子解的搜索.文献[12]将退火算法 (Quantumannealing, QA)、量子粒子群算法 (Quantum-behaved particle swarmoptimization, QPSO) 进行对比, 对15种不同的标准测试函数求最小值, 实验结果进一步说明了其具有更好的计算速度和计算精度.
MQHOA集中处理单峰优化问题, 主要根据量子隧道效应[12]保证全局采样点数目足够大以避免陷入局部较优解位置.算法收敛的依据是这些采样点都收敛到一个最优解上, 但对于多峰函数, 采样点将分散在各个峰上而得不到聚集, 导致算法无法收敛.为扩展MQHOA对多峰优化问题的应用, 本文提出了基于划分的多尺度量子谐振子算法 (Partition-based MQHOA, P-MQHOA), 它把求解域均匀划分成N个局部域, 把MQHOA的收敛特性应用在局部域内, 每个局部域中心内放置一个采样点进行谐振子迭代操作, 那么这些采样中心点的平衡位置即为所求的解.实验使用可调节参数的双势阱量子系统函数来验证算法对于等值极值点和不等极值点寻优准确性.对于全部极值点的搜寻效果, 本算法对比了蚁群算法对同一个函数的全极值搜寻, 实验结果超出了蚁群算法.此外, 通过对14个基准函数的全局多峰优化, 对比遗传算法、粒子群算法和差分进化算法效果都要好.当极值点极多的情况下也可以采用多级分辨率策略求解, 求解精度高, 查全率高, 速度快.
1. 多峰优化问题描述
工程上的多峰优化是把这个工程问题映射到一个多维函数上, 要得到多个最优方案就是在这个多维函数上求多个极大值 (或者是求这个多维函数取相反数的多个极小值).
为了全文描述方便, 本文给出以下四个定义:
定义 1.对于某 $D$ 维N极值函数 ( $N\ge2$ ), 第 $i$ 个极值点和第 $j$ 个极值点相邻, 他们坐标之间第 $k$ 维差的绝对值记作 $\lambda_{i, j, k}$ .取 $\lambda_{i, j}=\min\nolimits_{1\le k\le D}\lambda_{i, j, k}, \lambda_{\rm min}=\min\nolimits_{1\le i < j\le N}\lambda_{i, j}$ , 把 $\lambda_{\rm min}$ 称为该函数实际波长分辨率.
定义 2.工程上需要识别某 $D$ 维多极值函数上的极值点满足如下要求: $\lambda_{i, j}\ge \lambda_0\ge \lambda_{\rm min}$ .把 $\lambda_0$ 称为需求波长分辨率, 这些极值点称之为满足需求波长分辨率 $\lambda_0$ 的极值点.
定义 3.把算法求得的极值点坐标与对应的实际极值点坐标之间允许的各维最大差的绝对值称作该算法求解的位置精度, 记为 $A_l$ .
定义 4.某 $D$ 维函数定义域上界为 ${\pmb u}$ , 下界为 ${\pmb s}$ , 那么称 ${\pmb l}={\pmb u}-{\pmb s}=(l_1, l_2, \cdots, l_D)$ 为该函数的长度域, 其中 $l_k$ 为该函数在第 $k$ 维上定义域的长度.
多峰优化问题的数学描述即寻找 $D$ 维函数 $f (x_1, x_2, \cdots, x_D)$ 在定义域上某点 $X_i (x_{i, 1}, x_{i, 2}, \cdots, x_{i, D})$ 处取得极值. $X_i^-=(x_{i, 1}-dx_1, x_{i, 2}-dx_2, \cdots, x_{i, D}-dx_D), X_i^+=(x_{i, 1}+dx_1, x_{i, 2}+dx_2, \cdots, x_{i, D}+dx_D)$ , $f_i^-$ 为函数在点 $X_i^-$ 求得的值, $f_i^+$ 为函数在点 $X_i^+$ 求得的值, $f_i$ 为函数在点 $X_i$ 取得的函数值, $X_i^*$ 为实际极值点, 那么判断该函数在 $X_i$ 处取得极小值的依据:
$\begin{align}\left\{\begin{array}{ll} f_i^->f_i, f_i^+>f_i \\ X_i\in\{X\mid\left|X-X_i^*\right|\le A_l \} \\ \end{array} \right.\end{align}$
(1) 判断该函数在 $X_0$ 处取得极大值的依据:
$\begin{align}\left\{ \begin{array}{ll} f_i^-<f_i, f_i^+<f_i \\ X_i\in\{X\mid\left|X-X_i^*\right|\le A_l \} \\ \end{array} \right.\end{align}$
(2) 2. 基于划分的MQHOA原理
量子谐振子是物理上描述微粒动态振动的物理量, 可以近似描述分子的运动规律.MQHOA算法把量子谐振子群作为智能粒子群来解决优化问题, 振子的平衡位置即势能最低位置, 就是所求函数的最优解.量子谐振子从高能态到低能态是一个逐渐收敛的过程[13], 直到收敛到基能态: $\psi_0(x)=\frac{\sqrt{\alpha}}{\sqrt[4]{\pi}}{\rm e}^{-\alpha^2\frac{(x-x_i)^2}{2}}$ , 定义MQHOA算法在尺度 $\sigma_s$ 下的波函数为 $k$ 个以 $x_i$ 为中心的高斯概率密度函数的迭加:
$\begin{align}\psi_{QHO}(x)=\sum\limits_{i =1}^k{\frac{1}{\sqrt{2\pi}\sigma_s}}{\rm e}^{-\frac{(x-x_i)^2}{2\sigma_s^2}}\end{align}$
(3) $\psi_{QHO}(x)$ 近似代表了以 $k$ 个最优解为中心的目标函数可行域上的概率分布, 通过引入多尺度优化函数二进信息高斯采样模型, 二进小波尺度函数采用高斯函数按照不同尺度和不同精度要求进行聚焦搜索[9], 可以保证信息的不遗漏采样, 从而获得全局最优解, 详细算法描述参见文献[11].
P-MQHOA扩展了MQHOA对多峰优化问题的处理能力, 并且具有对所有最优解的记忆能力和更快的收敛速度.其主要思想是把MQHOA的全局优化特性局限在了一个较小划分区域内, 在这个区域内只设置一个高斯采样中心点, 并且让这个采样中心点在反复迭代过程中不断收敛到精度范围内.P-MQHOA对整个函数定义域划分的N个区域与N个采样中心点一一对应, 最终找到N个解 (包括重复解).由于这N个采样中心点之间的采样计算是独立进行的, 并不会由于某几个采样中心点收敛不了而影响其余采样中心点, 因此其收敛速度要比MQHOA要快得多, 并且算法通过对记忆的N个解进行整理可以得到全局最优解.
算法对函数定义域进行均匀划分, 划分区域内如何选择要生成的高斯采样点的标准差直接影响算法的性能和执行结果.各维度划分成的各均匀区段中心位置处随机生成 $M$ 个标准差为 $\sigma_s$ 的高斯分布点, 在同一个区段内, $\sigma_s$ 与这些分布点有如下关系: $\sigma_s$ 越大分布点越接近均匀分布, $\sigma_s$ 越小分布点越集中在中心点, 通过 $\sigma_s$ 的衰减可以使得这些分布点最终聚集在搜索精度范围内.算法初始的 $\sigma_s$ 大小很重要, 分布点太均匀的话很大概率会搜索到临近的其他区段, 分布点太集中的话很大概率会忽略区域边界附近的极值点.假如某个区段内初次生成的高斯分布点超出边界的概率为 $p$ , $M$ 次伯努利实验后可得到超出边界的分布点数目平均值为 $Mp$ .算法一般选择 $\sigma_s$ 为该区段长度的四分之一, 通过查找正态分布表可知, $p=0.0456$ , 即每个分布点超出区段的概率为0.0456.这样设置不仅可以保证对本区段内的极值搜索, 又不失高斯分布的中心点集中特性.
对定义域下界为 ${\pmb s}$ 长度域为 ${\pmb l}$ 的 $D$ 维函数在第 $j$ 维上进行均匀划分, 划分成 $k$ 个区段后高斯概率密度叠加可得第 $j$ 维上采样点的概率分布:
$\begin{align}\psi_{jQHO}(x)=\sum\limits_{i =1}^k{\frac{1}{\sqrt{2\pi}\sigma_s}}{\rm e}^{-\frac{\left(x-\frac{l_i(i-\frac{1}{2})}{k}-s_i\right)^2}{2\sigma_s^2}}\end{align}$
(4) 其对应的密度图如图 1所示, 虚线为原来的高斯概率分布, 实线为叠加后的概率分布.
叠加后某一区段的高斯分布由于临近区段分布点的 $p$ 概率跳入而弥补了本区段 $p$ 概率跳出流失的分布点, 这样就保证了 $M$ 个采样点在各区段数目相等, 并且有 $Mp$ 个点作为区段之间的交流媒介保证了各区段的联系.这种联系使得在某区段内找不到极值点的情况下可以跳出该区域协助临近区段寻找极值点.把各个维度上的均匀划分区段组成一个个均匀划分区域, 在各个划分区域内寻找极值是相互独立的操作, 故算法可以多线程并行执行.
为了保证在需求波长分辨率 $\lambda_0$ 下识别该分辨率下的所有极值点, 本算法需要在可行域上划分 $K$ 个区域并设置 $K$ 个采样中心点, $K$ 的大小决定算法能准确寻找极值的数目.以寻找一维函数极小值为例, 定义搜索区域长度为函数定义域上均匀划分的区段长度.如图 2所示, $\psi (x)$ 为初始采样点在区域内分布的概率密度函数.当搜索区域长度超出该一维函数的波长分辨率 $\lambda_0$ (两个波谷之间的水平距离) 时, 算法在当前搜索区域内至少能找到一个极值点, 如图中函数 $f_1 (x)$ 所示; 当搜索区域长度等于 $\lambda_0$ 时, 极值点或者在该区域被找到或者在邻近区域被找到, 如图中函数 $f_2 (x)$ 所示; 当搜索区域长度小于 $\lambda_0$ 时, 当前搜索区域无法找到极值点, 只能在邻近区域找到, 如图中函数 $f_3(x)$ 所示.因此, 本算法需要通过设置合理的采样中心点数使得搜索区域长度大于等于 $\lambda_0$ .
$K$ 的大小也影响算法的执行速度.以某 $D$ 维函数为例, 其可行域 ${\pmb l}=(l_1, l_2, \cdots, l_D)$ , 如果需要找到所有极值, 那么搜索区域长度正好为该函数的实际波长分辨率 $\lambda_{\rm min}$ , 以此算出合理的采样中心点数.在第 $i$ 维上置采样中心点个数 $k_i$ 为 $\lceil{ \frac{l_i}{\lambda_{\rm min}}}\rceil $ , 即在该维度上划分出 $k_i$ 个区段, 要找到 $D$ 维函数所有极值点则需初始划分区域数 $K=\prod\nolimits_{i=1}^Dk_i$ .当 $\lambda_{\rm min}$ 很小而 $D$ 很大时, 算法需要花费极多CPU时间在所有区域内进行计算.类似多分辨率的小波分析[14], 本文提出多级分辨率策略, 其基本思想是先在大波长分辨率下进行区域划分, 在找到极值点的区域内再进行小波长分辨率的划分, 而在不存在极值点的区域就不再进行划分计算, 如此重复可以使搜索区域长度逼近 $\lambda_{\rm min}$ , 从而找到满足需求分辨率极值点.这种策略避免了原来固定分辨率下对无极值区域的计算浪费, 而把更多的CPU计算集中在极值较多的区域.实际应用中需求波长分辨率 $\lambda_0$ 一般比 $\lambda_{\rm min}$ 大得多, $K$ 取 $\prod\nolimits_{i=1}^D\left\lceil\frac{\lambda_0}{l_i}\right\rceil$ .当 $K$ 不是很大时, 可以直接采用固定分辨率策略, 反之采用多级分辨率策略.
3. P-MQHOA算法过程描述
3.1 P-MQHOA算法的三个阶段
P-MQHOA主要有三个阶段.初始划分阶段:该阶段主要是对 $D$ 维度函数定义域进行区域均匀划分, 然后在区域中心放置一个采样中心点; 量子谐振子收敛阶段:该阶段采用高斯函数来进行采样最终定位中心点; 归并阶段:该阶段把位置精度范围内的中心点归为一类, 并从该类中选出最好的中心点.
3.1.1 初始划分阶段
输入需求波长分辨率 $\lambda_0$ 下所求 $D$ 维度函数定义域上界为 ${\pmb u}$ , 下界为 ${\pmb s}$ , 由定义4可求得长度域 ${\pmb l}$ .对于第 $j$ 维划分区段数记作 $k_j$ , 若算法采取固定分辨率策略, 初始第 $j$ 维上划分区段数 $k_j$ 为 $\lceil\frac{l_i}{\lambda_0}\rceil$ ; 若算法采取N级分辨率策略, 初始第 $j$ 维上划分区段数为 $\lceil\sqrt[N]{\frac{l_j}{\lambda_0}}\rceil$ .总采样点数 $K=\prod\limits_{j=1}^Dk_i$ , 初始化采样中心点, 某个中心点在第 $i$ 个区段第 $j$ 维的坐标值可以记作 $c_{i, j}=\frac{l_j (i-1/2)}{k_i}i+s_j, i\in[1, k_j], j\in[1, D]$ , $i$ 与 $j$ 均为整数.
3.1.2 量子谐振子收敛阶段
基于划分的MQHOA的采样中心点之间互不影响, 在每个划分区域内的收敛过程都是独立的, 因此各中心点可以并行执行.如图 3所示为寻找某一维函数 $f (x)$ 极大值的收敛过程.初始过程的采样点概率分布为6个高斯分布函数的叠加, 见式 (4).其区域长度正好为这个叠加函数的周期长度.按照前面说明, 该区域的高斯函数提供 $M (1-p)$ 个采样点, 由于在可行域外产生的采样点很少可以忽略不计, 那么很容易证明该区域内约有 $Mp$ 个采样点由其他区域的高斯采样提供, 而越靠近的区域提供越多.在图 3第一步中左起第二个区域内约有 $Mp/2$ 个采样点落在第三个区域, 而第三个区域内的采样点效果显然比原区域好, 因此下一轮采样将在第三个区域进行.在找到当前区域最大值采样点后, 以这个点为中心, 搜索区域长度衰减, 重新在新的搜索区域寻找最大值, 如图 3中第二第三步, 以此类推, 直到所有中心点要么收敛在极值点附近, 要么找不到极值点.图 3最后一步中所有中心点都收敛在极值点附近.
3.1.3 归并阶段
图 3最后收敛阶段有6个中心点分别收敛于4个波峰, 其中第2、3、4个中心点聚集在一个极值点附近, 需要对这些点进行归并操作.若任两个点各维坐标最大差的绝对值小于位置精度 $A_l$ , 按照定义3, 把这两个点归为一类, 最终从此类中取出对应函数值最大 (最小) 的点作为最终的极值点.
3.2 P-MQHOA算法流程
综上所述, 整个算法流程可以简要描述如下:
步骤 1.初始化:
步骤 1.1.输入 $\lambda_0, M, A_l, {\pmb u}, {\pmb s}$ .
步骤 1.2.根据公式初始化长度域 ${\pmb l}$ , 标准差 $\sigma_s=\frac{\lambda_0}{4}$ , 在各维度上均匀划分产生 $K$ 个中心点, 其中第 $i$ 个点 $center_i=(c_{i, 1}, c_{i, 2}, \cdots, c_{i, D})$ .
步骤 2.收敛过程:
步骤 2.1.在各维上以 $c_i$ 为中心点, $\sigma_s$ 为标准差进行 $M$ 次高斯采样, 当采样点超出可行域边界时以边界点为准.
步骤 2.2.从 $M$ 个采样点中选择对应函数值最大 (最小) 的一个点作为新的中心点, $\sigma_s$ 衰减为原来一半.如果满足 $\sigma_s\geq A_l$ , 重复执行步骤2.1, 否则继续执行步骤3.1.
步骤 3.归并:
步骤 3.1.对最后生成的 $K$ 个解根据式 (1) 或式 (2) 作极值判定, 去除不满足式子的点.
步骤 3.2.计算各相邻点之间各维度差值的最小值, 若小于 $A_l$ 则合并这些点并取其中最大 (最小) 的作为最终解.
3.3 P-MQHOA算法两种策略
针对维度较低、需求波长分辨率不高的函数进行多峰寻优, 采用固定分辨率策略即可解决.但维度较高或者需求波长分辨率很小情况下, 应当采用多级分辨率策略.
对于N级分辨率, 算法经历初始划分阶段和量子谐振子收敛阶段后, 找到输出中心点不为空所对应的区域, 对这些区域执行下一级的初始划分阶段的操作并对划分的子区域执行量子谐振子收敛阶段的操作, 递归搜寻到第N级后对所有输出的中心点执行归并阶段的操作, 可输出最终的极值点.
以图 4为例, 对 $f (x, y)$ 在需求分辨率为0.1下求极值, 长度域 $l_x=l_y=5$ , 采用三级分辨率策略初始划分区域数为 $\lceil\sqrt[3]{\frac{5}{0.1}}\times\sqrt[3]{\frac{5}{0.1}}\rceil=16$ , 如子图 $4 (a)$ 所示; 经过收敛和归并过程后确认极值点在6号与11号区域, 重新在这两个区域进行第二级划分, 如子图 4(b)所示; 重复以上过程第三级划分后实际分辨精度可达到 $\frac{5}{4^{3}}=0.078$ , 满足需求分辨率终止迭代, 如子图 4(c)所示.若每个中心点的采样次数为 $M$ , 可以计算总的采样次数为 $16(1+2+3) M=96M$ , 而采取固定分辨率需要次数为 $50\times50M=2 500M$ , 可见其省去很多计算时间, 但多级数情况下也容易丢失部分极值点.因为 $M$ 在粗分辨率下相对于细分辨率就要小, 就好比把黄豆随机撒在球场上和撒在球门里的区别.如子图 4(b)所示若11号区域未找到极值点那么在子图 4(c)中就会丢失很多极值点.因此 $M$ 的选择和级数的选择应该根据实际情况而定.
此外, 多级分辨率策略下不仅在同级区域内横向可并行执行, 不同级区域也可纵向并行执行, 计算过程不会相互干涉, 在保持原QHOA算法低耗时、快收敛、高精度的基础上该算法有很好的多峰寻优能力.
4. 算法实例分析
本节通过寻优精确性、全极值点寻优和全局多峰优化三个角度来测试该算法的综合性能, 实验使用的计算机采用Mac系统, 主频29 GHz, 在JDK7环境下的浮点计算精度为小数点后16位, 实验寻找的极值点默认为极小值点, 设置位置精度 $A_l$ 为0.000001.
4.1 寻优准确性
准确性是衡量一个算法好坏的重要标准, 为了测试本算法在需求波长分辨率下寻找极值点的准确度, 实验使用双势阱量子系统函数来验证.双势阱量子系统是研究微观物理学的一个重要物理模型, 其函数表达式为
$\begin{align}V_{sys}=V_0\frac{(x^2-a^2)^2}{a^4}+\delta x\end{align}$
(5) 该函数有两个极小值, 其中 $V_0$ 表示函数的值域长度, $2|a|$ 表示两个极小值点之间的水平距离, $2|\delta||a|$ 表示两个极小值点之间的垂直距离.置 $V_0=10$ , $a$ 分别取 $1, 0.5, 0.1, 0.05, \cdots, 0.0000001$ , 对应函数实际分辨率 $\lambda_{\rm min}$ 为 $2, 1, 0.5, 0.1, \cdots, 0.0000005$ .设置函数定义域为 $[-1.5, 1.5]$ , 初始采样中心点数 $K=30$ , 固定分辨率下可识别波长分辨率为 $\lceil\frac{1.5-(-1.5)}{30}\rceil=0.1$ , 实验设置需求波长分辨率 $\lambda_0=0.1$ , 测试算法在实际波长分辨率 $\lambda_{\rm min}$ 取不同值时识别极值的准确度.此外, 为了测试极值点间的垂直距离对寻找极值点准确度的影响, 实验分别取 $\delta=1$ 与0, 前者对应不等值极值双势阱函数, 后者对应等值极值双势阱函数.经过1 000次试验, 统计寻找极值点的次数, 不等值极值点函数统计结果如图 5所示, 等值极值点函数统计结果如图 6所示.其中横坐标为 $a$ , 三条曲线分别代表找到第一个极值点的次数、找到第二个极值点的次数和全部找到极值点次数.
从两张图可以发现本算法可以找出需求分辨率下几乎所有极值.令 $\lambda_0=2|a|$ , 解出 $a_0=0.05$ ; 令 $A_l=2|a|$ , 解出 $a_1=0.0000005. a_0$ 与 $a_1$ 正好是两张图曲线趋势变化的分界点.当 $a < a_0$ 时均能找到两个极值点; 当 $a=a_0$ 时极大概率找到两个极值点; 当 $a_0 < a\leq a_1$ 时, 找到两个极值点的概率趋于稳定; 当 $a>a_1$ 时找到两个极值点的概率逐渐减少.实验结果与前面图 2说明相符合, 证明了需求波长分辨率下寻优准确性.
此外可以发现:
1) 不管对于等极值函数和不等极值函数, 该算法在定义域内均能至少找到一个极值点, 并且随着 $\lambda_{\rm min}$ 减小, 三条曲线变化趋势几乎一致.
2) 图 4中 $a$ 取0.05与0.0005之间时算法更容易找到极值更小的第一个极值点, 也就是说对于不等值极值点, 该算法容易找到极值更明显的点, 这个特性使得该算法适合用于全局优化.
4.2 全极值点寻优
4.2.1 全极值点寻优性能
函数的全部极值点往往可以反映函数整体变化趋势, 大部分算法对全极值点的搜索消耗了很多时间, 而且精度不高, 例如文献[7]中的Ant算法与文献[8]中RePAMO算法全极值点的搜寻成为多峰优化问题上的一个难题, 这里采用两种策略的P-MQHOA作全极值寻优.
实验使用简单一维震荡函数 $y=x{\rm sin}\frac{1}{x}, x\in[-0.5, 0.5]$ , 如图 7.该函数在原点附近有无穷极值, 越靠近原点极值点分布越紧密, 直至波长分辨率趋近于0.故不可能找出理论上的所有极值点, 实验指定需求波长分辨率, 然后找出该波长分辨率下的极值点数.
指定需求分辨率下的理论极值点数目可以这样求.对 $y$ 求导, $y'=-\frac{1}{x}\cos{\frac{1}{x}}+\sin{\frac{1}{x}}$ , 令 $y'=0$ , 有 $\tan{\frac{1}{x}}=\frac{1}{x}$ , 当 $\frac{1}{x}\in[-0.5, 0.5]$ 时可以得到近似解, 其中正近似解 $x_+=\frac{2}{(2k+1)\pi}, k\in N^*$ .相邻两个波谷水平距离 $\lambda=\frac{2}{(2k+1)\pi}-\frac{2}{(2k+5)\pi}$ , 在 $\lambda_0 < \lambda$ 条件下可解得
$\begin{align}k<\sqrt{1+\frac{2}{\lambda_0\pi}}-1.5, k\in N^*\end{align}$
(6) 这是正近似解数, 由于函数对称性, 负近似解数等于正近似解数, 所以满足式 (6) 的最大正整数的两倍记作 $K_0$ , 即为满足需求分辨率的极值点个数.
为验证该算法对多维函数的多极值寻优性能分析, 实验使用Griewank函数, 函数表达形式: $f (x_i)=\sum\limits_{i=1}^D\frac{x_i^2}{4000}-\prod\limits_{i=1}^D\cos\left (\frac{x_i}{\sqrt{i}}\right)+1$ 这里选择维度为2, 如图 8所示, $x_i$ 在[ $-$ 100, 100]内, 该函数有固定等分辨率的极值点数记为 $K_0$ .
实验设置不同需求波长分辨率 $\lambda_0$ , 进行50次独立实验.对比固定分辨率和二级分辨率下算法的性能, 前者找到的极值点个数用 $K_1$ 表示, 后者用 $K_2$ 表示, 平均耗时 (Timeconsuming) 分别为 $TC_1$ 、 $TC_2$ , 单位为秒, 采样中心点数 (Center-point numbers) 分别用 $CN_1$ 、 $CN_2$ 表示.震荡函数结果见表 1.
表 1 一维震荡函数两种策略寻优性能比较Table 1 Comparison of two strategies for one dimensional shock functionK0 λ0 K1 CN1 TC1 K2 CN2 TC2 2 0.1 3.9 10 0.0030 4.0 4 0.0038 12 0.01 13.9 100 0.0180 11.9 10 0.0084 46 0.001 47.7 1 000 0.1499 41.6 31 0.0348 156 0.0001 157.0 10 000 1.312 144.2 100 0.167 500 0.00001 493.2 100 000 10.374 461.0 316 0.767 1592 0.000001 1515.0 1 000 000 95.614 1 419.6 1 000 3.067 $K_0$ 用来表示满足需求波长分辨率 $\lambda_0$ 的极值点数, 实际上算法可以在原点附近找到两个比需求波长分辨率还小的极值点, 如 $\lambda_0=0.1$ 时算法还可以找到分别在 ( $-$ 0.1, 0) 与 (0, 0.1) 域内的两个极值点, 但这两个极值点距离最近极值点的坐标距离远小于 $\lambda_0$ , 因此找到满足 $\lambda_0$ 的极值点数实际为 $K_1-2$ 和 $K_2-2$ .比较 $K_0$ 与 $K_1-2$ 和 $K_2-2 $ 的关系可知:随着 $\lambda_0$ 减小, 算法可以找到更多的极值点, 而且相比多级分辨率策略固定分辨率策略下找到满足 $\lambda_0$ 的极值点数更接近 $K_0$ , 但多级分辨率策略下算法花费时间较短.如取 $\lambda_0$ 为0.000001时固定分辨率与二级分辨率时间消耗相差约三十倍, 其原因是单个采样中心点迭代次数恒定的前提下两种策略算法所用的采样中心点数相差很大, 采用固定分辨率策略共需要 $CN_1=\lceil (\frac{L}{\lambda_0})^D\rceil=\lceil (\frac{(0.5-(-0.5))}{0.000001})^1\rceil=1 000 000$ 个采样中心点, 而采用二级分辨率策略初始只需 $CN_2=\lceil (\frac{L}{\lambda_0})^{D/N}\rceil=\lceil (\frac{(0.5-(-0.5))}{0.000001})^{1/2}\rceil=1 000$ 个采样中心点, 初始平均划分1 000个区域, 这1 000个区域中存在很多无极值的区域, 在第二级划分时不会对这些无极值区域划分, 这样总的采样中心点相对来说就变很少了.同样把函数换成二维Griewank函数进行50次独立实验, 对结果取平均值, 见表 2所示.
表 2 二维Griewank函数两种策略寻优性能比较Table 2 Comparison of two strategies for two dimensional Griewank functionK0 λ0 K1 CN1 TC1 K2 CN2 TC2 1 369 20 75.7 10 0.051 40.3 4 0.036 1 369 10 307.8 400 0.183 231.2 10 0.291 1 369 5 1 049.2 1 600 0.726 617.5 31 1.472 1 369 2 1 333.4 10 000 4.078 670.2 100 1.657 1 369 1 1 340.5 40 000 15.312 783.0 316 9.64 1 369 0.5 1 353.10 160 000 57.272 1 924.1 1 000 21.04 随着 $\lambda_0$ 减小算法找到的极值点变多, 到 $\lambda_0 < 5$ 后两种策略找到的极值点数均开始变得稳定.这是由于Griewank函数极值点坐标均匀分布在各维度上, 各个相邻极值点在各维度上的坐标差均相等, 等于实际波长分辨率, 当需求分辨率小于实际波长分辨率后算法找到的极值点数就趋于稳定了.此外, 采用固定分辨率策略划分的MQHOA几乎能够找到需求波长分辨率下所有极值点, 而二级分辨率策略下查全率较低.如上表所示, 二维Griewank函数在 $\lambda_0$ =2时采用固定分辨率策略几乎能够找到所有极值点, 而采用二级分辨率策略下只能找到实际数目的一半, 这是因为在粗分辨率下忽略的极值点所在区域不会参与到第二级的划分.但整体来说多级分辨率策略可以用较短的时间完成算法搜索, 尤其当需求波长分辨率很小的情况下这种时间差别更明显.
与其他智能算法比较, 本算法在时间和精度上都有很大优势.如文献[5]中利用蚁群算法对函数 $y=x_1^2+x_2^2-\cos (18x_1)-\cos (18x_2)$ , $x_1, x_2\in[-1, 1]$ 求所有极大值点.与利用P-MQHOA算法得出的结果比较如表 3所示.对36个极值点按照各维坐标进行基数排序, No.为找到极值的序号, $((x_1, x_2, f (x_1, x_2))$ 分别为坐标点和对应的值.
表 3 P-MQHOA与蚁群算法全极值搜索比较Table 3 Comparison of total-extremum-searching of P-MQHOA and ant colony algorithmNO. P-MQHOA
(x1, x2, f/(x1, x2))ANT 1 (-0.878093608, -0.878093599, 3.53255484) (-0.8781, -0.8781, 3.5326) 2 (-0.878093550, -0.526852857, 3.042136418) -0.8737, -0.5310, 3.0362 3 (-0.878093617, -0.17561706, 2.796928371) (-0.8725, -0.1810, 2.7873) 4 (-0.878093568, 0.175617049, 2.796928371) (-0.8725, -0.5310, 3.0362) 5 (-0.878093628, 0.526852825, 3.042136418) (-0.8700, 0.5175, 3.0176) 6 (-0.878093611, 0.878093569, 3.53255484) (-0.8675, 0.8675, 3.4966) 7 (-0.52685281, -0.878093573, 3.042136418) (-0.5310, -0.8737, 3.0362) 8 (-0.526852791, -0.526852812, 2.551717996) (-0.5268, -0.5268, 2.5517) 9 (-0.526852794, -0.175617069, 2.306509949) (-0.5352, -0.1773, 2.3056) 10 (-0.526852803, 0.17561705, 2.306509949) (-0.5215, 0.1700, 2.2969) 11 (-0.526852866, 0.526852858, 2.551717996) (-0.5200, 0.5200, 2.5367) 12 (-0.526852827, 0.878093628, 3.042136418) (-0.5175, 0.8700, 3.0176) 13 (-0.175617027, -0.878093592, 2.796928371) (-0.1810, -0.8725, 2.7873) 14 (-0.175617075, -0.526852805, 2.306509949) (-0.1773, -0.5252, 2.3056) 15 (-0.175617028, -0.175617046, 2.061301903) (-0.1756, -0.1756, 2.0613) 16 (-0.175617006, 0.175617028, 2.061301903) (-0.1756, 0.1756, 2.0559) 17 (-0.175617067, 0.526852847, 2.306509949) (-0.1700, 0.5215, 2.2969) 18 (-0.175617056, 0.878093620, 2.796928371) (-0.1675, 0.8700, 2.7759) 19 (0.175617041, -0.878093592, 2.796928371) (0.1675, -0.8700, 2.7759) 20 (0.175617068, -0.526852815, 2.306509949) (0.1700, -0.5215, 2.2969) 21 (0.175617061, -0.175617041, 2.061301903) (0.1715, -0.1715, 2.0559) 22 (0.175617051, 0.175617057, 2.061301903) (0.1756, 0.1756, 2.0613) 23 (0.175617054, 0.526852783, 2.306509949) (0.1773, 0.5252, 2.3056) 24 (0.175617039, 0.878093590, 2.796928371) (0.1810, 0.8725, 2.7873) 25 (0.526852762, -0.878093542, 3.042136418) (0.5175, -0.8700, 3.0176) 26 (0.526852876, -0.526852885, 2.551717996) (0.5200, -0.5200, 2.5367) 27 (0.526852785, -0.175617064, 2.306509949) (0.5215, -0.1700, 2.2969) 28 (0.526852823, 0.175617064, 2.306509949) (0.5252, 0.1773, 2.3056) 29 (0.52685278, 0.52685282, 2.551717996) (0.5268, 0.5268, 2.5517) 30 (0.526852781, 0.878093598, 3.042136418) (0.5310, 0.8737, 3.0362) 31 (0.878093552, -0.878093598, 3.53255484) (0.8675, -0.8675, 3.4966) 32 (0.87809359, -0.526852814, 3.042136418) (0.8700, -0.5175, 3.0176) 33 (0.878093528, -0.175617038, 2.796928371) (0.8700, -0.1675, 2.7759) 34 (0.878093604, 0.175617041, 2.796928371) (0.8725, 0.1810, 2.7873) 35 (0.878093617, 0.526852807, 3.042136418) (0.8737, 0.5310, 3.0362) 36 (0.878093597, 0.878093581, 3.53255484) (0.8781, 0.8781, 3.5326) 以上实验数据是在固定分辨率策略下设置中心点数为200, 可以算出可识别波长分辨率为0.142, 小于函数实际的波长分辨率, 因此可找到所有极值点.另外, 相对于蚁群算法求出所有极值点花费的3 203.297 s, P-MQHOA只需0.123 s, 并且其精度相对于蚁群算法的0.01精确到0.0000005, 无论是时间还是精度上P-MQHOA都占很大优势.实际上, 这里每个中心点只生成20个高斯采样点也能找出所有极值点, 时间花费可以更少, 因此P-MQHOA性能上还有进一步提高的空间.
综上讨论, P-MQHOA适合函数的全极值点寻优, 并且可以根据工程具体需求采用不同的策略:当需要查找所有极值点时可以采用固定分辨率策略, 当对时间要求较高时可以采用多级分辨率策略.实际上实验还发现这两种策略所找到的极值都包含全局最优解, 因此在不考虑求所有极值的情况下, 对于全局优化问题可以考虑用多分辨率策略.
4.3 全局多峰优化
在控制与决策领域中, 例如无人机的紧密编队飞行控制应用、航空发动机的燃烧主动控制应用、导弹控制参数的自动调节应用等[15], 系统的代价函数普遍地存在多个全局极值, 如何能够使系统的代价函数稳定地工作在全局极值点, 最大地发挥系统的性能, 是工业控制界一个亟待解决的问题.大部分多峰优化智能算法致力于解决这个问题, 而基于划分的MQHOA是适合全极值点寻优的算法, 可以从众多极值点中筛选出最佳极值点.
本节实验选用14个函数[6]作为测试集合, 其中 $f_1\sim f_7$ 和 $f_{12}$ 为一维函数, $f_{14}$ 为三维函数, 其余都是二维函数.比较P-MQHOA与文献[6]中其他智能算法的的区别.这些算法包括:基于双目标差分进化的多峰优化算法 (Multimodal optimization using abiobjective differential evolution algorithm, MOBi-DE)、拥挤度差分进化算法 (Crowding DE, CDE)、基于形态的差分进化算法 (Speciation-based DE, SDE)、欧氏适应度距离半径的粒子群算法 (Fitness-Euclidean distance ratio PSO, FER-PSO)、基于物种形成原理的粒子群算法 (Speciation-based PSO, SPSO)、协方差矩阵适应进化算法 (Covariance matrix adaptation evolutionstrategy, CMA-ES)、固定小生境半径的协方差矩阵适应进化算法 (CMA-ESwith fixed niche radius, CMA)、双目标多群体遗传算法 (Biobjective multipopulation genetic algorithm, BMPGA)、基于小生境的带有精英策略的非支配排序的遗传算法 (Niching-basedNSGA-Ⅱ) 和环型拓扑粒子群优化算法 (Particle swarm optimization using aring topology, rpso).其中r2pso中每个粒子只与环形拓扑中右边的粒子交互, r3pso与环形拓扑中的左右两边的粒子均交互.如表 4所示, 所有算法都经过50次独立重复实验, 统计找到的全局峰数, 对结果取平均值, 然后对每个算法进行排名记录到括号中去, 最终统计总排名.其中ε用来控制极值点的精度, No.表示对应函数应有全局峰数. P-MQHOA算法对一维函数采用固定分辨率策略, 二维函数采用二级分辨率策略, 三维函数采用三级分辨率策略.
表 4 二维Griewank函数极值分布拟合结果Table 4 Fitting results of extremum values0 distribution of Griewank functionFunc. No. ε P-MQHOA MOBi-DE CDE SDE S-CMA CMA SPSO PER-PSO r2pso r3pso Niching-based BNPGA
NSGA-Ⅱf1 1 0.05 1(4) 1(4) 1(4) 1(4) 1(4) 1(4) 0.44(12) 0.82(10) 0.76(11) 0.84(9) 1(4) 1(4) f2 1 0.05 1(5) 1(5) 1(5) 1(5) 1(5) 1(5) 0.40(12) 1(5) 0.88(11) 0.96(10) 1(5) 1(5) f3 2 1E-6 2(2.5) 2(2.5) 2(2.5) 1.96(5) 1.92(7) 1.95(6) 1.38(9) 0.8(10) 0.48(12) 0.6(11) 2(2.5) 1.5(8) f4 5 1E-6 5(1.5) 5(1.5) 3.84(10) 4.70(6) 0.04(12) 0.6(11) 4.88(3) 4.84(4) 4.68(7) 4.74(5) 4.37(9) 4.64(8) f5 1 1E-6 1(6) 1(6) 0.72(12) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) 1(6) f6 5 1E-6 5(1.5) 5(1.5) 3.96(9) 4.6(7) 0(12) 0.64(11) 4.92(4) 4.96(3) 4.88(5) 4.72(6) 4.52(8) 3.56(10) f7 1 1E-6 1(5.5) 1(5.5) 0.8(11) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 1(5.5) 0.6(12) 1(5.5) f8 4 5E-4 4(1.5) 4(1.5) 0.32(12) 3.72(4.5) 3.72(4.5) 3.43(7) 0.84(11) 3.68(6) 2.92(9) 2.76(10) 3.74(3) 3.34(8) f9 2 1E-6 2(3) 2(3) 0.04(12) 2(3) 1.6(7) 2(3) 0.08(11) 1.96(6) 1.44(9) 1.56(8) 2(3) 1.24(10) f10 1 1E-6 1(2) 1(2) 0.52(9) 0.32(10) 0.06(12) 0.18(11) 0.56(8) 1(2) 0.88(6) 0.76(7) 0.94(5) 0.96(4) f11 18 5 E-2 17.5(1) 17.4(2) 11.39(12) 12.78(9) 12.04(10) 12(11) 14.36(7) 15.61(6) 15.95(5) 16.45(4) 16.94(3) 13.72(8) f12 6 0.05 6(1.5) 6(1.5) 5.56(6.5) 4.88(12) 5.56(6.5) 5.81(3) 5.6(5) 5.28(10) 5.52(8) 5.16(11) 5.36(9) 5.71(4) f13 36 1E-3 35.6(1) 35.4(2) 33.8(3) 23.8(7) 24.6(6) 23.6(8) 25.72(5) 21.8(12) 22.4(9) 22.2(10) 22.05(11) 31.76(4) f14 218 1E-3 163.56(2) 175.88(1) 152(3) 50.6(8) 0.6(12) 32.5(11) 70.12(6) 68.6(7) 40.6(10) 45.4(9) 92.76(5) 121.54(4) Total ranks 38 39 111 92 109.5 102.5 104.5 92.5 113.5 111.5 85.5 88.5 参与比较的算法主要是对差分进化算法 (DE)、遗传算法 (GA) 和粒子群优化算法 (PSO) 的各种改进算法, 属于进化算法的分枝, 也都是基于群体智能的随机并行优化算法, 通过模仿生物群体内个体间的合作与竞争产生的启发式群体智能来指导优化搜索.这些算法使用小生境来提高种群多样性避免早熟, 可以同时使用多个子种群来定位和跟踪多个最优解.小生境的半径和小生境的数量影响整体的运算速度和精度, 算法性能对该参数很敏感, 需要调试多次, 其中大部分改进算法需要较多控制参数才能达到比较满意的结果, 而P-MQHOA只需要设置与需求波长分辨率相对应的采样中心点数即可.根据表 4统计的排名数据绘制热力图, 如图 9所示, 其中同一行颜色越淡表示找到的全局最优点越多, 反之越少.比如对于 $f_1$ 函数, 由颜色深浅可见SPSO算法的寻优效果最差, 其他PSO算法表现也不佳, 而P-MQHOA与其他遗传算法均能找到对应的全局最优点.
从上图可看出, P-MQHOA与MOBi-DE对14个函数测试下所找到的最优解数目差不多, 而其他算法相比之下较为逊色.单独比较两种算法对这14个函数测试所花费的平均时间, 如表 5所示, 时间单位为秒.
表 5 P-MQHOA算法与MOBi-DE算法时间消耗 (s)Table 5 Time consumption of P-MQHOA algorithm and MOBi-DE algorithm (s)Func. f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 f11 f12 f13 f14 P-MQHOA 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.1 0.1 0.1 1.76 0.28 4.32 12.14 MOBi-DE 1.56 1.46 2.82 2.64 1.54 2.04 1.68 18.34 15.28 41.24 46.28 36.49 51.74 72.36 可见采取固定分辨率策略的P-MQHOA的一维函数全局寻优时间消耗比MOBi-DE算法低两个数量级, 采取多级分辨率策略的P-MQHOA的多维函数全局寻优时间消耗比MOBi-DE算法也低很多.这是因为P-MQHOA局部收敛速度要比MOBi-DE等大部分智能算法要快的多.例如P-MQHOA对 $f_{14}$ 全局寻优过程中, 能够在13秒内找到约30 000个极值点, 然后筛选出全局极值点.
通过以上实验对比进行分析, 归纳如下:
GA采用二进制编码, 需要进行交叉变异选择遗传操作, 比较繁杂, 并且以前的知识随着种群的改变被改变, 没有记忆性, 有较多敏感的控制参数, 收敛速度慢; PSO一般采用实数编码, 粒子通过内部速度变化进行更新, 没有GA复杂的遗传操作, 此外通过当前搜索得到最优点进行信息共享, 较优解被所有粒子保存, 具有记忆能力, 收敛速度相对GA较快但很容易陷入局部最优解, 算法不稳定; DE采用实数编码、基于差分的简单变异操作和一对一的竞争生存策略, 降低了遗传操作的复杂性, 其独特的记忆能力使其可以动态跟踪当前的搜索情况以调整其搜索策略, 控制参数较少, 收敛速度相对于GA与PSO都要快. P-MQHOA无需编码, 拥有固定分辨率与多级分辨率两种策略, 每个区域的搜索相互独立, 每个采样中心点具有对所在区域极值的记忆能力, 操作简单, 控制参数较少, 收敛速度最快.
5. 结论
为解决多峰优化问题, 本文改进了MQHOA, 提出P-MQHOA算法, 把MQHOA的全局优化特性分散在多个均匀划分区域, 可以快速在各个区域找到极值点, 从而解决多峰优化问题.当实际问题除了最优解还需要备择解时, 该算法为全极值点优化中提供次优解提供了条件.此外针对具体问题本文提出了固定分辨率和多级分辨率两种策略, 前者有很好的全极值搜索能力, 后者有较强的实时性, 两种策略都能用于解决全局多峰优化问题.后面的研究工作将集中在该算法在多级分辨率策略下的全极值寻优性能研究, 提高其处理高维问题的能力, 扩展该算法在工业上的应用.
-
表 1 检测算法评价指标
Table 1 The evaluation index of detection method
指标名称 简要描述 使用场景 敏感度 即真阳性率, 表示被正确检测出的阳性比率 检测和分类评测 特异度 即真阴性率, 表示被正确识别为阴性的比率 分类评测 假阳性数 每个扫描数据集的检测结果中的假阳性个数, PFs/scan 检测和分类评测 假阴性数 每个扫描数据集的检测结果中的假阴性个数, FNs/scan 分类评测 $A_{z}$ 受试者操作特性曲线(ROC)下的面积 分类评测 Dice系数 一种集合相似度度量函数, 公式为: $s = \frac{2\left| X \bigcap Y \right |}{\left | X\right| + \left | Y\right |}$ 分割评测 Hausdorff距离[105] 一种基于对象边缘的、两对象间几何距离的度量方法 检测和分割评测 MAD距离[106] 两对象间的平均绝对距离, 同样基于对象边缘进行计算 检测和分割评测 表 2 肺部病变区域检测方法的性能比较
Table 2 The performance comparison of detection methods for lung lesion region
作者 年份 数据集 敏感度(%) 特异度(%) 假阳性 $A_z$ DSC 响应时间 研究对象 算法类型 张永强等[12] 2012 私有良性:75slices
恶性: 47slices- - - - 0.88
0.834.25s
4.36s结节分割 形态结构 王青竹等[92] 2013 私有scans: 150
Slices: 8250; LR: 109897.05 - 9.21 0.995 - - 结节检测 统计模
式识别Carvalho等[13] 2014 公开LIDC-IDRI; 800scans 85.91 97.70 1.82FP/scan
0.008FP/slice0.8062 - 13.56min/scan 结节检测 混合方法 Camarlinghi等[6] 2012 公开LIDC; scans: 69;
nodules: 11485 - 25FP/exam - - - 结节检测 混合方法 高婷等[7] 2014 私有scans:66; ST:2mm;
LR: 85; Size: 3 $\sim$ 20mm95.29 - 12.90 0.85 - 12.56s 结节检测 形态学 孙娟等[86] 2014 私有slices: 30 98.2 - 8.8 - - 3.76ms/slice 结节检测 聚类方法 王凯等[69] 2014 公开LIDC; scans: 17;
LR: 30; Size: 3 $\sim$ 30mm100 - 6.7/scan - - - 结节检测 形态结构 申正文等[75] 2016 私有4DCT;
两组10相位- - - - 0.868
0.830- 肿瘤分割 跟踪方法 Camarlinghi等[6] 2012 公开LIDC; scans: 69
ST $<$ 2mm; LR: 11485 - 25FP/scan - - - 结节检测 混合方法 Choi等[5] 2012 公开LIDC; scans:32;LR:76;
Slices:5453;Size:3 $\sim$ 30mm94.1 - 5.45/scan 0.967 - - 结节检测 遗传算法 Elizabeth等[17] 2012 私有slices: 1714 92.31 - 94.92 - - - 结节检测 区域增长 Netto等[18] 2012 私有exams: 29; LR: 48 86 91 0.138FP/exam - - - 结节检测 聚类方法 Song等[61] 2013 私有scans: 23 - - - - 0.81 - 结节/肿块 概率方法 Badura等[57] 2014 公开LIDC-IDRI; LR: 551 95.5 - 9.15 - - - 结节分割 聚类方法 Guo等[62] 2014 私有scans: 7 - - - - 0.85 - 肿瘤分割 概率方法 Jacobs等[48] 2014 公开scans: 318; LR:339;
Size:5 $\sim$ 33mm80 - 1.0/scan - - 122s/scan 结节检测 混合方法 Santos等[16] 2014 公开LIDC; scans:28;LR:252 90.6 85 1.17/scan - - - 结节检测 混合方法 Wang等[94] 2015 私有scans:103; nodules:127 - - 4FP/scan FROC:
0.88- - 结节检测 混合方法 Yong等[70] 2014 私有scans: 50 - - - - 91.2 - 肿瘤分割 形态结构 Demir等[91] 2015 公开LIDC-IDRI; scans: 200
Slices: 2775898.03 - 2.45/scan - - - 结节检测 统计模
式识别李阳等[79] 2013 私有scans: 20 100 95 - - - - 结节检测 统计模
式识别Lassen等[82] 2015 公开LIDC-IDRI; LR: 59 - - - - 0.52 结节分割 混合方法 Lu等[83] 2015 公开LIDC; scans:294;LR:631 85.2 - 3.31 - - - 结节检测 混合方法 Messay等[21] 2015 公开LIDC-IDRI; scans:456 - - - - 77.8 结节分割 混合方法 Qiang等[84] 2015 私有scans: 25; LR: 280 88.09 - 10.32 - 89.8 - 结节分割 形态结构 Setio等[77] 2015 LP: 238 98.3 - 4.0/scan - - - 结节检测 统计模
式识别Erdal等[88] 2015 公开LIDC; slices:138 - - - 0.9679 - - 结节检测 统计模
式识别Sheeraz等[22] 2016 公开LIDC; scans: 84 95.31 99.73 - - - - 结节检
测、分类统计模
式识别Senthilkumar等[72] 2016 私有LR: 25 88 84.05 2.05/scan - - - 肿瘤检测 区域增长 Manikandan等[49] 2016 私有scans: 106; LR: 801 100 93 0.38/patient - - - 肿瘤检测 数学
形态学Shin等[96] 2016 公开, 数据来自文献[107-108] - - TPR =0.85;
3FP/scan0.95 - - 淋巴结和
IDL检测深度学习 Dou等[98] 2016 公开LIDC-IDRI; scans:888 92.2 - 8FP/scan - - - 结节检测 深度学习 Cheng等[99] 2016 公开LIDC-IDRI;
良恶性结节:各70090.8±5.3 98.1±2.2 - 0.984±0.015 - - 结节分类 深度学习 Arnaud等[100] 2016 公开LIDC-IDRI; scans: 888 90.1 - 4/scan 0.996 - - 结节检测 深度学习 Bram等[101] 2015 公开LIDC-IDRI; scans: 865 77(上限78) - 8FP/scan - - - 结节分类 深度学习
拿来主义Francesco等[102] 2015 公开NELSON; nodules:4026 - - - 0.868 - - 结节分类 深度学习
拿来主义Thomas等[103] 2014 私有200002Dimagepatches;
5个类别- - 误分类率:
5.3- - - 结节分类 深度学习
半监督Song等[73] 2016 公开LIDC-IDRI; LR: 850;
私有LR: 12196.2 - 9.1 - - 8s 结节, 肿瘤 区域增长 刘峡壁等[109] 2015 公开LISS; scans: 271 70.2 97.2 - - - - 9类征
象分类混合方法 Guo等[46] 2016 公开LISS; 2D GGO: 45;
3D GGO 19100 33.13 - - - - GGO分类 阈值法 -
[1] World Health Organization. Cancer:fact sheet[Online], available:http://www.who.int/mediacentre/factsheets/fs297/en/, May 1, 2017 [2] Cascio D, Magro R, Fauci F, Iacomi M, Raso G. Automatic detection of lung nodules in CT datasets based on stable 3D mass-spring models. Computers in Biology and Medicine, 2012, 42(11):1098-1109 doi: 10.1016/j.compbiomed.2012.09.002 [3] Saito H, Oshima M, Kiuchi R, Watanabe K, Kitagawa K, Masuda S. Relationship between pleural indentation on computed tomography scans and pleural invasion in small peripheral lung cancer of 2, cm in size or less. Kyobu Geka, 2009, 62(9):767-770 http://www.ncbi.nlm.nih.gov/pmc/articles/instance/3590446/figure/fig1/ [4] Han G H, Liu X B, Han F F, Santika I, Zhao Y F, Zhao X M, Zhou C W. The LISS-a public database of common imaging signs of lung diseases for computer-aided detection and diagnosis research and medical education. IEEE Transactions on Biomedical Engineering, 2015, 62(2):648-656 doi: 10.1109/TBME.2014.2363131 [5] Choi W J, Choi T S. Genetic programming-based feature transform and classification for the automatic detection of pulmonary nodules on computed tomography images. Information Sciences, 2012, 212:57-78 doi: 10.1016/j.ins.2012.05.008 [6] Camarlinghi N, Gori I, Retico A, Bellotti R, Bosco P, Cerello P, Gargano G, Lopez T E, Megna R, Peccarisi M, Fantacci M E. Combination of computer-aided detection algorithms for automatic lung nodule identification. International Journal of Computer Assisted Radiology and Surgery, 2012, 7(3):455-464 doi: 10.1007/s11548-011-0637-6 [7] 高婷, 龚敬, 王远军, 聂生东, 孙希文.检测肺结节的3维自适应模板匹配.中国图象图形学报, 2014, 19(9):1384-1391 doi: 10.11834/jig.20140916Gao Ting, Gong Jing, Wang Yuan-Jun, Nie Sheng-Dong, Sun Xi-Wen. Three dimensional adaptive template matching algorithm for lung nodule detection. Journal of Image and Graphics, 2014, 19(9):1384-1391 doi: 10.11834/jig.20140916 [8] Hu S C, Bi K S, Ge Q X, Li M C, Xie X, Xiang X. Curvature-based correction algorithm for automatic lung segmentation on chest CT images. Journal of Biological Systems, 2014, 22(1):1-28 doi: 10.1142/S0218339014500016 [9] Lai J, Wei Q. Automatic lung fields segmentation in CT scans using morphological operation and anatomical information. Bio-medical Materials and Engineering, 2014, 24(1):335-340 http://content.iospress.com/articles/bio-medical-materials-and-engineering/bme815 [10] Li Q, Sone S, DoiK. Selective enhancement filters for nodules, vessels, and airway walls in two-and three-dimensional CT scans. Medical Physics, 2003, 30(8):2040-2051 doi: 10.1118/1.1581411 [11] 刘阳. 基于肺部高分辨率CT影像的肺结节识别方法研究[博士学位论文], 东北大学, 中国, 2010.Liu Yang. Lung Nodule Recognition Methods Based on Pulmonary High-Resolution CT Images[Ph., D. dissertation], Northeastern University, China, 2010. [12] 张永强, 夏顺仁, 陈兴灿, 郝欣.基于CT图像的肺结节分割算法研究.医疗卫生装备, 2012, 33(11):19-22 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ynws201211010&dbname=CJFD&dbcode=CJFQZhang Yong-Qiang, Xia Shun-Ren, Chen Xing-Can, Hao Xin. Study on segmentation algorithm of pulmonary nodule based on CT images. Chinese Medical Equipment Journal, 2012, 33(11):19-22 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ynws201211010&dbname=CJFD&dbcode=CJFQ [13] de Carvalho Filho A O, de Sampaio W B, Silva A C, de Paiva A C, Nunes R A, Gattass M. Automatic detection of solitary lung nodules using quality threshold clustering, genetic algorithm and diversity index. Artificial Intelligence in Medicine, 2014, 60(3):165-177 doi: 10.1016/j.artmed.2013.11.002 [14] Iqbal S, Iqbal K, Arif F, Shaukat A, Khanum A. Potential lung nodules identification for characterization by variable multistep threshold and shape indices from CT images. Computational and Mathematical Methods in Medicine, 2014, Article ID 241647, DOI:10.1155/2014/241647 [15] Zhou S J, Cheng Y Z, Tamura S. Automated lung segmentation and smoothing techniques for inclusion of juxtapleural nodules and pulmonary vessels on chest CT images. Biomedical Signal Processing and Control, 2014, 13:62-70 doi: 10.1016/j.bspc.2014.03.010 [16] Santos A M, de Carvalho Filho A O, Silva A C, de Paiva A C, Nunes R A, Gattass M. Automatic detection of small lung nodules in 3D CT data using Gaussian mixture models, Tsallis entropy and SVM. Engineering Applications of Artificial Intelligence, 2014, 36:27-39 doi: 10.1016/j.engappai.2014.07.007 [17] Elizabeth D S, Retmin Raj C S, Nehemiah H K, Kannan A. Computer-aided diagnosis of lung cancer based on analysis of the significant slice of chest computed tomography image. IET Image Processing, 2012, 6(6):697-705 doi: 10.1049/iet-ipr.2010.0521 [18] Netto S M B, Silva A C, Nunes R A, Gattass M. Automatic segmentation of lung nodules with growing neural gas and support vector machine. Computers in Biology and Medicine, 2012, 42(11):1110-1121 doi: 10.1016/j.compbiomed.2012.09.003 [19] Khan Z F, Kannan A. Intelligent segmentation of medical images using fuzzy bitplane thresholding. Measurement Science Review, 2014, 14(2):94-101 http://www.degruyter.com/dg/viewarticle.fullcontentlink:pdfeventlink/$002fj$002fmsr.2014.14.issue-2$002fmsr-2014-0013$002fmsr-2014-0013.pdf/msr-2014-0013.pdf?format=INT&t:ac=j$002fmsr.2014.14.issue-2$002fmsr-2014-0013$002fmsr-2014-0013.xml [20] Özkan H, Osman O, Sahin S, Boz A F. A novel method for pulmonary embolism detection in CTA images. Computer Methods and Programs in Biomedicine, 2014, 113(3):757-766 doi: 10.1016/j.cmpb.2013.12.014 [21] Messay T, Hardie R C, Tuinstra T R. Segmentation of pulmonary nodules in computed tomography using a regression neural network approach and its application to the lung image database consortium and image database resource initiative dataset. Medical Image Analysis, 2015, 22(1):48-62 doi: 10.1016/j.media.2015.02.002 [22] Akram S, Javed M Y, Akram M U, Qamar U, Hassan A. Pulmonary nodules detection and classification using hybrid features from computerized tomographic images. Journal of Medical Imaging and Health Informatics, 2016, 6(1):252-259 doi: 10.1166/jmihi.2016.1600 [23] Wang J H, Li F, Li Q. Automated segmentation of lungs with severe interstitial lung disease in CT. Medical Physics, 2009, 36(10):4592-4599 doi: 10.1118/1.3222872 [24] Wei Y, Shen G, Li J J. A fully automatic method for lung parenchyma segmentation and repairing. Journal of Digital Imaging, 2013, 26(3):483-495 doi: 10.1007/s10278-012-9528-9 [25] Mansoor A, Bagci U, Xu Z, Foster B, Olivier K N, Elinoff J M, Suffredini A F, Udupa J K, Mollura D J. A generic approach to pathological lung segmentation. IEEE Transactions on Medical Imaging, 2014, 33(12):2293-2310 doi: 10.1109/TMI.2014.2337057 [26] Park S, Lee S M, Kim N, Seo J B, Shin H. Automatic reconstruction of the arterial and venous trees on volumetric chest CT. Medical Physics, 2013, 40(7):071906-1-071906-13 doi: 10.1118/1.4811203 [27] Cheng Y Z, Hu X, Wang Y D, Wang J K, Tamura S. Automatic centerline detection of small three-dimensional vessel structures. Journal of Electronic Imaging, 2014, 23(1):013007, DOI: 10.1117/1.JEI.23.1.013007 [28] Koller T M, Gerig G, Szekely G, Dettwiler D. Multiscale detection of curvilinear structures in 2-D and 3-D image data. In:Proceedings of the 5th International Conference on Computer Vision. Cambridge, USA:IEEE, 1995. 864-869 [29] Lorenz C, Carlsen I-C, Buzug T M, Fassnacht C, Weese J. Multi-scale line segmentation with automatic estimation of width, contrast and tangential direction in 2D and 3D medical images. In:Proceedings of the 1st Joint Conference on Computer Vision, Virtual Reality and Robotics in Medicine and Medical Robotics and Computer-Assisted Surgery. Grenoble, France:Springer, 1997. 233-242 [30] Frangi A F, Niessen W J, Hoogeveen R M, Van Walsum T, Viergever M A. Model-based quantitation of 3-D magnetic resonance angiographic images. IEEE Transactions on Medical Imaging, 1999, 18(10):946-956 doi: 10.1109/42.811279 [31] Sato Y, Westin C-F, Bhalerao A, Nakajima S, Shiraga N, Tamura S, Kikinis R. Tissue classification based on 3D local intensity structures for volume rendering. IEEE Transactions on Visualization and Computer Graphics, 2000, 6(2):160-180 doi: 10.1109/2945.856997 [32] Malladi R, Sethian J A, Vemuri B C. Shape modeling with front propagation:a level set approach. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1995, 17(2):158-175 doi: 10.1109/34.368173 [33] 贾同, 魏颖, 吴成东.基于几何形变模型的三维肺血管图像分割方法.仪器仪表学报, 2010, 31(10):2296-2301 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yqxb201010024&dbname=CJFD&dbcode=CJFQJia Tong, Wei Ying, Wu Cheng-Dong. 3D lung vessel image segmentation scheme based on geometric active contour model. Chinese Journal of Scientific Instrument, 2010, 31(10):2296-2301 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yqxb201010024&dbname=CJFD&dbcode=CJFQ [34] Boskamp T, Rinck D, Link F, Kümmerlen B, Stamm G, Mildenberger P. New vessel analysis tool for morphometric quantification and visualization of vessels in CT and MR imaging data sets. Radiographics, 2004, 24(1):287-297 doi: 10.1148/rg.241035073 [35] 陈刚, 吕煊, 王志成, 陈宇飞.肺CT图像的血管骨架化方法.计算机科学, 2013, 40(5):274-278 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jsja201305067&dbname=CJFD&dbcode=CJFQChen Gang, Lv Xuan, Wang Zhi-Cheng, Chen Yu-Fei. Vessel skeletonization method for lung CT images. Computer Science, 40(5):274-278 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jsja201305067&dbname=CJFD&dbcode=CJFQ [36] Zhou C, Hadjiiski L M, Sahiner B, Chan H P, Patel S, Cascade P N, Kazerooni E A, Wei J. Computerized detection of pulmonary embolism in 3D computed tomographic (CT) images:vessel tracking and segmentation techniques. In:Proceedings of the 2003 Medical Imaging. San Diego, USA:SPIE, 2003. 1613-1620 [37] Kawata Y, Niki N, Umetani K, Nakano Y, Ohamatsu H, Moriyama N, Itoh H. Stochastic tracking of small pulmonary vessels in human lung alveolar walls using synchrotron radiation micro CT images. In:Proceedings of the 2013 Medical Imaging:Biomedical Applications in Molecular, Structural, and Functional Imaging. Lake Buena Vista, USA:SPIE, 2013. [38] Kitamura Y Z Y, Ito W, Ishikawa H. Data-dependent higher-order clique selection for Artery-Vein segmentation by energy minimization. International Journal of Computer Vision, 2016, 117(2):142-158 doi: 10.1007/s11263-015-0856-3 [39] Feulner J, Kevin Zhou S, Hammon M, Hornegger J, Comaniciu D. Lymph node detection and segmentation in chest CT data using discriminative learning and a spatial prior. Medical Image Analysis, 2013, 17(2):254-270 doi: 10.1016/j.media.2012.11.001 [40] Dornheim L, Dornheim J, Rössling I. Complete fully automatic model-based segmentation of normal and pathological lymph nodes in CT data. International Journal of Computer Assisted Radiology and Surgery, 2010, 5(6):565-581 doi: 10.1007/s11548-010-0530-8 [41] Steger S, Ebert D, Erdt M. Lymph node segmentation in CT slices using dynamic programming. In:Proceedings of the 2011 IEEE International Symposium on Biomedical Imaging:From Nano to Macro. Chicago, USA:IEEE, 2011. 1990-1993 [42] Wang Y, Beichel R. Graph-based segmentation of lymph nodes in CT data. In:Proceedings of the 6th International Conference on Advances in Visua Computing. Las Vegas, USA:Springer, 2010. 312-321 [43] Chen Q, Quan F, Xu J J, Rubin D L. Snake model-based lymphoma segmentation for sequential CT images. Computer Methods and Programs in Biomedicine, 2013, 111(2):366-375 doi: 10.1016/j.cmpb.2013.05.019 [44] 张艳玲, 何鑫驰, 李立.基于最大类间方差与形态学的淋巴结图像分割.计算机科学, 2013, 40(8):296-299 http://edu.wanfangdata.com.cn/Periodical/Detail/jsjkx201308064Zhang Yan-Ling, He Xin-Chi, Li Li. Lymph node image segmentation algorithm based on maximal variance between-class and morphology. Computer Science, 2013, 40(8):296-299 http://edu.wanfangdata.com.cn/Periodical/Detail/jsjkx201308064 [45] Yu P C, Poh C L. Region-based snake with edge constraint for segmentation of lymph nodes on CT images. Computers in Biology and Medicine, 2015, 60:86-91 doi: 10.1016/j.compbiomed.2015.02.011 [46] Guo K, Liu X B, Soomro N Q, Liu Y. A novel 2D ground-glass opacity detection method through local-to-global multilevel thresholding for segmentation and minimum bayes risk learning for classification. Journal of Medical Imaging and Health Informatics, 2016, 6(5):1193-1201 doi: 10.1166/jmihi.2016.1934 [47] Moltz J H, Bornemann L, Kuhnigk J M, Dicken V, Peitgen E, Meier S, Bolte H, Fabel M, Bauknecht H-C, Hittinger M, KieSSling A, Pusken M, Peitgen H O. Advanced segmentation techniques for lung nodules, liver metastases, and enlarged lymph nodes in CT scans. IEEE Journal of Selected Topics in Signal Processing, 2009, 3(1):122-134 doi: 10.1109/JSTSP.2008.2011107 [48] Jacobs C, Van Rikxoort E M, Twellmann T, Scholten E T, De Jong P A, Kuhnigk J M, Oudkerk M, de Koning H J, Prokop M, Schaefer-Prokop C, van Ginneken B. Automatic detection of subsolid pulmonary nodules in thoracic computed tomography images. Medical Image Analysis, 2014, 18(2):374-384 doi: 10.1016/j.media.2013.12.001 [49] Manikandan T, Bharathi N. Lung cancer detection using fuzzy auto-seed cluster means morphological segmentation and SVM classifier. Journal of Medical Systems, 2016, 40(7):1-9 doi: 10.1007/s10916-016-0539-9.pdf [50] Ge Z, Sahiner B, Chan H P, Hadjiiski L M, Cascade P N, Bogot N, Kazerooni E A, Wei J, Zhou C. Computer-aided detection of lung nodules:false positive reduction using a 3D gradient field method and 3D ellipsoid fitting. Medical Physics, 2005, 32(8):2443-2454 doi: 10.1118/1.1944667 [51] 魏颖, 徐陆, 李翔, 李锐.结合多尺度圆形滤波与MS聚类的疑似结节分割.仪器仪表学报, 2016, 37(1):192-199 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yqxb201601026&dbname=CJFD&dbcode=CJFQWei Ying, Xu Lu, Li Xiang, Li Rui. Segmentation of suspected nodules by multi-scale circular filtering combined with MS clustering. Chinese Journal of Scientific Instrument, 2016, 37(1):192-199 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yqxb201601026&dbname=CJFD&dbcode=CJFQ [52] Kanazawa K, Kawata Y, Niki N, Satoh H, Ohmatsu H, Kakinuma R, Kaneko M, Moriyama N, Eguchi K. Computer-aided diagnosis for pulmonary nodules based on helical CT images. Computerized Medical Imaging and Graphics, 1998, 22(2):157-167 doi: 10.1016/S0895-6111(98)00017-2 [53] Kubo M, Kubota K, Yamada N, Kawata Y, Niki N, Eguchi K, Ohmatsu H, Kakinuma R, Kaneko M, Kusumoto M, Mori K, Nishiyama H, Moriyama N. CAD system for lung cancer based on low-dose single-slice CT image. In:Proceedings of the 2002 Medical Imaging. San Diego, USA:SPIE, 2002. 1262-1269 [54] Yamada N, Kubo M, Kawata Y, Niki N, Eguchi K, Omatsu H, Kakinuma R, Kaneko M, Kusumoto M, Nishiyama H. ROI extraction of chest CT images using adaptive opening filter. In:Proceedings of the 2003 Medical Imaging. San Diego, USA:SPIE, 2003. 869-876 [55] Antonelli M, Frosini G, Lazzerini B, Marcelloni F. A CAD system for lung nodule detection based on an anatomical model and a fuzzy neural network. In:Proceedings of the 2006 Annual Meeting of the North American Fuzzy Information Processing Society. Montreal, Canada:IEEE, 2006. 448-453 [56] 李翠芳, 聂生东, 王远军, 孙希文, 郑斌.改进的模糊聚类亚实质肺结节3维分割.中国图形学报, 2013, 18(8):1019-1030 http://www.wenkuxiazai.com/doc/a77b8ef9a45177232f60a2e4-2.htmlLi Cui-Fang, Nie Sheng-Dong, Wang Yuan-Jun, Sun Xi-Wen, Zheng Bin. Segmentation of sub-solid pulmonary nodules based on improved fuzzy C-means clustering. Journal of Image and Graphics, 2013, 18(8):1019-1030 http://www.wenkuxiazai.com/doc/a77b8ef9a45177232f60a2e4-2.html [57] Badura P, Pietka E. Soft computing approach to 3D lung nodule segmentation in CT. Computers in Biology and Medicine, 2014, 53:230-243 doi: 10.1016/j.compbiomed.2014.08.005 [58] 宋佳, 聂生东, 王远军, 鲁雯.一种基于三维特征的肺结节概率分割算法.生物医学工程学杂志, 2014, 31(4):771-776 doi: 10.7507/1001-5515.20140144Song Jia, Nie Sheng-Dong, Wang Yuan-Jun, Lu Wen. A probability segmentation algorithm for lung nodules based on three-dimensional features. Journal of Biomedical Engineering, 2014, 31(4):771-776 doi: 10.7507/1001-5515.20140144 [59] Kouzani A Z, Lee S L A, Hu E J. Lung nodules detection by ensemble classification. In:Proceedings of the 2008 IEEE International Conference on Systems, Man, and Cybernetics. Singapore:IEEE, 2008. 324-329 [60] Bharath A C, Kumar D. An improved bayesian network model based image segmentation in detection of lung cancer. In:Proceedings of the 2014 International Conference on Recent Trends in Information Technology (ICRTIT). Chennai, India:IEEE, 2014. 1-7 [61] Song Q, Bai J J, Han D F, Bhatia S, Sun W Q, Rockey W, Bayouth J E, Buatti J M, Wu X D. Optimal co-segmentation of tumor in PET-CT images with context information. IEEE Transactions on Medical Imaging, 2013, 32(9):1685-1697 doi: 10.1109/TMI.2013.2263388 [62] Guo Y, Feng Y M, Sun J, Zhang N, Lin W, Sa Y, Wang P. Automatic lung tumor segmentation on PET/CT images using fuzzy Markov random field model. Computational and Mathematical Methods in Medicine, 2014, 401201, DOI:10.1155/2014/40120 [63] Farag A A, El-Baz A, Gimel'farb G G, Falk R, Hushek S G. Automatic detection and recognition of lung abnormalities in helical CT images using deformable templates. In:Proceedings of the 2004 International Conference on Medical Image Computing and Computer-Assisted Intervention. Saint-Malo, France:Springer, 2004. 856-864 [64] Guo L, Zhang Y T, Zhang Z W, Li D Y, Li Y. An improved random walk segmentation on the lung nodules. International Journal of Biomathematics, 2013, 6(6):1350043, DOI: 10.1142/S1793524513500435 [65] 孙文燕, 董恩清, 曹祝楼, 郑强.一种基于模糊主动轮廓的鲁棒局部分割方法.自动化学报, 2017, 43(4):611-621 http://www.aas.net.cn/CN/abstract/abstract19039.shtmlSun Wen-Yan, Dong En-Qing, Cao Zhu-Lou, Zheng Qiang. A robust local segmentation method based on fuzzy-energy based active contour. Acta Automatica Sinica, 2017, 43(4):611-621 http://www.aas.net.cn/CN/abstract/abstract19039.shtml [66] 陈侃, 李彬, 田联房.基于模糊速度函数的活动轮廓模型的肺结节分割.自动化学报, 2013, 39(8):1257-1264 http://www.oalib.com/paper/4417598Chen Kan, Li Bin, Tian Lian-Fang. A segmentation algorithm of pulmonary nodules using active contour model based on fuzzy speed function. Acta Automatica Sinica, 2013, 39(8):1257-1264 http://www.oalib.com/paper/4417598 [67] 郑强, 董恩清.窄带主动轮廓模型及在医学和纹理图像局部分割中的应用.自动化学报, 2013, 39(1):21-30 http://www.wenkuxiazai.com/doc/034f2351cfc789eb172dc8ce.htmlZheng Qiang, Dong En-Qing. Narrow band active contour model for local segmentation of medical and texture images. Acta Automatica Sinica, 2013, 39(1):21-30 http://www.wenkuxiazai.com/doc/034f2351cfc789eb172dc8ce.html [68] Li X P, Wang X, Dai Y X, Zhang P B. Supervised recursive segmentation of volumetric CT images for 3D reconstruction of lung and vessel tree. Computer Methods and Programs in Biomedicine, 2015, 122(3):316-329 doi: 10.1016/j.cmpb.2015.08.014 [69] 王凯, 张煜, 刘哲星, 林炳权, 吴志强, 曹蕾.基于自适应体窗结构分析的肺结节检测方法.南方医科大学学报, 2014, 34(6):759-765 http://www.cqvip.com/QK/91170X/201406/50044641.htmlWang Kai, Zhang Yu, Liu Zhe-Xing, Lin Bing-Quan, Wu Zhi-Qiang, Cao Lei. Structural analysis based on adaptive window for pulmonary nodule detection. Journal of Southern Medical University, 2014, 34(6):759-769 http://www.cqvip.com/QK/91170X/201406/50044641.html [70] Yong J R, Qi S, Van Triest H J, Kang Y, Qian W. Automatic segmentation of juxta-pleural tumors from CT images based on morphological feature analysis. Bio-Medical Materials and Engineering, 2014, 24(6):3137-3144 [71] 王秋萍, 冯筠, 强永乾, 于楠, 邓蕾, 郭佑民.空洞检测算法自动化检测肺结节内空泡和空洞的可行性.中国医学影像技术, 2014, 30(9):1309-1313 http://www.cnki.com.cn/Article/CJFDTotal-LCFS201411011.htmWang Qiu-Ping, Feng Jun, Qiang Yong-Qian, Yu Nan, Deng Lei, Guo You-Min. Feasibility of automatic detection of cavity and vacuole in pulmonary nodules using cavity detection algorithm. Chinese Journal of Medical Imaging Technology, 2014, 30(9):1309-1313 http://www.cnki.com.cn/Article/CJFDTotal-LCFS201411011.htm [72] Krishnamurthy S, Narasimhan G, Rengasamy U. Three-dimensional lung nodule segmentation and shape variance analysis to detect lung cancer with reduced false positives. Proceedings of the Institution of Mechanical Engineers, Part H:Journal of Engineering in Medicine, 2016, 230(1):58-70 doi: 10.1177/0954411915619951 [73] Song J D, Yang C Y, Fan L, Wang K, Yang F, Liu S Y, Tian J. Lung lesion extraction using a toboggan based growing automatic segmentation approach. IEEE Transactions on Medical Imaging, 2016, 35(1):337-353 doi: 10.1109/TMI.2015.2474119 [74] Iwao Y, Gotoh T, Kagei S, Iwasawa T, De Sales Guerra Tsuzuki M. Integrated lung field segmentation of injured region with anatomical structure analysis by failure-recovery algorithm from chest CT images. Biomedical Signal Processing and Control, 2014, 12:28-38 doi: 10.1016/j.bspc.2013.10.005 [75] 申正文, 高圆圆, 张煜.基于星形先验和图割的肺部四维计算机断层扫描肿瘤自动分割.生物医学工程学杂志, 2016, 33(2):295-302 doi: 10.7507/1001-5515.20160050Shen Zheng-Wen, Gao Yuan-Yuan, Zhang Yu. Automatic segmentation of four dimensional computed tomography of lung tumor based on star shape prior and graph cuts. Journal of Biomedical Engineering, 2016, 33(2):295-302 doi: 10.7507/1001-5515.20160050 [76] 范立南, 李道静, 孙申申, 常朝海. Top-hat和Gabor filter在肺结节ROI自动检测中的应用研究.计算机工程与应用, 2012, 48(4):190-192 http://www.wenkuxiazai.com/doc/5495a990524de518964b7dc7.htmlFan Li-Nan, Li Dao-Jing, Sun Shen-Shen, Chang Chao-Hai. Applied research on automatic detection of lung nodules ROI based on Top-hat and Gabor filter. Computer Engineering and Applications, 2012, 48(4):190-192 http://www.wenkuxiazai.com/doc/5495a990524de518964b7dc7.html [77] Setio A A, Jacobs C, Gelderblom J, Van Ginneken B. Automatic detection of large pulmonary solid nodules in thoracic CT images. Medical physics, 2015, 42(10):5642-5653 doi: 10.1118/1.4929562 [78] Cao L, Wang K, Xing Q Q, Lin B Q, Zhang Y. Auto detection of lung ground-glass opacity nodules based on high-pass filter and Gaussian mixture model. Journal of Medical Imaging and Health Informatics, 2016, 6(2):320-327 doi: 10.1166/jmihi.2016.1700 [79] 李阳, 史东承, 王珂, 王燕, 魏艳芳.基于图像模式的肺结节识别.吉林大学学报(工学版), 2013, 43(S1):463-467 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jlgy2013s1098&dbname=CJFD&dbcode=CJFQLi Yang, Shi Dong-Cheng, Wang Ke, Wang Yan, Wei Yan-Fang. Lung nodule recognition based on image pattern. Journal of Jilin University (Engineering and Technology Edition), 2013, 43(S1):463-467 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jlgy2013s1098&dbname=CJFD&dbcode=CJFQ [80] Tan Y Q, Schwartz L H, Zhao B S. Segmentation of lung lesions on CT scans using watershed, active contours, and Markov random field. Medical Physics, 2013, 40(4):043502, DOI: 10.1118/1.4793409 [81] Alilou M, Kovalev V, Snezhko E, Taimouri V. A comprehensive framework for automatic detection of pulmonary nodules in lung CT images. Image Analysis and Stereology, 2014, 33(1):13-27 doi: 10.5566/ias.v33.p13-27 [82] Lassen B C, Jacobs C, Kuhnigk J M, van Ginneken B, van Rikxoort E M. Robust semi-automatic segmentation of pulmonary subsolid nodules in chest computed tomography scans. Physics in Medicine and Biology, 2015, 60(3):Article ID 1307, DOI:10.1088/0031-9155/60/3/1307 [83] Lu L, Tan Y Q, Schwartz L H, Zhao B S. Hybrid detection of lung nodules on CT scan images. Medical Physics, 2015, 42(9):5042-5054 doi: 10.1118/1.4927573 [84] Qiang Y, Zhang X H, Ji G H, Zhao J J. Automated lung nodule segmentation using an active contour model based on PET/CT images. Journal of Computational and Theoretical Nanoscience, 2015, 12(8):1972-1976 doi: 10.1166/jctn.2015.4216 [85] Gonçalves L, Novo J, Campilho A. Hessian based approaches for 3D lung nodule segmentation. Expert Systems with Applications, 2016, 61:1-15 doi: 10.1016/j.eswa.2016.05.024 [86] 孙娟, 王兵, 杨颖, 田学东.聚类分析在肺结节识别中的应用.计算机应用, 2014, 34(7):2050-2053 doi: 10.11772/j.issn.1001-9081.2014.07.2050Sun Juan, Wang Bing, Yang Ying, Tian Xue-Dong. Research on cluster analysis in pulmonary nodule recognition. Journal of Computer Applications, 2014, 34(7):2050-2053 doi: 10.11772/j.issn.1001-9081.2014.07.2050 [87] 邱实, 汶德胜, 冯筠, 崔莹.一种新策略的肺结节检测算法.电子学报, 2016, 44(6):1413-1419 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxu201606023&dbname=CJFD&dbcode=CJFQQiu Shi, Wen De-Sheng, Feng Jun, Cui Ying. A new strategy of lung nodules detection algorithm. Acta Electronica Sinica, 2016, 44(6):1413-1419 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxu201606023&dbname=CJFD&dbcode=CJFQ [88] Taşci E, Uğur A U. Shape and texture based novel features for automated juxtapleural nodule detection in lung CTs. Journal of Medical Systems, 2015, 39(5):1-13 doi: 10.1007/s10916-015-0231-5.pdf [89] Wang Z, Chen S C. New least squares support vector machines based on matrix patterns. Neural Processing Letters, 2007, 26(1):41-56 doi: 10.1007/s11063-007-9041-1 [90] de Carvalho Filho A O, Silva A C, de Paiva A C, Nunes R A, Gattass M. 3D shape analysis to reduce false positives for lung nodule detection systems. Medical and Biological Engineering and Computing. Berlin:Springer-Verlag, 2016. 1-15 [91] Demir Ö, Yilmaz Ç A. Computer-aided detection of lung nodules using outer surface features. Bio-Medical Materials and Engineering, 2015, 26(S1):S1213-S1222 doi: 10.3233/BME-151418 [92] 王青竹, 康文炜, 王斌.三维隐SVM算法设计及在胸CT图像病灶检测中的应用.模式识别与人工智能, 2013, 26(5):460-466 doi: 10.3969/j.issn.1003-6059.2013.05.007Wang Qing-Zhu, Kang Wen-Wei, Wang Bin. Design of 3D latent-SVM and application to detection of lesions in chest CT. Pattern Recognition and Artificial Intelligence, 2013, 26(5):460-466 doi: 10.3969/j.issn.1003-6059.2013.05.007 [93] Brown M S, Goldin J G, Suh R D, McNitt Gray M F, Sayre J W, Aberle D R. Lung micronodules:automated method for detection at thin-section CT-initial experience. Radiology, 2003, 226(1):256-262 doi: 10.1148/radiol.2261011708 [94] Wang B, Tian X D, Wang Q, Yang Y, Xie H Z, Zhang S Y, Gu L X. Pulmonary nodule detection in CT images based on shape constraint CV model. Medical Physics, 2015, 42(3):1241-1254 doi: 10.1118/1.4907961 [95] 常亮, 邓小明, 周明全, 武仲科, 袁野, 杨硕, 王宏安.图像理解中的卷积神经网络.自动化学报, 2016, 42(9):1300-1312 doi: 10.16383/j.aas.2016.c150800.htmlChang Liang, Deng Xiao-Ming, Zhou Ming-Quan, Wu Zhong-Ke, Yuan Ye, Yang Shuo, Wang Hong-An. Convolutional neural networks in image understanding. Acta Automatica Sinica, 2016, 42(9):1300-1312 doi: 10.16383/j.aas.2016.c150800.html [96] Shin H C, Roth H R, Gao M C, Lu L, Xu Z Y, Nogues I, Summers R M. Deep convolutional neural networks for computer-aided detection:CNN architectures, dataset characteristics and transfer learning. IEEE Transactions on Medical Imaging, 2016, 35(5):1285-1298 doi: 10.1109/TMI.2016.2528162 [97] Tajbakhsh N, Shin J Y, Gurudu S R, Hurst R T, Kendall C B, Gotway M B, Liang J M. Convolutional neural networks for medical image analysis:full training or fine tuning? IEEE Transactions on Medical Imaging, 2016, 35(5):1299-1312 doi: 10.1109/TMI.2016.2535302 [98] Dou Q, Chen H, Yu L Q, Qin J, Heng P A. Multi-level contextual 3D CNNs for false positive reduction in pulmonary nodule detection. IEEE Transactions on Biomedical Engineering, 2016, 64(7):1558-1567 https://www.ncbi.nlm.nih.gov/pubmed/28113302 [99] Cheng J Z, Ni D, Chou Y H, Qin J, Tiu C M, Chang Y C, Huang S, Shen D G, Chen C M. Computer-aided diagnosis with deep learning architecture:applications to breast lesions in US images and pulmonary nodules in CT scans. Scientific Reports, 2016, 6:Article No. 24454 [100] Setio A A A, Ciompi F, Litjens G, Gerke P, Jacobs C, Van Riel S J, Wille M M W, Naqibullah M, Sanchez C I, van Ginneken B. Pulmonary nodule detection in CT images:false positive reduction using multi-view convolutional networks. IEEE Transactions on Medical Imaging, 2016, 35(5):1160-1169 doi: 10.1109/TMI.2016.2536809 [101] Van Ginneken B, Setio A A, Jacobs C, Ciompi F. Off-the-shelf convolutional neural network features for pulmonary nodule detection in computed tomography scans. In:Proceedings of the 12th IEEE International Symposium on Biomedical Imaging (ISBI). New York, USA:IEEE, 2015. 286-289 [102] Ciompi F, De Hoop B, Van Riel S J, Chung K, Scholten E T, Oudkerk M, de Jong P A, Prokop M, van Ginneken B. Automatic classification of pulmonary Peri-Fissural nodules in computed tomography using an ensemble of 2D views and a convolutional neural network out-of-the-box. Medical Image Analysis, 2015, 26(1):195-202 doi: 10.1016/j.media.2015.08.001 [103] Schlegl T, Ofner J, Langs G. Unsupervised pre-training across image domains improves lung tissue classification. In:Proceedings of the 2014 International MICCAI Workshop on Medical Computer Vision. Cambridge, USA:Springer, 2014. 82-93 [104] Armato S G I I I, Mclennan G, Bidaut L, McNitt-Gray M F, Meyer C R, Reeves A P, Zhao B S, Aberle D R, Henschke C I, Hoffman E A, Kazerooni E A, MacMahon H, van Beek Ed. The lung image database consortium (LIDC) and image database resource initiative (IDRI):a completed reference database of lung nodules on CT scans. Medical Physics, 2011, 38(2):915-931 doi: 10.1118/1.3528204 [105] Huttenlocher D P, Klanderman G A, Rucklidge W J. Comparing images using the Hausdorff distance. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1993, 15(9):850-863 doi: 10.1109/34.232073 [106] Chalana V, Linker D T, Haynor D R, Kim Y M. A multiple active contour model for cardiac boundary detection on echocardiographic sequences. IEEE Transactions on Medical Imaging, 1996, 15(3):290-298 doi: 10.1109/42.500138 [107] Depeursinge A, Vargas A, Platon A, Geissbuhler A, Poletti P A, Müller H. Building a reference multimedia database for interstitial lung diseases. Computerized Medical Imaging and Graphics, 2012, 36(3):227-238 doi: 10.1016/j.compmedimag.2011.07.003 [108] Seff A, Lu L, Cherry K M, Roth H R, Liu J M, Wang S J, Hoffman J, Turkbey E B, Summers R M. 2D view aggregation for lymph node detection using a shallow hierarchy of linear classifiers. Medical Image Computing and Computer-Assisted Intervention 2014, 17(Pt1):544-552 https://arxiv.org/pdf/1408.3337.pdf [109] Liu X B, Ma L, Song L, Zhao Y F, Zhao X M, Zhou C W. Recognizing common CT imaging signs of lung diseases through a new feature selection method based on fisher criterion and genetic optimization. IEEE Journal of Biomedical and Health Informatics, 2015, 19(2):635-647 doi: 10.1109/JBHI.2014.2327811 [110] Christ P F, Elshaer M E A, Ettlinger F, Tatavarty S, Bickel M, Bilic P, Rempfler M, Armbruster M, Hofmann F, D'Anastasi M, Sommer W H, Ahmadi S A, Menze B H. Automatic liver and lesion segmentation in CT using cascaded fully convolutional neural networks and 3D conditional random fields. In:Proceedings of the 2016 International Conference on Medical Image Computing and Computer-Assisted Intervention. Athens, Greece:Springer, 2016. 415-423 [111] Lipton Z C, Kale D C, Elkan C, Wetzell R. Learning to diagnose with LSTM recurrent neural networks. arXiv:1511. 03677, 2015. [112] Hofmanninger J, Langs G. Mapping visual features to semantic profiles for retrieval in medical imaging. In:Proceedings of the 2015 IEEE Conference on Computer Vision and Pattern Recognition. Boston, USA:IEEE, 2015. 457-465 期刊类型引用(9)
1. 张浩,解艺旋,杨晓楠,岳晨阳,马明,夏立坤. 融合迭代自组织数据分析的肺结节自动检测算法. 计算机与数字工程. 2024(05): 1336-1340+1358 . 百度学术
2. 崔冬梅,魏国辉,张喜科,马志庆. 肺部肿瘤影像辅助诊断分类算法研究进展. 中国胸心血管外科临床杂志. 2022(07): 892-898 . 百度学术
3. 张战成,张大龙,罗晓清. 带有语义特征得分的肺结节检测方法. 计算机应用. 2020(03): 925-930 . 百度学术
4. 于胡飞,温景熙,辛江,唐艳. 基于生成对抗网络的医学数据域适应研究. 大数据. 2020(05): 45-54 . 百度学术
5. 杨举平. CT与普通放射检查对肺部孤立性球形病变的诊断对比. 影像研究与医学应用. 2019(02): 85-86 . 百度学术
6. 霍其润,李建武,陆耀,秦明. 基于变分的CT图像环形伪影校正. 自动化学报. 2019(09): 1713-1726 . 本站查看
7. 顾广海. 螺旋CT扫描技术在肺结节患者临床诊断中的应用效果分析. 中国医药指南. 2019(35): 28-29 . 百度学术
8. 李祥霞,李彬,田联房,张莉,朱文博. 基于稀疏表示和随机游走的磨玻璃型肺结节分割. 自动化学报. 2018(09): 1637-1647 . 本站查看
9. 于超,王璐,程道文. 基于图像降阶的纹理特征挖掘方法在医学上的应用. 东北师大学报(自然科学版). 2018(03): 53-57 . 百度学术
其他类型引用(31)
-