WO2009107658A1 - Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置 - Google Patents

Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置 Download PDF

Info

Publication number
WO2009107658A1
WO2009107658A1 PCT/JP2009/053407 JP2009053407W WO2009107658A1 WO 2009107658 A1 WO2009107658 A1 WO 2009107658A1 JP 2009053407 W JP2009053407 W JP 2009053407W WO 2009107658 A1 WO2009107658 A1 WO 2009107658A1
Authority
WO
WIPO (PCT)
Prior art keywords
ray
image
lesion
simulated
simulated lesion
Prior art date
Application number
PCT/JP2009/053407
Other languages
English (en)
French (fr)
Inventor
浩一 廣川
大雅 後藤
嘉晃 菅谷
Original Assignee
株式会社 日立メディコ
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社 日立メディコ filed Critical 株式会社 日立メディコ
Priority to US12/866,484 priority Critical patent/US8270560B2/en
Priority to JP2010500716A priority patent/JP5588859B2/ja
Publication of WO2009107658A1 publication Critical patent/WO2009107658A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms

Definitions

  • the present invention relates to an X-ray computed tomography (CT), and more particularly to a technique for obtaining an image corresponding to imaging conditions by simulation.
  • CT computed tomography
  • X-ray CT systems have made it possible to collect a wide range of data in a short period of time by increasing the number of X-ray detectors and increasing the scanner rotation speed.
  • the exposure dose per examination is increasing, and reduction of exposure by low-dose imaging is regarded as important.
  • simple low-dose imaging alone causes an increase in image noise, and an image suitable for diagnosis may not be acquired.
  • Patent Document 1 a technique for simulating an image obtained under an arbitrary shooting condition designated by an operator before shooting has been proposed. According to this technique, an image obtained under the designated photographing condition can be presented to the operator before photographing. Therefore, the operator can determine whether or not an image suitable for diagnosis can be acquired according to the designated photographing condition.
  • An object of the present invention is to provide an X-ray CT scan simulator and an X-ray CT apparatus capable of displaying a lesion assumed by an operator on a simulation image and capable of highly accurate simulation including the assumption of a new lesion occurrence. There is.
  • an X-ray CT scan simulator assumes an image storage device that stores a reference image, a target noise value setting device that sets a noise target value of a desired image, and an operator.
  • a simulated lesion setting device for setting a simulated lesion condition; a simulation image generating device for generating a simulation image including a simulated lesion using the reference image based on the set target noise value and simulated lesion condition; and And an image display device for displaying a simulation image including a simulated lesion.
  • the X-ray CT apparatus an X-ray source that irradiates a subject with X-rays, an X-ray detector that is disposed opposite to the X-ray source and detects X-rays transmitted through the subject, A tomographic image of the subject is reproduced based on a rotating device that is mounted with the X-ray source and the X-ray detector and rotates around the subject, and transmitted X-ray doses in a plurality of directions detected by the X-ray detector.
  • an X-ray CT apparatus comprising: an image reconstructing apparatus to be configured; an imaging condition input apparatus for inputting the X-ray irradiation conditions and image reconstruction conditions; and an image display apparatus for displaying the tomographic image.
  • a target noise value setting device for setting a noise target value of a desired image
  • a simulated lesion setting device for setting conditions of a simulated lesion assumed by an operator
  • the set target noise value Include simulated lesions using the reference image based on simulated lesion conditions
  • simulation image generating device for generating a simulation image, further wherein the image display device and displaying a simulation image including the simulated lesion.
  • the lesion assumed by the operator can be displayed on the simulation image, and a highly accurate simulation including the assumption of occurrence of a new lesion is possible.
  • Configuration diagram of the X-ray CT apparatus of the present invention Configuration diagram of the X-ray CT scan simulator of the present invention The figure which shows the processing flow which shows the method of creating the simulation image The figure which shows an example of the condition setting screen for the simulated lesion Outline explanatory diagram of approximate model calculation method Explanatory drawing of how to create reference projection data Explanatory diagram of relationship between distribution of slice thickness direction of simulated lesion and X-ray path
  • 1 X-ray CT device 100 scan gantry section, 101 X-ray tube, 102 rotating disk, 103 collimator, 104 opening, 105 bed, 106 X-ray detector, 107 data collection device, 108 gantry control device, 109 bed control device , 110 X-ray control device, 120 console, 121 input device, 122 image calculation device, 123 storage device, 124 system control device, 125 display device, 2 X-ray CT scan simulator, 200 simulation image generation device, 201 object model Calculation means, 202 image noise calculation means, 203 addition noise amount calculation means, 204 reference projection data creation means, 205 noise projection data creation means, 206 simulated lesion projection data creation means, 207 projection data addition means, 208 image reconstruction means, 209 simulation image creation means, 210 image storage device, 220 condition setting device, 230 simulation image display device, 400 Simulated lesion condition setting screen, 410 reference image display area, 411 reference image, 412 simulated lesion coordinate setting cursor, 420 simulated lesion condition setting area,
  • FIG. 1 is an overall configuration diagram of an X-ray CT apparatus 1 to which the present invention is applied.
  • the X-ray CT apparatus 1 includes a scan gantry unit 100 and a console 120.
  • the scan gantry unit 100 includes an X-ray tube 101, a rotating disk 102, a collimator 103, an X-ray detector 106, a data collection device 107, a bed 105, a gantry control device 108, a bed control device 109, An X-ray control device 110.
  • the X-ray tube 101 is an apparatus that irradiates a subject placed on a bed 105 with X-rays.
  • the collimator 103 is a device that controls the radiation direction of X-rays emitted from the X-ray tube 101.
  • the X-ray detector 106 is a device that is arranged opposite to the X-ray tube 101 and detects X-rays transmitted through the subject, and generates an electrical signal corresponding to the transmitted X-ray dose.
  • the rotating disk 102 includes an opening 104 into which a subject placed on a bed 105 enters, and is equipped with an X-ray tube 101 and an X-ray detector 106, and rotates around the subject.
  • the data collection device 107 is a device that converts the X-rays detected by the X-ray detector 106 into a predetermined signal.
  • the gantry control device 108 is a device that controls the rotation of the rotary disk 102.
  • the bed control device 109 is a device that controls the vertical movement of the bed 105.
  • the X-ray control device 110 is a device that controls the output to the X-ray tube 101.
  • the console 120 includes an input device 121, an image calculation device 122, a display device 125, a storage device 123, and a system control device 124.
  • the input device 121 is a device for inputting a subject's name, examination date / time, imaging conditions, simulated lesion conditions, and the like, and is specifically a keyboard, a pointing device, or the like.
  • the image computation device 122 is a device that performs CT processing on the measurement data sent from the data collection device 107 and performs CT image reconstruction.
  • the display device 125 is a device that displays the CT image created by the image calculation device 122.
  • the storage device 123 is a device that stores data collected by the data collection device 107 and image data of a CT image created by the image calculation device 122.
  • the system control device 124 is a device that controls these devices, the gantry control device 108, the bed control device 109, and the X-ray control device 110.
  • the X-ray tube 101 is controlled by the X-ray control device 110, and emits X-rays based on the imaging conditions (X-ray tube voltage, X-ray tube current, etc.) input from the input device 121.
  • the X-ray detector 106 has a plurality of X-ray detection elements arranged in the circumferential direction of the rotating disk 102 (for example, 1000), or two-dimensionally in the circumferential direction of the rotating disk 102 and the rotational axis direction of the rotating disk 102. The X-rays irradiated from the X-ray tube 101 and transmitted through the subject are detected by each detection element.
  • the rotating disk 102 is controlled by the gantry control device 108 and rotates based on imaging conditions (scanning speed, etc.) input from the input device 121.
  • the bed 105 is controlled by the bed control device 109 and operates based on the imaging conditions (such as a helical pitch) input from the input device 121.
  • X-ray projection data from various directions is collected by the data collection device 107 by irradiating and detecting X-rays while rotating the rotating disk 102 around the subject.
  • the X-ray projection data collected by the data collection device 107 is sent to the image calculation device 122.
  • the image calculation device 122 reconstructs the X-ray projection data to obtain a CT image.
  • the reconstructed CT image is displayed on the display device 125, and is stored in the storage device 123 as image data together with the imaging conditions.
  • the X-ray CT scan simulator 2 Prior to imaging with the X-ray CT apparatus 1, a configuration diagram of an X-ray CT scan simulator 2 that simulates an image obtained by imaging is shown in FIG.
  • the X-ray CT scan simulator 2 includes an image storage device 210, a condition setting device 220, a simulation image generation device 200, and a simulation image display device 230.
  • the image storage device 210 is a device that stores a reference image serving as a reference for a simulation image created by the X-ray CT scan simulator 2, and is specifically an HD (Hard Disk).
  • the reference image is a past image obtained by CT imaging of the subject in the past or a human phantom image obtained by CT imaging of a human phantom that faithfully reproduces the internal tissue of the human body.
  • Sufficient X-ray dose should be used to capture human phantom images so that clear images can be obtained.
  • the human phantom image may be prepared by taking a representative FOV for each part.
  • the storage device 123 may be used as the image storage device 210.
  • the condition setting device 220 is a device for an operator to set a target noise value and a simulated lesion condition, and specifically, a keyboard, a pointing device, and the like.
  • the target noise value is the image SD (Standard Deviation) value input by the operator, or the imaging conditions (tube voltage, tube current, slice thickness, scan time, etc.) set by the operator in the X-ray CT apparatus and the subject to be described later
  • the image SD value calculated based on the approximate model may be used.
  • Examples of the conditions for the simulated lesion include the coordinates, shape, and size of the simulated lesion, and the difference in CT value from the surrounding tissue.
  • the input device 121 may be used as the condition setting device 220.
  • the simulation image generation device 200 is a device that generates a simulation image based on the reference image, the target noise value, and the condition of the simulated lesion.
  • the simulation image generation device 200 is a CPU that executes arithmetic processing described later, or a dedicated arithmetic circuit.
  • the simulation image generating apparatus 200 includes an object approximate model calculating unit 201, an image noise calculating unit 202, an added noise amount calculating unit 203, a reference projection data generating unit 204, a noise projection data generating unit 205, and a simulated lesion projection.
  • Data creation means 206, projection data addition means 207, image reconstruction means 208, and simulation image creation means 209 are provided.
  • the simulation image generation apparatus 200 may be mounted on the image calculation apparatus 122 or the system control apparatus 124 of the X-ray CT apparatus 1, or may be mounted on another apparatus other than the X-ray CT apparatus 1, and may be connected to a LAN (Local Area Area). Communication with the X-ray CT apparatus 1 may be performed via a network.
  • LAN Local Area Area
  • the subject approximate model calculation means 201 calculates the approximate transmission length of the homogeneous medium (for example, water) in the relevant cross section of the subject from the reference image, and calculates the subject approximate model having the calculated approximate transmission length. A method for calculating the approximate object model will be described later.
  • the homogeneous medium for example, water
  • the image noise calculation unit 202 calculates an image noise value included in the reference image. A method for calculating the image noise will be described later.
  • the addition noise amount calculation means 203 calculates the amount of noise to be added to the reference image in order to create an image that realizes a desired image noise value. A method for calculating the amount of noise to be added will be described later.
  • the reference projection data creation unit 204 creates reference projection data based on the reference image. A method for calculating the reference projection data will be described later.
  • the noise projection data creation means 205 creates noise projection data for a noise image to be added to the reference image while taking into account the value of the reference projection data. A method for calculating the noise projection data will be described later.
  • the simulated lesion projection data creation means 206 creates simulated lesion projection data based on the simulated lesion conditions (coordinates, shape, size, CT value difference with surrounding tissue) set via the condition setting device 220. Is.
  • the simulated lesion projection data is created as projection data that does not include noise. A method for calculating simulated lesion projection data will be described later.
  • Projection data adding means 207 adds noise projection data and simulated lesion projection data to create addition image projection data.
  • the image reconstruction unit 208 applies an reconstruction process (for example, the Filtered Back Projection method) to the addition image projection data, and creates an addition image to be added to the reference image.
  • an reconstruction process for example, the Filtered Back Projection method
  • the image calculation device 122 may be used as the image reconstruction unit 208.
  • the simulation image creation means 209 creates a simulation image by adding the addition image to the reference image.
  • the created simulation image is an image that satisfies the image noise value set via the condition setting device 220 and the condition of the simulated lesion.
  • the simulation image display device 230 is a device that displays a simulation image generated by the simulation image generation device 200, and specifically, a CRT (Cathode Ray Tube) or a liquid crystal display. Note that the display device 125 may be used as the simulation image display device 230.
  • FIG. 3 is a diagram showing a processing flow until a simulation image is generated according to the present embodiment.
  • the general flow is to set a reference image in steps S301 to S303, generate noise data in steps S304 to S308, generate simulated lesion data in step S309, and generate a simulation image in steps S310 to S312. .
  • FIG. 4 shows an example of a simulated lesion condition setting screen.
  • the simulated lesion condition setting screen 400 includes a reference image display area 410 and a simulated lesion condition setting area 420, and is displayed on the image display device 230.
  • the selected reference image 411 is displayed in the reference image display area 410, and various boxes 421 to 425 used for setting various conditions of the simulated lesion are displayed in the simulated lesion condition setting area 420.
  • Various boxes include x-coordinate setting box 421 and y-coordinate setting box 422 for setting the coordinate of the simulated lesion, 423 for setting the size of the simulated lesion, and for setting the CT value difference between the simulated lesion and surrounding tissue
  • a box 424 and a simulated lesion shape setting box 425 are displayed. Boxes 421 to 424 are boxes for numerical editing, and box 425 is a box for menu selection. The operator can set the conditions of the desired simulated lesion by operating these boxes. Alternatively, the coordinates of the simulated lesion may be set by a method such as clicking the desired position in the reference image 411 with the mouse and setting the simulated lesion coordinate setting cursor 412.
  • Step S301 The X-ray CT scan simulator 2 searches the image storage device 210 for a past image of the subject to be simulated, and determines whether there is a history of imaging the same part in the past. If it is determined that there is a past image, the process proceeds to step S302; otherwise, the process proceeds to step S303.
  • Step S302 The X-ray CT scan simulator 2 sets the past image as the reference image.
  • X-ray CT scan simulator 2 sets the human phantom image as the reference image.
  • Step S304 The subject approximate model calculating unit 201 calculates the subject approximate model based on the reference image set in S302 or S303 as follows.
  • Figure 5 outlines the approximate model calculation method.
  • the CT values of the reference image are integrated in the y direction as shown in Equation 1 to obtain profile1.
  • vch Virtual detector channel number
  • vch 0,1, ..., XMTX-1
  • XMTX Number of pixels in the x direction
  • ypt Pixel pitch in the y direction
  • YMTX Number of pixels in the y direction
  • CT of the reference image Accumulate the values in the x direction as shown in Equation 2 to obtain profile2.
  • the transmission length ypd in the y direction and the transmission length xpd in the x direction are calculated, and these are used as the diameters in the y direction and the x direction of the elliptic approximation model.
  • the approximate model is assumed to be composed of water, but is not necessarily limited to water.
  • nmaxave (prf, n) a function for calculating a simple addition average of n pieces of data from the larger value of the data string prf (step S305)
  • the image noise calculation means 202 calculates the image noise amount of the subject approximate model calculated in S304 as follows. In order to calculate the amount of noise to be added to the reference image in order to realize the target image SD, first, it is necessary to calculate the image noise variance in the reference image, that is, the square value of the image SD in the reference image.
  • the added noise amount calculating means 203 calculates the amount of noise to be added to the reference image based on the image noise amount of the reference image calculated in S305 and a preset target noise value. Also, a virtual scan condition for creating noise projection data corresponding to the amount of noise to be added is calculated. The calculation method will be described below.
  • the target image SD is SDtgt
  • the target image noise variance Vtgt is as shown in Equation 4.
  • a noise image having the image noise variance Vn_img represented by Equation 5 may be created.
  • the tube current / time product of the reference image is xmAs0
  • the tube current / time product xmAs_ni of the noise image to be added to the reference image is as shown in Formula 6.
  • the reference projection data creation means 204 creates reference projection data as shown in FIG. 6 by reprojecting the reference image.
  • the reprojection processing integration similar to Equation 1 or Equation 2 is performed.
  • the X-ray tube is set to 0 °, and the angle of the position of the X-ray tube increases clockwise. For convenience of explanation, only the reprojection result in two directions is shown. Yes. In actual processing, it is desirable to perform re-projection at the X-ray tube position based on the same view as the assumed scan, or at the X-ray tube position based on the half-round view.
  • Step S308 The noise projection data creation unit 205 creates noise projection data as follows based on the virtual scan condition calculated in S306 and the reference projection data created in S307.
  • the tube current / time product xmAs_ni of the noise image is shown in the above equation (6).
  • the projection data value Vr_ref (prf) when the projection data value is prf and scanning is performed with a certain reference tube current / time product xmAs_ref can be obtained in advance by an experiment, for example, expressed by Equation 7 it can.
  • bsprf (vch) The value of the reference projection data for the virtual detector channel number vch, that is, the projection data indicated by Vr_ni (vch) in which the average value is 0 and the variance is several 8 in each channel of the virtual detector is the noise projection data
  • the simulated lesion projection data creating means 206 creates simulated lesion projection data based on preset simulated lesion conditions (coordinates, size, and CT value difference with surrounding tissue) as follows.
  • ⁇ w Water X-ray attenuation coefficient
  • the distribution of the simulated lesion in the image slice thickness direction is considered. That is, as shown in Fig. 7 (a), the X-ray path in the cross-sectional image viewed from the body axis direction is a line segment, but it is necessary to consider that it is actually a plane having a depth corresponding to the slice thickness. is there.
  • the diameter ⁇ of the intersection Cfs between the simulated lesion and the X-ray path p shown in FIGS. 7 (b) and 7 (c) is expressed by equation (10).
  • Step S310 The projection data adding means 207 adds the noise projection data created in S308 and the simulated lesion projection data created in S309, and creates projection data of the image for addition.
  • Step S311 The image reconstruction unit 208 reconstructs an image by applying a filtered back projection method or the like to the projection data of the image for addition created in S310, and creates an image for addition.
  • Step S312 The simulation image creation means 209 adds the addition image created in S311 to the reference image set in S302 or S303 to create a simulation image. If the simulation image created in S312 is satisfactory for the operator in terms of visual effects of image noise and discriminability of the simulated lesion, the actual scan conditions for achieving the preset target image noise value are captured.
  • the scanning conditions may be used.
  • the tube current / time product xmAs_app to be applied as the scanning condition for the main imaging is as shown in Equation 12.
  • the operator may input the result of Formula 12 displayed on the simulation image display device by the simulation image generation apparatus 200 using the input device 121.
  • the simulation image generation apparatus 200 may set the result of Formula 12 in the X-ray control apparatus 110.
  • the simulation image generation apparatus 200 sets the result of Formula 12 in the X-ray control apparatus 110, the operation of the operator can be simplified.
  • the above-described functions may be realized by configuring a scan simulator program that causes a computer to execute the functions described in the above-described embodiment, and installing the program in a personal computer or a workstation.
  • the simulator image does not necessarily have to be displayed on the monitor of the personal computer or workstation on which the scan simulator program is installed, and the simulator image is transmitted to a terminal device connected via a network such as a LAN.
  • a simulator image may be displayed.

Abstract

　操作者が想定する病巣をシミュレーション画像上に表示でき、新たな病巣発生の想定も含めた精度の高いシミュレーションが可能なX線CTスキャンシミュレータ、X線CT装置を提供する。 　基準画像を格納する画像格納装置と、所望画像のノイズ目標値を設定する目標ノイズ値設定装置と、操作者が想定する模擬病巣の条件を設定する模擬病巣設定装置と、前記設定された目標ノイズ値と模擬病巣の条件に基づき前記基準画像を用いて模擬病巣を含んだシミュレーション画像を生成するシミュレーション画像生成装置と、前記模擬病巣を含んだシミュレーション画像を表示する画像表示装置と、を備える。

Description

Ｘ線ＣＴスキャンシミュレータ及びＸ線ＣＴ装置
　本発明は、X線コンピュータ断層装置(CT)に関し、特に撮影条件に対応する画像をシミュレーションにより得る技術に関するものである。
　近年、X線CT装置では、X線検出器の多列化とスキャナ回転速度の高速化により、短時間で広範囲のデータ収集が可能になった。その一方で、一検査あたりの被曝線量が増加しており、低線量撮影による被曝低減が重要視されている。しかし、単なる低線量撮影だけでは画像ノイズの増加を招き、診断に適した画像を取得できない場合がある。
　そこで、特許文献1のように、操作者が指定する任意の撮影条件において得られる画像を、撮影前にシミュレーションする技術が提案されている。この技術によれば、指定された撮影条件により得られる画像を操作者へ撮影前に提示できるため、操作者は指定された撮影条件により診断に適した画像を取得できるか否かを判断できる。
特開2004-057831号公報 特開2006-116137号公報
　しかしながら、特許文献1に開示されている技術では、操作者が想定する病巣がシミュレーション画像上でどのように見えるのか、すなわち腫瘍の発生・転移等を想定したシミュレーションを行なうことについて配慮がなされていなかった。
　本発明の目的は、操作者が想定する病巣をシミュレーション画像上に表示でき、新たな病巣発生の想定も含めた精度の高いシミュレーションが可能なX線CTスキャンシミュレータ、X線CT装置を提供することに有る。
　上記課題を解決するために、本発明に係るX線CTスキャンシミュレータは、基準画像を格納する画像格納装置と、所望画像のノイズ目標値を設定する目標ノイズ値設定装置と、操作者が想定する模擬病巣の条件を設定する模擬病巣設定装置と、前記設定された目標ノイズ値と模擬病巣の条件に基づき前記基準画像を用いて模擬病巣を含んだシミュレーション画像を生成するシミュレーション画像生成装置と、前記模擬病巣を含んだシミュレーション画像を表示する画像表示装置と、を備えることを特徴とする。
　また、本発明に係わるX線CT装置は、被検体にX線を照射するX線源と、前記X線源に対向配置され前記被検体を透過したX線を検出するX線検出器と、前記X線源と前記X線検出器を搭載し前記被検体の周囲を回転する回転装置と、前記X線検出器により検出された複数方向の透過Ｘ線量に基づき前記被検体の断層像を再構成する画像再構成装置と、前記X線の照射条件と画像再構成の条件を入力する撮影条件入力装置と、前記断層像を表示する画像表示装置と、を有するX線CT装置において、基準画像を格納する画像格納装置と、所望画像のノイズ目標値を設定する目標ノイズ値設定装置と、操作者が想定する模擬病巣の条件を設定する模擬病巣設定装置と、前記設定された目標ノイズ値と模擬病巣の条件に基づき前記基準画像を用いて模擬病巣を含んだシミュレーション画像を生成するシミュレーション画像生成装置と、をさらに備え、前記画像表示装置は前記模擬病巣を含んだシミュレーション画像を表示することを特徴とする。
　本発明によれば、操作者が想定する病巣をシミュレーション画像上に表示でき、新たな病巣発生の想定も含めた精度の高いシミュレーションが可能となる。
本発明のX線CT装置の構成図 本発明のX線CTスキャンシミュレータの構成図 シミュレーション画像を作成する方法を示した処理フローを示す図 模擬病巣の条件設定用画面の一例を示す図 近似モデル算出方法の概要説明図 基準投影データ作成方法の説明図 模擬病巣のスライス厚さ方向の分布とX線経路との関係説明図
符号の説明
　1　X線CT装置、100　スキャンガントリ部、101　Ｘ線管、102　回転円盤、103　コリメータ、104　開口部、105　寝台、106　X線検出器、107　データ収集装置、108　ガントリ制御装置、109　寝台制御装置、110　X線制御装置、120　操作卓、121　入力装置、122　画像演算装置、123　記憶装置、124　システム制御装置、125　表示装置、2　X線CTスキャンシミュレータ、200　シミュレーション画像生成装置、201　被検体モデル算出手段、202　画像ノイズ算出手段、203　加算ノイズ量算出手段、204　基準投影データ作成手段、205　ノイズ投影データ作成手段、206　模擬病巣投影データ作成手段、207　投影データ加算手段、208　画像再構成手段、209　シミュレーション画像作成手段、210　画像格納装置、220　条件設定装置、230　シミュレーション画像表示装置、400　模擬病巣の条件設定用画面、410　基準画像表示エリア、411　基準画像、412　模擬病巣座標設定用カーソル、420　模擬病巣条件設定エリア、421　x座標設定用ボックス、422　y座標設定用ボックス、423　大きさ設定用ボックス、424　CT値差設定用ボックス、425　形状設定用ボックス
発明を実施するための形態
　以下、添付図面に従って本発明の好ましい実施の形態について詳説する。
　《ハードウェア構成》
　図1は本発明を適用したX線CT装置1の全体構成図である。X線CT装置1はスキャンガントリ部100と操作卓120とを備える。
　スキャンガントリ部100は、X線管101と、回転円盤102と、コリメータ103と、X線検出器106と、データ収集装置107と、寝台105と、ガントリ制御装置108と、寝台制御装置109と、X線制御装置110と、を備えている。X線管101は寝台105上に載った被検体にX線を照射する装置である。コリメータ103はX線管101から照射されるX線の放射方向を制御する装置である。X線検出器106は、X線管101と対向配置され被検体を透過したX線を検出する装置であり、透過X線量に応じた電気信号を発生させる。回転円盤102は、寝台105上に載った被検体が入る開口部104を備えるとともに、X線管101とX線検出器106を搭載し、被検体の周囲を回転するものである。データ収集装置107は、X線検出器106で検出されたX線を所定の信号に変換する装置である。ガントリ制御装置108は回転円盤102の回転を制御する装置である。寝台制御装置109は、寝台105の上下前後動を制御する装置である。X線制御装置110はX線管101への出力を制御する装置である。
　操作卓120は、入力装置121と、画像演算装置122と、表示装置125と、記憶装置123と、システム制御装置124とを備えている。入力装置121は、被検体氏名、検査日時、撮影条件、模擬病巣の条件などを入力するための装置であり、具体的にはキーボードやポインティングデバイス等である。画像演算装置122は、データ収集装置107から送出される計測データを演算処理してCT画像再構成を行う装置である。表示装置125は、画像演算装置122で作成されたCT画像を表示する装置である。記憶装置123は、データ収集装置107で収集したデータ及び画像演算装置122で作成されたCT画像の画像データを記憶する装置である。システム制御装置124は、これらの装置及びガントリ制御装置108と寝台制御装置109とX線制御装置110を制御する装置である。
　X線管101はX線制御装置110によって制御され、入力装置121から入力された撮影条件(X線管電圧やX線管電流など)に基づいたX線を照射する。X線検出器106はX線検出素子を回転円盤102の円周方向に多数(例えば1000個)配列したもの、若しくは回転円盤102の円周方向と回転円盤102の回転軸方向とに2次元的に配列したもので、X線管101から照射され被検体を透過したＸ線を各検出素子で検出する。回転円盤102はガントリ制御装置108により制御され、入力装置121から入力された撮影条件(スキャン速度など)に基づいて回転する。寝台105は寝台制御装置109によって制御され、入力装置121から入力された撮影条件(らせんピッチなど)に基づいて動作する。
　回転円盤102を被検体の周りで回転させながらX線の照射と検出をすることにより、様々な方向(例えば1000方向)からのX線投影データがデータ収集装置107で収集される。データ収集装置107で収集されたX線投影データは、画像演算装置122へ送出される。画像演算装置122は、X線投影データを再構成演算してCT画像とする。再構成されたCT画像は表示装置125に表示され、また撮影条件とともに画像データとして記憶装置123に記憶される。
　《X線CTスキャンシミュレータの構成》
　X線CT装置1での撮影に先立ち、撮影により得られる画像をシミュレートするX線CTスキャンシミュレータ2の構成図を図2に示す。X線CTスキャンシミュレータ2は、画像格納装置210と、条件設定装置220と、シミュレーション画像生成装置200と、シミュレーション画像表示装置230を備える。
　画像格納装置210は、X線CTスキャンシミュレータ2が作成するシミュレーション画像の基準となる基準画像を格納する装置であり、具体的にはHD(Hard Disk)などである。基準画像は、被検体を過去にCT撮影して得られた過去画像、または人体の内部組織を忠実に再現した人体ファントムを予めCT撮影して得られた人体ファントム画像である。人体ファントム画像の撮影には充分なX線量を用い、鮮明な画像が得られるようにする。また人体ファントム画像には、部位ごとに代表的なFOVで撮影したものを用意しておけばよい。なお、画像格納装置210として、記憶装置123を使用しても良い。
　条件設定装置220は、目標ノイズ値や模擬病巣の条件の設定を操作者が行うための装置であり、具体的にはキーボードやポインティングデバイス等である。目標ノイズ値は、操作者が入力する画像SD(Standard Deviation)値、若しくは操作者がX線CT装置において設定した撮影条件(管電圧、管電流、スライス厚、スキャン時間等)と後述する被検体近似モデルに基づき算出される画像SD値でもよい。模擬病巣の条件としては、模擬病巣の座標、形状、大きさ、周囲組織とのCT値の差等が挙げられる。なお、条件設定装置220として、入力装置121を使用しても良い。
　シミュレーション画像生成装置200は、基準画像、目標ノイズ値、及び模擬病巣の条件に基づきシミュレーション画像を生成する装置であり、具体的には後述する演算処理を実行するCPU、若しくは専用の演算回路である。シミュレーション画像生成装置200は、被検体近似モデル算出手段201と、画像ノイズ算出手段202と、加算ノイズ量算出手段203と、基準投影データ作成手段204と、ノイズ投影データ作成手段205と、模擬病巣投影データ作成手段206と、投影データ加算手段207と、画像再構成手段208と、シミュレーション画像作成手段209を備える。なお、シミュレーション画像生成装置200は、X線CT装置1の画像演算装置122やシステム制御装置124に搭載されてもよいし、X線CT装置1以外の他の装置に搭載され、LAN(Local Area Network)を介してX線CT装置1と通信してもよい。
　被検体近似モデル算出手段201は、基準画像から被検体の該当横断面における均質媒体(例えば水)の近似透過長を求め、求めた近似透過長を持つ被検体近似モデルを算出するものである。被検体近似モデルの算出方法は後述する。
　画像ノイズ算出手段202は、基準画像が有する画像ノイズ値を算出するものである。画像ノイズの算出方法は後述する。
　加算ノイズ量算出手段203は、所望の画像ノイズ値を実現した画像を作成するために基準画像に加算すべきノイズ量を算出するものである。加算すべきノイズ量の算出方法は後述する。
　基準投影データ作成手段204は、基準画像に基づき基準投影データを作成するものである。基準投影データの算出方法は後述する。
　ノイズ投影データ作成手段205は、基準投影データの値を考慮しつつ基準画像に加算すべきノイズ画像のためのノイズ投影データを作成するものである。ノイズ投影データの算出方法は後述する。
　模擬病巣投影データ作成手段206は、条件設定装置220を介して設定された模擬病巣の条件(座標、形状、大きさ、周囲組織とのCT値の差)に基づき模擬病巣の投影データを作成するものである。模擬病巣投影データは、ノイズを含まない投影データとして作成される。模擬病巣投影データの算出方法は後述する。
　投影データ加算手段207は、ノイズ投影データと模擬病巣投影データを加算して加算画像用投影データを作成するものである。
　画像再構成手段208は、加算画像用投影データを再構成処理(例えば、Filtered Back Projection法)を適用して、基準画像に加算する加算用画像を作成するものである。なお、画像再構成手段208として画像演算装置122を使用しても良い。
　シミュレーション画像作成手段209は、基準画像に加算用画像を加算してシミュレーション画像を作成するものである。作成されたシミュレーション画像は、条件設定装置220を介して設定された画像ノイズ値と模擬病巣の条件を満たした画像となる。
　シミュレーション画像表示装置230は、シミュレーション画像生成装置200により生成されたシミュレーション画像を表示する装置であり、具体的にはCRT(Cathode Ray Tube)や液晶ディスプレイ等である。なお、シミュレーション画像表示装置230として、表示装置125を使用しても良い。
　《処理の流れ》
　図3は、本実施形態によりシミュレーション画像を生成するまでの処理の流れを示す図である。大まかな流れは、ステップS301～S303での基準画像の設定、ステップS304～S308でのノイズデータの生成、ステップS309での模擬病巣データの生成、ステップS310～S312でのシミュレーション画像の生成、である。
　なお、目標ノイズ値及び模擬病巣の条件(座標、大きさ、周囲組織とのCT値の差)は、条件設定装置220を介して予め設定されているものとする。図４に模擬病巣の条件設定用画面の一例を示す。模擬病巣の条件設定用画面400は、基準画像表示エリア410と模擬病巣条件設定エリア420から構成され、画像表示装置230に表示される。基準画像表示エリア410には選択された基準画像411が表示され、模擬病巣条件設定エリア420には模擬病巣の各種条件を設定するのに用いる各種ボックス421～425が表示される。各種ボックスとしては、模擬病巣の座標設定用のx座標設定用ボックス421とy座標設定用ボックス422と、模擬病巣の大きさ設定用ボックス423と、模擬病巣と周囲組織とのCT値差設定用ボックス424と、模擬病巣の形状設定用ボックス425が表示される。ボックス421～424は数値編集用のボックスであり、ボックス425はメニュー選択用のボックスである。操作者はこれらのボックス操作することにより、所望の模擬病巣の条件を設定できる。また基準画像411内の所望の位置をマウスでクリックして、模擬病巣座標設定用カーソル412を設定する等の方法により、模擬病巣の座標を設定しても良い。
　図3の各ステップについて以下で説明する。
　(ステップS301)
　X線CTスキャンシミュレータ2は、シミュレーション対象となる被検体の過去画像を画像格納装置210の中で検索し、過去に同一部位を撮影した履歴があるか否かを判断する。判断の結果、過去画像があればステップS302に、なければステップS303に進む。
　(ステップS302)
　X線CTスキャンシミュレータ2は、過去画像を基準画像に設定する。
　(ステップS303)。
　X線CTスキャンシミュレータ2は、人体ファントム画像を基準画像に設定する。
　(ステップS304)
　被検体近似モデル算出手段201は、S302若しくはS303で設定された基準画像に基づき、以下のように被検体近似モデルを算出する。
　図5に近似モデル算出方法の概要を示す。まず、基準画像のCT値を数1式のようにｙ方向に積算し、profile1を得る。
Figure JPOXMLDOC01-appb-M000001
 　　ただし、
　　vch：仮想的な検出器チャネル番号で、vch＝0,1,...,XMTX-1
　　XMTX：x方向の画素数
　　ypt：y方向の画素ピッチ
　　YMTX：y方向の画素数
　　org_img(vch,iy)：x座標=vch、y座標＝iyにおける基準画像のCT値
　次に、基準画像のCT値を数2式のようにx方向に積算し、profile2を得る。
Figure JPOXMLDOC01-appb-M000002
 　　ただし、
　　vch：仮想的な検出器チャネル番号で、vch＝0,1,...,YMTX-1
　　YMTX：y方向の画素数
　　xpt：x方向の画素ピッチ
　　XMTX：x方向の画素数
　　org_img(ix,vch)：x座標=ix、y座標＝vchにおける基準画像のCT値
　次に、数3式によりy方向の透過長ypdとx方向の透過長xpdを算出し、これらを楕円近似モデルのy方向とx方向の径とする。なお、本実施形態では近似モデルが水で構成されるものとしているが、必ずしも水でなくてもよい。
Figure JPOXMLDOC01-appb-M000003
 　　ただし、
　　nmaxave(prf,n)：データ列prfの値の大きい方からn個のデータによる単純加算平均を算出する関数
　(ステップS305)
　画像ノイズ算出手段202は、S304で算出された被検体近似モデルの画像ノイズ量を以下のように算出する。目標とする画像SDを実現するために基準画像に加算するノイズ量を算出するためには、まず基準画像における画像ノイズ分散、すなわち基準画像における画像SDの二乗値を算出する必要がある。そこで前記の楕円近似モデルに基づき、画像再構成に使用する各ビューの透過長や、X線管電圧、X線管電流、再構成時のビュー重み付け等を考慮して、たとえば特許文献2に開示された方法により、基準画像における画像ノイズ分散Vorg_imgを算出する。
　(ステップS306)
　加算ノイズ量算出手段203は、S305で算出した基準画像の画像ノイズ量と予め設定された目標ノイズ値に基づき、基準画像に加算すべきノイズ量を算出する。また、加算すべきノイズ量に対応するノイズ投影データ作成用の仮想的なスキャン条件も算出する。算出方法について以下で述べる。
　目標とする画像SDをSDtgtとすると、目標とする画像ノイズ分散Vtgtは数4式のとおりである。
Figure JPOXMLDOC01-appb-M000004
 　したがって、目標とする画像ノイズを実現するためには、数5式で示す画像ノイズ分散Vn_imgとなるノイズ画像を作成すればよい。
Figure JPOXMLDOC01-appb-M000005
 　また、基準画像の管電流・時間積をxmAs0とすれば、基準画像に加算するノイズ画像の
管電流・時間積xmAs_niは数6式で示すとおりとなる。
Figure JPOXMLDOC01-appb-M000006
 　(ステップS307)
　基準投影データ作成手段204は、基準画像を再投影処理することにより図6に示すような基準投影データを作成する。再投影処理には数1式や数2式と同様な積算を行う。なお、図6ではX線管が真上の状態を0°とし、X線管の位置をしめす角度は右回りに増大するものとしてあり、説明の便宜上、2方向の再投影結果のみを示している。実際の処理では想定するスキャンと同じビューによるＸ線管位置、または半周分のビューによるX線管位置において再投影を行うことが望ましい。
　(ステップS308)
　ノイズ投影データ作成手段205は、S306で算出された仮想的なスキャン条件とS307で作成した基準投影データに基づき、ノイズ投影データを以下のように作成する。　
　ノイズ画像の管電流・時間積xmAs_niは前記の数6式に示されている。また、投影データの値がprfであり、ある基準管電流・時間積xmAs_refでスキャンした場合の投影データノイズVr_ref(prf)は実験により予め求めておくことができ、たとえば数7式で表すことができる。
　　ただし、
　　cおよびa：予め実験によって求めておく定数
　　exp(x)：引数xに対し、自然対数の底eのx乗を算出する指数関数
　　prf：投影データの値
　よって、ある仮想チャネルvchにおける基準投影データの値に対応するノイズのみの投影データの分散値Vr_ni(vch)は数8式によって算出することができる。
Figure JPOXMLDOC01-appb-M000008
 　　ただし、
　　bsprf(vch)：仮想検出器チャネル番号vchにおける基準投影データの値
　すなわち、仮想検出器の各チャネルにおいて平均値が0で分散が数8式のVr_ni(vch)で示される投影データをノイズ投影データとする。
　このように基準画像を再投影して作成した基準投影データの値を考慮することにより、基準画像内での人体組織の分布の影響を考慮したノイズ投影データを作成することができる。
　(ステップS309)
　模擬病巣投影データ作成手段206は、予め設定された模擬病巣の条件(座標、大きさ、周囲組織とのCT値の差)に基づき、模擬病巣の投影データを以下のように作成する。
　図4に例示した模擬病巣の条件設定用画面等を介して、模擬病巣の中心座標、直径Dtum、模擬病巣と周囲組織とのCT値の差Htumが設定されているものとする。模擬病巣と周囲組織とのCT値の差Htumから、模擬病巣と周囲組織とのX線減弱係数差μcを数9式によって求める。
Figure JPOXMLDOC01-appb-M000009
 ただし、
μw：水のX線減弱係数
　また、模擬病巣の中心から距離dだけ隔たったX線経路における投影データを算出するため、模擬病巣の画像スライス厚方向の分布を考慮する。すなわち、図7(a)のように体軸方向からみた横断面画像でのX線経路は線分であるが、実際にはスライス厚さ分の奥行きを持つ面であることを考慮する必要がある。図7(b)及び(c)に示す模擬病巣とX線経路pとの交面Cfsの直径ξは数10式によって表される。
Figure JPOXMLDOC01-appb-M000010
 　図7(b)に示すように模擬病巣とX線経路との交面の直径ξが画像スライス厚Thkより小さい場合、もしくは図示しないがξ＝Thkである場合と、図7(c)に示すようにξ＞Thkである場合を考慮して、X線経路pにおける模擬病巣投影データの値prftum(p)は数11式によって表される。
Figure JPOXMLDOC01-appb-M000011
 　　ただし、
　　g：スケール合わせの定数。予め求めておく。
　(ステップS310)
　投影データ加算手段207は、S308で作成されたノイズ投影データとS309で作成された模擬病巣投影データを加算し、加算用画像の投影データを作成する。
　(ステップS311)
　画像再構成手段208は、S310で作成された加算用画像の投影データにFiltered Back Projection法等を適用して画像再構成し、加算用画像を作成する。
　(ステップS312)
　シミュレーション画像作成手段209は、S302若しくはS303で設定された基準画像に、S311で作成された加算用画像を加算し、シミュレーション画像を作成する。　
　S312で作成されたシミュレーション画像が画像ノイズの視覚的効果および模擬病巣の識別性の点で操作者にとって満足のいくものであれば、予め設定された目標画像ノイズ値を実現するスキャン条件を本撮影のスキャン条件として用いればよい。ここで、本撮影のスキャン条件として適用すべき管電流・時間積xmAs_appは数12式のとおりとなる。
Figure JPOXMLDOC01-appb-M000012
 　数12式の結果を適用して本撮影することにより診断に適した画像を取得することが可能となる。数12式の結果を本撮影に適用するには、シミュレーション画像生成装置200がシミュレーション画像表示装置に表示させた数12式の結果を操作者が入力装置121を用いて入力しても良いし、シミュレーション画像生成装置200が数12式の結果をX線制御装置110に設定しても良い。シミュレーション画像生成装置200が数12式の結果をX線制御装置110に設定した場合、操作者の操作を簡略化することができる。
　上記実施形態に記載の機能をコンピュータに実行させるスキャンシミュレータプログラムとして構成し、このプログラムをパーソナルコンピュータやワークステーションにインストールして上記機能を実現してもよい。その場合、シミュレータ画像は、必ずしもスキャンシミュレータプログラムをインストールしたパーソナルコンピュータやワークステーションのモニタに表示させる必要はなく、LANなどのネットワークを介して接続された端末装置にシミュレータ画像を送信し、その端末装置においてシミュレータ画像を表示させてもよい。このようにすることにより、例えば、X線CT装置の操作者である放射線技師が、X線CT装置の操作卓から離れた場所にいる放射線医に、撮影条件に対する意見をシミュレータ画像に基づきながら仰ぐことが可能になる。その結果、画像診断を行う放射線医の要望にそった画像が撮影されるようになり、より適切な画像診断がなされる。

Claims (7)

  1. 　基準画像を格納する画像格納装置と、
    　所望画像のノイズ目標値を設定する目標ノイズ値設定装置と、
    　操作者が想定する模擬病巣の条件を設定する模擬病巣設定装置と、
    　前記設定された目標ノイズ値と模擬病巣の条件に基づき前記基準画像を用いて模擬病巣を含んだシミュレーション画像を生成するシミュレーション画像生成装置と、
    　前記模擬病巣を含んだシミュレーション画像を表示する画像表示装置と、
    　を備えることを特徴とするX線CTスキャンシミュレータ。
  2. 　請求項1に記載のX線CTスキャンシミュレータにおいて、
    　前記模擬病巣設定装置で設定する病巣の条件として、病巣の位置、病巣の大きさ、病巣と周囲組織とのコントラストを含むことを特徴とするX線CTスキャンシミュレータ。
  3. 　請求項2に記載のX線CTスキャンシミュレータにおいて、
    　前記模擬病巣設定装置は、前記基準画像上で指定された点を前記病巣の位置として取得することを特徴とするX線CTスキャンシミュレータ。
  4. 　被検体にX線を照射するX線源と、前記X線源に対向配置され前記被検体を透過したX線を検出するX線検出器と、前記X線源と前記X線検出器を搭載し前記被検体の周囲を回転する回転装置と、前記X線検出器により検出された複数方向の透過X線量に基づき前記被検体の断層像を再構成する画像再構成装置と、前記X線の照射条件と画像再構成の条件を入力する操作に用いられる撮影条件入力装置と、前記断層像を表示する画像表示装置と、を有するX線CT装置において、
    　基準画像を格納する画像格納装置と、
    　所望画像のノイズ目標値を設定する目標ノイズ値設定装置と、
    　操作者が想定する模擬病巣の条件を設定する模擬病巣設定装置と、
    　前記設定された目標ノイズ値と模擬病巣の条件に基づき前記基準画像を用いて模擬病巣を含んだシミュレーション画像を生成するシミュレーション画像生成装置と、
    　をさらに備え、前記画像表示装置は前記模擬病巣を含んだシミュレーション画像を表示することを特徴とするX線CT装置。
  5. 　請求項4に記載のX線CT装置において、
    　前記シミュレーション画像生成装置が前記シミュレーション画像を生成する際に算出したX線照射条件に基づいて前記X線源を制御するX線制御装置をさらに備えることを特徴とするX線CT装置。
  6. 請求項4に記載のX線CT装置において、
    　前記模擬病巣設定装置で設定する病巣の条件として、病巣の位置、病巣の大きさ、病巣と周囲組織とのコントラストを含むことを特徴とするX線CT装置。
  7. 　請求項6に記載のX線CT装置において、
    　前記模擬病巣設定装置は、前記基準画像上で指定された点を前記病巣の位置として取得することを特徴とするX線CT装置。
PCT/JP2009/053407 2008-02-25 2009-02-25 Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置 WO2009107658A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US12/866,484 US8270560B2 (en) 2008-02-25 2009-02-25 X-ray CT scan simulator and X-ray CT apparatus
JP2010500716A JP5588859B2 (ja) 2008-02-25 2009-02-25 Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2008042308 2008-02-25
JP2008-042308 2008-02-25

Publications (1)

Publication Number Publication Date
WO2009107658A1 true WO2009107658A1 (ja) 2009-09-03

Family

ID=41016046

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/053407 WO2009107658A1 (ja) 2008-02-25 2009-02-25 Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置

Country Status (3)

Country Link
US (1) US8270560B2 (ja)
JP (1) JP5588859B2 (ja)
WO (1) WO2009107658A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017029461A (ja) * 2015-08-03 2017-02-09 株式会社日立製作所 医用画像処理装置および医用画像処理方法
JP2020519863A (ja) * 2017-05-01 2020-07-02 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 定量的分子撮像のための正確なハイブリッドデータセットの生成
KR20200104997A (ko) * 2019-02-28 2020-09-07 (주)토탈소프트뱅크 심혈관 형상변형 시뮬레이션 시스템

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2521489A2 (en) * 2010-01-06 2012-11-14 Koninklijke Philips Electronics N.V. Method for simulating reduction of acquisition dosage of an x-ray system, computer system and x-ray system
US9492136B2 (en) 2012-01-06 2016-11-15 Indiana University Research And Technology Corporation Method and apparatus that automates tube current and voltage selection for CT scans
US9050056B2 (en) 2012-12-26 2015-06-09 Biosense Webster (Israel) Ltd. Reduced X-ray exposure by simulating images
CN106233357B (zh) * 2013-10-07 2020-09-04 曼帝斯有限公司 基于医疗程序模拟的辐射估算和防护
JP2017051395A (ja) * 2015-09-09 2017-03-16 富士フイルム株式会社 放射線画像処理装置、方法およびプログラム
AU2018225165B2 (en) * 2017-02-24 2023-07-06 Bayer Healthcare Llc Systems and methods for generating simulated computed tomography (CT) images
EP3384847B1 (de) * 2017-04-05 2019-06-12 Siemens Healthcare GmbH Bestimmen eines referenzdosisparameters einer computertomographiebildgebung
CN108511043B (zh) * 2018-02-27 2022-06-03 华东师范大学 基于数值模拟的x-ct虚拟数据采集及图像重建方法及系统
US10896503B2 (en) 2018-03-23 2021-01-19 International Business Machines Corporation Identification of areas of interest in imaging applications
JP7282487B2 (ja) * 2018-06-07 2023-05-29 キヤノンメディカルシステムズ株式会社 医用画像診断装置
EP4092684A1 (en) * 2021-05-18 2022-11-23 Koninklijke Philips N.V. Method for providing feedback data in a medical imaging system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007074772A1 (ja) * 2005-12-28 2007-07-05 Hitachi Medical Corporation Ｘ線ｃｔ装置
WO2007114470A1 (ja) * 2006-04-04 2007-10-11 Hitachi Medical Corporation Ｘ線ｃｔスキャンシミュレータ装置、ｘ線ｃｔ装置、及びｘ線ｃｔスキャンシミュレータプログラム
JP2008048880A (ja) * 2006-08-24 2008-03-06 Ge Medical Systems Global Technology Co Llc 画像作成方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3732568B2 (ja) * 1996-04-03 2006-01-05 株式会社東芝 Ｘ線コンピュータ断層撮影装置
US6668073B1 (en) * 1998-11-12 2003-12-23 The University Of British Columbia Anthropomorphic film phantom for three-dimensional dosimetry
US6829323B2 (en) * 2002-07-29 2004-12-07 Ge Medical Systems Global Technology Company, Llc Method and system for low dose image simulation for imaging systems
JP4439202B2 (ja) * 2003-05-09 2010-03-24 株式会社東芝 Ｘ線コンピュータ断層撮影装置及び画像ノイズシミュレーション装置
JP4509903B2 (ja) * 2005-09-27 2010-07-21 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007074772A1 (ja) * 2005-12-28 2007-07-05 Hitachi Medical Corporation Ｘ線ｃｔ装置
WO2007114470A1 (ja) * 2006-04-04 2007-10-11 Hitachi Medical Corporation Ｘ線ｃｔスキャンシミュレータ装置、ｘ線ｃｔ装置、及びｘ線ｃｔスキャンシミュレータプログラム
JP2008048880A (ja) * 2006-08-24 2008-03-06 Ge Medical Systems Global Technology Co Llc 画像作成方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017029461A (ja) * 2015-08-03 2017-02-09 株式会社日立製作所 医用画像処理装置および医用画像処理方法
JP2020519863A (ja) * 2017-05-01 2020-07-02 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 定量的分子撮像のための正確なハイブリッドデータセットの生成
JP7326160B2 (ja) 2017-05-01 2023-08-15 コーニンクレッカ フィリップス エヌ ヴェ 定量的分子撮像のための正確なハイブリッドデータセットの生成
KR20200104997A (ko) * 2019-02-28 2020-09-07 (주)토탈소프트뱅크 심혈관 형상변형 시뮬레이션 시스템
KR102228314B1 (ko) * 2019-02-28 2021-03-16 (주) 토탈소프트뱅크 심혈관 형상변형 시뮬레이션 시스템

Also Published As

Publication number Publication date
JP5588859B2 (ja) 2014-09-10
US8270560B2 (en) 2012-09-18
US20100322375A1 (en) 2010-12-23
JPWO2009107658A1 (ja) 2011-07-07

Similar Documents

Publication Publication Date Title
JP5588859B2 (ja) Ｘ線ｃｔスキャンシミュレータ及びｘ線ｃｔ装置
US7142633B2 (en) Enhanced X-ray imaging system and method
JP5047945B2 (ja) Ｘ線ｃｔスキャンシミュレータ装置、ｘ線ｃｔ装置、及びｘ線ｃｔスキャンシミュレータプログラム
US9036771B2 (en) System and method for denoising medical images adaptive to local noise
JP6066596B2 (ja) Ｘ線撮像における散乱補正の方法及びシステム
JP5209979B2 (ja) 無較正の幾何学的構成における三次元撮像の方法及びシステム
WO2015105132A1 (ja) Ｘ線ｃｔ装置及びｃｔ画像表示方法
JP5918374B2 (ja) Ｘ線ｃｔ装置およびｘ線ｃｔ画像の処理方法
US20100208870A1 (en) Scattered radiation correction method and scattered radiation correction apparatus
JP2010269048A (ja) Ｘ線ｃｔ装置
JP2004057831A (ja) 画像形成システムの低線量画像シミュレーションのための方法及びシステム
US20080118128A1 (en) Methods and systems for enhanced accuracy image noise addition
Neubauer et al. Comparing the image quality of a mobile flat-panel computed tomography and a multidetector computed tomography: a phantom study
JP6767997B2 (ja) 歯科用画像生成システムの画像データの画像向上のための方法
CN104665854A (zh) 图像处理器、治疗系统和图像处理方法
WO2019116619A1 (ja) 放射線断層画像処理装置および放射線断層撮影装置
US20140049542A1 (en) Layer Display of Volume Data
EP3809376A2 (en) Systems and methods for visualizing anatomical structures
Duan et al. Deep-learning convolutional neural network-based scatter correction for contrast enhanced digital breast tomosynthesis in both cranio-caudal and mediolateral-oblique views
CN107809954B (zh) 计算机断层切片图像相对于要被成像的对象的深度位置的显示
US20230309935A1 (en) Methods and systems for contrast-to-noise evaluation of computed tomography systems
WO2018116791A1 (ja) 医用画像処理装置及びそれを備えたｘ線ｃｔ装置、医用画像処理方法
Shaheen et al. The influence of position within the breast on microcalcification detectability in continuous tube motion digital breast tomosynthesis
Burrell Computed Tomography Physics and Instrumentation
You New techniques for chest digital tomosynthesis imaging

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09715754

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2010500716

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 12866484

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 09715754

Country of ref document: EP

Kind code of ref document: A1