JP2012512672A - 医用画像内病変自動検出方法およびシステム - Google Patents
医用画像内病変自動検出方法およびシステム Download PDFInfo
- Publication number
- JP2012512672A JP2012512672A JP2011541052A JP2011541052A JP2012512672A JP 2012512672 A JP2012512672 A JP 2012512672A JP 2011541052 A JP2011541052 A JP 2011541052A JP 2011541052 A JP2011541052 A JP 2011541052A JP 2012512672 A JP2012512672 A JP 2012512672A
- Authority
- JP
- Japan
- Prior art keywords
- intensity
- image
- pixel
- value
- normalized
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 116
- 230000003902 lesion Effects 0.000 title claims abstract description 87
- 210000001519 tissue Anatomy 0.000 claims abstract description 102
- 230000003211 malignant effect Effects 0.000 claims abstract description 58
- 210000004003 subcutaneous fat Anatomy 0.000 claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 25
- 210000000038 chest Anatomy 0.000 claims description 54
- 238000002604 ultrasonography Methods 0.000 claims description 54
- 210000000481 breast Anatomy 0.000 claims description 52
- 238000013507 mapping Methods 0.000 claims description 36
- 238000004422 calculation algorithm Methods 0.000 claims description 35
- 206010028980 Neoplasm Diseases 0.000 claims description 22
- 238000001514 detection method Methods 0.000 claims description 22
- 201000011510 cancer Diseases 0.000 claims description 16
- 230000036210 malignancy Effects 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 9
- 238000003064 k means clustering Methods 0.000 claims description 8
- 238000009499 grossing Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 208000031513 cyst Diseases 0.000 claims description 4
- 238000009877 rendering Methods 0.000 claims description 4
- 206010011732 Cyst Diseases 0.000 claims description 3
- 230000008520 organization Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 description 32
- 230000008569 process Effects 0.000 description 26
- 238000004458 analytical method Methods 0.000 description 21
- 210000002445 nipple Anatomy 0.000 description 13
- 230000009467 reduction Effects 0.000 description 13
- 238000004195 computer-aided diagnosis Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 9
- 239000000523 sample Substances 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 8
- 238000009792 diffusion process Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 210000003484 anatomy Anatomy 0.000 description 6
- 238000013459 approach Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000012549 training Methods 0.000 description 6
- 210000000577 adipose tissue Anatomy 0.000 description 5
- 210000002976 pectoralis muscle Anatomy 0.000 description 5
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 4
- 229910052791 calcium Inorganic materials 0.000 description 4
- 239000011575 calcium Substances 0.000 description 4
- 210000003195 fascia Anatomy 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000012552 review Methods 0.000 description 4
- 241001135902 Peanut clump virus Species 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 3
- 230000003176 fibrotic effect Effects 0.000 description 3
- 210000004907 gland Anatomy 0.000 description 3
- 210000003041 ligament Anatomy 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 206010006187 Breast cancer Diseases 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 2
- 238000007792 addition Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000010339 dilation Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000003706 image smoothing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 210000002307 prostate Anatomy 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000007920 subcutaneous administration Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000004434 Calcinosis Diseases 0.000 description 1
- 206010033128 Ovarian cancer Diseases 0.000 description 1
- 206010061535 Ovarian neoplasm Diseases 0.000 description 1
- 206010060862 Prostate cancer Diseases 0.000 description 1
- 208000000236 Prostatic Neoplasms Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007435 diagnostic evaluation Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 230000000762 glandular Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 210000001672 ovary Anatomy 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 210000000779 thoracic wall Anatomy 0.000 description 1
- 210000001685 thyroid gland Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4312—Breast evaluation or disorder diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
-
- 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/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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/30096—Tumor; Lesion
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Gynecology & Obstetrics (AREA)
- Mathematical Physics (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Quality & Reliability (AREA)
- Reproductive Health (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【解決手段】本発明は、医用画像を処理するシステムおよび方法を提供する。入力医用画像は初めに、皮下脂肪を含むコントロールポイント組織のピクセル強度を用いて正規化される。クラスタ化密度マップと悪性確率マップは、正規化された画像から生成され、病変を示す可能性のある、共通内部特性の部位、つまりブロブを同定するように、さらに分析される。これらのブロブは、真性の可能性のあるブロブを、医用画像でよく見られる他のタイプの非悪性腫瘤から識別するために、分析および分類される。
Description
本発明は、概して、医用画像のコンピュータ処理の分野に関するものであり、特に、医用画像における組織層の同定、医用画像における病変自動検出および医用画像におけるピクセル強度の正規化に関する。
多くの国において、癌は死亡の主な原因として認識されている。一般に、癌の早期発見および診断、それに伴う癌の早期治療は死亡率を減少させることに役立つと考えられている。乳癌、卵巣癌、前立線癌等の癌の検出および診断のための様々な画像技術が開発されている。例えば、乳がんの検出および診断のための現在の画像技術は、特に、他の技術の中で、マンモグラフィー、MRIおよびソノグラフィを含む。
ソノグラフィは超音波に基づく画像技術であり、一般には身体の軟組織の画像化に使用される。通常、トランスデューサは患者の体をスキャンするために使用される。生体組織および器官の超音波(US)画像は、トランスデューサにより受信された超音波エコーから生成される。撮像された腫瘤の形状、輪郭および縁の特徴的記述子ならびにエコー輝度は、通常、医用超音波画像の診断に使用される。ソノグラフィは、良性と悪性腫瘤の分類における効果的な撮像モダリティであることが証明されている。
しかし、超音波画像の正確な診断においては、放射線医の経験が重要な役割を担うことが多い。経験豊かな専門家により得られる感度と陰性適中率は、経験の浅い放射線医により必ずしも得られるものではない。さらに、スキャン技術は、悪性病変および良性病変の特徴を識別する定量化および能力に強い影響を及ぼす。このような強い影響はまた、経験値の異なる放射線医の間における超音波画像の一貫性のない診断の一因となる。
また、超音波画像の一貫性を有する分析は、絶対強度の変動によりさらに困難になる。組織タイプの絶対強度は、主に、画像取得中にハードウエア操作者により設定される利得係数などのオペレータ依存変数によって、異なる超音波画像の間で大いに異なる。利得係数は、組織のエコー輝度からグレーピクセル強度まで、マッピング決定における重要な役割を果たす。異なるオペレータによる利得係数の設定はスキャンの間で大きく異なることがあり、その結果、超音波画像の一貫性を有する分析をさらに困難にする。
オペレータに依存する他の設定としては、タイム・ゲイン・コントロール(TGC)設定があり、超音波画像の全体的な利得係数に密接に関係する。TGCは、組織深さに応じて強度マッピングに対するエコー輝度を調節する。組織深さは通常、ピクセルのY座標により示される。一貫性に欠けるTGC設定、または一貫性のないTGC設定の一貫した補償は、一貫性を有する統一された画像分析への別の困難をもたらす。
オペレータ依存の効果を克服し、胸部ソノグラフィの診断能を改善するために、コンピュータ支援診断(CAD)システムが開発されている。CADシステムの機能の1つは、コンピューターベース画像処理アルゴリズムを画像に適用することによる、超音波画像における疑わしい領域の自動検出および画定である。これは、ソノグラムにおける多くの検鏡ノイズおよび構造的作為により困難な作業となる。異なる画像取得条件は、一貫性のある画像分析をより困難にする。更なる困難としては、超音波画像における正常解剖構造の腫瘍状外観を含む。クーパー靭帯、腺組織および皮下脂肪は特に、真性病変と同一の多くのエコー輝度および形態学的特徴を共有する、正常な胸部解剖構造である。
本発明の目的は、上記困難の少なくとも1つを軽減または除去することである。
本発明は、医用画像における、組織層の同定、自動病変検出およびピクセル強度正規化に関する。
本発明は、医用画像を処理するシステムおよび方法を提供する。入力医用画像は初めに、皮下脂肪を含むコントロールポイント組織のピクセル強度を用いて正規化される。クラスタ化された密度マップと悪性確率地図は正規化された画像から生成され、病変を示す可能性のある共通内部特性の部位、つまりブロブを同定するためにさらに分析される。これらのブロブは、真性の可能性のある病変を、医用画像でよく見られる他のタイプの非悪性腫瘤から識別するために、さらに分析および分類される。
本発明の一態様において、超音波医用画像における疑われる病変の同定方法を提供する。この方法は、前記医用画像内の皮下脂肪ピクセルの推定代表脂肪強度値を計算するステップと、正規化画像を取得するために、正規化脂肪強度値と前記代表脂肪強度値の間のマッピング関係を利用して、前記医用画像のピクセル値から正規化グレーピクセル値を計算するステップと、個別領域を形成する前記正規化画像内のピクセルを同定し、前記個別領域の各々は一貫性のある内部特性を有するステップと、前記各個別領域から、記述的特徴を抽出するステップと、前記個別領域の前記抽出された記述的特徴を分析し、前記各個別領域が病変である尤度値を前記各個別領域に割当てるステップと、所定の基準を満たす尤度値を有する全ての個別領域を病変候補として同定するステップとを備える。
本発明の別の態様において、病変に対応する可能性のある医用画像における領域を自動同定するシステムを提供する。このシステムは、前記医用画像内のコントロールポイント組織の推定強度を、前記医用画像のピクセル値から計算するよう構成される強度部と、入力ピクセルと正規化ピクセルとの間のマッピング関係を生成し、グレーピクセル値を正規化ピクセル値に変換して、前記マッピング関係に従って正規化画像を取得するよう構成される正規化部と、パラメータ値を入力画像の各ピクセルに割り当てて、パラメータマップを生成するマップ生成部と、前記パラメータマップにあるブロブを、検出および画定するよう構成されるブロブ検出部と、前記検出ブロブの記述的特徴を、検出および計算するよう構成される特徴抽出部と、ブロブが悪性である推定尤度値をブロブの記述的特徴から計算し、前記尤度値を前記ブロブに割り当てるブロブ分析部と、を備える。
別の態様において、デジタルまたはデジタル化された医用画像における組織のグレースケール強度推定方法を提供する。この方法は、クラスタリング操作を医用画像のピクセルの強度値に適用して、前記強度値を個別強度クラスタに分類するステップと、前記個別強度クラスタの1つを、前記組織に対応する強度クラスタとして同定するステップと、前記強度クラスタの代表グレースケール強度値を、前記強度クラスタのピクセルのグレースケール強度から推定するステップと、前記代表グレースケール強度値を前記組織に割り当てるステップと、を備える。
本発明の別の態様において、超音波胸部画像処理方法を提供する。この方法は、胸部の層モデルを構築し、前記モデルの隣接層の各対は前記隣接層の各対間の境界面を定義するステップと、複数のサンプル超音波胸部画像でモデルをキャリブレーションし、前記複数のサンプル超音波胸部画像はそれぞれ、前記サンプル超音波胸部画像の境界面を同定するために手動でセグメント化され、キャリブレーションされた前記モデルはパラメータ化表面モデルを含み、各パラメータ化表面モデルは、サイズパラメータの離散値に対応する一組の境界面参照テーブル(LUT)を含むステップと、前記超音波胸部画像の前記サイズパラメータの推定値を受信するステップと、前記サイズパラメータの前記推定値に対応する新規表面モデルを、前記パラメータ化された表面モデルから確立し、前記新規表面モデルは、前記サイズパラメータの前記推定値に対応する一組の計算された境界面LUTを含むステップと、前記超音波胸部画像内の一次層のピクセルを同定するために、境界面の推定位置を、前記新規表面モデルの前記一組の計算された境界面LUTから計算するステップと、を備える。
本発明のさらに別の態様において、超音波胸部画像における病変同定方法を提供する。この方法は、一次層組織を分離する表面の推定位置を計算し、前記一次層組織は***領域の組織を含むステップと、前記***領域のピクセルを同定するステップと、前記***領域の各ピクセルに対してピクセル特性ベクター(PCV)を構築し、前記PCVは、少なくとも前記ピクセルの各々の近隣の特性を含むステップと、前記***領域の各ピクセルに対して、前記各ピクセルのPCVから悪性確率値を計算し、前記ピクセルのそれぞれに悪性確率値を割り当てて、前記割り当てられた悪性確率値が閾値を上回る場合は、ピクセルを病変の可能性のあるピクセルとして同定するステップと、全ての病変の可能性のあるピクセルの隣接領域を、潜在病変として報告するステップと、を備える。
別の態様において、本発明は、上記の態様の様々な組み合わせおよびサブセットを提供する。
限定ではなく説明の目的で、本発明の上記及び他の態様を、添付図面を参照して詳細に説明する。
下記の説明およびその中の実施例は、本発明の原理の特定の実施例の一つまたは複数の例を説明するために与えられたものである。これらの例は、これらの原理と発明を限定することなく説明を目的として与えられたものである。下記説明において、同様の要素は明細書および図面の全体にわたり同じ符号を付す。
本発明は、概して、医用画像を処理するシステムおよび方法に関する。具体的には、本発明は、超音波医用画像における病変候補の検出に関する。
一実施例において、一連の画像処理ルーチンは、より一層の診断的再検討を要することがある各病変候補の検出および分類のために、1つの胸部超音波画像(またはボリュームデータセット)等の入力画像に適用される。図1は、処理100の概略を提供するフローチャートである。
準備ステップ110として、初めに、入力画像(またはボリューム)が処理のため受信される。入力画像は、画像システムにより取得した画像を保存および提供する、医用デジタル画像と通信(DICOM)アーカイブ等の画像保存装置から読み出すことができる。また、入力画像は超音波プローブ等の画像取得装置から直接受信することができる。取得パラメータは画像データと共に読み出しが可能である。取得パラメータは、深度およびトランスデューサ周波数を含む、異なる業者の超音波トランスデューサ機器の間の差異によるハードウエアパラメータと、例えばトランスデューサ圧力や時間利得補償等の技術者の機器または取得設定等のオペレータパラメータとを含む。これらのハードウエアパラメータおよびオペレータパラメータは、DICOM基準に定義されるデータヘッダから抽出、または、画像システムにより取得された画像がリアルタイムで処理される時に、画像データと共に直接送信できる。
この準備ステップ110に続いて、処理100はノイズ除去、すなわちノイズ除去またはノイズ低減ステップ120を有し、入力画像からノイズを除去または低減する。ノイズは、例えばエッジ保存拡散アルゴリズムを入力画像に適用することにより、除去または低減することができる。このようなアルゴリズムは、画像内の物体のエッジを保存および好ましくは強調しながら、画像からノイズを除去または低減するために使用することができる。
次のステップは、画像強度の正規化(ステップ130)である。当業者にはよく知られているように、ほとんどの医用画像モダリティのハードウエア装置により生成されたピクセルの強度は、通常、画像取得ハードウエアの差異によりもたらされる非一貫性に見舞われる。これらはまた、超音波画像取得中にハードウエアオペレータに選択される利得係数等、オペレータにより適用される取得技術の不一致に見舞われる。これらの不一致を減少させることが望ましい。非一貫性を減少させる1つの手法は、例えば、一貫して目に見える皮下脂肪等、よく特徴づけられる共通の組織タイプをコントロールポイントとして選択し、このコントロールポイントに対して画像強度を正規化することである。理想的には、強度は、入力画像から動的に判定または測定されたコントロールポイント組織の代表強度に対して正規化される。コントロールポイント、つまりコントロールポイント組織の代表強度は、入力組織タイプから出力強度まで、マッピング関数を確立する。コントロールポイントは、例えば、皮下脂肪の計算された平均強度値であり、出力ダイナミックレンジの中間点にマッピングされる。皮下脂肪は、一貫して、胸部超音波(BUS)画像の上部に非常に近接した皮膚線の下に表れる。BUS画像では、皮下脂肪は強度正規化のための信頼できるコントロールポイント組織であると考えられている。コントロールポイントとして通常選択することができる他の画像要素には、例えば、無響嚢胞等の器官組織、皮膚組織、線維腺組織、カルシウムがあるが、必ずしもこれらに限定されない。前立腺または甲状腺等、相当量の皮下脂肪が必ずしも存在しない器官においては、別の組織タイプが正規化処理のために持続的にみられることがある。これらの別のコントロールポイント組織は、これらの他の撮像された器官に通常見られる典型的な解剖構造を説明するために使用することができる。
正規化は、ピクセル値をそれぞれの初期値から正規化値へマッピングすることであり、すなわち、初期値と正規化数との間のマッピング関係に従ってそれぞれの初期値から正規化値への変換することである。画像は、平均脂肪強度に基づいたマッピング関係に従って正規化することができる。皮下脂肪が位置することが予想される画像領域が選択され、その後、例えば強度クラスタリングアルゴリズムを選択領域に適用することにより皮下脂肪の強度値が計算される。k平均アルゴリズム等のロバストクラスタリングアルゴリズムを、皮下脂肪の強度推定値の計算に使用することができる。他のロバストクラスタリングアルゴリズムとして、ファジーc平均またはR.O. Dudaによる"Pattern Classification", John Wiley & Sons Inc., 2001に記載された期待値最大化クラスタリング技術などがある。クラスタリングアルゴリズムは、ピクセル強度によりグループ化されたいくつかのクラスタを生成する。皮下脂肪に対応するクラスタが同定され、同定されたクラスタ内のピクセルのピクセル値から代表脂肪強度としての平均脂肪強度が計算できるようになる。入力画像を正規化するためのマッピング関係は、脂肪組織のグレーレベルが中間レベルのグレーに表れるように確立することができる。これは、クラスタリングアルゴリズムにより計算された代表脂肪強度を変換するために、マッピング関係を確立する。他の組織または撮像対象を示すピクセルの強度値は、マッピング関係に従って、それぞれの入力値から正規化値に変換される。このステップでは、「脂肪強度正規化画像」が出力される。他のコントロールポイント組織を含むこともできる。便宜上、コントロールポイント組織の検出強度からそれらの各正規化強度へのマッピングを確立することができる。適切な補間を有するマッピング関係は、画像を正規化し、またより一貫性のある正規化画像を生成するために使用することができる。正規化画像の画像強度は、強度と組織のエコー輝度との間のより一貫性のあるマッピングを提供する。
正規化画像は、次に、一致するまたは同様の内部強度特性を有する隣接ピクセルの個別エリア、つまり「ブロブ」を検出するよう処理される(ステップ140)。異なるブロブ検出方法を採用することができる。一般に、まず、選択パラメータに対して、パラメータマップ、すなわち画像の各ピクセルにおけるパラメータ値の空間的変動を生成する。次に、閾値を上回る、閾値を下回るまたは所定の範囲内等、一定の基準を満たすパラメータ値を有し、且つ、個別のエリアを形成する隣接ピクセルは、ブロブに属すると同定され、各個別エリアは一つの検出ブロブとなる。パラメータは、結果として生じるマップ、特に、検出されたブロブが病変または他の下層組織構造の検出を支援するよう選択される。このようなマップは、脂肪正規化画像のグレーピクセル強度値に基づいた密度マップである。このようなマップ内のブロブはピクセルの強度クラスタに対応し、これは一般に認められている相対エコー輝度に基づいて、胸部を構成する胸部組織の対応クラスに分類することができる。他のパラメータマップを生成するために他のパラメータを選択することも可能であり、これらについては後述する。
便宜上、潜在的病変のエコー輝度をいくつかのカテゴリの1つとして分類するためにBI-RADS地図を使用することができる。
1.無響:内部エコーがなく暗穴に類似する
2.低エコー:脂肪と関連性があると定義され、領域全体にわたる低レベルのエコーを特徴とする
3.等エコー:脂肪と同一のエコー輝度を有する
4.高エコー:脂肪に比べ増加されたエコー輝度、または線維腺組織に等しい
5.複合:低エコー(嚢胞性)とエコー輝度(固形体)の両要素を含む
1.無響:内部エコーがなく暗穴に類似する
2.低エコー:脂肪と関連性があると定義され、領域全体にわたる低レベルのエコーを特徴とする
3.等エコー:脂肪と同一のエコー輝度を有する
4.高エコー:脂肪に比べ増加されたエコー輝度、または線維腺組織に等しい
5.複合:低エコー(嚢胞性)とエコー輝度(固形体)の両要素を含む
クラスタリングアルゴリズムは、それらのグレーピクセル値に基づいてピクセルをクラスタリングする密度マップに適用可能である。例えば、結果として生じるクラスタは、上術の様々なBI-RADSカテゴリに対応する密度マップ内の領域を描写することができる。この処理は通常、多くのクラスタまたは隣接ピクセルの領域、すなわち、個別の画像領域を生成し、それらのうちいくつかの強度は潜在的病変の範囲にある。これらの領域または形状のそれぞれは「密度ブロブ」として特定される。
上述の通り、グレーピクセル密度値に基づく密度マップに加えて、他のパラメータを使用することもできる。そのうちの1つはピクセルの悪性である確率値である。各ピクセルに対して該ピクセルの悪性確率値が計算され、各ピクセルに割り当てられる。その結果、悪性確率マップが生成される。悪性確率マップの生成方法については、後述する。75%または他の適切な値の閾値を上回る悪性確率を有するピクセルを、複数の個別の領域にグループ化することができる。これらの個別領域は、それぞれ「悪性ブロブ」として同定される。
同定されたブロブの全てが真性の病変ではない。それらのいくつかは、真性の病変ではなく、偽陽性病変候補でありうる。偽陽性病変候補の数を減らすために、ステップ150では、ブロブの特徴に基づく分析が行われる。この特徴に基づく分析の詳細は後述する。簡潔に言えば、各ブロブの記述子が、各ブロブの特徴を数値化するために推定される。これらの記述子は通常、形状、方位、深度および局部的背景に対するブロブコントラスト等の特徴や大域的皮下脂肪の強度に関連する。
本処理の次のステップ160では、ステップ150で推定された記述子を使って、専門家が同定した病変に一致する真性病変候補と、誤って報告された候補との間の微妙な相違を識別する。偽陽性は、除去される。真性病変の可能性のあるものと偽陽性候補とを識別する1つの手法は、各ブロブの記述子をCART(Classification And Regression Trees)アルゴリズムを介して供給することである。CARTアルゴリズムは、初めにトレーニング画像の代表的な一組でトレーニングされる。CARTアルゴリズムをトレーニングするために、トレーニング画像の各画像から抽出または計算されたブロブの特徴は、それらの各記述子と、対応する専門家による分類とに関連付けられる。ステップ160では、ステップ150で推定された記述子が、トレーニングされたCARTアルゴリズムに供給される。その結果、あるブロブが病変であるという推定確率が生成され、その値がブロブに割り当てられる。ステップ160において、閾値を下回る、すなわち所定の基準を満たさない推定確率を有するブロブは、偽陽性として扱われ除去される。ステップ170で、残ったブロブのみが、さらなる再検討と研究のための病変候補として同定および報告される。
図2は医用画像の処理および診断のためのCADシステム200の概略図である。CADシステムは医用画像の1つまたは複数のソースと通信する。ソース202は、患者の超音波画像がリアルタイムで取得される超音波画像システム等の医用画像取得システムで構成されることができる。ソースは、画像システムにより取得した画像を1つまたは複数のコンピュータ可読メディアに保存したDICOMアーカイブ等の画像アーカイブであっても良い。また、ソースは、医師により既に読み出された、医師のコンピュータシステムの局部的な保存メディアに保存された画像データであっても良い。画像ソース202と連動する画像読出部204は、入力画像データを受信する。当該分野の当業者にとって明らかなように、画像読出部204はまた、CADシステム200により提供される画像読出し機能でもあり、必ずしも1つまたは複数の特定のモジュールに属するわけではない。取得パラメータ部206は、保存された、もしくは医用画像データと共に転送された取得パラメータを抽出する。取得パラメータ部206は、DICOMデータを処理し、これらのパラメータを画像データ内のDICOMデータヘッダから抽出する。また、非標準データ形式に対応するように構成され、任意の独自形式で保存された画像データからこれらのパラメータを抽出しても良い。
ノイズ除去等の前処理部、つまりノイズ低減部208は、ノイズレベル低減のために提供される。この目的のため、任意の適切なノイズ低減または除去アルゴリズムを組み入れることができる。画像読出部204は、受信画像を前処理部へ送る。前処理部は、ノイズ低減のために組み入れられたノイズ低減又は除去アルゴリズムを、受信画像に適用して、ほとんどのUS画像に見られる周知のスペックルノイズアーチファクト等のノイズを低減する。
システム200は強度測定部210、つまり強度部を含む。強度測定部はノイズ低減画像等の画像をノイズ低減部208から受信し、脂肪または皮膚の平均強度または中間強度等の選択されたコントロールポイント組織の代表強度を測定する。組織強度を測定するために異なる方法を実施することも可能である。例えば、使用者は画像内のある領域をコントロールポイント組織に属するとして同定することができ、次に強度測定部が使用者に同定されたその領域内の各ピクセルの強度値を推定し、そこからコントロールポイント組織の平均強度を計算する。クラスタリングにより強度値を抽出する等、より精緻な方法を使用することも可能である。平均強度値を計算するためのk平均クラスタリングを使用する例は、後述する。
また、システム200は強度正規化部212、つまり正規化部を含む。強度正規化部212は、コントロールポイント組織の代表強度に基づいて、画像のピクセル強度を正規化し、これにより、正規化後、異なるハードウエアユニットまたは異なるオペレータにより取得された画像は、同タイプの組織の対しより一貫性のある強度を有する。強度正規化部212は通常、ノイズ低減部208により前処理されたノイズ低減画像を入力とするが、画像読出部204から直接転送された画像を正規化することもできる。強度正規化部212は強度測定部210からの出力、すなわち、コントロールポイント組織の代表強度値を使って画像を正規化する。脂肪の強度のみを正規化処理に組み込む場合、強度正規化部は脂肪正規化画像を出力する。通常、一組のコントロールポイント組織の強度が強度正規化部212により取り入れられ、その結果、一般強度正規化画像が得られる。一組のコントロールポイント組織に基づいた画像の正規化に採用される方法の詳細については後述する。
また、システム200は、密度マップ生成部216および悪性マップ生成部218等の1つ以上の異なるマップ生成部を含むマップ生成モジュール214を含む。これらのマップ生成部は、密度値または悪性確率値等のパラメータ値を各ピクセルに割当てる。その結果、密度マップまたは悪性確率マップが生成される。密度マップ部216は、各ピクセルに正規化グレーピクセル値、すなわち対応する密度値に割り当てることにより、密度マップを生成する。悪性マップ生成部218は、通常はノイズ除去され強度が正規化された画像内の各ピクセルに、ピクセルの悪性確率値、すなわち組織の悪性領域に属するピクセルの確率値を割当て、それにより悪性確率マップを生成する。また、胸部解剖マップ(BAM)部218'を有することも可能である。BAM部218'は、正規化画像またはノイズ除去画像等の入力医用画像を受信し、画像の各ピクセルを可能性のある組織タイプに分類する。分類されたピクセルを有する画像は、悪性マップ生成部218によってさらに処理することができる。ピクセルの悪性である確率値は各ピクセルに割り当てることができ、これによっても悪性マップを生成できる。これらのマップ生成部で実行可能な処理は、詳細に後述する。
これらのマップは、ブロブにクラスタ化される。ブロブ検出部220は、マップ内のピクセルをクラスタして、各ブロブを含む関心領域(ROI、以下「ブロブROI」)を検出するために、提供される。ブロブは、上述のように、クラスタリングにより、または所定の基準を満たすパラメータを有するピクセルをグループ化することにより検出することができる。ブロブの境界を追跡する、もしくは境界を判定することにより、ブロブ検出部220はまた、検出されたブロブを画定する。特徴抽出部222は、検出されたブロブから特徴または特性を抽出するために提供される。特徴または記述子の異なるカテゴリを定義し分類することができる。例えば、ブロブの形状に関する特徴や、グレーレベルの変動に関する特徴、または撮像された領域内の解剖構造に関するブロブの空間的位置に関する特徴などが挙げられる。特徴抽出部222は、特徴の定義された各カテゴリに従って特徴又は記述子を抽出する、すなわち検出および/または計算をするよう実行される。当然、定義された特徴のカテゴリが増すとともに、特徴抽出部222の機能性は、特徴の幅広い領域に対応するよう拡張される。
検出されたブロブは、さらに分析される。例えば、詳細に後述するように、形状、稠密度および伸長等のブロブの形態学的特徴が計算および分析される。分析の前に、ブロブは、ブロブROI内のいかなる「穴」も埋めるなどの形態的変更、または形状学的フィルタリングを用いたROIの境界輪郭の平滑化等を受けることができる。これらのブロブは、ブロブ分析部224により、抽出された特徴と、適応可能な場合は各特徴に割り当てられた数値とを考慮して分析される。この分析の結果は、ブロブが悪性、すなわち真性病変の可能性があるという推定尤度値を計算するために組み合わされる。ブロブ分析部224はまた、一旦その値が計算または決定されると、推定された尤度値をブロブに割り当てる。所定の閾値を上回る尤度値を有する全てのブロブは、放射線医による更なる検討のために報告されることができる、もしくはさらなる自動診断的評価を受けることができる。報告する病変候補(または疑わしいブロブ)の同定およびこれらの病変候補の報告は、病変報告部226により実行される。
さらなる改善として、システム200はまた粗胸部解剖マップ(CBAM)モデリング部228を含むことができる。詳細は後述するため簡単に述べるなら、CBAMモデルは胸部の層状モデルであり、胸部の一般的解剖構造とマッチングさせるために、胸部画像をいくつかの一次層に分割したものである。CBAMモデリングは、皮下脂肪または***領域等の、一次層の位置を推定することに対する自動的なアプローチを提供する。境界面の推定位置は、例えば、脂肪の強度を推定する強度検出部210により、またはBAM’部218により、***領域のピクセルのみを分類するために使用することができる。CBAMモデリング部228により実行することができる処理の詳細は、後述する。
図3は、超音波画像等の医用画像を自動区画し、医用画像内で検出された区画された複数の腫瘤を病変候補と偽陽性とに自動分類する処理を示す。この方法は図2に示すCADシステム、またはCADの一部と、数ある可能性の中でも特に、画像取得ハードウエアに埋め込まれた画像取得システムとを用いて実行することができる。
最初のステップは、医用画像データと、それに関連する画像取得パラメータとを含む入力画像データを受信することである(ステップ302)。これは、例えば、画像読出部204と取得部206とにより実行することができる。次に、入力医用画像は、ノイズ低減のために前処理される(ステップ304)。超音波入力画像内の通常のハイレベルノイズを低減するために、このステップは通常、ノイズ低減およびノイズ除去を含む。ノイズ除去技術はノイズを除去するのみでなく、画像エッジの位置を不鮮明または変更することなく除去するべきである。例えば、入力画像は、画像エッジを維持および強調する一方で、超音波画像からノイズを除去するエッジ保存拡散アルゴリズムにより強調することができ、これにより、極めて限局させることを確実にする。従って、このようなノイズ除去ステップは、ノイズ除去、画像平滑化およびエッジ強調の目的を同時に達成することができる。ノイズ低減部288(図2参照)は、ノイズ除去ステップを行うために、任意の適切なノイズ除去、すなわちノイズ低減および除去アルゴリズムを実施することができる。画像エッジの位置を不鮮明または変更せずに画像に収められた特徴のエッジを強調する目的で、ノイズ低減またはノイズ除去アルゴリズムが選択されることが好ましい。さらに、エッジ強調またはノイズ除去処理は、オペレータによる設定またはハードウエアフィルタリングによる画像の特徴の内在的相違を計上するために、画像取得パラメータの関数として構成される。
多くの異なる適切なエッジ保存画像フィルタリングアルゴリズムが使用可能であるが、以下の説明では、これが唯一の適切または利用可能な方法ではないという理解のもと、非線形拡散法を使用する。非線形拡散法は、関心対象のエッジを保存する一方で、入力画像内の無関係又は誤った細部を除去するためによく使用される、周知の画像処理強調技術である。非線形拡散平滑化は、インターリージョン平滑化よりもイントラリージョン平滑化を優先して促す選択的フィルタリングであり、エッジの鮮明さを保存する。この方法は、非線形偏微分方程式の反復的解決
の形式を有し、kはコントラストパラメータの役割を果たし、すなわちkよりも大きい勾配値の構造はエッジであるとみなされ、拡散率は0に近く、kよりも小さい勾配値は内部領域に属するとみなされる。アルゴリズムは、Weickert, Jによる"Anisotropic Diffusion in Image Processing", ECMI Series, Verlag, 1998に記載されている。エッジ保存拡散の適用例は、図4に示す。図4aはノイズ除去アルゴリズム適用前の入力画像を示し、図4bは平滑化画像を示す。
病変候補の検出を単純化するため、入力画像Iは画像ピクセル強度と組織エコー輝度との間の一貫性のあるマッピングを確保するために正規化されて、画像間における利得係数設定の変動性を軽減させる。正規化ステップ306は、ノイズ除去画像に適用される。正規化されたピクセル値は、画像に示された組織に対し、より一貫したピクセル値の領域を有する。皮下脂肪のエコー輝度は、正規化画像の中間レベルグレー強度により示されることが好ましい。典型的な8ビットの超音波画像では、強度値の中間点は強度127に相当し、グレーレベルの範囲は0から255の間となる。強度正規化を適用するために、まず脂肪の強度がステップ308にて検出される。正規化ステップ306の結果、脂肪正規化画像が生成される。より一般的な手法では、皮下脂肪に加えていくつかのコントロールポイント組織の強度が、例えば強度測定部210により測定または検出される(ステップ310)。マッピング関係は、マッピング参照テーブル(LUT)により示される。マッピングLUTは、コントロールポイント組織の代表強度およびそれらに割り当てられた各値から計算される(ステップ312)。次に、マッピングLUTに従って、画像を正規化する(ステップ306)。以下は、脂肪に基づく正規化を最初に説明する。
図5は、皮下脂肪の強度を判定し、その平均強度を脂肪強度の代表値として使用して原入力画像を正規化する自動処理500を示す。超音波画像取得プロトコルは、解釈を容易にするために、一般に、強度に対するエコー輝度の均一なマッピングを確保するようにソノグラファが、時間利得補償設定を構成するよう促す。このことは、任意の空間的変動性を考慮して、ここに記載する方法の一層の改良を行うことができることは当然であるが、マッピングの空間的変動性が最小であるという仮定を容認する。
この方法は、超音波入力画像等の原入力画像を受信することにより開始する(ステップ502)。次に、主に皮下脂肪ピクセルから構成される入力画像内のROIを選択または同定する(ステップ504)。実験的に、皮下脂肪領域の前面の標準的深度はおよそ1.5mmであり、皮下脂肪領域の後面の標準的深度はおよそ6mmである。皮下脂肪領域の正確な位置における変動に関わらず、深度1.5mmと6mmの間の画像領域のかなりの部分は、脂肪により構成される傾向がある。
脂肪領域は、ターゲット深度の間のノイズ除去画像をトリミングすることにより選択できる。従って、ROIは、原入力画像の脂肪領域を示すために取得される。乳頭領域の周辺等の一部の領域では、皮下脂肪領域は、画像の最上部におけるより通常の位置よりも後方となる。脂肪領域の選択は、画像領域内の乳頭、乳頭パッドまたは他の特徴の存在を検出することにより、一層改良することができ、必要な場合は、これらのケースに対応するために、皮下脂肪画像ストリップの推定位置の輪郭を移動または変更する。胸部超音波画像の組織タイプの深度のモデリング等の他の方法も、皮下脂肪領域の境界線を描写するために採用することができる。そのようなモデリング方法の一つとして、いわゆるCBAM方法を簡潔に後述する。
次に、ステップ506にて、k平均クラスタリングアルゴリズム等のロバストクラスタリングアルゴリズムが、皮下脂肪領域の強度に適用される。k平均クラスタリングアルゴリズムの使用につきここで説明するが、前述のように、k平均アルゴリズムの代わりに、ファジーc平均または期待値最大化クラスタリング技術等の他のロバストクラスタリングアルゴリズムを使用することもできる。k平均クラスタリングアルゴリズムは、k=3である場合、皮下脂肪画像ストリップ等の脂肪領域を3つのピクセル強度クラスタに分割するよう構成される。ストリップの無響および低エコー領域は最低ピクセル強度のクラスタにより同定され、ストリップの等エコー領域は中間レベル強度のクラスタにより同定され、高エコー領域は高強度クラスタにより同定される。
ステップ508において、中間レベル強度クラスタの強度を計算または抽出する。ステップ510において、代表強度、すなわち本例においては平均脂肪強度が、抽出された強度から計算される。このクラスタリング操作の結果、画像ストリップ内の脂肪の強度のロバスト推定値が得られる。脂肪組織は同一の強度を有することが期待されるため、脂肪組織が皮下ストリップ内または他所に位置するかに関わらず、この推定値は画像全体にわたる脂肪の代表強度として使用することができる。最後に、ステップ512にて、ステップ502で受信された原入力画像が、推定された脂肪強度を使用して正規化され、その結果、脂肪正規化画像が生成される。正規化処理は、本明細書においてさらに詳細に説明する。
前述のように、皮下脂肪領域はまた、胸部超音波画像の様々な組織タイプの深度のモデリングを通して同定されることができる。図6は、4つの一次層を含む胸部超音波画像の通常の構造を示す。皮膚層602は、画像最上部近くで明るい水平領域として表れ、次に、***後方筋膜608により腺領域または***領域606から分離された皮下脂肪604の均一領域(多くの場合、水平画像ストリップの形状である)が続く。画像の下端には、***後方領域610、すなわち胸壁領域、通常は胸筋および肋骨が、***領域606から筋膜608により分離された暗い領域として、再度示される。
皮膚線612、すなわち皮膚層602の外面は、様々な一次層の深度を測定する堅固な基準表面を提供する。各一次層始まりと終わりの深度は、皮膚線612から一次層とその隣接する一次層との間の境界面までの距離により便宜的に測定することができる。例えば、皮膚層602と皮下脂肪604とを分離する第一境界面614の深度により、皮膚層602の厚さの測定値が得られる。同様に、皮下脂肪604の厚さは、皮下脂肪604と***領域606の間の第二境界面616の深度を測定し、第一境界面614と第二境界面616の深度の差を計算することにより得ることができる。また、***領域606と***後方領域610との間の第三境界面618の深度が分かる場合、***領域606の厚さは、第二境界面616と第三境界面618の深度の差から計算することができる。
周知の通り、一次層の厚さは胸部にわたり、胸部のサイズにより変化する。有利なことに、粗胸部解剖マップ(CBAM)モデルは、胸部内の一次層の近似位置のモデルを作る、すなわち推定するために確立することができる(図6)。一次層の位置は、隣接層から分離する境界面により示すことができる。CBAMモデルは、一組の境界面参照テーブル(LUT)を使用して、境界面の位置を示す。胸部のサイズのばらつきを考慮するため、異なる組の境界面LUTをパラメータ化するために胸部のサイズを反映したパラメータが選択され、従ってパラメータ化されたCBAMモデルとなる。このようなサイズパラメータのひとつは、皮膚表面から測定された境界面の最大深度であってもよい。
下記において、最大***領域深度MZDmaxは、パラメータ化された胸部組織の層状モデル生成、キャリブレーションおよび適用を説明ためのサイズパラメータとして使用される。***部深度MZDは、皮膚からの***部の深度、すなわち皮膚線612と、***領域606を後方***領域610から分離する第三境界面618との間の距離である。通常、MZDmaxは乳頭の位置としばしば一致する領域である胸部の中央領域にあり、これは顕著な解剖学的特徴である。撮像された胸部のサイズを反映する他の任意の適切なサイズパラメータを、モデリングのために選択できることに注意されたい。
図7は、胸部の一次組織層のパラメータ化モデルを確立し、一時組織層の位置、つまりBUS画像内の皮膚の深度、皮下脂肪、***組織および後方***組織層の位置を推定する処理700のフローチャートである。なお、本処理は二次元(2D)および三次元(3D)胸部超音波画像のいずれに対しても等しく適用可能である。
処理700は、大きく分けて3つの段階を含む。第1段階は、一次層のパラメータ化モデルを構築することである。パラメータ化モデルは、多くの入力画像により構築され、サイズパラメータの選択範囲にわたり生成される。第2段階では、新規BUSデータセットを受信すると、モデルのサイズパラメータの値であるMZDmaxが取得パラメータから推定または判定される。新規BUSデータセットに対するサイズパラメータの推定値は、新規BUS画像に対する新規CBAMモデルを動的に生成するために、パラメータ化モデルに渡される。第3段階では、一次層の境界面の位置が、新規CBAMモデルから計算される。これらの各段階における、各段階のステップおよびそれらのいくつかの変動について、下記で詳述する。
第一段階では、モデルを構築する段階であり、例えば、胸部を4つの一次層に分けて、その後、パラメータ化モデルを多くのサンプルBUS画像でトレーニング、すなわちキャリブレーションすることにより、4つの一次層の推定位置を計算することによる、物理的胸部解剖の表現を生成する。各サンプルBUS画像は、手動でセグメント化される、すなわち専門家によりマークされた一次層の間の境界面を有する。モデルはまた、胸部の全体厚み領域に広がるように(例えば、2cmから5cmの間)、多くのサイズパラメータ値、すなわち最大領域深度を用いてキャリブレーションされる。
図7を参照すると、第一ステップ710において、多くのサンプルBUS画像データが読み出され、各サンプル画像が手動でセグメント化される。各サンプル画像と共に、ピクセルサイズ、トランスデューサ角度、最大***領域深度値MZDmaxおよび胸部クロックマップの胸部上のトランスデューサの位置等の、スキャン取得パラメータが受信される。次に、ステップ720において、スキャン取得パラメータを使って各セグメント画像に対する各境界面の位置が計算される。例えば、ピクセルサイズパラメータは深度値とピクセル値の間で変換するために使用することができる。同様に、トランスデューサスキャン角度は(利用可能であれば)、トランスデューサの配向が皮膚表面に対して垂直でない場合はいつでも、三角形状に基づいて深度値を再計算するために使用できる。次に、ステップ730において、各境界面に対して、多くの表面参照テーブル(LUT)が生成され、各LUTはMZDmax値、つまりトレーニングBUS画像のMZDmax値のビン(bin)に対応する。表面LUTにより境界面の位置を計算することができる。各表面LUT(およびその対応する境界面)は、乳頭を参照して測定されるような、胸部上の位置の関数である。乳頭は、MZD値が最も高いMZDmaxの位置であることが多い。
図8および図9は、MZD表面、つまり第3境界面618の例であり、***領域606と後方***筋膜608との間の境界面である。図8は、通常は***領域組織の層が最も厚くなる乳頭804から始まる、胸部クロックマップ802の位置に対するMZD値の変化例を示す。この例は、MZDmax値が3cmの中央領域804を示し、乳頭の位置と一致する。従って、乳頭の位置におけるMZD値は、MZDmaxパラメータの100%として示される。遠位細胞810は、乳頭804の位置から離れるに従って、MZD値が小さくなる。この例は、ほぼ対称なMZD表面を示す。例えば、中央領域804までほぼ同じ距離にあるセル806、808は、ほぼ同じMZD値を有するが、実際には、表面の形状は非対称であることが多い。MZD表面910の3Dビューを図9aに示し、2Dプロファイル920を図9bに示す。両者とも、MZDmaxのポイント、つまり乳頭の位置930から離れるにつれてMZD値が減少する一般的な傾向、およびMZDmaxのポイントの周囲におけるほぼ対称な性質を明示している。上述のように、MZD表面に対して表面LUTが生成される。同様に、他の一次層境界、例えば図8および9に示すような、皮膚層602と皮下脂肪604との間の境界面、そして皮下脂肪604と***領域606との間の境界面を計算することができ、これらの境界面に対して類似のLUTを生成することができる。
離散MZDmax値の3D(または2D)LUTが、保存される。これらはその後、新たな一組の新規LUTが新規胸部USスキャン、すなわち新規MZDmax値に対して計算される時に、利用される。その後、ピクセルサイズパラメータが、ピクセル値に対するMZD値を図るために使用され、トランスデユーサスキャン角度は、トランスデューサが皮膚表面に対して垂直でない場合に、三角形状に基づいてMZD値の再計算に使用されることができる。
次の段階では、新規胸部US画像、すなわちステップ302で受信した医用画像の粗胸部解剖マップ(CBAM)モデルを生成する。図7を参照すると、画像データと関連するスキャン取得パラメータがステップ740で読み出される。または、既にステップ302で受信している場合は、単にCBAMモジュール228を処理のため転送する。上述のように、受信した画像は、それに関連付けられた推定MZDmax値を有する。受信した医用画像のMZDmax値に対して、表面LUTの一組がステップ750にて計算される。便宜上、新規BUS画像の推定MZDmax値が、パラメータ化層モデルの生成およびキャリブレーションに用いられた離散MZDmax値の1つに一致する場合、対応する表面LUTは単に読み出され、次のステップで直接使用することができる。そうでない場合は、画像のMZDmax値に対応する表面LUTの新たな一組が、パラメータ化層モデルから計算されなければならない。1つの手法では、MZDmax値が一括するまたは画像のMZDmax値に最も近い二組の表面LUTを単純に選択し、表面LUTの新たな一組を、これらの表面LUTの2つの組から、内挿または単純にこれら2つのモデルの間の加重算術平均により計算する。特定のMZDmaxがモデル内の2つの表面により一括されない場合は、表面LUTの新たな一組は、最も類似するMZDmax値を有するモデルから外挿することができる。当然、表面LUTの新たな一組を計算するために、二組以上のLUTを使用するさらに改良された方法を使用することも可能である。計算された表面LUTの新たな一組は、新規CBAMモデルを構成する。
LUTの新たな一組が計算されると、最終段階として推定境界表面、すなわち一次組織層の位置を、計算されたLUTの新たな一組を使用して計算する(ステップ760)。必要な場合は、異なるスキャンパラメータにより生じた偏差を補正、すなわち補償するために、スキャン取得パラメータを使用することもできる。
本処理はCBAMモデルを利用するため、ここではCBAM方法700と言及される。CBAM方法は、例えば任意のBUS画像における一次層の位置推定等、一般的な複数のアプリケーションを有するが、一つのアプリケーションは皮下脂肪領域604を同定することである。第一境界面614と第二境界面616との間のピクセル(図6参照)は、主に脂肪組織を構成すると考えられる。平均脂肪強度は、前述のように、複数の方法の中で、kクラスタリング技術の適用、または単純平均化により、これらのピクセルの強度から抽出することができる。
CBAM方法はいくつかの一般的アプリケーションを有することを理解されたい。例えば、グレーレベル強度は、一次層から一次層の間で著しく変化する傾向がある。撮像された胸部を可視化するため、CBAMモデルから同定された各層は別々にレンダーリングされることができ、これにより、大きく異なる強度により引き起こされる困難を軽減する。または、可視化のために画像をレンダーリングする時に、撮像された胸部の一部のみがレンダーリングされる。例えば、皮膚層または後部***層、またはその両方は、***層のレンダーリングから除外することができる。また、後述するように、病変のより一層の処理および同定は、特定の一次層、もしくは一定タイプの組織または病変がより発生する見込みのある層にフィルタの適用を限定することにより、異なる一次層の位置における知識を利用することができ、フィルタはこれらのタイプの組織または病変を検出するよう設計される。
皮下脂肪に対応する画像領域を同定し(例えばCBAM方法を使用する)、代表脂肪強度を推定したため、図3に戻り、ステップ306で代表脂肪強度を入力画像の正規化に使用する。正規化ステップでは、入力画像の入力範囲[MinGreyLevel, MaxGrayLevel]を、0から2Nrbits−1までに広がるダイナミックレンジにマッピングする。ここで、Nrbitsが出力画像のピクセル毎のビット数を示す。皮下脂肪の推定強度は、出力ダイナミックレンジ(2Nrbits−1)/2の中間強度、例えばNrbits=8の場合127のように、通常は強度の中間レベルに近い点にマッピングされる。正規化は、単なる例であるNrbits=8に限定されないことに注意されたい。
図10は、強度正規化および画像平滑化の前(図10a)と後(図10b)の超音波画像の例であり、画像の上端領域の皮下脂肪領域(1002aおよび1002b)の中間レベルグレーへのマッピングを示す。正規化画像において、低エコー領域と無響領域は、より簡単に同定される。病変はこれらの2つのエコー輝度分類のいずれかに当てはまることが非常に多いため、正規化処理によりそれらの検出と描写が容易になる。
画像強度を正規化した後、次の段階は、病変候補ブロブ検出である。図3を参照すると、様々なタイプの組織エコー輝度に対する自動病変候補検出に対応するために、最初にパラメータマップが生成される。下記では、ステップ314の密度マップ生成、すなわち正規化グレーピクセル値を使用した密度マップの生成を想定する。他のタイプのパラメータマップの生成と処理については後述する。正規化グレーピクセル値は、胸部の組織組成に類似する胸部組織に対するいくつかのエコー輝度領域にクラスタ化される。無響クラス、低エコークラスまたは等エコークラスのいずれの1つに属する領域は、潜在的病変候補として個別に追跡および保存され、後述のように、それらの分類を支援するためのさらなる分析を受ける。
密度マップからブロブを検出するために(ステップ318)、密度マップは初めにクラスタ化されてクラスタリング密度マップが生成される。密度マップからのクラスタリング密度マップの計算は、画像の無響領域、低エコー領域および等エコー領域を検出するためにロバストクラスタリングまたは正規化画像の強度への分類アルゴリズム適用により構成される。この方法に続いて、入力画像のピクセルが初めにピクセル強度に基づいて5つの領域カテゴリにグループ化される。第1クラスタは、画像の無響の暗領域を含む。第2クラスタは、BUS内の低エコー領域をとらえる。第3クラスタは、等エコー脂肪領域を含む。第4クラスタは、わずかに低エコーの腺領域を含み、最後に、皮膚、クーパー靭帯、スペックル雑音および微小石灰化が、高エコーの第5クラスタを構成する。ピクセルをクラスタ化するために、k=5であるk平均クラスタリングアルゴリズムが正規化画像に適用される。このことは、正規化入力画像のダイナミックレンジを、上記の5つのカテゴリに対応する5つの強度クラスタに区分化する。クラスタ内の各隣接領域は、一貫したまたは類似した内部強度特性を有し、そのため「密度ブロブ」である。これらの隣接した個別領域はそれぞれ、潜在的病変として別々に分析することができる。その後、ブロブの検出は、クラスタ化された隣接する領域、すなわち密度ブロブの概略的輪郭の生成のみを必要とする。
真性病変の可能性のあるものを、医用画像に見られる他のタイプの腫瘤から区別するために、ブロブに関連する異なる特徴または特性が抽出され、ブロブを分類するために分析される(ステップ320)。通常、これらの特徴又は特性は、ブロブが何を示すかを説明するため、記述子と称される。以下は、特徴に基づくブロブ分析の1つの手法を説明する。この分析は、最初に一組の関連記述子を選択することにより開始する。次に、これらの記述子は分析され、ブロブの悪性である推定確率との関連性が評価される。これらの記述子を評価し、ブロブが悪性であるという推定確率を生成するために、例えば、CARTツリーを採用することができる。次に、ブロブが悪性である尤度を示す値が計算され、ブロブに割り当てられる。高推定確率を有するブロブは、可能性のある病変としてマークされる。これらの各ステップは後ほど詳述する。
前述のように、これらの特徴、つまり記述子は、CART操作に入力される。CARTアルゴリズムは最初に、これらのブロブ記述子、およびこれらに対応する、トレーニング画像の代表的な一組に対する専門家による分類により展開される。次に、CARTアルゴリズムは入力画像から同定された各ブロブの一組の記述子に適用される。CARTツリー操作の出力は、ブロブが病変であることの尤度値または推定確率を取得するために使用される。偽陽性は、分析されたブロブに割り当てられた尤度値に基づいて除去される。(ステップ322)。
CART操作を準備するために、まず、固体小結節と疑わしくない病変候補とを区別するために、いくつかのブロブ特徴が定義される。特徴又は記述子は、形状、グレーレベル変動および空間的位置の3つの主なカテゴリに分類することができる。記述子の各カテゴリは、さらにサブクラスに分類することができる。例えば、形状記述子は、セグメント化されたブロブ候補から生じた特徴と、楕円をブロブの輪郭に合わせた結果生じる特徴との2つのサブクラスに分類することができる。稠密度指標と伸長値は、第1サブクラスに属する。稠密度指標は、下記の様に計算することができる。
ここで、円は稠密度1を有し、正方形は稠密度値π/4を有する。上記の数式では、BlobAreaはブロブ領域であり、BlobPerimeterはブロブの概略的輪郭の全長を示す。
伸長指標は、幅と高さの割合を使用して定義され、下記の様に計算することができる。
上記の数式において、WidthBlobBoundingBoxは密接にブロブの境界をつける長方形ボックスの幅であり、HeightBlobBoundingBoxは長方形境界ボックスの高さである。伸長度の値は常に0よりも大きい。1の値は、ほぼ正方形または円形のオブジェクトを表す。伸長度の値が無限である傾向に従って、オブジェクトは水平方向により伸長される一方で、その伸長度の値が0に近づくに従ってオブジェクトは垂直方向により延伸される。
同様に、楕円をブロブの輪郭に合わせた結果生成された特徴は、主軸の偏心および配向の2つのサブクラスにさらに分けられる。偏心は、下記の様に定義される。
上記の数式において、ShortEllipseAxisLengthは合わせた楕円の短軸の長さであり、LongEllipseAxisLengthは合わせた楕円の長軸の長さである。偏心値は、厳密に0より大きく、1以下である。
主軸の配向は、以下から計算される。
μ11、μ20、μ02はオブジェクト内でピクセルが腫瘤の中心からどれだけ分散しているかを測定する二次モーメントである。より一般的には、中央モーメントμmnは下記の様に定義される。
ここで、(Xmean, Ymean)は腫瘤の中央の座標である。
グレーレベル変動記述子も、病変内のグレーレベル変動を記述するカテゴリと、病変の局部的および全体的背景に対して病変グレーレベルコントラストを記述するカテゴリとの、2つのカテゴリに分けられる。病変内のピクセルのグレーレベル変動を記述する特徴は、さらに分割され下記の4つのサブカテゴリに分類される。
a.二次中心モーメントである、分散
ここで、ブロブはNピクセルを含み、ブロブ内のithピクセルのグレーレベルはxiにより示され、μはブロブ内の全てのピクセルの平均グレーレベルを示す。
b.三次中心モーメントであり確率変数の非対称を表す、歪度
ここで、σは標準偏差または変数の平方根である。負の歪度は左に歪んだデータを示し、正の歪度は右に歪んだデータを示す。
c.第4モーメントである、尖度
正の尖度は「尖頂を有する」分布を示し、負の尖度は「平らな」分布を示す。
d.セグメント化された病変に含まれる全てのピクセルのL2勾配ノルム
病変の背景と皮下脂肪領域に対して病変内のピクセルのグレーレベルの変動を表す特徴は、下記の2つのサブカテゴリに分類することができる。
a.任意の背景上の病変の可視性。背景は、ブロブの境界輪郭の外縁から始まる領域として定義され、輪郭を超えていくつかのピクセルを伸長する。ブロブ領域を囲むこのような背景領域を描出する1つの方法は、ブロブ領域を示す二値画像に形態学的膨張フィルタを適用し、その後、元のブロブ領域を膨張領域から差引いて、ブロブ領域を直接包囲する管状背景領域を残すことである。さらに単純な背景領域は、ブロブ領域と境を接する長方形のボックスを当てて、背景領域を示すために、ブロブ領域外で当てられた長方形内の全てのピクセルを考慮することにより抽出することができる。可視性指標は、以下から計算することができる。
この数式において、MeanLesionは病変内のピクセルの平均グレー値であり、MeanBackgroundは背景領域のピクセルの平均グレー値である。
b.病変内のピクセルの平均グレーピクセル値の正規化値と皮下領域からの脂肪値:
上述のように画像から計算または抽出された特徴の分析および偽陽性除去の後、最終ステップ(ステップ324)にて、閾値を超える尤度値を有する全てのブロブが、例えば、放射線医によるより一層の検討のために表示装置に表示されたり、またはより一層の自動分析のために追加CADモジュールに転送されることにより、報告される。これらは尤度値の降順に表示されるよう分類されてから後に報告することができる。これにより、自動病変検出処理が完了する。
病変候補の自動検出のパフォーマンス向上のため、上記の処理にバリエーションを取り入れることができる。例えば、上記の説明は、概して脂肪正規化画像の処理に関連している。通常、入力画像の正規化は、皮下脂肪のみの代わりに、いくつかのコントロールポイント組織を含むことができる。このことは、脂肪の強度に加え、入力画像強度範囲内にいくつかの基準強度値(「コントロールポイント値」)を確立することを要する。各コントロールポイントは、入力またはノイズ除去された画像から測定されたコントロールポイント組織の代表強度と、コントロールポイント組織の期待される、または割り当てられた強度との2つの値とを含む。ステップ308(処理500として詳述されるステップ)で判定される脂肪の強度と他のコントロールポイントの強度は、画像正規化参照テーブルを生成するために使用される。
一実施例において、ステップ130(図3参照)で、正規化LUTに対するコントロールポイント値を計算するために、ロバストクラスタリング操作(例えば、k平均クラスタリング)が、皮下脂肪領域のみの代わりに全体画像に適用される。最初のクラスタリング実行で取得された脂肪および高エコーの強度値は、全体画像を分類する一方、脂肪および皮膚の後続するクラスタ中心を設定するために使用される。この第2クラスタリング操作の結果、無響、低エコー輝度、および高エコー輝度を示す強度をとらえるいくつかの追加の特異的なクラスタ中心が推定される。これらの幅広いクラス(例えば高エコー輝度)は、クラスタリングにより個別に検出される必要があるエコー輝度のいくつかのサブクラスを含むことが多い。例として、類似しているが統計的に異なる、皮膚と線維腺組織のエコー輝度が挙げられる。スペックルノイズのレベルを減らすために、平滑化画像を使用することができる。あるいは、クラスタリングアルゴリズムを適用する前に、読み出した画像にメジアンフィルターを適用することができる。
全てのコントロールポイントの強度が推定または測定された後、前処理画像のピクセルの強度と、結果として生じる画像との間の一貫性のあるマッピングを確立する必要がある。マッピング関係(ステップ312)は、様々な方法で確立することができる。例えば、コントロールポイントは所定の値または割当てられた値にマッピングされ、隣接するコントロールポイントの間のピクセル値が補間される。また、最小または最大ピクセル値とその隣接するコントロールポイントとの間の強度を有するピクセル値を、補間することができる。例えば、スプライン当てはめのような、任意の周知の曲線当てはめ方法を使用して、最小または最大ピクセル値を含む滑らかな曲線をこれらのコントロールポイントに当てはめ、これらのコントロールポイントの間に、値を補間することができる。
図11は、測定されたコントロールポイント組織の、それらのそれぞれ割り当てられたグレーピクセル強度値へのマッピングを表す。本実施例において、これらのコントロールポイント1104(そのうちの1つは脂肪1106)を連結する滑らかな曲線1102が最初に求める。例えば、スプライン補間はいくつかのコントロールポイントを入力し、コントロールポイントを連結するように滑らかな曲線を当てはめる。次に、滑らかな曲線1102に基づく参照テーブルが生成され、すなわち、滑らかな曲線1102に示される近似関数を使用して計算され、入力画像の最初のピクセル値1108から正規化画像の出力ピクセル値1110までの速いマッピングを容易にする。
マッピングLUTまたは正規化LUTにおける各コントロールポイント位置は、入力画像の基準強度と、あらかじめ確立された相対的エコー輝度関係に密接に従う各基準組織に対する、所定または割り当てられた出力強度との関数である。このような相対エコー輝度関係の1つは、A.T.Stavrosによる"Breast Ultrasound", Lippincott Williams and Wilkins, 2004により提案されている。この方法を説明するために、下記では一組のコントロールポイントの割当てを説明する。
1.最小入力グレーピクセル値は0にマッピングされる。
2.無響領域の強度は出力強度PCyst*(2NrBits-1)にマッピングされ、PCystは最大出力強度の所定の割合値(通常5%)である。
3.皮下脂肪の推定強度は、出力ダイナミックレンジ(2NrBits-1)/2の中強度、例えばNrBits=8である場合127に位置づけされる。
4.皮膚強度は、入力画像の高エコーピクセル強度の範囲の底部で認識され、新規強度PSkin*(2NrBits-1)にマッピングされ、PSkinは90%など所定の高い割合値である。
5.線維腺組織のエコー輝度の強度は、入力画像の高エコーピクセル強度の中程度で同定され、出力強度PFibroglandular*(2NrBits-1)に割り当てられ、PFibroglandularはPSkinよりも大きい所定の割合値で、例えば95%である。
6.カルシウム強度は、入力画像の最高強度クラスタにおける平均グレーピクセル値であると推定される。これらの入力強度はPCalcium*(2NrBits-1)に割り当てられ、PCalciumはPFibroglandularよりもさらに大きい所定の割合値であり、例えば98%である。
7.最後に、最大入力グレーピクセル値は(2NrBits-1)に割り当てられ、例えばNrBits=8の場合255である。
NrBitsという用語は、入力ピクセル強度値を示すために使用されるビット数を示し、本例では通常の値8が示され、入力ピクセル強度のダイナミックレンジが0から255となる。より大きなデータタイプ(NrBitsの値が16、32または64)、または非整数強度値をサポートする浮動小数点タイプもまた、これらの目的のため使用することができる。
1.最小入力グレーピクセル値は0にマッピングされる。
2.無響領域の強度は出力強度PCyst*(2NrBits-1)にマッピングされ、PCystは最大出力強度の所定の割合値(通常5%)である。
3.皮下脂肪の推定強度は、出力ダイナミックレンジ(2NrBits-1)/2の中強度、例えばNrBits=8である場合127に位置づけされる。
4.皮膚強度は、入力画像の高エコーピクセル強度の範囲の底部で認識され、新規強度PSkin*(2NrBits-1)にマッピングされ、PSkinは90%など所定の高い割合値である。
5.線維腺組織のエコー輝度の強度は、入力画像の高エコーピクセル強度の中程度で同定され、出力強度PFibroglandular*(2NrBits-1)に割り当てられ、PFibroglandularはPSkinよりも大きい所定の割合値で、例えば95%である。
6.カルシウム強度は、入力画像の最高強度クラスタにおける平均グレーピクセル値であると推定される。これらの入力強度はPCalcium*(2NrBits-1)に割り当てられ、PCalciumはPFibroglandularよりもさらに大きい所定の割合値であり、例えば98%である。
7.最後に、最大入力グレーピクセル値は(2NrBits-1)に割り当てられ、例えばNrBits=8の場合255である。
NrBitsという用語は、入力ピクセル強度値を示すために使用されるビット数を示し、本例では通常の値8が示され、入力ピクセル強度のダイナミックレンジが0から255となる。より大きなデータタイプ(NrBitsの値が16、32または64)、または非整数強度値をサポートする浮動小数点タイプもまた、これらの目的のため使用することができる。
図11は、上記の各コントロールポイント組織へのピクセル強度値のグレーレベル割り当てを示す。本例は一組の7つのコントロールポイントを説明するが、当然、他の適切な組のコントロールポイントを選択および定義することも可能である。コントロールポイントの選択は、画像診断法(例えば、MRI画像は異なる一組のコントロールポイントを必要とする)および撮像された解剖領域(例えば、肺、前立腺、または卵巣の画像は異なる一組のコントロールポイントを使用することがより好ましく標準化される)に左右されることが多い。
正規化LUTは、ステップ306にて、ノイズ除去画像の正規化に使用することができる。その結果、全般強度正規化画像が生成される。強度正規化画像を処理するさらなるステップは、脂肪正規化画像の処理に関する上記の処理、つまり、ステップ314、318乃至324と基本的に同様である。これらのさらなるステップの説明はここでは繰り返さない。
他の変形態様は、ハードウエアおよびオペレータパラメータにおける変動を補償するために使用できる悪性マップの使用に関する。前述の様に、画像処理に関与する様々な取得パラメータは、画像データの一貫性に影響を及ぼす。これらのパラメータは2つのグループに分類することができる。第1クラスは、異なるメーカーの超音波トランスデューサ機材における変動によるパラメータであり、深度およびトランスデューサ周波数を含む。第2クラスは、トランスデューサ圧力およびTGC等、技術者による手動設定に関連する要素を含む。
ステップ316では、悪性可能性マップが生成される。このステップは、密度マップの生成(ステップ314)の代替または追加となりえる。悪性可能性マップは、入力画像の各ピクセルに、ピクセルの悪性である確率値を割り当てる。確率値は良性の0.0と悪性の1.0の間に広がる。
当業者にとって周知のように、確率は0から1の間の任意の値を取ることができる。一方、専門家によるマークはバイナリ方式で病変領域を示す。つまり、病変領域内のピクセルは1の値を有し、画像背景の悪性値は0に設定される。バイナリ変数0と1を扱うためのモデルを生成するために、ロジスティック回帰が使用される。モデルは放射線医によりマークされた病変領域を含む多くの画像でトレーニングされる。ロジスティックモデルは、下記のような画像ピクセル値、ハードウエアおよびオペレータ変数を組み込んで生成される。
1.正規化グレーピクセル値
2.クラスタ化された密度マップ値
3.皮膚表面からの検査領域の深度を示すmmで表されるピクセルサイズ
4.トランスデューサ周波数
5.TGC値
1.正規化グレーピクセル値
2.クラスタ化された密度マップ値
3.皮膚表面からの検査領域の深度を示すmmで表されるピクセルサイズ
4.トランスデューサ周波数
5.TGC値
従って、悪性確率マップは、正規化グレーピクセル値、密度マップ値および取得パラメータを考慮して生成されるため、これらのパラメータの不一致を最小限に抑える。
ロジスティックモデルの1つの主な利点は、皮膚表面からの物理領域深度を考慮することであり、その結果得られる悪性確率マップは、後部特徴として陰影または部分的な陰影を有する病変の上部を示すことができることである。一方、通常は深度に関する情報を持たない密度マップは、後部特徴を示すことができない。従って、密度マップのみの検査では見逃される特定の病変候補を、悪性可能性マップから特定することができる。
悪性確率マップは、所定の確率閾値を確率マップに適用することにより、ステップ318にてクラスタ化される。0.75等の所定の閾値よりも大きい確率値を有する全ての領域は、ブロブ、つまり「悪性ブロブ」に分類される。検出された悪性ブロブの分析およびそれに関する報告を行うさらなるステップは、密度ブロブの分析および報告(ステップ320および322)に関連する説明と同様であるため、ここでは説明を繰り返さない。
悪性確率マップの生成は、上記のロジスティックモデル手法に限定されない。以下は、悪性確率マップを生成する他の方法を詳細に説明する。この方法は、ピクセルを異なる組織タイプに分類するために、ピクセルに特有の特徴およびピクセルの近隣に関連する特徴を分析する。ピクセルが一組の組織タイプのそれぞれに属する確率値は、ピクセルが病変に属する確率を含めて、ピクセルに割り当てられ、そこから悪性確率マップが生成される。この方法は、後節で詳細に説明する。
図12を参照すると、まず、入力画像が前処理される(ステップ1210)。これはノイズ除去(サブステップ1212)、エッジ保存フィルタリング(サブステップ1214)および画像正規化(サブステップ1216)を含むが、必ずしもこの順序とは限らない。入力画像の強度は、動的に検出されたコントロールポイント組織強度(皮下脂肪等)を使用して正規化される(サブステップ1216)。異方性エッジ保存フィルタリング(サブステップ1214)が、超音波画像から典型的なスペックルノイズを除去するために適用される。フィルターパラメータは、特定のスキャナの低いネイティブ解像度、または他のスキャナに適用される増強された「固有平滑化」を説明するメーカー特有の画像特性に対応するためにチューニングされる。他の前処理ステップを含むこともできる。例えば、フィルタチューニングは、前処理に対するフィルタ構成に影響することが多いトランスデューサ周波数等の要素を含むことができる。
次に、この方法は、ベクター、すなわちピクセルおよびその近隣を示す一組のパラメータの計算または構築に関与する(ステップ1220)。一組のパラメータは、ピクセル特性ベクター(PCV)と称される。当然、3D画像の場合、一組のパラメータはボクセルを表し、ボクセル特性ベクター(VCV)と称される。ここで説明する方法は、2Dまたは3Dに関わらず、等しく適用可能である。下記において、PCVとVCVを区別しない。両者とも、PCVと称される。ピクセル固有特性値は、例えば、ピクセルの正規化グレーレベルと、皮膚ラインからのピクセルの物理的深度とを含む。入手できる限り、乳頭からのピクセルの垂直距離や、乳頭(胸部クロック位置)を横断するSI(上下位)ラインからのピクセル位置の角度のようなピクセルの概位を含むこともできる。これらのピクセル固有特性は、サブステップ1222にて画像データから抽出され、その後ピクセルのPCVに含められる。
各ピクセルの近隣の特性もまた、多重解像度ピクセル近隣解析で測定される(サブステップ1224)。いくつかの近隣サイズおよびスケールは分析に使用され、分析は各スケールの各ピクセルの物理的サイズ(例えば、mm/ピクセル)を具体的に説明する。いくつかの多重解析フィルタが、ピクセルの近隣に適用される。多重解析フィルタに対する反応は、近隣における、質感、エッジ構造、線状パターンまたはグレーレベル強度等の、対象ピクセルの近隣特性を提供する。これらのフィルタの反応は、カテゴリに分類され、PCVの近隣部に含められる。以下の例は、組織タイプの同定を支援するために適用可能なフィルタータイプを説明する。
1.線状パターンフィルタは、画像内の線形構造に対する高反応を実現する。胸筋は、その特徴的短い線状「小体」により同定されることが多いため、適切なスティックフィルタが高反応を起こす独特の質感を与える。
2.多重スケールで高反応を起こすエッジフィルタは、長い構造を示すことが多い。画像の上半分における長エコー輝度ラインは、画像上部において皮膚を示す可能性のある水平状でない限り、乳筋膜またはクーパー靭帯に一致することがある。
3.直径4mm未満の円形低エコー領域は、通常、健康な管または終末乳管小葉単位(TDLUs)を示し、テンプレート円形コンボリューションカーネルは、これらの領域に対する高反応を起こす傾向がある。
4.大きな近隣サイズにおけるフーリエ分析は、シャープなエッジの細部の量に対して徐々に変化する領域内のシェーディングを示し、これは脂肪の見込みがある、低周波数の質感変動を有する、等エコー組織の領域を同定するために使用することができる。
次に、ピクセル特性ベクター(PCV)が、サブステップ1222と1224で獲得されたピクセル固有値とピクセル近隣値から各ピクセルに対して構成される(サブステップ1226)。
各ピクセルに対するPCVを獲得した後、次のステップ(1230)では、各ピクセルのPCVを分類する。PCVを分類するために、任意の適切な分類子を使用することができる。一般に、多次元分類コンポーネントはPCVを入力とし、問題となっているピクセルが各特定出力クラスに属する確率を表す出力ベクターを生成する。一組の出力クラスは、以下のうち1つ以上を含む。
・脂肪組織
・管
・TDLU
・筋膜
・クーパー靭帯
・病変組織
・胸筋(筋)組織
・リンパ節
・肋骨
・脂肪組織
・管
・TDLU
・筋膜
・クーパー靭帯
・病変組織
・胸筋(筋)組織
・リンパ節
・肋骨
一実施例において、PCVを分類するためのCARTツリーは、多次元分類子を構成、すなわちキャリブレーションするために専門家にマークされたトレーニングデータのセットを使用して生成される。各出力タイプに対する最大識別反応を発生させるようにフィルタを整調するために、PCVを獲得するために使用されるフィルタセットと類似した、もしくは同一の多重解像度フィルタセットが、これらのBUS画像を背景に適用される。CARTツリーでのPCVの処理は、PCVが可能性のある出力カテゴリのいずれかひとつ、特に、病変組織カテゴリに当てはまる確率を返す。各ピクセル、すなわち各PCVが分類された後(ステップ1230)、BUS画像内の各ピクセルが属する組織および解剖タイプを表すBAMが生成される。
これらの2つのステップ1220と1230およびそれらのサブステップ、すなわち各ピクセルに対する一組のPCVを抽出するために、一組の多重解像度フィルタをBUS画像に適用することと、PCVを分類するために続けて多次元分類子を適用することは、BAM部218’で実行されることができる。次に、BAM部218’は、悪性確率マップを抽出または生成するために、分類されたPCVを悪性マップ部218に渡すことができる。
病変組織カテゴリは、ピクセルが疑わしく、さらなる調査を必要とするBUS画像の領域に属することを示す。PCV分類時に、病変組織に属する対応するピクセルの確率が取得されて、ピクセルに割り当てられる。悪性確率値をピクセルにマッピングすることにより、悪性確率マップが生成される(ステップ1240)。
ニューラルネットワークと混合モデル(クラスタに基づく)技術を含む、別の分類システムを使用することも可能である。これらの他の分類子内の内部処理は異なるが、全ての手法は同一タイプの入力を取り、ピクセルを分類し、同一タイプの出力を生成すると考えることができる。
図3に示す処理の、各ステップの代替手段により紹介された変形に加え、ここに記載した異なる代替手段を組み合わせることにより、図3に示す処理の変形も可能である。図13は図3に示す例の1つの代替手段の例を示す。
この変形形態によると、前述のように、その取得パラメータと共に画像は最初に読み出しおよびノイズ除去される(ステップ1302)。次に、エッジ保存フィルタリングが、ノイズ除去画像に適用される(ステップ1304)。一方、MZDmaxでパラメータ化されたCBAMモデルが生成される(ステップ1306)。CBAMモデル構築の詳細なステップは、図6から9を参照に既に説明済みであるため、ここでは繰り返さない。画像の推定MZDmax値に基づいて、画像の新規CBAMモデルが生成され、そこから一次層境界面の位置を推定することができる(ステップ1308)。次に、この画像は、脂肪のみの代表強度またはいくつかのコントロールポイント組織の代表強度を使用して正規化される(ステップ1310)。
前述のように、CBAMモデル方法の出力により提供された一次層の位置に対する知識は、後に続くステップの精度と効果を改善するために使用することができる。特定タイプの組織または病変に反応するフィルタは、このタイプの組織または病変が発生すると最も考えられる一次層のみにおいて選択的に適用可能であり、それらの存在が期待されない1つまたは複数の層には適用されない。通常、例えば、胸筋は、1mm2よりも大きな近隣で明白に可視できる特徴的質感を有するが、通常は後方***領域のみに予期される。効率を上げ偽陽性を低減するために、胸筋を検出するよう設計されたフィルタを、後方***領域で同定されるピクセルのみに適用することができる。フィルタ反応は、他の全ての一次層に対して0に設定される。このことは、皮膚層の胸筋を誤って報告する可能性をなくす。同様に、***領域にもっともよく起こる病変を検出するよう設計されたフィルタは、他の層ではなく、***領域内のピクセルのみに適用することができる。
次のステップでは、画像内の各ピクセルに対してPCVを構築(計算)する時(ステップ1312)、上記の様に効率と精度を改善するために、***領域606内のみのピクセルに計算を限定することができる。次に、***領域のピクセルは、計算された各タイプに属する確率値を有する異なる組織タイプに分類される(ステップ1314)。ピクセルの悪性である確率値から、悪性マップが生成される(ステップ1316)。次のステップでは、例えば悪性確率の閾値を悪性マップに適用することにより、ブロブを分離する(ステップ1318)。さらなる分析ステップ、すなわち、ブロブ検出(ステップ1320)、ブロブ分析(ステップ1322)および偽陽性の除去(ステップ1324)を実行することができる。最後に、病変の候補が報告される(ステップ1326)。これらのステップは、全て詳細に説明済みであるため、ここではその説明を繰り返さない。
本発明の様々な実施例を、詳細に説明した。当業者であれば、本発明の範囲から逸脱することなく、これらの実施例に対し種々の変更、適応および変形を行うことが可能であることを理解されたい。上記のベストモードにおける変更または追加を、本発明の本質、精神または範囲から逸脱することなく行うことができるため、本発明はこれらの詳細にではなく、添付された特許請求の範囲によってのみ限定される。
Claims (30)
- 超音波医用画像内における疑われる病変を同定する方法において、
前記医用画像内の皮下脂肪ピクセルの推定代表脂肪強度値を計算するステップと、
正規化画像を取得するために、正規化脂肪強度値と前記代表脂肪強度値の間のマッピング関係を利用して、前記医用画像のピクセル値から正規化グレーピクセル値を計算するステップと、
複数の個別領域を形成する前記正規化画像内のピクセルを同定し、前記個別領域の各々は一貫性のある内部特性を有するステップと、
前記各個別領域から、記述的特徴を抽出するステップと、
前記各個別領域の前記抽出された記述的特徴を分析し、前記各個別領域が病変である尤度値を前記各個別領域に割当てるステップと、
所定の基準を満たす尤度値を有する全ての個別領域を病変候補として同定するステップと、
を備えることを特徴とする方法。 - 前記医用画像の正規化グレーピクセル値を計算する前に、ノイズ除去画像を得るために前記医用画像のノイズを除去するステップをさらに備え、前記正規化グレーピクセル値は、前記ノイズ除去画像のピクセル値から計算されることを特徴とする、請求項1記載の方法。
- 前記ノイズを除去するステップは、エッジ保存拡散画像フィルタリングを前記医用画像に適用することをさらに備えることを特徴とする、請求項2記載の方法。
- 前記ノイズを除去するステップは、エッジの鮮明さを保ちイントラリージョン平滑化を促進する選択的フィルタリングを前記医用画像に適用することをさらに備えることを特徴とする、請求項2記載の方法。
- 前記超音波画像内の皮下脂肪領域を選択するステップをさらに備え、
前記推定代表脂肪強度値は、前記選択された皮下脂肪領域内のピクセルから計算されることを特徴とする、請求項1記載の方法。 - 前記推定代表脂肪強度値を計算するステップは、k平均クラスタリングを前記選択された皮下脂肪領域に適用することさらに備えることを特徴とする、請求項5記載の方法。
- 前記推定代表脂肪強度値を計算する前記ステップは、前記k平均クラスタリングから取得した中間レベル強度クラスタ内のピクセルの平均ピクセル強度値を計算することをさらに備えることを特徴とする、請求項6記載の方法。
- 皮下脂肪を含む複数のコントロールポイント組織を選択するステップと、
強度クラスタを生成するためにピクセル強度をクラスタリングし、前記複数のコントロールポイント組織はそれぞれ少なくとも1つの前記強度クラスタにより示されるステップと、
前記複数のコントロールポイント組織の各々に対して、
少なくとも1つの強度クラスタに対応するピクセルのグレー強度値から、前記各コントロールポイント組織に対する代表強度値を判定するステップと、
正規化強度値を前記各コントロールポイント組織に割り当てるステップと、をさらに備え、
前記マッピング関係は、前記正規化強度値を前記複数のコントロールポイント組織のそれぞれの代表強度値に関連付けることを特徴とする、請求項1記載の方法。 - 悪性確率マップを生成するために、前記正規化画像の各ピクセルに対する悪性確率値を推定することをさらに備え、
前記個別領域を同定するステップは、
前記個別領域を形成するために、閾値を上回る悪性確率値を有するピクセルを隣接領域に分類することを含むことを特徴とする、請求項1記載の方法。 - 前記悪性確率値を推定するステップは、
サンプル画像のコレクションでキャリブレーションされたロジスティックモデルを確立し、前記ロジスティックモデルは、画像ピクセル値と、ハードウエアおよびオペレータ取得パラメータとを取り組むステップと、
前記悪性確率値を推定するために、前記ロジスティックモデルを前記正規化画像に適用するステップと、
をさらに備えることを特徴とする、請求項9記載の方法。 - 病変に対応する可能性のある医用画像における領域を自動的に同定するシステムにおいて、
前記医用画像内のコントロールポイント組織の推定強度を、前記医用画像のピクセル値から計算するよう構成される強度部と、
入力ピクセルと正規化ピクセルとの間のマッピング関係を生成し、グレーピクセル値を正規化ピクセル値に変換して、前記マッピング関係に従って正規化画像を取得するよう構成される正規化部と、
パラメータ値を入力画像の各ピクセルに割り当てて、パラメータマップを生成するマップ生成部と、
前記パラメータマップにあるブロブを、検出および画定するよう構成されるブロブ検出部と、
前記検出ブロブの記述的特徴を、検出および計算するよう構成される特徴抽出部と、
ブロブが悪性である推定尤度値をブロブの記述的特徴から計算し、前記尤度値を前記ブロブに割り当てるブロブ分析部と、
を備えることを特徴とするシステム。 - 前記コントロールポイント組織は、皮下脂肪、無響嚢胞、皮膚、および線維腺組織の少なくとも1つを含むことを特徴とする、請求項11記載のシステム。
- 前記マップ生成部は密度マップ部を含み、前記パラメータ値は密度であり、前記パラメータマップは密度マップであることを特徴とする、請求項11記載のシステム。
- 前記マップ生成部は悪性確率マップ部をさらに含み、前記悪性確率マップ部は、各ピクセルに対する悪性確率を計算し、前記悪性確率を前記各ピクセルに割り当てて悪性確率マップを生成し、前記検出部により検出および画定された前記悪性確率マップ内のブロブは、前記特徴抽出部および前記ブロブ分析部により処理された前記ブロブに含まれることを特徴とする、請求項11記載のシステム。
- 前処理部をさらに備え、前記前処理部は、前記医用画像のノイズ除去を行いノイズ除去画像を生成することを特徴とする、請求項11記載のシステム。
- 報告部をさらに備え、前記報告部は、所定の閾値を上回る尤度値を有する全てのブロブを同定することを特徴とする、請求項11記載のシステム。
- デジタルまたはデジタル化された医用画像における組織のグレースケール強度を推定する方法において、
クラスタリング操作を医用画像のピクセルの強度値に適用して、前記強度値を個別強度クラスタに分類するステップと、
前記個別強度クラスタの1つを、前記組織に対応する強度クラスタとして同定するステップと、
前記強度クラスタの代表グレースケール強度値を、前記強度クラスタのピクセルのグレースケール強度から推定するステップと、
前記代表グレースケール強度値を前記組織に割り当てるステップと、
を備えることを特徴とする方法。 - 前記代表グレースケール強度値は、前記グレースケール強度の平均強度であることを特徴とする、請求項17記載の方法。
- 前記組織は皮下脂肪であることを特徴とする、請求項17記載の方法。
- 前記クラスタリングアルゴリズムは、k=3を有するk平均アルゴリズムであり、前記同定された個別クラスタは中強度クラスタであることを特徴とする、請求項17記載の方法。
- 前記組織が位置することが見込まれる医用画像内の領域を画定することをさらに備え、
前記クラスタリング操作を前記画定された領域に適用することを特徴とする、請求項17記載の方法。 - 胸部超音波画像を処理する方法において、
胸部の層モデルを構築し、前記モデルの隣接層の各対は前記隣接層の各対間の境界面を定義するステップと、
複数のサンプル超音波胸部画像でモデルをキャリブレーションし、前記複数のサンプル超音波胸部画像はそれぞれ、前記サンプル超音波胸部画像の境界面を同定するために手動でセグメント化され、キャリブレーションされた前記モデルはパラメータ化表面モデルを含み、各パラメータ化表面モデルは、サイズパラメータの離散値に対応する一組の境界面参照テーブル(LUT)を含むステップと、
前記超音波胸部画像の前記サイズパラメータの推定値を受信するステップと、
前記サイズパラメータの前記推定値に対応する新規表面モデルを、前記パラメータ化された表面モデルから計算し、前記新規表面モデルは、前記サイズパラメータの前記推定値に対応する一組の計算された境界面LUTを含むステップと、
前記超音波胸部画像内の一次層のピクセルを同定するために、境界面の推定位置を、前記新規表面モデルの前記一組の計算された境界面LUTから計算するステップと、
を備えることを特徴とする方法。 - 複数のコントロールポイント組織の代表強度値を獲得し、前記複数のコントロールポイント組織の前記代表強度値は、皮下脂肪層のグレーピクセル値の代表強度値を含むステップと、
前記複数のコントロールポイント組織の前記代表強度値と前記複数のコントロールポイント組織の正規化値との間のマッピング関係に基づいて、前記複数のコントロールポイント組織に対して前記超音波胸部画像を正規化することにより正規化画像を得るステップと、
をさらに備えることを特徴とする、請求項22記載の方法。 - 前記正規化画像を可視化するためにレンダーリングするステップをさらに備えることを特徴とする、請求項23記載の方法。
- 前記層モデルは皮膚層と後部***層とを含み、前記方法は、
前記皮膚層を画定する第1境界面と、前記後部***層を画定する第2境界面の位置を推定するステップと、
前記皮膚層と前記後部***層とを除外する、前記超音波胸部画像の一部を可視化するためにレンダーリングするステップと、
をさらに備えることを特徴とする、請求項22記載の方法。 - 前記一次層は***領域であり、前記方法は、
前記***領域の各ピクセルに対して、各ピクセルの悪性である確率値を推定するステップと、
閾値を上回る確率値を有するピクセルを、疑わしい病変として隣接領域に分類するステップ、
をさらに備えることを特徴とする、請求項22記載の方法。 - 胸部超音波画像内の病変を同定する方法において、
一次層組織を分離する表面の推定位置を計算し、前記一次層組織は***領域の組織を含むステップと、
前記***領域のピクセルを同定するステップと、
前記***領域の各ピクセルに対してピクセル特性ベクター(PCV)を構築し、前記PCVは、少なくとも前記ピクセルの各々の近隣の特性を含むステップと、
前記***領域の各ピクセルに対して、前記各ピクセルのPCVから悪性確率値を計算し、前記ピクセルのそれぞれに悪性確率値を割り当てて、前記割り当てられた悪性確率値が閾値を上回る場合は、ピクセルを病変の可能性のあるピクセルとして同定するステップと、
全ての病変の可能性のあるピクセルの隣接領域を、潜在病変として報告するステップと、
を備えることを特徴とする方法。 - 皮下脂肪を含む複数のコントロールポイント組織の代表強度値を計算するステップと、
割当てられた強度値と前記複数のコントロールポイント組織の各代表強度値との間のマッピング関係に従って、前記***領域のピクセルのグレーピクセル値を正規化するステップと、
をさらに備え、
前記ピクセルに対する前記PCVは、前記***領域のピクセルのグレーピクセル値から構築されることを特徴とする、請求項27記載の方法。 - 前記各ピクセルに対する前記PCVは、前記各ピクセルのみから判定されたピクセル固有特性を含むことを特徴とする、請求項27記載の方法。
- 手動でマークされた胸部画像のコレクションで分類子をキャリブレーションし、前記手動でマークされた各胸部画像はマークされた病変ピクセルを含み、前記マークされた病変ピクセルはそれぞれPCVを有し、前記分類子の出力は病変組織クラスを含むステップ、
をさらに備え、
前記PCVから前記悪性確率値を計算するステップは、悪性確率値を取得するために、前記分類子を前記PCVに適用することを含むこと特徴とする、請求項27記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13972308P | 2008-12-22 | 2008-12-22 | |
US61/139,723 | 2008-12-22 | ||
PCT/CA2009/001881 WO2010071999A1 (en) | 2008-12-22 | 2009-12-22 | Method and system of automated detection of lesions in medical images |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2012512672A true JP2012512672A (ja) | 2012-06-07 |
Family
ID=42266183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011541052A Pending JP2012512672A (ja) | 2008-12-22 | 2009-12-22 | 医用画像内病変自動検出方法およびシステム |
Country Status (7)
Country | Link |
---|---|
US (2) | US20100158332A1 (ja) |
EP (2) | EP2378978B1 (ja) |
JP (1) | JP2012512672A (ja) |
CN (3) | CN102438529B (ja) |
CA (1) | CA2783867A1 (ja) |
HK (1) | HK1199127A1 (ja) |
WO (1) | WO2010071999A1 (ja) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140057130A (ko) * | 2012-10-25 | 2014-05-12 | 삼성전자주식회사 | 이미지 처리 시스템 |
WO2016051764A1 (ja) * | 2014-09-29 | 2016-04-07 | 富士フイルム株式会社 | 光音響画像生成装置 |
JP2016187748A (ja) * | 2016-08-12 | 2016-11-04 | コニカミノルタ株式会社 | 超音波画像診断装置及び超音波画像診断方法 |
US10383602B2 (en) | 2014-03-18 | 2019-08-20 | Samsung Electronics Co., Ltd. | Apparatus and method for visualizing anatomical elements in a medical image |
KR20200092466A (ko) * | 2019-01-07 | 2020-08-04 | 재단법인대구경북과학기술원 | 의료 영상의 분석 모델을 학습시키는 학습 장치 및 그 학습 방법 |
JP2021118840A (ja) * | 2020-01-30 | 2021-08-12 | 株式会社日立製作所 | 医用画像処理装置及び医用画像処理方法 |
JP2021527473A (ja) * | 2018-06-14 | 2021-10-14 | ケイロン メディカル テクノロジーズ リミテッド | 即時精密検査 |
KR20220040507A (ko) * | 2017-05-11 | 2022-03-30 | 베라톤 인코포레이티드 | 확률 맵 기반의 초음파 스캐닝 |
US11348228B2 (en) | 2017-06-26 | 2022-05-31 | The Research Foundation For The State University Of New York | System, method, and computer-accessible medium for virtual pancreatography |
Families Citing this family (117)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8238637B2 (en) * | 2006-10-25 | 2012-08-07 | Siemens Computer Aided Diagnosis Ltd. | Computer-aided diagnosis of malignancies of suspect regions and false positives in images |
US8077927B1 (en) * | 2006-11-17 | 2011-12-13 | Corelogic Real Estate Solutions, Llc | Updating a database with determined change identifiers |
JP5523891B2 (ja) * | 2009-09-30 | 2014-06-18 | 富士フイルム株式会社 | 病変領域抽出装置、その作動方法およびプログラム |
US8799013B2 (en) | 2009-11-24 | 2014-08-05 | Penrad Technologies, Inc. | Mammography information system |
US8687860B2 (en) | 2009-11-24 | 2014-04-01 | Penrad Technologies, Inc. | Mammography statistical diagnostic profiler and prediction system |
US20120014578A1 (en) * | 2010-07-19 | 2012-01-19 | Qview Medical, Inc. | Computer Aided Detection Of Abnormalities In Volumetric Breast Ultrasound Scans And User Interface |
US20110238446A1 (en) * | 2010-03-27 | 2011-09-29 | Chaudhry Mundeep | Medical record entry systems and methods |
US9392986B2 (en) * | 2011-02-14 | 2016-07-19 | University Of Rochester | Method and apparatus for cone beam breast CT image-based computer-aided detection and diagnosis |
JP5788230B2 (ja) * | 2011-06-09 | 2015-09-30 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、超音波画像処理プログラム |
US9119559B2 (en) | 2011-06-16 | 2015-09-01 | Salient Imaging, Inc. | Method and system of generating a 3D visualization from 2D images |
KR20130012297A (ko) * | 2011-07-25 | 2013-02-04 | 삼성전자주식회사 | 병변 검출 장치, 병변 검출 방법 및 병변 진단 장치 |
KR20130049638A (ko) * | 2011-11-04 | 2013-05-14 | 삼성전자주식회사 | 병변 진단 장치 및 방법 |
KR20130062153A (ko) * | 2011-12-02 | 2013-06-12 | 삼성전자주식회사 | 병변 주변 정보를 이용한 병변 진단 장치 및 방법 |
CN102512247B (zh) * | 2011-12-16 | 2014-04-16 | 高磊 | 超声医师给乳腺小结节三维定位的装置 |
KR101331225B1 (ko) | 2012-05-17 | 2013-11-19 | 동국대학교 산학협력단 | 악성 흑색종 판별 장치 및 방법 |
TWI483711B (zh) * | 2012-07-10 | 2015-05-11 | Univ Nat Taiwan | Tumor detection system and method of breast ultrasound image |
US9378462B2 (en) | 2012-08-15 | 2016-06-28 | Questor Capital Holdings Ltd. | Probability mapping system |
WO2014034948A1 (ja) * | 2012-09-03 | 2014-03-06 | 株式会社東芝 | 超音波診断装置及び画像処理方法 |
KR101993716B1 (ko) * | 2012-09-28 | 2019-06-27 | 삼성전자주식회사 | 카테고리별 진단 모델을 이용한 병변 진단 장치 및 방법 |
JP6533363B2 (ja) * | 2012-12-14 | 2019-06-19 | テクトロニクス・インコーポレイテッドTektronix,Inc. | アーチファクト検出方法及びビデオアーチファクト検出器 |
US8976190B1 (en) * | 2013-03-15 | 2015-03-10 | Pme Ip Australia Pty Ltd | Method and system for rule based display of sets of images |
KR20150002284A (ko) * | 2013-06-28 | 2015-01-07 | 삼성전자주식회사 | 병변 검출 장치 및 방법 |
CN105683973B (zh) * | 2013-07-29 | 2019-06-25 | 皇家飞利浦有限公司 | 具有集成病变分级器的报告工具、方法和计算机存储介质 |
US10238368B2 (en) | 2013-09-21 | 2019-03-26 | General Electric Company | Method and system for lesion detection in ultrasound images |
JP6182045B2 (ja) * | 2013-10-11 | 2017-08-16 | キヤノン株式会社 | 画像処理装置およびその方法 |
KR20150098119A (ko) * | 2014-02-19 | 2015-08-27 | 삼성전자주식회사 | 의료 영상 내 거짓양성 병변후보 제거 시스템 및 방법 |
DE102014208411A1 (de) * | 2014-05-06 | 2015-11-12 | Siemens Aktiengesellschaft | Auswertung eines während einer Mammographie erzeugten Röntgenbildes einer Brust |
WO2015196300A1 (en) * | 2014-06-27 | 2015-12-30 | Sunnybrook Research Institute | Systems and methods for generating an imaging biomarker that indicates detectability or conspicuity of lesions in a mammographic image |
US9999402B2 (en) * | 2014-07-21 | 2018-06-19 | International Business Machines Corporation | Automatic image segmentation |
US11213220B2 (en) | 2014-08-11 | 2022-01-04 | Cubisme, Inc. | Method for determining in vivo tissue biomarker characteristics using multiparameter MRI matrix creation and big data analytics |
KR101599891B1 (ko) | 2014-09-19 | 2016-03-04 | 삼성전자주식회사 | 초음파 진단장치, 방법 및 컴퓨터 판독 가능한 저장매체 |
US9996935B2 (en) * | 2014-10-10 | 2018-06-12 | Edan Instruments, Inc. | Systems and methods of dynamic image segmentation |
US20160113546A1 (en) * | 2014-10-23 | 2016-04-28 | Khalifa University of Science, Technology & Research | Methods and systems for processing mri images to detect cancer |
EP3258890B1 (en) | 2015-02-17 | 2023-08-23 | Siemens Healthcare GmbH | Method and system for personalizing a vessel stent |
CN104739447A (zh) * | 2015-03-25 | 2015-07-01 | 井晓燕 | 集成超声成像诊断分析仪 |
US9922433B2 (en) | 2015-05-29 | 2018-03-20 | Moira F. Schieke | Method and system for identifying biomarkers using a probability map |
WO2016198990A1 (en) * | 2015-06-12 | 2016-12-15 | Koninklijke Philips N.V. | System and method for identifying cancerous tissue |
US20160377717A1 (en) * | 2015-06-29 | 2016-12-29 | Edan Instruments, Inc. | Systems and methods for adaptive sampling of doppler spectrum |
CN107613881B (zh) * | 2015-07-31 | 2020-05-15 | 安德拉生命科学公司 | 用于校正脂肪引起的像差的方法和*** |
US10548528B2 (en) * | 2015-08-07 | 2020-02-04 | Ryan James Appleby | Smartphone device for body analysis |
CN108603922A (zh) * | 2015-11-29 | 2018-09-28 | 阿特瑞斯公司 | 自动心脏体积分割 |
US10290109B2 (en) * | 2015-12-22 | 2019-05-14 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for cardiac image segmentation |
CN105869155B (zh) * | 2016-03-24 | 2018-11-16 | 天津大学 | 面部皮肤发红区域的提取方法 |
CN105931224A (zh) * | 2016-04-14 | 2016-09-07 | 浙江大学 | 基于随机森林算法的肝脏平扫ct图像病变识别方法 |
US11137462B2 (en) * | 2016-06-10 | 2021-10-05 | Board Of Trustees Of Michigan State University | System and method for quantifying cell numbers in magnetic resonance imaging (MRI) |
CN106127738B (zh) * | 2016-06-16 | 2018-12-11 | 上海荣盛生物药业有限公司 | 凝集试验判读方法 |
TWI574671B (zh) * | 2016-06-27 | 2017-03-21 | 太豪生醫股份有限公司 | ***影像的分析方法以及其電子裝置 |
US11182901B2 (en) * | 2016-06-29 | 2021-11-23 | Koninklijke Philips N.V. | Change detection in medical images |
US10776963B2 (en) | 2016-07-01 | 2020-09-15 | Cubismi, Inc. | System and method for forming a super-resolution biomarker map image |
CN109640830B (zh) * | 2016-07-14 | 2021-10-19 | 医视特有限公司 | 基于先例的超声聚焦 |
CN106326856A (zh) * | 2016-08-18 | 2017-01-11 | 厚凯(天津)医疗科技有限公司 | 一种手术图像处理方法及装置 |
US10276265B2 (en) | 2016-08-31 | 2019-04-30 | International Business Machines Corporation | Automated anatomically-based reporting of medical images via image annotation |
US10729396B2 (en) * | 2016-08-31 | 2020-08-04 | International Business Machines Corporation | Tracking anatomical findings within medical images |
CN106558045B (zh) * | 2016-10-20 | 2019-07-19 | 上海联影医疗科技有限公司 | 一种肺组织分割方法、装置,医学图像处理*** |
CN106530290A (zh) * | 2016-10-27 | 2017-03-22 | 朱育盼 | 医学图像分析方法和装置 |
US10192295B2 (en) * | 2016-11-09 | 2019-01-29 | AI Analysis, Inc. | Methods and systems for normalizing images |
FR3060169B1 (fr) * | 2016-12-09 | 2019-05-17 | Universite D'orleans | Detection de nerfs dans une serie d'images echographiques |
JP7054787B2 (ja) * | 2016-12-22 | 2022-04-15 | パナソニックIpマネジメント株式会社 | 制御方法、情報端末、及びプログラム |
JP2020510463A (ja) | 2017-01-27 | 2020-04-09 | アーテリーズ インコーポレイテッド | 全層畳み込みネットワークを利用する自動化されたセグメンテーション |
JP6583316B2 (ja) * | 2017-03-08 | 2019-10-02 | カシオ計算機株式会社 | 識別装置、識別方法及びプログラム |
GB201705876D0 (en) | 2017-04-11 | 2017-05-24 | Kheiron Medical Tech Ltd | Recist |
GB201705911D0 (en) * | 2017-04-12 | 2017-05-24 | Kheiron Medical Tech Ltd | Abstracts |
US11232853B2 (en) | 2017-04-21 | 2022-01-25 | Cubisme, Inc. | System and method for creating, querying, and displaying a MIBA master file |
USD855651S1 (en) | 2017-05-12 | 2019-08-06 | International Business Machines Corporation | Display screen with a graphical user interface for image-annotation classification |
US10380739B2 (en) * | 2017-08-15 | 2019-08-13 | International Business Machines Corporation | Breast cancer detection |
US11426075B1 (en) * | 2017-08-23 | 2022-08-30 | Lumicell, Inc. | System and method for residual cancer cell detection |
CN109684901B (zh) * | 2017-10-19 | 2023-06-06 | 富士通株式会社 | 图像处理装置和图像处理方法 |
TWI683276B (zh) | 2017-11-10 | 2020-01-21 | 太豪生醫股份有限公司 | 病灶偵測裝置及其方法 |
WO2019103912A2 (en) | 2017-11-22 | 2019-05-31 | Arterys Inc. | Content based image retrieval for lesion analysis |
CA3085441A1 (en) * | 2018-01-08 | 2019-07-11 | Progenics Pharmaceuticals, Inc. | Systems and methods for rapid neural network-based image segmentation and radiopharmaceutical uptake determination |
EP3511866A1 (en) * | 2018-01-16 | 2019-07-17 | Koninklijke Philips N.V. | Tissue classification using image intensities and anatomical positions |
US11058390B1 (en) * | 2018-02-23 | 2021-07-13 | Robert Edwin Douglas | Image processing via a modified segmented structure |
KR102122767B1 (ko) * | 2018-03-05 | 2020-06-15 | 고려대학교산학협력단 | 초음파 검사 지원 장치 및 방법 |
KR102072476B1 (ko) * | 2018-05-09 | 2020-02-03 | 재단법인 대구경북첨단의료산업진흥재단 | 영상 유도 수술을 위한 종양 검출 장치 및 그 방법 |
CN117438053A (zh) * | 2018-08-27 | 2024-01-23 | 卡西欧计算机株式会社 | 显示控制装置、显示控制***、显示控制方法及记录介质 |
WO2020142131A2 (en) | 2018-10-19 | 2020-07-09 | Inkbit, LLC | High-speed metrology |
CN109447975A (zh) * | 2018-11-01 | 2019-03-08 | 肖湘江 | 基于图像的手掌红色斑点检测方法 |
US11354466B1 (en) | 2018-11-02 | 2022-06-07 | Inkbit, LLC | Machine learning for additive manufacturing |
AU2019374148A1 (en) | 2018-11-02 | 2021-05-27 | Inkbit, LLC | Intelligent additive manufacturing |
EP3856481A2 (en) | 2018-11-16 | 2021-08-04 | Inkbit, LLC | Inkjet 3d printing of multi-component resins |
CN109785296B (zh) * | 2018-12-25 | 2023-07-04 | 西安电子科技大学 | 一种基于cta图像的三维球形指数测定方法 |
US11793594B2 (en) | 2018-12-31 | 2023-10-24 | Lumicell, Inc. | System and method for thresholding for residual cancer cell detection |
WO2020146490A1 (en) * | 2019-01-08 | 2020-07-16 | Inkbit, LLC | Depth reconstruction in additive fabrication |
WO2020146481A1 (en) | 2019-01-08 | 2020-07-16 | Inkbit, LLC | Reconstruction of surfaces for additive manufacturing |
CN111476794B (zh) * | 2019-01-24 | 2023-10-20 | 武汉兰丁智能医学股份有限公司 | 一种基于unet的宫颈病理组织分割方法 |
CN110009600A (zh) * | 2019-02-14 | 2019-07-12 | 腾讯科技(深圳)有限公司 | 一种医学图像区域过滤方法、装置及存储介质 |
KR20200107555A (ko) * | 2019-03-08 | 2020-09-16 | 에스케이텔레콤 주식회사 | 영상 분석 장치, 방법 및 이에 이용되는 영상 분석 모델의 생성 방법 |
CN110334722B (zh) * | 2019-03-29 | 2022-07-05 | 上海联影智能医疗科技有限公司 | 图像分类方法、装置、计算机设备和存储介质 |
TWI728369B (zh) * | 2019-05-24 | 2021-05-21 | 臺北醫學大學 | 人工智慧雲端膚質與皮膚病灶辨識方法及其系統 |
US11334994B2 (en) * | 2019-05-24 | 2022-05-17 | Lunit Inc. | Method for discriminating suspicious lesion in medical image, method for interpreting medical image, and computing device implementing the methods |
KR102245219B1 (ko) * | 2019-05-24 | 2021-04-27 | 주식회사 루닛 | 의료 영상에서 악성의심 병변을 구별하는 방법, 이를 이용한 의료 영상 판독 방법 및 컴퓨팅 장치 |
WO2021006991A2 (en) * | 2019-06-12 | 2021-01-14 | Carnegie Mellon University | System and method for vessel segmentation |
KR102338018B1 (ko) * | 2019-07-30 | 2021-12-10 | 주식회사 힐세리온 | 초음파 이미지상의 특징점들을 이용한 초음파 지방간 자동 진단 장치 및 이를 이용한 원격 의료 진단 방법 |
US11712837B2 (en) | 2019-11-01 | 2023-08-01 | Inkbit, LLC | Optical scanning for industrial metrology |
US10994477B1 (en) | 2019-11-01 | 2021-05-04 | Inkbit, LLC | Optical scanning for industrial metrology |
CN110833395B (zh) * | 2019-11-15 | 2022-06-28 | 广州七喜医疗设备有限公司 | 乳腺脂肪水平确定方法和装置 |
JP7294996B2 (ja) * | 2019-11-28 | 2023-06-20 | 富士フイルムヘルスケア株式会社 | 超音波診断装置及び表示方法 |
CN111275719B (zh) * | 2020-01-19 | 2023-04-07 | 推想医疗科技股份有限公司 | 钙化假阳识别及其模型训练的方法、装置、终端及介质 |
KR102216697B1 (ko) * | 2020-02-28 | 2021-02-17 | 주식회사 루닛 | 의료 영상 장치 및 의료 영상 처리 방법 |
CN111444971A (zh) * | 2020-03-31 | 2020-07-24 | 联想(北京)有限公司 | 信息处理方法、装置及电子设备和存储介质 |
CN111603199B (zh) * | 2020-04-24 | 2023-03-14 | 中国人民解放军总医院第二医学中心 | 一种基于体表定位测量仪的三维重建超声诊断*** |
CN111832574B (zh) * | 2020-07-13 | 2023-06-30 | 福建省妇幼保健院 | 用于检测人***瘤状病毒感染性病变的图像识别方法 |
US10994490B1 (en) | 2020-07-31 | 2021-05-04 | Inkbit, LLC | Calibration for additive manufacturing by compensating for geometric misalignments and distortions between components of a 3D printer |
US20220067919A1 (en) * | 2020-08-26 | 2022-03-03 | GE Precision Healthcare LLC | System and method for identifying a tumor or lesion in a probabilty map |
CN112116027A (zh) * | 2020-09-29 | 2020-12-22 | 宁波工程学院 | 一种基于oct成像图像的光学强度及梯度的皮肤癌分类方法 |
CN112184683A (zh) * | 2020-10-09 | 2021-01-05 | 深圳度影医疗科技有限公司 | 一种超声图像的识别方法、终端设备及存储介质 |
CN112508850B (zh) * | 2020-11-10 | 2021-07-20 | 广州柏视医疗科技有限公司 | 基于深度学习的甲状腺细胞病理切片恶性区域的检测方法 |
CN112785609B (zh) * | 2021-02-07 | 2022-06-03 | 重庆邮电大学 | 一种基于深度学习的cbct牙齿分割方法 |
CN113034460B (zh) * | 2021-03-21 | 2022-09-09 | 湖南科迈森医疗科技有限公司 | 子宫内膜腺体密度预估方法 |
CN113408596B (zh) * | 2021-06-09 | 2022-09-30 | 北京小白世纪网络科技有限公司 | 病理图像的处理方法、装置、电子设备及可读存储介质 |
CN113408595B (zh) * | 2021-06-09 | 2022-12-13 | 北京小白世纪网络科技有限公司 | 病理图像处理的方法、装置、电子设备及可读存储介质 |
US20230027734A1 (en) * | 2021-07-16 | 2023-01-26 | Johnson & Johnson Enterprise Innovation Inc. | System and Method for Predicting the Risk of Future Lung Cancer |
CN113657553B (zh) * | 2021-09-01 | 2023-12-26 | 什维新智医疗科技(上海)有限公司 | 一种结节回声类型判断装置 |
CN113952580A (zh) * | 2021-11-18 | 2022-01-21 | 苏州大学 | 一种正念冥想训练的方法、装置、设备以及存储介质 |
CN115272313B (zh) * | 2022-09-27 | 2022-12-30 | 广州辉博信息技术有限公司 | 基于深度图像的肌肉均衡程度分析方法、***及设备 |
CN116958147B (zh) * | 2023-09-21 | 2023-12-22 | 青岛美迪康数字工程有限公司 | 基于深度图像特征的目标区域确定方法、装置和设备 |
CN117640932B (zh) * | 2024-01-25 | 2024-04-26 | 陕西通达伟业医疗供应链管理有限公司 | 用于远程医疗的神经内科图像压缩传输方法 |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5003979A (en) * | 1989-02-21 | 1991-04-02 | University Of Virginia | System and method for the noninvasive identification and display of breast lesions and the like |
IL119767A (en) * | 1993-08-13 | 1998-02-08 | Sophis View Tech Ltd | System and method for diagnosis of living tissue diseases |
US5657362A (en) * | 1995-02-24 | 1997-08-12 | Arch Development Corporation | Automated method and system for computerized detection of masses and parenchymal distortions in medical images |
DE69630935T2 (de) * | 1995-09-29 | 2004-11-04 | Koninklijke Philips Electronics N.V. | Bilverarbeitungsverfahren und -vorrichtung zur automatischen Erfassung von Bereichen eines vorbestimmten Krebstyps in einem Intensitätsbild |
US6317617B1 (en) * | 1997-07-25 | 2001-11-13 | Arch Development Corporation | Method, computer program product, and system for the automated analysis of lesions in magnetic resonance, mammogram and ultrasound images |
US5984870A (en) * | 1997-07-25 | 1999-11-16 | Arch Development Corporation | Method and system for the automated analysis of lesions in ultrasound images |
US6173034B1 (en) * | 1999-01-25 | 2001-01-09 | Advanced Optical Technologies, Inc. | Method for improved breast x-ray imaging |
US6901156B2 (en) * | 2000-02-04 | 2005-05-31 | Arch Development Corporation | Method, system and computer readable medium for an intelligent search workstation for computer assisted interpretation of medical images |
IT1320956B1 (it) * | 2000-03-24 | 2003-12-18 | Univ Bologna | Metodo, e relativa apparecchiatura, per la rilevazione automatica dimicrocalcificazioni in segnali digitali di tessuto mammario. |
US6956975B2 (en) * | 2001-04-02 | 2005-10-18 | Eastman Kodak Company | Method for improving breast cancer diagnosis using mountain-view and contrast-enhancement presentation of mammography |
FR2825816B1 (fr) * | 2001-06-08 | 2004-07-23 | Ge Med Sys Global Tech Co Llc | Procede d'affichage d'images d'organes |
US7054473B1 (en) * | 2001-11-21 | 2006-05-30 | R2 Technology, Inc. | Method and apparatus for an improved computer aided diagnosis system |
AU2003219689A1 (en) * | 2002-02-08 | 2003-09-02 | Maryellen L. Giger | Method and system for risk-modulated diagnosis of disease |
GB2395263A (en) * | 2002-11-12 | 2004-05-19 | Qinetiq Ltd | Image analysis |
US6999549B2 (en) * | 2002-11-27 | 2006-02-14 | Ge Medical Systems Global Technology, Llc | Method and apparatus for quantifying tissue fat content |
US7298883B2 (en) * | 2002-11-29 | 2007-11-20 | University Of Chicago | Automated method and system for advanced non-parametric classification of medical images and lesions |
US7769215B2 (en) * | 2003-11-28 | 2010-08-03 | The University Of Chicago | Method, system, and medium for prevalence-based computerized analysis of medical images and information |
US7534210B2 (en) * | 2004-02-03 | 2009-05-19 | Siemens Medical Solutions Usa, Inc. | Methods for adaptively varying gain during ultrasound agent quantification |
US7736313B2 (en) * | 2004-11-22 | 2010-06-15 | Carestream Health, Inc. | Detecting and classifying lesions in ultrasound images |
US7840046B2 (en) * | 2006-06-27 | 2010-11-23 | Siemens Medical Solutions Usa, Inc. | System and method for detection of breast masses and calcifications using the tomosynthesis projection and reconstructed images |
US8238637B2 (en) * | 2006-10-25 | 2012-08-07 | Siemens Computer Aided Diagnosis Ltd. | Computer-aided diagnosis of malignancies of suspect regions and false positives in images |
US20090226057A1 (en) * | 2008-03-04 | 2009-09-10 | Adi Mashiach | Segmentation device and method |
US20100124364A1 (en) * | 2008-11-19 | 2010-05-20 | Zhimin Huo | Assessment of breast density and related cancer risk |
US20100142790A1 (en) * | 2008-12-04 | 2010-06-10 | New Medical Co., Ltd. | Image processing method capable of enhancing contrast and reducing noise of digital image and image processing device using same |
-
2009
- 2009-12-21 US US12/643,337 patent/US20100158332A1/en not_active Abandoned
- 2009-12-22 CA CA2783867A patent/CA2783867A1/en not_active Abandoned
- 2009-12-22 CN CN200980152219.0A patent/CN102438529B/zh not_active Expired - Fee Related
- 2009-12-22 EP EP09833984.9A patent/EP2378978B1/en not_active Not-in-force
- 2009-12-22 JP JP2011541052A patent/JP2012512672A/ja active Pending
- 2009-12-22 WO PCT/CA2009/001881 patent/WO2010071999A1/en active Application Filing
- 2009-12-22 EP EP13193330.1A patent/EP2712554A1/en not_active Withdrawn
- 2009-12-22 CN CN201410055212.XA patent/CN103854028A/zh active Pending
- 2009-12-22 CN CN201410055240.1A patent/CN103824290A/zh active Pending
-
2013
- 2013-03-20 US US13/847,789 patent/US20130343626A1/en not_active Abandoned
-
2014
- 2014-12-10 HK HK14112431.6A patent/HK1199127A1/xx unknown
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140057130A (ko) * | 2012-10-25 | 2014-05-12 | 삼성전자주식회사 | 이미지 처리 시스템 |
KR102128325B1 (ko) | 2012-10-25 | 2020-07-01 | 삼성전자주식회사 | 이미지 처리 시스템 |
US10383602B2 (en) | 2014-03-18 | 2019-08-20 | Samsung Electronics Co., Ltd. | Apparatus and method for visualizing anatomical elements in a medical image |
WO2016051764A1 (ja) * | 2014-09-29 | 2016-04-07 | 富士フイルム株式会社 | 光音響画像生成装置 |
JP2016187748A (ja) * | 2016-08-12 | 2016-11-04 | コニカミノルタ株式会社 | 超音波画像診断装置及び超音波画像診断方法 |
KR20220040507A (ko) * | 2017-05-11 | 2022-03-30 | 베라톤 인코포레이티드 | 확률 맵 기반의 초음파 스캐닝 |
KR102409090B1 (ko) | 2017-05-11 | 2022-06-15 | 베라톤 인코포레이티드 | 확률 맵 기반의 초음파 스캐닝 |
US11348228B2 (en) | 2017-06-26 | 2022-05-31 | The Research Foundation For The State University Of New York | System, method, and computer-accessible medium for virtual pancreatography |
JP2021527473A (ja) * | 2018-06-14 | 2021-10-14 | ケイロン メディカル テクノロジーズ リミテッド | 即時精密検査 |
KR102186632B1 (ko) | 2019-01-07 | 2020-12-02 | 재단법인대구경북과학기술원 | 의료 영상의 분석 모델을 학습시키는 학습 장치 및 그 학습 방법 |
KR20200092466A (ko) * | 2019-01-07 | 2020-08-04 | 재단법인대구경북과학기술원 | 의료 영상의 분석 모델을 학습시키는 학습 장치 및 그 학습 방법 |
JP2021118840A (ja) * | 2020-01-30 | 2021-08-12 | 株式会社日立製作所 | 医用画像処理装置及び医用画像処理方法 |
JP7382240B2 (ja) | 2020-01-30 | 2023-11-16 | 富士フイルムヘルスケア株式会社 | 医用画像処理装置及び医用画像処理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102438529B (zh) | 2014-10-22 |
CN103824290A (zh) | 2014-05-28 |
HK1199127A1 (en) | 2015-06-19 |
US20100158332A1 (en) | 2010-06-24 |
EP2712554A1 (en) | 2014-04-02 |
EP2378978B1 (en) | 2014-11-12 |
US20130343626A1 (en) | 2013-12-26 |
CA2783867A1 (en) | 2010-07-01 |
CN103854028A (zh) | 2014-06-11 |
EP2378978A1 (en) | 2011-10-26 |
CN102438529A (zh) | 2012-05-02 |
WO2010071999A1 (en) | 2010-07-01 |
EP2378978A4 (en) | 2013-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2378978B1 (en) | Method and system for automated generation of surface models in medical images | |
WO2021129323A1 (zh) | 超声图像病灶描述方法、装置、计算机设备和存储介质 | |
US7903861B2 (en) | Method for classifying breast tissue density using computed image features | |
US5657362A (en) | Automated method and system for computerized detection of masses and parenchymal distortions in medical images | |
US7466848B2 (en) | Method and apparatus for automatically detecting breast lesions and tumors in images | |
US6795521B2 (en) | Computer-aided diagnosis system for thoracic computer tomography images | |
US6309353B1 (en) | Methods and apparatus for tumor diagnosis | |
US8340388B2 (en) | Systems, computer-readable media, methods, and medical imaging apparatus for the automated detection of suspicious regions of interest in noise normalized X-ray medical imagery | |
JP2010504129A (ja) | 肺結節の高度コンピュータ支援診断 | |
US10376230B2 (en) | Obtaining breast density measurements and classifications | |
CN109919254B (zh) | 乳腺密度分类方法、***、可读存储介质和计算机设备 | |
Ganvir et al. | Filtering method for pre-processing mammogram images for breast cancer detection | |
CN116580068A (zh) | 一种基于点云配准的多模态医学配准方法 | |
US20110222752A1 (en) | Microcalcification enhancement from digital mammograms | |
CN107767386B (zh) | 超声图像处理方法及装置 | |
Edwin et al. | Liver and tumour segmentation from abdominal CT images using adaptive threshold method | |
EP3977400A1 (en) | A computer based method for classifying a mass of an organ as a cyst | |
CN111383323B (zh) | 一种超声成像方法和***以及超声图像处理方法和*** | |
Joda et al. | Digital mammogram enhancement based on automatic histogram clipping | |
CN112541907A (zh) | 一种图像识别方法、装置、服务器及介质 | |
Ivković et al. | Segmentation and classification of Bi-Rads medical images with the imaging biomarkers according to level of detail | |
Sample | Computer assisted screening of digital mammogram images | |
Liu et al. | Segmentation of Mammography Images Based on Spectrum Clustering Method | |
CN114651273A (zh) | 用于图像标准化的方法和*** | |
BHATTACHARJEE et al. | Identification of Abnormal Masses in Digital Mammogram Using Statistical Decision Making |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20120215 |