一种图像对称结构的提取方法
【技术领域】
[0001] 本发明涉及医学图像处理领域，尤其涉及一种医学图像对称结构的提取方法。
【背景技术】
[0002] 在医学图像处理中，图像的线性对称结构（如二维图像的对称线，三位图像的对 称面）的提取具有广泛的应用。
[0003] 美国专利申请文件US5533143提出了一种检测图像是否对称的方法，判断准则主 要依靠图像两边像素值之差是否为0。但是所述专利申请中，假设前提是图像的对称轴是 和通过图像中心的一条垂直的中心线一致的，这显然应用是很狭隘的。而且两边在差值为 〇也是很苛刻的条件，很容易受噪声、病变等影响。
[0004] 美国专利申请文件US6263097提出以图像两边像素灰度值的均方差为判断准则， 不断优化对称线的角度来得到最后的对称线。另一篇美国专利申请文件US7409085首先计 算出图像质心（xc，yc)，选定一个搜索区域（_15mm〈 = [x-xc]〈 = 15mm，-2mm〈 = [y-yc]〈 =2mm，-20。〈 = β〈 = 20。），计算此区域内可能性最大的一条对称线（x0，y0, β0)，再重 定义搜索区域（_3mm〈 = [x_x0]〈 = 3mm，_3mm〈 = [y-y0]〈 = 3mm，-3。〈= [β-β0]〈 = 3。），重复计算可能性最大的对称线。显然上述方法是基于遍历寻找最优解的方法，与基于 迭代的最优化方法相比，在速度和准确性等方面都有着明显的不足。
[0005] 文献 Babak A. Ardekani, etc, 'Automatic Detection of the Mid-Sagittal Plane in3_D Brain Images'，IEEE transactions on medical imaging,vol.,16, N).6, Decemberl997.采用单纯形下降法优化求解了 3D头部图像的对称面，代价函数选取互相关 函数。考虑到医学图像的复杂性，并且代价函数一般为对称面方程参数的非线性函数，单纯 形法并不是非常适合的。
[0006] 美国专利申请文件US20070276219用优化方法寻找相互正交的三个面，其中第一 个平面为所求的对称面。代价函数选取三个面与各自相平行的临近数个面对应点的灰度熵 值，采用基于单纯形的Nelder-Mead优化方法。美国专利US20117986823目的是为了保证从 多次扫描检测出的对称面的一致性（Consistent),把医生（MR设备operator)标定的MSP 作为金标准，然后以互信息为代价函数衡量自动检测的MSP同金标准的差异。实现上需要 专家检测MSP，如果不同病人的MSP第一次扫描（localizer)均需人工检测时，操作上会非 常麻烦。
【发明内容】
[0007] 本发明解决的问题是提供一种医学图像的对称结构的提取方法，用以对医学图像 的对称结构进行提取。
[0008] 为了解决上述问题，本发明提供了一种医学图像对称结构的提取方法，包括：提 供医学图像及初始对称，根据所述初始对称将所述医学图像划分为位于所述初始对称两侧 的两部分；基于所述两侧的像素信息建立基于对称参数的目标函数；迭代优化计算对称参 数，进而提取所述对称参数对应的所述医学图像对称结构。
[0009] 可选的，所述初始对称包括初始中心线和初始中心面：若所述医学图像为二维图 像，则所述初始中心线为对称线方程x+k*y+c = 0,其中，k，c为待定的对称参数；若所述医 学图像为三维图形，则所述初始中心面为对称面方程X = kiy+k2z+b，其中，Kl，K2，b为待定 的对称参数。
[0010] 可选的，所述迭代优化包括：利用lucas-kanade算法、梯度下降法或者BFGS方法 对所述目标函数优化迭代。
[0011] 可选的，建立所述目标函数包括：基于所述两侧像素的均方差、互相关或者互信息 建立所述目标函数。
[0012] 可选的，基于所述目标函数计算对称参数还包括：对所述目标函数进行泰勒展开， 以获取增量目标函数。
[0013] 可选的，包括：通过迭代方法求解所述对称参数，每次迭代中，获取令所述增量目 标函数值最小对应的增量，并将所述增量后的对称参数替换当前对称参数，进行迭代优化。
[0014] 可选的，所述目标函数为C=E Plu[I(R(Pl|Q))-I(Pl)]2,其中，其中，I(P 1)为一 侧的点P1的像素值，L为遍历所述点P1所在一侧的像素值，I (R(P1IQ))为与P1对称的点 R(P11Q)的像素值，所述对称为以当前对称参数Q的对称。
[0015] 可选的，通过迭代方法求解对称参数Q，及每次迭代使得所述目标函数最小的增量 Λ Q，并令增量后的对称参数Q+ △ Q替换当前对称参数Q进行迭代优化。
[0016] 可选的，还包括采用泰勒展开对所述目标函数进行 展开：若所述医学图像为二维图像，展开的目标函数表达式为
t其中，
b医学图像 在点R(P1IQ)的偏导数；若所述医学图像为三维图像，则所述展开的目标函数为
^其图像I在点 R(P1Iq)的偏导数。
[0017] 可选的，所述目标函数在AQ处具有极小值，即所述目标函数对AQ的偏导为0,若 所述医学图像为二维图像，所述目标函数对AQ的偏导的表达式如下：
[0018]
突取 AQ = H 1S,其 中
?
[0019] 若所述医学图像为三维图像，则所述目标函数对AQ的偏导的表达式如下
获取AQ = H 1S,其中
[0020] 可选的，若所述医学图像为二维图像，设像素沿对称线x+k*y+c =〇的对称点的变换为R(RX，Ry)，对于点P1U, y)则其对称变换为：
x-k,*v-k"*z-h = 0的对称占的亦拖为RiR . R . R ) ·对干占 njx，y，y，）则其对称变换
i中 F = x-k^y-k^z-b，KF =1+V+k/，并可求得对称变换Rx，Ry，Rz对参数kp k2, b的偏导数。
[0022] 可选的，所述迭代的终止条件为：当所述对称参数的增量的模小于最小迭代步长 时或者迭代次数不小于最大迭代次数时，所述迭代终止。
[0023] 本发明对二维图像选取对称线方程x+k*y+c = 0,对于三维图像选取对称面方程X =kiy+k2z+b，借鉴类似图像配准的迭代优化方法，利用对称坐标方程将图像分成位于对称 线或对称面两侧的左、右两部分，基于两边像素点的目标函数作为代价函数进行优化迭代， 上述优化方法具有更加准确、快速的优点，并且鲁棒性更好。
【附图说明】
[0024] 图1所示为本发明一个实施例的医学图像的对称结构提取方法的方法流程示意 图。
【具体实施方式】
[0025] 在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以 很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况 下做类似推广，因此本发明不受下面公开的具体实施的限制。
[0026] 其次，本发明利用示意图进行详细描述，在详述本发明实施例时，为便于说明，所 述示意图只是实例，其在此不应限制本发明保护的范围。
[0027] 为了解决上述问题，本发明提供了一种医学图像对称结构的提取方法，包括：提 供医学图像及初始对称，根据所述初始对称将所述医学图像划分为位于所述初始对称两侧 的两部分；基于所述两侧的像素信息建立基于对称参数的目标函数；迭代优化计算对称参 数，进而提取所述对称参数对应的所述医学图像对称结构。建立所述目标函数包括：基 于所述两侧像素的均方差、互相关或者互信息建立所述目标函数。所述迭代优化包括利用 lucas-kanade算法、梯度下降法或者BFGS方法对所述目标函数优化迭代。
[0028] 进一步地，通过迭代方法求解所述对称参数，每次迭代中，获取令所述增量目标函 数的值最小对应的增量，并将所述增量后的对称参数替换当前对称参数，进行迭代优化。基 于所述目标函数计算对称参数包括：对所述目标函数进行泰勒展开，以获取增量目标函数。
[0029] 其中，所述初始对称包括初始中心线和初始中心面：若所述医学图像为二维图像， 则所述初始中心线为对称线方程x+k*y+c = 0,其中，k，c为待定的对称参数；若所述医学 图像为三维图形，则所述初始中心面为对称面方程X = kiy+k2z+b，其中，Kl，K2，b为待定的 对称参数。
[0030] 所述目标函数为C=E P1eJI(R(P1Iq))-I (P1)]2,其中，其中，I(P1)为一侧的点P1 的像素值，L为遍历所述两侧中的其中一侧的像素值，I (R(P1 Iq))为与P1对称的点R(P1 |q) 的像素值，所述对称为以当前对称参数Q的对称。
[0031] 具体地，通过迭代方法求解对称参数Q包括：每次迭代获取使得C = Σ Pl ^ [I (R(P11Q+ Λ Q)) -I (P1) ]2最小的增量Λ Q，通过求解获取的增量后的对称参数Q+ Λ Q 替换当前对称参数Q，进行迭代优化。
[0032] 采用泰勒展开对所述目标函数进行优化：若所述医学图像为二维图像，优化的目 标函数表达式为
医学图像在点R(P1IQ)的偏导数；若所述医学图像为三维图像，则所述优化的目标函数 为
I图像I在点 R(P1Iq)的偏导数。
[0033] 如图1所示为本发明一个实施例的医学图像的对称结构提取方法，包括：执行步 骤S1，输入图像，计算灰度质心并初始化参数；执行步骤S2,判断像素点是否位于左侧。所 述左侧为位于当前对称参数对应的对称线或者对称面的左侧。
[0034] 上述判断步骤S2中若为否，则执行步骤S3,输入下一点像素点，并继续后续的判 断和方法流程，即执行步骤S2,判断像素点是否位于左侧。
[0035] 若所述判断步骤S2为是，则执行步骤S4,计算原坐标点与对称点灰度差，及图像 在对称点梯度，并将对称变换对系数求偏导，累加 S及H矩阵，具体地，S，H矩阵将下述具体 实施例详述。
[0036] 接着执行步骤S5,判断像素是否遍历完毕。若判断为是，则执行步骤S6,求H逆矩 阵，并更新系数。若判断为否，则执行步骤S3,输入下一点像素点。并继续执行步骤S3的后 续步骤，