CN102549618B - 用于从图像分析组织的方法和系统 - Google Patents
用于从图像分析组织的方法和系统 Download PDFInfo
- Publication number
- CN102549618B CN102549618B CN201080044509.6A CN201080044509A CN102549618B CN 102549618 B CN102549618 B CN 102549618B CN 201080044509 A CN201080044509 A CN 201080044509A CN 102549618 B CN102549618 B CN 102549618B
- Authority
- CN
- China
- Prior art keywords
- breast
- tissue
- image
- reference point
- edge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/44—Analysis of texture based on statistical description of texture using image operators, e.g. filters, edge density metrics or local histograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30068—Mammography; Breast
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
公开了一种从图像分析组织的方法,其包括:提供组织(100,400,450,600,800,1100)的电子图像,从图像确定参考值(1070,1170,1270),建立图像的hint表示(500,700),并且在组织的分析中使用hint表示,以量化乳房以及计算校准误差。还公开了一种系统,其在电子图像上运行内部乳房边缘检测算法(1310)以检测所述图像上的内部乳房边缘(1315),并且如果校准误差是不可接受(1324)则改进内部乳房边缘位置(1340)。还公开了乳房成分的自动估计以及图像的时间分析。
Description
发明领域
本发明涉及图像的分析。建立了图像上的参考值以及使图像可以自校准并因此量化。
发明背景
本领域的很多人员为了不同的目的而试图量化图像,但是在成像物理数据中遇到与误差或者未知值有关的问题。例如,试图基于绝对物理模型连同组织的假设属性来量化图像需要所有的成像物理数据是已知并且精确的。甚至量化梯形楔块包括在每一个图像中仍然需要大量的成像物理数据并且假定它们是精确的。因此,需要使用图像内参考值。
例如,通过向乳房发送X射线光子并然后检测多少X射线光子通过它来创建乳房X射线照片。通过的X射线光子的数量越小,乳房组织就越密
乳房密度通过很多研究与发展乳癌的概率相关联,并且这些研究中的大部分通过使用视觉或者半自动式方法来评估乳房密度。一般综述例如在乳癌研究的评论丛书中给出,包括Vachon等人的“Mammographicdensity,breastcancerriskandriskprediction”(乳腺X射线摄像密度,乳腺癌的风险以及风险预测)(BreastCancerResearch,2007,vol9:217),Martin和Boyd的“Potentialmechanismsofbreastcancerriskassociatedwithmammographicdensity”(与乳腺X射线摄像密度有关联的乳腺癌风险的潜在机制)(BreastCancerResearch,2008,VoI10:201)以及Yaffe的“Measurementofmammographicdensity”(乳腺X射线摄像密度的测量)(BreastCancerResearch,2008,VoI10:209)。
现有技术
特别是,在这里参考:Highnam(1992)“Model-basedmammographicimageprocessing”(基于乳腺X射线摄像图像处理的模型),(PhDThesis,UniversityofOxford1992),Highnam和Brady1999“MammographicImageAnalysis”(KluwerAcademicPublishers,1999),Highnam2006“Breastcompositionmeasurementsusingretrospectivestandardmammogramform(SMF)”(使用过去标准乳腺X射线照片形式的乳腺构图测量)(PhysicsinMedicine&Biology,2006),Highnam2007“Comparingmeasurementsofbreastdensity”(对比乳腺密度的测量)(PhysicsinMedicine&Biology,2007),VanEngeland2006“Volumetricbreastdensityestimation”(乳腺体积密度评估)(IEEEMedicalImaging,2006),Hartman2008“Volumetricassessmentofbreasttissuecomposition”(乳腺组织组成的体积核计)(InternationalWorkshoponDigitalMammography2008),Kaufhold2002“Calibrationapproachtoglandulartissuecomposition”(腺性组织构图的校准方法)(MedicalPhysics,2002),Blot2005“Volumetricapproachtoriskassessment”(风险核计的体积的途径)(PhysicsinMedicine&Biology,2005),Heine2006“Effectivex-rayattenuationmeasures”(有效地X射线衰减测量)(MedicalPhysics,VoI33,2006),Heine2008“Effectivex-rayattenuationcoefficientmeasurementsfromtwofull-fielddigitalmammographysystems”(两个全领域的数字化乳腺X射线相片的有效地X射线衰减系数测量)(BioMedicalEngineeringOnLine,2008,VoI7:13)
一般量化肌理并将它们应用于乳房X线照片的标准的图像处理技术的使用还在下列文献中被描述:Miller1991,“Classificationofbreasttissuebytextureanalysis”(肌理分析的乳腺组织的分类)(BritishMachineVisionAssociationConference,1991),Karssemeijer1998,“Automatedclassificationofparenchymalpatternsinmammograms”(乳腺X射线照片中实质模式的自动分类)(PhysicsinMedicine&Biology,1998),Olsen2007“Automaticsegmentationoffibroglandulartissue”(腺性纤维组织的自动切分)(ImageAnalysis2007),Raundahl2008,“Automatedeffect-specificmammographicpatternmeasures”(自动乳腺X射线摄像模式测量的具体效果)(IEEEMedicalImaging2008),Hadley2007“Riskclassificationofmammogramsusinganatomicallinearstructuresanddensity”(使用解剖的乳腺X射线照片的线性结构和密度的风险分类)(PatternRecognitionandImageAnalysis,2007),Freixent2008“Eigendetectionofmassesconsideringfalsepositivereductionandbreastdensityinformation”(大众考虑的假阳性衰减以及乳腺密度信息的特征检测)(MedicalPhysics,2008),Raundahl2007“Breasttissuedensitymeasure”(乳腺组织密度测量:患上)(WO/2007/090892),Magnin1986“Mammographictextureanalysis:anevaluationofriskfordevelopingbreastcancer”(乳腺X射线摄像肌理分析:患上乳腺癌的风险的评估)(OpticalEngineering1986),Caldwell1990。
此外:“Characterizationofmammographicparenchymalpatternbyfractaldimension”(通过分型维数的乳腺X射线摄像的实质模式的特征化)(VoI35,PhysicsinMedicineandBiology1990),Luo2007“Methodforclassifyingbreasttissuedensity”(乳腺组织密度的分类方法)(WO/2008/088478),Giger2003“Risk-modulateddiagnosisofdisease”(疾病的风险收调诊断)(WO/2003/067371),Hui1999“Computerizedassessmentofbreastcancerrisk”(乳腺癌风险的计算机化评估)(WO/1999/063480),Thornton1998“Earlydetectionandaidtodiagnosis”(早期发现及诊断援助)(WO/1999/008225),Li2006“Mammaryglandquantification”(乳腺量化)(CN01234026),Zhou2001“Estimationofbreastdensityonmammograms”(在乳腺X射线照片上的乳腺密度的评估)(MedicalPhysics28(6),2001),Saha2001“Breasttissuedensityquantificationviadigitizedmammograms”(通过数字化的乳腺X射线照片的乳腺密度量化)(IEEEMedicalImaging,VoI20,8,2001),Bovis2002“Classificationofmammographicbreastdensityusingacombinedclassifierparadigm”(使用综合分类范式的乳腺X射线摄像乳腺密度的分类)(IWDM2002,p77-80,2002),Petroudi2003“Classificationofmammographicparenchymalpatterns”(乳腺X射线摄像实质模式的分类)(IEEEMedicine&Biology,VoI2,p416-423,2003),Lui2004“Mammographicparenchymalpatternsforassessingrisk”(用于风险评估的乳腺X射线摄像实质模式)(MedicalPhysics,VoI31,no3,2004),Karssemeijer1998“Automatedclassificationofparenchymalpatternsinmammograms”(在乳腺X射线照片中实质模式的自动分类)(PhysicsinMedicine&Biology,1998),Sivaramakrishna2001“Automaticsegmentationofmammographicdensity”(乳腺X射线摄像的自动切分)(AcademicRadiology,),Wei2004“CorrelationbetweenmammographicdensityandvolumetricfibroglandulartissueestimatedonbreastMRimages”(在乳腺MR图像上的乳腺X射线摄像密度与腺性纤维组织评估之间的相关性)(MedicalPhysics,31(4),2004),Boehm2008“Automatedclassificationofbreastparenchymaldensity:topologicanalysisofx-rayattenuationpatterns”(乳腺实质密度的自动分类:X射线衰减模式的拓扑分析)(AJR,191,p275-282,2008),Jamal2006“BreastdensityintoTabar'spatterns”(塔巴尔模式中的乳腺密度)(PhysicsinMedicine&Biology,2006),Torres-Majia2005“ComparisonofQualitativeandQuantitativeEvaluations”(定性与定量的评估的比较)(CancerEpidemiologyBiomarkers&Prevention,2005),Heine2008。
其他来源包括:“Automatedapproachforestimationofbreastdensity”(用于乳腺密度的评估的自动方法)(CancerEpidemiology,Biomarkers&Prevention,VoI17,Nov2008),Oliver2008“Novelbreasttissuedensityclassificationmethodology”(创新式的乳腺组织密度分类方法)(IEEEITinBiomedicine,2008),Li2005“Computerizedtextureanalysisofmammographicparenchymalpatterns”(乳腺X射线摄像实质模式的计算机化肌理分析)(AcademicRadiology2005),Megalooikionomou2007“Analysisoftextureinmedicalimageswithanapplicationtobreastimaging”(使用乳腺成像应用的医学图像中的肌理分析)(SPIEMedicalImaging2007),Chang2002“Computerizedassessmentoftissuecomposition”(组织构成的计算机化评估)(AcademicRadiology,2007),Glide-Hurst2007“Anewmethodforquantitativeanalysisofmammographicdensity”(用于乳腺X射线摄像密度量化分析的一种新方法)(MedicalPhysics),Taylor1994“Measuringimagetexturetoseparatedifficultfromeasymammograms”(检测图像肌理从简单乳腺X射线照片中分离困难乳腺X射线照片),以及Heine2000“Statisticalmethodologyformammographicdensitydetection”(用于乳腺X射线摄像密度检测的统计方法)(MedicalPhysics,27,2000)。
这些参考文献还指出了尝试仿真半自动化阈值技术的已知手段,例如积云,其被理解为对乳房密度测量是典范的。
乳房组织的量化的一个解决方案需要梯形楔块插入每个图像中。例如:Diffey2006“Newstepwedgeforvolumetricmeasurementofdensity”(用于密度的体积测量的新式步进式光楔)(IWDM2006),Augustine2006“Volumetricbreastdensityestimation”(IWDM2006),Patel2006“Automatedbreasttissuemeasurementofwomen”(女性的自动乳腺组织检测)(IWDM2006)以及Shepherd2005“Noveluseofsinglex-rayabsorptiometryformeasuringbreastdensity”(为测量乳腺密度的单X射线吸收仪的创新式使用)(TechnologyinCancerResearch&Treatment,2005,VoI4:173-182)。
在本领域中已知的可选途径包括:使用简单的物理过程将乳房分为密度类,即使这样的方法并不表现为使用图像内参考数据(Merelmeier2007“Displayinganx-rayimagerecordedonmammography”(展示记录在乳腺X射线相片上的X射线图像),WO/2008/052854);不使用任何图像校准的简单组织密度方法(Tasaki2007“Imagingsystem”,WO/2008/038525);通过侦察X射线曝光的组织密度(Jing2003“Mammographywithtissueexposurecontrol,Tomosynthesis”(组织曝光控制的乳腺X射线相片,层析X射线照相组合),WO/2004/049949);使用在压缩板顶部上的栅格以找到乳房厚度(Yang2003“Measuringthethicknessofcompressedobjects”(测量压迫物件的厚度),WO/2004/000121);使用幻影图像计算出乳房边缘厚度(Rico2003“Determiningperipheralbreastthickness”(确定周围的乳腺厚度),WO/2004/000110);用于密度的人工评估的物理栅格的使用(Bershtejn2008“Adiposetissuepercentageevaluationinmammarygland”(乳腺中脂肪组织百分比的估算),RU2325852);用于预后目的的密度的使用(Giger2003“Computerisedimageanalysisprognosis”(计算机化图像的预后分析),US2004101181);用于重建目的的密度图像的使用(Clause2008“Quantitativeimagereconstructionmethod”(量化图像的分析方法),US20080292217);涵盖胶片曲线的胶片工作(Highnam2000“X-rayimageprocessing”(X射线图像过程),WO/2000/052641);使用胸肌以标准化数据(Kotsuma2008“QuantitativeassessmentofmammographicdensityandbreastcancerriskforJapanesewomen”(对于日本女性的乳腺X射线摄像以及乳腺癌风险的量化核计),Breast,2008);以及改进的半自动阈值技术(Byng1994“Quantitativeanalysisofmammographicdensities”(乳腺X射线摄像密度的量化分析),PhysicsinMedicine&Biology,1994);依照骨骼密度使用双能量(Shepherd2002“Measurementofbreastdensitywithdualx-rayabsorptiometry:feasibility”(使用双X摄像吸光测量的乳腺密度的测量:可行性),Radiology2002);从MRI计算密度(Klife2004"QuantificationofbreasttissueindexfromMRdatausingfuzzyclustering"(使用模糊簇聚的核磁共振数据的乳腺组织系数的量化),IEEEEngineeringinMedicine&BiologySociety,2004);以及从超声波计算密度(Glide2007“Novelapproachtoevaluatingbreastdensityutilizingultrasoundtomography”(利用超声波层析X射线照相法得到的乳腺密度的估算的创新式方法),MedicalPhysics,2007)。
然而,图像物理中的误差以及未知值使得这样的方法固有地不可靠。
Highnam和Brady从事的研究导致对如何使用图像处理和X射线物理学与相关成像物理数据的混合从乳房X射线照片自动计算乳房组织的密度的理解,该相关成像物理数据通常涉及在乳房X射线照片中的坐标(x,y)处测量的像素值P(x,y)到脂肪厚度hfat(x,y)厘米以及“所关注的组织”的厚度hint(x,y)厘米的转换,其中“所关注的组织”可以是纤维组织、腺组织、水或者癌组织。
此外,近来对通过一种称为层析X射线照相组合的技术来创建乳房的三维图像有兴趣,根据该技术以在乳房周围的变化的角度拍摄一系列乳房X射线照片,并重建三维图像。
如Highnam和Brady(1999)提出的,到此为止Hint表示被计算的过程包括如在图2示出的8个步骤。相似的过程由VanEngeland2006(虽然不包括步骤S3,S6,S7,S8)和Hartman2008提出,Hartman注意到S4对数字乳房X射线相片并不是必要的,并且报告了对S6的修改。
Kaufhold(2002)以及Heine(2006,2008)以更系统的方式考虑了生成hint的过程,在数据收集之前在一个X射线系统上生成被校准的数据,而不使用用于自校准目的的图像数据。
更详细地,并且参考图2:来源于胶片乳房X射线照片202、数字乳房X射线照片204和/或数字断层投影206的数据须经以下列步骤开始的多个步骤:
步骤1210分割–某种形式的图像处理应用来描绘内部乳房边缘310以及外部乳房边缘360以及胸肌380。
步骤2215像素值至能量–将乳房X射线照片中的在位置(x,y)处的像素值P(x,y)转换至传给在(x,y)处的检测器的能量。这涉及到对特定检测器的物理过程的理解。
步骤3220检测器模糊校正–使用特定检测器的物理过程的知识来移除它可能引入的任何模糊。对于数字系统,这个步骤通常是不必要的。
步骤4225校正阳极足跟–考虑到在乳房上给出在空间上变化的入射辐射强度的X射线管的阳极足跟效应,并且校正它,以便给出出现得好像入射辐射强度是恒定的图像。在数字系统中,这通常由生产商自动执行作为他们“平场处理”过程的一部分。
步骤5230校正压缩板倾斜–确认当乳房在准备乳房X射线相片中被压迫时,顶部压缩板倾斜;该校正(一般)通过添加脂肪的虚拟楔块来补偿该倾斜。
步骤6235估计乳房厚度–认识到由X射线机记录的乳房厚度可以是不精确的,并且通过找到图像的大部分是脂肪的区域并且然后计算出必须有多少脂肪可以从图像本身进行另一个估计。
步骤7240移除分散分量;以及
步骤8245移除焦外分量。传给检测器的能量拥有主要分量(不受阻碍地直通过乳房的X射线光子)、分散分量(X射线光子被偏转)以及在乳房边缘周围与焦外辐射有关的分量(远离焦斑从X射线管出去的X射线光子)。在这些步骤中,估计分散分量以及焦外分量以仅仅获得主要分量。
步骤9250从主要分量生成hint255。从主要分量中,可以计算出必要的X射线衰减,并且因此以数学方式得出脂肪以及所关注的组织的厚度:
方程1:
其中φ是光子通量,Vt是X射线管电压,Ap是像素面积,ts是曝光时间,N0 rel是该管电压在能量ε处的入射X射线光子的相对数量,S(ε)是检测器吸收,G(ε)是电网传输,且最后的两项是归因于有机玻璃压缩板的衰减和归因于乳房的衰减。
从前面的步骤可知方程的左手侧,由于值与图像一起被自动存储,成像物理数据(X射线管电压等)对数字图像是已知的,因此唯一未知的是归因于乳房的衰减,hμ(x,y,ξ)。归因于乳房的衰减可使用迭代技术来确定,该迭代技术需要尝试各种所关注的组织的厚度以及脂肪的厚度,直到所需的衰减被看到:
hμ(x,y,ξ)=hint(x,y)μint(ξ)+hfat(x,y)μfat(ξ)
这也可使用线性近似来解决。显然,这有两个未知量,但是它们不是独立的:
hint(x,y)+hfat(x,y)=H(压缩的乳房厚度,从步骤6)。
从以上的方程应当清楚,所有的成像物理数据都很重要,但是特别是,压缩的乳房厚度H是关键参数,正如Highnam以及Brady1992、1999、2006等所证实的。有问题的是,胶片乳房X射线照片很少有所记录的乳房厚度,并且这就是为什么来自图像的估计通常对胶片图像作出。
存在所描述的方法可被使用的很多出版物:例如由Ding、McCormack和Jefferies,并且特别是对于胶片,该方法表现为适度地好得起作用。
就整个乳房的体积实际上是乳房内的像素数目乘以乳房厚度来说,乳房厚度特别是相关的:
乳房体积(厘米3)=像素数目×像素面积×乳房厚度
并且乳房密度是所关注的组织与乳房体积的比值。密度以及体积均实质上依赖于精确的乳房厚度。事实上,学术工作(Highnam与Brady1999,Blot2005)提出乳房厚度中的1mm的误差可导致乳房密度中的2-3%的误差。
最近,Hartman2008报告了在数字乳房X射线照片上使用基于Highnam与Brady方法的软件,由此Highnam与Brady方法被调整为结合所估计和记录的乳房厚度,而不是仅仅所估计的乳房厚度。
Highnam(2006)证实了乳房厚度的重要性,并且描述了乳房厚度的估计如何需要成像物理数据补偿的形式。也就是说,通过使用图像本身并通过调整乳房厚度来校正在曝光时间、管电压以及其他成像物理数据中的误差。这在能够找到表示完全多脂组织的图像的区域的基础上来完成。
虽然在对照视觉评估生成有效的hint值以及合理的乳房密度方面是成功的,该方法还是有一些缺陷,特别是对于所记录的乳房厚度从数字机器本身知道的数字图像。例如,所记录的乳房厚度可能不是主要误差,特别是在量控制的位置处,在该位置上,所记录的乳房厚度在定期的基础上被仔细地校准及检查。因此,通过乳房厚度的估计执行的误差校正可能不以正确的方式校正误差。如果误差处于曝光时间(ts),则误差校正应该是倍增的;相反,误差以指数方式在积分内被有效地校正。如果被认为是脂肪组织的图像的区域不是脂肪组织而是致密组织,则有潜在的进一步误差。
VanEngeland2006采取了这样的观点:所记录的乳房厚度是可靠的,并且简化了方程以得出像素值P与hint之间的关系。
方程2:
hint(x,y)=ln(P(x,y)/Pfat)/(μfat-μdense)
其中Pfat是一块组织的像素值,该组织在其内只有脂肪(而不是皮肤),并且μfat以及μdense是有效的X射线线性衰减系数,该衰减系数由管电压、靶、过滤器、以及所记录的乳房厚度计算出。
大体上,这样的区域可使用图像处理找到,但是对于致密的乳房,该图像处理常常失败。在图4中,左侧的乳房X射线照片是多脂乳房400,并且对确定全脂肪的区域相对简单。右侧的乳房是非常致密的450,并且找到这样的全脂肪的区域更困难。
本发明通过总是找到在图像中的可靠参考点并由此计算明显的校准误差来解决计算和计算底数中的错误的指示以及误差,这意味着,精确的成像物理数据是不必要的,并且事实上毫安(mAs)以及很多检测器信息可被忽略,而不是假设检测器随着已知偏移是线性的
本方法因此对成像数据中的误差以及未知数据更加强大,并且具有相关的校准误差因子,该校准误差因子可被用于在必要时警告用户。
找到好的参考点的以前尝试依赖于找到在内部乳房边缘内的抗噪声最大值或者最小值(以指示乳房的脂肪最多的部分)(见图3),其中内部乳房边缘310通过找到乳房边缘周围的最里面的均匀曲线来确定(HighnamandBrady1999)。
其他解决方案试图使用局部相位及方向来找到内部乳房边缘(“Theuseofmulti-scalemonogenicsignalonstructureorientationidentificationandsegmentation”,Pan,Brady,Highnam,Declerck,InternationalWorkshoponDigitalMammography2004);并且使用所记录的乳房厚度以及球形乳房边缘模型来从外部乳房边缘移动一定的距离(VanEngeland2006)。使用均匀曲线方法很慢;由于噪声问题,使用局部相位是不可靠的;并且虽然使用球形模型很快并对多脂乳房很好地起作用,然而是它对致密的乳房不能很好地起作用。
本发明以一种方式确定内部乳房边缘,该方式对于所有乳房都是可靠的,并且因而导致对所有图像的可靠的Pfat值,以及由此可计算出成像物理数据误差。
发明概述
根据本发明的第一个方面,提供了一种从图像分析组织的方法,该方法包括如下步骤:
提供组织的电子图像;
确定参考像素;
建立图像的hint表示;以及
在组织的分析中使用hint表示,包括提供校准误差。
根据本发明的第二个方面,提供了一种方法和系统,由此在图像数据中的坐标(x,y)处测量的像素值P(x,y)被转换成“组织A”(例如脂肪)的厚度hfat(x,y)厘米以及“组织B”的厚度hint(x,y)厘米,其中“组织B”包括“所关注的组织”,例如纤维组织、腺组织、水或者癌组织。可选地,图像数据可被转换成单个组织厚度。
本发明因此涉及对图像的分析,由此根据在内部乳房边缘的内部的乳房组织或者在内部乳房边缘和外部乳房边缘之间的乳房组织的分析来选择参考值,其中该内部乳腺边缘是使用相位一致性或者由乳房组织图像统计建立的阈值来最初确定的。一旦找到参考值,图像就可被自校准并因此量化。
本发明通过总是找到在图像中的可靠的参考点解决了在计算和计算底数中的错误指示和误差,该参考点随后允许明显校准误差的计算。
本方法因此对成像数据中的误差以及未知数据更加可靠,并且具有相关的校准误差因子,该校准误差因子用于在必要时警告用户。
本方法优选地包括用于建立参考点的步骤,其中,
生成图像中乳房组织的直方图;
检测位于直方图的两个峰值之间的间隙;以及
该间隙位置对应于内部乳房边缘,该内部乳房边缘用于建立搜索参考点的地方。通过使用相位一致性或者通过确定阈值也可以检测该间隙,该阈值通过查看乳房图像统计例如最小和最大像素值、平均像素值或者中间像素值来计算。
可选地,通过从乳房边缘中的数据外推参考像素值来建立参考点。
有时候,对于某些乳房来说,参考点并不理想,且可以通过将内部乳房边缘向内移动例如通过下列方式中的一个或多个来实现改进:基于乳房区域;基于乳房边缘区域;和/或采用迭代的方式,直至沿着所述内部乳房边缘的足够长度找到强度边缘,或者找到平滑的均匀边缘。
一旦得到好的参考像素,明显校准误差(κ)就可被引入方程1中。该误差可被调整,直至图像数据(在该方程左侧)匹配参考像素(如该方程右侧所示)的理论数据:
认为1的校准误差是示例性的;在1之上或者1之下的运动指示误差。系统参数可被设置在任意期望值,然而接近1的校准误差会提供最精确的数据。
根据进一步的方面,提供了用于从乳房的图像分析组织的系统,该系统包括:
(a)提供组织的电子图像;
(b)在电子图像上运行内部乳房边缘检测算法,以检测图像上的内部乳房边缘以产生结果;
(c)检查结果的校准误差是否是可接受的;
(d)根据校准误差的偏差,将结果通知或者不通知给用户。
因此,到现在为止讨论的在本发明的两个方面中的组织的分析通过计算hint与hfat的值、计算hint与hfat的值的和并随后计算乳房密度来导致对乳房成分的自动估计。
另一个分析选项包括:优化显示准备让放射科医生观看的图像;追踪对象的随着时间的过去的发展以确定它是否增长;计算机辅助检测,由此,计算机在图像的可疑区域做标记以供放射科医生考虑;以及大致分析多个图像(时间的、双能量、视图等)。
优选地,根据本发明的方法和系统,在一段时间内收集的多个图像被分析并被自动比较。优选地通过使标准化图像对齐来实现对图像的这个时间比较。优选地使用参考点例如乳头和胸壁来便于这个对齐。报告可被生成,或者如果检测到变化,该报告被标记给放射科医生。
在本方法和系统的优选实施方式中,组织的钙化被量化并且从随后的分析中移除。这很有用,因为钙化是致密组织,且它的存在使得组织看起来更致密,因而影响任何结果。
优选地,在本方法和系统中,一系列参考点或者参考值被找到,以能够从不同成像模式区分多于两种类型的组织。
在所描述的方法和系统中,电子图像的像素值是立体像素值,该立体像素值相应于空间中的三维位置(x,y,z)。
该描述参照了一个实施方式:乳房图像且特别是乳房X射线照片(乳房X射线)的量化。
密度与乳房组织的成分紧密相关:多脂组织的密度是最小的;纤维组织以及癌组织是适度的密度;钙化是非常致密的。为了总的乳房成分的目的,钙化一般可被忽略,因为它以这样微小的量出现。然而,钙化的量化是我们的发明的一个重要的副产品,并且会在下面被提到。
令:
乳房中的多脂组织的体积=Vfat厘米3
乳房中的所关注的(无脂肪)组织=Vint厘米3
以及:
乳房密度=100*Vint/(Vint+Vfat)
本发明的一个实施方式涉及从图像例如乳房X射线自动估计乳房成分(体积以及乳房密度)的方法和系统。为了解释的清楚,本说明书将聚焦于X射线成像,虽然也可以对乳房MRI以及其它方式进行乳房成分估计(其中乳房密度是1)。
在优选的实施方式中,图像包括了X射线、MRI以及其它方式。
此外,本发明的特性在分析通过X射线、MRI以及其他方式生成的图像——包括随着时间过去的乳房X射线相片、双能量成像以及层析X射线照相组合以及其他——中被使用。理想地,使用序列中的每个乳房X射线照片的hint表示来执行多个图像到层析X射线照相组合图像的三维重建。
在另一实施方式中,本文描述的方法和系统反复应用于在一段时间内收集的图像,因此实现自动时间比较,从而实现预测分析。例如,癌组织的概率及发展通过原始数据的时间交叉引用以及优化监控来指示。
术语“乳房X射线照片”在本文件中用于指胶片乳房X射线照片、数字乳房X射线照片或者数字层析X射线相片投影之一。
附图简述
图1示出了使用正好在乳房上方并向下发送X射线的X射线源拍摄的乳房X射线照片(左侧);在技术上,头尾向的乳房X射线照片以及在处理右侧的hint曲线之后;
图2示出了用于计算所关注的组织的现有技术过程;
图3示出了将图像分成与组织类型有关的不同的部分的现有技术分割步骤;
图4示出了具有不同成分的不同乳房的两张乳房X射线照片;
图5示出了图6的乳房组织的例证性直方图;
图6示出了乳房X射线照片,该照片示出了具有多脂边缘部分的致密乳房;
图7示出了另一个乳房边缘;
图8示出了图7的乳房边缘的直方图;
图9示出了另外的乳房边缘;
图10示出了典型乳房的图像,该图像示出了参考区域;
图11示出了具有参考组织的区域的图像,该参考组织进入乳房边缘很远;
图12示出了不是脂肪的参考组织的区域的图像;以及
图13示出了根据本发明的系统。
所示实施方式的详细描述
在图1中,乳房X射线照片示出了相对多脂的乳房100(如图像的相对黑暗所指示的),一些更致密的组织明显地朝向乳头110(较白部分)。在称为hint表示150(在图片右侧)的部分中,朝着乳头160的峰值指示hint的高值,即,较致密的组织。
一旦生成了hint表示,许多临床上有用的任务就可以被执行了,包括但不限于:通过对hint以及hfat的值求和并随后计算乳房密度来自动估计乳房成分;优化准备让放射科医生查看的图像的优化显示;追踪对象随着时间的过去的发展以确定它是否增长;计算机辅助检测,由此,计算机在图像的可疑区域做标记以供放射科医生考虑;大致分析多个图像(时间的、双能量、视图等)。
第一个实施方式指出,致密乳房600具有多脂乳房边缘610、510连同致密内部620、520,例如图6所示的图像,其例证性直方图500在图5中提供。内部乳房边缘检测算法通过例如使用相位一致性或者从在乳房组织上计算的图像统计获得阈值设置来检测乳房直方图的两个峰值之间的间隙530。这样的技术实现起来快速且容易,并且为图6中示出的致密乳房提供示例性内部乳房边缘640。
如图6中所示的,所检测的内部乳房边缘640全部位于多脂乳房组织上(在多脂乳房边缘610中);它并不覆盖致密组织620并且因此产生高乳房密度,而不是错误的低乳房密度。如现在确定的内部乳房边缘640对致密乳房很好地起作用。然而,对于多脂乳房800它是次最优的,该多脂乳房800具有更相似于图7所示的直方图的直方图700,图7提供了乳房边缘710,如图8所示的810。
内部乳房边缘检测算法检测乳房直方图的两个峰值之间的间隙730,然而,乳房边缘太向外,并且由此生成过高的乳房密度。根据本发明,如果乳房是非致密的,即,由多脂组织850组成,上面确定的内部乳房边缘740、840向内移动。
这可以用很多方式实现。首先,当投影在图像上时,多脂乳房倾向于比致密乳房更大,所以基于乳房区域,内部乳房边缘可被带向内。可选地,内部乳房边缘可被带向内,直至你到达强度边缘或者直至你达到某一平滑度。
通过经由内部乳房边缘确定参考值,可靠的参考值被实现并且减小或者消除对一些成像物理因素的依赖。
以下结果是对单个图像在不同曝光时间下实现的:
并且通过改变在不同乳房密度处的四张图像的kVp:
下面表示从在具有四个不同乳房密度的乳房图像上运行具有乳房厚度(H)中的误差的软件得到的结果。
标准值 | 乳房厚度(H)的变化 |
Hmm | kVp | -20% | -10% | 0% | 10% | 20% |
86 | 30 | 2.5% | 2.4% | 2.3% | 2.2% | 2.1% |
32 | 25 | 15.1% | 13.6% | 12.4% | 11.4% | 10.6% |
20 | 25 | 23.7% | 21.4% | 19.6% | 18.0% | 16.8% |
40 | 28 | 33.6% | 30.5% | 28.1% | 26.2% | 24.4% |
找到好的参考值的另一方法是使用乳房边缘中的数据来外推参考像素值。数字乳房X射线相片具有宽的动态范围,这意味着乳房边缘信息相对精确地被捕获,即使它在没有周围的增强的情况下没有很好地显示。在外部和内部乳房边缘是已知(例如,从图像处理,或者从假设半圆形乳房边缘)并且所记录的乳房厚度是已知的场合,图9中的线A与线B之间的组织可被识别为皮肤和纯脂肪(乳房边缘按照定义几乎是脂肪)。如果假设在点A只有小量的脂肪,并且如果B在哪里是已知,并且已知半圆形乳房边缘模型是良好的估计,然后在点B处的脂肪量的估计可被进行,并且因此参考像素值可被计算出以用于直接从图像确认参考点值,或者可被替代地使用。
可选地,代替校准多脂组织,以上方法或者在如上面引用的Highnam中公开的方法可被用于确定它是致密乳房并且因此新的参考点被需要。因此获取主要乳房组织内的最大X射线衰减,并且假设皮下脂肪层顶部和底部都是致密组织。由此,可计算出量化表示并且接着向前到在每个像素处的值。
通过上面途径的任一个找到参考脂肪值Pfat,看来似乎可靠地提供脂肪参考值。从这个已知的值中,方程1以及方程2可被调整为包括显式线性置信度量(κ):
方程3:
κ被调整直至以上方程的左手侧以及右手侧在脂肪参考值处匹配。
因此,如图10所示的,在内部乳房边缘1040与外部乳房边缘1060被标记的场合,如圆圈1070指示的看来大部分是脂肪的区域指示具有典型的压缩的乳房厚度:即,远离在边缘中的乳房滚降。
作为另一例子,图11示出了太向外的参考组织1170的区域。因此,如果大部分是脂肪的Hcm的区域太向外而进入内部乳房边缘1140中并且事实上小于脂肪的Hcm,那么过高并且κ>1.0。
图12示出了不是脂肪的参考组织1270的区域。因此如果是大部分是脂肪的Hcm的区域相反是致密组织的区域,那么:
过低并且K<1.0。
其中有许多的校准误差,包括可能在所记录的乳房厚度中的校准误差,接着κ从1.0被向上或者向下调整。
因此,如果所有的校准数据接近正确,并且找到了图像合适的部分作为大部分是脂肪的区域,κ会接近1.0。所以,远离1.0的κ意味着在模型中有误差,并且需要校正行为(虽然很难知道误差是从哪里出现的)。总之,找到了一种识别误差以及生成显式置信测量的方法。
另外,可在软件中提供切断,使得一旦获得例如在0.95和1.05之间的某一κ值,数据就被认为是可接受的并且结果被传递给用户。
图13是根据本发明的系统的流程图。提供了电子图像以及在图像数据上运行内部乳房边缘检测算法(1300)。当产生结果(1315)时,查询询问κ值或者结果的校准误差是否是可接受的(1320);如果是(1322),那么组织在结果的基础上被分析(1330);如果否(1324),那么用户被告知有校准问题[改变图片]。
如果已知多脂乳房被检查,并且发现κ距离1.0很远,很明显使用X射线机器校准设置有重大误差。由此,用户可被警告以执行质量控制或者系统可开始生成校正参数。
参考乳房图像,随着体积更加可靠并且没有与乳房压缩有关的问题,比较可优先于乳房密度而基于区域中纤维腺组织的体积。结果被显示为减影图像,放射科医生可从经验中更容易理解该减影图像。
可通过例如在医院中的服务器来便于图像的比较。可选地,图像可以全部或者部分地被存储,但是在对以前运行的DICOM结构报告(其存储在PACS系统中)内足以使以前的信息能够在区域基础上被比较。该历史信息可能是加密的,使得用户无法看到它。由于MLO图像还包括复杂化比较的胸肌,可仅对CC进行时间比较。
另一实施方式包括通过“云”或者特别的方式的图像的存储以及时间比较,由此对可配置的计算资源(例如,网络、服务器、存储器、应用以及服务)的共享池的方便的点播网络访问可被迅速提供并以最小的管理努力或者服务提供商互动而释放。此外,部署为“作为服务的软件”实现容易的更新,使得所有图像使用相同的算法来处理。
在优选的实施方式中,乳房X射线照片的时间比较通过标准图像的对齐而被影响,该标准图像也许使用参考点例如乳头和胸壁如上所述而产生,从参考点创建报告内的图像,该图像展示在一段时间内在乳房中更多的组织在哪里发展。如果更多的组织发展(或者组织变得“更致密”),那么癌症可能会生长,并且放射科医生被警告更细致地检查此。
本发明进一步的优势在于用户能够上传原始数据用于分析和存储,以为比较时间数据作准备,避免对用户例如医院的需求,以储存原始数据。中等体积随着时间的过去用作规范器,由此不是匹配hint,而是使用hint/中等hint,移除了影响,例如更年期问题。
进一步的优势是,乳房的成功量化通过例如钙化的自动检测来考虑钙化的量化、通过我们发明的背景元素的相减、以及接着其余信号的量化。
应当认识到,这些图示只是用于说明目的,并且其他配置是可能的。
本发明通过具有修改和变更的一些实施方式被描述,但是在阅读并且理解了本描述后,另外的实施方式以及修改对于本领域的技术人员是明显的。所有这样的实施方式以及修改被规定为落在由所附权利要求限定的本发明的范围内。
Claims (11)
1.一种从乳房图像分析组织的方法,包括以下步骤:
提供乳房组织(100,400,450,600,800,1100)的电子乳房图像;
根据所述电子乳房图像中的像素值生成表示所述乳房图像的直方图;
分析所述直方图以确定称为内部乳房边缘(310,640)的所述乳房组织的边缘;
在所述乳房组织中,使用所述内部乳房边缘来确定所述乳房图像上的参考点(1070,1170,1270);
建立纤维组织、腺组织或者癌组织的图像的厚度hint表示;
确定所述hint表示的校准误差,且如果所述校准误差高于或低于示例性值(1324),则重新定位所述参考点并计算新的校准误差;
重复重新定位所述参考点并计算新的校准误差的步骤,直到获得可接受的校准误差(1322);
其中,通过以迭代方式向内移动所述内部乳房边缘的图像的边界并计算置信度量κ,直到获得所述置信度量κ的预定值,由此来改善所述内部乳房边缘的位置,所述置信度量κ是根据以下方程来定义和计算的:
其中:x是所述电子乳房图像中的像素的x坐标,以及y是所述电子乳房图像中的像素的y坐标;φ是光子通量,Vt是X射线管电压,Ap是像素面积,ts是曝光时间,N0 rel是该管电压在能量ε处的入射X射线光子的相对数量,S(ε)是检测器吸收,G(ε)是电网传输,hμ(x,y,ε)是归因于乳房的衰减,μluc(ε)hplate归因于有机玻璃压缩板的衰减,ε是管能量。
2.根据权利要求1所述的从乳房图像分析组织的方法,包括用于建立所述参考点的步骤,其中,
检测位于所述直方图的两个峰值之间的间隙(530,730);以及
所述间隙的位置的像素对应于所述内部乳房边缘(640,840,1040,1140)上的像素,所述内部乳房边缘用于为搜索所述参考点(1070,1170,1270)确定界限。
3.根据权利要求1和2中任一项所述的从乳房图像分析组织的方法,包括建立所述参考点,由此:
使用相位一致性检测对应于所述内部乳房边缘的位置;以及,
所述位置用于为搜索所述参考点(1070,1170,1270)确定界限。
4.根据权利要求1和2中任一项所述的从乳房图像分析组织的方法,包括建立所述参考点,由此:
通过确定阈值来检测对应于内部乳房边缘的位置,所述阈值通过查看像素值的乳房图像统计值来计算;以及
所述位置用于为搜索所述参考点(1070,1170,1270)确定界限。
5.根据权利要求4所述的从乳房图像分析组织的方法,其中像素值的乳房图像统计值包括:最小像素值,最大像素值以及平均像素值。
6.根据权利要求1、2和5中的任一项所述的从乳房图像分析组织的方法,其中,通过从所述乳房边缘中的数据外推参考像素值来建立参考点。
7.根据权利要求1、2和5中的任一项所述的从乳房图像分析组织的方法,包括:
建立所述乳房图像的厚度表示,其中所述厚度对应于脂肪组织,所述厚度表示可称为hfat表示;以及
通过以下操作来自动估计乳房成分:
计算厚度hint与脂肪组织的厚度hfat的值;
对hint与hfat的值求和;以及
随后计算乳房密度。
8.根据权利要求1、2和5中的任一项所述的从乳房图像分析组织的方法,其中,对在一段时间内收集的多个图像进行分析和自动比较。
9.根据权利要求8所述的从乳房图像分析组织的方法,其中通过产生对齐的标准化图像来便于对所述多个图像的比较。
10.根据权利要求9所述的从乳房图像分析组织的方法,其中,使用参考点乳头和胸壁来便于所述对齐。
11.根据权利要求1、2、5、9和10中的任一项所述的从乳房图像分析组织的方法,其中,一系列参考点被找到,以能够从不同的成像模式区分多于两种类型的组织。
Applications Claiming Priority (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NZ57876309 | 2009-08-03 | ||
NZ57876409 | 2009-08-03 | ||
NZ578764 | 2009-08-03 | ||
NZ578763 | 2009-08-03 | ||
NZ584397 | 2010-04-01 | ||
NZ58439710 | 2010-04-01 | ||
PCT/GB2010/001472 WO2011015818A1 (en) | 2009-08-03 | 2010-08-03 | A method and system for analysing tissue from images |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102549618A CN102549618A (zh) | 2012-07-04 |
CN102549618B true CN102549618B (zh) | 2015-11-25 |
Family
ID=42829518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080044509.6A Active CN102549618B (zh) | 2009-08-03 | 2010-08-03 | 用于从图像分析组织的方法和系统 |
Country Status (7)
Country | Link |
---|---|
US (1) | US9008382B2 (zh) |
EP (1) | EP2462561B1 (zh) |
CN (1) | CN102549618B (zh) |
DK (1) | DK2462561T3 (zh) |
ES (1) | ES2699471T3 (zh) |
PL (1) | PL2462561T3 (zh) |
WO (1) | WO2011015818A1 (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9146663B2 (en) * | 2008-12-08 | 2015-09-29 | Hologic, Inc. | Displaying computer-aided detection information with associated breast tomosynthesis image information |
CN102549618B (zh) | 2009-08-03 | 2015-11-25 | 马塔基纳科技有限公司 | 用于从图像分析组织的方法和系统 |
WO2011044295A2 (en) * | 2009-10-07 | 2011-04-14 | Hologic, Inc. | Processing and displaying computer-aided detection information associated with breast x-ray images |
EP2646980B1 (en) * | 2010-11-30 | 2017-01-11 | Volpara Health Technologies Limited | An imaging technique and imaging system |
US9304973B2 (en) * | 2010-12-15 | 2016-04-05 | H. Lee Moffitt Cancer Center And Research Institute, Inc. | Method for assessing breast density |
EP2675360B1 (en) * | 2011-02-14 | 2022-08-10 | University Of Rochester | Method and apparatus of cone beam breast ct image-based computer-aided detection and diagnosis |
US8417005B1 (en) * | 2011-09-15 | 2013-04-09 | Sunnybrook Health Sciences Centre | Method for automatic three-dimensional segmentation of magnetic resonance images |
CN103886571B (zh) * | 2012-11-13 | 2016-03-02 | 上海联影医疗科技有限公司 | 一种乳房三维图像的分割方法和装置 |
WO2014195669A1 (en) * | 2013-06-06 | 2014-12-11 | Matakina Technology Ltd | A method of reconstruction of an object from projection views |
KR20150018665A (ko) * | 2013-08-08 | 2015-02-24 | 삼성전자주식회사 | 엑스선 영상 장치 및 그 제어 방법 |
EP3035850B1 (en) | 2013-08-20 | 2020-05-13 | Densitas Incorporated | Methods and systems for determining breast density |
JP6200839B2 (ja) * | 2014-03-19 | 2017-09-20 | 富士フイルム株式会社 | 乳房厚測定装置、乳房厚測定方法及び放射線撮影システム |
DE102014208411A1 (de) * | 2014-05-06 | 2015-11-12 | Siemens Aktiengesellschaft | Auswertung eines während einer Mammographie erzeugten Röntgenbildes einer Brust |
US9256939B1 (en) * | 2014-07-17 | 2016-02-09 | Agfa Healthcare | System and method for aligning mammography images |
WO2016200983A1 (en) * | 2015-06-09 | 2016-12-15 | The Board of Trustees of the Leand Stanford Junior University | System for determining tissue density values using polychromatic x-ray absorptiometry |
US9858385B2 (en) | 2015-07-23 | 2018-01-02 | International Business Machines Corporation | Identifying errors in medical data |
DK3371774T3 (da) * | 2015-11-05 | 2021-10-04 | Volpara Health Tech Limited | Fremgangsmåde til kvantificeringen af billeder |
US9959617B2 (en) * | 2016-01-28 | 2018-05-01 | Taihao Medical Inc. | Medical image processing apparatus and breast image processing method thereof |
US9943280B2 (en) * | 2016-03-07 | 2018-04-17 | General Electric Company | Breast tomosynthesis with flexible compression paddle |
CN109478325A (zh) * | 2016-07-19 | 2019-03-15 | 唯盼健康科技有限公司 | 用于临床决策优化的系统和装置 |
US10282838B2 (en) | 2017-01-09 | 2019-05-07 | General Electric Company | Image analysis for assessing image data |
CN110892448B (zh) * | 2017-06-16 | 2023-10-03 | 唯盼健康科技有限公司 | 用于动脉钙化的检测和量化的方法 |
AU2020354898A1 (en) * | 2019-09-27 | 2022-02-17 | Hologic, Inc. | Motion detection for internal breast tissue in tomosynthesis |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1273516A (zh) * | 1997-08-28 | 2000-11-15 | 特性计算公司 | 从数字乳房x射线照片中自动检测成簇的微钙化的方法和系统 |
CN101292261A (zh) * | 2005-10-21 | 2008-10-22 | 卡尔斯特里姆保健公司 | 用于医用图像的增强可视化的方法 |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5452367A (en) * | 1993-11-29 | 1995-09-19 | Arch Development Corporation | Automated method and system for the segmentation of medical images |
AUPO848197A0 (en) | 1997-08-08 | 1997-09-04 | Breast Screening Decision Support System R&D Syndicate | Breast screening - early detection and aid to diagnosis |
US6282305B1 (en) | 1998-06-05 | 2001-08-28 | Arch Development Corporation | Method and system for the computerized assessment of breast cancer risk |
GB9904692D0 (en) * | 1999-03-01 | 1999-04-21 | Isis Innovation | X-ray image processing |
WO2003067371A2 (en) | 2002-02-08 | 2003-08-14 | Giger Maryellen L | Method and system for risk-modulated diagnosis of disease |
CA2391133C (en) | 2002-06-21 | 2011-02-15 | Jiwei Yang | Method and apparatus for measuring the thickness of compressed objects |
CA2391132A1 (en) * | 2002-06-21 | 2003-12-21 | Dan Rico | Method and apparatus for determining peripheral breast thickness |
US7418123B2 (en) | 2002-07-12 | 2008-08-26 | University Of Chicago | Automated method and system for computerized image analysis for prognosis |
GB2395263A (en) * | 2002-11-12 | 2004-05-19 | Qinetiq Ltd | Image analysis |
US7123684B2 (en) | 2002-11-27 | 2006-10-17 | Hologic, Inc. | Full field mammography with tissue exposure control, tomosynthesis, and dynamic field of view processing |
US8340373B2 (en) | 2003-12-23 | 2012-12-25 | General Electric Company | Quantitative image reconstruction method and system |
WO2006119426A2 (en) * | 2005-05-03 | 2006-11-09 | Regents Of The University Of California | Biopsy systems for breast computed tomography |
EP1938754A4 (en) * | 2005-10-19 | 2011-01-05 | Hitachi Medical Corp | ECHOGRAPH FOR CREATING AN ELASTIC IMAGE |
GB0602739D0 (en) | 2006-02-10 | 2006-03-22 | Ccbr As | Breast tissue density measure |
RU2325852C1 (ru) | 2006-08-08 | 2008-06-10 | ГУН НИИ онкологии им. проф. Н.Н. Петрова Росздрава | Способ оценки доли жировой ткани в молочной железе |
JP5056760B2 (ja) | 2006-09-26 | 2012-10-24 | コニカミノルタエムジー株式会社 | 画像撮影システム |
DE102006051778A1 (de) | 2006-11-02 | 2008-05-15 | Siemens Ag | Verfahren und Einrichtung zur Anzeige eines im Rahmen einer Mammographie aufgenommenen Röntgenbildes |
US7903861B2 (en) | 2006-12-28 | 2011-03-08 | Carestream Health, Inc. | Method for classifying breast tissue density using computed image features |
CN101234026B (zh) | 2008-03-07 | 2010-04-14 | 李立 | 一种乳腺病变定量化影像评价系统 |
CN102549618B (zh) | 2009-08-03 | 2015-11-25 | 马塔基纳科技有限公司 | 用于从图像分析组织的方法和系统 |
US8675933B2 (en) * | 2010-04-30 | 2014-03-18 | Vucomp, Inc. | Breast segmentation in radiographic images |
-
2010
- 2010-08-03 CN CN201080044509.6A patent/CN102549618B/zh active Active
- 2010-08-03 EP EP10754771.3A patent/EP2462561B1/en active Active
- 2010-08-03 ES ES10754771T patent/ES2699471T3/es active Active
- 2010-08-03 WO PCT/GB2010/001472 patent/WO2011015818A1/en active Application Filing
- 2010-08-03 US US13/261,165 patent/US9008382B2/en active Active
- 2010-08-03 PL PL10754771T patent/PL2462561T3/pl unknown
- 2010-08-03 DK DK10754771.3T patent/DK2462561T3/en active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1273516A (zh) * | 1997-08-28 | 2000-11-15 | 特性计算公司 | 从数字乳房x射线照片中自动检测成簇的微钙化的方法和系统 |
CN101292261A (zh) * | 2005-10-21 | 2008-10-22 | 卡尔斯特里姆保健公司 | 用于医用图像的增强可视化的方法 |
Non-Patent Citations (3)
Title |
---|
乳腺X线检查的临床应用及进展;姜洪;《黑龙江医学》;20010731;第25卷(第7期);493-495 * |
基于一致性直方图的超声乳腺图像分割;于瀛星;《中国优秀硕士学位论文全文数据库》;20061215;摘要,第3章 * |
基于特征点对齐度的图像配准方法;葛永新;《电子与信息学报》;20070228;第29卷(第2期);425-428 * |
Also Published As
Publication number | Publication date |
---|---|
ES2699471T3 (es) | 2019-02-11 |
EP2462561B1 (en) | 2018-10-31 |
WO2011015818A1 (en) | 2011-02-10 |
US20120189175A1 (en) | 2012-07-26 |
EP2462561A1 (en) | 2012-06-13 |
DK2462561T3 (en) | 2018-12-10 |
US9008382B2 (en) | 2015-04-14 |
PL2462561T3 (pl) | 2019-02-28 |
CN102549618A (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102549618B (zh) | 用于从图像分析组织的方法和系统 | |
US8571290B2 (en) | Automated quantification of digital radiographic image quality | |
US8634610B2 (en) | System and method for assessing cancer risk | |
Szczykutowicz et al. | Variation in CT number and image noise uniformity according to patient positioning in MDCT | |
EP2976746B1 (en) | A method and x-ray system for computer aided detection of structures in x-ray images | |
US11594005B2 (en) | System, method and apparatus for assisting a determination of medical images | |
US9361683B2 (en) | Imaging technique and imaging system | |
Badal et al. | Reproducing two-dimensional mammograms with three-dimensional printed phantoms | |
Kim et al. | Impact of radiation dose and iterative reconstruction on pulmonary nodule measurements at chest CT: a phantom study | |
Siegelman et al. | Pulmonary nodules with ground-glass opacity can be reliably measured with low-dose techniques regardless of iterative reconstruction: results of a phantom study | |
Lau et al. | Towards visual-search model observers for mass detection in breast tomosynthesis | |
Liu et al. | Artificial intelligence (AI) for lung nodules, from the AJR special series on AI applications | |
Mainprize et al. | The relationship between anatomic noise and volumetric breast density for digital mammography | |
US20100321404A1 (en) | Method and device for displaying an x-ray image recorded in the course of mammography | |
US11783476B2 (en) | System and method for analyzing three-dimensional image data | |
Clark et al. | Reproducing 2D breast mammography images with 3D printed phantoms | |
GB2474319A (en) | Analysing breast tissue image using reference spot and calibration error | |
Pinto Pereira et al. | Automated registration of diagnostic to prediagnostic x‐ray mammograms: Evaluation and comparison to radiologists’ accuracy | |
US10755454B2 (en) | Clinical task-based processing of images | |
EP4105880A1 (en) | An x-ray image analysis method and apparatus | |
CN108475423B (zh) | 量化图像的方法 | |
Miguel et al. | A Limited PACS Dedicated to Computer-Aided Diagnosis in Chest and Breast Radiography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C56 | Change in the name or address of the patentee | ||
CP01 | Change in the name or title of a patent holder |
Address after: Wellington, New Zealand Patentee after: MATAKINA TECHNOLOGY LTD Address before: Wellington, New Zealand Patentee before: Matakina Technology Ltd. |