跟踪区域组织体积和肺功能变化使用图像配准

文摘

之前我们已经演示了24小时再分配和再吸收的支气管肺泡灌洗(BAL)流体传递到肺支气管镜的过程中正常的志愿者。在这个工作我们利用影像匹配过程关联液再分配和再吸收区域肺功能的变化。肺CT数据集从六个人类受试者被用于这项研究。每个主题在四个时间点扫描前后的平衡过程。执行图像配准对齐图像在不同的时间点和不同的通货膨胀水平。生成的致密的位移场进行跟踪组织体积变化和当地的实质组织定量揭示变形模式。登记的准确性评估通过测量具有里程碑意义的匹配错误,1毫米的。结果表明,quantitative-assessed体液同意与bronchoscopist-reported unretrieved BAL体积在整个肺(平方线性相关系数是0.81)。肺组织体积的平均差异在基线和24小时后是2%左右,这表明BAL肺部液体几乎是24小时后吸收。地区肺功能变化与BAL流体的存在,和区域功能回到基线的流体被重吸收。

1。介绍

支气管肺泡灌洗(BAL)具有重要的临床应用和通常用于诊断肺部疾病,如感染(1)、肺癌和肺间质疾病。在球过程中,液体喷到的一小部分通过支气管镜肺,然后想起考试。

极大的兴趣了解分布和解决的进度落下帷幕。凯利et al。2)用数字减法技术可视化的解剖分布包含低浓度的盐水radio-opaque染料。格尼等人,陈等人。3,4]表明,倾向于减少缺陷的程度和频率随着时间的推移,和清除后约24小时。加布等。5)观察到流体运动最终决议之前叶和肺之间,表明叶回到基线后24小时。此外,肺功能的变化由于BAL过程研究unretrieved BAL的影响是很重要的。克莱恩et al。6]表明,肺力学可以显著改变一个小时或更长时间后落下帷幕。然而,一些研究已经完成跟踪unretrieved BAL和量化的解决过程的区域水平的对肺通气功能的影响。

多排探头CT (MDCT)检测可用于获得多个静态屏息肺CT图像在不同的肺容积。将图像配准技术应用于CT数据,我们能够发现密集的变形场变换不同肺容积之间的肺。转换可以分析计算voxel-by-voxel密度变化,估计当地的肺组织扩张,使其他生物力学测量。当结合图像分割结果,功能性和生物力学测量可以报道肺癌、叶,任何任意形状的次区域基础(7]。

霍夫曼和Ritman [8)用CT密度测量计算区域组织和肺部的空气含量。这种方法反映了组织密度准确、区域通风有效估计(5,9]。克里斯腾森等人,莱因哈特等人估计的局部组织变形使用Jacobian-based通风措施7,10]。这些测量的组织内容和通风功能可以利用图像配准跟踪组织体积变化和力学性能随时间变化区域。

本文描述了一种基于图像配准的方法定量跟踪unretrieved BAL的解决过程和测量区域肺通气功能在24小时后BAL过程。我们评估登记通过跟踪具有里程碑意义的动作。我们表明,unretrieved落下帷幕逐渐吸收和产品内容返回到基线后24小时。另外,我们观察到局部组织通风函数返回基线状态后24小时。

本研究是最接近5]。他们都观察到该决议进展后落下帷幕。主要的区别在于,在5]在叶级BAL研究决议进展需要叶细分为每个数据集,虽然提出了一个方法来研究球决议在区域层面上利用图像配准技术。同时,当地的肺功能变化由于BAL中观察到这个工作。

2。材料和方法

2.1。图像数据集

综述了协议和爱荷华大学的机构审查委员会批准。CT数据集从六个健康人体。每个主题在四个时间点:扫描基线,立即post-lavage(30分钟)内,post-lavage post-lavage 4小时,24小时。在每个时间点的功能余气量(FRC)和肺活量(TLC)扫描。因此,每个主题经历了完全8 CT扫描在24小时内从四个阶段。每一对扫描收购西门子感觉64多探头行CT扫描仪(Forchheim、德国)在屏息扫描阶段相同。每个体积数据集收购一段间距0.5毫米和512×512的重建矩阵。平面像素间距大约是0.6毫米×0.6毫米。扫描重建使用B31f重建的内核。表1列出了数据集获得每个主题在四个阶段和每个扫描使用的名称。图1显示FRC的肺容积和TLC压力,和体积差异FRC在每个阶段和薄层色谱扫描了六科目。

受试者洗胃的中部叶和唇舌。每个主题总共收到了300毫升的整除;总量来标示。BAL液体也被检索和衡量bronchoscopist在这个过程;检索到的数量表示。然后受试者接受post-lavage扫描在30分钟内灌洗(post0相位扫描)。unretrieved BAL流体的体积可以通过减去计算体积从体积,检索所示(1)。计算为每个主题表所示2。注意主题2观察咳嗽大量BAL流体的肺。

2.2。方法概述

我们的目标是利用图像配准跟踪产品体积变化和组织通风后24小时内BAL过程变化区域。图2显示了注册用于分析。两种类型的登记进行CT数据集。Intra-phase登记注册FRC图像到TLC图像在一个阶段。使用这些结果来估计局部肺功能在每一个阶段,不同阶段之间的比较来衡量变化的函数。Inter-phase注册登记薄层色谱图像post-lavage阶段基线。这些结果用来跟踪局部组织(或真空)含量变化在四个不同的阶段。

2.3。组织容量评估

我们假设肺是空气和组织/血液的混合物(真空)。所以Hounsfield单位(胡)肺部CT图像中组织和空气含量的函数。胡的肺部CT图像、组织体积和风量air-tissue后可以估计混合模型由霍夫曼和Ritman [8]。组织体积在体素位置可以被估计为 在哪里体素的体积吗。同样,风量体素可以估计 在本文中,我们假设和。和介绍了符号的简单性和。区域组织体积计算

2.4。图像配准

登记的目的是找到使两幅图像的空间映射到对齐。让和代表两个三维图像注册。转换定义了如何从模板映像点映射到相应的点目标图像。在三维空间,让定义一个立体像素坐标在图像域。转换是一个向量值函数定义在目标图像的立体像素点阵给予其相应的位置在模板映像。

对于每一个主题,进行了两种不同的注册:intra-phase注册匹配FRC数据显示其薄层色谱数据在同一阶段和inter-phase注册匹配TLC post-lavage阶段基线TLC数据中的数据。不同的注册算法被用来执行intra-phase和inter-phase注册。

2.4.1。Intra-Phase注册:RTVP

正规化的组织体积和vesselness措施保护非刚性的登记(RTVP)算法(11,12从FRC)被用来估计转换TLC在同一阶段。RTVP算法的平方和最小化组织体积差异(SSTVD) [9,13,14和vesselness测量差异(SSVMD),利用提供的丰富的强度和形状信息。这种方法已被证明是有效的注册在肺部CT图像精度高(11,12]。

对于一对FRC和薄层色谱扫描在同一阶段,并购可以忽略的时间间隔和假设组织体积不变。因此,登记可由组织体积保存在两个图像。的平方的总和组织体积差异(SSTVD)相似成本函数(9)强度的变异占肺CT图像在呼吸。这种相似性准则最小化肺内的当地组织体积差扫描在不同压力水平。组织的体积CT扫描可以估计(2)。然后强度相似度量SSTVD被定义为(9] 在哪里和是模板和目标图像强度函数,分别。表示肺区域联盟的目标图像和变形模板映像。雅可比矩阵的变换估计当地的体积变化导致映射图像变形。因此,组织量的形象和是相关的。

作为小直径血管分支,从容器体素提供原始灰度信息很少贡献指导灰度相似性度量。为了更好地利用血管位置的信息,我们使用vesselness测量基于海赛矩阵的特征值的图像强度。海赛矩阵的特征值的大小顺序,Frangi vesselness函数(15)被定义为 与 在哪里区分平板状和管状结构,占偏离状的结构,管状结构和噪音之间的区别。,,控制vesselness测量的敏感性。本文实验使用,,。vesselness测量计算多尺度和选为最大响应。vesselness形象是新,可以被视为一种像估计的vesselness特性。了基于相似性度量,平方之和vesselness测量差异(SSVMD),是为了匹配相似vesselness模式在两个图像。鉴于和vesselness措施的图像和在位置分别形成新的代价函数

执行限制转换有助于产生生理上注册更有意义的结果。连续介质力学模型,如线性弹性可以用来调整转换。本文使用拉普拉斯算符来调整位移场在哪里。正则化项组成 在哪里和。使用线性弹性微分算子可以帮助平滑转换,并帮助消除突然位移场的变化。

最后,总成本被定义为一个线性组合的灰度度规,vesselness措施保护指标和拉普拉斯算子约束 常量和是权重调整的意义三个方面。在这部作品中,权重常数选择通过最小化三个独立的成本方面同时登记实验的基础上的数据集从六个科目。这些参数被设置为和对于所有intra-phase RTVP注册。

转换参数化使用立方b样代表变换: 在哪里描述了控制节点的位移和是一个三维的三次b样条基函数的张量积。空间多分辨率从粗到细的过程是用于注册提高速度、准确性和鲁棒性。总成本(10)是优化使用内存有限,拟牛顿最小化方法和边界(L-BFGS-B) [16)算法。基于函数的充分条件,保证当地吸水参数化均匀三次b样崔和李(提出的17),b样系数约束的转换维护两个图像的拓扑结构。

2.4.2。Inter-Phase注册:RIVP

因为落下帷幕的影响液体灌洗后24小时内解决,组织体积(真空卷)不同阶段之间的不同。因此,组织体积保持假设是无效的情况下从post-lavage阶段登记薄层色谱扫描基线薄层色谱扫描。差的平方和(SSD)定义的 这是用于inter-phase登记。SSD的潜在假设是,图像强度在两个图像之间对应点应该是相似的。这是真的当注册图像相同的形态。在这种情况下,如果图像完全映射,相应的强度应该是相同的,这意味着每个点相同的底层结构具有相同的强度值在两个图像注册。然而,考虑到CT变化强度由于不同的空气含量和落下帷幕的影响,肺内的灰度范围是不同地区两个薄层色谱图像获得在不同的阶段。平衡这些灰度范围差异,正常化的强度是必需的。直方图匹配过程之前使用固态硬盘驱动注册修改模板图像的直方图,以便它是类似于目标图像。

除了相似性度量是改变从SSTVD到SSD(直方图匹配之后),其他组件和inter-phase登记注册制度是相同的与intra-phase登记。inter-phase薄层色谱扫描匹配使用一个正规化的强度和vesselness保留非刚性的登记(RIVP)。变换估计通过最小化的平方的总和强度差异,vesselness测量差的平方和,拉普拉斯算子约束。总成本函数 摘要加权常数选择使用类似的标准,在intra-phase RTVP注册,设置为和对于所有inter-phase RIVP注册。

2.5。评估登记通过跟踪具有里程碑意义的运动精度

对于每一个主题,20解剖地标手动选择和跟踪所有八个图片在四个阶段评估登记的准确性。这些地标被选为可辨认的气道分歧点,如图3。转换决定从肺注册可以用来预测不同坐标系之间的具有里程碑意义的运动。具有里程碑意义的错误被定义为之间的欧几里得距离registration-predicted地标位置和其真实位置相同的图像坐标测量匹配的准确性。

对于每个扫描对基线数据,一套分布式的地标选为船树枝点也被定义。的地标性建筑形象第一次被选为分段血管树的分支,18]。半自动系统(19)是用于指导观察者发现的地标图像与相应的体素的形象。地标是选择在肺部。一个例子的点分布如图4。专家为每个基线扫描选定的100多个里程碑式的双一双六个科目。

2.6。跟踪肺组织体积变化

对于每一个主题,交叉注册所有post-lavage TLC图像映射到基准坐标系,并提供voxel-wise通讯post-lavage从基线阶段三个阶段。这些映射使评估组织体积变化对于一个给定的体素位置在四个阶段。

使用拉格朗日坐标系,图5显示了一个示例的一个地区的位置在post0TLC baseTLC变形的不同形状,post4TLC, post24TLC转换,,,分别。假设该地区baseTLC对应的红色矩形区域内红色曲线包围的后阶段在图的图像5。同一地区的卷四个阶段,,,。这些卷可以分解为组织体积分数和空气体积分数基于立方体内的意思是立体像素强度。组织体积计算,,和,分别。在这里来标示的雅可比矩阵的转换;,,,强度函数的四个图片。随着组织空气的比率增加,体素的CT强度增加(亮)。通过这种方式,我们能够跟踪组织体积的次区域穿过四个阶段使用inter-phase注册结果。

组织总额可以集成在一个地区使用(4)。在同一地区定义在基线坐标。让,,,代表组织总量从baseTLC post0TLC, post4TLC和post24TLC分别。他们计算 与此同时,我们定义,,的组织体积变化指数三个post-lavage阶段相比,基础阶段,并定义,和组织体积变化指数的比率。他们制定如下:

2.7。评估肺功能的雅可比矩阵

肺组织变形模式是一个索引来评估肺功能。在这部作品中,雅可比矩阵的行列式转换领域得到了图像配准是用来估计局部组织变形(7]。

雅可比矩阵的行列式(通常是简单地称为雅可比矩阵)(20.- - - - - -22)是一种测量估计逐点的扩张和收缩变形。在三维空间,让的向量值转换变形模板映像到目标图像。雅可比矩阵的变换在位置被定义为 使用拉格朗日坐标系的雅可比矩阵值对应于零扩张或收缩。局部组织扩张对应于一个大于1的雅可比矩阵和局部组织收缩小于1对应于一个雅可比矩阵。

在不同的阶段,雅可比矩阵估计与两个因素有关。第一个因素是unretrieved落下帷幕。根据(6),肺力学可以显著改变一个小时或更长时间后落下帷幕。第二个因素是呼吸的努力。如图1、盐液保留落下帷幕后的肺手术可能会导致降低呼吸努力,也就是说,小体积变化FRC TLC。更小的体积变化导致较低的平均肺功能。为了我们的分析集中在如何unretrieved BAL影响肺功能,我们消除呼吸努力因素的累积分布函数(CDF)通过计算雅可比矩阵的估计。

在我们的例子中,提供的雅可比矩阵描述了雅可比矩阵值将被发现的概率小于或等于给定值。直观地说,这是“面积曲线”概率分布的函数。它可以被视为“排名”信息:提供1意味着等级最高(最大的雅可比矩阵),同时提供值意味着等级下的(最小的雅可比矩阵)。在以下描述,我们称之为CDF的雅可比矩阵的秩的便利。无论不同的总量变化,肺功能对应高的地区,更高的排名。通过这种方式,我们消除呼吸努力因素虽然比较unretrieved BAL流体引起的肺功能变化。让,,,代表雅克比的秩估计从基地,post0, post4,分别和post24阶段。然后排名变化被定义为

2.8。预处理

预处理开始通过识别所有图片使用肺肺地区工作站2.0(维达诊断公司,爱荷华州的城市,IA)。baseTLC图像的每一个主题,一个自动叶分割算法(23)是用于薄壁组织区域划分为五个不同的叶。图片和面具downsampled倍在每个维度减少图像配准的计算时间。

预处理后,三个inter-phase注册映射TLC图像三个post-lavage阶段为每个主题进行薄层色谱基线图像使用RIVP登记。这些生成的转换是用来跟踪组织体积变化在四个阶段。然后四个intra-phase注册FRC图像映射到TLC图像在每个阶段为每个主题使用RTVP进行登记。由此产生的变换被用来估计区域肺变形模式在每个阶段。结合inter-phase转换,该地区肺功能是跟踪四个阶段。

RTVP和RIVP登记、多分辨率策略是用于优化过程。从低到高收益的空间分辨率和图像分辨率从八分之一提高两倍,直到完全解决。与此同时,一个网格的层次b样条空间使用由大变小。最好的b样条网格空间中使用的实验是8毫米。图片和图像网格空间的细化或者。

3所示。结果

3.1。评估登记的准确性

注册结果可以定性判断观察员通过视觉评估。可视化图像的强度匹配是一个直观的方法来确定该地区边界和相应的结构是一致的。的例子inter-phase和intra-phase登记结果如图6。数据6(一)- - - - - -6 (e)显示的匹配结果inter-phase登记:(a)和(b)显示冠状切片baseTLC和post0TLC数据集相同的主题,分别;FRC的畸形的切片图像TLC图像所示(c);(d)显示了彩色融合片之间(a)和(b)之前登记;和(e)节目之间的彩色融合片(a)和(c)后登记。匹配结果intra-phase登记baseTLC和baseFRC之间同样在图所示6 (f)- - - - - -6 (j)。很明显看到,除了肺边界,重要的结构,如气道,船,注册后和裂缝的位置是一致的。

对于每一个主题,登记精度定量评估通过跟踪运动20 airway-branch地标在TLC图像四个阶段,和FRC和TLC图像之间在同一个阶段。箱线图的里程碑式的跟踪误差inter-phase登记和intra-phase登记在六主题图所示7。具有里程碑意义的错误超过100的箱线图广泛分布vessel-branch地标基线intra-phase登记(翘曲baseFRC baseTLC)如图8。这些每个主题的意思是错误的1毫米。

3.2。比较Quantitative-Assessed和Bronchoscopist-Reported Unretrieved落下帷幕

根据实验协议,组织之间的体积变化的定量评估基线阶段和post0阶段应该与unretrieved BAL流体从bronchoscopist报告数据(5]。从(1),计算从bronchoscopist-reported卷灌输和检索在平衡过程。从(15),组织体积变化计算baseTLC和变形post0TLC基线坐标定义为。实验发现的关系进行了定量评估和临床数据。之间的相关系数和通过线性回归计算。

表2列表unretrieved BAL流体的体积和组织体积变化为每个主题。线性回归和相关系数图所示9(一个)。如上所述,主题2观察咳嗽大量BAL流体的肺,洗胃过程中并没有收集。因此,它是合理的把这一主题从线性回归分析。线性回归和相关系数从主题删除数据后2如图9 (b)。

3.3。跟踪全球和区域组织体积变化

使用(14),我们在相同的赫芝电地区组织体积计算在不同的阶段(定义在基线阶段)。替换与基线整个肺分割,我们追踪全球组织体积变化使用(15)。测量整个肺组织体积变化比率在不同阶段表中列出3。

对于每个扫描对基线数据,自动叶分割算法(23)是用于薄壁组织区域划分为五个不同的叶。替换在(14)和五个不同叶分割,我们跟踪lobe-based组织体积随时间变化的使用(15)。图10显示了lobe-based组织体积变化比率在不同阶段平均在6个科目。

3.4。跟踪区域肺功能变化

对于每个intra-phase登记,我们计算voxel-wise雅可比矩阵的变换来估计肺扩张和收缩,它反映了当地在呼吸过程中肺功能。为了观察区域肺功能的趋势,我们将肺分成30个长方形板从背侧到腹侧肺,肺,从顶点到基地。雅可比矩阵计算平均在每个板,并绘制在图11。每个颜色线显示了雅可比矩阵估计从一个intra-phase登记。

intra-phase和inter-phase登记相结合在一起,我们雅可比矩阵估计在不同阶段映射到相同的基准坐标进行比较。图12说明了雅可比矩阵及其等级(CDF)分布在baseTLC横向切片定义坐标。左列显示强度模式在四个不同的阶段;中间列显示相应的雅可比矩阵地图在每个阶段;和正确的列显示相应的雅可比矩阵的秩地图估计。

如图12 (d),该地区下unretrieved BAL(灌洗区域)CT强度远高于其他地区。这种强度差异使我们能够部分手动清洗区域。区域功能随时间变化可以通过跟踪观察雅可比矩阵估计的排名变化灌洗区域和non-lavage区域内四个阶段,如图13。

4所示。讨论

我们已经描述了一个方法研究区域的进展BAL决议须就产权证和肺功能变化。inter-phase和intra-phase登记取得了良好的准确性通过目视检查如图6,通过跟踪具有里程碑意义的运动。均值地标(气管分支点)跟踪误差在6个科目mm inter-phase登记毫米intra-phase登记。这些错误是subvoxel范围内,这表明注册结果能够描述肺变形在同一阶段,不同阶段之间可容忍的错误。

Inter-phase登记是用来跟踪任何赫芝肺组织体积变化地区通过不同阶段。在整个肺层面,bronchoscopist-reported unretrieved落下帷幕和quantitative-assessed组织体积变化从基线到post0阶段表中列出2为每个主题。我们注意到总是低于决议,这可能由于液体灌洗的时间间隔期间,第一个post-lavage扫描(30分钟)。线性回归分析证明好bronchoscopist-reported之间的相关性和quantitative-assessed,如图9(一个)。相关系数是斜率为0.83。删除数据后从主题2由于咳嗽unretrieved BAL测量不准确,相关系数显著增加斜率为0.98,如图9 (b)。这一分析表明,quantitative-assessed组织体积变化与bronchoscopist-measured数据高度相关,这表明该方法用于追踪组织体积变化通过inter-phase注册全球意义上的有意义。

inter-phase注册,我们还可以通过四个阶段跟踪组织体积变化。如表所示3,意味着组织体积变化比基线状态在post0阶段,那么它滴post4阶段,最后下降在post24阶段。这表明流体不断灌洗后在24小时内解决。24小时后,大多数灌洗液吸收和肺组织卷回到基线状态。

除了评估组织体积变化在全球层面,我们可以跟踪组织体积变化在基准数据集上定义的条件与inter-phase登记。例如,大叶性基线图像分割是用来跟踪每个叶的组织体积变化通过四个阶段。如图10,右中肺和左肺上有最大的组织灌洗后体积变化。这是因为灌洗液是交付给两叶。然而,其他三个叶也经验组织体积增加,这表明流体不同叶之间的重新分配。24小时后,大多数灌洗液吸收和组织体积的所有叶回到基线状态。

Intra-phase登记结果提供当地肺功能测量评估可比在每个阶段。将整个肺地区划分为30板从背侧到腹侧肺的高度,或者从顶点到基础肺长度、平均雅可比矩阵跟踪与inter-registration跨四个阶段每个板区域的结果。如图11雅可比矩阵,平均在每个板高基线,然后在post0急剧下降阶段,并增加在post4,最后返回到基线状态post24阶段。这是因为肺的流体,所有科目除了主题6灌洗后立即有较小的灵感努力,如图1。作为流体的影响解析,肺部呼吸恢复正常状态和有类似的努力,基线状态。因此,肺功能恢复具有相同的灵感的努力,显示为蓝线之间的重叠(基线)和紫色的线(post24阶段)在图11。我们也注意到,有每个意味着雅可比矩阵梯度曲线从背侧到腹侧肺的高度,从顶点到基地以及肺长度。肺体积变化不同,地区较大的肺功能(背侧肺和肺基地)比其他地区所受影响更大。

图12在横向切片显示当地的雅可比矩阵图。扭曲所有post-lavage图像基线后,相同的片上的变化在24小时跟踪。流体输送和解决强度图像中清晰可见。雅可比矩阵地图显示功能可以减少在post0和post4阶段,在post24阶段和恢复。雅可比矩阵的累积密度函数(CDF)给了我们一个雅可比矩阵的秩信息,消除了影响不同的呼吸的努力。它表明在post0和post4阶段,军衔non-lavage地区几乎不变,但灌洗地区明显减少。24小时后,每个地区的排名恢复。这种模式的变化也反映在图13。灌洗区域有6%左右功能等级灌洗后立即减少,而non-lavage地区有一个增加的对比。这些排名变化24小时后消失。

这项工作的结果同意先前的研究在3- - - - - -6]。参考文献(3- - - - - -5)表明,流体吸收逐渐在24小时内,6)展示了肺力学变化后落下帷幕。这些研究大多在全球层面,或取决于该地区为每个数据集提供的利益分割。我们的工作使用图像配准图像之间产生密集的位移场和能够跟踪流体内容改变随着时间的推移在任何给定地区定量。肺功能改变也分析了从注册产生的转换。

5。结论

本文描述了一种方法来研究流体决议进展和肺功能变化区域利用图像配准技术。结果表明注册能够描述肺运动在同一阶段,不同阶段。Inter-phase注册使我们能够跟踪在任何赫芝肺组织体积变化地区通过不同阶段。同时,intra-phase注册提供了本地肺功能映射在同一个阶段。结合inter-phase intra-phase注册,肺功能变化是通过跟踪不同阶段在当地地区水平。这种方法也可以应用于放疗设计跟踪疾病进展和帮助。

信息披露

e·a·霍夫曼和j·m·莱因哈特维达诊断的股东,公司。

确认

作者要感谢k·墨菲博士和b . van Ginneken提供软件生成和注释的iX地标。这项工作是HL079406赠款支持,HL64368 HL080285, EB004126, CA129022。

引用

1. a·j·w·亨德森,支气管肺泡灌洗。”儿童疾病档案,卷70,不。3、167 - 169年,1994页。视图:谷歌学术搜索
2. c·a·凯利,c . j . Kotre c·沃德et al .,“解剖分布的支气管肺泡灌洗液的评估数字减放射显影术,”胸腔,42卷,不。8,624 - 628年,1987页。视图:谷歌学术搜索
3. k·j·w·格尼w·c·哈里森西尔斯et al .,“支气管肺泡灌洗:影像学表现,”放射学,卷163,不。1,第74 - 71页,1987。视图:谷歌学术搜索
4. m·p·c . c . Chen Andrich, j . Shelhamer”异常在肺通气/灌注扫描引起的支气管肺泡灌洗,”核医学杂志》,34卷,不。11日,第1858 - 1854页,1993年。视图:谷歌学术搜索
5. l·m·加布k·m·贝克·e·j·r·范发现,g . w . Hunninghake j·m·莱因哈特和e·a·霍夫曼”效应的节段支气管肺泡灌洗的定量ct肺,“放射学学术,18卷,不。7,876 - 884年,2011页。视图:出版商的网站|谷歌学术搜索
6. 美国克莱因,w·卡尔扎伊,p .齐默尔曼et al .,“肺力学变化后纤维支气管肺泡灌洗在机械通气患者中,“重症监护医学,24卷,不。12日,第1293 - 1289页,1998年。视图:出版商的网站|谷歌学术搜索
7. j·m·莱因哈特k .叮,k .曹g·e·克里斯滕森e·a·霍夫曼和美国诉博达,“须就产权证估计当地的肺组织扩张氙CT措施具体的通风相比,“医学图像分析,12卷,不。6,752 - 763年,2008页。视图:出版商的网站|谷歌学术搜索
8. e·a·霍夫曼和e . l . Ritman”身体取向对地区的影响在狗和懒惰,肺扩张”应用生理学杂志卷,59号2、481 - 491年,1985页。视图:谷歌学术搜索
9. y阴、e·a·霍夫曼和c·l·林”质量保护非刚性的肺部CT图像,采用三次b样条登记”医学物理学,36卷,不。9日,第4222 - 4213页,2009年。视图:出版商的网站|谷歌学术搜索
10. g·e·克里斯坦森,j . h的歌,w,即Naqa, a和d低,“跟踪肺组织运动和扩张/压缩与逆一致的图像配准和肺量测定法,“医学物理学,34卷,不。6,2155 - 2163年,2007页。视图:出版商的网站|谷歌学术搜索
11. 曹k, k叮,g·e·克里斯坦森和j·m·莱因哈特”组织体积和vesselness测量保护肺CT图像的非刚性的登记,”医学成像2010:图像处理卷,7623学报学报p。762309年,圣地亚哥,加利福尼亚州,美国,2010年2月。视图:出版商的网站|谷歌学术搜索
12. 曹k, k叮,g·e·克里斯滕森m . l . Raghavan r . e . Amelon和j·m·莱因哈特”统一在灰度非刚性的肺部CT血管信息登记”诉讼的第四届国际研讨会在生物医学图像配准(WBIR 10)卷,6204在计算机科学的课堂讲稿施普林格,页1 - 12,2010。视图:出版商的网站|谷歌学术搜索
13. k .曹g·e·克里斯滕森k .丁和j·m·莱因哈特”Intensity-and-landmarkdriven,逆一致,b样条登记和分析肺图像”第二届国际研讨会上肺图像分析,页137 - 148,伦敦,英国,2009年9月。视图:谷歌学术搜索
14. 崔y阴,j . e·a·霍夫曼m . h . Tawhai和c·l·林,“模拟肺气流与科目的边界条件,”生物力学杂志,43卷,不。11日,第2163 - 2159页,2010年。视图:出版商的网站|谷歌学术搜索
15. a . f . Frangi w . j . Niessen k . l . Vincken和m . a . Viergever“多尺度船舶增强过滤”第11届国际研讨会论文集医学影像计算和计算机辅助介入(MICCAI ' 98)卷,1496在计算机科学的课堂讲稿,页130 - 137年,剑桥,质量,美国,1998年10月。视图:出版商的网站|谷歌学术搜索
16. r·h·伯德·卢j . Nocedal和c·朱,“有限记忆绑定约束优化的算法,”暹罗期刊在科学计算,16卷,不。5,1190 - 1208年,1995页。视图:出版商的网站|谷歌学术搜索
17. y . Choi和s·李”,吸水条件2 d和3 d均匀三次b样条函数,“图形化的模型,卷62,不。6,411 - 427年,2000页。视图:出版商的网站|谷歌学术搜索
18. h . Shikata说,e·a·霍夫曼和m . Sonka”自动分割肺血管树的3 d CT图像,”2004年医学成像:生理学、功能,从医学图像和结构,a . a . Amini和a . Manduca Eds。卷,5369学报学报,页107 - 116,圣地亚哥,加利福尼亚州,美国,2004年2月。视图:出版商的网站|谷歌学术搜索
19. k·墨菲,b . van Ginneken j·p·w·Pluim s克莱因和m .盯着“半自动定量参考标准建设evaluatoin肺CT的登记,”学报》第11届国际会议上医学影像计算和计算机辅助干预,第二部分(MICCAI ' 08)卷,5242在计算机科学的课堂讲稿,页1006 - 1013,纽约,纽约,美国,2008年9月。视图:出版商的网站|谷歌学术搜索
20. r·c·巴克先进的计算,麦格劳-希尔,圣路易斯,密苏里州,美国,第3版,1978年版。
21. g·e·克里斯滕森r·d·洛夫·m·米勒et al .,“光滑的拓扑属性解剖地图,”在医学成像信息处理c . Braillot, y Bizais, r·d·Paola Eds。第三卷计算成像和愿景,第112 - 101页,1995年。视图:谷歌学术搜索
22. g·e·克里斯坦森和h·j·约翰逊,“一致的图像配准,”IEEE医学成像,20卷,不。7,568 - 582年,2001页。视图:出版商的网站|谷歌学术搜索
23. 美国Ukil和j·m·莱因哈特Anatomy-guided肺叶在x射线CT图像分割,“IEEE医学成像,28卷,不。2、202 - 214年,2009页。视图:出版商的网站|谷歌学术搜索

版权

版权©2012 Kunlin曹等。这是一个开放的分布式下文章知识共享归属许可,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。