Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7526503B2 - Image processing device and image processing method - Google Patents
[go: Go Back, main page]

JP7526503B2 - Image processing device and image processing method - Google Patents

Image processing device and image processing method Download PDF

Info

Publication number
JP7526503B2
JP7526503B2 JP2022172459A JP2022172459A JP7526503B2 JP 7526503 B2 JP7526503 B2 JP 7526503B2 JP 2022172459 A JP2022172459 A JP 2022172459A JP 2022172459 A JP2022172459 A JP 2022172459A JP 7526503 B2 JP7526503 B2 JP 7526503B2
Authority
JP
Japan
Prior art keywords
pixel
luminance value
unit
interest
image processing
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
JP2022172459A
Other languages
Japanese (ja)
Other versions
JP2024064108A (en
Inventor
聖 正岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Asahi Roentgen Industries Co Ltd
Original Assignee
Asahi Roentgen Industries Co Ltd
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 Asahi Roentgen Industries Co Ltd filed Critical Asahi Roentgen Industries Co Ltd
Priority to JP2022172459A priority Critical patent/JP7526503B2/en
Publication of JP2024064108A publication Critical patent/JP2024064108A/en
Application granted granted Critical
Publication of JP7526503B2 publication Critical patent/JP7526503B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Description

本発明は、X線撮影画像に写っている特定の物質を抽出する技術に関する。 The present invention relates to a technology for extracting specific substances that appear in X-ray images.

X線撮影画像を処理する画像処理装置では、金属アーチファクトの低減などのために、X線撮影画像に写っている金属を正確に抽出することが求められている。 Image processing devices that process X-ray images are required to accurately extract metals that appear in X-ray images in order to reduce metal artifacts, etc.

特開2013-240584号公報JP 2013-240584 A

特許文献1で提案されている画像処理装置は、投影データを基礎として、所定のX線吸収係数を持つ仮想金属体がX線CT撮影領域に擬似的に配置されたときの擬似投影データを演算により取得し、その擬似投影データを利用してX線撮影画像に写っている金属を抽出している。 The image processing device proposed in Patent Document 1 uses projection data as a basis to calculate and obtain pseudo projection data when a virtual metal object with a specified X-ray absorption coefficient is placed in the X-ray CT imaging area, and uses the pseudo projection data to extract the metal that appears in the X-ray image.

しかしながら、特許文献1で提案されている画像処理装置は、擬似投影データを再構成するという処理を行っているため計算時間がかかる。 However, the image processing device proposed in Patent Document 1 requires a long calculation time because it performs a process of reconstructing pseudo-projection data.

また、X線撮影画像に写っている金属以外の物質を抽出することが求められる場合もある。 In addition, there may be cases where it is necessary to extract non-metallic substances that appear in X-ray images.

本発明は、上記の状況に鑑み、X線撮影画像に写っている特定の物質を抽出するときに、抽出精度の向上及び処理時間の短縮を見込むことができる画像処理装置及び画像処理方法を提供することを目的とするものである。 In view of the above, the present invention aims to provide an image processing device and an image processing method that can improve extraction accuracy and reduce processing time when extracting a specific substance that appears in an X-ray image.

上記目的を達成するために本発明に係る画像処理装置は、X線撮影画像において着目ピクセルを囲む所定の領域を設定する設定部と、前記所定の領域を構成する複数の単位領域それぞれの平均輝度値に基づき、暫定的に第1の物質が写っている第1領域とした単位領域全体の代表的輝度値を算出する算出部と、前記代表的輝度値に対する前記着目ピクセルの輝度値の比と閾値との大小関係に基づき、前記着目ピクセルに前記第1の物質が写っているか前記第1の物質以外の物質が写っているかを判定する判定部と、を備える構成(第1の構成)である。 In order to achieve the above object, the image processing device according to the present invention is configured (first configuration) to include a setting unit that sets a predetermined region surrounding a pixel of interest in an X-ray image, a calculation unit that calculates a representative luminance value of the entire unit region that is provisionally set as a first region in which a first substance is captured, based on the average luminance value of each of a plurality of unit regions that make up the predetermined region, and a determination unit that determines whether the pixel of interest captures the first substance or a substance other than the first substance, based on the magnitude relationship between the ratio of the luminance value of the pixel of interest to the representative luminance value and a threshold value.

上記第1の構成の画像処理装置において、前記閾値は、前記代表的輝度値に応じて可変する構成(第2の構成)であってもよい。 In the image processing device of the first configuration, the threshold value may be variable depending on the representative luminance value (second configuration).

上記第1又は第2の構成の画像処理装置において、前記閾値は、前記X線撮影画像が撮影されるときの被写体に対するX線の照射角度に応じて可変する構成(第3の構成)であってもよい。 In the image processing device of the first or second configuration described above, the threshold value may be configured to vary depending on the irradiation angle of X-rays with respect to the subject when the X-ray image is captured (third configuration).

上記目的を達成するために本発明に係る画像処理方法は、X線撮影画像において着目ピクセルを囲む所定の領域を設定する設定ステップと、前記所定の領域を構成する複数の単位領域それぞれの平均輝度値に基づき、暫定的に第1の物質が写っている第1領域とした単位領域全体の代表的輝度値を算出する算出ステップと、前記代表的輝度値に対する前記着目ピクセルの輝度値の比と閾値との大小関係に基づき、前記着目ピクセルに前記第1の物質が写っているか前記第1の物質以外の物質が写っているかを判定する判定ステップと、を備える構成(第4の構成)である。 In order to achieve the above object, the image processing method according to the present invention is configured (fourth configuration) to include a setting step of setting a predetermined region surrounding a pixel of interest in an X-ray image, a calculation step of calculating a representative luminance value of the entire unit region provisionally set as a first region in which a first material is captured, based on the average luminance value of each of a plurality of unit regions constituting the predetermined region, and a determination step of determining whether the pixel of interest captures the first material or a material other than the first material, based on the magnitude relationship between the ratio of the luminance value of the pixel of interest to the representative luminance value and a threshold value.

本発明によると、X線撮影画像に写っている特定の物質を抽出するときに、抽出精度の向上及び処理時間の短縮を見込むことができる。 The present invention is expected to improve extraction accuracy and reduce processing time when extracting specific substances seen in X-ray images.

本発明の一実施形態に係る情報処理装置の構成を示す図FIG. 1 is a diagram showing a configuration of an information processing device according to an embodiment of the present invention. 図1に示す情報処理装置の機能の一例を示す機能ブロック図FIG. 2 is a functional block diagram showing an example of a function of the information processing device shown in FIG. 1 . 着目ピクセル及び所定の領域の一例を示す図FIG. 1 is a diagram showing an example of a pixel of interest and a predetermined region; 着目ピクセル及び所定の領域の他の例を示す図FIG. 13 is a diagram showing another example of a pixel of interest and a predetermined region; 所定の領域における単位領域の平均輝度値分布を示す図FIG. 1 is a diagram showing the average luminance distribution of unit areas in a given area. 暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値と着目ピクセルの輝度値との関係の一例を示す図FIG. 13 is a diagram showing an example of the relationship between the average luminance value of an entire unit region, which is provisionally set as a first region in which a non-metallic object is captured, and the luminance value of a pixel of interest; 金属が写っているピクセルと非金属が写っているピクセルとが異なる輝度値である二値化画像を示す図A binarized image in which pixels representing metals and pixels representing non-metals have different brightness values. アキシャル画像を示す図Axial image diagram 二値化アキシャル画像を示す図A diagram showing a binarized axial image. 順投影画像を示す図Diagram showing forward projection image

本発明の実施形態について図面を参照して以下に説明する。 Embodiments of the present invention are described below with reference to the drawings.

図1は、本発明の一実施形態に係る情報処理装置の構成を示す図である。本発明の一実施形態に係る情報処理装置1(以下、情報処理装置1という)は、制御部2、記憶部3、通信部4、表示部5、及び操作部6を備える。 Figure 1 is a diagram showing the configuration of an information processing device according to one embodiment of the present invention. An information processing device 1 according to one embodiment of the present invention (hereinafter referred to as information processing device 1) includes a control unit 2, a storage unit 3, a communication unit 4, a display unit 5, and an operation unit 6.

制御部2は、例えばマイクロコンピュータである。制御部2は、情報処理装置1の全体を統括的に制御する。制御部2は、不図示のCPU(Central Processing Unit)、RAM(Random Access Memory)、及びROM(Read Only Memory)を含む。 The control unit 2 is, for example, a microcomputer. The control unit 2 performs overall control of the information processing device 1. The control unit 2 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and a ROM (Read Only Memory), all of which are not shown.

記憶部3は、例えばフラッシュメモリ、ハードディスクドライブ等である。記憶部3は、再構成画像を生成するための画像再構成処理プログラム、X線撮影画像に写っている特定の物質を抽出する画像処理プログラム等の各種プログラム、各種プログラムを実行する際に用いられる各種パラメータの設定値や画像データ等の各種データを記憶する。 The storage unit 3 is, for example, a flash memory, a hard disk drive, etc. The storage unit 3 stores various programs such as an image reconstruction processing program for generating a reconstructed image, an image processing program for extracting a specific substance shown in an X-ray image, and various data such as setting values of various parameters and image data used when executing the various programs.

通信部4は、外部装置との通信を行うための通信インターフェースである。通信部4と外部装置との通信方法は、有線通信でもよく、無線通信でもよく、有線と無線とを組み合わせた通信であってもよい。外部装置としては、例えばX線撮影画像を撮影するX線撮影装置、X線撮影画像を記憶している記憶装置等を挙げることができる。 The communication unit 4 is a communication interface for communicating with an external device. The communication method between the communication unit 4 and the external device may be wired communication, wireless communication, or a combination of wired and wireless communication. Examples of the external device include an X-ray imaging device that captures X-ray images, a storage device that stores X-ray images, etc.

表示部5は、例えば液晶表示装置、有機EL(Electro Luminescence)表示装置等である。表示部5は、制御部2の制御に基づいて各種の画像を表示する。 The display unit 5 is, for example, a liquid crystal display device, an organic EL (Electro Luminescence) display device, etc. The display unit 5 displays various images based on the control of the control unit 2.

操作部6は、例えばキーボード、ポインティングデバイス等である。操作部6は、ユーザの操作内容に応じた信号を制御部2に出力する。 The operation unit 6 is, for example, a keyboard, a pointing device, etc. The operation unit 6 outputs a signal according to the user's operation to the control unit 2.

図2は、情報処理装置1の機能の一例を示す機能ブロック図である。情報処理装置1は、画像処理プログラムを実行することによって画像処理装置10として機能する。画像処理装置10は、設定部11と、算出部12と、判定部13と、を備える。 FIG. 2 is a functional block diagram showing an example of the functions of the information processing device 1. The information processing device 1 functions as an image processing device 10 by executing an image processing program. The image processing device 10 includes a setting unit 11, a calculation unit 12, and a determination unit 13.

設定部11は、X線撮影画像において着目ピクセルを囲む所定の領域を設定する。例えば、設定部11は、図3に示すようにX線撮影画像の一例である透過X線画像20において着目ピクセル21を囲む水平方向mピクセル×垂直方向nピクセルである所定の領域22を設定する。m、nはそれぞれ複数の自然数である。mとnとは、同じ値であってもよく、互いに異なる値であってもよい。ここでは、所定の領域22は着目ピクセル21を含む領域であるが、所定の領域22は着目ピクセル21を含まない(m×n-1)ピクセルの領域であってもよい。また、設定部11は、X線撮影画像を分割する単位領域も設定する。単位領域は、水平方向aピクセルであり、垂直方向bピクセルである。a、bはそれぞれ自然数である。aとbとは、同じ値であってもよく、互いに異なる値であってもよい。そして、aはmより小さく、bはnより小さい。また、mはaの倍数であり、nはbの倍数である。m、n、a、bは、例えばm=240、n=420、a=12、b=12に設定することができる。なお、a=1、b=1である場合、単位領域は、着目ピクセル21そのものとなる。 The setting unit 11 sets a predetermined region surrounding a pixel of interest in an X-ray image. For example, as shown in FIG. 3, the setting unit 11 sets a predetermined region 22 of m pixels in the horizontal direction by n pixels in the vertical direction surrounding a pixel of interest 21 in a transmitted X-ray image 20, which is an example of an X-ray image. m and n are each a plurality of natural numbers. m and n may be the same value or different values. Here, the predetermined region 22 is a region that includes the pixel of interest 21, but the predetermined region 22 may be a region of (m×n-1) pixels that does not include the pixel of interest 21. The setting unit 11 also sets a unit region into which the X-ray image is divided. The unit region is a pixels in the horizontal direction and b pixels in the vertical direction. a and b are each natural numbers. a and b may be the same value or different values. a is smaller than m, and b is smaller than n. m is a multiple of a, and n is a multiple of b. For example, m, n, a, and b can be set to m = 240, n = 420, a = 12, and b = 12. Note that when a = 1 and b = 1, the unit region is the pixel of interest 21 itself.

図3に示すように、着目ピクセル21が所定の領域22の水平方向の略中心及び所定の領域22の垂直方向の中央部分に位置するように、設定部11は所定の領域22を設定することが望ましい。ただし、着目ピクセル21が透過X線画像20の外縁に近い場所に位置する場合には、例えば図4に示すように、着目ピクセル21の位置が所定の領域22の水平方向の中央部分からの偏り及び所定の領域22の垂直方向の中央部分からの偏りができるだけ小さくなるように、設定部11は所定の領域22を設定してもよい。 As shown in FIG. 3, it is desirable for the setting unit 11 to set the specified region 22 so that the pixel of interest 21 is located approximately at the horizontal center of the specified region 22 and at the vertical center of the specified region 22. However, when the pixel of interest 21 is located near the outer edge of the transmitted X-ray image 20, the setting unit 11 may set the specified region 22 so that the deviation of the position of the pixel of interest 21 from the horizontal center of the specified region 22 and the deviation from the vertical center of the specified region 22 are as small as possible, for example, as shown in FIG. 4.

算出部12は、X線撮影画像を分割した複数の単位領域それぞれの平均輝度を算出する。算出部12は、所定の領域22の水平方向mピクセル×垂直方向nピクセルの平均輝度値L1(図5参照)を求める。具体的には、算出部12は、着目ピクセル21が属している単位領域を中心として、着目ピクセル21が属している単位領域を含む所定の領域22の水平方向mピクセル×垂直方向nピクセルの平均輝度値L1(図5参照)を求める。平均輝度値L1を求める計算及び以下の計算では、ピクセルの輝度値の代わりに単位領域の平均輝度値を使用する。 The calculation unit 12 calculates the average luminance of each of a plurality of unit areas obtained by dividing the X-ray image. The calculation unit 12 calculates an average luminance value L1 (see FIG. 5) of m pixels in the horizontal direction by n pixels in the vertical direction of a predetermined area 22. Specifically, the calculation unit 12 calculates an average luminance value L1 (see FIG. 5) of m pixels in the horizontal direction by n pixels in the vertical direction of a predetermined area 22 including the unit area to which the target pixel 21 belongs, with the unit area to which the target pixel 21 belongs being the center. In the calculation to calculate the average luminance value L1 and the following calculations, the average luminance value of the unit area is used instead of the luminance value of the pixel.

次に、算出部12は、平均輝度値L1より高い平均輝度値である単位領域を暫定的に非金属が写っている第1領域とする。 Next, the calculation unit 12 provisionally determines the unit area having an average brightness value higher than the average brightness value L1 as the first area in which a non-metallic object is depicted.

そして、算出部12は、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2(図5参照)を求める。 Then, the calculation unit 12 calculates the average brightness value L2 (see FIG. 5) of the entire unit area, which is provisionally determined to be the first area in which non-metallic material is captured.

判定部13は、暫定的に非金属が写っている第1領域とした単位領域全体の代表的輝度値に対する着目ピクセル21の輝度値の比と閾値との大小関係に基づき、着目ピクセル21に非金属が写っているか非金属以外の物質が写っているかを判定する。 The determination unit 13 determines whether the pixel of interest 21 contains a non-metal or a material other than a non-metal based on the magnitude relationship between the ratio of the luminance value of the pixel of interest 21 to the representative luminance value of the entire unit area, which is provisionally determined to be the first area containing a non-metal, and a threshold value.

図6は、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2と着目ピクセル21の輝度値との関係の一例を示す図である。図6のように、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2を横軸にとり、着目ピクセル21の輝度値を縦軸にとると、特性線T1によって金属が写っているピクセル(図6中の丸印で示されるピクセル)と非金属が写っているピクセル(図6中の菱形印で示されるピクセル)とは、特性線T1によって分離することができる。 Figure 6 is a diagram showing an example of the relationship between the average brightness value L2 of the entire unit region provisionally designated as the first region in which a non-metal is depicted and the brightness value of the pixel of interest 21. As shown in Figure 6, when the average brightness value L2 of the entire unit region provisionally designated as the first region in which a non-metal is depicted is plotted on the horizontal axis and the brightness value of the pixel of interest 21 is plotted on the vertical axis, the pixels in which a metal is depicted (pixels indicated by circles in Figure 6) and the pixels in which a non-metal is depicted (pixels indicated by diamonds in Figure 6) can be separated by the characteristic line T1.

したがって、例えば、判定部13は、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2に対する着目ピクセル21の輝度値の比が閾値よりも小さい場合に着目ピクセル21に金属が写っていると判定する。閾値は、特性線T1に対応する値に設定される。特性線T1が略直線である場合、閾値を固定値にすることができる。一方、特性線T1が曲線である場合、判定部13が、閾値を平均輝度値L2に応じて可変することが望ましい。閾値を平均輝度値L2に応じて可変することによって、判定精度が向上する。 Therefore, for example, the determination unit 13 determines that metal is included in the pixel of interest 21 when the ratio of the luminance value of the pixel of interest 21 to the average luminance value L2 of the entire unit area, which is provisionally determined to be a first area in which a non-metal is included, is smaller than a threshold value. The threshold value is set to a value corresponding to the characteristic line T1. When the characteristic line T1 is an approximately straight line, the threshold value can be set to a fixed value. On the other hand, when the characteristic line T1 is a curve, it is desirable for the determination unit 13 to vary the threshold value according to the average luminance value L2. By varying the threshold value according to the average luminance value L2, the determination accuracy is improved.

画像処理装置10は、X線撮影画像内の全てのピクセルを1つずつ順次「着目ピクセル」として、上述した処理を繰り返す。画像処理装置10は、判定結果を反映した画像、例えば、金属が写っているピクセルと非金属が写っているピクセルとが異なる輝度値である図7に示す二値化画像を生成する。図7に示す二値化画像では、金属が写っているピクセル(図7中の黒色部分)の輝度値は、非金属が写っているピクセル(図7中の白色部分)の輝度値より小さい。 The image processing device 10 repeats the above-mentioned process by sequentially treating all pixels in the X-ray image as "pixels of interest" one by one. The image processing device 10 generates an image that reflects the judgment result, for example, the binary image shown in FIG. 7 in which pixels depicting metals and pixels depicting non-metals have different brightness values. In the binary image shown in FIG. 7, the brightness value of pixels depicting metals (black parts in FIG. 7) is smaller than the brightness value of pixels depicting non-metals (white parts in FIG. 7).

例えば歯科用X線撮影装置のCT撮影では、被写体に対するX線の照射角度を変えながら、複数のX線撮影画像が撮影される。被写体に対するX線の照射角度が変わった場合、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2と着目ピクセル21の輝度値との関係も変わる。したがって、判定部13が、閾値をX線撮影画像が撮影されるときの被写体に対するX線の照射角度に応じて可変することが望ましい。閾値をX線撮影画像が撮影されるときの被写体に対するX線の照射角度に応じて可変することによって、判定精度が向上する。 For example, in CT imaging using a dental X-ray imaging device, multiple X-ray images are taken while changing the X-ray irradiation angle of the subject. When the X-ray irradiation angle of the subject changes, the relationship between the average luminance value L2 of the entire unit area, which is provisionally set as the first area in which non-metal is captured, and the luminance value of the pixel of interest 21 also changes. Therefore, it is desirable for the determination unit 13 to vary the threshold value depending on the X-ray irradiation angle of the subject when the X-ray image is taken. By varying the threshold value depending on the X-ray irradiation angle of the subject when the X-ray image is taken, the determination accuracy is improved.

画像処理装置10は、被写体に対するX線の照射角度が異なる複数の上記二値化画像に対して再構成処理を行うことで3D画像を生成し、当該3D画像から図8に示すアキシャル画像を生成する。次に、画像処理装置10は、図8に示すアキシャル画像に対して閾値処理を行い、図9に示す二値化アキシャル画像を生成する。その後、画像処理装置10は、図9に示す二値化アキシャル画像に対して順投影処理を行い、図10に示す順投影画像を生成する。図10に示す順投影画像では、金属が写っているピクセルが白色になっており、金属が抽出されている。 The image processing device 10 generates a 3D image by performing reconstruction processing on the above-mentioned multiple binary images with different X-ray irradiation angles to the subject, and generates the axial image shown in FIG. 8 from the 3D image. Next, the image processing device 10 performs threshold processing on the axial image shown in FIG. 8 to generate the binary axial image shown in FIG. 9. After that, the image processing device 10 performs forward projection processing on the binary axial image shown in FIG. 9 to generate the forward projection image shown in FIG. 10. In the forward projection image shown in FIG. 10, pixels that show metal are white, and the metal has been extracted.

画像処理装置10は、特許文献1で提案されている画像処理装置と比較して、X線撮影画像に写っている特定の物質(例えば金属)を抽出するときに、抽出精度の向上及び処理時間の短縮を見込むことができる。 Compared to the image processing device proposed in Patent Document 1, the image processing device 10 is expected to improve extraction accuracy and reduce processing time when extracting a specific substance (e.g., metal) that appears in an X-ray image.

発明の構成は、上記実施形態のほか、発明の主旨を逸脱しない範囲で種々の変更を加えることが可能である。上記実施形態は、全ての点で例示であって、制限的なものではないと考えられるべきであり、本発明の技術的範囲は、上記実施形態の説明ではなく、特許請求の範囲によって示されるものであり、特許請求の範囲と均等の意味及び範囲内に属する全ての変更が含まれると理解されるべきである。 In addition to the above-described embodiment, the configuration of the invention can be modified in various ways without departing from the spirit of the invention. The above-described embodiment should be considered to be illustrative in all respects and not restrictive, and the technical scope of the present invention is indicated by the claims, not by the description of the above-described embodiment, and should be understood to include all modifications that are equivalent in meaning and scope to the claims.

例えば、上述した実施形態では、図10に示す順投影画像において金属が抽出されているが、図7に示す二値化画像金属が抽出されてもよい。 For example, in the embodiment described above, metal is extracted in the forward projection image shown in FIG. 10, but metal may be extracted from the binarized image shown in FIG. 7.

例えば、上述した実施形態では、X線減弱が大きい物質が抽出されているが、逆にX線減弱が小さい物質が抽出されてもよい。なお、上述した実施形態では、X線減弱が大きい物質の一例として、金属が抽出されているが、X線減弱が大きい物質は、金属以外の物質であってもよい。 For example, in the above-described embodiment, a material with high X-ray attenuation is extracted, but conversely, a material with low X-ray attenuation may be extracted. Note that in the above-described embodiment, metal is extracted as an example of a material with high X-ray attenuation, but the material with high X-ray attenuation may be a material other than metal.

例えば、上述した実施形態では、透過X線画像(測定画像)が処理されているが、透過X線画像(測定画像)に対して白黒反転している投影画像が処理されてもよい。被写体を設置して撮影したときの透過X線画像(測定画像)におけるピクセルの輝度値をBv、被写体を設置しないで撮影したときの透過X線画像におけるピクセルの輝度値をWとしたとき、下記式(1)の対数変換で得られる値xが投影画像におけるピクセルの輝度値である。この対数変換をここでは白黒反転と定義することにする。

Figure 0007526503000001
For example, in the above-mentioned embodiment, a transmitted X-ray image (measurement image) is processed, but a projected image in which black and white are inverted with respect to the transmitted X-ray image (measurement image) may be processed. If the luminance value of a pixel in a transmitted X-ray image (measurement image) captured with a subject installed is Bv, and the luminance value of a pixel in a transmitted X-ray image captured without a subject installed is W, the value x obtained by logarithmic transformation of the following formula (1) is the luminance value of a pixel in the projected image. This logarithmic transformation is defined here as black and white inversion.
Figure 0007526503000001

例えば、上述した実施形態では、暫定的に非金属が写っている第1領域とした単位領域全体の代表的輝度値として、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2が用いられているが、これはあくまで一例である。例えば、暫定的に非金属が写っている第1領域とした単位領域全体の代表的輝度値として、暫定的に非金属が写っている第1領域とした単位領域全体の平均輝度値L2に所定の定数を乗じて得られる値が用いられてもよい。 For example, in the above-described embodiment, the average brightness value L2 of the entire unit area provisionally designated as the first area in which a non-metal is temporarily depicted is used as the representative brightness value of the entire unit area provisionally designated as the first area in which a non-metal is temporarily depicted, but this is merely one example. For example, a value obtained by multiplying the average brightness value L2 of the entire unit area provisionally designated as the first area in which a non-metal is temporarily depicted by a predetermined constant may be used as the representative brightness value of the entire unit area provisionally designated as the first area in which a non-metal is temporarily depicted.

閾値が複数設けられ、複数の閾値それぞれを用いた画像処理が並行して実行されてもよい。例えば、AuAgPa合金が写っているか、AuAgPa合金以外の物質が写っているかを判定するための第1閾値と、AuAgPa合金又はTiが写っているか、AuAgPa合金及びTi以外の物質が写っているかを判定するための第2閾値とが設けられてもよい。 Multiple thresholds may be provided, and image processing using each of the multiple thresholds may be performed in parallel. For example, a first threshold may be provided for determining whether an AuAgPa alloy or a material other than an AuAgPa alloy is shown, and a second threshold may be provided for determining whether an AuAgPa alloy or Ti is shown, or a material other than an AuAgPa alloy and Ti is shown.

例えば、上述した実施形態では、歯科用X線撮影装置で撮影されたX線撮影画像が用いられているが、これはあくまで一例である。例えば、非破壊検査装置で撮影されたX線撮影画像に対して本発明を適用してもよい。 For example, in the above-described embodiment, X-ray images taken by a dental X-ray device are used, but this is merely one example. For example, the present invention may also be applied to X-ray images taken by a non-destructive testing device.

1 本発明の一実施形態に係る情報処理装置
2 制御部
3 記憶部
4 通信部
5 表示部
6 操作部
10 画像処理装置
11 設定部
12 算出部
13 判定部
20 透過X線画像
21 着目ピクセル
22 所定の領域
REFERENCE SIGNS LIST 1 Information processing device according to one embodiment of the present invention 2 Control unit 3 Storage unit 4 Communication unit 5 Display unit 6 Operation unit 10 Image processing device 11 Setting unit 12 Calculation unit 13 Determination unit 20 Transmission X-ray image 21 Pixel of interest 22 Predetermined area

Claims (4)

X線撮影画像において着目ピクセルを囲む所定の領域を設定する設定部と、
前記所定の領域を構成する複数の単位領域それぞれの平均輝度値に基づき、暫定的に第1の物質が写っている第1領域とした単位領域全体の代表的輝度値を算出する算出部と、
前記代表的輝度値に対する前記着目ピクセルの輝度値の比と閾値との大小関係に基づき、前記着目ピクセルに前記第1の物質が写っているか前記第1の物質以外の物質が写っているかを判定する判定部と、
を備える、画像処理装置。
A setting unit that sets a predetermined region surrounding a pixel of interest in an X-ray image;
a calculation unit that calculates a representative luminance value of an entire unit region that is provisionally determined as a first region in which a first material is captured, based on an average luminance value of each of a plurality of unit regions that constitute the predetermined region;
a determination unit that determines whether the first material or a material other than the first material is included in the pixel of interest based on a magnitude relationship between a ratio of the luminance value of the pixel of interest to the representative luminance value and a threshold value;
An image processing device comprising:
前記閾値は、前記代表的輝度値に応じて可変する、請求項1に記載の画像処理装置。 The image processing device according to claim 1, wherein the threshold value varies depending on the representative luminance value. 前記閾値は、前記X線撮影画像が撮影されるときの被写体に対するX線の照射角度に応じて可変する、請求項1又は請求項2に記載の画像処理装置。 The image processing device according to claim 1 or 2, wherein the threshold value varies depending on the angle of X-ray irradiation with respect to the subject when the X-ray image is captured. X線撮影画像において着目ピクセルを囲む所定の領域を設定する設定ステップと、
前記所定の領域を構成する複数の単位領域それぞれの平均輝度値に基づき、暫定的に第1の物質が写っている第1領域とした単位領域全体の代表的輝度値を算出する算出ステップと、
前記代表的輝度値に対する前記着目ピクセルの輝度値の比と閾値との大小関係に基づき、前記着目ピクセルに前記第1の物質が写っているか前記第1の物質以外の物質が写っているかを判定する判定ステップと、
を備える、画像処理方法。
A setting step of setting a predetermined region surrounding a pixel of interest in the X-ray image;
a calculation step of calculating a representative luminance value of the entire unit region provisionally determined as a first region in which a first material is captured, based on an average luminance value of each of a plurality of unit regions constituting the predetermined region;
a determining step of determining whether the first material or a material other than the first material is depicted in the pixel of interest based on a magnitude relationship between a ratio of the luminance value of the pixel of interest to the representative luminance value and a threshold value;
An image processing method comprising:
JP2022172459A 2022-10-27 2022-10-27 Image processing device and image processing method Active JP7526503B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022172459A JP7526503B2 (en) 2022-10-27 2022-10-27 Image processing device and image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022172459A JP7526503B2 (en) 2022-10-27 2022-10-27 Image processing device and image processing method

Publications (2)

Publication Number Publication Date
JP2024064108A JP2024064108A (en) 2024-05-14
JP7526503B2 true JP7526503B2 (en) 2024-08-01

Family

ID=91034454

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022172459A Active JP7526503B2 (en) 2022-10-27 2022-10-27 Image processing device and image processing method

Country Status (1)

Country Link
JP (1) JP7526503B2 (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002153454A (en) 2000-11-20 2002-05-28 Hitachi Medical Corp X-ray ct device
JP2011172926A (en) 2010-02-23 2011-09-08 Carestream Health Inc Method and system for cone beam computed tomography high density object artifact reduction
US20120002855A1 (en) 2010-06-30 2012-01-05 Fujifilm Corporation Stent localization in 3d cardiac images
CN105118030A (en) 2015-08-11 2015-12-02 上海联影医疗科技有限公司 Medical image metal artifact correction method and device
CN111223156A (en) 2019-11-06 2020-06-02 深圳市深图医学影像设备有限公司 A metal artifact removal method for dental cone beam CT system
JP2021173716A (en) 2020-04-30 2021-11-01 朝日レントゲン工業株式会社 Inspection device and inspection method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6375575B2 (en) * 2014-10-22 2018-08-22 朝日レントゲン工業株式会社 Image processing apparatus, image processing method, and X-ray imaging apparatus

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002153454A (en) 2000-11-20 2002-05-28 Hitachi Medical Corp X-ray ct device
JP2011172926A (en) 2010-02-23 2011-09-08 Carestream Health Inc Method and system for cone beam computed tomography high density object artifact reduction
US20120002855A1 (en) 2010-06-30 2012-01-05 Fujifilm Corporation Stent localization in 3d cardiac images
CN105118030A (en) 2015-08-11 2015-12-02 上海联影医疗科技有限公司 Medical image metal artifact correction method and device
CN111223156A (en) 2019-11-06 2020-06-02 深圳市深图医学影像设备有限公司 A metal artifact removal method for dental cone beam CT system
JP2021173716A (en) 2020-04-30 2021-11-01 朝日レントゲン工業株式会社 Inspection device and inspection method

Also Published As

Publication number Publication date
JP2024064108A (en) 2024-05-14

Similar Documents

Publication Publication Date Title
JP6416582B2 (en) Method and apparatus for metal artifact removal in medical images
JP5828649B2 (en) Image processing apparatus, image processing method, and computer program
JP5224898B2 (en) Radiation imaging apparatus and driving method thereof
JP5722414B1 (en) Osteoporosis diagnosis support device
JP5416377B2 (en) Image processing apparatus, X-ray foreign object detection apparatus including the same, and image processing method
CN110176010B (en) Image detection method, device, equipment and storage medium
JP2011125462A (en) X-ray image processing apparatus, x-ray image processing method, and computer program
JP2011255033A (en) Radiation image region-extracting apparatus, radiation image region-extracting program, and radiographic apparatus
JP5534411B2 (en) Image processing device
JP6492553B2 (en) Image processing apparatus and program
WO2005110232A1 (en) Image processing device and method thereof
WO2013146283A1 (en) Image processing device and image processing method
US9619893B2 (en) Body motion detection device and method
JP2021036969A (en) Machine learning equipment, machine learning methods and programs
KR101804699B1 (en) Method and apparatus for correcting artifact of computed tomography image
US11151760B2 (en) Image reconstruction processing method, image reconstruction processing program, and tomography device equipped with same
JP4416823B2 (en) Image processing apparatus, image processing method, and computer program
CN115439423B (en) CT image-based identification method, device, equipment and storage medium
JP7526503B2 (en) Image processing device and image processing method
US20130308841A1 (en) Method and apparatus for image processing
JP2023553918A5 (en)
JP5784408B2 (en) Medical image processing apparatus and medical image processing method
JP4675133B2 (en) Density analysis method from X-ray CT data, computer program for executing density analysis method, and density analysis system
JP2017070590A (en) Image processing device, method for controlling the same, and computer program
JP6120913B2 (en) Image processing apparatus, image processing method, and computer program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231109

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: 20240625

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240712

R150 Certificate of patent or registration of utility model

Ref document number: 7526503

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150