用于执行彩色多普勒超声成像的方法和系统.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112334072A(43)申请公布日2021.02.05(21)申请号201980043599.8(74)专利代理机构永新专利商标代理有限公司(22)申请日2019.06.2472002代理人刘兆君(30)优先权数据18290072.02018.06.28EP(51)Int.Cl.A61B8/06(2006.01)(85)PCT国际申请进入国家阶段日A61B8/08(2006.01)2020.12.28A61B8/14(2006.01)(86)PCT国际申请的申请数据PCT/EP2019/0665672019.06.24(87)PCT国际申请的公布数据WO2020/002171EN2020.01.02(71)申请人皇家飞利浦有限公司地址荷兰艾恩德霍芬(72)发明人B·张O·博纳富权利要求书2页说明书9页附图4页(54)发明名称用于执行彩色多普勒超声成像的方法和系统(57)摘要本发明提供一种彩色多普勒超声成像方法。所述方法包括采集多个稀疏彩色多普勒超声图像帧和B模式超声图像帧。所述彩色多普勒图像帧和B模式图像帧然后在借助于通过优化框架来求解所述图像帧的全视场中的流向量场而基于所述帧来估计向量流之前被预处理，在所述优化框架中，掩膜适于定义所述稀疏彩色多普勒图像帧内的观察到的彩色多普勒测量结果的稀疏区域。新彩色多普勒图像帧基于所述估计来生成并且被包括在输出图像中。CN112334072ACN112334072A权利要求书1/2页1.一种彩色多普勒超声成像方法，所述方法包括：采集多个稀疏彩色多普勒超声图像帧(110)和B模式超声图像帧(120)；对所述多个稀疏彩色多普勒超声图像帧和B模式超声图像帧中的当前稀疏彩色多普勒超声图像帧和B模式超声图像帧进行预处理(130)；基于经预处理的当前稀疏彩色多普勒超声图像帧和B模式超声图像帧来估计向量流(180)，其中，对向量流的所述估计包括通过将优化框架应用到所述当前稀疏彩色多普勒图像帧来求解所述图像帧的全视场中的流向量场(u)，其中，所述优化框架包括掩膜(M)，所述掩膜适于定义所述当前稀疏彩色多普勒图像帧内的观察到的彩色多普勒测量结果(m)的稀疏区域，使得所求解的向量场(u)与所述掩膜(M)内的所述观察到的彩色多普勒测量结果(m)兼容；基于所估计的向量流来生成新彩色多普勒图像帧(190)；并且生成包括所述新彩色多普勒图像帧的输出图像(200)。2.根据权利要求1所述的方法，其中，所述优化框架包括平滑项。3.根据权利要求1或2所述的方法，其中，所述优化框架包括零散度约束。4.根据权利要求1至3中的任一项所述的方法，其中，所述优化框架包括估计的边界速度。5.根据权利要求1至4中的任一项所述的方法，其中，所述掩膜(M)是逐像素加权掩膜，其中，所述加权掩膜针对所述当前稀疏彩色多普勒图像帧中的无观察到的彩色多普勒测量结果(m)存在的那些像素中的每一个像素具有零值。6.根据权利要求5所述的方法，其中，所述加权掩膜针对所述当前稀疏彩色多普勒图像帧中的存在观察到的彩色多普勒测量结果(m)的那些像素中的每一个像素具有在零和一之间的范围内的值，所述值反映关于每个相应像素处的所述观察到的彩色多普勒测量结果(m)的置信度水平。7.根据权利要求6所述的方法，其中，所述置信度水平是每个相应像素处的噪声水平的函数。8.根据权利要求1至7中的任一项所述的方法，其中，对输出图像的所述生成包括将所述新彩色多普勒图像帧叠加在所述当前B模式超声图像帧上。9.根据权利要求1至8中的任一项所述的方法，其中，对新彩色多普勒图像帧的所述生成包括将在所述图像帧的所述全视场中求解的所述流向量场(u)重新投影在沿着其已经采集了所述观察到的彩色多普勒测量结果(m)的超声波束取向(s)上。10.根据权利要求1至9中的任一项所述的方法，其中，所述预处理包括在所述当前稀疏彩色多普勒超声图像帧上执行相位展开(140)。11.根据权利要求1至10中的任一项所述的方法，其中，所述预处理包括在所述当前B模式超声图像帧上执行流边界分割(150)。12.根据权利要求11所述的方法，其中，所述预处理包括基于所述流边界分割来导出先前B模式图像帧(170)与所述当前B模式帧之间的边界(160)的运动。13.根据权利要求12所述的方法，其中，对所述边界(160)的所述运动的所述导出基于斑点跟踪。14.根据权利要求1至13中的任一项所述的方法，其中，所述多个稀疏彩色多普勒超声2CN112334072A权利要求书2/2页图像帧包括3D彩色多普勒超声图像体积。15.根据权利要求14所述的方法，其中，所述3D彩色多普勒超声体积通过以下各项来生成：采集多个空间位置处的稀疏彩色多普勒超声图像帧的集合；在所述稀