JP6987342B2 - Image processing equipment, methods and programs - Google Patents

Image processing equipment, methods and programs Download PDF

Info

Publication number
JP6987342B2
JP6987342B2 JP2018030552A JP2018030552A JP6987342B2 JP 6987342 B2 JP6987342 B2 JP 6987342B2 JP 2018030552 A JP2018030552 A JP 2018030552A JP 2018030552 A JP2018030552 A JP 2018030552A JP 6987342 B2 JP6987342 B2 JP 6987342B2
Authority
JP
Japan
Prior art keywords
lesion
suspected
image processing
area
processing device
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
Application number
JP2018030552A
Other languages
Japanese (ja)
Other versions
JP2019141478A (en
Inventor
英恵 吉田
昌広 影山
昌宏 荻野
Original Assignee
富士フイルムヘルスケア株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 富士フイルムヘルスケア株式会社 filed Critical 富士フイルムヘルスケア株式会社
Priority to JP2018030552A priority Critical patent/JP6987342B2/en
Publication of JP2019141478A publication Critical patent/JP2019141478A/en
Application granted granted Critical
Publication of JP6987342B2 publication Critical patent/JP6987342B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Description

本発明は、医用画像を処理する画像処理装置、その方法及びプログラムに関する。 The present invention relates to an image processing apparatus for processing a medical image, a method and a program thereof.

X線CT(X−ray Computed Tomography)装置やMRI(Magnetic Resonance Imaging)装置等に代表される医用画像検査装置を用いた診断では、撮影された三次元医用画像を、連続した二次元断面として再構成し、その二次元断面画像を観察して読影を行うことが一般的である。 In diagnosis using a medical image inspection device represented by an X-ray CT (X-ray Computed Tomography) device or an MRI (Magnetic Response Imaging) device, the captured three-dimensional medical image is reproduced as a continuous two-dimensional cross section. It is common to construct and observe the two-dimensional cross-sectional image for interpretation.

撮影装置の高度化により、生成される三次元医用画像の三次元分解能も向上しており、データサイズは増加する傾向にある。特に、先に述べた二次元断面の生成間隔はより細かくすることが可能となり、医用画像上に現れる病変のより詳細な観察が可能となってきているが、結果的に三次元医用画像あたりの枚数も増加している。 With the sophistication of imaging equipment, the three-dimensional resolution of the generated three-dimensional medical image is also improved, and the data size tends to increase. In particular, it has become possible to make the generation interval of the two-dimensional cross section described above finer, and it has become possible to observe the lesions appearing on the medical image in more detail. The number of sheets is also increasing.

また特にCT装置においては、低線量で高画質な三次元医用画像の撮影が可能になってきたこともあり、CT画像の撮影機会も増加傾向にある。 Further, in particular, in a CT apparatus, it has become possible to take a three-dimensional medical image with a low dose and high image quality, and the opportunity to take a CT image is increasing.

これらの理由により、膨大な医用三次元医用画像を読影する際に医師や技師にかかる負担を軽減し、主に病変の見落としを防ぐためにCAD(Computer Aided Detection)という技術の開発が進められている。CADは、コンピュータにより陰影の検出やサイズ計測、陰影の正常/異常の識別や異常陰影の病変種類の区別等を、画像処理技術を応用して自動あるいは半自動で行うことを目指したものである。 For these reasons, a technology called CAD (Computer Aided Detection) is being developed to reduce the burden on doctors and technicians when interpreting a huge amount of medical three-dimensional medical images and mainly to prevent oversight of lesions. .. CAD aims to automatically or semi-automatically detect and measure the size of shadows, identify normal / abnormal shadows, and distinguish lesion types of abnormal shadows by applying image processing technology.

ここでは、医用画像の特徴から病変の疑いが高い陰影を提示することを目的とするCADについて述べる。画像処理分野の特に画像認識技術の発展に伴い,CADは精度を向上させてきた。一方で,医用画像における病変画像の画像特徴はおなじ病名であったとしても多岐にわたり,主に機械学習を用いて異常陰影を識別する場合に,学習データが不十分であるために求められる精度を達成できないという課題があった。 Here, CAD for the purpose of presenting a shadow with a high suspicion of a lesion due to the characteristics of a medical image will be described. With the development of image processing technology, especially in the field of image processing, CAD has improved the accuracy. On the other hand, the image features of lesion images in medical images are diverse even if they have the same disease name, and when identifying abnormal shadows mainly using machine learning, the accuracy required due to insufficient learning data is obtained. There was a problem that it could not be achieved.

これらの課題を解決するため、特許文献1には、顔画像の識別精度を向上させる方法として、入力データを、互いに異なる画像特徴を基にして2クラス以上に分類する複数の識別器を用いてクラスに分類し、複数の属性情報に基づいて、識別器のパラメータを調整する方法などが提案されている。 In order to solve these problems, Patent Document 1 uses a plurality of classifiers that classify input data into two or more classes based on different image features as a method for improving the identification accuracy of facial images. A method of classifying into classes and adjusting the parameters of the classifier based on multiple attribute information has been proposed.

特開2017-102671公報Japanese Unexamined Patent Publication No. 2017-102671

CADによる病変疑い陰影の自動検出では一般的には、画像から得られる特徴量を用いて病変の疑わしさを定義する定義式を設定し、その病変の疑わしさが高い陰影を提示する。ここで定義式を設定する際に,病変であると既に診断を受けている過去の画像の情報を参考にする場合,おなじ病名であっても画像の特徴が多岐にわたることや、診断済の過去の画像の量が十分でないことにより,定義式の精度が低くなるという課題がある。また、複数の画像データに対してラベリングをする技術において、同じラベルの学習データであっても見た目の変化が大きいという課題に対する対策としてこれまで提案されている技術においては,互いに異なる画像特徴を基にした複数の属性情報に基づいて、識別器のパラメータを調整することは可能ではあるが,分類の基にする属性情報については、全ての場合について最適な設定となっているとは限らないという課題がある。 In the automatic detection of a suspicious lesion shadow by CAD, a definition formula that defines the suspicion of a lesion is generally set using a feature amount obtained from an image, and a shadow with a high suspicion of the lesion is presented. When setting the definition formula here, when referring to the information of the past image that has already been diagnosed as a lesion, the characteristics of the image are diverse even if the disease name is the same, and the past that has been diagnosed. There is a problem that the accuracy of the definition formula becomes low because the amount of images in the image is not sufficient. Further, in the technique of labeling a plurality of image data, the techniques proposed so far as a countermeasure against the problem that the appearance changes greatly even if the learning data has the same label are based on different image features. Although it is possible to adjust the parameters of the classifier based on the multiple attribute information set in, the attribute information used as the basis for classification may not be the optimum setting in all cases. There are challenges.

本発明の目的は、参考にする過去の画像の量が不十分である場合にも高精度な病変の疑わしさの定義式の設定が可能な医用画像処理装置、医用画像表示方法、及びプログラムを提供することにある。 An object of the present invention is to provide a medical image processing device, a medical image display method, and a program capable of setting a highly accurate definition formula of suspicion of a lesion even when the amount of past images to be referred to is insufficient. To provide.

特に、従来の技術による定義式の設定では病変の検出率が低い場合、例えば放射線科医が読影した場合の結節検出率6割程度を下回る場合に有効な技術を提供する。 In particular, we provide an effective technique when the detection rate of lesions is low in the setting of the definition formula by the conventional technique, for example, when the nodule detection rate is less than about 60% when the radiologist interprets the image.

上記の課題を解決するために、本発明においては、医用画像を処理する画像処理装置であって、医用画像内の複数の病変疑い領域と、病変疑い領域の複数の特徴量を表示する病変疑い情報表示部と、表示された複数の特徴量から選択された特徴量により、病変疑い領域を分類する病変疑い領域分類部と、を備える構成の画像処理装置を提供する。 In order to solve the above problems, the present invention is an image processing device that processes a medical image, and displays a plurality of suspected lesion regions and a plurality of feature quantities of the suspected lesion region in the medical image. Provided is an image processing apparatus having a configuration including an information display unit and a lesion suspected area classification unit that classifies a lesion suspected area according to a feature amount selected from a plurality of displayed feature amounts.

また、上記の課題を解決するため、本発明においては、画像処理装置により医用画像を処理する画像処理方法であって、画像処理装置は、医用画像内の複数の病変疑い領域と、病変疑い領域の複数の特徴量をモニタに表示し、表示された複数の特徴量から選択された特徴量により、病変疑い領域を分類する画像処理方法を提供する。 Further, in order to solve the above-mentioned problems, the present invention is an image processing method for processing a medical image by an image processing apparatus, wherein the image processing apparatus has a plurality of suspected lesion regions and suspected lesion regions in the medical image. Provided is an image processing method for displaying a plurality of feature quantities of the above on a monitor and classifying the suspected lesion region according to the feature quantities selected from the plurality of displayed feature quantities.

更に、上記の課題を解決するため、本発明においては、画像処理装置に医用画像を処理させる画像処理プログラムであって、画像処理装置が、医用画像内の複数の病変疑い領域と、病変疑い領域の複数の特徴量をモニタに表示し、表示された複数の特徴量から選択された特徴量により、病変疑い領域を分類する、よう動作させる画像処理プログラムを提供する。 Further, in order to solve the above-mentioned problems, in the present invention, an image processing program for causing an image processing device to process a medical image, wherein the image processing device has a plurality of suspected lesion regions and suspected lesion regions in the medical image. Provided is an image processing program that displays a plurality of feature quantities of the above on a monitor and classifies the suspected lesion region according to the feature quantities selected from the displayed plurality of feature quantities.

本発明によれば、参考にする過去の医用画像の量が不十分である場合にも、高精度な病変の疑わしさの定義式の設定が可能になる。 According to the present invention, it is possible to set a highly accurate definition formula for suspicion of a lesion even when the amount of past medical images to be referred to is insufficient.

実施例1に係る医用画像処理装置を含むシステム構成例を示す図である。It is a figure which shows the system configuration example including the medical image processing apparatus which concerns on Example 1. FIG. 実施例1に係る、病変疑い領域分類部と病変尤度設定部の構成例を示す図である。It is a figure which shows the structural example of the lesion suspected area classification part and the lesion likelihood setting part which concerns on Example 1. FIG. 実施例1に係る、病変疑い領域分類式定義処理フローの一例を示す図である。It is a figure which shows an example of the lesion suspected area classification type definition processing flow which concerns on Example 1. FIG. 実施例1に係る、各病変疑い領域の特徴量と、対応する病変疑い領域の画像を重畳表示した図である。It is a figure which superposed and displayed the feature amount of each lesion suspected area and the image of the corresponding lesion suspected area which concerns on Example 1. FIG. 実施例1に係る、各病変疑い領域の他の特徴量を、人体モデルの画像上に重畳表示した図である。It is a figure which superposed and displayed the other feature amount of each lesion suspected area which concerns on Example 1. on the image of a human body model. 実施例1に係る、病変疑い情報表示部による表示から、分類の必要性を判断するための画面表示例を示す図である。It is a figure which shows the screen display example for determining the necessity of classification from the display by the lesion suspicion information display part which concerns on Example 1. FIG. 実施例1に係る、病変尤度算出式の一例を示す図である。It is a figure which shows an example of the lesion likelihood calculation formula which concerns on Example 1. FIG. 実施例1に係る、病変尤度算出式定義処理フローの一例を示す図である。It is a figure which shows an example of the lesion likelihood calculation formula definition processing flow which concerns on Example 1. FIG. 実施例1に係る、病変尤度算出式の一例を示す図である。It is a figure which shows an example of the lesion likelihood calculation formula which concerns on Example 1. FIG. 実施例1に係る、病変尤度算出処理フローの一例を示す図である。It is a figure which shows an example of the lesion likelihood calculation processing flow which concerns on Example 1. FIG. 現状の医用画像処理技術の課題を説明するための図である。It is a figure for demonstrating the problem of the present medical image processing technology. 実施例1に係る、医用画像処理装置を利用した場合の評価データの評価を説明するための図である。It is a figure for demonstrating the evaluation of the evaluation data at the time of using the medical image processing apparatus which concerns on Example 1. FIG.

以下、本発明の画像処理装置の実施の形態を図面に従い順次説明する。本発明の画像処理装置は、医用画像を処理して、医用画像上の病変疑い領域を分類するが、病変疑い領域とは、読影医の医学的知識や疾病診断に対する医学的根拠(エビデンス)等に基づいて判断される、病変の疑いが高い点、および領域を意味する。ここで診断対象となる病変とは、医用画像上に現れた場合に、周囲即ち病変の疑いが低い領域との輝度の違いや、その分布の違いから判断できる可能性が高いものとする。 Hereinafter, embodiments of the image processing apparatus of the present invention will be sequentially described with reference to the drawings. The image processing apparatus of the present invention processes a medical image to classify a suspected lesion area on the medical image, and the suspected lesion area is the medical knowledge of an image interpreter and the medical basis (evidence) for disease diagnosis. It means a point and area with a high suspicion of lesion, which is judged based on. Here, it is highly possible that the lesion to be diagnosed can be determined from the difference in brightness and the difference in its distribution from the surrounding area, that is, the area where the suspicion of the lesion is low, when it appears on the medical image.

実施例1は、医用画像を処理する画像処理装置であって、医用画像内の複数の病変疑い領域と、病変疑い領域の複数の特徴量を表示する病変疑い情報表示部と、表示された複数の特徴量から選択された特徴量により、病変疑い領域を分類する病変疑い領域分類部と、を備える構成の画像処理装置の実施例であり、参考にする過去の医用画像の量が不十分である場合にも、高精度な病変の疑わしさの定義式の設定が可能とし、高精度な病変疑い箇所を検出することができる医用画像処理装置、及びその方法、プログラムの実施例である。 The first embodiment is an image processing apparatus for processing a medical image, which includes a plurality of suspected lesion regions in the medical image, a suspected lesion information display unit for displaying a plurality of feature quantities of the suspected lesion region, and a plurality of displayed parts. This is an example of an image processing device having a configuration including a suspected lesion area classification unit for classifying suspected lesion areas according to the feature amount selected from the feature amounts of the above, and the amount of past medical images to be referred to is insufficient. This is an example of a medical image processing apparatus, a method, and a program capable of setting a highly accurate definition formula of suspicion of a lesion and detecting a suspicious part of a lesion with high accuracy.

本実施例においては、CT医用画像撮影装置により得られる再構成三次元医用画像について述べるが、本技術は他の医用画像撮影装置により得られるデータについても応用可能である。例えばMRI撮影装置等により得られるデータであっても、複数の二次元断面の積み重ねとして表現できる三次元画像を得るもので、画素分布に病変特徴が現れるとされているものであれば適用することができる。 In this embodiment, a reconstructed three-dimensional medical image obtained by a CT medical imaging device will be described, but this technique can also be applied to data obtained by another medical imaging device. For example, even if the data is obtained by an MRI imaging device or the like, a three-dimensional image that can be expressed as a stack of a plurality of two-dimensional cross sections is obtained, and if it is said that lesion features appear in the pixel distribution, it should be applied. Can be done.

上述したように、本実施例における病変疑い領域とは、読影医の医学的知識や当該疾病診断に対する医学的根拠(エビデンス)等に基づいて判断される、病変の疑いが高い点、および領域、部位を指す。ここで対象となる病変とは、医用画像上に現れた場合に、周囲即ち病変の疑いが低い領域との輝度の違いや、分布の違いから判断できる可能性が高いものとする。例えば肺結節の場合は一般に、そのCT値が周辺の空気領域よりも高い画素を多く含む領域としてCT画像上に現れることが知られている。胸部CT画像上で高輝度画素を多く含む他のオブジェクトとしては、血管や骨があるが、高輝度値の分布形状に応じて、それらと区別し、病変すなわちこの例における肺結節である疑いの高さを判別できると言われている。 As described above, the lesion suspicious area in this example is a point and area where there is a high suspicion of a lesion, which is determined based on the medical knowledge of the image interpreter and the medical basis (evidence) for the diagnosis of the disease. Refers to the part. Here, it is highly possible that the target lesion can be determined from the difference in brightness and the difference in distribution from the surrounding area, that is, the area where the suspicion of the lesion is low, when it appears on the medical image. For example, in the case of a lung nodule, it is generally known that its CT value appears on a CT image as a region containing many pixels higher than the surrounding air region. Other objects that contain a lot of high-intensity pixels on the chest CT image are blood vessels and bones, which are distinguished from them according to the distribution shape of the high-intensity values, and are suspected to be lesions, that is, lung nodules in this example. It is said that the height can be determined.

図1に、本実施例に係る医用画像処理装置を含むシステム構成の一例を示す。同図に示すように本実施例のシステムは、医用画像処理装置11と、操作者の入力等を受診し医用画像処理装置11に送信する入力装置10と、医用画像処表示装置11から得られる医用画像と病変疑い情報を表示するモニタ12とから構成されている。本システム構成は、通常のコンピュータで構成可能であり、医用画像処理装置は種々の処理プログラムを実行する中央処理部(CPU)とそれらの処理プログラムやデータを記憶するメモリで実現できる。 FIG. 1 shows an example of a system configuration including the medical image processing apparatus according to the present embodiment. As shown in the figure, the system of this embodiment is obtained from a medical image processing device 11, an input device 10 that receives an operator's input and the like and transmits the input to the medical image processing device 11, and a medical image processing display device 11. It is composed of a medical image and a monitor 12 for displaying suspected lesion information. This system configuration can be configured by an ordinary computer, and the medical image processing apparatus can be realized by a central processing unit (CPU) that executes various processing programs and a memory that stores those processing programs and data.

医用画像処理装置11は、医用画像を記憶する医用画像記憶部20と、各医用画像に対応する病変疑い情報を記憶する病変疑い情報記憶部21と、複数の医用画像と病変疑い情報を一度にモニタ12に表示する病変疑い情報表示部22と、病変疑い領域を分類する病変疑い領域分類部23と、分類した病変疑い領域について病変尤度を設定する病変尤度設定部24とから構成されている。 The medical image processing device 11 includes a medical image storage unit 20 that stores medical images, a lesion suspicious information storage unit 21 that stores lesion suspicious information corresponding to each medical image, and a plurality of medical images and lesion suspicious information at once. It is composed of a lesion suspicion information display unit 22 displayed on the monitor 12, a lesion suspicion area classification unit 23 for classifying the lesion suspicious area, and a lesion likelihood setting unit 24 for setting the lesion likelihood for the classified lesion suspicion area. There is.

図2は、病変疑い領域分類部23と病変尤度設定部24のシステム構成の一例を示す。同図の上部に示すように病変疑い領域分類部23は、入力装置10からの入力を受けて病変疑い領域分類式を定義する病変疑い領域分類式定義部30と、その式を記憶する病変疑い領域分類式記憶部31、病変疑い領域分類式に従って入力病変疑い領域を分類する病変疑い領域分類実行部32とから構成されている。 FIG. 2 shows an example of the system configuration of the lesion suspected region classification unit 23 and the lesion likelihood setting unit 24. As shown in the upper part of the figure, the lesion suspected region classification unit 23 receives the input from the input device 10 and defines the lesion suspected region classification formula, and the lesion suspected region classification formula definition unit 30 and the lesion suspicion that stores the formula. It is composed of a region classification type storage unit 31 and a lesion suspected area classification execution unit 32 that classifies the input lesion suspected area according to the lesion suspected area classification formula.

図2の下部に示すように病変尤度設定部24は、分類された病変疑い領域の情報から、分類された群ごとに病変尤度算出式を定義する病変尤度算出式定義部33と、その式を記憶する病変尤度算出式記憶部34と、病変尤度算出式に従って入力病変疑い領域の病変尤度を算出する病変尤度算出部35とから構成されている。 As shown in the lower part of FIG. 2, the lesion likelihood setting unit 24 includes a lesion likelihood calculation formula definition unit 33 that defines a lesion likelihood calculation formula for each classified group from information on the classified suspected lesion region. It is composed of a lesion likelihood calculation formula storage unit 34 that stores the formula and a lesion likelihood calculation unit 35 that calculates the lesion likelihood of the input lesion suspected region according to the lesion likelihood calculation formula.

次に、図1、図2に示したシステム構成を持つ本実施例の医用画像処理装置による処理の流れを説明する。ここでは、システムからの入力もしくはユーザからの指示を受けた医用画像処理装置11が、病変疑い領域分類式記憶部31に病変疑い領域分類式を記憶する病変疑い領域分類式定義処理と、システムからの入力もしくはユーザからの指示を受けた医用画像処理装置11が、病変尤度算出式記憶部34に病変尤度算出式を記憶する病変尤度算出式定義処理と、に分割してそれぞれ説明する。 Next, the flow of processing by the medical image processing apparatus of this embodiment having the system configuration shown in FIGS. 1 and 2 will be described. Here, the medical image processing device 11 that receives the input from the system or the instruction from the user stores the lesion suspected area classification formula in the lesion suspected area classification formula storage unit 31 and the lesion suspected region classification formula definition processing from the system. The medical image processing device 11 that receives the input of the above or the instruction from the user divides into the lesion likelihood calculation formula definition process that stores the lesion likelihood calculation formula in the lesion likelihood calculation formula storage unit 34, and each of them will be described separately. ..

まず、図3の処理フローを用いて、病変疑い領域分類式定義処理について説明する。病変疑い情報表示部22は、システムからの入力もしくはユーザからの指示により、医用画像記憶部20と病変疑い情報記憶部21から、複数の医用画像とそれらに対応する病変疑い情報を取得し、医用画像中の病変疑い領域の複数の特徴量を各々算出し、それを直感的にわかりやすい形でモニタ12に表示する(ステップ101)。 First, the lesion suspected region classification formula definition process will be described using the process flow of FIG. The lesion suspicion information display unit 22 acquires a plurality of medical images and corresponding lesion suspicion information from the medical image storage unit 20 and the lesion suspicion information storage unit 21 by input from the system or instructions from the user, and is used for medical purposes. Each of the plurality of feature quantities of the suspected lesion region in the image is calculated and displayed on the monitor 12 in an intuitive and easy-to-understand format (step 101).

次に病変疑い領域分類部23の病変疑い領域分類式定義部30は、入力装置10に対するユーザからの入力により、複数の特徴量各々に対する分類の必要性の有無を取得し(ステップ102)、それを元に病変疑い領域分類式を定義し(ステップ103)、病変疑い領域分類式記憶部31に記憶する(ステップ104)。 Next, the lesion suspicious area classification formula definition unit 30 of the lesion suspicious area classification unit 23 acquires whether or not there is a need for classification for each of the plurality of feature quantities by input from the user to the input device 10 (step 102). A lesion suspicious area classification formula is defined based on the above (step 103), and is stored in the lesion suspicious area classification formula storage unit 31 (step 104).

次に、図4、図5を用いて、病変疑い情報表示部22による病変疑い領域、及び特徴量の表示の方法の一例について説明する。図4の(a)は、各病変疑い領域の輝度、同図の(b)は各病変疑い領域の大きさを特徴量として計算し、これら特徴量である輝度もしくは大きさのヒストグラムと、ヒストグラムの値に対応する代表的な病変疑い領域の画像を重畳表示したものである。 Next, an example of a method of displaying the suspected lesion region and the feature amount by the lesion suspicion information display unit 22 will be described with reference to FIGS. 4 and 5. FIG. 4A shows the brightness of each suspected lesion area, and FIG. 4B calculates the size of each suspected lesion area as a feature amount, and a histogram of the brightness or size of these feature amounts and a histogram. The image of the representative suspected lesion area corresponding to the value of is superimposed and displayed.

また図5の(a)は、各病変疑い領域と画像内の二次元位置、同図の(b)は各病変疑い領域と被検体である人体のボリュームデータ内の体軸方向位置を特徴量として計算し、人体モデルの断面画像上にそれぞれ重畳表示したものである。図4、図5のような表示画面により、ユーザは特徴量ごとに、どのような特徴を持つ病変疑い領域が、どのくらいの割合で存在するのかを把握することが可能となる。なお、図4、図5中の点線は、後で説明する病変疑い領域を分類するための閾値(thresold)を示している。 Further, (a) of FIG. 5 features each suspected lesion region and a two-dimensional position in the image, and (b) of FIG. 5 features each suspected lesion region and the axial position in the volume data of the human body as a subject. It is calculated as, and superimposed and displayed on the cross-sectional image of the human body model. With the display screens as shown in FIGS. 4 and 5, the user can grasp what kind of characteristic the suspected lesion region exists for each feature amount and at what ratio. The dotted line in FIGS. 4 and 5 indicates a threshold value for classifying the suspected lesion region, which will be described later.

図6は、医用画像処理装置11の病変疑い情報表示部22によるモニタ12への表示から、ユーザが分類の必要性を判断するための表示画面例を示す図である。図6の(a)、(b)、(c)、(d)に示すように、モニタ12の表示画面上で、図4、図5で図示した重畳表示画像の横に、例えばチェックボックス34を表示する。そして、分類が必要と判断される特徴量には、ユーザによりユーザ入力35がなされる。病変疑い領域分類部23の病変疑い領域分類式定義部30は特徴量のチェックの有無を取得し、図6の(e)にその一例を示す分類必要性判断表にその情報を反映することで、分類の必要の有無の情報を取得して病変疑い領域分類式を定義する。 FIG. 6 is a diagram showing an example of a display screen for the user to determine the necessity of classification from the display on the monitor 12 by the lesion suspicion information display unit 22 of the medical image processing device 11. As shown in (a), (b), (c), and (d) of FIG. 6, on the display screen of the monitor 12, for example, a check box 34 beside the superimposed display image shown in FIGS. 4 and 5. Is displayed. Then, the user inputs 35 to the feature amount that is determined to require classification. The lesion suspicious area classification formula definition unit 30 of the lesion suspicious area classification unit 23 acquires the presence or absence of the feature quantity check, and reflects the information in the classification necessity judgment table showing an example in (e) of FIG. , Acquire information on the necessity of classification and define the lesion suspected area classification formula.

ここで、病変疑い領域分類式定義部30の病変疑い領域分類式について、具体例を用いて説明する。この病変疑い領域分類式は、複数の病変疑い領域を、その特徴量によって複数の群に分類するための式である。例えば、複数の病変疑い領域ROI[n](n=0〜N)が、領域の中心位置(Cx[n]、Cy[n])と半径R[n]で表すことができる場合で、図3のステップ102において分類の必要性があると判断された特徴量が、図6の(e) の分類必要性判断表に示すように画像内の二次元位置である場合で、病変疑い領域をG0とG1の二つの群に分類する場合について述べる。 Here, the lesion suspected region classification formula of the lesion suspected region classification formula definition unit 30 will be described with reference to a specific example. This lesion suspected region classification formula is a formula for classifying a plurality of lesion suspected regions into a plurality of groups according to their feature amounts. For example, when a plurality of suspected lesion regions ROI [n] (n = 0 to N) can be represented by the center position (Cx [n], Cy [n]) of the region and the radius R [n]. When the feature amount determined to be in need of classification in step 102 of 3 is a two-dimensional position in the image as shown in the classification necessity determination table of FIG. 6 (e), the lesion suspected region is defined. The case of classifying into two groups of G0 and G1 will be described.

図7にこの場合に定義され、病変疑い領域分類式記憶部31に記憶される病変疑い領域分類式の一例を式1として示した。この式1は、二次元画像内の縦軸が画像下に向かって増加する場合で、病変疑い領域分類実行部32が、画像内で中心位置が閾値(threshold)より上にある病変疑い領域はG0に、閾値(thresold)と同じまたは下にある病変疑い領域はG1に分類することを意味する。この閾値は事前に決めておいてもよいし、病変疑い領域から算出する特徴量(ここではCx[n])の平均値や中央値、もしくは判別分析法などの手法で自動的に求めてもよいし、ユーザが調整して決定するものとしてもよい。図4、5では閾値を点線で示しているが、この点線をユーザが画面上で動かせるようなユーザインタフェース(UI)にしておけば、ユーザがモニタ12の表示画面を使って簡便に調整することができる。 FIG. 7 shows an example of the lesion suspected region classification formula defined in this case and stored in the lesion suspected region classification formula storage unit 31 as Equation 1. In this equation 1, the vertical axis in the two-dimensional image increases toward the bottom of the image, and the lesion suspected area classification execution unit 32 has the center position in the image above the threshold value (threshold). At G0, a suspected lesion area equal to or below the threshold is meant to be classified as G1. This threshold value may be determined in advance, or may be automatically obtained by the average value or median value of the feature amount (here, Cx [n]) calculated from the suspected lesion region, or by a method such as a discriminant analysis method. Alternatively, it may be adjusted and decided by the user. In FIGS. 4 and 5, the threshold value is shown by a dotted line, but if the user interface (UI) is set so that the user can move the dotted line on the screen, the user can easily adjust the threshold value using the display screen of the monitor 12. Can be done.

次に、図8を用いて、病変尤度設定部24の病変尤度算出式定義部33における病変尤度算出式定義処理について説明する。医用画像処理装置11は、システムからの入力もしくはユーザからの指示により、病変疑い情報記憶部21と医用画像記憶部20とから、一つ以上の病変疑い領域とその領域に対応する医用画像を取得して、画像情報付き病変疑い領域を生成する(ステップ201)。 Next, with reference to FIG. 8, the lesion likelihood calculation formula definition process in the lesion likelihood calculation formula definition unit 33 of the lesion likelihood setting unit 24 will be described. The medical image processing device 11 acquires one or more suspected lesion regions and medical images corresponding to the regions from the suspected lesion information storage unit 21 and the medical image storage unit 20 by input from the system or instructions from the user. Then, a suspected lesion area with image information is generated (step 201).

病変疑い領域分類実行部32は、図7に示したような病変疑い領域分類式に従って、これらの医用画像情報付き病変疑い領域を複数の群に分類する(ステップ202)。 The lesion suspected region classification execution unit 32 classifies these lesion suspected region with medical image information into a plurality of groups according to the lesion suspected region classification formula as shown in FIG. 7 (step 202).

次に病変尤度算出式定義部33は、群に分類された医用画像情報付き病変疑い領域について画像特徴を示す特徴量を把握し, 分類された群ごとに病変尤度算出式を定義し(ステップ203),病変尤度算出式記憶部32に記憶する(ステップ204)。 Next, the lesion likelihood calculation formula definition unit 33 grasps the feature amount indicating the image feature for the lesion suspected region with medical image information classified into groups, and defines the lesion likelihood calculation formula for each classified group ( Step 203), store in the lesion likelihood calculation formula storage unit 32 (step 204).

病変尤度算出式定義部33が定義する病変尤度算出式について、具体例を用いて説明する。病変尤度算出式は、予め病変尤度が判明している病変疑い領域を用いて、病変疑い領域分類処理にて分類された複数の群ごとに、疑い領域の病変尤度Lprev[n]を算出する式である。複数の病変疑い領域ROI[n]が、G0、G1の二つの群に分類されている場合で、群G0について説明する。 The lesion likelihood calculation formula defined by the lesion likelihood calculation formula definition unit 33 will be described with reference to a specific example. In the lesion likelihood calculation formula, the lesion likelihood Lprev [n] of the suspected region is calculated for each of a plurality of groups classified by the lesion suspect region classification process using the lesion likelihood region whose lesion likelihood is known in advance. It is a formula to calculate. In the case where a plurality of suspected lesion regions ROI [n] are classified into two groups of G0 and G1, the group G0 will be described.

ここでG0に分類された病変疑い領域をROI_0[n(0〜N_0)]とし、ROI_0[n]は予め病変尤度が判明しており、その確定病変尤度L[n]は0か1かの2値で表されている場合とする。確定病変尤度L[n]は、L[n]=1が医師の診断により病変であった病変疑い領域、L[n]=0がそれ以外、即ち病変ではなかった病変疑い領域を表す。 Here, the suspected lesion region classified as G0 is set to ROI_0 [n (0 to N_0)], the lesion likelihood of ROI_0 [n] is known in advance, and the definite lesion likelihood L [n] is 0 or 1. It is assumed that it is represented by the two values. The definite lesion likelihood L [n] represents a lesion suspected region where L [n] = 1 was a lesion according to a doctor's diagnosis, and L [n] = 0 represents another, that is, a lesion suspected region which was not a lesion.

図9にこの場合の病変尤度算出式の一例を式2として示す。式2でπは円周率、R[n]は病変疑い領域ROI[n]の半径、P[n]は病変疑い領域ROI[n]の平均輝度、AmaxはROI[n]がとり得る面積の最大値、PmaxはROI[n]がとり得る平均輝度の最大値を示す。この式2に従うと、病変疑い領域は、面積が大きく、輝度が高い方が病変らしいという結果となる。 FIG. 9 shows an example of the lesion likelihood calculation formula in this case as Equation 2. In Equation 2, π is the pi, R [n] is the radius of the suspected lesion region ROI [n], P [n] is the average brightness of the suspected lesion region ROI [n], and Amax is the area that ROI [n] can take. The maximum value of, Pmax, indicates the maximum value of the average brightness that ROI [n] can take. According to this equation 2, the suspected lesion area has a larger area and higher brightness, which is more likely to be a lesion.

この式2では病変尤度Lprev[n]は連続値を取るが、閾値を用いて多値化したものを最終的な病変尤度とすることもできる。ここでは説明のため、単純な病変尤度算出式を例に挙げたが、この式の設計は所謂ルールベースにより、真に病変であった病変尤度疑い領域とそうでなかった病変尤度疑い領域の特徴量等の違いを判断し、それぞれのとり得る範囲から人手で設計することもできるし、SVM(Support Vector Machine)やDeepLearning等、機械学習を用いて設計することもできる。 In this equation 2, the lesion likelihood Lprev [n] takes a continuous value, but a multi-valued one using a threshold value can be used as the final lesion likelihood. Here, for the sake of explanation, a simple lesion likelihood calculation formula is given as an example, but the design of this formula is based on the so-called rule base, and the lesion likelihood region that was truly a lesion and the lesion likelihood suspicion that was not. It is possible to judge the difference in the feature amount of the region and manually design from each possible range, or it is possible to design by using machine learning such as SVM (Support Vector Machine) or Deep Learning.

以上説明した医用画像処理装置11の処理により、参考にする過去の画像の量が不十分である場合にも、高精度な病変の疑わしさの定義式すなわち、病変尤度算出式の設定が可能になる。 By the processing of the medical image processing apparatus 11 described above, it is possible to set a highly accurate definition formula of suspicion of a lesion, that is, a lesion likelihood calculation formula even when the amount of past images to be referred to is insufficient. become.

次に、図10の病変尤度算出処理フローを用いて、こうして設定され、病変尤度算出式記憶部34に記憶された病変尤度算出式により病変尤度を算出する場合の処理、すなわちシステムからの入力もしくはユーザからの指示を受けた医用画像処理装置11が、病変尤度設定部24により病変尤度を算出する病変尤度算出処理について説明する。 Next, using the lesion likelihood calculation processing flow of FIG. 10, the process in which the lesion likelihood is calculated by the lesion likelihood calculation formula set in this way and stored in the lesion likelihood calculation formula storage unit 34, that is, the system. The lesion likelihood calculation process for calculating the lesion likelihood by the lesion likelihood setting unit 24 by the medical image processing apparatus 11 receiving the input from or instructed by the user will be described.

医用画像処理装置11は、システムからの入力もしくはユーザからの指示により、病変疑い情報記憶部21と医用画像記憶部20とから、処理対象となる病変疑い領域とその領域に対応する医用画像を取得して、画像情報付き病変疑い領域とする(ステップ301)。 The medical image processing device 11 acquires a suspected lesion area to be processed and a medical image corresponding to the area from the suspected lesion information storage unit 21 and the medical image storage unit 20 by input from the system or an instruction from the user. Then, it is set as a suspected lesion area with image information (step 301).

次に病変疑い領域分類実行部32が、あらかじめ設定された医用画像分類式に従って、処理対象の画像情報付き病変疑い領域が属する郡を決定する(ステップ302)。 Next, the lesion suspicious area classification execution unit 32 determines the group to which the lesion suspicious area with image information to be processed belongs according to a preset medical image classification formula (step 302).

次に病変尤度算出部35は、画像情報付き病変疑い領域が所属する群に対応する病変尤度算出式に従い、病変尤度を算出する(ステップ303)。なお、図9に例示した病変尤度算出式は、病変疑い領域分類部23で分類する群の数だけ存在する。ステップ303では、これら複数の算出式により得られる複数の病変尤度を、論理和もしくは論理積等により1つの病変尤度とした後に出力するものとする。 Next, the lesion likelihood calculation unit 35 calculates the lesion likelihood according to the lesion likelihood calculation formula corresponding to the group to which the lesion suspicious region with image information belongs (step 303). The lesion likelihood calculation formula illustrated in FIG. 9 exists as many as the number of groups classified by the lesion suspicious region classification unit 23. In step 303, the plurality of lesion likelihoods obtained from these plurality of calculation formulas are converted into one lesion likelihood by logical sum or logical product, and then output.

ここで、上述した本実施例の医用画像処理装置11を利用することで解決される課題の具体例について説明する。まず、図11を用いて、現状の医用画像処理装置に発生する課題を説明する。 Here, a specific example of the problem solved by using the medical image processing apparatus 11 of the present embodiment described above will be described. First, problems that occur in the current medical image processing apparatus will be described with reference to FIG.

同図において、丸や三角、四角などの図形は、学習データや評価データの病変疑い領域を模式的に示している。図11の(a)に示す学習データにおいて、点線の丸で示す病変疑い領域が,真に病変である領域とする。この学習データにおいて点線の丸で示された領域とそうでない領域の違いを考えると、“形状は丸で、色は黒”である領域のみが、真に病変であると捉えることができる。従って、形状は丸で、色は黒である領域の結節尤度が高い領域であるとする式が設計されると考えられる。 In the figure, figures such as circles, triangles, and squares schematically show lesion-suspicious areas of learning data and evaluation data. In the learning data shown in FIG. 11A, the suspected lesion region indicated by the dotted circle is defined as the region that is truly a lesion. Considering the difference between the area indicated by the dotted circle and the area not indicated by the dotted line in this training data, it can be considered that only the area where the shape is a circle and the color is black is truly a lesion. Therefore, it is considered that an equation is designed in which the region having a round shape and a black color has a high nodule likelihood.

この病変尤度算出式を用い、図11の(b)に示す評価データの病変尤度を算出した場合、“形状は丸で色は黒”となる領域の病変尤度が高いという結果となる。この結果を、識別器による評価として示す。一方で、医師による評価として示した図は、同じ評価データを医師が評価した場合を想定した結果である。このように、識別器による評価とは結果が異なり、“形状は丸で色は黒”以外の領域にも点線の丸で示される領域がある。このように、現状の技術では、学習データのバリエーションが十分でない場合、結節の特徴を正しく捉えた病変尤度算出式が設計できない可能性が高い。 When the lesion likelihood of the evaluation data shown in FIG. 11 (b) is calculated using this lesion likelihood calculation formula, the result is that the lesion likelihood in the region where the shape is round and the color is black is high. .. This result is shown as an evaluation by a classifier. On the other hand, the figure shown as an evaluation by a doctor is a result assuming that the same evaluation data is evaluated by a doctor. As described above, the result is different from the evaluation by the discriminator, and there is a region indicated by a dotted circle in the region other than "the shape is a circle and the color is black". As described above, with the current technology, if the variation of the learning data is not sufficient, there is a high possibility that a lesion likelihood calculation formula that correctly captures the characteristics of the nodule cannot be designed.

次に図12を用い、本実施例の医用画像処理装置11を利用した場合について説明する。図12の(a)、(b)は、図11の場合と同じ学習データ、評価データを用いているものとする。また、ここで説明する例は、図7の式1で説明した病変疑い領域分類式を用いるものとする。図12の(a)の学習データにある点線は、図7で説明した閾値(thresold)を示す。従って、G0は閾値より上の病変疑い領域の集合、G1は閾値より下の病変疑い領域の集合となる。G0、G1それぞれに病変尤度算出式を設計すると、G0の場合は“形状が丸”、G1の場合は“色は黒”の病変疑い領域の病変尤度が高いという結果となると考えられる。 Next, a case where the medical image processing apparatus 11 of this embodiment is used will be described with reference to FIG. It is assumed that (a) and (b) of FIG. 12 use the same learning data and evaluation data as in the case of FIG. Further, in the example described here, the lesion suspected region classification formula described in the formula 1 of FIG. 7 is used. The dotted line in the training data of FIG. 12A indicates the threshold value described in FIG. 7. Therefore, G0 is a set of suspected lesion regions above the threshold value, and G1 is a set of suspected lesion regions below the threshold value. When the lesion likelihood calculation formulas are designed for each of G0 and G1, it is considered that the result is that the lesion likelihood of the suspected lesion region is high in the case of G0, which is “round shape” and in the case of G1, which is “black in color”.

ここで、G0により設計された病変尤度算出式を用いる病変尤度算出部35を識別器0、G1により設計された病変尤度算出式を用いる病変尤度算出部35を識別器1とする。図12の(b)の左部に識別部0、識別器1それぞれで想定される評価データに対する評価結果を示す。ここでは、複数の識別器による結果の論理和を、最終的な出力とする。最終的な識別結果を図12の(b)の右部に示す。こうして、本実施例の医用画像処理装置により、病変疑い領域の画像内二次元位置、ここでは画像の上部と下部とで別の病変尤度算出式を設定すべきというユーザの知見を取り入れることで、ユーザによる識別結果と一致する精度の高い識別結果を取得することができる。 Here, the lesion likelihood calculation unit 35 using the lesion likelihood calculation formula designed by G0 is referred to as the discriminator 0, and the lesion likelihood calculation unit 35 using the lesion likelihood calculation formula designed by G1 is referred to as the discriminator 1. .. The left part of FIG. 12B shows the evaluation results for the evaluation data assumed by the discriminator 0 and the discriminator 1. Here, the logical sum of the results obtained by a plurality of classifiers is used as the final output. The final identification result is shown in the right part of FIG. 12 (b). In this way, by incorporating the user's knowledge that the medical image processing apparatus of this embodiment should set different lesion likelihood calculation formulas for the two-dimensional position in the image of the suspected lesion region, here, the upper part and the lower part of the image. , It is possible to obtain a highly accurate identification result that matches the identification result by the user.

以上詳述した実施例1の医用画像処理装置によれば、病変疑い陰影の自動検出において、病変であると既に診断を受けている過去の画像の情報を参考にする場合に,参考にする過去の画像の量が不十分である場合にも、高精度な病変の疑わしさの定義式の設定が可能な医用画像処理装置、方法、及びプログラムを提供することができる。特に、従来の技術による定義式の設定では病変の検出率が低い場合、例えば放射線科医が読影した場合の結節検出率6割程度を下回る場合に有効な技術を提供することができる。 According to the medical image processing apparatus of Example 1 described in detail above, the past to be referred to when referring to the information of the past image that has already been diagnosed as a lesion in the automatic detection of the suspected lesion shadow. It is possible to provide a medical image processing device, a method, and a program capable of setting a highly accurate definition formula of suspicion of a lesion even when the amount of images in the image is insufficient. In particular, it is possible to provide an effective technique when the detection rate of lesions is low in the setting of the definition formula by the conventional technique, for example, when the nodule detection rate is less than about 60% when the radiologist interprets the image.

実施例1の医用画像処理装置11の処理技術は、CADの稼動前に精度の初期調整として利用することもできる。この場合、病変疑い領域分類式および病変尤度算出式は、経験的に妥当と考えられる精度に予め設定しておき、CAD納入時に、そのときの利用者のニーズや施設の読影方針等に応じて調整するという使い方を想定している。初期値となる病変疑い領域分類式および病変尤度算出式の設定方法は、手動で式を設定することもできるし、利用可能な医用画像データベースを用いて、病変疑い領域分類式定義処理および病変尤度算出式定義処理を実行して設定することもできる。 The processing technique of the medical image processing apparatus 11 of the first embodiment can also be used as an initial adjustment of accuracy before the operation of CAD. In this case, the lesion suspected area classification formula and the lesion likelihood calculation formula are set in advance to the accuracy considered to be empirically appropriate, and at the time of CAD delivery, according to the needs of the user at that time, the interpretation policy of the facility, etc. It is supposed to be used for adjustment. The initial values of the suspected lesion area classification formula and the lesion likelihood calculation formula can be set manually, or the available medical image database can be used to define the suspected lesion area classification formula and the lesion. It can also be set by executing the likelihood calculation formula definition process.

ここで調整前、即ち出荷時に設定する病変疑い領域分類式定義処理および病変尤度算出式定義処理は、納入先の施設にて取得された医用画像データベースを用いて実施することもできるし、もしくは納入先の施設に近い読影方針を持つ別の施設にて取得された医用画像データベースの利用が可能であれば、それを利用してもよい。もしくは、納入先の施設に近い読影方針を持つ別の施設にて調整された結果、病変疑い領域分類式およびの病変尤度算出式を利用することとしてもよい。 Here, the lesion suspicious area classification formula definition process and the lesion likelihood calculation formula definition process set before adjustment, that is, at the time of shipment, can be performed using the medical image database acquired at the delivery destination facility, or can be performed. If it is possible to use a medical image database acquired at another facility that has an image interpretation policy close to that of the delivery destination, it may be used. Alternatively, as a result of adjustment at another facility having an image interpretation policy close to the delivery destination facility, the lesion suspected area classification formula and the lesion likelihood calculation formula may be used.

このように上述してきた処理技術を初期調整として利用する場合、確信度の提示や調整に用いられる医用画像としては、過去に撮影された臨床データを利用することができる。この場合、病変疑い領域情報の種類によって予め分類しておき、確信度算出式に用いる特徴量が顕著に異なるような病変疑い陰影を示すことで、より効率的な調整が行うことも可能である。もしくは、ここで用いる医用画像は、必ずしも実際の臨床画像でなくてもよく、例えば人体に近い形状と素材を持つファントムを撮影したファントムデータを用いることや、臨床画像から模擬的に病変疑い位置を複製・生成した模擬データを利用することもできる。 When the processing technique described above is used as the initial adjustment, clinical data taken in the past can be used as the medical image used for presenting and adjusting the certainty. In this case, it is possible to perform more efficient adjustment by classifying in advance according to the type of lesion suspected area information and showing a lesion suspected shadow such that the feature amount used in the certainty calculation formula is significantly different. .. Alternatively, the medical image used here does not necessarily have to be an actual clinical image. For example, phantom data obtained by photographing a phantom having a shape and material close to the human body can be used, or a lesion suspected position can be simulated from the clinical image. It is also possible to use the simulated data duplicated and generated.

また、CADの稼動前の初期調整のみではなく、CAD稼働中に、随時実施する調整に利用することも可能である。上記のように初期調整を実施する場合、利用するデータとしては過去画像、ファントム画像、模擬画像等を利用するが、稼働中に随時実施する調整に利用する場合、当該施設の撮影装置や撮影方針に従って撮影した画像を用いた調整が可能となり、よりユーザの希望に近い精度での検出が可能になるという効果がある。 Further, it can be used not only for the initial adjustment before the operation of the CAD but also for the adjustment to be performed at any time during the operation of the CAD. When performing initial adjustment as described above, past images, phantom images, simulated images, etc. are used as the data to be used, but when using for adjustment to be performed at any time during operation, the shooting device and shooting policy of the facility concerned. It is possible to make adjustments using the images taken according to the above, and there is an effect that detection can be performed with an accuracy closer to the user's desire.

なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。先に述べたように、本発明はCT医用画像撮影装置により得られる再構成三次元医用画像のみならず、MRI撮影装置等により得られ他の医用画像データに対しても応用可能である。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。実施例1には、医用画像処理装置11に医用画像撮影装置(図示せず)を含まなかったが、医用画像処理装置11は医用画像撮影装置を含んでもよく、また医用画像処理装置11は医用画像撮影装置の一部として機能してもよい。 The present invention is not limited to the above-described embodiment, and includes various modifications. For example, the above-described embodiment has been described in detail in order to explain the present invention in an easy-to-understand manner, and is not necessarily limited to the one including all the described configurations. As described above, the present invention can be applied not only to reconstructed three-dimensional medical images obtained by a CT medical imaging device but also to other medical image data obtained by an MRI imaging device or the like. Further, it is possible to replace a part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment. Further, it is possible to add / delete / replace a part of the configuration of each embodiment with another configuration. In Example 1, the medical image processing device 11 does not include a medical image taking device (not shown), but the medical image processing device 11 may include a medical image taking device, and the medical image processing device 11 is for medical use. It may function as part of an imaging device.

更に、上述した各構成、機能、医用画像処理装置等は、それらの一部又は全部を実現する処理プログラムを作成する例を説明したが、それらの一部又は全部を例えば集積回路で設計する等によりハードウェアで実現しても良いことは言うまでもない。すなわち、画像処理装置の全部または一部の機能は、プログラムに代え、例えば、ASIC(Application Specific Integrated Circuit)、FPGA(Field Programmable Gate Array)などの集積回路などにより実現してもよい。 Further, for each of the above-mentioned configurations, functions, medical image processing devices, etc., an example of creating a processing program that realizes a part or all of them has been described, but a part or all of them is designed, for example, by an integrated circuit. Needless to say, it may be realized by hardware. That is, all or part of the functions of the image processing apparatus may be realized by, for example, an integrated circuit such as an ASIC (Application Specific Integrated Circuit) or an FPGA (Field Programmable Gate Array) instead of the program.

以上詳述した明細書の記載から明らかなように、本明細書中には特許請求の範囲に規定した発明以外に、数多くの発明が含まれている。それらの一部を下記に列記する。 As is clear from the description of the specification described in detail above, a large number of inventions are included in the present specification in addition to the inventions defined in the claims. Some of them are listed below.

[列記1]
画像処理装置により医用画像を処理する画像処理方法であって、
前記画像処理装置は、
前記医用画像内の複数の病変疑い領域と、前記病変疑い領域の複数の特徴量をモニタに表示し、
表示された前記複数の特徴量から選択された特徴量により、前記病変疑い領域を分類し、
前記特徴量は、前記病変疑い領域の輝度と大きさを含む、
ことを特徴とする画像処理方法。
[List 1]
An image processing method for processing medical images with an image processing device.
The image processing device is
A plurality of suspected lesion areas in the medical image and a plurality of feature quantities of the suspected lesion area are displayed on the monitor.
The suspected lesion area is classified according to the feature amount selected from the displayed plurality of feature amounts.
The feature amount includes the brightness and size of the suspected lesion area.
An image processing method characterized by that.

[列記2]
列記1に記載の画像処理方法であって、
前記画像処理装置は、
前記モニタに前記特徴量のヒストグラムを表示する、
ことを特徴とする画像処理方法。
[List 2]
The image processing method according to the list 1.
The image processing device is
A histogram of the feature amount is displayed on the monitor.
An image processing method characterized by that.

[列記3]
列記1に記載の画像処理方法であって、
前記特徴量は、前記病変疑い領域の前記医用画像内の二次元位置を含む、
ことを特徴とする画像処理方法。
[List 3]
The image processing method according to the list 1.
The feature amount includes the two-dimensional position of the suspected lesion region in the medical image.
An image processing method characterized by that.

[列記4]
列記1に記載の画像処理方法であって、
前記特徴量は、前記病変疑い領域のボリュームデータ内の被検体の体軸方向位置を含む、
ことを特徴とする画像処理方法。
[List 4]
The image processing method according to the list 1.
The feature amount includes the body axis position of the subject in the volume data of the suspected lesion region.
An image processing method characterized by that.

[列記5]
画像処理装置に医用画像を処理させる画像処理プログラムであって、
前記画像処理装置が、
前記医用画像内の複数の病変疑い領域と、前記病変疑い領域の複数の特徴量をモニタに表示し、
表示された前記複数の特徴量から選択された特徴量により、前記病変疑い領域を分類する、よう動作させ、
前記特徴量に、前記病変疑い領域の輝度と大きさを含む、
ことを特徴とする画像処理プログラム。
[List 5]
An image processing program that causes an image processing device to process medical images.
The image processing device
A plurality of suspected lesion areas in the medical image and a plurality of feature quantities of the suspected lesion area are displayed on the monitor.
The suspected lesion area is classified according to the feature amount selected from the displayed plurality of feature amounts.
The feature amount includes the brightness and size of the suspected lesion area.
An image processing program characterized by this.

[列記6]
列記5に記載の画像処理プログラムであって、
前記画像処理装置が、
前記モニタに前記特徴量のヒストグラムを表示する、よう動作させる、
ことを特徴とする画像処理プログラム。
[List 6]
The image processing program described in Listing 5
The image processing device
Operate to display the histogram of the feature amount on the monitor.
An image processing program characterized by this.

[列記7]
列記5に記載の画像処理プログラムであって、
前記特徴量は、前記病変疑い領域の前記医用画像内の二次元位置を含む、
ことを特徴とする画像処理プログラム。
[List 7]
The image processing program described in Listing 5
The feature amount includes the two-dimensional position of the suspected lesion region in the medical image.
An image processing program characterized by this.

[列記8]
列記5に記載の画像処理プログラムであって、
前記特徴量は、前記病変疑い領域のボリュームデータ内の被検体の体軸方向位置を含む、
ことを特徴とする画像処理プログラム。
[List 8]
The image processing program described in Listing 5
The feature amount includes the body axis position of the subject in the volume data of the suspected lesion region.
An image processing program characterized by this.

[列記8]
列記5に記載の画像処理プログラムであって、
前記画像処理装置が、
分類された前記病変疑い領域の群ごとに、当該群に属する前記病変疑い領域の病変尤度を算出する、よう動作させる
ことを特徴とする画像処理プログラム。
[List 8]
The image processing program described in Listing 5
The image processing device
An image processing program characterized in that the lesion likelihood of the suspected lesion region belonging to the group is calculated for each group of the suspected lesion region classified.

[列記9]
列記8に記載の画像処理プログラムであって、
前記画像処理装置が、
選択された前記特徴量を用いて病変疑い領域分類式を定義し、当該病変疑い領域分類式に基づき、複数の前記病変疑い領域を前記群に分類する、よう動作させる
ことを特徴とする画像処理プログラム。
[List 9]
The image processing program described in Listing 8
The image processing device
Image processing characterized in that a lesion suspected region classification formula is defined using the selected feature amount, and a plurality of the lesion suspected regions are classified into the group based on the lesion suspected region classification formula. program.

[列記10]
列記8に記載の画像処理プログラムであって、
前記画像処理装置が、、
前記群に分類された前記病変疑い領域の前記特徴量に基づき、前記群ごとに病変尤度算出式を定義し、当該病変尤度算出式に基づき、前記病変疑い領域の病変尤度を算出する、よう動作させる
ことを特徴とする画像処理プログラム。
[List 10]
The image processing program described in Listing 8
The image processing device ...
A lesion likelihood calculation formula is defined for each group based on the characteristic amount of the lesion suspected region classified into the group, and the lesion likelihood of the lesion suspected region is calculated based on the lesion likelihood calculation formula. An image processing program characterized by operating like,.

10 入力装置
11 医用画像処理装置
12 モニタ
20 医用画像記憶部
21 病変疑い情報記憶部
22 病変疑い情報表示部
23 病変疑い領域分類部
24 病変尤度算出部
30 病変疑い領域分類式定義部
31 病変疑い領域分類式記憶部
32 病変尤度算出式定義部
33 病変尤度算出式記憶部
34 チェックボックス
35 ユーザ入力
10 Input device 11 Medical image processing device 12 Monitor 20 Medical image storage unit 21 Lesion suspicious information storage unit 22 Lesion suspicious information display unit 23 Lesion suspicious area classification unit 24 Lesion likelihood calculation unit 30 Lesion suspicious area classification formula definition unit 31 Lesion suspicion Area classification formula Storage unit 32 Lesion likelihood calculation formula Definition unit 33 Lesion likelihood calculation formula Storage unit 34 Check box 35 User input

Claims (7)

医用画像を処理する画像処理装置であって、
前記医用画像内の複数の病変疑い領域と、前記病変疑い領域の複数の特徴量を表示する病変疑い情報表示部と、
表示された前記複数の特徴量から選択された特徴量により、前記病変疑い領域を分類する病変疑い領域分類部と、
前記病変疑い領域分類部により分類された群ごとに、当該群に属する前記病変疑い領域の病変尤度を算出する病変尤度設定部と、を備える、
ことを特徴とする画像処理装置。
An image processing device that processes medical images.
A plurality of suspected lesion areas in the medical image, a lesion suspected information display unit for displaying a plurality of feature quantities of the suspected lesion area, and a lesion information display unit.
A lesion suspected area classification unit that classifies the lesion suspected area according to the feature amount selected from the displayed plurality of feature amounts, and a lesion suspected area classification unit.
For each group classified by the suspected lesion region classification unit, a lesion likelihood setting unit for calculating the lesion likelihood of the suspected lesion region belonging to the group is provided.
An image processing device characterized by this.
請求項1に記載の画像処理装置であって、
前記特徴量は、前記病変疑い領域の輝度と大きさを含む、
ことを特徴とする画像処理装置。
The image processing apparatus according to claim 1.
The feature amount includes the brightness and size of the suspected lesion area.
An image processing device characterized by this.
請求項2に記載の画像処理装置であって、
前記病変疑い情報表示部は、前記特徴量のヒストグラムを表示する、
ことを特徴とする画像処理装置。
The image processing apparatus according to claim 2.
The lesion suspicion information display unit displays a histogram of the feature amount.
An image processing device characterized by this.
請求項2に記載の画像処理装置であって、
前記特徴量は、前記病変疑い領域の前記医用画像内の二次元位置を含む、
ことを特徴とする画像処理装置。
The image processing apparatus according to claim 2.
The feature amount includes the two-dimensional position of the suspected lesion region in the medical image.
An image processing device characterized by this.
請求項2に記載の画像処理装置であって、
前記特徴量は、前記病変疑い領域のボリュームデータ内の被検体の体軸方向位置を含む、
ことを特徴とする画像処理装置。
The image processing apparatus according to claim 2.
The feature amount includes the body axis position of the subject in the volume data of the suspected lesion region.
An image processing device characterized by this.
画像処理装置により医用画像を処理する画像処理方法であって、
前記画像処理装置は、
前記医用画像内の複数の病変疑い領域と、前記病変疑い領域の複数の特徴量をモニタに表示し、
表示された前記複数の特徴量から選択された特徴量により、前記病変疑い領域を分類し、
分類された前記病変疑い領域の群ごとに、当該群に属する前記病変疑い領域の病変尤度を算出する、
ことを特徴とする画像処理方法。
An image processing method for processing medical images with an image processing device.
The image processing device is
A plurality of suspected lesion areas in the medical image and a plurality of feature quantities of the suspected lesion area are displayed on the monitor.
The suspected lesion area is classified according to the feature amount selected from the displayed plurality of feature amounts.
For each group of the suspected lesion regions classified, the lesion likelihood of the suspected lesion region belonging to the group is calculated.
An image processing method characterized by that.
画像処理装置に医用画像を処理させる画像処理プログラムであって、
前記画像処理装置が、
前記医用画像内の複数の病変疑い領域と、前記病変疑い領域の複数の特徴量をモニタに表示し、
表示された前記複数の特徴量から選択された特徴量により、前記病変疑い領域を分類し、
分類された前記病変疑い領域の群ごとに、当該群に属する前記病変疑い領域の病変尤度を算出する、よう動作させる、
ことを特徴とする画像処理プログラム。
An image processing program that causes an image processing device to process medical images.
The image processing device
A plurality of suspected lesion areas in the medical image and a plurality of feature quantities of the suspected lesion area are displayed on the monitor.
The suspected lesion area is classified according to the feature amount selected from the displayed plurality of feature amounts.
For each group of the suspected lesion regions classified, the lesion likelihood of the suspected lesion region belonging to the group is calculated .
An image processing program characterized by this.
JP2018030552A 2018-02-23 2018-02-23 Image processing equipment, methods and programs Active JP6987342B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018030552A JP6987342B2 (en) 2018-02-23 2018-02-23 Image processing equipment, methods and programs

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018030552A JP6987342B2 (en) 2018-02-23 2018-02-23 Image processing equipment, methods and programs

Publications (2)

Publication Number Publication Date
JP2019141478A JP2019141478A (en) 2019-08-29
JP6987342B2 true JP6987342B2 (en) 2021-12-22

Family

ID=67773070

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018030552A Active JP6987342B2 (en) 2018-02-23 2018-02-23 Image processing equipment, methods and programs

Country Status (1)

Country Link
JP (1) JP6987342B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110969613B (en) * 2019-12-06 2023-04-14 广州柏视医疗科技有限公司 Intelligent pulmonary tuberculosis identification method and system with image symptom explanation
CN111709932B (en) * 2020-06-17 2023-04-07 广州麦仑信息科技有限公司 Method for detecting paper defects from coarse to fine

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4253055B2 (en) * 1998-04-27 2009-04-08 東芝医用システムエンジニアリング株式会社 Computer-aided diagnosis device
US7783089B2 (en) * 2002-04-15 2010-08-24 General Electric Company Method and apparatus for providing mammographic image metrics to a clinician
JP2008253292A (en) * 2007-03-30 2008-10-23 Fujifilm Corp Case image search apparatus and system
JP5536987B2 (en) * 2008-05-16 2014-07-02 株式会社東芝 Image processing device

Also Published As

Publication number Publication date
JP2019141478A (en) 2019-08-29

Similar Documents

Publication Publication Date Title
EP2710958B1 (en) Method and system for intelligent qualitative and quantitative analysis of digital radiography softcopy reading
Niemeijer et al. Information fusion for diabetic retinopathy CAD in digital color fundus photographs
KR101874348B1 (en) Method for facilitating dignosis of subject based on chest posteroanterior view thereof, and apparatus using the same
JP5197029B2 (en) Medical image processing device
US7072435B2 (en) Methods and apparatus for anomaly detection
Hogeweg et al. Fusion of local and global detection systems to detect tuberculosis in chest radiographs
JP2004531315A (en) Computer-aided diagnosis method and system for assisting diagnosis of pulmonary nodules in digital volumetric medical images
JP6824845B2 (en) Image processing systems, equipment, methods and programs
JP2016214323A (en) Medical diagnosis support device, information processing method, and program
JP2006034585A (en) Picture display device and picture display method, and program thereof
JP2007151645A (en) Medical diagnostic imaging support system
JP2017534316A (en) Image report annotation identification
JP6987342B2 (en) Image processing equipment, methods and programs
WO2021037790A1 (en) Image processing for stroke characterization
US11058383B2 (en) Apparatus for the detection of opacities in X-ray images
US11423554B2 (en) Registering a two-dimensional image with a three-dimensional image
WO2017221537A1 (en) Image processing device and method
JP2006181025A (en) Abnormal shadow detecting method, device and program
WO2010143015A2 (en) Breast tissue density measure
Bandyopadhyay A Method for Lung Boundary Detection
JP2003310587A (en) Display device for abnormal shadow detected result
Park et al. Separation of left and right lungs using 3D information of sequential CT images and a guided dynamic programming algorithm
JP2010075330A (en) Medical image processor and program
JP2023173115A (en) Image recognition program, image recognition method and image recognition device
WO2017018230A1 (en) Image processing device, method, and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200421

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210303

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210323

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210520

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210928

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20211027

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211028

R150 Certificate of patent or registration of utility model

Ref document number: 6987342

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150