专利名称：一种早期检测特发性脊柱侧凸的智能增强系统及其方法
技术领域：
本发明涉及一种利用胸部放射成像早期检测特发性脊柱侧凸的智能增强系统及其方法。
背景技术：
随着医学影像管理与通讯系统(PACS)的出现，逐渐积累的大量数字图像数据库和数字地图更突出了获取图像内容的问题。在传统的医学图像数据库中，大多数索引和检索方法都基于文本信息如研究身份、患者名字、日期等。但是，在很多图像数据库中，信息本身不是字母和数字的，而是大量的图像特征诸如颜色、结构和形状特征，这些特征都属于图像内容而不能基于传统的数据库的检索方法而找到。
事实上，在很多研究、临床和教学实践中，如果图像数据能够基于图像内容而检索到，则更多的信息能够被利用从而促进某些方法学的图像的回顾性和展望性研究。这些图像内容的检索系统如果能很好的开发将有助于临床大夫和教育者以及科研工作者通过自动匹配所查询图像与PACS图像数据库中属于相近类别或相似特征的图像而找到具有类似特征的图像。
医学图像，如胸部影像图像的分类是研究人员公认的难题，因为这些医学图像极其相似，它们通常只在一些微小的细节上有差别。目前的基于特定图像特征如结构、颜色、形状等进行基于内容的索引技术，应用于医学图像上还不够准确。
2000年以来，眼科视光学及放射科学系自主研发了一种影像管理与通讯系统。从2001年以来，该PACS应用于临床上的放射成像临床学中。从此，该图像数据库得到了快速发展。PACS中的现有索引系统是为日常系统的运行而设计的。然而，当我们的关注点为具体图像特征或形状时，该索引系统则难以满足研究与教学目的。
在香港，学生的脊柱侧凸的发病率约为2.7％。在放射成像临床中，通过大学生的胸部X片也检测到特发脊柱侧凸。特发脊柱变形的检测似乎只是整形外科医生和社区健康医生关心的问题。在大多数情况下，脊柱侧凸的影像不会在第一次就诊时就被分出。这将使青少年的脊柱侧凸的随后研究和回顾性研究复杂化。
现有PACS基本结构和脊柱侧凸的现有病例促进了智能型图像内容检索和分类系统的开发。该系统的开发有助于社区医疗，因为目前仍缺乏一种较好的整合了PACS的图像内容检索系统。因此，一种全新的基于图像匹配算法和计算机智能方法的图像内容检索系统显示了开发的良好前景。
发明内容
本发明的目的在于提供一种利用胸部放射成像早期检测特发性脊柱侧凸的智能增强系统，该智能增强系统包括a.图像处理�？椋靡曰竦眯夭坑跋窆庹な荩籦.向量化�？椋靡越枷裰械募怪蛄炕玫郊怪兰浦行牡悖籧.代价函数计算�？椋靡约扑阃枷竦募怪蛄课呛现毕叩某潭龋玫骄轿蟛睿籨.判定�？椋靡耘卸暇轿蟛钣脬兄礡的大小，从而将图像分类为正常脊柱和脊柱侧凸。
其中，向量化�？榭梢园ㄏ喙胤椒ǖピ�、�：兰频ピ把≡竦ピ鲅≡竦ピ挥谙喙胤椒ǖピ湍：兰频ピ埃靡匝≡窠邢喙胤椒ɑ故悄：兰�。
本发明还提供了利用上述智能增强系统进行早期检测特发性脊柱侧凸的方法，该方法包括如下步骤(1)处理图像，获得胸部影像的光栅数据；(2)将图像的脊柱向量化，得到脊柱估计中心点；(3)计算代价函数，得到均方误差；(4)将图像分类通过比较均方误差与阈值R的大小而将图像分为正常脊柱和脊柱侧凸。
其中，步骤(2)的向量化可以采用相关方法和�：兰疲⑶以谙喙胤椒ê湍：兰浦安捎醚≡癫街瑁貌街栌靡耘卸鲜遣捎孟喙胤椒ɑ故悄：兰�。
本发明的智能增强系统是用于检测脊柱侧凸的发病和严重程度，其依据数字化胸部放射影像，该放射影像是基于�：呒璠1，2，3，5]，神经网络[4，6]和最小二乘法[7]。不像采用活动形状模型(Active ShapeModeling(ASM))方法那样对脊柱的轮廓一点一点地研究，采用本发明的系统时，只用一条垂直线表示脊柱，找到胸部放射影像的每条直线所匹配的脊柱估计中心点。该方法不仅快速简便而且易于自动定位胸部放射影像中的脊柱。此外还定义了代价函数，以便量化脊柱侧凸指数评价所估计的脊柱位点与垂直直线之间的偏差。所有的病例最初都被认为是正常的，并且代价函数都收敛到零。在系统的训练模式期间，通过专家输入数据而获得误差水平R值。在运行模式期间，当代价函数值低于误差水平时，病例被归为正常脊柱；相反，当代价函数值高于误差水平时，病例被归为不正常类，即脊柱侧凸，见图1智能系统工作概图。系统的实施包括两个计算过程(a)脊柱中心点的向量化；(b)代价函数的计算。
图1为智能系统工作概图；图2为智能复合方法的流程图；图3表示模板为胸部放射影像中的一条线段的像素值曲线图，3a为所选的线段，3b为该线段对应的曲线图；图4为传统神经模糊网络的结构；图5为表示一阶B-样条(B-Spline)函数的基函数；图6为智能方法的运行全图。
具体实施例方式
下面结合附图和具体实施例对本发明作进一步详细说明。
(a)脊柱中心点的向量化在检测脊柱侧凸的发病及评价脊柱侧凸的严重程度时，并不需要描绘脊骨的轮廓。找到数字胸部放射影像的每条线对应的脊柱中心点就足以满足实施的要求。但是，胸部放射影像是以光栅格式将器官和组织进行迭加而成，而其本身并不提供以向量格式存在的器官或组织的形状和位置。本发明利用一种智能复合方法，将胸部影像的脊柱形状进行向量化，该智能复合方法基于预先设定的标准便于将相关方法转变为模糊估计。在介绍本方法之前，需要先说明一下表1中的符号。
表1符号的定义
图2为智能复合方法的流程图。该方法处理了胸部影像中除底部边缘的每一条线，在此将底部边缘定义为位于图像底部的、高度为总高度十分之一的部分。流程图的第二个判定框将该底部边缘纳入该方法中。该边缘包括遮盖肝脏的脊柱区，是否属于脊柱侧凸并不根据该部分脊柱而确定。在第一个判定框中，通过判断三个条件以确定来自相关方法或�：兰频氖涑鍪菽囊桓鲇Ω檬亲钪盏墓兰萍怪行牡�。第一个条件是用来确认相关方法的可靠性，即Cj＞0.94。第二个条件是用于确保脊柱的连续性，即|Xcj-Xj-1|＜0.0284*W。第三个条件是j＞0.341*H，当评价上纵膈处脊柱时，此时不符合该条件，则使用�：兰疲蛭：兰圃诟们蛑杏τ眯Ч芎谩６逼兰巯伦蓦醮怪保耸碧跫闪ⅲ捎谛脑嘁跤�、大动脉和脊柱的叠加而需要更复杂的�：兰疲壹扑闶奔湟哺ぃ耸备屎侠孟喙胤椒ā�0.94、0.0284和0.341是在本方法的实验测试中通过一定数量的细算得出的参数。
(i)相关方法由样品图像获得下纵膈处的脊柱水平剖面模板。该模板只是一系列n像素数值，PT＝[P(1)，P(2)，...，P(n)]，如图3所示。
设图像第j行和第i 列的像素值为P(i，j)。对下纵膈处的每一条线(如j＞0.341*H)，考虑连续像素窗口，即，Pw(x，j)＝[P(x+1，j)，P(x+2，j)，...，P(x+n，j)] (1)其中，x＝0，1，...，W-n。相关系数C(x，j)被定义为 其中，cov(a，b)为向量a和b的协方差，var(a)是向量a的方差。第j行的C(x，j)最大值为Cj＝max C(x，j)，x＝0，1，...，W-n (3)因此，估计的脊柱中心点为xcj＝xmax+n/2(4)其中xmax是所有x的平均值，因而C(x，j)＝Cj。这样，该相关方法对于图像的每条线得出两个数值(xcj，Cj)，如表1所示。
(ii)�：兰撇煌钠鞴�、组织或物质依据其各自密度和厚度对X射线曝光具有不同程度的衰减。骨骼、软组织和空气应该由放射影像中不同范围的像素值来表示。我们可以把像素数据看作只表示一个单色图像平面，其中样本值意于在进行兴趣体积区(VOI)灰度转化后显为黑色，即DICOM标准中的光度解释等于‘单色2’。骨骼、软组织和空气应分别用高、中、低范围的灰度像素值表示。然而，找到脊柱的中心点并不是如此直接，因为该点不一定在高范围内具有最大像素值。因为神经�：缡且恢种悄苣Ｐ停浒ǘㄎ�(localized)基函数和训练权重[1，2]，脊柱图像的图案可以很容易地被融入该模型中。可以将基于已经广泛用于系统识别的神经�：缃峁沟哪：兰疲糜谘罢页晌怪行牡愕南袼氐乃迫唬缓蠡谡庑┧迫唤泄兰�。
传统神经�：绲慕峁瓜允居谕�4中。假设神经�：缇哂械ヒ皇涑觯淙肟占溆烧じ衩枋觯闹行脑谡じ竦母鞯闵�。B-样条(B-Spline)函数通常作为神经�：绲幕齕1-6]。
神经�：绲氖涑黾扑闳缦聐(t)=Q(x(t))=Σi=1nθi'μi(x(t))----(5)]]>其中，μi和θ’i分别为第i基函数和第i权重。数据可通过N对输入-输出对表示，如[x(1)，y(1)]，[x(2)，y(2)]，...，[x(N)，y(N)]。输入-输出的映射关系表示为y＝Q(x)。对于m个多输出，映射关系为Qk(x)=Σi=1nθik'μik(x)]]>其中k＝1，2，...，m (6)为了估计脊柱的中心点，将输入x定义为一个变量，其代表沿胸部影像上的一条线的水平位置，输出Qk(x)表示作为基于输入x的脊柱中心点的水平位置(k+n/2)的似然。假定n为偶数，m＝W-n。权重θ’ik为坐标(k-n/2+i-1，j)上的像素值，即P(k-n/2+i-1，j)，基函数μik是一阶B-样条(B-Spline)函数 图5为基函数的图示。
从图5可以看出，这种基函数具有一个重要特性，∫0Wμik(x)dx=1]]>i，k (8)通过沿着图像第j条线对Qk(x)进行积分，并将(8)代入积分中，(k-n/2+i-1，j)的总似然Lkj为，
Lkj=∫0WΣi=1nP(k-n2+i-1,j)μik(x)dx]]>=Σi=1nP(k-n2+i-1,j)∫0Wμik(x)dx]]>=Σi=1nP(k-n2+i-1,j)]]>其中k＝1，2，..，W-n注意，根据待检测特征的稀少度，可以利用高阶B-样条(B-Spline)函数。然而，如果采取高阶B-样条(B-Spline)函数，则总似然将更复杂，因此，计算时间更长。由于脊柱的像素值在灰度的高范围内，�：兰频氖涑鑫轿恢胟的平均值，其中总似然为最大值，即xjf=ave(argmaxkLjk)----(10)]]>注意，对于选择总似然的算子，如我们这里选用的‘arg max’，依我们想检测的特征而定，也可以随不同目标对象而改变。
因此，表示脊柱第j条线中心点的向量(u(j)，v(j))为(xfj，j)或(xcj，j)，后者视相关方法得到的Cj值而定。
(b)代价函数的计算脊柱侧凸的发病可在获得胸部影像的每条线的向量(u(j)，v(j))之后而确定。通过归谬法进行确定，其中假设脊柱是正常脊柱，脊柱为近乎一条直线。在这种假设下，向量(u(j)，v(j))能够符合下式的直线u＝m·v+c (11)下式给出零模型误差e(j)＝u(j)-m·v(j)-c≈0 (12)然而，数据不可避免的包括一定程度的噪音量，因此，这里采用了最小二乘法[7]以获得最吻合的直线。设胸部影像有N条线。根据式(12)，模型误差为
ϵ=e(1)e(2)...e(N)=u(1)u(2)...u(N)-v(1)1v(2)1......v(N)1mc=U-Vw----(13)]]>其中，ε＝[e(1)，e(2)，...，e(N)]T，U＝[u(1)，u(2)，...，u(N)]T，V＝[v(1)，v(2)，...，v(N)]T，且w＝[m，c]T。根据最小二乘法，w的最佳估计值为w^=(VTV)-1VTV----(14)]]>从(13)和(14)可以得出，相应模型误差为ϵ^=U-Vw^----(15)]]>为了评价数据吻合直线的程度，采用均方误差(M.S.E)M.S.E=1NΣj=1Ne^(j)2=1Nϵ^Tϵ^----(16)]]>因此，通过代价函数将脊柱侧凸指数量化为均方误差。
M.S.E.的值可被分成S和R两个部分，即M.S.E.＝S+R。R表示脊柱侧凸和噪音的可接受程度，S表示脊柱侧凸的显著程度。R值即误差水平是由专家确定的。R值的确定方法如下利用交互式用户界面来接收放射学家输入的数据。本智能系统此时处于训练模式，系统将呈现胸部放射影像集合及其相应的计算后的误差水平给放射学专家。放射学专家则评价该病例是否属于脊柱侧凸，因此在确定R值中放射学专家用其评价来训练该系统。放射学家利用现有的胸部放射影像的训练子集确定R值。例如，如果一名患者并不患有脊柱侧凸，从系统得到的R值为70，放射学家在用户界面中选中“no”，则R＝70将被认为是正常的脊柱。或者，如果一张胸片显示脊柱侧凸，此时R＝90，专家确定存在脊柱侧凸，并选中“yes”，因此，当R＝90时，胸片被归类为脊柱侧凸。然后，系统根据一系列病例计算R值的平均值和标准偏差，从而获得R值，此时R值为平均值。阈值R，即平均值R就是指示器，只要大于该值的图像，被分类为脊柱侧凸。
如果M.S.E.＜R，假设有效，则脊柱正常。如果M.S.E.＞R，脊柱不接近直线，根据归谬法，脊柱则异常。本智能系统的整个运行如图6所示。
本发明的智能系统的特征有1.早期检测特发脊柱侧凸2.监测脊柱侧凸的发展状况和侧凸程度3.基于X-射线胸部影像，检测脊柱侧凸的严重程度，无需进一步对整个脊柱进行X-照射，减少了辐射剂量4.没有人为干预，从图像数据库中自动提取可能的脊柱侧凸图像本发明智能系统的应用范围应用本系统是为了形成一种图像内容索引方法以从PACS中检索图像，从而实现1.开发基于内容智能图像检索方法2.基于图像特征提供参照例以助临床诊断3.为图像分析的图像数据研究提供一种方法应用实施例从香港理工大学的临床放射学中的PACS患者数据库中选出100例病例，其中一半为正常，另一半为脊柱侧凸。智能增强系统能准确地自动检测出具有脊柱侧凸的病例。结果与Cobb’s角相关，Cobb’s角是一种表示脊柱侧凸弯曲度的指标，其通过直接测量放射影像胸片而得，正如解剖学家所测定的。
智能增强系统能够从PACS图像数据库中直接读出DICOM图像。胸部影像可被传送到应用服务器上，通过我们的系统进行筛选，或者我们的系统能直接读出图像服务器内的胸部影像，从而检测出脊柱侧凸。
在香港理工大学的装置中，胸部影像被传送到智能增强系统中的应用服务器上，从而检测胸片是否为脊柱侧凸。
权利要求
1.一种早期检测特发性脊柱侧凸的智能增强系统，其特征在于，所述系统包括a.图像处理�？椋靡曰竦眯夭坑跋窆庹な荩籦.向量化模块，用以将图像中的脊柱向量化，得到脊柱估计中心点；c.代价函数计算模块，用以计算图像的脊柱向量吻合直线的程度，得到均方误差；d.判定�？椋靡耘卸暇轿蟛钣脬兄礡的大�。佣枷穹掷辔＜怪图怪嗤�。
2.根据权利要求1所述的智能增强系统，其特征在于，所述向量化�？榘ㄏ喙胤椒ǖピ⒛：兰频ピ把≡竦ピ鲅≡竦ピ挥谙喙胤椒ǖピ湍：兰频ピ埃靡匝≡窠邢喙胤椒ɑ故悄：兰�。
3.一种利用权利要求1或2所述的智能增强系统进行早期检测特发性脊柱侧凸的方法，该方法包括如下步骤(1)处理图像，获得胸部影像的光栅数据；(2)将图像的脊柱向量化，得到脊柱估计中心点；(3)计算代价函数，得到均方误差；(4)将图像分类通过比较均方误差与阈值R的大小而将图像分为正常脊柱和脊柱侧凸。
4.根据权利要求3所述的方法，其特征在于，步骤(2)的向量化是采用相关方法和�：兰疲⑶以谙喙胤椒ê湍：兰浦笆茄≡癫街瑁貌街栌靡耘卸鲜遣捎孟喙胤椒ɑ故悄：兰�。
全文摘要
本发明涉及一种利用胸部放射成像早期检测特发性脊柱侧凸的智能增强系统，该智能增强系统包括a.图像处理�？椋靡曰竦眯夭坑跋窆庹な荩籦.向量化�？椋靡越枷裰械募怪蛄炕玫郊怪兰浦行牡悖籧.代价函数计算�？椋靡约扑阃枷竦募怪蛄课呛现毕叩某潭龋玫骄轿蟛睿籨.判定�？椋靡耘卸暇轿蟛钣脬兄礡的大�。佣枷穹掷辔＜怪图怪嗤�。本发明还涉及利用上述智能增强系统进行早期检测特发性脊柱侧凸的方法。
文档编号G01T1/00GK1795823SQ20041010209
公开日2006年7月5日 申请日期2004年12月23日 优先权日2004年12月23日
发明者邓福禧, 陈颖志 申请人:香港理工大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
