JP2012110399A - Dynamic diagnosis support information generation system - Google Patents
Dynamic diagnosis support information generation system Download PDFInfo
- Publication number
- JP2012110399A JP2012110399A JP2010259650A JP2010259650A JP2012110399A JP 2012110399 A JP2012110399 A JP 2012110399A JP 2010259650 A JP2010259650 A JP 2010259650A JP 2010259650 A JP2010259650 A JP 2010259650A JP 2012110399 A JP2012110399 A JP 2012110399A
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- image
- frame
- imaging
- dynamic
- 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.)
- Granted
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 59
- 238000004458 analytical method Methods 0.000 claims abstract description 280
- 238000012545 processing Methods 0.000 claims abstract description 145
- 238000003384 imaging method Methods 0.000 claims description 165
- 238000000034 method Methods 0.000 claims description 145
- 230000008569 process Effects 0.000 claims description 135
- 230000005855 radiation Effects 0.000 claims description 94
- 238000001514 detection method Methods 0.000 claims description 37
- 230000033001 locomotion Effects 0.000 claims description 27
- 230000000241 respiratory effect Effects 0.000 claims description 7
- 238000004891 communication Methods 0.000 description 50
- 238000003860 storage Methods 0.000 description 46
- 210000004072 lung Anatomy 0.000 description 42
- 230000008859 change Effects 0.000 description 34
- 238000009423 ventilation Methods 0.000 description 26
- 238000012937 correction Methods 0.000 description 23
- 238000004364 calculation method Methods 0.000 description 22
- 230000017531 blood circulation Effects 0.000 description 19
- 238000007726 management method Methods 0.000 description 16
- 230000029058 respiratory gaseous exchange Effects 0.000 description 14
- 210000000038 chest Anatomy 0.000 description 13
- 238000011156 evaluation Methods 0.000 description 12
- 238000012546 transfer Methods 0.000 description 9
- 230000007423 decrease Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000000284 resting effect Effects 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000010191 image analysis Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 230000003434 inspiratory effect Effects 0.000 description 5
- 230000010349 pulsation Effects 0.000 description 5
- 239000004065 semiconductor Substances 0.000 description 5
- 238000004904 shortening Methods 0.000 description 5
- 230000002861 ventricular Effects 0.000 description 5
- 230000000747 cardiac effect Effects 0.000 description 4
- 230000036387 respiratory rate Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 206010047571 Visual impairment Diseases 0.000 description 3
- 238000012790 confirmation Methods 0.000 description 3
- 230000002950 deficient Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000010408 film Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000002438 flame photometric detection Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 210000000115 thoracic cavity Anatomy 0.000 description 2
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 1
- 208000001871 Tachycardia Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- OJIJEKBXJYRIBZ-UHFFFAOYSA-N cadmium nickel Chemical compound [Ni].[Cd] OJIJEKBXJYRIBZ-UHFFFAOYSA-N 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000011976 chest X-ray Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052987 metal hydride Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- PXHVJJICTQNCMI-UHFFFAOYSA-N nickel Substances [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 1
- -1 nickel metal hydride Chemical class 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 125000003003 spiro group Chemical group 0.000 description 1
- 230000006794 tachycardia Effects 0.000 description 1
- 208000008203 tachypnea Diseases 0.000 description 1
- 206010043089 tachypnoea Diseases 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- YDLQKLWVKKFPII-UHFFFAOYSA-N timiperone Chemical compound C1=CC(F)=CC=C1C(=O)CCCN1CCC(N2C(NC3=CC=CC=C32)=S)CC1 YDLQKLWVKKFPII-UHFFFAOYSA-N 0.000 description 1
- 229950000809 timiperone Drugs 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
本発明は、動態診断支援情報生成システムに関する。 The present invention relates to a dynamic diagnosis support information generation system.
近年、FPD(flat panel detector)等の半導体イメージセンサを利用して胸部等の動態画像を撮影し、診断に応用する試みがなされるようになってきている。例えば、特許文献1には、胸部の動態画像を撮影する放射線画像取得装置が記載されている。この特許文献1には、動態画像を医師が観察した際に呼吸運動の安定した追跡認識を可能とするために、動態画像の撮影速度は1秒あたり3フレーム以上とすることが記載されている。 In recent years, attempts have been made to take a dynamic image of the chest using a semiconductor image sensor such as a flat panel detector (FPD) and apply it to diagnosis. For example, Patent Document 1 describes a radiological image acquisition apparatus that captures a dynamic image of a chest. This Patent Document 1 describes that the imaging speed of a dynamic image is 3 frames or more per second in order to enable stable tracking recognition of respiratory motion when a dynamic image is observed by a doctor. .
また、動態撮影で得られた一連のフレーム画像に基づいて当該動態に関する特徴量を算出し、早期診断に向けて診断支援情報として医師に提供することが提案されている。例えば、特許文献2には、動態画像を構成する複数のフレーム画像間で信号値の差を示す差分画像を生成し、生成された差分画像から各信号値の最大値を求めて表示する技術が記載されている。 In addition, it has been proposed to calculate a feature amount related to the dynamics based on a series of frame images obtained by dynamic imaging and to provide the doctor with diagnosis support information for early diagnosis. For example, Patent Document 2 discloses a technique for generating a difference image indicating a difference in signal values between a plurality of frame images constituting a dynamic image, and obtaining and displaying the maximum value of each signal value from the generated difference image. Are listed.
しかしながら、特許文献1及び特許文献2のいずれにおいても、動態画像を解析して特徴量を算出する際の最適な撮影速度(フレームレート)については記載されていない。
一般的には、フレームレートを増加させるほど精度の良い解析結果が得られるように思われる。しかし、例えば、動態撮影の各撮影における照射線量を一定とした場合、1秒当たりの撮影枚数であるフレームレートが大きくなるほど患者の被曝線量が増加し、好ましくない。一方、動態撮影におけるトータルの照射線量を一定とすると、フレームレートが大きくなるほど個々のフレーム画像におけるノイズが増え、解析結果に影響する。また、フレームレートが増大するほどデータ量が増え、解析に時間がかかってしまう。
そこで、本願発明者らは、鋭意検討の結果、解析用の動態画像を撮影する際に、診断に使用し得る解析結果が出力でき、データ転送や解析演算の処理時間及び患者の被曝を最低限に抑え、ハードウエアの構成をいたずらに増大させることのないフレームレートを見出した。
However, neither Patent Document 1 nor Patent Document 2 describes an optimum shooting speed (frame rate) when a dynamic image is analyzed to calculate a feature amount.
In general, it seems that more accurate analysis results are obtained as the frame rate is increased. However, for example, if the irradiation dose in each dynamic imaging is constant, the patient exposure dose increases as the frame rate, which is the number of images taken per second, increases. On the other hand, if the total irradiation dose in dynamic imaging is constant, the noise in each frame image increases as the frame rate increases, which affects the analysis result. Further, as the frame rate increases, the amount of data increases and analysis takes time.
Therefore, as a result of intensive studies, the inventors of the present application can output an analysis result that can be used for diagnosis when capturing a dynamic image for analysis, minimizing the processing time of data transfer and analysis operation and patient exposure. And a frame rate that does not unnecessarily increase the hardware configuration.
本発明の課題は、診断に使用し得る解析結果が出力でき、処理時間及び患者の被曝を最低限に抑え、ハードウエアの構成をいたずらに増大させることのない動態診断支援情報生成システムを提供することである。 An object of the present invention is to provide a dynamic diagnosis support information generation system that can output analysis results that can be used for diagnosis, minimizes processing time and patient exposure, and does not unnecessarily increase the hardware configuration. That is.
上記課題を解決するため、請求項1に記載の発明は、
パルス照射可能な放射線源と、
2次元状に配置された複数の検出素子を有し、前記複数の検出素子のそれぞれにおいて前記放射線源によりパルス照射された放射線を検出して順次フレーム画像を生成する放射線検出器と、
前記放射線源及び前記放射線検出器により前記被写体の動態を撮影することにより得られた複数のフレーム画像に基づいて前記被写体の動態の特徴量を算出する解析手段と、
を備える動態診断支援情報生成システムであって、
前記放射線検出器におけるフレームレートが3.75枚/秒以上である。
In order to solve the above-mentioned problem, the invention described in claim 1
A radiation source capable of pulse irradiation;
A radiation detector that has a plurality of detection elements arranged in a two-dimensional shape, detects radiation irradiated by the radiation source at each of the plurality of detection elements, and sequentially generates a frame image;
Analyzing means for calculating feature quantities of the dynamics of the subject based on a plurality of frame images obtained by photographing the dynamics of the subject by the radiation source and the radiation detector;
A dynamic diagnosis support information generation system comprising:
The frame rate in the radiation detector is 3.75 sheets / second or more.
請求項2に記載の発明は、請求項1に記載の発明において、
前記複数のフレーム画像の画素数を低減するための間引き処理を施す間引き処理手段を備え、
前記解析手段は、前記間引き処理手段により間引き処理された前記複数のフレーム画像に基づいて特徴量の算出を行う。
The invention according to claim 2 is the invention according to claim 1,
Comprising a thinning processing means for performing a thinning process for reducing the number of pixels of the plurality of frame images;
The analysis unit calculates a feature amount based on the plurality of frame images thinned by the thinning processing unit.
請求項3に記載の発明は、請求項2に記載の発明において、
前記間引き処理手段は、前記放射線検出器に設けられている。
The invention according to claim 3 is the invention according to claim 2,
The thinning processing means is provided in the radiation detector.
請求項4に記載の発明は、請求項1〜3の何れか一項に記載の発明において、
前記解析手段は、前記被写体の動態1サイクル以上のフレーム画像群に基づいて前記被写体の動態の特徴量を算出する。
The invention according to claim 4 is the invention according to any one of claims 1 to 3,
The analysis means calculates a feature quantity of the subject's dynamics based on a frame image group of one cycle or more of the subject's dynamics.
請求項5に記載の発明は、請求項1〜4の何れか一項に記載の発明において、
前記動態診断支援情報生成システムは、前記放射線源と前記放射線検出器により胸部の動態撮影を行い、
前記解析手段は、前記被写体の呼吸運動に関する特徴量を算出する。
The invention according to claim 5 is the invention according to any one of claims 1 to 4,
The dynamic diagnosis support information generation system performs dynamic imaging of the chest with the radiation source and the radiation detector,
The analysis unit calculates a feature amount related to the respiratory motion of the subject.
請求項6に記載の発明は、請求項5に記載の発明において、
前記解析手段は、呼気期間、吸気期間のそれぞれにおける気流速度の指標となる特徴量を算出する。
The invention according to claim 6 is the invention according to claim 5,
The analysis means calculates a feature quantity that serves as an index of airflow velocity in each of the expiration period and the inspiration period.
本発明によれば、診断に使用し得る解析結果が出力でき、処理時間及び患者の被曝を最低限に抑え、ハードウエアの構成をいたずらに増大させることのない動態診断支援情報生成システムを提供することができる。 According to the present invention, a dynamic diagnosis support information generation system that can output analysis results that can be used for diagnosis, minimizes processing time and patient exposure, and does not unnecessarily increase the hardware configuration is provided. be able to.
以下、本発明に係る診断支援情報生成システムの実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。 Hereinafter, embodiments of a diagnosis support information generation system according to the present invention will be described with reference to the drawings. However, the present invention is not limited to the following illustrated examples.
[第1の実施の形態]
図1は、本発明の第1の実施形態に係る診断支援情報生成システム100の全体構成例を示す図である。
図1に示す撮影室R1〜R3は、患者の身体の一部である被写体(すなわち患者の撮影部位)に放射線を照射して被写体の動態撮影又は静止画撮影を行うための室である。
動態撮影とは、被写体に対し、X線等の放射線をパルス的に連続照射して複数の画像を取得(即ち、連続撮影)することをいう。動態撮影では、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性(サイクル)を持つ被写体の動態を撮影する。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
静止画撮影とは、従来のフィルム方式やCR方式と同様に撮影部位の濃度分解能に基づく診断に使用されるもので、被写体に対し、X線等の放射線を1回照射して一枚の静止画像を取得する ことをいう。
[First Embodiment]
FIG. 1 is a diagram showing an example of the overall configuration of a diagnosis support information generation system 100 according to the first embodiment of the present invention.
The imaging rooms R1 to R3 shown in FIG. 1 are rooms for performing dynamic imaging or still image imaging of a subject by irradiating a subject that is a part of the patient's body (that is, a patient's imaging region) with radiation.
Dynamic imaging refers to acquiring a plurality of images (ie, continuous imaging) by continuously irradiating a subject with radiation such as X-rays in a pulsed manner. In dynamic imaging, for example, the dynamics of a subject having a periodicity (cycle), such as changes in the shape of lung expansion and contraction associated with respiratory motion, heart pulsation, and the like are captured. A series of images obtained by this continuous shooting is called a dynamic image. Each of the plurality of images constituting the dynamic image is called a frame image.
Still image photography is used for diagnosis based on the density resolution of the imaged part in the same way as in the conventional film method and CR method. A single piece of still image is irradiated by irradiating the subject with radiation such as X-rays once. This refers to acquiring an image.
撮影室R1は、単射及び連射が可能な放射線源3aが設けられた、被写体の動態撮影又は静止画撮影をするための室である。
撮影室R1には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3aと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が備えられている。
The imaging room R1 is a room for performing dynamic imaging or still image imaging of a subject provided with a radiation source 3a capable of single shot and continuous shot.
In the radiographing room R1, for example, a bucky device 1 for standing photography, a bucky device 2 for standing photography, a radiation source 3a, a cradle 4, a console 5, an operation console 6, and an access point AP are provided. , Is provided.
撮影室R2は、単射のみが可能な放射線源3b及びポータブル撮影用の放射線源3cが設けられた、被写体の静止画撮影をするための室である。
撮影室R2には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3b、3cと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が設けられている。
The imaging room R2 is a room for taking a still image of a subject, provided with a radiation source 3b capable of only a single shot and a radiation source 3c for portable imaging.
In the radiographing room R2, for example, a bucky device 1 for standing photography, a bucky device 2 for standing photography, radiation sources 3b and 3c, a cradle 4, a console 5, a console 6, and an access point are provided. AP is provided.
撮影室R3は、放射線源3bが設けられた、被写体の静止画撮影をするための室である。
撮影室R3には、例えば、立位撮影用のブッキー装置1と、臥位撮影用のブッキー装置2と、放射線源3bと、クレードル4と、コンソール5と、操作卓6と、アクセスポイントAPと、が設けられている。
The imaging room R3 is a room in which a radiation source 3b is provided for taking a still image of a subject.
In the photographing room R3, for example, a bucky device 1 for standing position photographing, a bucky device 2 for standing position photographing, a radiation source 3b, a cradle 4, a console 5, a console 6, and an access point AP are provided. , Is provided.
なお、各撮影室R1〜R3には前室Raと撮影実施室Rbが設けられ、前室Raにコンソール5及び操作卓6が備えられることで、撮影技師等の操作者の被曝を防止するようになっている。 Each of the photographing rooms R1 to R3 is provided with a front room Ra and a photographing room Rb, and the front room Ra is provided with a console 5 and a console 6 so as to prevent exposure of an operator such as a photographing engineer. It has become.
ブッキー装置1は、立位での撮影時にFPD9a又は9bを保持して撮影を行うための装置である。
図2に、ブッキー装置1の機能構成例を示す。図2に示すように、ブッキー装置1は、
制御部11と、検出器装着部12と、通信I/F13と、駆動部14とを備えて構成されている。
The bucky device 1 is a device for holding the FPD 9a or 9b when shooting in a standing position.
FIG. 2 shows a functional configuration example of the bucky device 1. As shown in FIG.
A control unit 11, a detector mounting unit 12, a communication I / F 13, and a drive unit 14 are provided.
制御部11は、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)により構成される。制御部11のROMには、ブッキー装置1の各部を制御するための各種処理プログラムや処理に必要なデータ、当該ブッキー装置1の識別情報であるブッキーID等が記憶されている。CPUは、ROMに記憶されているプログラムとの協働によりブッキー装置1各部の動作を統括的に制御する。
例えば、制御部11は、検出器装着部12にFPD9a又は9bが装着されると、装着されたFPDにコネクタ12bを介してFPDID(FPDの識別情報)の送信要求を行い、FPDIDが受信されると、自己の識別番号であるブッキーIDをFPDIDに対応付けて通信I/F13を介してコンソール5に送信する。また、受信されたFPDIDを一時的にRAMに記憶する。
また、例えば、制御部11は、検出器装着部12からFPDが抜き取られると、通信I/F13を介してコンソール5にFPDIDを送信し、当該FPDIDの(撮影管理テーブル521からの)消去要求を行う。
The control unit 11 includes a CPU (Central Processing Unit), a ROM (Read Only Memory), and a RAM (Random Access Memory). The ROM of the control unit 11 stores various processing programs for controlling each unit of the bucky device 1, data necessary for the processing, a bucky ID that is identification information of the bucky device 1, and the like. The CPU comprehensively controls the operation of each part of the bucky device 1 in cooperation with a program stored in the ROM.
For example, when the FPD 9a or 9b is attached to the detector attachment unit 12, the control unit 11 sends an FPDID (FPD identification information) transmission request to the attached FPD via the connector 12b, and the FPDID is received. Then, the Bucky ID which is its own identification number is associated with the FPDID and transmitted to the console 5 via the communication I / F 13. Also, the received FPDID is temporarily stored in the RAM.
Further, for example, when the FPD is extracted from the detector mounting unit 12, the control unit 11 transmits the FPDID to the console 5 via the communication I / F 13, and issues a request for deleting the FPDID (from the imaging management table 521). Do.
検出器装着部12は、FPD(FPD9a又はFPD9b)を保持するための保持部12aと、保持部12aに装着されたFPDのコネクタ94を接続するためのコネクタ12bとを有する。コネクタ12bは、保持部12aに装着されたFPDとの間でデータ送受信を行ったり、FPDに電力を供給したりする。 The detector mounting unit 12 includes a holding unit 12a for holding an FPD (FPD 9a or FPD 9b) and a connector 12b for connecting a connector 94 of the FPD mounted on the holding unit 12a. The connector 12b transmits / receives data to / from the FPD attached to the holding unit 12a and supplies power to the FPD.
通信I/F13は、アクセスポイントAPを介してコンソール5等の外部機器と通信ケールを介してデータ送受信を行うためのインターフェースである。 The communication I / F 13 is an interface for performing data transmission / reception via an access point AP and an external device such as the console 5 via a communication cable.
駆動部14は、図示しないフットスイッチ等の操作に応じて検出器装着部12を垂直方向又は水平方向に移動させる。 The drive unit 14 moves the detector mounting unit 12 in the vertical direction or the horizontal direction in accordance with an operation of a foot switch or the like (not shown).
ブッキー装置2は、臥位での撮影時にFPD9a又は9bを保持して撮影を行うための装置である。
ブッキー装置2は、制御部21と、検出器装着部22と、通信I/F23と、駆動部24とを備えて構成されている。制御部21、検出器装着部22、通信I/F23、駆動部24は、それぞれ上述の制御部11、検出器装着部12、通信I/F13、駆動部14と同様の構成であるので説明を援用する。更に、ブッキー装置2は、被写体を載置するための被写体台26を備える。
The bucky device 2 is a device for holding the FPD 9a or 9b during shooting in the supine position.
The bucky device 2 includes a control unit 21, a detector mounting unit 22, a communication I / F 23, and a drive unit 24. The control unit 21, the detector mounting unit 22, the communication I / F 23, and the driving unit 24 have the same configurations as the control unit 11, the detector mounting unit 12, the communication I / F 13, and the driving unit 14, respectively. Incorporate. Further, the bucky device 2 includes a subject table 26 for placing a subject.
放射線源3aは、単射及び連射(パルス照射)が可能な放射線発生装置である。放射線源3aは、例えば、撮影室R1、R3の天井から吊り下げられており、撮影時にはコンソール5からの指示に基づいて起動され、図示しない駆動機構によりにより所定の位置、向きに調整されるようになっている。そして、放射線の照射方向を変えることで、立位用のブッキー装置1又は臥位用のブッキー装置2に装着されたFPD9a又は9bに対して放射線を照射することができるようになっている。また、放射線源3aは、コンソール5からの指示に従って放射線を1回又は連続して照射し、静止画撮影又は動態撮影を行う。 The radiation source 3a is a radiation generator capable of single shot and continuous shot (pulse irradiation). The radiation source 3a is suspended from the ceilings of the imaging rooms R1 and R3, for example, and is activated based on an instruction from the console 5 at the time of imaging, and is adjusted to a predetermined position and orientation by a driving mechanism (not shown). It has become. Then, by changing the irradiation direction of the radiation, the FPD 9a or 9b mounted on the standing-side bucky device 1 or the recumbent bucky device 2 can be irradiated with the radiation. The radiation source 3a irradiates the radiation once or continuously in accordance with an instruction from the console 5 to perform still image capturing or dynamic image capturing.
放射線源3bは、単射のみが可能な放射線発生装置である。放射線源3bは、例えば、撮影室R2の天井から吊り下げられており、撮影時にはコンソール5からの指示に基づいて起動され、図示しない駆動機構によりにより所定の位置、向きに調整されるようになっている。そして、放射線の照射方向を変えることで、立位用のブッキー装置1又は臥位用のブッキー装置2に装着されたFPD9a又は9bに対して放射線を照射することができるようになっている。また、放射線源3bは、コンソール5からの指示に従って放射線を1回照射し、静止画撮影を行う。 The radiation source 3b is a radiation generating device capable of only a single shot. For example, the radiation source 3b is suspended from the ceiling of the imaging room R2, and is activated based on an instruction from the console 5 at the time of imaging, and is adjusted to a predetermined position and orientation by a drive mechanism (not shown). ing. Then, by changing the irradiation direction of the radiation, the FPD 9a or 9b mounted on the standing-side bucky device 1 or the recumbent bucky device 2 can be irradiated with the radiation. The radiation source 3b emits radiation once in accordance with an instruction from the console 5 and performs still image shooting.
放射線源3cは、移動可能なポータブル撮影用の放射線源である。放射線源3cは、単射のみが可能である。 The radiation source 3c is a movable portable radiation source. The radiation source 3c can only be a single shot.
クレードル4は、装着されたFPDと接続するための図示しないコネクタを有し、FPDが装着されると、コネクタを介して装着されたFPDからFPDIDを取得して、コンソール5に通知する。
コンソール5は、当該FPDIDを取得すると、以後、当該FPDを自己の制御下におき、起動やスリープ遷移等を制御する。
なお、本実施の形態においては、FPDを撮影室内に持ち込んだ際及び持ち出す際にクレードル4に装着することで、FPDの撮影室への侵入及び持ち出しをクレードル4を介してコンソール5で検知できるようになっている。
なお、FPDの撮影室への侵入及び持ち出しは、上記のクレードル方式以外も採用可能で、例えば、WO2008/111355号公報に開示されているRFID方式等を用いることができる。
The cradle 4 has a connector (not shown) for connecting to the mounted FPD. When the FPD is mounted, the cradle 4 acquires the FPDID from the FPD mounted via the connector and notifies the console 5 of the FPDID.
After acquiring the FPDID, the console 5 thereafter puts the FPD under its own control and controls activation, sleep transition, and the like.
In the present embodiment, the FPD can be detected by the console 5 via the cradle 4 by attaching the FPD to the cradle 4 when the FPD is brought into or taken out of the photographic room. It has become.
In addition, the FPD can be entered and taken out of the photographic room other than the above cradle method, and for example, the RFID method disclosed in WO2008 / 111355 can be used.
コンソール5は、放射線源3a、3bやFPD9a、9bを制御することにより撮影を制御するための装置である。コンソール5は、LAN(Local Area Network)を介してHIS/RIS(Hospital Information System/ Radiology Information System)7、解析用WS8、PACS(Picture Archiving and Communication System)10等に接続されており、HIS/RIS7から送信された撮影オーダ情報に基づいて、コンソール5が対応つけられている(設置されている)撮影室でそのオーダの撮影が可能であるか否かを判断し、判断結果を表示する。そして、撮影が可能である場合、コンソール5は撮影に用いられる放射線源及びFPDを起動させる等の制御をして撮影を行わせる。 The console 5 is a device for controlling imaging by controlling the radiation sources 3a and 3b and the FPDs 9a and 9b. The console 5 is connected to a HIS / RIS (Hospital Information System / Radiology Information System) 7, an analysis WS 8, a PACS (Picture Archiving and Communication System) 10, etc. via a LAN (Local Area Network). Based on the imaging order information transmitted from, it is determined whether or not imaging of the order is possible in the imaging room to which the console 5 is associated (installed), and the determination result is displayed. If imaging is possible, the console 5 controls the radiation source and FPD used for imaging to perform imaging.
図3に、コンソール5の要部構成例を示す。図3に示すように、コンソール5は、制御部51、記憶部52、入力部53、表示部54、通信I/F55、ネットワーク通信部56等を備えて構成されており、各部はバス57により接続されている。 In FIG. 3, the example of a principal part structure of the console 5 is shown. As shown in FIG. 3, the console 5 includes a control unit 51, a storage unit 52, an input unit 53, a display unit 54, a communication I / F 55, a network communication unit 56, and the like. It is connected.
制御部51は、CPU、RAM等により構成される。制御部51のCPUは、記憶部52に記憶されているシステムプログラムや処理プログラム等の各種プログラムを読み出してRAMに展開し、展開されたプログラムに従って各種処理を実行する。
例えば、制御部51は、通信I/F55を介してFPDID及びブッキーIDが受信されると、記憶部52の撮影管理テーブル521(図4参照)の受信されたブッキーIDに対応する領域にFPDIDを書き込む。また、制御部51は、通信I/F55を介してクレードル4からFPDIDが受信されると、記憶部52の撮影管理テーブル521のブッキーIDが対応付けられていない領域にFPDIDを書き込む。また、制御部51は、ブッキー装置1又は2を介して画像データが受信されると、画像受信時刻を撮影管理テーブル521の送信元のブッキー装置のブッキーIDに対応する領域に記憶させる。
The control unit 51 includes a CPU, a RAM, and the like. The CPU of the control unit 51 reads out various programs such as system programs and processing programs stored in the storage unit 52, expands them in the RAM, and executes various processes according to the expanded programs.
For example, when the FPDID and the Bucky ID are received via the communication I / F 55, the control unit 51 assigns the FPDID to an area corresponding to the received Bucky ID in the shooting management table 521 (see FIG. 4) of the storage unit 52. Write. In addition, when the FPDID is received from the cradle 4 via the communication I / F 55, the control unit 51 writes the FPDID in an area in the imaging management table 521 of the storage unit 52 that is not associated with a bucky ID. Further, when the image data is received via the bucky device 1 or 2, the control unit 51 stores the image reception time in an area corresponding to the bucky ID of the source bucky device in the photographing management table 521.
また、例えば、制御部51は、所定時間毎にネットワーク通信部56を介してHIS/RIS7に問い合わせを行い、新たにHIS/RIS7で登録された撮影オーダ情報を取得する。
また、例えば、制御部51は、後述する撮影・解析処理を実行し、HIS/RIS7から取得した撮影オーダ情報に基づいて、このコンソール5が設置されている撮影室でそのオーダの撮影が可能であるか否かを判断し、判断結果を表示する。そして、撮影が可能である場合、コンソール5は撮影に用いられる放射線源及び撮影に用いられるFPDを制御して撮影を行わせる。
Further, for example, the control unit 51 makes an inquiry to the HIS / RIS 7 via the network communication unit 56 every predetermined time, and acquires imaging order information newly registered in the HIS / RIS 7.
In addition, for example, the control unit 51 executes imaging / analysis processing described later, and based on the imaging order information acquired from the HIS / RIS 7, the order can be captured in the imaging room in which the console 5 is installed. It is judged whether there is, and the judgment result is displayed. And when imaging | photography is possible, the console 5 controls the radiation source used for imaging | photography, and FPD used for imaging | photography, and performs imaging | photography.
記憶部52は、例えばHDD(Hard Disk Drive)や半導体の不揮発性メモリ等で構成
されている。
The storage unit 52 includes, for example, an HDD (Hard Disk Drive), a semiconductor nonvolatile memory, or the like.
記憶部52には、各種のプログラム及びデータが記憶されている。
例えば、記憶部52には、各撮影室における撮影を管理するための撮影管理テーブル521が記憶されている。
図4に、撮影管理テーブル521のデータ格納例を示す。図4に示すように、撮影管理テーブル521には、「ブッキーID」、「管球タイプ」、「FPDID」、「画像受信時刻」等の項目が設けられている。「ブッキーID」、「管球タイプ」の領域には、当該コンソール5が設置されている撮影室に設けられているブッキー装置及び放射線源のタイプの情報が予め設定されている。ブッキーIDに対応付けられている「FPDID」の領域は、そのブッキーIDに装着されているFPDを管理するための領域であり、ブッキー装置からFPDID及びブッキーIDが受信された際に、受信されたFPDIDがブッキーIDに対応付けて格納される。また、ブッキーIDが対応付けられていない「FPDID」の領域は、撮影室に存在するFPDを管理するための領域であり、クレードル4からFPDIDが受信された際に、受信されたFPDIDが格納される。なお、ブッキー装置からFPDが抜き取られてFPDIDの消去要求が受信された場合には、制御部51により消去要求されたFPDIDのうち当該ブッキーIDに対応して格納されているFPDIDは消去される。また、ブッキーIDが対応付けられていない「FPDID」の領域に格納されているFPDIDがクレードル4から受信されると、制御部51によりそのFPDIDのFPDは撮影室から持ち出された(即ち、撮影室に存在しなくなった)と判断され、そのFPDIDは撮影管理テーブル521から消去される。また、「画像受信時刻」には、通信I/F55から画像データが受信された際に、その時刻が格納される。
The storage unit 52 stores various programs and data.
For example, the storage unit 52 stores a shooting management table 521 for managing shooting in each shooting room.
FIG. 4 shows a data storage example of the imaging management table 521. As shown in FIG. 4, the shooting management table 521 includes items such as “Bucky ID”, “tube type”, “FPDID”, and “image reception time”. In the “Bucky ID” and “tube type” areas, information on the type of the Bucky device and the radiation source provided in the radiographing room in which the console 5 is installed is set in advance. The “FPDID” area associated with the Bucky ID is an area for managing the FPD attached to the Bucky ID, and is received when the FPDID and the Bucky ID are received from the Bucky device. The FPDID is stored in association with the Bucky ID. The “FPDID” area not associated with the Bucky ID is an area for managing the FPD existing in the photographing room, and when the FPDID is received from the cradle 4, the received FPDID is stored. The When an FPDID is deleted from the Bucky device and an FPDID deletion request is received, the FPDID stored in correspondence with the Bucky ID among the FPDIDs requested to be deleted by the control unit 51 is deleted. When an FPDID stored in the “FPDID” area not associated with a Bucky ID is received from the cradle 4, the FPDID of the FPDID is taken out of the shooting room by the control unit 51 (that is, the shooting room). The FPDID is deleted from the imaging management table 521. The “image reception time” stores the time when image data is received from the communication I / F 55.
また、記憶部52には、画像データから患部を検出するための自動部位認識に基づく階調処理・周波数処理等の画像処理を行うためのプログラム等、各種のプログラムが記憶されているほか、撮影画像の画像データを部位毎の診断に適した画質に調整するための画像処理パラメータ(階調処理に用いる階調曲線を定義したルックアップテーブル、周波数処理の強調度等)等が記憶されている。 In addition, the storage unit 52 stores various programs such as a program for performing image processing such as gradation processing and frequency processing based on automatic part recognition for detecting an affected part from image data. Stores image processing parameters (look-up table defining tone curve used for tone processing, enhancement degree of frequency processing, etc.) for adjusting image data of image to image quality suitable for diagnosis for each part .
また、記憶部52には、撮影の種類(動態又は静止画)と撮影部位との組み合わせに対応付けて放射線照射条件及び画像読取条件が記憶されている。放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、フィルタ種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、連続撮影において、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ(マトリックスサイズ)等である。フレームレートは、1秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、連続撮影において、1回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。 The storage unit 52 stores radiation irradiation conditions and image reading conditions in association with combinations of imaging types (dynamics or still images) and imaging regions. The radiation irradiation conditions are, for example, pulse rate, pulse width, pulse interval during continuous irradiation, the number of imaging frames per imaging, X-ray tube current value, X-ray tube voltage value, filter type, and the like. The pulse rate is the number of times of radiation irradiation per second, and matches the frame rate described later. The pulse width is a radiation irradiation time per one irradiation. The pulse interval is the time from the start of one radiation irradiation to the start of the next radiation irradiation in continuous imaging, and coincides with a frame interval described later. The image reading conditions are, for example, a frame rate, a frame interval, a pixel size, an image size (matrix size), and the like. The frame rate is the number of frame images acquired per second and matches the pulse rate. The frame interval is the time from the start of one frame image acquisition operation to the start of the next frame image acquisition operation in continuous shooting, and coincides with the pulse interval.
また、記憶部52には、診断支援情報生成システム100に登録されている各FPD9a、9bのFPDIDと、そのFPDで撮影可能な撮影の種類(動態、静止画)とが対応付けて記憶されている。また、記憶部52には、診断支援情報生成システム100に登録されている各ブッキー装置1、2のブッキーIDと、そのブッキー装置で撮影可能な***(立位又は臥位)とが対応付けて記憶されている。 Further, the storage unit 52 stores the FPDID of each FPD 9a, 9b registered in the diagnosis support information generation system 100 and the type of shooting (dynamic, still image) that can be shot with the FPD in association with each other. Yes. Further, the storage unit 52 associates the bucky ID of each of the bucky devices 1 and 2 registered in the diagnosis support information generation system 100 with the body posture (standing or standing) that can be photographed by the bucky device. It is remembered.
また、記憶部52は、所定時間毎にHIS/RIS7から送信される撮影オーダ情報を記憶する。 The storage unit 52 stores imaging order information transmitted from the HIS / RIS 7 every predetermined time.
入力部53は、文字入力キー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを、入力信号として制御部51に出力する。 The input unit 53 includes a keyboard having character input keys, numeric input keys, various function keys, and the like, and a pointing device such as a mouse, and a key pressing signal pressed by the keyboard and an operation signal by the mouse. Are output to the control unit 51 as an input signal.
表示部54は、例えば、CRT(Cathode Ray Tube)やLCD(Liquid Crystal Display)等のモニタを備えて構成されており、制御部51から入力される表示信号の指
示に従って、各種画面を表示する。
なお、表示部54の画面上に、透明電極を格子状に配置した感圧式(抵抗膜圧式)のタッチパネル(図示せず)を形成し、表示部54と入力部53とが一体に構成されるタッチスクリーンとしてもよい。この場合、タッチパネルは、手指やタッチペン等で押下された力点のXY座標を電圧値で検出し、検出された位置信号が操作信号として制御部51に出力されるように構成される。なお、表示部54は、一般的なPC(Personal Computer)に用いられるモニタよりも高精細のものであってもよい。
The display unit 54 includes a monitor such as a CRT (Cathode Ray Tube) or an LCD (Liquid Crystal Display), for example, and displays various screens according to instructions of display signals input from the control unit 51.
A pressure-sensitive (resistive film pressure type) touch panel (not shown) in which transparent electrodes are arranged in a grid pattern is formed on the screen of the display unit 54, and the display unit 54 and the input unit 53 are configured integrally. It may be a touch screen. In this case, the touch panel is configured to detect the XY coordinates of the power point pressed by a finger, a touch pen, or the like as a voltage value, and output the detected position signal to the control unit 51 as an operation signal. The display unit 54 may have a higher definition than a monitor used in a general PC (Personal Computer).
通信I/F55は、ブッキー装置1、ブッキー装置2、放射線源3a〜3c、FPD9a又は9bとアクセスポイントAPを介して接続し、無線、または有線によりデータ送受信を行うためのインターフェースである。本実施形態において、通信I/F55はアクセスポイントAPを介して必要に応じてFPD9a、9bに対してポーリング信号を送信する。 The communication I / F 55 is an interface for connecting to the Bucky device 1, the Bucky device 2, the radiation sources 3a to 3c, the FPD 9a or 9b via the access point AP, and performing data transmission / reception wirelessly or by wire. In the present embodiment, the communication I / F 55 transmits a polling signal to the FPDs 9a and 9b as necessary via the access point AP.
ネットワーク通信部56は、ネットワークインターフェース等により構成され、スイッチングハブを介して通信ネットワークNに接続された外部機器との間でデータの送受信を行う。 The network communication unit 56 includes a network interface or the like, and transmits / receives data to / from an external device connected to the communication network N via a switching hub.
操作卓6は、撮影室内の放射線源に接続され、放射線照射指示を入力するための入力装置である。 The console 6 is an input device that is connected to a radiation source in the imaging room and inputs a radiation irradiation instruction.
HIS/RIS7は、問診結果等に基づくオペレータによる登録操作に応じて撮影オーダ情報を生成する。撮影オーダ情報は、例えば被写体となる患者の氏名等の患者情報や、撮影部位、撮影方向、***(立位、臥位)、撮影方法、解析要否、解析項目等の撮影予約に関する情報等を含んでいる。なお、撮影オーダ情報はここに例示したものに限定されず、これ以外の情報を含んでいてもよいし、上記に例示した情報のうちの一部でもよい。 The HIS / RIS 7 generates imaging order information in response to a registration operation by an operator based on an inquiry result or the like. The imaging order information includes patient information such as the name of the patient who is the subject, information on the imaging reservation such as the imaging part, imaging direction, body position (standing position, standing position), imaging method, necessity of analysis, analysis items, etc. Contains. Note that the imaging order information is not limited to that exemplified here, but may include other information, or may be a part of the information exemplified above.
解析用WS8は、CPU、RAM等により構成される制御部と、解析用プログラムを記憶する記憶部と、入力部と、表示部と、通信ネットワークNを介してコンソール5等の外部機器とデータ送受信を行うための通信部とを備えて構成されるワークステーションである。解析用WS8は、制御部と記憶部に記憶されている解析用プログラムとの協働により、コンソール5から送信された一連のフレーム画像に基づいて解析処理を行い、解析結果をコンソール5に送信する。 The analysis WS 8 transmits / receives data to / from an external device such as the console 5 through a communication network N, a control unit including a CPU, a RAM, a storage unit for storing an analysis program, an input unit, a display unit, and a communication network N. A workstation configured to include a communication unit. The analysis WS 8 performs analysis processing based on a series of frame images transmitted from the console 5 in cooperation with the control unit and the analysis program stored in the storage unit, and transmits the analysis result to the console 5. .
FPD9aは、パルス照射撮影可能な動態撮影及び静止画撮影対応の放射線検出器である。
図5に、FPD9aの機能構成例を示す。図5に示すように、FPD9aは、制御部91、検出部92、記憶部93、コネクタ94、バッテリ95、無線通信部96等を備えて構成され、各部はバス97により接続されている。
The FPD 9a is a radiation detector compatible with dynamic imaging and still image imaging capable of pulse irradiation imaging.
FIG. 5 shows a functional configuration example of the FPD 9a. As shown in FIG. 5, the FPD 9 a includes a control unit 91, a detection unit 92, a storage unit 93, a connector 94, a battery 95, a wireless communication unit 96, and the like, and each unit is connected by a bus 97.
制御部91は、CPU、RAM等により構成される。制御部91のCPUは、記憶部93に記憶されているシステムプログラムや処理プログラム等の各種プログラムを読み出してRAMに展開し、展開されたプログラムに従って各種処理を実行する。
例えば、制御部91は、コネクタ94を介して接続されたブッキー装置1又は2からの要求に応じてFPD9aの識別情報であるFPDIDを記憶部93から読み出して要求元のブッキー装置に送信する。
また、例えば、制御部91は、コンソール5から入力された画像読取条件に基づいて検出部92のスイッチング部を制御して、各放射線検出素子(以下、検出素子)に蓄積された電気信号の読み取りをスイッチングしていき、検出部92に蓄積された電気信号を読み取ることにより、画像データ(静止画像又はフレーム画像)を生成する。そして、制御部91は、生成した画像データを、順次コネクタ94及びブッキー装置1又は2を介してコンソール5に出力する。尚、撮影により取得された各フレーム画像は、一旦FPD9aの記憶部93に記憶され、全撮影終了後に、纏めてFPD9aからコンソール5に出力することとしても良い。
なお、FPD9aは、ブッキーに装填されない単体使用時には、バッテリ駆動及び無線通信する構成であるが、動態撮影の場合には、登録第4,561,730号公報に開示されているように、ブッキーを介した外部電力供給及び有線通信する構成に変更することが好ましい。これは、静止画撮影に比べ、データ転送容量や転送時間が圧倒的に増えるので、一のフレーム画像の転送中に他のフレーム画像の撮影(読取り)に対しノイズを与えないため、且つ、転送時間自体を短くするためであり、さらに、一連の撮影途中のバッテリ切れを防止するためである。
The control unit 91 is configured by a CPU, a RAM, and the like. The CPU of the control unit 91 reads out various programs such as system programs and processing programs stored in the storage unit 93, expands them in the RAM, and executes various processes according to the expanded programs.
For example, the control unit 91 reads FPDID, which is identification information of the FPD 9a, from the storage unit 93 in response to a request from the bucky device 1 or 2 connected via the connector 94, and transmits the FPDID to the requesting bucky device.
Further, for example, the control unit 91 controls the switching unit of the detection unit 92 based on the image reading condition input from the console 5 to read the electrical signal accumulated in each radiation detection element (hereinafter, detection element). The image data (still image or frame image) is generated by reading the electrical signal stored in the detection unit 92. Then, the control unit 91 sequentially outputs the generated image data to the console 5 via the connector 94 and the bucky device 1 or 2. Each frame image acquired by shooting may be temporarily stored in the storage unit 93 of the FPD 9a, and may be output from the FPD 9a to the console 5 after completion of all shooting.
Note that the FPD 9a is configured to be battery-driven and wirelessly communicated when used alone without being attached to the bucky, but in the case of dynamic shooting, as disclosed in Japanese Patent No. 4,561,730, It is preferable to change to a configuration for external power supply and wired communication. This is because data transfer capacity and transfer time are overwhelmingly increased compared to still image shooting, so that no noise is given to the shooting (reading) of another frame image during transfer of one frame image and transfer. This is to shorten the time itself, and to prevent the battery from running out during a series of photographing.
検出部92は、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源3a〜3cの何れかから照射されて少なくとも被写体を透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の検出素子が二次元状に配列されている。検出素子は、フォトダイオード等の半導体イメージセンサにより構成される。各検出素子は、例えばTFT(Thin Film Transistor)等のスイッチング部に接続され、スイッチング部により電気信号の蓄積及び読み出しが制御される。
なお、生成された静止画像又はフレーム画像を構成する各画素は、検出部92の各検出素子のそれぞれから出力された信号値(ここでは、濃度値と呼称する)を示す。
The detection unit 92 includes, for example, a glass substrate, and detects radiation that has been irradiated from any of the radiation sources 3a to 3c and transmitted through at least the subject at a predetermined position on the substrate according to the intensity thereof. A plurality of detection elements that convert detected radiation into electrical signals and store them are two-dimensionally arranged. The detection element is configured by a semiconductor image sensor such as a photodiode. Each detection element is connected to a switching unit such as a TFT (Thin Film Transistor), for example, and the switching unit controls the accumulation and reading of electric signals.
Each pixel constituting the generated still image or frame image indicates a signal value (referred to herein as a density value) output from each detection element of the detection unit 92.
記憶部93は、例えば半導体の不揮発性メモリ等で構成されている。記憶部93には、検出部92を制御するための各種プログラムや自己の識別情報であるFPDID等が記憶されている。また、記憶部93には、検出部92から出力された画像データが一時的に記憶される。 The storage unit 93 is configured by, for example, a semiconductor nonvolatile memory. The storage unit 93 stores various programs for controlling the detection unit 92, FPDID that is identification information of itself, and the like. The storage unit 93 temporarily stores the image data output from the detection unit 92.
コネクタ94は、ブッキー装置1、2側のコネクタと接続し、ブッキー装置1又は2とのデータ送受信を行う。また、コネクタ94は、ブッキー装置1又は2のコネクタから供給される電力を各機能部へ供給する。なお、バッテリ95を充電する構成としても良い。 The connector 94 is connected to the connector on the Bucky device 1 or 2 side and performs data transmission / reception with the Bucky device 1 or 2. The connector 94 supplies power supplied from the connector of the bucky device 1 or 2 to each function unit. Note that the battery 95 may be charged.
バッテリ95は、制御部91の制御に基づいて、FPD9aの各部に電力を供給する。バッテリ95としては、例えばニッカド電池、ニッケル水素電池、リチウムイオン電池等の充電自在な電池等を適用することができる。 The battery 95 supplies power to each unit of the FPD 9a based on the control of the control unit 91. As the battery 95, for example, a rechargeable battery such as a nickel cadmium battery, a nickel metal hydride battery, or a lithium ion battery can be used.
FPD9bは、FPD9aと同様に、制御部91、検出部92、記憶部93、コネクタ94、バッテリ95を備えて構成されるが、フレームレートの設定はできない。即ち、FPD9bは、静止画撮影のみが可能である。
また、単体での使用は勿論のこと、ブッキー装置に装填しても使用可能で、ブッキー装置に装填時には、コネクタ接続により、バッテリ/無線方式から、有線/電力供給方式に切り替えることができる。従い、複数の患者を連続的に静止画撮影する場合に於いても、バッテリ切れを気にする必要がなくなる。
Like the FPD 9a, the FPD 9b includes a control unit 91, a detection unit 92, a storage unit 93, a connector 94, and a battery 95, but the frame rate cannot be set. That is, the FPD 9b can only take still images.
In addition, it can be used as a stand-alone unit, or can be used by being loaded into a Bucky device. When the Bucky device is loaded, it can be switched from a battery / wireless system to a wired / power supply system by connecting a connector. Accordingly, even when a plurality of patients are continuously photographed, it is not necessary to worry about running out of the battery.
PACS10は、画像データ等を保存するサーバ装置と、当該サーバ装置から診断用の画像を取得して表示するための読影端末とを備えて構成される。PACS10のサーバ装置は、コンソール5から送信された画像データや解析結果のデータを撮影オーダ情報に対応付けて記憶する。 The PACS 10 includes a server device that stores image data and the like, and an image interpretation terminal that acquires and displays a diagnostic image from the server device. The server device of the PACS 10 stores image data and analysis result data transmitted from the console 5 in association with imaging order information.
次に、診断支援情報生成システム100における撮影動作について説明する。
図6に、診断支援情報生成システム100において実行される撮影・解析処理の流れを示す。図6のコンソール5側の処理は、コンソール5の制御部51と記憶部52に記憶されているプログラムとの協働により実行される。解析用WS8側の処理は、解析用WS8の制御部と記憶部に記憶されている解析プログラムとの協働により実行される。
Next, the imaging operation in the diagnosis support information generation system 100 will be described.
FIG. 6 shows a flow of imaging / analysis processing executed in the diagnostic support information generation system 100. The processing on the console 5 side in FIG. 6 is executed in cooperation with the control unit 51 of the console 5 and the program stored in the storage unit 52. The processing on the analysis WS8 side is executed in cooperation with the control unit of the analysis WS8 and the analysis program stored in the storage unit.
まず、撮影技師等の操作者は、何れかの撮影室のコンソール5の入力部53を操作して撮影オーダ情報の一覧を表示する撮影オーダリスト画面を表示部54に表示させる。そして、入力部53を操作することにより撮影オーダリスト画面から撮影対象の撮影オーダ情報を指定する。 First, an operator such as a photographing engineer operates the input unit 53 of the console 5 in any of the photographing rooms to display a photographing order list screen for displaying a list of photographing order information on the display unit 54. Then, by operating the input unit 53, the imaging order information of the imaging target is designated from the imaging order list screen.
コンソール5においては、入力部53により撮影対象の撮影オーダ情報が指定されると(ステップS1)、記憶部52の撮影管理テーブル521が参照され、当該コンソール5が設置されている撮影室で選択された撮影オーダ情報に基づく撮影が可能であるか否かが判断される(ステップS2)。例えば、撮影オーダ情報により動態撮影が指示されていた場合、撮影管理テーブル521が参照され、その撮影室に、連射が可能な管球タイプの放射線源及び動態撮影対応のFPDが存在し、かつ、そのFPDが使用中ではない場合(画像受信時刻から予め定められた時間経過している場合)に、動態撮影が可能であると判断される。 In the console 5, when shooting order information to be shot is specified by the input unit 53 (step S1), the shooting management table 521 of the storage unit 52 is referred to and selected in the shooting room in which the console 5 is installed. It is determined whether or not shooting based on the shooting order information is possible (step S2). For example, when dynamic imaging is instructed by the imaging order information, the imaging management table 521 is referred to, and there is a tube-type radiation source capable of continuous shooting and an FPD corresponding to dynamic imaging in the imaging room, and When the FPD is not in use (when a predetermined time has elapsed from the image reception time), it is determined that dynamic imaging is possible.
当該コンソール5が設置されている撮影室において、選択された撮影オーダ情報に基づく撮影が可能であると判断されると(ステップS2;YES)、処理はステップS4に移行する。 If it is determined that shooting based on the selected shooting order information is possible in the shooting room in which the console 5 is installed (step S2; YES), the process proceeds to step S4.
当該撮影室では選択された撮影オーダ情報に基づく撮影が不可能であると判断されると(ステップS2:NO)、表示部54に警告が表示される(ステップS3)。例えば、撮影オーダ情報により動態撮影が指示されているにもかかわらず、連射可能な管球が撮影室に存在していないと判断されると(図1の撮影室R2のような場合)、「この撮影室では撮影できません」等の警告が表示される。また、例えば、撮影オーダ情報により指定された***を動態撮影するためのブッキー装置に動態撮影対応のFPD9aが装着されていないと判断されると、「立位(臥位)用のブッキー装置にFPD9aを装着してください」等の警告が表示される。そして、撮影・解析処理は終了する。なお、撮影オーダ情報により指定された***を動態撮影するためのブッキー装置に動態撮影対応のFPD9aが装着されていない場合、動態撮影対応のFPDを該当するブッキー装置に装着すれば、撮影管理テーブル521の内容は更新される。従って、撮影可能と判断され、処理はステップS4に移行する。
なお、この場合、一旦、撮影室内でFPDの入替操作を行って、再度コンソールまで戻り、再度ステップS1の処理から行わせ、確実性を高めるフローとしても良い。
If it is determined that shooting based on the selected shooting order information is impossible in the shooting room (step S2: NO), a warning is displayed on the display unit 54 (step S3). For example, when dynamic shooting is instructed by the shooting order information, it is determined that there is no continuous shooting tube in the shooting room (such as the shooting room R2 in FIG. 1). A warning such as “Cannot shoot in this room” is displayed. Further, for example, if it is determined that the FPD 9a corresponding to dynamic imaging is not attached to the Bucky device for performing dynamic imaging of the body position specified by the imaging order information, the FPD 9a Please put on "warning. Then, the photographing / analysis process ends. If the FPD 9a corresponding to dynamic imaging is not attached to the Bucky device for dynamic imaging of the posture specified by the imaging order information, the imaging management table 521 can be provided by attaching an FPD corresponding to dynamic imaging to the corresponding Bucky device. The contents of are updated. Therefore, it is determined that photographing is possible, and the process proceeds to step S4.
In this case, the FPD replacement operation may be performed once in the photographing room, and the process may return to the console again, and the process may be performed again from step S1 to increase the certainty.
ステップS4においては、指定された撮影オーダ情報の撮影が可能な放射線源及びFPDが起動され、使用されるブッキー装置に応じて放射線源の向き及び位置が調整される。撮影技師により被写体に応じてFPDやブッキー装置の位置等が調整されると、それに応じて放射線源の向き及び位置が調整される(ステップS4)。また、記憶部52から撮影する部位や動態撮影か静止画撮影かの別に応じた放射線照射条件及び画像読取条件が読み出され、放射線源に放射線照射条件が設定されるとともに、ブッキー装置を介してFPDに画像読取条件が設定される(ステップS5)。動態撮影の結果を用いて解析を行う場合は、診断に使用し得る解析精度を確保するため、フレームレートが3.75枚/秒以上に設定される。ここで、操作者は、肺野の動態撮影の場合、安静呼吸の動態を撮影するために被験者に楽にするように指示し、安静呼吸を続けるよう促す。撮影準備が整った時点で、技師は前室へ移動し、操作卓6を操作して放射線照射指示を入力する。 In step S4, the radiation source and FPD capable of imaging the designated imaging order information are activated, and the orientation and position of the radiation source are adjusted according to the used Bucky device. When the position of the FPD or the bucky device is adjusted according to the subject by the imaging engineer, the direction and position of the radiation source are adjusted accordingly (step S4). In addition, the radiation irradiation condition and the image reading condition corresponding to the part to be imaged, dynamic imaging or still image imaging are read out from the storage unit 52, the radiation irradiation condition is set in the radiation source, and via the Bucky device Image reading conditions are set in the FPD (step S5). When analysis is performed using the result of dynamic imaging, the frame rate is set to 3.75 frames / second or more in order to ensure analysis accuracy that can be used for diagnosis. Here, in the case of lung field dynamic imaging, the operator instructs the subject to make the rest breathing easier in order to capture the rest breathing dynamics, and urges the patient to continue the rest breathing. When preparation for imaging is completed, the technician moves to the front room and operates the console 6 to input a radiation irradiation instruction.
操作卓6からの放射線照射指示が入力されると(ステップS6;YES)、撮影に使用される放射線源及びFPDが制御され、撮影が行われる(ステップS7)。
動態撮影の場合は、ステップS5で設定されたパルス間隔で放射線源3aにより放射線が照射され、ステップS5で設定されたフレームレートでFPD9aによりフレーム画像が取得される。予め定められたフレーム画像数の撮影が終了すると、制御部51により放射線源3a及びFPD9aに撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム画像数は、少なくとも動態1サイクルが撮影できる枚数である。撮影により取得されたフレーム画像はFPD9aからブッキー装置を介して順次コンソール5に入力される。
なお、撮影オーダ情報により解析が指示されていない場合には、オフセット補正用のダーク画像を読取り、コンソール5に入力することとしても良い。
静止画撮影の場合は、ステップS5で設定された条件で被写体の1枚の静止画像及びオフセット補正用の1または数枚のダーク画像が撮影される。撮影により取得された静止画像及びダーク画像は、FPDからブッキー装置を介してコンソール5に入力される。
When a radiation irradiation instruction is input from the console 6 (step S6; YES), the radiation source and FPD used for imaging are controlled and imaging is performed (step S7).
In the case of dynamic imaging, radiation is emitted from the radiation source 3a at the pulse interval set in step S5, and a frame image is acquired by the FPD 9a at the frame rate set in step S5. When photographing of a predetermined number of frame images is completed, the control unit 51 outputs a photographing end instruction to the radiation source 3a and the FPD 9a, and the photographing operation is stopped. The number of frame images to be photographed is the number that can be photographed in at least one dynamic cycle. Frame images obtained by photographing are sequentially input from the FPD 9a to the console 5 via the bucky device.
If analysis is not instructed by the shooting order information, a dark image for offset correction may be read and input to the console 5.
In the case of still image shooting, one still image of the subject and one or several dark images for offset correction are shot under the conditions set in step S5. The still image and the dark image acquired by photographing are input to the console 5 from the FPD via the bucky device.
次いで、解析用WS8で解析を行うか否かが判断される(ステップS8)。ここで、解析用WS8において解析を行うか否かの判断は、例えば、ステップS1で指定された撮影オーダ情報に基づいて判断される。撮影オーダ情報により静止画撮影がオーダされていれば、解析は行わないと判断される。動態撮影がオーダされている場合は、撮影オーダ情報により解析が必要である旨を示す情報が含まれている場合に、解析が必要であると判断される。 Next, it is determined whether or not the analysis is performed by the analysis WS 8 (step S8). Here, the determination as to whether or not to perform analysis in the analysis WS 8 is made based on, for example, the imaging order information specified in step S1. If still image shooting is ordered based on the shooting order information, it is determined that analysis is not performed. When dynamic imaging is ordered, it is determined that analysis is necessary when the imaging order information includes information indicating that analysis is necessary.
解析用WS8で解析を行わないと判断されると(ステップS8;NO)、撮影により得られた画像に対して補正処理が行われ(ステップS9)、処理はステップS10に移行する。ステップS9の補正処理においては、上述のダーク画像を用いたオフセット補正処理、ゲイン補正処理、欠陥画素補正処理、ラグ(残像)補正処理等の補正処理が必要に応じて行われる。解析を行う場合は、処理時間の短縮化優先するためにこれらの補正処理を省略しても良く、処理はステップS10に移行する。
本件発明者等は、動態解析においては、静止画のような個々の画素の絶対的出力値はあまり重要ではなく、個々の画素におけるフレーム間の相対的出力値(変動成分)に基づく特徴量の算出が基本となっているため、上記の補正処理の一部或いは全部を省略しても、補正処理を行った場合と略同等の解析結果を得ることが可能である知見を得た。従い、解析結果を得るまでの時間を短縮するため、一部或いは全部の補正処理を省略することできる。
If it is determined that analysis is not performed by the analysis WS 8 (step S8; NO), a correction process is performed on the image obtained by photographing (step S9), and the process proceeds to step S10. In the correction process in step S9, correction processes such as the offset correction process using the dark image, the gain correction process, the defective pixel correction process, and the lag (afterimage) correction process are performed as necessary. When performing analysis, these correction processes may be omitted in order to give priority to shortening the processing time, and the process proceeds to step S10.
In the dynamic analysis, the inventors of the present invention are not so important in the absolute output value of each pixel such as a still image, and the feature amount based on the relative output value (variation component) between frames in each pixel. Since the calculation is fundamental, it has been found that even if a part or all of the correction process is omitted, an analysis result substantially equivalent to that obtained when the correction process is performed can be obtained. Accordingly, some or all of the correction processing can be omitted in order to shorten the time until the analysis result is obtained.
ステップS10においては、撮影により得られたフレーム画像又は静止画像が撮影オーダ情報と対応付けて記憶部52に記憶される(ステップS10)。なお、撮影に得られたフレーム画像には、撮影順を示す番号が付与され、各画像のヘッダ情報等に記憶される。 In step S10, a frame image or a still image obtained by shooting is stored in the storage unit 52 in association with shooting order information (step S10). It should be noted that the frame image obtained for shooting is given a number indicating the shooting order and stored in the header information of each image.
次いで、順次入力された画像が間引き処理され、表示部54に表示される(ステップS11)。ここでいう間引き処理は、各フレーム画像や静止画像の画素数を減らす処理を指す。例えば、所定の画素間隔毎の画素で構成された間引き画像を作成する処理(単純間引き処理と呼ぶ)と、フレーム画像を所定サイズの画素ブロック単位、例えば、2mm×2mm角単位の小領域に分割し、各小領域内の画素の信号値の代表値(ここでは、平均信号値とする)を算出して小領域内の画素の信号値を算出された代表値に置き換えるビニング処理が含まれる。ビニング処理では、各小領域単位を一つの画素とみなして取り扱うことで、処理対象の画素数を低減することができる。ビニング処理では、動態画像の場合、各フレーム画像間の対応する各小領域は検出素子の同じ位置の出力を示す画素群からなるように分割が行われる。例えば、フレーム画像上の同一画素位置(0、0)を基点として2mm×2mm角で分割が行われる。なお、ビニング処理における各画素ブロックのサイズは、診断対象(即ち、解析対象)となる撮影部位に応じたサイズとすることが好ましい。また、後段で解析処理を行う場合には、解析により算出される特徴量に応じたサイズとすることが好ましい。 Next, the sequentially input images are thinned out and displayed on the display unit 54 (step S11). The thinning-out processing here refers to processing for reducing the number of pixels of each frame image or still image. For example, a process for creating a thinned image composed of pixels at a predetermined pixel interval (referred to as simple thinning process) and a frame image are divided into pixel blocks of a predetermined size, for example, a small area of 2 mm × 2 mm square unit In addition, a binning process is included in which a representative value (here, an average signal value) of the signal values of the pixels in each small area is calculated and the signal values of the pixels in the small area are replaced with the calculated representative values. In the binning process, the number of pixels to be processed can be reduced by treating each small area unit as one pixel. In the binning process, in the case of a dynamic image, division is performed so that each corresponding small region between each frame image is composed of a pixel group indicating an output at the same position of the detection element. For example, the division is performed by 2 mm × 2 mm square with the same pixel position (0, 0) on the frame image as a base point. In addition, it is preferable that the size of each pixel block in the binning process is a size corresponding to an imaging region that is a diagnosis target (that is, an analysis target). Moreover, when performing an analysis process in the latter stage, it is preferable to set the size according to the feature amount calculated by the analysis.
次いで、間引きされたフレーム画像が表示部54に表示される(ステップS12)。
撮影技師は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された(撮影OK)か、再撮影が必要(撮影NG)か、を判断する。そして、入力部53を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、一旦FPD9aの記憶部93に記憶され、全撮影の終了後に纏めてFPD9aからコンソール5へ出力されるようにしても良い。
Next, the thinned frame image is displayed on the display unit 54 (step S12).
The imaging engineer confirms the positioning and the like based on the displayed dynamic image, and determines whether an image suitable for diagnosis is acquired by imaging (imaging OK) or re-imaging is necessary (imaging NG). Then, the input unit 53 is operated to input the determination result. Each frame image acquired by photographing may be temporarily stored in the storage unit 93 of the FPD 9a, and may be output from the FPD 9a to the console 5 after completion of all photographing.
入力部53の所定の操作により撮影NGを示す判断結果が入力されると(ステップS13;NO)、記憶部52に記憶された一連のフレーム画像が削除され(ステップS14)、本処理は終了する。なお、この場合、再撮影が行われることとなる。
入力部53の所定の操作により撮影OKを示す判断結果が入力されると(ステップS13;YES)、解析を行うか否かが判断される(ステップS15)。解析を行うか否かの判断は、例えば、ステップS8で説明したのと同様の判断により行われる。解析を行わないと判断されると(ステップS15;NO)、撮影された静止画像又はフレーム画像が必要に応じて画像処理され、ネットワーク通信部56を介してPACS10のサーバ装置に送信される(ステップS16)。なお、PACS10のサーバ装置においては、受信された静止画像又はフレーム画像が撮影オーダ情報に対応付けて保存される。
When a determination result indicating imaging NG is input by a predetermined operation of the input unit 53 (step S13; NO), a series of frame images stored in the storage unit 52 is deleted (step S14), and this process ends. . In this case, re-photographing is performed.
When a determination result indicating photographing OK is input by a predetermined operation of the input unit 53 (step S13; YES), it is determined whether or not to perform analysis (step S15). The determination as to whether or not to perform the analysis is performed based on the same determination as described in step S8, for example. If it is determined not to perform analysis (step S15; NO), the captured still image or frame image is subjected to image processing as necessary and transmitted to the server device of the PACS 10 via the network communication unit 56 (step S15). S16). In the server device of PACS 10, the received still image or frame image is stored in association with the shooting order information.
一方、解析を行うと判断されると(ステップS15;YES)、撮影により得られた一連のフレーム画像群から解析に使用するフレーム画像群を選択する処理が行われる(ステップS17)。 On the other hand, if it is determined that analysis is to be performed (step S15; YES), a process of selecting a frame image group to be used for analysis from a series of frame image groups obtained by photographing is performed (step S17).
ここで、ステップS17の処理について説明する。ステップS17においては、まず、表示部54に、撮影により得られた一連のフレーム画像のサムネイル画像が撮影順に並べて表示された選択用画面541が表示される。図7に、選択用画面541の一例を示す。図7に示すように、選択用画面541には、一連のフレーム画像が並べて表示されているとともに、選択枠541aが表示されている。 Here, the process of step S17 will be described. In step S <b> 17, first, a selection screen 541 on which thumbnail images of a series of frame images obtained by shooting are arranged and displayed in the shooting order is displayed on the display unit 54. FIG. 7 shows an example of the selection screen 541. As shown in FIG. 7, the selection screen 541 displays a series of frame images side by side and a selection frame 541a.
解析に必要となる枚数以上のフレーム画像を用いて解析を行うと、データ転送時間や解析処理時間がかかり、好ましくない。処理時間短縮のためには、撮影された一連のフレーム画像の中から解析に使用する枚数のフレーム画像を選択することが必要となる。一方、動態画像の解析は、後述するように、主として動態周期や隣接するフレーム画像間での信号値の差分等に基づき行われるので、動態画像の解析結果を被写体の真の特徴量を表す結果とするには、被写体の動態一周期以上の枚数の、連続する一連のフレーム画像が必要である。そこで、ステップS17においては、解析に必要な枚数分の連続するフレーム画像を操作者が正しく選択するためのGUIとして選択用画面541を提供する。具体的には、選択枠541aを用いて選択させるようにすることで、連続しない1コマ毎の画像選択や、一周期に満たない選択を防止する。 If analysis is performed using more frame images than necessary for analysis, it takes data transfer time and analysis processing time, which is not preferable. In order to shorten the processing time, it is necessary to select the number of frame images to be used for analysis from a series of photographed frame images. On the other hand, as will be described later, dynamic image analysis is performed mainly based on a dynamic cycle, a difference in signal values between adjacent frame images, and the like. In order to achieve this, a continuous series of frame images having a number equal to or greater than one cycle of the subject's dynamics is required. Therefore, in step S17, a selection screen 541 is provided as a GUI for the operator to correctly select the continuous frame images for the number of sheets necessary for analysis. Specifically, by using the selection frame 541a to select, it is possible to prevent image selection for each frame that is not continuous or selection that is less than one cycle.
選択枠541aを表示するために、まず、解析に必要な画像枚数nが算出される。例えば、記憶部52に平均的な成人等の動態周期を部位毎に記憶しておき、撮影に利用された画像読取条件に含まれるフレームレートと、被写体部位の動態周期とに基づいて、解析に必要な画像枚数nが算出される。例えば、肺の換気についての解析を行うのであれば、平均的な成人であれば呼吸周期が3.3秒前後であるので、フレームレートが5枚/秒であれば20枚、フレームレートが7.5枚/秒であれば30枚程度必要となる。解析に必要な枚数nが算出されると、連続するn枚のフレーム画像のサムネイル画像を囲う選択枠541aが表示される。デフォルトでは、撮影順番が1番〜n番までが選択枠541aに囲まれて表示される。選択用画面541には、次ページボタン541a、最終ページボタン541c、確定ボタン541d等が設けられている。次ページボタン541aが押下されると、次のページが表示される。最終ページボタン541cが押下されると、最終ページが表示される。操作者により一のサムネイル画像が選択されると、選択されたフレーム画像を先頭としてn枚の画像を囲った選択枠541aが表示される。確定ボタン541dが押下されると、選択枠541aで囲まれたサムネイル画像のフレーム画像群が解析に用いる画像として選択される。
なお、選択されたフレーム画像を先頭とすると必要な枚数が確保できない場合がある。このような場合は、選択用画面541上に警告をポップアップ表示する。このポップアップ画面でフレーム画像数不足のまま解析を行うか、選択しなおすかをユーザに選択させるようにしてもよい。また、本実施の形態においては、解析に使用する始点となるフレーム画像を選択させることとしているが、最終画像をユーザに選択させる構成としてもよい。
In order to display the selection frame 541a, first, the number n of images necessary for analysis is calculated. For example, the dynamic cycle of an average adult or the like is stored for each region in the storage unit 52, and the analysis is performed based on the frame rate included in the image reading conditions used for imaging and the dynamic cycle of the subject region. The required number of images n is calculated. For example, when analyzing the ventilation of the lung, the breathing cycle is around 3.3 seconds for an average adult, so 20 frames when the frame rate is 5 frames / second, and a frame rate of 7 If it is 5 sheets / second, about 30 sheets are required. When the number n required for analysis is calculated, a selection frame 541a surrounding thumbnail images of consecutive n frame images is displayed. By default, shooting order numbers 1 to n are displayed surrounded by a selection frame 541a. The selection screen 541 is provided with a next page button 541a, a last page button 541c, a confirmation button 541d, and the like. When the next page button 541a is pressed, the next page is displayed. When the last page button 541c is pressed, the last page is displayed. When one thumbnail image is selected by the operator, a selection frame 541a surrounding n images starting from the selected frame image is displayed. When the confirm button 541d is pressed, a frame image group of thumbnail images surrounded by the selection frame 541a is selected as an image used for analysis.
If the selected frame image is the head, there may be a case where a necessary number cannot be secured. In such a case, a warning is displayed in a pop-up on the selection screen 541. On this pop-up screen, the user may be allowed to select whether to perform analysis while the number of frame images is insufficient or to select again. In the present embodiment, the frame image that is the starting point used for the analysis is selected. However, the final image may be selected by the user.
なお、解析項目によっては、連続する動態周期以上の一連のフレーム画像とともに、例えば肺野であれば最大呼気、最大吸気の画像が必要となる場合もある。このような場合には、選択枠541aによる選択に加えて、単発でフレーム画像を選択することを可能とする。例えば、入力部53のマウスによりフレーム画像のサムネイル画像をクリックすると、その画像が始点となるように選択枠541aが移動し、ダブルクリックすると単発でフレーム画像が選択される。 Depending on the analysis item, for example, in the lung field, an image of maximum exhalation and maximum inspiration may be required together with a series of frame images having a continuous dynamic period or more. In such a case, in addition to selection by the selection frame 541a, it is possible to select a frame image in a single shot. For example, when a thumbnail image of a frame image is clicked with the mouse of the input unit 53, the selection frame 541a is moved so that the image becomes the start point, and when a double-click is performed, the frame image is selected in a single shot.
一連のフレーム画像群から解析に使用するフレーム画像群が選択されると、選択されたフレーム画像の間引きデータが撮影オーダ情報と対応付けられ、ネットワーク通信部56を介して解析用WS8に送信される(ステップS18)。例えば、選択された一連のフレーム画像の間引きデータのそれぞれに、動態画像を識別するための識別IDや、患者情報、撮影部位、放射線照射条件、画像読取条件(フレーム間隔等)、撮影順を示す番号、フレーム枚数、撮影日時等の情報が付帯され(例えば、DICOM(Digital Imaging and Communications in Medicine)マルチフォーマットのファイルフォーマットで画像データのヘッダ領域に書き込まれ)、ネットワーク通信部56を介して解析用WS8に送信される。また、解析項目も併せて解析用WS8に通知される。 When a frame image group to be used for analysis is selected from a series of frame image groups, the thinned data of the selected frame image is associated with the imaging order information and transmitted to the analysis WS 8 via the network communication unit 56. (Step S18). For example, identification data for identifying a dynamic image, patient information, imaging region, radiation irradiation condition, image reading condition (frame interval, etc.), and imaging order are indicated for each of a series of selected frame image thinning data. Information such as the number, number of frames, shooting date and time is attached (for example, written in the header area of the image data in DICOM (Digital Imaging and Communications in Medicine) multi-format file format), and for analysis via the network communication unit 56 Sent to WS8. The analysis items are also notified to the analysis WS 8.
解析用WS8においては、受信したフレーム画像群に基づいて解析処理が行われる(ステップS19)。
解析処理の内容は、撮影部位毎に異なる。ここでは、肺野の解析を行う場合を例にとり説明する。
In the analysis WS 8, analysis processing is performed based on the received frame image group (step S 19).
The content of the analysis process differs for each imaging region. Here, a case where lung field analysis is performed will be described as an example.
肺野の解析としては、肺野の局所的な動きを示す特徴量を算出する解析と、肺野全体の動きを示す特徴量を算出する解析とがある。また、換気機能を対象とする解析と、血流機能を対象とする解析とがある。 The lung field analysis includes an analysis for calculating a feature value indicating local movement of the lung field and an analysis for calculating a feature value indicating the movement of the entire lung field. Further, there are an analysis for the ventilation function and an analysis for the blood flow function.
肺野の換気機能による局所的な信号変化を示す特徴量を算出する解析としては、例えば、下記の(1)〜(6)の項目が挙げられる。以下、各特徴量毎に、その算出方法の手順について簡単に説明する。なお、下記においては、解析に必要な処理を明確にするため、撮影されたフレーム画像の生データ(RAWデータ)から特徴量を算出する処理内容について説明するが、本実施の形態においては、解析用WS8へのフレーム画像送信時のデータ量の削減及び処理時間の短縮化の観点から、すでに各フレーム画像に、間引き処理の1種であるビニング処理が施されている(即ち、所定サイズの小領域に分割され、小領域毎に信号値の平均化が行われている)。 For example, the following items (1) to (6) are listed as the analysis for calculating the feature quantity indicating the local signal change due to the ventilation function of the lung field. Hereinafter, the procedure of the calculation method for each feature amount will be briefly described. In the following, in order to clarify the processing necessary for the analysis, the processing content for calculating the feature amount from the raw data (RAW data) of the captured frame image will be described, but in the present embodiment, the analysis is performed. From the viewpoint of reducing the amount of data when transmitting a frame image to the WS 8 and shortening the processing time, each frame image has already been subjected to a binning process which is a kind of thinning process (that is, a small size of a predetermined size). It is divided into regions, and signal values are averaged for each small region).
解析用WS8の解析処理においては、従来のようにワーピング処理を施してフレーム画像間の画素を対応付けて特徴量を算出するのではなく、ワーピング処理は施さずにFPD9aにおける同一位置の検出素子の出力を示す画素を互いに対応付けて特徴量の算出を行い、特徴量の精度を維持したまま処理時間を大幅に短縮した点に特徴がある。 In the analysis process of the analysis WS 8, the warping process is not performed and the feature amount is calculated by associating the pixels between the frame images as in the related art. The feature amount is calculated by associating pixels indicating outputs with each other, and the processing time is significantly shortened while maintaining the accuracy of the feature amount.
(1)換気−フレーム間差分画像
フレーム間差分画像は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→フレーム間差分処理→ノイズ除去処理
ビニング処理は、上述のように、各フレーム画像において、画像領域を所定サイズの画素ブロック単位の小領域に分割し、小領域毎に領域内画素の信号値の代表値、例えば、平均値を算出する(平均化する)処理である。なお、代表値としては、平均値に限らず、中央値、平均値、最頻値としてもよい。画素ブロックのサイズは、解析対象となる部位、及び/又は解析により算出される特徴量に応じたものとすることが解析精度向上の点から好ましい。
時間軸方向のローパスフィルタ処理は、換気による信号値の時間変化を抽出するための処理であり、例えば、カットオフ周波数0.5Hzでフィルタリングする。
フレーム間差分処理は、一連のフレーム画像の同じ画素位置の小領域(FPDの同じ位置の検出素子から出力された領域)を互いに対応付け、各小領域毎に、隣接するフレーム画像間で信号値の差分値を算出し、フレーム間差分画像を作成する処理である。
なお、フレーム間差分画像の静止画像を作成する場合は、肺野全体の濃度変化又は横隔膜の位置の変化を解析することにより一連のフレーム画像における吸気期間と呼気期間を算出し、小領域毎に吸気期間については正のフレーム間差分値の絶対値を積算し、呼気期間については負のフレーム間差分値の絶対値を積算した画像を作成する。
(1) Ventilation-interframe difference image An interframe difference image is calculated by performing the following processing on a series of frame images.
Binning process → Low-pass filter process in the time axis direction → Inter-frame difference process → Noise removal process As described above, the binning process divides the image area into small areas in units of pixel blocks of a predetermined size in each frame image. This is a process of calculating (averaging) a representative value, for example, an average value, of signal values of pixels in the area for each area. The representative value is not limited to the average value, and may be a median value, an average value, or a mode value. The size of the pixel block is preferably in accordance with the part to be analyzed and / or the feature amount calculated by the analysis from the viewpoint of improving the analysis accuracy.
The low-pass filter process in the time axis direction is a process for extracting a time change of a signal value due to ventilation, and for example, filtering with a cutoff frequency of 0.5 Hz.
The inter-frame difference processing associates small regions at the same pixel position (regions output from detection elements at the same position of the FPD) of a series of frame images with each other, and signal values between adjacent frame images for each small region. The difference value is calculated, and an inter-frame difference image is created.
When creating a still image of the inter-frame difference image, the inspiratory period and expiratory period in a series of frame images are calculated by analyzing the change in the concentration of the entire lung field or the change in the position of the diaphragm. For the inhalation period, an absolute value of the positive inter-frame difference value is integrated, and for the expiration period, an image of the absolute value of the negative inter-frame difference value is integrated.
(2)換気−波形描画
換気−波形描画は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→波形描画
ビニング処理及び時間軸方向のローパスフィルタ処理は上述のとおりである(以下同じ)。
波形描画処理は、一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、横軸を撮影開始からの経過時間、縦軸を画素の平均信号値とした座標平面を作成して、各フレーム画像の撮影開始からの経過時間とその小領域について算出された平均信号値が交わる点をプロットすることより、換気量を示す信号値の時間変化を示す波形を描画する処理である。
(2) Ventilation-waveform drawing Ventilation-waveform drawing is calculated by performing the following processing on a series of frame images.
Binning process → low-pass filter process in time axis direction → waveform drawing The binning process and the low-pass filter process in time axis direction are as described above (hereinafter the same).
The waveform drawing process associates small regions at the same pixel position (regions of pixel blocks output from detection elements at the same position of the FPD) of a series of frame images with each other, and the horizontal axis for each small region starts from the start of imaging. By creating a coordinate plane with the vertical axis representing the average signal value of the pixel and plotting the points where the elapsed time from the start of capturing each frame image and the average signal value calculated for that small region intersect This is a process of drawing a waveform indicating a time change of a signal value indicating a ventilation amount.
(3)換気−気流速度
気流速度は、各小領域の肺のやわらかさ(肺コンプライアンス)を示す特徴量である。気流速度は、一連のフレーム画像に、以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→フレーム間差分処理→フレーム間差分値の代表値(最大値もしくは平均値)を算出
(3) Ventilation-Airflow velocity The airflow velocity is a feature amount indicating the softness of the lungs (lung compliance) in each small area. The air velocity is calculated by performing the following processing on a series of frame images.
Binning processing → Low-pass filter processing in the time axis direction → Difference processing between frames → Calculate the representative value (maximum value or average value) of difference values between frames
なお、最大値を代表値とする場合、各小領域毎に、呼気期間と吸気期間のそれぞれの期間における気流速度の最大値を表す指標として、呼気期間、吸気期間における信号変化(フレーム間差分値)の最大値をそれぞれ算出し、その比(最大流速比)の分布を示すヒストグラムを作成するとともに、何れか1つのフレーム画像上の各小領域を最大流速比に応じた輝度もしくは色で示す画像を作成し、両者を並べて解析結果として提供してもよい。この解析を最大流速比のヒストグラム解析と呼ぶ。
最大流速比のヒストグラム解析では、図11A〜図11C等に示すように、各小領域の吸気気流速度の最大値(絶対値)と呼気気流速度の最大値(絶対値)との比の値がヒストグラム表示されるとともに、COPDであるか否かの指標となる肺野全体での平均値や標準偏差が表示された画像が生成される。また、併せて静止画像上の各小領域を比の値に応じた輝度もしくは色で示すことで、異常個所の分布を医師が容易に把握できるような診断情報提供する。
When the maximum value is used as a representative value, the signal change (difference value between frames) is used as an index representing the maximum value of the airflow velocity in each of the expiration period and the inspiration period for each small area. ) Respectively, and a histogram showing the distribution of the ratio (maximum flow rate ratio) is created, and each small region on any one of the frame images is displayed in luminance or color according to the maximum flow rate ratio. May be prepared, and both may be arranged and provided as an analysis result. This analysis is called histogram analysis of the maximum flow rate ratio.
In the histogram analysis of the maximum flow rate ratio, as shown in FIGS. 11A to 11C, the value of the ratio between the maximum value (absolute value) of the inspiratory airflow velocity and the maximum value (absolute value) of the expiratory airflow velocity in each small region is In addition to the histogram display, an image is generated in which the average value and standard deviation of the entire lung field, which is an indicator of whether or not the image is COPD, are displayed. In addition, by displaying each small region on the still image with luminance or color according to the ratio value, diagnostic information is provided so that the doctor can easily understand the distribution of abnormal parts.
(4)換気量の振幅
換気量の振幅は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、呼吸1サイクル中の最大信号値(極大値)−最小信号値(極小値)の算出
(4) Ventilation Amplitude The ventilation volume amplitude is calculated by performing the following processing on a series of frame images.
Binning processing → Low-pass filter processing in the time axis direction → Small regions at the same pixel position (regions of pixel blocks output from detection elements at the same position in the FPD) of a series of frame images are associated with each other, Calculation of maximum signal value (maximum value)-minimum signal value (minimum value) during one breath cycle
(5)吸気遅延時間
吸気遅延時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、肺野全体の濃度変化又は横隔膜位置の変化を解析することにより安静呼気位のフレーム画像を抽出し、各小領域毎に、安静呼気位のフレーム画像から、吸気時において安静呼気位の信号値との差が所定の閾値以上となるまでの時間の算出
(5) Intake delay time The intake delay time is calculated by performing the following processing on a series of frame images.
Binning processing → Low-pass filter processing in the time axis direction → Small regions at the same pixel position (regions of pixel blocks output from detection elements at the same position in the FPD) in a series of frame images are associated with each other, and density changes in the entire lung field Alternatively, a frame image of a resting expiratory position is extracted by analyzing a change in the diaphragm position, and a difference between the signal value of the resting expiratory position at the time of inspiration from a frame image of the resting expiratory position for each small region is a predetermined threshold value. Calculation of time until
(6)吸気時間、呼気時間
吸気時間、呼気時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のローパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、呼吸1サイクル中の最大信号値(極大値)、最小信号値(極小値)の算出→最大信号値から最小信号値までの時間を呼気時間、最小信号値から最大信号値までの時間を呼気時間として算出
(6) Inspiration time and expiration time The inspiration time and expiration time are calculated by performing the following processing on a series of frame images.
Binning processing → Low-pass filter processing in the time axis direction → Small regions at the same pixel position (regions of pixel blocks output from detection elements at the same position in the FPD) of a series of frame images are associated with each other, Calculation of maximum signal value (maximum value) and minimum signal value (minimum value) during one breath cycle → time from maximum signal value to minimum signal value is expiration time, time from minimum signal value to maximum signal value is expiration time Calculated as
肺野内の血流による局所的な信号変化を示す特徴量を算出する解析としては、例えば、下記の(7)〜(10)が挙げられる。以下、各特徴量毎に、その算出方法の手順について簡単に説明する。なお、下記においては、撮影されたフレーム画像の生データ(RAWデータ)から特徴量を算出する処理内容について説明するが、本実施の形態においては、解析用WS8へのフレーム画像送信時のデータ量の削減及び処理時間の短縮化の観点から、すでに各画像にビニング処理が施されている(即ち、所定サイズの小領域に分割され、小領域毎に信号値の平均化が行われている)。 For example, the following (7) to (10) may be mentioned as the analysis for calculating the feature amount indicating the local signal change due to the blood flow in the lung field. Hereinafter, the procedure of the calculation method for each feature amount will be briefly described. In the following, processing contents for calculating a feature amount from raw data (RAW data) of a captured frame image will be described. In this embodiment, the data amount at the time of transmitting a frame image to the analysis WS 8 is described. Each image has already been subjected to binning processing (that is, divided into small areas of a predetermined size, and signal values are averaged for each small area) from the viewpoint of reduction of image quality and processing time. .
(7)血流−フレーム間差分画像
フレーム間差分画像は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→フレーム間差分処理→ノイズ除去
時間軸方向のハイパスフィルタ処理は、血流による信号値の時間変化を抽出するための処理であり、例えば、カットオフ周波数0.7Hzでフィルタリングする。その他は、上述の(1)換気−フレーム間差分画像で説明した処理と同様である。
(7) Blood flow-interframe difference image The interframe difference image is calculated by performing the following processing on a series of frame images.
Binning process → High-pass filter process in the time axis direction → Difference process between frames → Noise removal The high-pass filter process in the time axis direction is a process for extracting a temporal change in the signal value due to blood flow. Filter at 7 Hz. Others are the same as the processing described in the above (1) Ventilation-frame difference image.
(8)血流−波形描画
血流−波形描画は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→各小領域毎の波形描画
ビニング処理及び時間軸方向のハイパスフィルタ処理は上述のとおりである(以下同じ)。
波形描画処理は、上述の(2)換気−波形描画で説明した処理と同様の処理である。
(8) Blood flow-waveform drawing Blood flow-waveform drawing is calculated by performing the following processing on a series of frame images.
Binning process → High-pass filter process in the time axis direction → Waveform drawing for each small region The binning process and the high-pass filter process in the time axis direction are as described above (the same applies hereinafter).
The waveform drawing process is the same as the process described in (2) Ventilation-waveform drawing described above.
(9)血流量の振幅
血流量の振幅は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、各小領域毎に、心拍1サイクル中の最大信号値(極大値)−最小信号値(極小値)の算出
(9) Blood Flow Amplitude The blood flow amplitude is calculated by performing the following processing on a series of frame images.
Binning process → High-pass filter process in time axis direction → Small area at the same pixel position (area of the pixel block output from the detection element at the same position of the FPD) of a series of frame images is associated with each other, Calculation of maximum signal value (maximum value)-minimum signal value (minimum value) during one heartbeat cycle
(10)心室収縮遅延時間
吸気遅延時間は、一連のフレーム画像に以下の処理を施すことによって算出される。
ビニング処理→時間軸方向のハイパスフィルタ処理→一連のフレーム画像の同じ画素位置の小領域(FPDの同一位置の検出素子から出力された画素ブロックの領域)を互いに対応付け、心室領域の濃度変化又は心壁位置の変化を解析することにより心室拡張期の終わりに相当するフレーム画像を抽出し、各小領域毎に、心室拡張期の終わりに相当するフレーム画像から、心室収縮期において心室拡張期の終わりの信号値との差が所定の閾値以上となるまでの時間の算出
(10) Ventricular systolic delay time The inspiratory delay time is calculated by performing the following processing on a series of frame images.
Binning processing → High-pass filter processing in the time axis direction → Small regions at the same pixel position (regions of pixel blocks output from detection elements at the same position in the FPD) of a series of frame images are associated with each other, A frame image corresponding to the end of the ventricular diastole is extracted by analyzing the change in the ventricular wall position, and for each subregion, the frame image corresponding to the end of the ventricular diastole is extracted from the frame image corresponding to the end of the ventricular diastole. Calculating the time until the difference from the end signal value exceeds a predetermined threshold
肺野全体の動きを示す特徴量を算出する解析としては、例えば、下記の(11)〜(15)が挙げられる。以下、各特徴量毎に、その算出方法の手順について簡単に説明する。なお、本実施の形態においては、解析用WS8へのフレーム画像送信時のデータ量の削減及び処理時間の短縮化の観点から、各フレーム画像にビニング処理が施されているが、ビニング処理は行っても行わなくても下記特徴量は算出可能である。 For example, the following (11) to (15) may be mentioned as the analysis for calculating the feature amount indicating the movement of the entire lung field. Hereinafter, the procedure of the calculation method for each feature amount will be briefly described. In this embodiment, binning processing is performed on each frame image from the viewpoint of reducing the amount of data when transmitting the frame image to the analysis WS 8 and shortening the processing time. However, the binning processing is performed. Even if it is not performed, the following feature amount can be calculated.
(11)横隔膜移動量解析
横隔膜移動量は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
画像解析により各フレーム画像から横隔膜の位置を抽出→各フレーム画像の横隔膜の位置を追跡し、移動量を算出
(12)胸郭移動量解析
胸郭移動量は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
画像解析により各フレーム画像から上部胸郭(上肋骨(第2〜第6肋骨))、下部胸郭(下肋骨(第7〜第10肋骨))の位置を抽出→各フレーム画像の上胸郭、下胸郭の位置を追跡し、移動量を算出する。
(13)呼吸数、呼吸周期
呼吸数、呼吸周期は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
各フレーム画像から画像解析により求めた横隔膜位置(肺尖から横隔膜までの距離)の変化若しくはローパスフィルタ処理後の肺野全体の信号変化(信号値(平均信号値)の極大値から次の極小値までの時間間隔)より呼吸周期を求め、呼吸周期の逆数から単位時間当たりの呼吸数を算出する。
(14)心拍数、心周期
心拍数、心周期は、撮影された各フレーム画像に、以下の処理を施すことによって算出される。
各フレーム画像から画像解析により求めた心壁位置の変化若しくはハイパスフィルタ処理後の肺野全体の信号変化(信号値(平均信号値)の極大値から次の極小値までの時間間隔)より心周期を求め、心周期の逆数から単位時間当たりの呼吸数を算出する。
(15)スパイロ検査相当の値の算出
横隔膜位置の時間変化波形を算出し、FEV1.0%(一秒率)相当の値を算出する。
胸郭と横隔膜位置の変化から肺野面積の変化量を算出し、別途測定した胸厚の変化量を乗ずることで、VC(肺活量)相当の値を算出する。
(11) Diaphragm movement amount analysis The diaphragm movement amount is calculated by performing the following processing on each captured frame image.
Extract the position of the diaphragm from each frame image by image analysis → track the position of the diaphragm in each frame image and calculate the movement amount (12) Thoracic movement amount analysis The thoracic movement amount is Calculated by performing processing.
Extract the positions of the upper rib cage (upper ribs (second to sixth ribs)) and lower rib cage (lower ribs (seventh to tenth ribs)) from each frame image by image analysis → upper and lower rib cages of each frame image Is tracked and the amount of movement is calculated.
(13) Respiration rate and respiration cycle The respiration rate and respiration cycle are calculated by performing the following processing on each captured frame image.
Changes in the diaphragm position (distance from the apex to the diaphragm) determined by image analysis from each frame image, or the signal change (signal value (average signal value)) from the local maximum to the next local minimum after low-pass filtering The respiratory cycle is obtained from the time interval until the respiratory rate is obtained, and the respiratory rate per unit time is calculated from the reciprocal of the respiratory cycle.
(14) Heart rate and cardiac cycle The heart rate and cardiac cycle are calculated by performing the following processing on each captured frame image.
Cardiac cycle based on changes in heart wall position obtained by image analysis from each frame image or signal changes in the entire lung field after high-pass filter processing (time interval from the maximum value of the signal value (average signal value) to the next minimum value) And the respiration rate per unit time is calculated from the reciprocal of the cardiac cycle.
(15) Calculation of a value corresponding to a spiro test A time change waveform of the diaphragm position is calculated, and a value corresponding to FEV 1.0% (one second rate) is calculated.
A change in lung field area is calculated from changes in the rib cage and diaphragm positions, and a value corresponding to VC (pulmonary vital capacity) is calculated by multiplying the change in chest thickness measured separately.
解析処理が終了すると、解析用WS8において、解析結果のデータが通信ネットワークNを介してコンソール5に送信される(ステップS20)。コンソール5においては、ネットワーク通信部56により解析結果のデータが受信されると、受信された解析結果のデータが撮影オーダ情報と対応付けられてPACS10に送信される(ステップS21)。そして、撮影・解析処理は終了する。なお、特徴量の解析に使用した各フレーム画像データに関しては、静止画のような濃度諧調ベースの病変部読影には使用できないものであるので、保存データ容量削減の観点から、保存せず、算出された特徴量データのみを保存することが好ましい。
また、解析用WS8にて、解析結果のデータと、コンソール5か受信した撮影オーダ情報を対応付け、これらのデータを解析用WS8からPACS10に送信する構成としても良い。
When the analysis process ends, analysis result data is transmitted to the console 5 via the communication network N in the analysis WS 8 (step S20). In the console 5, when the analysis result data is received by the network communication unit 56, the received analysis result data is transmitted to the PACS 10 in association with the imaging order information (step S21). Then, the photographing / analysis process ends. Note that each frame image data used for feature amount analysis cannot be used for density-gradation-based lesion interpretation like a still image, so it is calculated without saving from the viewpoint of saving data capacity. It is preferable to store only the feature amount data.
Alternatively, the analysis WS 8 may associate the analysis result data with the imaging order information received from the console 5, and transmit these data from the analysis WS 8 to the PACS 10.
PACS10においては、受信された解析結果データが撮影オーダ情報に対応付けてサ
ーバ装置のHDD等に記憶される。そして、読影用端末からの要求に応じて読影用端末表示部に解析結果が表示される。
In the PACS 10, the received analysis result data is stored in the HDD or the like of the server device in association with the imaging order information. The analysis result is displayed on the interpretation terminal display unit in response to a request from the interpretation terminal.
ここで、従来、動態撮影で得られた一連のフレーム画像に基づき当該動態に関する特徴量を算出し、診断支援情報として提供するシステムにおいては、診断精度の向上のため、複数のフレーム画像において被写体の同一部分が描画された領域を互いに対応付ける、所謂ワーピング処理を行うことが必要であると考えられていた(例えば、特許文献1、2)。
このワーピング処理を行うには、一のフレーム画像を複数の小領域に分割し、当該一のフレーム画像の各小領域に描画された構造物の部分と同一の部分を描画した小領域を各フレーム画像毎に抽出しなければならない。ワーピング処理では、一般的には肺野内構造物による空間的な濃度変化をもとに位置合わせを行うため、各フレーム画像にわたり当該構造物の濃度が忠実に(均一に)再現される必要があり、従って、検出器の各画素の出力バラツキ等は極力抑える必要があり(従い、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理、ラグ補正処理等の各種補正処理を行って、バラツキを補正する必要がある)、補正処理に時間を要し、更に、高精度のワーピング処理を行うには、それに従い分解能の細かい画像を必要とするので画素サイズの小さな検出器が必要となり、各フレーム画像毎のデータ容量が増え、処理対象となる全体データ容量が大幅増加となる。そのため、大容量メモリや高速処理CPU等のハードウエアが必要であり、かつ処理時間も必要となっていた。
Here, conventionally, in a system that calculates feature values related to dynamics based on a series of frame images obtained by dynamic imaging and provides them as diagnosis support information, in order to improve diagnosis accuracy, It has been thought that it is necessary to perform a so-called warping process in which areas where the same part is drawn are associated with each other (for example, Patent Documents 1 and 2).
In order to perform this warping process, one frame image is divided into a plurality of small areas, and a small area in which the same part of the structure drawn in each small area of the one frame image is rendered as each frame. It must be extracted for each image. In the warping process, alignment is generally performed based on spatial density changes caused by structures in the lung field, so the density of the structure must be faithfully (uniformly) reproduced over each frame image. Therefore, it is necessary to suppress the output variation of each pixel of the detector as much as possible (accordingly, various correction processes such as an offset correction process, a gain correction process, a defective pixel correction process, and a lag correction process are performed to correct the variation. Correction process takes time, and in order to perform high-accuracy warping, an image with a fine resolution is required accordingly, so a detector with a small pixel size is required. The total data capacity to be processed is greatly increased. For this reason, hardware such as a large-capacity memory and a high-speed processing CPU is required, and processing time is also required.
しかし、本願発明者等が鋭意検討を行った結果、ワーピング処理を行わずとも、FPD個々の検出素子単位、又は複数画素をひとまとまりとした画素ブロック単位に動態撮影にかかる一連のフレーム画像間の差異を比較することで、同等の解析結果を得られることを見出した。 However, as a result of intensive studies by the inventors of the present application, even without performing warping processing, between FPD individual detection element units or a series of frame images related to dynamic imaging in a pixel block unit including a plurality of pixels as a group. It was found that equivalent analysis results can be obtained by comparing the differences.
以下、肺野を例にとり、ワーピング処理を行なわずとも同等の解析結果が得られる理由について説明する。
まず、体厚方向(z方向;側面)における信号値の変化ついて説明する。
図8は、安静呼気位である時刻T1における肺野の体厚方向(z方向)を模式的に示す図、時刻T1から吸気していくことにより安静吸気位となった時刻T2における肺野の体厚方向を模式的に示す図、及びFPDの検出素子位置(体軸方向(y方向))を模式的に示す図である。図8においては、吸気することにより肺胞a及び肺胞bのy方向の位置は下方へ移動し、時刻T1における肺胞bのy方向の位置と時刻T2における肺胞aのy方向の位置が一致した例を示している。
Hereinafter, the reason why an equivalent analysis result can be obtained without performing the warping process will be described using the lung field as an example.
First, changes in the signal value in the body thickness direction (z direction; side surface) will be described.
FIG. 8 is a diagram schematically showing the body thickness direction (z direction) of the lung field at time T1 in the resting expiratory position, and the lung field at time T2 at which the inspiratory position is reached by inhaling from time T1. It is a figure which shows a body thickness direction typically, and a figure which shows the detection element position (body-axis direction (y direction)) of FPD typically. In FIG. 8, the positions of alveoli a and b in the y direction move downward by inhaling, and the position of alveoli b in the y direction at time T1 and the position of alveolar a in the y direction at time T2. Shows an example of matching.
吸気によって肺野内における肺胞の位置は移動する。そこで、フレーム画像間で同一の肺胞に対して位置合わせを行い、ワーピングしてから信号値の差分をとると、肺野外の部分でのz方向のX線減衰量は肺野のy方向の位置によって異なるため、肺胞の位置合わせを行うことにより、逆に、呼吸による肺胞の密度変化による信号増加分に対し、肺野外部分でのX線減衰量の差が誤差因子として加わってしまうことになる。 The position of the alveoli in the lung field is moved by inspiration. Therefore, when the same alveoli between the frame images is aligned and the difference between the signal values is taken after warping, the X-ray attenuation amount in the z-direction in the portion outside the lung field is in the y-direction of the lung field. As the alveoli are aligned, the difference in X-ray attenuation in the lung field is added as an error factor to the signal increase due to the change in alveolar density due to respiration. It will be.
例えば、図8において、時刻T1におけるフレーム画像とT2におけるフレーム画像とで、肺胞bを位置合わせ及びワーピングをしてから信号値の差分をとると、この差分値の中には、時刻T1の実線の矢印で示す肺野外でのX線減衰量と、時刻T2の点線矢印で示す肺野外でのX線減衰量との差分も含んでいることとなり、この肺野外でのX線減衰量の差分が、同一肺胞間での呼吸による密度変化による信号値の変化に、誤差として加わることとなる。これにより、肺胞の密度変化による信号変化量の計算精度が下がる。 For example, in FIG. 8, when the difference between the signal values is obtained after aligning and warping the alveoli b between the frame image at the time T1 and the frame image at the time T2, the difference value includes the difference at the time T1. This also includes the difference between the X-ray attenuation outside the lung field indicated by the solid arrow and the X-ray attenuation outside the lung field indicated by the dotted arrow at time T2, and the X-ray attenuation outside the lung field is also included. The difference is added as an error to the change in the signal value due to the density change due to respiration between the same alveoli. Thereby, the calculation accuracy of the signal change amount due to the change in alveolar density is lowered.
ここで、肺胞の位置合わせ及びワーピングを行わず、時刻T1におけるフレーム画像の肺胞bと、時刻T2におけるフレーム画像の肺胞aを描画した画素(画素ブロック)、つまり、FPDの同一位置の検出素子(検出素子群)から出力された信号値の差分を計算することとする。このとき、この画素(画素ブロック)に描画されている肺胞は異なるが、肺野の同一のy方向位置に対して差分値を計算しているため、図8に示すように、肺野外におけるX線減衰量は変化しない。従って、異なる肺胞間で信号値の差分値を計算した場合、同一肺胞間での呼吸による密度変化の差による信号変化に対して、異なる肺胞間での密度の差による信号変化(図8の同一タイミングにおけるaとbの密度の差)が誤差として加わることとなる。 Here, without performing alignment and warping of the alveoli, pixels (pixel blocks) in which the alveoli b of the frame image at time T1 and the alveoli a of the frame image at time T2 are drawn, that is, at the same position of the FPD. The difference between the signal values output from the detection elements (detection element group) is calculated. At this time, alveoli drawn in this pixel (pixel block) are different, but since the difference value is calculated for the same y-direction position of the lung field, as shown in FIG. The amount of X-ray attenuation does not change. Therefore, when the difference value of the signal value between different alveoli is calculated, the signal change due to the difference in density between different alveoli (Fig. The difference between the density of a and b at the same timing of 8) is added as an error.
「肺野内のy方向位置が異なる肺胞間での密度の違いによる信号変化」は「y方向位置が異なる肺野外におけるX線減衰量の違いによる信号変化」と同等又はそれ以下である。そのため、フレーム画像間での肺胞の位置合わせ及びワーピング処理は行わず、そのままFPDの画素単位で差分をとる方が、処理の手間が省け、同等レベルの誤差で肺胞の密度変化による信号変化量が算出できる。 “Signal change due to difference in density between alveoli having different y-direction positions in the lung field” is equal to or less than “signal change due to difference in X-ray attenuation in the outside of the lung field having different y-direction positions”. For this reason, the alveolar alignment and the warping process between frame images are not performed, and it is possible to save the processing effort by taking the difference for each pixel of the FPD as it is, and the signal change due to the alveolar density change with the same level of error. The amount can be calculated.
特に、個々の画素や小領域毎に内在する誤差成分は、肺野全体の換気量情報を算出する際に加算されて相殺されることになり、肺野全体の換気や血流に関する特徴量を算出する場合には、ワーピング処理を行うと、ワーピング処理による処理時間延長というマイナス効果しか得られないことになる。 In particular, the error component inherent in each pixel or small area is added and canceled when calculating the ventilation information of the entire lung field, and the characteristic quantity related to ventilation and blood flow of the entire lung field is calculated. In the case of calculation, if the warping process is performed, only a negative effect of extending the processing time by the warping process can be obtained.
次に、x−y方向について検討する。図9に、肺野を正面(x−y方向)から見た図を示す。図9における実線は、図8における時刻T1のフレーム画像をx−y方向(正面)から見た図を示し、図9における点線は、図8における時刻T2のフレーム画像をx−y方向から見た図を示している。
図9に示すように、通常、吸気においては、肺胞は左肺野の場合左下方向に、右肺野の場合右下方向に移動する。これを鉛直方向(y方向)への移動と水平方向(x方向)への移動に分解する。y方向への肺胞の移動に対するワーピング処理は前述のとおりである。以下に、安静換気時の、肺胞のx方向への移動について説明する。
Next, the xy direction will be considered. FIG. 9 shows a view of the lung field viewed from the front (xy direction). The solid line in FIG. 9 shows a view of the frame image at time T1 in FIG. 8 viewed from the xy direction (front), and the dotted line in FIG. 9 shows the frame image at time T2 in FIG. The figure is shown.
As shown in FIG. 9, normally, in inspiration, the alveoli move in the lower left direction in the case of the left lung field and in the lower right direction in the case of the right lung field. This is broken down into movement in the vertical direction (y direction) and movement in the horizontal direction (x direction). The warping process for the movement of the alveoli in the y direction is as described above. Hereinafter, the movement of the alveoli in the x direction during rest ventilation will be described.
安静換気時の場合、胸郭の変化幅は最大でも10mm程度である。このとき、安静呼気位から安静吸気位に至る変化を考えた場合、x方向の移動量は、胸郭すぐ内側に位置する肺胞が最も大きく、5mm程度である。3.75枚/秒のフレームレートで動態画像を取得し、フレーム間差分値を算出する場合を考えると、隣接フレーム画像間での肺胞の移動量は更に小さくなり、無視しうる。安静換気時においてはx方向の移動量は小さく、ワーピング処理を施したときと施さないときの信号値の変化量は、ほぼ同等である。従って、x−Y方向においてもワーピング処理は不要である。 In the case of rest ventilation, the change width of the rib cage is about 10 mm at the maximum. At this time, when considering the change from the resting expiratory position to the resting inspiratory position, the amount of movement in the x direction is the largest in the alveoli located just inside the rib cage and is about 5 mm. Considering the case where dynamic images are acquired at a frame rate of 3.75 frames / second and the difference value between frames is calculated, the amount of movement of the alveoli between adjacent frame images is further reduced and can be ignored. During rest ventilation, the amount of movement in the x direction is small, and the amount of change in signal value when the warping process is performed and when the warping process is not performed is substantially the same. Accordingly, no warping process is required in the x-Y direction.
以上の知見によれば、ワーピング処理を施さず、FPDの個々の検出素子の出力を示す画素単位、或いは小領域単位の演算のみを行えばよいので、解析用WS8は、FPDからの出力信号形態に対する制約を行わず、各社の動態撮影対応のFPD9aで生成されるフレーム画像、例えば、ビニング処理を施されたフレーム画像を処理することができ、いわばオープンシステム化を可能とする。 According to the above knowledge, since the warping process is not performed, it is only necessary to perform the calculation in units of pixels indicating the outputs of the individual detection elements of the FPD or small area units. Therefore, it is possible to process a frame image generated by the FPD 9a corresponding to dynamic photographing of each company, for example, a frame image subjected to binning processing, which enables an open system.
ところで、動態撮影は、解析に使用するフレーム画像数が静止画に比べて飛躍的に増加し、処理時間が飛躍的に増大する難点を有している。本願発明者らは、この観点からも鋭意検討を行った。 By the way, dynamic imaging has a problem that the number of frame images used for analysis increases dramatically compared to still images, and the processing time increases dramatically. The inventors of the present application have also made extensive studies from this viewpoint.
従来の胸部や***の異常陰影候補を検出して医師へ診断支援情報として提供するシステム(CAD)では、原画像データへ検出アルゴリズムを適用する方式であるので、アルゴリズムを適用する画像自体の濃度分解能や画素サイズが検出アルゴリズムにマッチすること(画素サイズや濃度分解能が高精細であること)が重要視されていた。
また、従来のシステムにおいて、原画像から間引き画像が作成されていたが、これは被写体の関心領域が解析可能な範囲内に収まっているかの確認や、検出アルゴリズムに適合する濃度範囲に検出対象画像の濃度範囲を収めるための階調処理条件の算出等、検出処理の前準備段階で使用されるだけであり、CADによる検出処理においては使用されることなく、廃棄されるものであった。
In the conventional system (CAD) that detects abnormal shadow candidates of the breast and breast and provides them to the doctor as diagnosis support information, the detection algorithm is applied to the original image data, so the density resolution of the image itself to which the algorithm is applied And that the pixel size matches the detection algorithm (the pixel size and density resolution are high-definition) have been emphasized.
Also, in the conventional system, a thinned image is created from the original image. This is done by confirming that the region of interest of the subject is within the analyzable range, and detecting the target image within the density range that matches the detection algorithm. It is only used in the preparatory stage of detection processing, such as calculation of gradation processing conditions for keeping the density range, and is not used in the detection processing by CAD, but is discarded.
一方、上述のように、動態画像の特徴量解析においては、主として隣接するフレーム間での差分値処理を基本としている。そのため、動態画像の特徴量解析においては、各画素の絶対的出力値そのものを閾値と比較したり、微細な構造物を抽出したりする従来のCADにおける解析処理のようには各画素の厳密な出力値を必要とせず、また個々の画像の画素サイズに影響は受けないことを本願発明者らは見出した。 On the other hand, as described above, the feature value analysis of dynamic images is mainly based on the difference value processing between adjacent frames. Therefore, in the dynamic image feature quantity analysis, the absolute output value of each pixel itself is compared with a threshold value, or a strict analysis of each pixel is performed as in analysis processing in a conventional CAD in which a fine structure is extracted. The present inventors have found that no output value is required and the pixel size of each image is not affected.
また、本願発明者らは、更に鋭意検討を行った結果、画素単位の演算(解析)から、特定の画素ブロックサイズ(ブロックサイズ)の小領域単位の演算(各小領域内の一の画素値を使用した演算、小領域内の画素値の代表値(平均値等)を使用した演算)に変換しても同様の結果が得られることを見出した。そして、これにより、解析処理の演算に要するデータ量を大幅に削減し、処理時間を大幅に短縮可能であることを見出した。 In addition, as a result of further intensive studies, the inventors of the present application have calculated from a calculation (analysis) in units of pixels to a calculation in units of small areas of a specific pixel block size (block size) (one pixel value in each small area) It has been found that the same result can be obtained even if the calculation is performed by using the above-described method, and the calculation using a representative value (average value or the like) of pixel values in a small region. As a result, the present inventors have found that the amount of data required for the calculation of the analysis processing can be greatly reduced, and the processing time can be greatly shortened.
図10は、小領域の画素ブロックサイズを0.5mm角〜10mm角の間で変化させたときの上記解析処理の代表的な項目の解析結果及び処理時間の評価結果を示す図である。
ここでは、換気量のフレーム間差分画像、最大流速比のヒストグラム解析、血流量のフレーム間差分画像を対象項目とした場合を例として示している。なお、各フレーム画像は、照射線量及びフレームレート等の放射線照射条件及び画像読取条件を一定として撮影された画像を用いている。図10では、フレームレート7.5枚/秒、10秒間での入射表面線量は0.2mGy、FPD9aの画素サイズは200μmである。なお、図10(及び後述する図11A及び図11B)において、評価結果△以上は診断に使用できる精度で特徴量の算出が可能であることを示し、△→○→◎の順に解析精度が低い→高いとなっていることを示す。また、処理時間×は処理時間が長く、あまり実用化には向かないことを示し、△以上は処理時間が実用上許容できる範囲であることを示す。なお、この処理時間の場合、Xレベルであっても、全画素データを使う場合にくらべると、かなりの処理時間短縮は可能にはなっているので、使用不可というわけでは無い。
FIG. 10 is a diagram illustrating analysis results of typical items of the above analysis processing and evaluation results of processing time when the pixel block size of the small region is changed between 0.5 mm square to 10 mm square.
Here, an example is shown in which a ventilation amount inter-frame difference image, a maximum flow rate ratio histogram analysis, and a blood flow amount inter-frame difference image are set as target items. Note that each frame image uses an image taken with radiation irradiation conditions such as an irradiation dose and a frame rate, and image reading conditions constant. In FIG. 10, the incident surface dose at a frame rate of 7.5 frames / second and 10 seconds is 0.2 mGy, and the pixel size of the FPD 9a is 200 μm. In FIG. 10 (and FIGS. 11A and 11B described later), the evaluation result Δ or higher indicates that the feature quantity can be calculated with accuracy that can be used for diagnosis, and the analysis accuracy is low in the order of Δ → ○ → ◎. → Indicates that the price is high. Further, the processing time x indicates that the processing time is long and is not suitable for practical use, and Δ or more indicates that the processing time is within a practically allowable range. In the case of this processing time, even at the X level, the processing time can be considerably reduced as compared with the case where all pixel data is used.
図10に示すように、ブロックサイズ0.5mm角〜1mm角である場合、換気量のフレーム間差分値算出、最大流速比のヒストグラム解析、血流量のフレーム間差分値の全ての項目において解析結果は高い精度で得ることができた。ブロックサイズ2mm角では、ブロックサイズ0.5mm角〜1mm角より若干解析精度が落ちるものの、診断には十分な結果が得られた。ブロックサイズ5mm角では、各項目について解析精度に低下が見られた。ただし、診断に使用可能な解析精度は確保されている。ブロックサイズ10mm角である場合、診断に使用可能な精度の特徴量を得ることができなかった。一方、処理時間は、図10に示すように、ブロックサイズが小さいほど処理時間は長くなった。 As shown in FIG. 10, when the block size is 0.5 mm square to 1 mm square, the analysis results are obtained for all items of calculation of the difference value of the ventilation volume between frames, histogram analysis of the maximum flow velocity ratio, and difference value of the blood flow volume between frames. Could be obtained with high accuracy. With a block size of 2 mm square, the analysis accuracy was slightly lower than the block size of 0.5 mm square to 1 mm square, but sufficient results were obtained for diagnosis. When the block size was 5 mm square, the analysis accuracy decreased for each item. However, analysis accuracy that can be used for diagnosis is ensured. When the block size was 10 mm square, it was not possible to obtain feature quantities with accuracy that could be used for diagnosis. On the other hand, as shown in FIG. 10, the processing time is longer as the block size is smaller.
ここで解析結果の評価について一例を挙げて説明する。
図11A〜図11Cに、ブロックサイズをそれぞれ2mm角、5mm角、10mm角としたときの最大流速比のヒストグラム解析の解析結果の一例を示す。
なお、ブロックサイズを0.5mm角、1mm角とした場合は、図11Aに示すヒストグラムとほとんど変化はみられなかった。
図11Aに示すブロックサイズ2mm角の解析結果と図11Bに示すブロックサイズ5mm角の解析結果を比較すると、若干ヒストグラムの形状に変化が見られる。図11Aに示すブロックサイズ2mm角の解析結果と図11Cに示すブロックサイズ10mm角の解析結果を比較すると、ヒストグラムの形状に大幅な変化が見られ、解析精度が大幅に低下することがわかる。よって、解析精度の点からいえば、ブロックサイズ5mm角以下、特に、2mm角以下であることが必要である。一方、処理時間の点からすると、上述の図10に示すように、ブロックサイズが小さいほど処理時間は長くなる。
よって、解析精度の維持及び処理時間短縮化の観点から、肺野の換気の解析ではブロックサイズは2mm角から5mm角程度とすることが好ましい。
Here, an example of the analysis result evaluation will be described.
11A to 11C show examples of analysis results of histogram analysis of the maximum flow rate ratio when the block sizes are 2 mm square, 5 mm square, and 10 mm square, respectively.
When the block size was 0.5 mm square and 1 mm square, there was almost no change from the histogram shown in FIG. 11A.
When the analysis result of the block size of 2 mm square shown in FIG. 11A is compared with the analysis result of the block size of 5 mm square shown in FIG. 11B, the shape of the histogram is slightly changed. Comparing the analysis result of the block size of 2 mm square shown in FIG. 11A and the analysis result of the block size of 10 mm square shown in FIG. 11C, it can be seen that a significant change is seen in the shape of the histogram and the analysis accuracy is greatly reduced. Therefore, in terms of analysis accuracy, the block size needs to be 5 mm square or less, particularly 2 mm square or less. On the other hand, from the viewpoint of processing time, as shown in FIG. 10 described above, the processing time increases as the block size decreases.
Therefore, from the viewpoint of maintaining the analysis accuracy and shortening the processing time, the block size is preferably about 2 mm square to 5 mm square in the analysis of lung field ventilation.
更に、肋骨の移動による信号変化の除去の観点からも、ブロックサイズは2mm程度が好ましい。例えば、肋骨の移動による信号変化の影響を除去すべく、複数の肋骨を含むようにブロックサイズを広げると、肋骨の太さが10mm〜20mmなので、複数の肋骨が含まれるブロックサイズは50mm角〜100mm角である。これでは分解能が低すぎて、局所的な解析には不向きである。肋骨一本+安静呼吸時の肋骨移動量が含まれるようにブロックサイズを広げると、肋骨の太さが上述のように10mm〜20mm、安静呼吸時の肋骨の移動量は数mmであるので、ブロックサイズは15mm角〜20mm角となる。これでも分解能が低すぎて、局所的な解析には不向きである。また、他の肋骨が含まれないように小領域を区切るのは困難である。そこで、ある程度ブロックサイズを小さくして肋骨による信号変化が周囲より大きい領域を検出し、その部分を解析から省く、又は、周期の信号値により置き換える等を行うことで、肋骨による影響を除去する。この場合のブロックサイズは2mm角〜5mm角程度が好ましい。
以上の観点からも、肺野の換気に関する解析処理を行う場合のブロックサイズは2mm角〜5mm角程度が好ましい。
Furthermore, the block size is preferably about 2 mm from the viewpoint of removing signal changes due to movement of the ribs. For example, when the block size is expanded to include a plurality of ribs in order to remove the influence of signal changes due to the movement of the ribs, the thickness of the ribs is 10 mm to 20 mm, so the block size including the plurality of ribs is 50 mm square to 100 mm square. This is too low resolution and is not suitable for local analysis. When the block size is expanded to include one rib + the amount of rib movement at rest breathing, the rib thickness is 10 mm to 20 mm as described above, and the amount of rib movement at rest breathing is several mm, The block size is 15 mm square to 20 mm square. Even so, the resolution is too low to be suitable for local analysis. In addition, it is difficult to divide the small area so that other ribs are not included. Therefore, the influence of the ribs is removed by detecting a region where the signal change due to the rib is larger than the surroundings by reducing the block size to some extent and omitting that portion from the analysis or replacing it with a signal value of the period. In this case, the block size is preferably about 2 mm square to 5 mm square.
Also from the above viewpoint, the block size in the case of performing the analysis processing related to the ventilation of the lung field is preferably about 2 mm square to 5 mm square.
また、本願発明者らは、FPD9aのフレームレートによる解析への影響についても検討を行った。
図12A〜図12Bは、FPD9aのフレームレートを2枚/秒〜20枚/秒と変化させたときの上記解析処理の代表的な項目の解析結果、処理時間、間引きによる影響、患者被曝の程度、の評価結果を示す図である。図12Aは、各フレーム画像における照射線量を一定にした場合の評価結果を示す。図12Bは、一回の動態撮影のトータルの照射線量を一定(撮影10秒間での入射照射線量0.2mGy相当)にした場合の評価結果を示す。図12A、図12Bとも画素サイズは200μm、ブロックサイズは2mm角である。
In addition, the inventors of the present application also examined the influence on the analysis by the frame rate of the FPD 9a.
12A to 12B show the analysis results of typical items of the above analysis processing, the processing time, the influence of thinning, and the degree of patient exposure when the frame rate of the FPD 9a is changed from 2 frames / second to 20 frames / second. It is a figure which shows the evaluation result of. FIG. 12A shows an evaluation result when the irradiation dose in each frame image is constant. FIG. 12B shows an evaluation result when the total irradiation dose of one dynamic imaging is constant (equivalent to an incident irradiation dose of 0.2 mGy in 10 seconds of imaging). 12A and 12B, the pixel size is 200 μm, and the block size is 2 mm square.
各フレーム画像の照射線量を一定とした場合、図12Aに示すように、換気量、最大流速比のヒストグラム解析については3.75枚/秒以上のフレームレートであれば診断に利用できる精度の特徴量を得ることができた。なお、換気遅延時間については診断に利用できる精度の特徴量を得るには10枚/秒以上、血流遅延時間については30フレーム以上が必要であった。一方、処理時間についてはフレームレートが高くなるほど長くなり、15枚/秒以上は実用に耐えなかった。なお、画素間隔を1/8に間引いた場合は、処理時間は全てのフレームレートで良好であった。患者の被曝線量については、フレームレートが高くなるほど増加し、30枚/秒以上となると、静止画の胸部X線画像の撮影時の通常の入射表面線量の2倍以上となり、許容できる範囲を超えた。 Assuming that the irradiation dose of each frame image is constant, as shown in FIG. 12A, for the histogram analysis of the ventilation volume and the maximum flow rate ratio, a feature of accuracy that can be used for diagnosis if the frame rate is 3.75 frames / second or more. I was able to get the amount. In addition, about the ventilation delay time, in order to obtain the feature quantity of the precision which can be used for a diagnosis, 10 sheets / second or more and the blood flow delay time required 30 frames or more. On the other hand, the processing time was longer as the frame rate was higher, and the processing time of 15 sheets / second or more was not practical. When the pixel interval was thinned down to 1/8, the processing time was good at all frame rates. The patient exposure dose increases as the frame rate increases, and when it exceeds 30 frames / second, it is more than twice the normal incident surface dose at the time of taking a chest X-ray image of a still image, exceeding the allowable range. It was.
また、動態撮影においてトータルの照射線量を一定とし、フレームレートが増えるにつれて1フレーム画像あたりの照射線量が少なくなるようにした場合、図12Bに示すように、各解析項目において、フレームレートが15枚/秒以上では図12Aの評価結果に比べて低い評価となっている。これは、フレームレートが高くなるにつれて各フレーム画像のS/N比が悪くなり、画像が劣化したためである。 In addition, when the total irradiation dose is constant in dynamic imaging and the irradiation dose per frame image decreases as the frame rate increases, as shown in FIG. / Second or more, the evaluation is lower than the evaluation result of FIG. 12A. This is because as the frame rate increases, the S / N ratio of each frame image deteriorates and the image deteriorates.
ここで、図12Bの解析結果の評価について一例を挙げて説明する。
図13A〜図13Eに、動態撮影においてトータルの照射線量を一定とし、フレームレートを2枚/秒〜30枚/秒と変化させたときの最大流速比のヒストグラム解析の解析結果の一例を示す。
図13Aに示すように、フレームレート2枚/秒とした場合、フレームレート3.75枚/秒とした場合に比べ、平均値は1.15→1.32(15%増)、分散値は0.38→0.52(37%増)と解析結果に差異が生じている。また、異常値となる小領域が増えており(図13A〜図13Eの肺野領域の中で静止画上に最大流速比に応じた輝度が重畳されていない領域)、解析精度は低下している。ここで、吸気、呼気のそれぞれにおいてフレーム間差分値の最大値を算出する際、ノイズの影響を抑制する目的で、フレーム間差分値の最大値が所定閾値より小さい場合は異常値としてその小領域は解析対象から除外している。より具体的には、吸気時は信号値が増加するため、フレーム間差分値の最大値が正の所定閾値以上である場合を正常値、呼気時は信号値が減少するため、フレーム間差分値の最大値が負の所定閾値以下である場合を正常値、それ以外を異常値と判定している。
Here, an example of the evaluation of the analysis result in FIG. 12B will be described.
FIGS. 13A to 13E show examples of analysis results of histogram analysis of the maximum flow rate ratio when the total irradiation dose is constant in dynamic imaging and the frame rate is changed from 2 frames / second to 30 frames / second.
As shown in FIG. 13A, when the frame rate is 2 frames / second, the average value is 1.15 → 1.32 (15% increase) and the dispersion value is compared with the case where the frame rate is 3.75 frames / second. There is a difference in analysis results from 0.38 to 0.52 (37% increase). In addition, small areas that become abnormal values are increasing (areas in which the luminance corresponding to the maximum flow rate ratio is not superimposed on the still image in the lung field areas of FIGS. 13A to 13E), and the analysis accuracy decreases. Yes. Here, when calculating the maximum value of the inter-frame difference value in each of inspiration and expiration, if the maximum value of the inter-frame difference value is smaller than a predetermined threshold for the purpose of suppressing the influence of noise, the small area is set as an abnormal value. Are excluded from the analysis. More specifically, since the signal value increases during inspiration, the maximum value of the interframe difference value is equal to or greater than a predetermined positive threshold value, and the signal value decreases when exhalation. When the maximum value is equal to or less than a predetermined negative threshold value, it is determined as a normal value, and other values are determined as abnormal values.
フレームレート3.75枚/秒以上の場合、図13B〜図13Eに示すように、平均的な呼吸数(安静時15回〜20回/分。18回/分とすると呼吸周期は3.33秒。)の成人の解析結果であれば略同じ解析結果が得られる。ただし、レアケースな頻呼吸患者(単位時間当たりの呼吸数が多い患者。24回〜40回/分。40回/分とすると呼吸周期は1.5秒。)を考慮すると、7.5枚/秒以上がより好ましい。高フレームレートでは、1フレーム画像あたりの線量が減ったとしてもフレーム画像数が増えた分、時間軸方向のローパスフィルタに使用するタップ数(積算回数)を増やすことで、ノイズが低減され、解析結果の劣化を抑制することができる。 When the frame rate is 3.75 frames / second or more, as shown in FIG. 13B to FIG. 13E, the average respiratory rate (15 to 20 times / minute at rest. If it is the analysis result of an adult in seconds.), Almost the same analysis result is obtained. However, in consideration of rare cases of tachypnea patients (patients with high respiratory rate per unit time. 24 to 40 times / minute. If 40 times / minute, the respiratory cycle is 1.5 seconds), 7.5 sheets / Second or more is more preferable. At high frame rates, even if the dose per frame image decreases, the number of frame images increases, so by increasing the number of taps (number of integrations) used for the low-pass filter in the time axis direction, noise is reduced and analysis is performed. The deterioration of the result can be suppressed.
しかし、高フレームレートを実現するには、FPD9aのデータ読み出し及び転送時間の高速化、短パルスを出力するためのX線発生装置の複雑化が必要となり、ハードウエアが高コストとなる。そのため、トータル線量を一定とする場合、換気機能の解析において高フレームレートとするメリットはあまりない。
例えば、SID=200cm(被写体厚を20cmと想定)、フレームレート7.5枚/秒における胸部動態撮影の撮影条件の一例としては、管電圧100kV、管電流50mA、パルス幅2ms、付加フィルタAl0.5mm+Cu0.1mmが挙げられる。このとき、1フレーム画像あたりのX線パルスの照射線量は0.1mAs(=50mA×0.002s)であり、トータル線量が一定となるようにフレームレートに応じて1フレーム画像あたりの照射線量を変えた場合、フレームレート15枚/秒における1フレーム画像あたりのX線パルスの照射線量は0.05mAs、フレームレート30枚/秒における1フレーム画像あたりのX線パルスの照射線量は0.025mAsとなる。しかし、X線パルス幅を1msに制御することは難しく、また、管電流を下げようとすると、X線管のアノード・カソード間に存在するコンデンサに蓄えられた電荷を放電する時定数が長くなり、その結果、管電圧の下降が穏やかになるため、上記電荷の放電を急峻に行うための回路が必要となり、装置コストが増大する。従って、X線パルスの制御の観点では、低フレームレートの方が有利となる。
また、一般的に、FPDにおける残像(ラグ)は、時間(撮影間隔)に対して指数関数的に減少するため、トータル線量が一定となるようにフレームレートに応じて1フレーム画像当たりの照射線量を変えた場合でも、即ち、フレームレートに反比例させて線量を減少させた場合でも、残像の視認性の観点では、フレームレートが小さいほうが有利となる。
従って、換気の解析では、フレームレートを3.75枚/秒〜7.5枚/秒とすることが好ましい。
However, in order to realize a high frame rate, it is necessary to increase the data reading and transfer time of the FPD 9a and to make the X-ray generator complicated for outputting a short pulse, which increases the cost of hardware. Therefore, when the total dose is constant, there is not much merit in setting a high frame rate in the analysis of the ventilation function.
For example, as an example of imaging conditions for chest dynamic imaging at SID = 200 cm (assuming the subject thickness is 20 cm) and a frame rate of 7.5 frames / second, tube voltage 100 kV, tube current 50 mA, pulse width 2 ms, additional filter Al0. 5mm + Cu0.1mm is mentioned. At this time, the irradiation dose of the X-ray pulse per frame image is 0.1 mAs (= 50 mA × 0.002 s), and the irradiation dose per frame image is set according to the frame rate so that the total dose is constant. When changed, the irradiation dose of X-ray pulses per frame image at a frame rate of 15 frames / second is 0.05 mAs, and the irradiation dose of X-ray pulses per frame image at a frame rate of 30 frames / second is 0.025 mAs. Become. However, it is difficult to control the X-ray pulse width to 1 ms, and when trying to lower the tube current, the time constant for discharging the charge stored in the capacitor between the anode and cathode of the X-ray tube becomes longer. As a result, the decrease in the tube voltage becomes gentle, so that a circuit for sharply discharging the electric charge is required, which increases the device cost. Therefore, a low frame rate is more advantageous from the viewpoint of X-ray pulse control.
In general, the afterimage (lag) in the FPD decreases exponentially with respect to time (imaging interval), so that the irradiation dose per frame image according to the frame rate so that the total dose is constant. Even when the value is changed, that is, when the dose is decreased in inverse proportion to the frame rate, the smaller frame rate is advantageous from the viewpoint of the afterimage visibility.
Therefore, in the ventilation analysis, the frame rate is preferably 3.75 frames / second to 7.5 frames / second.
図14A〜図14Dに、動態撮影においてトータルの照射線量を一定とし、フレームレートを3.75枚/秒〜30枚/秒と変化させたときのフレーム間差分画像(血液の拍出タイミング)の解析結果の一例を示す。
フレームレート3.75枚/秒とした場合、フレーム画像間隔が離れすぎているため、場合によっては心臓からの拍出タイミングを捉えられない場合がある。フレームレート7.5枚/秒、15枚/秒とした場合、図14B、図14Cに示すように、平均的な脈拍数(安静時50回〜100回/分。60回/分とすると心拍周期は1.0秒。)の成人の画像であれば血液の拍出タイミングを捉えることができる。ただし、レアケースな頻脈患者(単位時間当たりの脈拍数が多い患者。100回〜120回/分。120回/分とすると心拍周期は0.5秒。)の場合、15枚/秒以上が好ましい。フレームレート30枚/秒とした場合は、1フレーム画像当たりの線量が少ないためにノイズが多く、フレーム間差分画像の画質が低下している。
従って、血流の解析では、フレームレート7.5枚/秒〜15枚/秒とすることが好ましい。
14A to 14D show inter-frame difference images (blood pumping timings) when the total irradiation dose is constant in dynamic imaging and the frame rate is changed from 3.75 frames / second to 30 frames / second. An example of an analysis result is shown.
When the frame rate is set to 3.75 frames / second, the frame image interval is too far away, and in some cases, the timing to beat from the heart may not be captured. When the frame rate is 7.5 frames / second and 15 frames / second, as shown in FIGS. 14B and 14C, the average pulse rate (resting 50 to 100 times / minute. If it is an adult image with a period of 1.0 second, the blood ejection timing can be captured. However, in rare case tachycardia patients (patients with high pulse rate per unit time. 100 to 120 times / minute. If 120 times / minute, heart rate cycle is 0.5 seconds), 15 / second or more Is preferred. When the frame rate is 30 frames / second, the dose per frame image is small, so there is a lot of noise and the image quality of the inter-frame difference image is degraded.
Therefore, in the analysis of blood flow, it is preferable to set the frame rate to 7.5 / 15 to 15 / second.
次に、ビニングと単純間引きの解析への影響について検討する。ここでは、ビニング処理及び単純間引き処理を施した場合の最大流速比のヒストグラム解析、換気のフレーム間差分画像、血流フレーム間差分画像を比較して、両者の影響を検討する。 Next, the effect of binning and simple decimation on the analysis will be examined. Here, the histogram analysis of the maximum flow rate ratio when the binning process and the simple decimation process are performed, the difference image between ventilation frames, and the difference image between blood flow frames are compared, and the effects of both are examined.
図15A〜図15Cに、ビニング処理及び単純間引き処理を施した場合の最大流速比のヒストグラム解析結果、換気フレーム間差分画像、血流フレーム間差分画像の一例を示す。図15Aに示す解析は、フレームレート7.5秒/枚、トータル線量(撮影10秒間での入射表面線量)が0.2mGy相当の動態画像を使用している。図15Aの上段は、2mm角でのビニング処理を施した場合、図15Aの下段は、2mm間隔での単純間引き処理を施した場合を示す。図15Bの上段は、2mm角でのビニング処理を施した場合、図15Bの下段は、2mm間隔での単純間引き処理を施した場合を示す。図15Cの上段は、2mm角でのビニング処理を施した場合、図15Cの下段は、2mm間隔での単純間引き処理を施した場合を示す。 FIGS. 15A to 15C show examples of histogram analysis results of the maximum flow velocity ratio, ventilation frame difference image, and blood flow frame difference image when binning processing and simple thinning processing are performed. The analysis shown in FIG. 15A uses a dynamic image corresponding to a frame rate of 7.5 seconds / frame and a total dose (incident surface dose in 10 seconds of imaging) equivalent to 0.2 mGy. The upper part of FIG. 15A shows a case where binning processing is performed at 2 mm square, and the lower part of FIG. 15A shows a case where simple thinning processing is performed at intervals of 2 mm. The upper part of FIG. 15B shows a case where binning processing is performed at 2 mm square, and the lower part of FIG. 15B shows a case where simple thinning processing is performed at intervals of 2 mm. The upper part of FIG. 15C shows the case where binning processing is performed at 2 mm square, and the lower part of FIG. 15C shows the case where simple thinning processing is performed at intervals of 2 mm.
図15A、図15Bに示すように、最大流速比のヒストグラム解析、換気フレーム間差分画像においては、解析結果に対するビニング処理と単純間引き処理は大差ない。即ち、換気の解析においては、ビニング処理と単純間引き処理の影響は大差ない。一方、図15Cに示すように、血流フレーム間差分画像においては、単純間引き処理ではビニング処理に比べて画像上のノイズが大きく、血流情報が検出できなくなる小領域(図15Cの肺野領域の中で静止画上にフレーム間差分値が加算されていない領域)が多くなる。よって、換気機能の解析の場合は、ビニング処理と単純間引きの何れを施してもよい(上記解析処理の個々の項目において、ビニング処理と単純間引きの何れを行ってもよい)が、血流機能の解析の場合は、ビニング処理を施すことが好ましい。 As shown in FIGS. 15A and 15B, in the histogram analysis of the maximum flow rate ratio and the difference image between ventilation frames, the binning process and the simple thinning process for the analysis result are not significantly different. That is, in the analysis of ventilation, the influence of the binning process and the simple thinning process is not much different. On the other hand, as shown in FIG. 15C, in the difference image between blood flow frames, the simple thinning process has a larger noise on the image than the binning process, and blood flow information cannot be detected (the lung field area in FIG. 15C). In this case, there are many areas in which no inter-frame difference value is added on a still image. Therefore, in the case of analysis of the ventilation function, either binning processing or simple thinning may be performed (in each item of the above analysis processing, either binning processing or simple thinning may be performed). In the case of this analysis, it is preferable to perform a binning process.
また、ビニング処理を施せば、小領域内における信号値の平均化により個々の画素のバラツキ影響が緩和されるので、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理等の補正処理を省略できる可能性が高くなり、処理時間を大幅に短縮化することができる。
また、ポジショニング確認用にビニング処理及び単純間引き処理された画像を解析用に再利用することができるので、処理時間を大幅に短縮化することができる。
これは、既存の動画撮影対応、例えば透視用のFPD等に於いては、FPD自体が、例えば、2×2画素等のビニング処理を行ったデータを表示装置に出力してリアルタイム表示を行い、手術具等の位置確認等に使用するものがあるが、このような型式の出力信号を、そのまま特徴量の解析に使用することが可能となり、既存の撮影装置の出力信号を使用しても動態に係る特徴量の解析が行えることを示している。
Also, if binning processing is performed, the influence of individual pixel variations is mitigated by averaging signal values within a small region, so that correction processing such as offset correction processing, gain correction processing, and defective pixel correction processing can be omitted. And the processing time can be greatly shortened.
Further, since the image subjected to binning processing and simple thinning processing for positioning confirmation can be reused for analysis, processing time can be greatly shortened.
For example, in the case of existing video shooting support, for example, FPD for fluoroscopy, the FPD itself outputs data subjected to binning processing such as 2 × 2 pixels to a display device for real-time display, Some types are used for confirming the position of surgical tools, etc., but this type of output signal can be used as it is for analysis of features, and even if the output signal of an existing imaging device is used It is shown that the feature amount related to can be analyzed.
なお、上記第1の実施の形態においては、各撮影室にコンソール5を配置し、各撮影室のコンソール5で撮影室内における撮影を制御することとして説明したが、図16に示す診断支援情報生成システム200のように、撮影室の外に1又は複数のコンソール5を設置して各撮影室の操作卓6及びアクセスポイントAPと接続可能な構成とし、各コンソール5が撮影室R1〜R3の何れの撮影についても制御できるようにしてもよい。 In the first embodiment, it has been described that the console 5 is arranged in each photographing room and the photographing in the photographing room is controlled by the console 5 in each photographing room. However, the diagnosis support information generation shown in FIG. As in the system 200, one or a plurality of consoles 5 are installed outside the shooting room so that the console 6 and the access point AP of each shooting room can be connected to each of the shooting rooms R1 to R3. It may be possible to control the shooting.
図16に示す診断支援情報生成システム200においては、例えば、各コンソール5の撮影管理テーブル521が「撮影室NO.」、「管球タイプ」、「ブッキーID」、「FPDID」、「画像受信時刻」等の項目を有し、各撮影室NO.に対応付けて当該撮影室NO.の撮影室に設けられている放射線源の管球タイプ及びブッキー装置のブッキーID、そのブッキー装置に現在装着されているFPDのFPDID、現在その撮影室に存在するFPDのFPDID、画像受信時刻等を記憶可能とする。各撮影室のブッキー装置においては、FPDが装着された際に、装着されたFPDからFPDIDを取得し、撮影室NO.、FPDID及び自己のブッキーIDを各コンソール5に送信する。同様に、各撮影室のクレードル4においては、FPDが装着された際に、装着されたFPDからFPDIDを取得し、撮影室NO.及びFPDIDを各コンソール5に送信する。このように構成すれば、各コンソール5において各撮影室R1〜R3の状況を把握し、指定された撮影オーダ情報に基づく撮影が可能な撮影室を判断することが可能となる。
例えば、各コンソール5において、動態撮影を指示する撮影オーダ情報が指定された際に、制御部51により撮影管理テーブル521に記憶されている情報、例えば、撮影室に存在する管球のタイプやFPDの種類の情報に基づいて各撮影室で動態撮影が可能であるか否かを事前に判断し、撮影が可能な撮影室の選択画面を表示部54に表示することが可能となる。又は、選択画面に各撮影室に設けられている放射線源の管球タイプ、ブッキー装置、FPDの種類等を表示し、操作者が撮影オーダ情報に基づく撮影が可能な撮影室を容易に選択できるようにしてもよい。或いは、操作者が撮影室を選択するための選択画面を表示部54に表示し、選択された撮影室での撮影が可能か否かを判断して撮影が不可能な場合に警告を表示することとしてもよい。このような制御により、コンソール5と撮影室とが1:1対応でなくm:nであっても、操作者は、間違いなく撮影オーダ情報に応じた撮影を行うことができる。また、操作者が動態撮影が不可能な放射線源やFPDで誤って撮影を開始してしまうことを防止することができる。撮影室が選択されると、制御部51は選択された撮影室の放射線源及びFPDを起動させる。
なお、診断支援情報生成システム200におけるその他の動作については診断支援情報生成システム100で説明したものと同様である。
In the diagnosis support information generation system 200 shown in FIG. 16, for example, the imaging management table 521 of each console 5 includes “imaging room No.”, “tube type”, “bookie ID”, “FPDID”, “image reception time”. ”And the like. In association with the shooting room NO. The tube type of the radiation source provided in the radiographing room, the bucky ID of the bucky device, the FPDID of the FPD currently installed in the bucky device, the FPDID of the FPD currently in the radiographing room, the image reception time, etc. Make memorable. When the FPD is installed in the shooting device in each shooting room, the FPDID is acquired from the mounted FPD, and the shooting room NO. , The FPDID and its own Bucky ID are transmitted to each console 5. Similarly, in the cradle 4 of each photographing room, when the FPD is attached, the FPDID is acquired from the attached FPD, and the photographing room NO. And FPDID are transmitted to each console 5. If comprised in this way, it will become possible to grasp | ascertain the condition of each imaging | photography room R1-R3 in each console 5, and to judge the imaging | photography room which can image | photograph based on the designated imaging | photography order information.
For example, the information stored in the imaging management table 521 by the control unit 51 when the imaging order information for instructing dynamic imaging is specified in each console 5, for example, the type of tube existing in the imaging room or the FPD Based on this type of information, it is possible to determine in advance whether or not dynamic shooting is possible in each shooting room, and to display a selection screen for a shooting room in which shooting can be performed on the display unit 54. Alternatively, the tube type of the radiation source provided in each imaging room, the bucky device, the type of FPD, etc. are displayed on the selection screen, and the operator can easily select an imaging room in which imaging can be performed based on imaging order information. You may do it. Alternatively, a selection screen for the operator to select a shooting room is displayed on the display unit 54, and it is determined whether or not shooting in the selected shooting room is possible, and a warning is displayed when shooting is impossible. It is good as well. With such control, even if the console 5 and the photographing room are not 1: 1 and m: n, the operator can definitely perform photographing according to the photographing order information. In addition, it is possible to prevent the operator from erroneously starting imaging with a radiation source or FPD for which dynamic imaging is impossible. When the imaging room is selected, the control unit 51 activates the radiation source and FPD of the selected imaging room.
Other operations in the diagnosis support information generation system 200 are the same as those described in the diagnosis support information generation system 100.
[第2の実施の形態]
次に、本発明の第2の実施の形態について説明する。
まず、構成について説明する。
図17に、第2の実施の形態における診断支援情報生成システム300の全体構成を示す。診断支援情報生成システム300は、開業医やクリニック等の小規模施設において撮影された動態画像を解析センタの解析サーバ30で解析し、小規模施設に解析結果を提供するシステムである。第1の実施の形態で説明したように、ワーピング処理を施さなくてよい場合、解析処理を行う側でFPDの画素サイズや個々の画素のダイナミックレンジ等を考慮する必要がない。そのため、各施設で使用しているFPDの型式を問わず、以下に説明するようにオープンシステムで解析処理を提供することができる。
[Second Embodiment]
Next, a second embodiment of the present invention will be described.
First, the configuration will be described.
FIG. 17 shows an overall configuration of a diagnosis support information generation system 300 according to the second embodiment. The diagnosis support information generation system 300 is a system that analyzes a dynamic image taken at a small-scale facility such as a medical practitioner or a clinic by the analysis server 30 of the analysis center and provides an analysis result to the small-scale facility. As described in the first embodiment, when it is not necessary to perform the warping process, it is not necessary to consider the pixel size of the FPD, the dynamic range of each pixel, or the like on the analysis process side. Therefore, regardless of the type of FPD used in each facility, analysis processing can be provided by an open system as described below.
図17に示すように、小規模施設内の受付には受付装置20が設けられ、診察室にはコンソール50が設けられ、撮影室には第1の実施の形態で説明した撮影室内の装置(ブッキー装置1、2、放射線源3a、FPD9a等)が設けられている。コンソール50は、受付装置20及びブッキー装置とLAN等を介して接続されている。また、コンソール50は、インターネットPN等を介して解析サーバ30と接続可能に構成されている。 As shown in FIG. 17, a reception device 20 is provided for reception in a small-scale facility, a console 50 is provided in an examination room, and an apparatus (in the imaging room described in the first embodiment) ( Bucky devices 1 and 2, radiation source 3a, FPD 9a, etc.) are provided. The console 50 is connected to the accepting device 20 and the bucky device via a LAN or the like. The console 50 is configured to be connectable to the analysis server 30 via the Internet PN or the like.
受付装置20は、来院した患者の受付登録、会計計算、保険点数計算等を行うためのコンピュータ装置である。受付装置20は、受付番号、患者情報(「患者ID」、「氏名」、「生年月日」、「性別」、「年齢」、「住所」、「電話番号」、「被保険者番号」等)が入力されると、入力された受付番号及び患者情報をコンソール50に送信する。 The reception apparatus 20 is a computer apparatus for performing reception registration, accounting calculation, insurance score calculation, and the like of a patient who has visited the hospital. The reception device 20 includes a reception number, patient information (“patient ID”, “name”, “birth date”, “gender”, “age”, “address”, “phone number”, “insured person number”, etc. ) Is input, the input reception number and patient information are transmitted to the console 50.
コンソール50は、撮影室内の各装置を制御して撮影を行う。また、解析サーバ30と通信接続し、動態画像の解析結果を取得し、表示する。
コンソール50の構成は、図3に示すコンソール5と同様に、制御部51、記憶部52、入力部53、表示部54、通信I/F55、ネットワーク通信部56等を備えて構成され、各部はバス57により接続されている。
コンソール50の記憶部52には、コンソール50に対応する各種プログラムが記憶されており、制御部51は、当該プログラムに従って後述する撮影・解析処理Bをはじめとする各種処理を実行する。なお、コンソール50においては、撮影管理テーブル521は不要である。また、ネットワーク通信部56は、スイッチングハブを介して接続された施設内の装置のほか、インターネットPNを介して解析サーバ30等の外部機器と通信接続可能である。その他コンソール50の各部の構成は、コンソール5と同様であるので説明を援用する。
The console 50 performs imaging by controlling each device in the imaging room. Further, the communication server 30 is connected to the analysis server 30 to acquire and display the dynamic image analysis result.
The configuration of the console 50 includes a control unit 51, a storage unit 52, an input unit 53, a display unit 54, a communication I / F 55, a network communication unit 56, etc., similar to the console 5 shown in FIG. They are connected by a bus 57.
Various programs corresponding to the console 50 are stored in the storage unit 52 of the console 50, and the control unit 51 executes various processes including an imaging / analysis process B described later according to the program. In the console 50, the imaging management table 521 is not necessary. Further, the network communication unit 56 can be communicably connected to an external device such as the analysis server 30 via the Internet PN in addition to the devices in the facility connected via the switching hub. The configuration of each part of the console 50 is the same as that of the console 5, so the description is cited.
撮影室内の各装置の構成は、第1の実施の形態で説明したものと同様であるので説明を援用する。 Since the configuration of each apparatus in the photographing room is the same as that described in the first embodiment, the description will be cited.
解析サーバ30は、外部の解析センタに設けられたサーバ装置である。解析サーバ30は、CPU等により構成される制御部、解析用プログラムが記憶された記憶部、入力部、記憶部、通信部等を備えて構成されている。解析サーバ30は、コンソール50からの要求に応じて制御部と記憶部に記憶されている解析用プログラムとの協働により撮影・解析処理Bを行い、解析結果をコンソール50に送信する。なお、解析サーバ30は、上述の解析用WS8と同様に、ワーピング処理を施さない。ワーピング処理を施す場合、解析サーバ30側で各社のFPDの画素サイズや個々の画素のダイナミックレンジ、照射量等が必要となるが、ワーピング処理を省略するので各社のFPDに対応することができる。また、解析処理の処理時間を大幅に低減することができる。 The analysis server 30 is a server device provided in an external analysis center. The analysis server 30 includes a control unit configured by a CPU and the like, a storage unit in which an analysis program is stored, an input unit, a storage unit, a communication unit, and the like. In response to a request from the console 50, the analysis server 30 performs imaging / analysis processing B in cooperation with the control unit and the analysis program stored in the storage unit, and transmits the analysis result to the console 50. Note that the analysis server 30 does not perform the warping process, similar to the analysis WS 8 described above. When performing the warping process, the analysis server 30 side needs the pixel size of each company's FPD, the dynamic range of each pixel, the irradiation amount, etc., but since the warping process is omitted, it can cope with the FPD of each company. Moreover, the processing time of the analysis process can be significantly reduced.
次に、診断支援情報生成システム300における動作について説明する。
上述のように、患者が来院すると、受付担当により患者に受付番号が付与され、受付装置20において受付番号、患者情報等が入力される。受付装置20の制御部においては、患者の受付番号及び患者情報等が入力部により入力されると、入力された情報が記憶部に記憶されて受付登録が行われるとともに、入力された情報(受付リスト情報)が通信部によりコンソール50に送信される。コンソール50の制御部においては、受付装置20からネットワーク通信部56により受付リスト情報が受信されると、受信された受付リスト情報が記憶部52に記憶される。また、入力部53からの操作に応じて、本日の受付リスト情報の一覧が表示部54に表示される。
Next, the operation in the diagnosis support information generation system 300 will be described.
As described above, when the patient comes to the hospital, a reception number is given to the patient by the receptionist, and the reception number, patient information, and the like are input in the reception device 20. In the control unit of the reception device 20, when the patient reception number, patient information, and the like are input by the input unit, the input information is stored in the storage unit, reception registration is performed, and the input information (reception List information) is transmitted to the console 50 by the communication unit. In the control unit of the console 50, when reception list information is received from the reception device 20 by the network communication unit 56, the received reception list information is stored in the storage unit 52. Further, according to an operation from the input unit 53, a list of today's reception list information is displayed on the display unit 54.
受付番号が付与された患者が診察室に移動すると、医師は、入力部53を操作して、表示部54に表示された受付リスト情報から診察対象の患者の受付リスト情報を選択する。コンソール50においては、受付リスト情報から一の受付リストが選択されると、選択された受付リストに対応する患者の診断用画面が表示部54に表示される。診断用画面は、例えば、当該患者について撮影された画像の表示欄、診療情報の入力欄等を有する画面である。次いで、医師は、患者に対し問診を行い、実施する撮影を決定する。問診により動態画像の撮影及び解析が必要であると判断した場合には、医師は入力部53を操作して、コンソール50に以下の撮影・解析処理Bを実行させる。 When the patient to whom the reception number is assigned moves to the examination room, the doctor operates the input unit 53 to select the reception list information of the patient to be examined from the reception list information displayed on the display unit 54. In the console 50, when one reception list is selected from the reception list information, a patient diagnosis screen corresponding to the selected reception list is displayed on the display unit 54. The diagnosis screen is, for example, a screen having a display field for images taken for the patient, an input field for medical information, and the like. Next, the doctor interrogates the patient and decides the imaging to be performed. When it is determined that the dynamic image needs to be captured and analyzed through the inquiry, the doctor operates the input unit 53 to cause the console 50 to execute the following imaging / analysis process B.
以下、撮影・解析処理Bについて説明する。
図18に、診断支援情報生成システム300において実行される動態解析処理の流れを示す。撮影・解析処理Bのコンソール50側の処理は、コンソール50の制御部51と記憶部52に記憶されているプログラムとの協働により実行される。解析サーバ30側の処理は、解析サーバ30の制御部と記憶部に記憶されているプログラムとの協働より実行される。
The imaging / analysis process B will be described below.
FIG. 18 shows a flow of dynamic analysis processing executed in the diagnosis support information generation system 300. The process on the console 50 side of the imaging / analysis process B is executed in cooperation with the control unit 51 of the console 50 and the program stored in the storage unit 52. The processing on the analysis server 30 side is executed by the cooperation of the control unit of the analysis server 30 and the program stored in the storage unit.
まず、コンソール50において、ネットワーク通信部56により解析サーバ30へのアクセスが行われ、動態解析開始が要求される(ステップT1)。
解析サーバ30においては、動態解析開始要求が受信されると、コンソール50に対し、解析対象の部位、特徴量、及び使用するFPDのフレームレートの問い合わせが行われる(ステップT2)。
コンソール50においては、解析対象の部位、特徴量(解析対象項目)、及び使用するFPDのフレームレートの情報が取得され、ネットワーク通信部56により解析サーバ30に送信される(ステップT3)。ステップT3では、解析サーバ30からの問い合わせ内容が表示部54に表示され、医師が解析対象の部位及び特徴量、使用するFPD9aのフレームレートを入力部53により入力すると、入力した情報がネットワーク通信部56により解析サーバ30に送信される。特徴量としては、例えば、解析対象が肺野である場合、第1の実施の形態の図6のステップS19において説明した(1)〜(15)の項目特徴量が挙げられる。
First, in the console 50, the network communication unit 56 accesses the analysis server 30, and a dynamic analysis start is requested (step T1).
In the analysis server 30, when the dynamic analysis start request is received, the console 50 is inquired about the part to be analyzed, the feature quantity, and the frame rate of the FPD to be used (step T2).
In the console 50, information on a part to be analyzed, a feature amount (analysis target item), and an FPD frame rate to be used are acquired and transmitted to the analysis server 30 by the network communication unit 56 (step T3). In step T3, the content of the inquiry from the analysis server 30 is displayed on the display unit 54. When the doctor inputs the part to be analyzed, the feature amount, and the frame rate of the FPD 9a to be used by the input unit 53, the input information is the network communication unit. 56 to the analysis server 30. As the feature amount, for example, when the analysis target is a lung field, the item feature amounts (1) to (15) described in step S19 of FIG. 6 of the first embodiment may be mentioned.
解析サーバ30においては、コンソール50から解析対象の部位、特徴量及びフレームレートが受信されると、解析に必要なフレーム画像数が算出され、コンソール50に通知される(ステップT4)。解析に必要なフレーム画像数の算出については、例えば、図6のステップS17で説明したのと同様である。即ち、解析サーバ30の記憶部に平均的な成人等の動態周期を部位毎に記憶しておき、フレームレートと解析対象の部位の動態周期に基づいて解析に必要なフレーム画像数が算出される。 In the analysis server 30, when the part to be analyzed, the feature amount, and the frame rate are received from the console 50, the number of frame images necessary for the analysis is calculated and notified to the console 50 (step T4). The calculation of the number of frame images necessary for the analysis is the same as that described in step S17 of FIG. 6, for example. In other words, the dynamic cycle of an average adult or the like is stored for each region in the storage unit of the analysis server 30, and the number of frame images necessary for analysis is calculated based on the frame rate and the dynamic cycle of the region to be analyzed. .
コンソール50においては、撮影に必要なフレーム画像数が通知されると、通知されたフレーム画像数、撮影部位、***等が入力部53により入力される(ステップT5)。そして、入力された情報に基づいて、撮影処理が行われる(ステップT6)。
撮影処理は、図6のステップS4〜S13(又はS14)の処理と略同様の処理である。ここでは、医師や撮影技師等の撮影実施者は、患者を撮影室に連れて行き、撮影する***のブッキー装置(1又は2)にFPD9aを装着する。コンソール50においては、入力されたフレーム画像数、撮影部位、***に基づいて、放射線源3a及びブッキー装置(1又は2)が起動され、放射線源3aの位置及び向きが調整される。また、入力されたフレーム画像数、撮影部位、***に基づいて、放射線源3aに放射線照射条件が設定され、FPD9aに画像読取条件が設定される。動態撮影の結果を用いて解析を行う場合は、診断に使用し得る解析精度を確保するため、フレームレートが3.75枚/秒以上に設定される。操作卓6から放射線照射指示が入力されると、コンソール50により放射線源3a及びFPD9aが制御され、動態撮影が行われる。入力されたフレーム画像数+αの撮影が終了すると、撮影動作が停止される。
撮影により取得されたフレーム画像は順次FPD9aのコネクタ94によりブッキー装置を介してコンソール50に入力される。入力されたフレーム画像は、記憶部52に記憶され、間引き処理が行われる。間引き処理は、ビニング処理及び/又は単純間引き処理である。何れの間引き処理を行うか、及び、ビニングする画素数(ブロックサイズ)もしくは単純間引きの画素間隔は、解析サーバ30が受信した解析対象の部位及び特徴量に応じた値が、フレーム画像数と併せて、撮影前に解析サーバ30からコンソール50へ通知されるようにしてもよい。間引き後の画像は表示部54に表示される。撮影実施者は、表示部54に表示されたフレーム画像を見てポジショニング等の確認を行い、撮影により診断に適した画像が取得された(撮影OK)か、再撮影が必要(撮影NG)かの判断結果を入力部53により入力する。尚、撮影により取得された各フレーム画像は、一旦FPD9aの記憶部93に記憶され、全撮影の終了後に纏めてFPD9aからコンソール50へ出力されるようにしても良い。入力部53により撮影NGの判断結果が入力されると、入力されたフレーム画像が記憶部52から削除される。この場合、再撮影が行われる。入力部53の所定の操作により撮影OKを示す判断結果が入力されると、処理はステップT7へ移行する。
In the console 50, when the number of frame images necessary for imaging is notified, the notified number of frame images, imaging region, body position, and the like are input by the input unit 53 (step T5). Then, a photographing process is performed based on the input information (step T6).
The photographing process is substantially the same as the process in steps S4 to S13 (or S14) in FIG. Here, an imaging operator such as a doctor or an imaging engineer takes the patient to the imaging room and attaches the FPD 9a to the bucky device (1 or 2) in a posture for imaging. In the console 50, the radiation source 3a and the bucky device (1 or 2) are activated based on the input number of frame images, imaging region, and body position, and the position and orientation of the radiation source 3a are adjusted. Further, based on the input number of frame images, imaging region, and body position, radiation irradiation conditions are set for the radiation source 3a, and image reading conditions are set for the FPD 9a. When analysis is performed using the result of dynamic imaging, the frame rate is set to 3.75 frames / second or more in order to ensure analysis accuracy that can be used for diagnosis. When a radiation irradiation instruction is input from the console 6, the radiation source 3a and the FPD 9a are controlled by the console 50, and dynamic imaging is performed. When shooting of the number of input frame images + α is completed, the shooting operation is stopped.
The frame images obtained by photographing are sequentially input to the console 50 via the bucky device by the connector 94 of the FPD 9a. The input frame image is stored in the storage unit 52, and a thinning process is performed. The thinning process is a binning process and / or a simple thinning process. Which decimation process is performed, and the number of pixels to be binned (block size) or the simple decimation pixel interval, the value according to the analysis target part and the feature quantity received by the analysis server 30 are combined with the number of frame images. Thus, the console may be notified from the analysis server 30 before shooting. The thinned image is displayed on the display unit 54. The photographer confirms positioning and the like by looking at the frame image displayed on the display unit 54, and whether an image suitable for diagnosis has been acquired by photographing (photographing OK) or re-photographing is necessary (photographing NG). The determination result is input by the input unit 53. Each frame image acquired by shooting may be temporarily stored in the storage unit 93 of the FPD 9a, and may be output from the FPD 9a to the console 50 after completion of all shooting. When the determination result of shooting NG is input by the input unit 53, the input frame image is deleted from the storage unit 52. In this case, re-photographing is performed. When a determination result indicating photographing OK is input by a predetermined operation of the input unit 53, the process proceeds to step T7.
図18のステップT7においては、コンソール50において、解析に使用するフレーム画像の選択が行われる(ステップT7)。ここで、上述のように解析に必要となる枚数以上のフレーム画像を用いて解析を行うと、データ転送時間や解析処理時間がかかり、好ましくない。処理時間短縮のためには、撮影された一連のフレーム画像の中から解析に使用する枚数のフレーム画像を選択することが必要となる。一方、動態画像の解析は、動態周期や隣接するフレーム画像間での信号値の差分等に基づき行われるので、動態画像の解析結果を被写体の真の特徴量を表す結果とするには、被写体の動態1サイクル以上の枚数の、連続する一連のフレーム画像が必要である。そこで、ステップT7においては、解析に必要な枚数分の連続するフレーム画像を操作者が正しく選択するためのGUIとして選択用画面541を提供する。具体的には、選択枠541aを用いて選択させるようにすることで、連続しない1コマ毎の画像選択や、一周期に満たない選択を防止する。選択用画面541及びその動作については、第1の実施形態で説明したものと同様であるので説明を援用する。 In step T7 of FIG. 18, the console 50 selects a frame image to be used for analysis (step T7). Here, if the analysis is performed using the frame images more than the number necessary for the analysis as described above, it takes a data transfer time and an analysis processing time, which is not preferable. In order to shorten the processing time, it is necessary to select the number of frame images to be used for analysis from a series of photographed frame images. On the other hand, analysis of dynamic images is performed based on a dynamic cycle, a difference in signal values between adjacent frame images, and the like. A continuous series of frame images of at least one cycle is required. Therefore, in step T7, a selection screen 541 is provided as a GUI for the operator to correctly select the continuous frame images as many as necessary for analysis. Specifically, by using the selection frame 541a to select, it is possible to prevent image selection for each frame that is not continuous or selection that is less than one cycle. Since the selection screen 541 and the operation thereof are the same as those described in the first embodiment, the description is incorporated.
解析に使用するフレーム画像群が選択されると、送信する画像数がネットワーク通信部56により解析サーバ30に通知される(ステップT8)。解析サーバ30においては、通知された画像数の画像データの受信が可能な状態となると、送信許可通知がコンソール50に送信される(ステップT9)。コンソール50においては、ネットワーク通信部56により解析サーバ30から送信許可通知が受信されると、選択されたフレーム画像群の間引きデータが解析サーバ30へ送信される(ステップT10)。なお、選択された一連のフレーム画像群の間引きデータに、患者情報、撮影順を示す番号等の情報が付帯される(例えば、DICOMマルチフォーマットのファイルフォーマットで画像データのヘッダ領域に書き込まれる)。 When a group of frame images to be used for analysis is selected, the number of images to be transmitted is notified to the analysis server 30 by the network communication unit 56 (step T8). When the analysis server 30 is ready to receive image data for the notified number of images, a transmission permission notification is transmitted to the console 50 (step T9). In the console 50, when the transmission permission notification is received from the analysis server 30 by the network communication unit 56, the thinned data of the selected frame image group is transmitted to the analysis server 30 (step T10). In addition, information such as patient information and a number indicating the imaging order is attached to the thinned data of the selected series of frame image groups (for example, written in the header area of the image data in the DICOM multi-format file format).
解析サーバ30においては、コンソール50より通知された枚数の画像の取得が終了すると、取得完了通知がコンソール50に送信される(ステップT11)。そして、解析処理が実行される(ステップT11)。解析処理においては、第1の実施の形態で説明したように、ワーピング処理を行わないことが処理時間の短縮化の点から好ましい。解析処理の例は、第1の実施の形態において図6のステップS19で説明したものと同様であるので説明を援用する。 In the analysis server 30, when the acquisition of the number of images notified from the console 50 is completed, an acquisition completion notification is transmitted to the console 50 (step T11). Then, analysis processing is executed (step T11). In the analysis process, as described in the first embodiment, it is preferable not to perform the warping process from the viewpoint of shortening the processing time. An example of the analysis process is the same as that described in step S19 of FIG. 6 in the first embodiment, and thus description is incorporated.
解析が終了すると、解析結果のデータがコンソール50に送信される。コンソール50においては、解析結果のデータが診断対象として選択されている患者情報と対応付けて記憶部52に記憶され、撮影・解析処理Bは終了する。 When the analysis is completed, the analysis result data is transmitted to the console 50. In the console 50, the analysis result data is stored in the storage unit 52 in association with the patient information selected as the diagnosis target, and the imaging / analysis process B ends.
第2の実施の形態の診断支援情報生成システム300によれば、間引き処理後のフレーム画像が解析サーバ30に送信されるので、撮影された画像データ全体を送信する場合に比べて通信時間を抑え、短時間で解析サーバ30に解析に必要な画像データを送信することができる。間引きされたフレーム画像としては、ポジショニング確認用にコンソール5で作成されたものを再使用することができるので、わざわざ再度間引き処理を行う必要がない。また、解析サーバ30において、ワーピング処理を施さずに間引きデータを用いて解析処理を行うことで、画像データ全体を用いてワーピング処理を行って解析処理を行う場合に比べて大幅に解析処理の処理時間を短縮することが可能となり、応答性のよい解析サービスを提供することが可能となる。更に、コンソール50において、送信するフレーム画像からROI領域を抽出し、抽出されたROI領域の画像データのみを解析サーバ30に送信することとすると、通信時間及び解析サーバ30での処理時間を短縮することができ、より早期に解析結果を取得することが可能となる。 According to the diagnosis support information generation system 300 of the second embodiment, the frame image after the thinning process is transmitted to the analysis server 30, so that the communication time is reduced compared to the case where the entire captured image data is transmitted. The image data necessary for the analysis can be transmitted to the analysis server 30 in a short time. As the thinned frame image, a frame image created by the console 5 for positioning confirmation can be reused, so that it is not necessary to perform the thinning process again. Further, the analysis server 30 performs the analysis process using the thinned data without performing the warping process, so that the process of the analysis process is greatly compared with the case where the analysis process is performed by performing the warping process using the entire image data. The time can be shortened, and an analysis service with good responsiveness can be provided. Furthermore, when the console 50 extracts the ROI area from the frame image to be transmitted and transmits only the image data of the extracted ROI area to the analysis server 30, the communication time and the processing time in the analysis server 30 are shortened. It is possible to obtain the analysis result earlier.
以上説明したように、診断支援情報生成システムによれば、動態撮影を行う場合、コンソール5は、FPD9aにフレームレート3.75枚/秒の画像読取条件を設定する。FPD9aは、設定されたフレームレートで撮影を行う。
従って、撮影対象の部位に応じて、診断に使用し得る解析結果が出力でき、処理時間及び患者の被曝を最低限に抑え、ハードウエアの構成をいたずらに増大させることのない動態診断支援情報生成システムを提供することができる。
As described above, according to the diagnosis support information generation system, when performing dynamic imaging, the console 5 sets an image reading condition with a frame rate of 3.75 frames / second in the FPD 9a. The FPD 9a performs shooting at a set frame rate.
Therefore, analysis results that can be used for diagnosis can be output according to the region to be imaged, generation of dynamic diagnosis support information that minimizes processing time and patient exposure, and does not unnecessarily increase the hardware configuration. A system can be provided.
また、コンソール5は、撮影により得られた複数のフレーム画像の画素数を低減するための間引き処理を施して解析用WS8に送信し、解析用WS8は間引きにより画素数が低減されたフレーム画像を用いて特徴量を算出するので、演算に要するデータ量を大幅に 削減し、解析に要する処理時間を大幅に短縮することができる。 Further, the console 5 performs a thinning process for reducing the number of pixels of a plurality of frame images obtained by photographing and transmits the result to the analysis WS8. The analysis WS8 receives the frame image with the number of pixels reduced by the thinning. Since the feature value is calculated using this method, the amount of data required for computation can be greatly reduced, and the processing time required for analysis can be greatly reduced.
また、解析用WS8及び解析サーバ30は、被写体の動態1サイクル以上のフレーム画像群に基づいて被写体の動態の特徴量を算出するので、動態の正しい解析結果を得ることができる。 In addition, the analysis WS 8 and the analysis server 30 calculate the feature amount of the subject's dynamics based on the frame image group of one or more cycles of the subject's dynamics, so that a correct analysis result of the dynamics can be obtained.
被写体の呼吸運動に関する特徴量、例えば、被写体がCOPDであるか否かの指標となる特徴量である最大気流速度の比(呼気期間、吸気期間のそれぞれにおける気流速度の比)等を算出する場合は、フレームレートを3.75枚/秒とすれば、ハードウエアの構成をいたずらに増大させることなく、診断に使用し得る解析結果が出力でき、処理時間及び患者の被曝を最低限に抑えることが可能となる。 When calculating a feature quantity related to the breathing motion of the subject, for example, a ratio of the maximum airflow velocity (a ratio of airflow velocity in each of the expiration period and the inspiration period), which is a feature quantity indicating whether the subject is COPD or not. If the frame rate is 3.75 frames / second, analysis results that can be used for diagnosis can be output without unnecessarily increasing the hardware configuration, minimizing processing time and patient exposure. Is possible.
なお、上記実施の形態は本発明の好適な一例であり、これに限定されない。
例えば、上記実施の形態においては、コンソール5やコンソール50において間引き処理を行うこととして説明したが、例えば、特許第4,546,174号公報のように、動態画像を撮影する撮影装置側(本実施の形態ではFPD9a)でビニング処理や単純間引き処理を行い、処理済みのフレーム画像をコンソール5に送信することとしてもよい。このようにすれば、FPDとコンソールとの間での画像データの転送時間も短縮することができるので、更に好ましい。
In addition, the said embodiment is a suitable example of this invention, and is not limited to this.
For example, in the above-described embodiment, the thinning process is performed in the console 5 or the console 50. However, for example, as in Japanese Patent No. 4,546,174, the photographing apparatus side (this book) In the embodiment, binning processing or simple thinning processing may be performed by the FPD 9 a), and the processed frame image may be transmitted to the console 5. This is more preferable because the transfer time of the image data between the FPD and the console can be shortened.
また、上記第1の実施の形態においては、解析用WS8をコンソール5とは別体で設ける構成としたが、コンソール5が解析用プログラムを備えて解析処理を行う構成としてもよい。このようにすれば、コンソール5から解析用WS8へ画像データを送信する時間を省くことが可能となるとともに、他のコンソール5による解析によって処理が遅れるといった事態を防止することができる。 In the first embodiment, the analysis WS 8 is provided separately from the console 5. However, the console 5 may include an analysis program and perform analysis processing. In this way, it is possible to save time for transmitting image data from the console 5 to the analysis WS 8 and to prevent a situation in which processing is delayed due to analysis by another console 5.
また、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてHDDや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、CD-ROM等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。 In the above description, an example in which an HDD, a semiconductor nonvolatile memory, or the like is used as a computer-readable medium for the program according to the present invention is disclosed, but the present invention is not limited to this example. As another computer-readable medium, a portable recording medium such as a CD-ROM can be applied. Further, a carrier wave is also applied as a medium for providing program data according to the present invention via a communication line.
その他、診断支援情報生成システムを構成する各装置の細部構成及び細部動作に関しても、発明の趣旨を逸脱することのない範囲で適宜変更可能である。 In addition, the detailed configuration and detailed operation of each device constituting the diagnosis support information generation system can be changed as appropriate without departing from the spirit of the invention.
100 診断支援情報生成システム
1 ブッキー装置
11 制御部
12 検出器装着部
13 通信I/F
14 駆動部
15 バス
2 ブッキー装置
21 制御部
22 検出器装着部
23 通信I/F
24 駆動部
25 バス
3a 放射線源
3b 放射線源
3c 放射線源
4 クレードル
5 コンソール
51 制御部
52 記憶部
521 撮影管理テーブル
53 入力部
54 表示部
541 選択用画面
55 通信I/F
56 ネットワーク通信部
57 バス
9a FPD
9b FPD
91 制御部
92 検出部
93 記憶部
94 コネクタ
95 バッテリ
96 無線通信部
97 バス
6 操作卓
7 HIS/RIS
8 解析用WS
10 PACS
20 受付装置
50 コンソール
30 解析サーバ
200 診断支援情報生成システム
300 診断支援情報生成システム
DESCRIPTION OF SYMBOLS 100 Diagnosis support information generation system 1 Bucky device 11 Control part 12 Detector mounting part 13 Communication I / F
14 Drive unit 15 Bus 2 Bucky device 21 Control unit 22 Detector mounting unit 23 Communication I / F
24 driving unit 25 bus 3a radiation source 3b radiation source 3c radiation source 4 cradle 5 console 51 control unit 52 storage unit 521 imaging management table 53 input unit 54 display unit 541 selection screen 55 communication I / F
56 Network communication unit 57 Bus 9a FPD
9b FPD
91 control unit 92 detection unit 93 storage unit 94 connector 95 battery 96 wireless communication unit 97 bus 6 console 7 HIS / RIS
8 WS for analysis
10 PACS
20 reception device 50 console 30 analysis server 200 diagnosis support information generation system 300 diagnosis support information generation system
Claims (6)
2次元状に配置された複数の検出素子を有し、前記複数の検出素子のそれぞれにおいて前記放射線源によりパルス照射された放射線を検出して順次フレーム画像を生成する放射線検出器と、
前記放射線源及び前記放射線検出器により前記被写体の動態を撮影することにより得られた複数のフレーム画像に基づいて前記被写体の動態の特徴量を算出する解析手段と、
を備える動態診断支援情報生成システムであって、
前記放射線検出器におけるフレームレートが3.75枚/秒以上である動態診断支援情報生成システム。 A radiation source capable of pulse irradiation;
A radiation detector that has a plurality of detection elements arranged in a two-dimensional shape, detects radiation irradiated by the radiation source at each of the plurality of detection elements, and sequentially generates a frame image;
Analyzing means for calculating feature quantities of the dynamics of the subject based on a plurality of frame images obtained by photographing the dynamics of the subject by the radiation source and the radiation detector;
A dynamic diagnosis support information generation system comprising:
A dynamic diagnosis support information generation system in which a frame rate in the radiation detector is 3.75 frames / second or more.
前記解析手段は、前記間引き処理手段により間引き処理された前記複数のフレーム画像に基づいて特徴量の算出を行う請求項1に記載の動態診断支援情報生成システム。 Comprising a thinning processing means for performing a thinning process for reducing the number of pixels of the plurality of frame images;
The dynamic analysis support information generation system according to claim 1, wherein the analysis unit calculates a feature amount based on the plurality of frame images thinned by the thinning processing unit.
前記解析手段は、前記被写体の呼吸運動に関する特徴量を算出する請求項1〜4の何れか一項に記載の動態診断支援情報生成システム。 The dynamic diagnosis support information generation system performs dynamic imaging of the chest with the radiation source and the radiation detector,
The dynamic analysis support information generation system according to any one of claims 1 to 4, wherein the analysis unit calculates a feature amount related to respiratory motion of the subject.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010259650A JP5625800B2 (en) | 2010-11-22 | 2010-11-22 | Dynamic diagnosis support information generation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010259650A JP5625800B2 (en) | 2010-11-22 | 2010-11-22 | Dynamic diagnosis support information generation system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012110399A true JP2012110399A (en) | 2012-06-14 |
JP5625800B2 JP5625800B2 (en) | 2014-11-19 |
Family
ID=46495313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010259650A Active JP5625800B2 (en) | 2010-11-22 | 2010-11-22 | Dynamic diagnosis support information generation system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5625800B2 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014091977A1 (en) * | 2012-12-12 | 2014-06-19 | コニカミノルタ株式会社 | Image-processing device and program |
US9901317B2 (en) | 2013-05-31 | 2018-02-27 | Konica Minolta, Inc. | Image processing apparatus for acquiring a dynamic image and storage medium |
JP2018110637A (en) * | 2017-01-10 | 2018-07-19 | コニカミノルタ株式会社 | Dynamic image processing system |
JP2020171475A (en) * | 2019-04-10 | 2020-10-22 | コニカミノルタ株式会社 | Dynamic image analysis apparatus, dynamic image analysis method, and program |
JP7078098B1 (en) | 2020-12-14 | 2022-05-31 | コニカミノルタ株式会社 | Mobile radiography equipment, programs, methods, controls and radiography systems |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004000412A (en) * | 2002-04-03 | 2004-01-08 | Canon Inc | Method, device, and system for radiography, program, computer-readable storage medium, and method and system for supporting image diagnosis |
JP2004312434A (en) * | 2003-04-08 | 2004-11-04 | Canon Inc | Image processing apparatus and method, and system |
JP2009050529A (en) * | 2007-08-28 | 2009-03-12 | Konica Minolta Medical & Graphic Inc | Radiological imaging system and support apparatus for radiological imaging |
JP2009153678A (en) * | 2007-12-26 | 2009-07-16 | Konica Minolta Medical & Graphic Inc | Kinetic image processing system |
JP2009172190A (en) * | 2008-01-25 | 2009-08-06 | Konica Minolta Holdings Inc | Image generating apparatus, program, and image generating method |
JP2009273605A (en) * | 2008-05-14 | 2009-11-26 | Konica Minolta Medical & Graphic Inc | Dynamic image diagnosing support system |
JP4546174B2 (en) * | 2004-07-15 | 2010-09-15 | キヤノン株式会社 | Imaging apparatus, control method therefor, and program |
-
2010
- 2010-11-22 JP JP2010259650A patent/JP5625800B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004000412A (en) * | 2002-04-03 | 2004-01-08 | Canon Inc | Method, device, and system for radiography, program, computer-readable storage medium, and method and system for supporting image diagnosis |
JP2004312434A (en) * | 2003-04-08 | 2004-11-04 | Canon Inc | Image processing apparatus and method, and system |
JP4546174B2 (en) * | 2004-07-15 | 2010-09-15 | キヤノン株式会社 | Imaging apparatus, control method therefor, and program |
JP2009050529A (en) * | 2007-08-28 | 2009-03-12 | Konica Minolta Medical & Graphic Inc | Radiological imaging system and support apparatus for radiological imaging |
JP2009153678A (en) * | 2007-12-26 | 2009-07-16 | Konica Minolta Medical & Graphic Inc | Kinetic image processing system |
JP2009172190A (en) * | 2008-01-25 | 2009-08-06 | Konica Minolta Holdings Inc | Image generating apparatus, program, and image generating method |
JP2009273605A (en) * | 2008-05-14 | 2009-11-26 | Konica Minolta Medical & Graphic Inc | Dynamic image diagnosing support system |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014091977A1 (en) * | 2012-12-12 | 2014-06-19 | コニカミノルタ株式会社 | Image-processing device and program |
JP5578297B1 (en) * | 2012-12-12 | 2014-08-27 | コニカミノルタ株式会社 | Image processing apparatus and program |
CN104853677A (en) * | 2012-12-12 | 2015-08-19 | 柯尼卡美能达株式会社 | Image-processing device and program |
US9639952B2 (en) | 2012-12-12 | 2017-05-02 | Konica Minolta, Inc. | Image-processing apparatus and storage medium |
US9901317B2 (en) | 2013-05-31 | 2018-02-27 | Konica Minolta, Inc. | Image processing apparatus for acquiring a dynamic image and storage medium |
JP2018110637A (en) * | 2017-01-10 | 2018-07-19 | コニカミノルタ株式会社 | Dynamic image processing system |
US11238590B2 (en) | 2017-01-10 | 2022-02-01 | Konica Minolta, Inc. | Dynamic image processing apparatus |
JP2020171475A (en) * | 2019-04-10 | 2020-10-22 | コニカミノルタ株式会社 | Dynamic image analysis apparatus, dynamic image analysis method, and program |
JP7078098B1 (en) | 2020-12-14 | 2022-05-31 | コニカミノルタ株式会社 | Mobile radiography equipment, programs, methods, controls and radiography systems |
JP2022093888A (en) * | 2020-12-14 | 2022-06-24 | コニカミノルタ株式会社 | Movable radiographic apparatus, program, method, control device, and radiographic system |
Also Published As
Publication number | Publication date |
---|---|
JP5625800B2 (en) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2012110400A (en) | Dynamic diagnosis support information generation system | |
JP5919717B2 (en) | Dynamic medical image generation system | |
JP5672147B2 (en) | Chest diagnosis support information generation system | |
JP5874636B2 (en) | Diagnosis support system and program | |
JP5948950B2 (en) | Dynamic analysis system | |
JP5617559B2 (en) | Radiation imaging system | |
JPWO2012026146A1 (en) | Chest diagnosis support system and program | |
JP2009153677A (en) | Kinetic image processing system | |
JP6540807B2 (en) | Shooting console | |
JP5625800B2 (en) | Dynamic diagnosis support information generation system | |
JP6507581B2 (en) | Chest dynamic imaging support system | |
JP5625799B2 (en) | Dynamic diagnosis support information generation system | |
JP5617577B2 (en) | Chest diagnosis support information generation method | |
JP5772653B2 (en) | Dynamic diagnosis support information generation system, dynamic diagnosis support information generation method, and dynamic analysis device | |
JP5125750B2 (en) | Dynamic image diagnosis support system | |
JP2014147844A (en) | Dynamic diagnosis support information generation system | |
JP2012110398A (en) | Diagnosis support information generation system | |
JP5772523B2 (en) | Medical image photographing system, medical image processing apparatus, and program | |
JP2020000475A (en) | Dynamic image processing apparatus and program | |
JP2012115582A (en) | Thoracic part diagnosis aiding system | |
JP6950483B2 (en) | Dynamic photography system | |
JP2020089612A (en) | Image display device, image display method, and image display program | |
JP2009273603A (en) | Dynamic image capturing system | |
JP2017217047A (en) | Image display system | |
JP2003290202A (en) | Radiographing apparatus, radiographing method, system for radiographic image, program, and computer readable storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20130416 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20131008 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140318 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140603 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140707 |
|
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: 20140902 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140915 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5625800 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |