JP6770121B2 - 画像表示方法 - Google Patents

画像表示方法 Download PDF

Info

Publication number
JP6770121B2
JP6770121B2 JP2019038445A JP2019038445A JP6770121B2 JP 6770121 B2 JP6770121 B2 JP 6770121B2 JP 2019038445 A JP2019038445 A JP 2019038445A JP 2019038445 A JP2019038445 A JP 2019038445A JP 6770121 B2 JP6770121 B2 JP 6770121B2
Authority
JP
Japan
Prior art keywords
data
image
dimension
dimensional
column
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2019038445A
Other languages
English (en)
Other versions
JP2019144250A (ja
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 同方威視技術股▲分▼有限公司
Publication of JP2019144250A publication Critical patent/JP2019144250A/ja
Application granted granted Critical
Publication of JP6770121B2 publication Critical patent/JP6770121B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • 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/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • 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/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/003—Navigation within 3D models or images
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5223—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G06T2207/10124—Digitally reconstructed radiograph [DRR]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30196—Human being; Person
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/436—Limited angle
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/016—Exploded view
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/028—Multiple view windows (top-side-front-sagittal-orthogonal)
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20—Indexing scheme for editing of 3D models
    • G06T2219/2004—Aligning objects, relative positioning of parts

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Public Health (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Optics & Photonics (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Pulmonology (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Computer Graphics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Geometry (AREA)
  • Remote Sensing (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Description

本願は、コンピュータ断層イメージング（ＣＴ）技術に関し、具体にはＣＴシステムにおける画像表示方法に関する。
米国の９・１１事件が発生した後で、航空分野のセキュリティーチェック作業はますます重視されている。ＣＴ技術に基づくセキュリティーチェック製品はたぐいまれな先進性を有しており、性能指標の優勢は明らかである。米国ＴＳＡは、ＣＴ機器及びその検出方法が空港のセキュリティーチェックを実現する重要な方法であると考える。普通の物品機器であるＸ線放射器は透視してイメージングし、明瞭なＤＲ画像を得ることができる。このようなＤＲ（Ｄｉｇｉｔａｌ Ｒａｄｉｏｇｒａｐｈｙ）画像は、ユーザに固定的な印象が形成されている。ユーザはＤＲ画像をよく知っている。ＤＲシステムにおいて、放射線源および探測器は固定で、回転しないものである。これらのＤＲ画像は、通常、頂視角画像または底視角画像であり、直観的で、便利に見られることができる。画像の読取及び画像の判定から形成される豊富な経験は、いずれも直観的なＤＲ画像に基づくものである。ある優秀な画像判定員は、必ず数万点以上の実際なＤＲ画像を見たことがある。通常に、９０日以上の業務を行うことによって優秀な画像判定員を育成することができる。現在、多数の優秀な画像判定員は、医療、セキュリティーチェックなどの分野で作業している。ＤＲ画像が存在しない場合、これらの経験豊富な画像判定員は、正確な画像判定を行うことができない。
現在のＣＴシステムには既製の３次元データがあるが、ＤＲ画像が存在しない。画像判定員が突然に３次元画像を見たとき、慣れないと感じている。例えば、観察角度の不適宜により、３次元画像における危険物（例えば、ピストル）は小さい集塊であり、形態が不明瞭である。画像判定員が突然に３次元画像を見たとき、如何にしてピストルが小さい集塊となるかを分かることができない。画像判定員が３次元画像を単独に観察しても、旅客が何を持っているかを明瞭にすることができない。このため、３次元画像が存在しても、画像判定員も３次元画像を以前によく知っているＤＲ画像と比較し、この２つの画像を対比し、既存のＤＲ画像の豊富な画像判定経験を利用してほしい。すなわち、画像判定員は、左辺にＤＲ画像を表示するディスプレイが置かれ、右辺にＣＴ画像を表示するディスプレイが置かれるような環境を望んでいる。既存のＤＲ画像の判定経験によって、多数の画像判定員は、旅客が何を持っているかをより迅速、簡便、正確に把握することができる。このように、多くの手間を省くことができる。
また、ＤＲシステムを有するＣＴ機器において、ユーザが特定の位置を選択する場合には、特定の位置のスライスを表示させるためにＤＲ−ＣＴの位置合わせが必要となる。ハードウェア装置のベルト番号に基づいて、特定の位置が荷物のヘッド、中部、又はテールであるかが分かる。しかしながら、搬送過程においてベルト番号が紛失又はなくなりやすく、あるいは、マッチングできない場合があるため、この方法によって正確に位置合わせることができない。新たな効果的な位置合わせ方法が期待されている。
従来技術の１つまたは複数の技術課題を考慮し、ＣＴシステムにＤＲイメージング機器を備えない場合に、ＤＲ画像と被検体の３次元画像とを同時に表示できるＣＴシステムの画像表示方法を提供する。
本発明の一局面において、ＣＴシステムの画像表示方法を提供し、被検体に対してＣＴ走査を行い、ＣＴ投影データを取得するステップと、ＣＴ投影データを所定の間隔で整理するステップと、ある固定角度を最初角度として３６０度を間隔とし、整理されたＣＴ投影データから基本データを抽出するステップと、抽出された基本データに基づいてＤＲ画像を形成するステップと、ＣＴ投影データから被検体の３次元画像を再構成するステップと、ＤＲ画像と再構成された３次元画像とを同時にスクリーンに表示するステップと、を含む。
いくつかの実施例によれば、前記画像表示方法は、前記固定角度に１８０度を加えて最初角度とし、３６０度を間隔として、整理されたＣＴ投影データから補充データを抽出するステップと、抽出された補充データを行列に記憶するステップと、前記行列に対して左右鏡像処理を行うステップと、行を単位として、鏡像処理された行列中の補充データを前記基本データに間隔挿入し、補充された基本データを取得するステップと、補充された基本データに基づいてＤＲ画像を形成するステップと、を含む。
いくつかの実施例によれば、前記画像表示方法は、別の固定角度を最初角度として３６０度を間隔とし、整理されたＣＴ投影データから基本データを抽出するステップと、抽出された基本データに基づいて別のＤＲ画像を形成するステップと、前記ＤＲ画像、前記別のＤＲ画像と再構成された３次元画像とを同時にスクリーンに表示するステップと、を含む。
いくつかの実施例によれば、ＣＴ投影データを所定の間隔で整理するステップは、サンプリング周波数が所定の間隔あたりにデータを１回投影することより高い場合に、サブサンプリング規則に従って、ＣＴ投影データを所定の間隔あたりにデータを１回投影することに調整するステップと、サンプリング周波数が所定の間隔あたりにデータを１回投影することより低い場合に、線形補間に従って、ＣＴ投影データを所定の間隔あたりにデータを１回投影することに調整するステップと、を含む。
いくつかの実施例によれば、前記基本データが高いエネルギーデータ及び低いエネルギーデータである場合に、物質識別アルゴリズムによって前記基本データを処理して階調データ及び材料データを取得する。
いくつかの実施例によれば、ＤＲ画像をユーザが要求する解像度に補間して、スクリーンに表示する。
いくつかの実施例によれば、再構成された３次元画像のデータは、物理密度データ、電子密度データ、原子番号データ、高いエネルギーデータ、低いエネルギーデータ、階調データ、材料データの１つである。
本発明の他の局面において、ＣＴシステムの画像表示方法を提供し、被検体に対してＣＴ走査を行い、ＣＴ投影データを取得するステップと、ＣＴ投影データから被検体の３次元画像を再構成するステップと、被検体の３次元画像データＨ（ｘ，ｙ，ｚ）をある視角で投影し、この視角でのＤＲデータを取得するステップと、この視角でのＤＲデータに基づいてＤＲ画像を形成するステップと、ＤＲ画像と３次元画像とを同時にスクリーンに表示するステップと、を含み、その中、３次元画像データＨ（ｘ，ｙ，ｚ）の各次元のサイズはＸ×Ｙ×Ｚであり、次元ｘは１からＸまで変化し、この変化は水平面内におけるベルト運動に垂直する方向であり、次元ｙは１からＹまで変化し、この変化は垂直方向の上向き方向であり、次元ｚは１からＺまで変化し、この変化は水平面内におけるベルト運動に沿う方向である。
いくつかの実施例によれば、前記画像表示方法は、被検体の３次元画像データＨ（ｘ，ｙ，ｚ）を別の視角で投影し、この視角でのＤＲデータを取得するステップと、この別の視角でのＤＲデータに基づいてＤＲ画像を形成するステップと、２つのＤＲ画像と３次元画像とを同時にスクリーンに表示するステップと、を含む。
いくつかの実施例によれば、前記視角及び前記別の視角は、水平面内におけるベルト運動に垂直する方向、垂直方向及び水平面内におけるベルト運動に沿う方向から選択されるものである。
いくつかの実施例によれば、３次元データＨ（ｘ，ｙ，ｚ）を次元ｘの方向に沿って投影し、２次元データＩ(ｙ,ｚ)を取得し、投影式は以下の通りである：
当該２次元データは側視角のＤＲデータであり、データの次元サイズはＹ×Ｚである。
いくつかの実施例によれば、３次元データＨ（ｘ，ｙ，ｚ）を次元ｙの方向に沿って投影し、２次元データＪ(ｘ,ｚ)を取得し、投影式は以下の通りである：
当該２次元データは底視角または頂視角のＤＲデータであり、データの次元サイズはＸ×Ｚである。
いくつかの実施例によれば、３次元データＨ（ｘ，ｙ，ｚ）を次元ｚであるベルト運動の方向に沿って投影し、２次元データＫ(ｘ,ｙ)を取得し、投影式は以下の通りである：
当該２次元データは前視角のＤＲデータであり、データの次元サイズはＸ×Ｙである。
いくつかの実施例によれば、再構成された３次元画像のデータは、物理密度データ、電子密度データ、原子番号データ、高いエネルギーデータ、低いエネルギーデータ、階調データ、材料データの１つである。
上記技術案によって、ＣＴ投影データを処理してＤＲデータを取得し、ＤＲデータを取得した後で、ＤＲデータ処理アルゴリズムによって、ＤＲ画像を取得する。これにより、画像判定員は既存のＤＲ画像の画像判定経験によって、旅客が持っている物品をより正確、快速に検査することができる。
本発明のまた他の局面において、ＤＲ−ＣＴシステムにおいて画像の位置合わせを行う方法を提供し、前記ＤＲ−ＣＴイメージングシステムはＤＲ機器およびＣＴ機器を含み、前記方法は、ＣＴ機器によって被検体を走査してＣＴ投影データを取得するステップと、前記ＣＴ投影データから被検体の３次元画像を再構成するステップと、ＤＲ機器によって被検体を走査してＤＲ画像を取得するステップと、前記ＤＲ画像及び前記３次元画像をスクリーンに表示するステップと、ユーザのＤＲ画像中における特定の位置に対する選択に応答して、ＤＲ画像のデータから対応する列を抽出するステップと、前記３次元画像から前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップと、前記ＤＲ画像の前記列と補助ＤＲデータの各列との関連性を算出するステップと、補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示するステップと、を含む。
いくつかの実施例によれば、前記３次元画像から前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは、次元ｙ方向に沿って被検体の３次元画像データＨ（ｘ，ｙ，ｚ）を投影し、この視角でのＤＲデータを取得するステップを含み、その中、３次元画像データＨ（ｘ，ｙ，ｚ）の各次元のサイズはＸ×Ｙ×Ｚであり、次元ｘは１からＸまで変化し、この変化は水平面内におけるベルト運動に垂直する方向であり、次元ｙは１からＹまで変化し、この変化は垂直方向の上向き方向であり、次元ｚは１からＺまで変化し、この変化は水平面内におけるベルト運動に沿う方向である。
いくつかの実施例によれば、３次元データＨ（ｘ，ｙ，ｚ）を次元ｙの方向に沿って投影し、２次元データＪ(ｘ,ｚ)を取得し、投影式は以下の通りである：
この２次元データは底視角または頂視角のＤＲデータであり、データの次元サイズはＸ×Ｚである。
本発明のまた一局面において、ＤＲ−ＣＴシステムにおいて画像の位置合わせを行う方法を提供し、前記ＤＲ−ＣＴイメージングシステムはＤＲ機器およびＣＴ機器を含み、前記方法は、ＣＴ機器によって被検体を走査してＣＴ投影データを取得するステップと、ＤＲ機器によって被検体を走査してＤＲ画像を取得するステップと、前記ＤＲ画像をスクリーンに表示するステップと、ユーザのＤＲ画像中における特定の位置に対する選択に応答して、ＤＲ画像のデータから対応する列を抽出するステップと、前記ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップと、前記ＤＲ画像の前記列と補助ＤＲデータの各列との関連性を算出するステップと、前記ＣＴ投影データから被検体の３次元画像を再構成するステップと、補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示するステップと、を含む。
いくつかの実施例によれば、前記ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは、ＣＴ投影データを所定の間隔で整理するステップと、ある固定角度を最初角度として３６０度を間隔とし、整理されたＣＴ投影データから前記補助データを抽出して取得するステップと、を含む。
いくつかの実施例によれば、前記方法は、前記固定角度に１８０度を加えて最初角度とし、３６０度を間隔として、整理されたＣＴ投影データから補充データを抽出するステップと、抽出された補充データを行列に記憶するステップと、前記行列に対して左右鏡像処理を行うステップと、行を単位とし、鏡像処理された行列中の補充データを前記補助ＤＲデータに間隔挿入し、補充された補助ＤＲデータを取得するステップと、を含む。
以下は、本発明をより良く理解できるように、図面に基づいて本発明を詳細に説明する。
本発明の実施例によるＣＴ機器の構造模式図を示す。 図１に示すようなコンピュータデータプロセッサの構造ブロック図を示す。 本発明の実施形態によるコントローラの構造ブロック図を示す。 本発明の一実施例によるＣＴシステムにおいて画像を表示する方法を説明するフローチャートである。 ＣＴ投影データを整理する過程を説明する模式図である。 ＣＴ投影データから得られたＤＲ画像の模式図を示す。 本発明の別の実施例によるＣＴシステムにおいて画像を表示する方法を説明するフローチャートである。 本発明のまた一実施例によるＣＴシステムにおいてＤＲ画像とＣＴ画像を位置合わせる方法を説明するフローチャートである。 本発明のまた別の一実施例によるＣＴシステムにおいてＤＲ画像とＣＴ画像とを位置合わせる方法を説明するフローチャートである。
以下、本発明の具体的な実施例を詳細に説明する。ここで説明した実施例は、例として説明するためのものであり、本発明は、これに限らないと理解すべきである。以下の説明において、本発明に対する透徹した理解をさせるため、大量の特定の細部を描写した。しかし、必ずこれらの特定の細部を採用して本発明を実現することではないことが当業者にとって明らかになっている。その他の実例においては、本発明との混同を避けるために、周知の回路、材料または方法に対する具体的な説明を省略した。
本明細書の全体において、言及した「一実施例」、「実施例」、「一示例」または「示例」は、該実施例または示例に結合して描写した特定の特徴、構造または特性が、本発明の少なくとも一実施例に含まれていることを意味する。従って、明細書全体の各箇所に現れた短文「一実施例において」、「実施例において」、「一示例」または「示例」は、必ず同一の実施例または示例を指したものではない。また、何らかの適宜な組み合わせ及び/またはサブ組み合わせによって、特定の特徴、構造または特性を一つまたは複数の実施例または示例に組み合わせることができる。また、当業者は、ここの「及び／又は」という用語が一つ又は複数の関連するアイテムの任意及び全ての組合を含むことを理解すべきである。
従来技術のＣＴ機器にＤＲイメージングモジュールを備えないが画像判定の過程において画像判定員がＤＲ画像を表示する必要があるという課題に対して、本発明のいくつかの実施例は、ＣＴ投影データから被検体のＤＲ画像を取得し、それと再構成された３次元画像とを同時にスクリーンに表示し、これにより、画像判定員が既存の画像判定の経験によってセキュリティーチェックを行うことを提供する。上記技術案によって、旅客が持っている物品をより正確、快速に検査することができる。いくつかの実施例において、再構成された３次元画像のデータから特定の視角での投影データを抽出し、これらの投影データを組み合わせてこの視角でのＤＲデータとなる。他の実施例において、再構成された３次元画像のデータを異なる視角で投影し、この視角でのＤＲデータを取得し、ＤＲ画像を形成する。
他の実施例において、ＤＲ機器が配置されたＣＴ機器においてベルト番号がなくなり又は紛失したことによるＤＲ画像とＣＴ画像とを位置合わせることができない課題に対して、ＣＴ投影データから、または再構成された３次元画像のデータから、特定の視角でのＤＲ画像を取得し、ＤＲ機器により取得されたＤＲ画像の選択列と、３次元画像のデータから得られたＤＲ画像とに対して関連性分析を行い、関連性が最大である列が位置するスライスを確定して、ＤＲ画像と同時に表示することを提供する。
図１は本発明の実施形態によるＣＴ機器の構造模式図である。図１に示すように、本実施形態によるＣＴ機器は、フレーム２０、搭載機構４０、コントローラ５０、コンピュータデータプロセッサ６０などを含む。フレーム２０は、例えばＸ線放射器という検査用Ｘ放射線を放射する放射線源１０と、探測・採集装置３０とを含む。搭載機構４０は、被検査荷物７０をフレーム２０の放射線源１０と探測・採集装置３０との間の走査領域を通すように搭載するとともに、フレーム２０が被検査荷物７０の直進方向を取り込むように回し、これにより、放射線源１０からの放射線が被検査荷物７０を透過し、被検査荷物７０に対してＣＴ走査を行うことができる。
探測・採集装置３０は、例えば全体モジュール構造を有する探測器・データ採集器であり、例えば平板探測器であり、検出される液体物品を透過する放射線を探測し、アナログ信号を取得し、アナログ信号をデジタル信号に変換し、これにより、被検査荷物７０のＸ放射線に対する投影データを出力する。コントローラ５０は、システム全体の各部分を同期して動作するように制御する。コンピュータデータプロセッサ６０は、データ採集器により採集されたデータを処理し、データを処理して再構成し、結果を出力する。
図１に示すように、放射線源１０は、被検体が放置可能な側に置かれ、探測・採集装置３０は被検査荷物７０の他側に置かれ、探測器・データ採集器を含み、被検査荷物７０の多角度の投影データを取得するためである。データ採集器は、データ拡大成形回路を含み、（電流）積分方式又はパルス（カウント）方式で動作することができる。探測・採集装置３０のデータ出力ケーブルは、コントローラ５０及びコンピュータデータプロセッサ６０に接続され、トリガーコマンドに基づいて、採集されたデータをコンピュータデータプロセッサ６０に記憶する。
図２は図１に示すようなコンピュータデータプロセッサ６０の構造ブロック図を示す。図２に示すように、データ採集器が採集したデータは、インタフェース手段６８及びバス６４によってメモリ６１に記憶される。リードオンリーメモリ（ＲＯＭ）６２にコンピュータデータプロセッサの配置情報及びプログラムが記憶される。ランダムアクセスメモリ（ＲＡＭ）６３は、プロセッサ６６の動作過程に各種のデータを一時に記憶する。なお、メモリ６１に、データ処理を行うためのコンピュータプログラムが更に記憶される。内部バス６４は、上記したメモリ６１、リードオンリーメモリ６２、ランダムアクセスメモリ６３、入力装置６５、プロセッサ６６、表示装置６７及びインタフェース手段６８を接続している。
ユーザが例えばキーボード及びマウスなどの入力装置６５によって操作コマンドを入力した後で、コンピュータプログラムの指令コードコマンドプロセッサ６６は、所定のデータ処理アルゴリズムを実行し、データ処理の結果を取得した後で、それを例えばＬＣＤディスプレイなどの表示装置６７に表示し、あるいは、例えば印刷などのハードコピーの形式によって処理結果を出力する。
図３は本発明の実施形態によるコントローラの構造ブロック図を示す。図３に示すように、コントローラ５０は、コンピュータ６０からの指令によって、放射線源１０、搭載機構４０、探測・採集装置３０を制御する制御手段５１と、制御手段の制御によって、放射線源１０、探測・採集装置３０、搭載機構４０の動作をトリガーするためのトリガーコマンドを生成するトリガー信号生成手段５２と、トリガー信号生成手段５２が制御手段５１の制御によって生成したトリガーコマンドに基づいて、搭載機構４０を被検査荷物７０を搬送するように駆動する第１の駆動機器５３と、トリガー信号生成手段５２が制御手段５１の制御によって生成したトリガーコマンドに基づいて、フレーム２０を回転するように駆動する第２の駆動機器５４と、を含む。探測・採集装置３０により取得された投影データがコンピュータ６０に記憶され、ＣＴ断層画像の再構成を行って、被検査荷物７０の断層画像データを取得する。そして、コンピュータ６０は、例えばソフトウェアを実行することにより、断層画像データから被検査荷物７０の少なくとも１つの視角でのＤＲ画像を取得し、再構成された３次元画像と同時に表示し、これにより、画像判定員がセキュリティーチェックを便利に行うことができる。他の実施例によって、上記ＣＴイメージングシステムは、ダブルエネルギーＣＴシステムであってもよく、すなわち、フレーム２０の放射線源１０は、高いエネルギー及び低いエネルギーという２つの放射線を放射することができ、探測・採集装置３０が異なるエネルギーレベルの投影データを探測した後で、コンピュータデータプロセッサ６０によってダブルエネルギーＣＴの再構成を行って、被検査荷物７０の各断層の同価原子番号及び同価電子密度データを取得する。
図４は本発明の一実施例によるＣＴシステムにおいて画像を表示する方法を説明するフローチャートである。図４に示すように、ステップＳ４１において、ＣＴシステムによって被検体７０に対してＣＴ走査を行い、ＣＴ投影データを取得する。ここのＣＴ投影データは、物理密度又は電子密度データ、原子番号データ、高いエネルギーデータ、低いエネルギーデータ、階調データ、材料データであってもよいが、これらに限定されるものではない。この実施例において、探測器の列数は１列を例とする。探測器が複数列である場合、複数列を全体１列としてデータ処理を行う。
ステップＳ４２において、ＣＴ投影データを所定の間隔で整理する。例えば、ＣＴ投影データが１度を間隔として整理される。すなわち、１度あたりに、１回のデータ投影がある。いくつかの実施例において、サンプリング周波数が１度あたりにデータを１回投影することより高い場合に、サブサンプリング規則に従って、１度あたりにデータを１回投影することに調整する。サンプリング周波数が低すぎる場合に、線形補間に従って、１度あたりにデータを１回投影することに調整する。
例えば、ＣＴ投影データであるサイノグラムを行列Ａに整理して並ぶことができる。そのエレメントをＡ(ｍ,ｎ)にする。ｎは探測器のアレイ要素を示し、ｎは探測器アレイの１〜Ｎに変化し、Ｎは探測器アレイにおけるアレイ要素の個数である。ｍは角度を示し、ｍの変化は角度数の変化であり、すなわち、行列Ａの第１行目は第１度の探測器アレイの１〜Ｎデータであり、第２行目は第２度の探測器アレイの１〜Ｎデータである。全ての角度は３６０を周期とし、図５のように周期的に配列される：１度〜３６０度、３６１度〜７２０度、７２１度〜１０８０度、１０８１度〜１４４０度…。
ステップＳ４３において、固定角度を最初角度とし、３６０度を間隔として、基本データを整理されたＣＴ投影データから抽出する。例えば、ＣＴ投影データから、円周期３６０度を間隔として、固定の最初角度Ｓで、例えば、１度、３６１度、７２１度、１０８１度、１４４１度という若干の角度のデータを抽出する。仮に、固定の最初角度がＳであれば、第ｉ回目の抽出で対応する角度数Ａｎｇｌｅは以下の通りである：
Ａｎｇｌｅ＝Ｓ＋ｉ×３６０
ただし、ｉ＝０,１,２, … （１）
抽出された基本データを行列Ｂに改めて組み合わせる。
ステップＳ４４において、図６に示すように、抽出された基本データに基づいてＤＲ画像を形成する。例えば、組み合わせ直された行列Ｂを既存のＤＲデータ処理アルゴリズムに従って処理してＤＲ画像を形成する。データが高いエネルギーデータ及び低いエネルギーデータである場合に、物質識別アルゴリズムによって階調データ及び材料データを取得する。
ステップＳ４５において、コンピュータデータプロセッサ６０は、ソフトウェアを運行してＣＴ投影データから被検体の３次元画像を再構成する。
ステップＳ４６において、ＤＲ画像と再構成された３次元画像とを同時にスクリーンに表示する。いくつかの状況において、ＤＲ画像をユーザが要求する解像度に補間して、ディスプレイに表示する。
他の実施例において、行列Ｂのデータ量が不十分であり、データを増加する必要がある場合、以下の式に基づいて抽出を続ける。仮に、固定の最初角度がＳであれば、第ｉ回目の抽出で対応する角度数Ａｎｇｌｅは、以下の通りである：
Ａｎｇｌｅ＝ Ｓ＋ｉ×３６０＋１８０
ただし、ｉ＝０,１,２, … （２）
抽出されたデータを行列Ｂ１に改めて組み合わせる。
そして、行列Ｂ１に対して左右鏡像処理を行い、すなわち、左右に対称するエレメントを互いに交換する。式は以下の通りである：
Ｂ１(ｍ,１)<＝＝>Ｂ１(ｍ,Ｎ)
Ｂ１(ｍ,２)<＝＝>Ｂ１(ｍ,Ｎ−１) （３）
……
次に、行を単位として、得られた行列Ｂ１を行列Ｂに間隔挿入し、行列Ｂと組み合わせて、行列Ｃを形成する。データの次元サイズはＭ×Ｎである。組み合わせの式は以下の通りである：
（４）
最後に、行列Ｃを既存のＤＲデータ処理アルゴリズムに従って処理してＤＲ画像を形成する。データが高いエネルギーデータ及び低いエネルギーデータである場合に、物質識別アルゴリズムによって階調データ及び材料データを取得する。
上記実施例において、固定の最初角度ＳはＤＲ画像の視角方向を決定した。例えば、Ｓ＝１である場合に、頂視角を得ることができる；Ｓ＝１８０である場合に、底視角を得ることができる；Ｓ＝９０である場合に、左側視角を得ることができる；Ｓ＝２７０である場合に、右側視角を得ることができる。
以上の実施例において、１つの視角でのＤＲ画像を取得する。必要であれば、上記角度Ｓと異なる別の固定角度を最初角度として当該角度のＤＲデータを形成することができ、対応するＤＲ画像を取得した後で上記ＤＲ画像と再構成された３次元画像とを同時にスクリーンに表示し、画像判定員がセキュリティーチェックを行う。
また、いくつかの実施例において、ＣＴシステムの探測器が一列又は複数列であってもよく、そして、ＣＴデータによっていずれか１つの視角（頂視角、底視角、側視角を含む）の画像を取得してもよい。あるいは、ＣＴデータによっていずれか２つの、すなわち、ダブル視角の画像を取得してもよい。そして、ＣＴデータによっていずれか複数の、すなわち、複数視角の画像を取得してもよい。
図７は本発明の別の実施例によるＣＴシステムにおいて画像を表示する方法を説明するフローチャートである。図７に示すように、ステップＳ７１において、ＣＴシステムによって被検体に対してＣＴ走査を行い、ＣＴ投影データを取得し、そして、ステップＳ７２において、コンピュータ６０によってＣＴ投影データから被検体の３次元画像を再構成する。ここの再構成された３次元データは、物理密度又は電子密度データ、原子番号データ、高いエネルギーデータ、低いエネルギーデータ、階調データ、材料データであることができるが、これらに限定されるものではない。
ステップＳ７３において、被検体の３次元画像データＨ（ｘ，ｙ，ｚ）を一視角で投影し、当該視角でのＤＲデータを取得し、その中、３次元画像データＨ（ｘ，ｙ，ｚ）の各次元のサイズはＸ×Ｙ×Ｚであり、次元ｘは１からＸまで変化し、この変化は水平面内におけるベルト運動に垂直する方向であり、次元ｙは１からＹまで変化し、この変化は垂直方向の上向き方向であり、次元ｚは１からＺまで変化し、この変化は水平面内におけるベルト運動に沿う方向である；
例えば、必要な視角で投影すると、具体的には３つの状況に分かられる：
(１)３次元データＨ（ｘ，ｙ，ｚ）を次元ｘの方向に沿って投影し、２次元データＩ(ｙ,ｚ)を取得し、投影式は以下の通りである：
（５）
当該２次元データは側視角のＤＲデータである。データの次元サイズはＹ×Ｚである。
(２)３次元データＨ（ｘ，ｙ，ｚ）を次元ｙの方向に沿って投影し、２次元データＪ(ｘ,ｚ)を取得し、投影式は以下の通りである：
（６）
当該２次元データは底視角または頂視角のＤＲデータである。データの次元サイズはＸ×Ｚである。
(３)３次元データＨ（ｘ，ｙ，ｚ）を次元ｚであるベルト運動の方向に沿って投影し、２次元データＫ(ｘ,ｙ)を取得し、投影式は以下の通りである：
（７）
この２次元データは前視角のＤＲデータである。データの次元サイズはＸ×Ｙである。
ステップＳ７４において、この視角でのＤＲデータに基づいてＤＲ画像を形成する。例えば、ＤＲデータを取得した後で、既存のＤＲデータ処理アルゴリズムによって処理してＤＲ画像を形成する。データが高いエネルギーデータ及び低いエネルギーデータである場合、物質識別アルゴリズムによって、階調データ及び材料データを取得する。
ステップＳ７５において、ＤＲ画像と３次元画像とを同時にスクリーンに表示する。例えば、ＤＲ画像をユーザが要求する解像度に補間して、ディスプレイに表示する。
上記にＣＴシステムにＤＲ機器を備えない状況が記述されるが、ＤＲ機器を備えるいくつかのＣＴシステムにおいて、ＤＲ画像とスライス画像とを位置合わせる必要がある。例えば、ＤＲシステムを備えるＣＴ機器には、ユーザがＤＲ画像のある所をクリックした後で、３次元データからスライスの正確な位置が荷物のヘッド、中部又はテールに位置するかを取得して、当該スライスを表示することが要求される。３次元データに若干のスライスがあり、３次元データから正確な位置を取得してスライスを表示するために、ハードウェア装置のベルト番号に基づいて荷物のヘッド、中部又はテールであるかを分かることができる。しかしながら、ベルト番号が搬送過程において紛失し又はなくなり、あるいは、マッチングできない場合に、当該方法によって正確に位置合わせることができない。
上記課題に対して、本発明の実施例は、ＤＲ機器により得られたＤＲ画像の一列と再構成された３次元画像から抽出されたＤＲデータの各列に対して関連性算出を行い、関連性が最大である列に対応するスライス画像とＤＲ画像とを同時にスクリーンに表示することを提供する。ここのＤＲデータは、物理密度又は電子密度データ、原子番号データ、高いエネルギーデータ、低いエネルギーデータ、階調データ、材料データなどであることができるが、これらに限定されるものではない。
図８は本発明のまた一実施例によるＣＴシステムにおいてＤＲ画像とＣＴ画像を位置合わせる方法を説明するフローチャートである。
ステップＳ８１において、ＣＴ機器によって被検体を走査してＣＴ投影データを取得する。ステップＳ８２において、コンピュータ６０がソフトウェアを実行することによってＣＴ投影データから被検体の３次元画像を再構成する。
ステップＳ８３において、ＤＲ機器によって被検体を走査してＤＲ画像を取得し、そして、ステップＳ８４において、ＤＲ機器により得られたＤＲ画像と再構成された３次元画像をスクリーンに表示し、画像判定員が便利にセキュリティーチェックを行う。
画像判定員は、より正確な検査を行うため、ＤＲ画像のある列に対応するスライス画像を見たい場合には、例えばマウスなどの入力装置６５を操作して、ＤＲ画像のその列をクリックする。ステップＳ８５において、ユーザのＤＲ画像中における特定の位置に対する選択に応答して、ＤＲ画像のデータから対応する列を抽出する。次に、ステップＳ８６において、コンピュータ６０は、３次元画像から、ＤＲ画像と同じ視角での補助ＤＲデータを取得する。
例えば、ＤＲシステムはＤＲデータを生成し、ユーザは、当該ＤＲデータにおいて、ある位置をマウスで指定する。この位置が位置するＤＲデータを探測器アレイ要素の方向に沿って広げて、Ａｒｒａｙ０(ｎ)として記述し、ただし、ｎは１からＮまで変化する。Ｎは探測器アレイのアレイ要素の個数である。当該データは、実際に配列であり、全ての探測器アレイ要素のデータが記憶されている。
図７に示すような方法によって、ＣＴデータからＤＲデータを取得し、行列Ｊとして記述し、そのエレメントはＪ(ｘ,ｚ)である。行列Ｊを入り替えてＪ１を取得し、そのエレメントはＪ１(ｚ,ｘ)であり、Ｊ１を要求された次元サイズＭ×Ｎに従って補間して行列Ｃを取得し、データの次元サイズはＭ×Ｎである。
いくつかの実施例において、３次元画像から前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは以下のことを含む：被検体の３次元画像データＨ（ｘ，ｙ，ｚ）を次元ｙ方向に沿って投影し、この視角でのＤＲデータを取得し、その中、３次元画像データＨ（ｘ，ｙ，ｚ）の各次元のサイズはＸ×Ｙ×Ｚであり、次元ｘは１からＸまで変化し、この変化は水平面内におけるベルト運動に垂直する方向であり、次元ｙは１からＹまで変化し、この変化は垂直方向の上向き方向であり、次元ｚは１からＺまで変化し、この変化は水平面内におけるベルト運動に沿う方向である。例えば、３次元データＨ（ｘ，ｙ，ｚ）を次元ｙの方向に沿って投影し、上記式（６）を参照して２次元データＪ(ｘ,ｚ)を取得する。
そして、ステップＳ８７において、コンピュータ６０は前記ＤＲ画像の前記列と補助ＤＲデータの各列との関連性を算出し、ステップＳ８８において、補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示する。
例えば、行列Ｃのある行を抽出し、例えばｍ＝１である場合、ｎは１からＮまで変化し、配列Ａｒｒａｙ１(ｎ)に記入する。下記の式に基づいて誤差Ｅを算出する。
（８）
行列Ｃの全ての行が走査されるまでに、上記ステップを繰返し、次の行を抽出し、例えばｍ＝２である場合に、ｎは１からＮまで変化し、配列Ａｒｒａｙ１(ｎ)に記入する。最後に、全てのＥから逐一に探して最小値に対応するｍを見つけて、この位置をマッチングする。３次元データＨ（ｘ，ｙ，ｚ）のスライスデータＨ(ｘ,ｙ,ｚ＝ｍ)を抽出して、ＤＲ画像と同時にディスプレイに表示する。
上記には、３次元画像を再構成した後で補助ＤＲデータを取得して位置合わせる実施例が記述されている。類似的には、３次元画像を再構成する前に、補助ＤＲデータを取得し、画像の位置合わせに用いてもよい。図９は、本発明のまた別の実施例によるＣＴシステムにおいてＤＲ画像とＣＴ画像を位置合わせる方法を説明するフローチャートである。
図９に示すようなＤＲ−ＣＴシステムにおいて画像の位置合わせを行う方法では、ＤＲ−ＣＴイメージングシステムは、ＤＲ機器およびＣＴ機器を含む。ステップＳ９１において、ＣＴ機器によって被検体を走査してＣＴ投影データを取得し、ステップＳ９２において、ＤＲ機器によって被検体を走査してＤＲ画像を取得する。そして、ステップＳ９３において、前記ＤＲ画像をスクリーンに表示する。
画像判定員は、より正確な検査を行うため、ＤＲ画像のある列に対応するスライス画像を見たい場合には、例えばマウスなどの入力装置６５を操作して、ＤＲ画像のある列をクリックする。ステップＳ９４において、ユーザのＤＲ画像中における特定の位置に対する選択に応答して、ＤＲ画像のデータから対応する列を抽出する。
ステップＳ９５において、前記ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得する。
ステップＳ９６において、前記ＤＲ画像の前記列と補助ＤＲデータの各列との関連性を算出する。
ステップＳ９７において、前記ＣＴ投影データから被検体の３次元画像を再構成する。
ステップＳ９８において、補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示する。
いくつかの実施例において、ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは以下のことを含む：ＣＴ投影データを所定の間隔で整理する；ある固定角度を最初角度として３６０度を間隔として、整理されたＣＴ投影データから前記補助データを抽出して取得する。
他の実施例において、上記方法は以下のステップを更に含んでいる：この固定角度に１８０度を加えて最初角度とし、３６０度を間隔とし、整理されたＣＴ投影データから補充データを抽出する。抽出された補充データを行列に記憶する；前記行列に対して左右鏡像処理を行う；行を単位として、鏡像処理された行列中の補充データを前記補助ＤＲデータに間隔挿入して、補充された補助ＤＲデータを取得する。
以上、詳細の記載は、模式図、フローチャット及び/または例を使用することによって、画像表示方法及び画像を位置合わせる方法に係る数多くの実施例を説明した。このような模式図、フローチャット及び/または例が、機能及び/または操作が一つまたは複数含まれた場合には、当業者は、このような模式図、フローチャットまたは例における各機能及び/または操作が、各種の構造、ハードウェア、ソフトウェア、ファームウェアまたは実質上のこれらの任意の組み合わせで、個別及び/または共同で実現できると理解すべきである。一つの実施例において、本発明の実施例の前記主題のいつかの部品は、専用集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、デジタル信号プロセッサ（ＤＳＰ）、またはその他の集積フォーマットで実現できる。しかしながら、当業者は、ここで公開した実施例の一部が、全体または部分的に集積回路で同じく実現することができる。例えば、一台または複数台のコンピュータ上で運行する一つまたは複数のコンピュータプログラム（例えば、一台または複数台のコンピュータシステム上で運行する一つまたは複数のプログラム）によって実現させても良いし、一つまたは複数のプロセッサ上で運行する一つまたは複数のプログラム（例えば、一つまたは複数のマイクロプロセッサ上で運行する一つまたは複数のプログラム）によって実現させても良いし、ファームウェアまたは実質上に上記形態の任意組み合わせによって実現させても良いと理解すべきである。また、当業者は、本開示を元に、回路の設計及び/またはソフトウェアの書き込み及び/またはファームウェアのコーディングの能力を備える。また、当業者には理解されるように、本開示の前記主題のメカニズムは、複数の形態のプログラム製品として配分できると共に、実際に配分の信号載置媒体の具体的な類型が何かであろうか、本開示の前記主題の実施例は何れも適用できる。信号載置媒体の例示は、例えば、ソフトディスク、ハートディスクドライブ、コンパクトディスク（ＣＤ）、デジタル汎用ディスク（ＤＶＤ）、デジタルテープ、コンピュータメモリ等の記録可能な記録型媒体、及び例えば、デジタル及び/またはアナログ通信媒体（例えば、光ファイバ、導波管、有線通信リング、無線通信リングなど）の搬送型媒体を含むが、それらに限らない。
以上、本発明の典型的な実施例に基づいて本発明を説明したが、当業者は、使用された用語が、説明するための例であって、本発明を限定する用語ではないと理解すべきである。また、本発明は、精神及び主旨を逸脱しない限りに、種々の形態で具体的に実施できるので、上記の実施例は、前述の詳細に限らず、特許請求の範囲によって限定されるものとして、広く解釈できると理解すべきである。特許請求の範囲または等価の範囲内での全ての変化や改良は、特許請求の範囲内に含まれていることを理解すべきである。

Claims (6)

  1. ＤＲ−ＣＴイメージングシステムにおいて画像の位置合わせを行う方法であって、前記ＤＲ−ＣＴイメージングシステムはＤＲ機器およびＣＴ機器を含み、前記方法は、
    前記ＣＴ機器によって被検体を走査してＣＴ投影データを取得するステップと、
    前記ＣＴ投影データから前記被検体の３次元画像を再構成するステップと、
    前記ＤＲ機器によって前記被検体を走査してＤＲ画像を取得するステップと、
    前記ＤＲ画像及び前記３次元画像をスクリーンに表示するステップと、
    ユーザの前記ＤＲ画像における特定の位置に対する選択に応答して、前記ＤＲ画像のデータから対応する列を抽出するステップと、
    前記３次元画像から前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップと、
    前記ＤＲ画像の前記列と前記補助ＤＲデータの各列との関連性を算出するステップと、
    前記補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示するステップと、を含む、画像の位置合わせを行う方法。
  2. 前記３次元画像から前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは、
    次元ｙの方向に沿って前記被検体の前記３次元画像のデータＨ（ｘ，ｙ，ｚ）を投影し、前記視角でのＤＲデータを取得するステップを含み、その中、前記３次元画像の前記データＨ（ｘ，ｙ，ｚ）の各次元のサイズはＸ×Ｙ×Ｚであり、次元ｘは１からＸまで変化し、前記次元ｘの方向は水平面内の方向であって前記被検体を搬送するためのベルト運動に垂直する方向であり、次元ｙは１からＹまで変化し、前記次元ｙの方向は前記水平面に垂直する方向の上向き方向であり、次元ｚは１からＺまで変化し、前記次元ｚの方向は前記水平面内における前記ベルト運動に沿う方向である、請求項１に記載の方法。
  3. 前記３次元画像の前記データＨ（ｘ，ｙ，ｚ）を前記次元ｙの方向に沿って投影し、２次元データＪ(ｘ,ｚ)を取得し、投影式は以下の通りである：
    この２次元データは前記水平面内に垂直する方向における底面視の視角または上面視の視角でのＤＲデータであり、データの次元サイズはＸ×Ｚである、請求項２に記載の方法。
  4. ＤＲ−ＣＴイメージングシステムにおいて画像の位置合わせを行う方法であって、前記ＤＲ−ＣＴイメージングシステムはＤＲ機器およびＣＴ機器を含み、前記方法は、
    前記ＣＴ機器によって被検体を走査してＣＴ投影データを取得するステップと、
    前記ＤＲ機器によって前記被検体を走査してＤＲ画像を取得するステップと、
    前記ＤＲ画像をスクリーンに表示するステップと、
    ユーザの前記ＤＲ画像における特定の位置に対する選択に応答して、前記ＤＲ画像のデータから対応する列を抽出するステップと、
    前記ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップと、
    前記ＤＲ画像の前記列と前記補助ＤＲデータの各列との関連性を算出するステップと、
    前記ＣＴ投影データから前記被検体の３次元画像を再構成するステップと、
    前記補助ＤＲデータの、前記ＤＲ画像の前記列との関連性が最大である列に対応するスライス画像と前記ＤＲ画像とを同時にスクリーンに表示するステップと、を含む、画像の位置合わせを行う方法。
  5. 前記ＣＴ投影データから前記ＤＲ画像と同じ視角での補助ＤＲデータを取得するステップは、
    前記ＣＴ投影データを所定の間隔で整理するステップと、
    ある固定角度を第１最初角度として３６０度を間隔とし、整理されたＣＴ投影データから前記補助ＤＲデータを抽出して取得するステップと、を含む、請求項４に記載の方法。
  6. 前記固定角度に１８０度を加えて第２最初角度とし、３６０度を間隔として、整理されたＣＴ投影データから補充データを抽出するステップと、
    抽出された補充データを行列に記憶するステップと、
    前記行列に対して左右鏡像処理を行うステップと、
    行を単位とし、鏡像処理された行列中の補充データを前記補助ＤＲデータに間隔挿入し、補充された補助ＤＲデータを取得するステップと、を含む、請求項５に記載の方法。
JP2019038445A 2014-05-14 2019-03-04 画像表示方法 Active JP6770121B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201410202944.7 2014-05-14
CN201410202944.7A CN105094725B (zh) 2014-05-14 2014-05-14 图像显示方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2018035383A Division JP2018128458A (ja) 2014-05-14 2018-02-28 画像表示方法

Publications (2)

Publication Number Publication Date
JP2019144250A JP2019144250A (ja) 2019-08-29
JP6770121B2 true JP6770121B2 (ja) 2020-10-14

Family

ID=53189681

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2015099504A Pending JP2015227872A (ja) 2014-05-14 2015-05-14 画像表示方法
JP2018035383A Pending JP2018128458A (ja) 2014-05-14 2018-02-28 画像表示方法
JP2019038445A Active JP6770121B2 (ja) 2014-05-14 2019-03-04 画像表示方法

Family Applications Before (2)

Application Number Title Priority Date Filing Date
JP2015099504A Pending JP2015227872A (ja) 2014-05-14 2015-05-14 画像表示方法
JP2018035383A Pending JP2018128458A (ja) 2014-05-14 2018-02-28 画像表示方法

Country Status (7)

Country Link
US (1) US9805504B2 (ja)
EP (1) EP3021283B1 (ja)
JP (3) JP2015227872A (ja)
CN (1) CN105094725B (ja)
HK (1) HK1217796A1 (ja)
RU (1) RU2599596C1 (ja)
WO (1) WO2015172726A1 (ja)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2536650A (en) 2015-03-24 2016-09-28 Augmedics Ltd Method and system for combining video-based and optic-based augmented reality in a near eye display
KR101772327B1 (ko) * 2016-04-29 2017-08-29 (의료)길의료재단 부메랑 형상의 커널을 사용하는 3차원 병렬영상복원기법
DE102016209674B4 (de) * 2016-06-02 2023-10-26 Siemens Healthcare Gmbh Ermittlung einer räumlichen Verteilung von Materialeigenschaftswerten auf Basis einer Einzelenergie-CT-Bildaufnahme mit Hilfe eines iterativen Optimierungsverfahrens
CN106383132B (zh) * 2016-10-17 2019-05-14 北京君和信达科技有限公司 辐射检查系统和方法
US12067650B2 (en) * 2018-03-20 2024-08-20 Nec Corporation Imaging apparatus and imaging method
WO2019211741A1 (en) 2018-05-02 2019-11-07 Augmedics Ltd. Registration of a fiducial marker for an augmented reality system
US10740983B2 (en) 2018-06-01 2020-08-11 Ebay Korea Co. Ltd. Colored three-dimensional digital model generation
US11766296B2 (en) 2018-11-26 2023-09-26 Augmedics Ltd. Tracking system for image-guided surgery
CN109949899B (zh) * 2019-02-28 2021-05-28 未艾医疗技术(深圳)有限公司 图像三维测量方法、电子设备、存储介质及程序产品
CN110349151B (zh) * 2019-07-16 2021-12-03 科大讯飞华南人工智能研究院（广州）有限公司 一种目标识别方法及装置
US11980506B2 (en) 2019-07-29 2024-05-14 Augmedics Ltd. Fiducial marker
CN112834968A (zh) * 2019-11-25 2021-05-25 通用电气精准医疗有限责任公司 用于磁共振成像系统的扫描控制系统及方法
US11382712B2 (en) 2019-12-22 2022-07-12 Augmedics Ltd. Mirroring in image guided surgery
CN113643361A (zh) * 2020-05-11 2021-11-12 同方威视技术股份有限公司 目标区域定位方法、装置、设备、介质和程序产品
CN112652054B (zh) * 2020-12-22 2022-12-20 上海索骥信息科技有限公司 一种dr/ct数字物体构建方法和系统
US11896445B2 (en) 2021-07-07 2024-02-13 Augmedics Ltd. Iliac pin and adapter
WO2023245505A1 (en) * 2022-06-22 2023-12-28 Syngular Technology Limited A system and a method for 3d image processing, and a method for rendering a 3d image
WO2024057210A1 (en) 2022-09-13 2024-03-21 Augmedics Ltd. Augmented reality eyewear for image-guided medical intervention

Family Cites Families (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02205986A (ja) * 1989-02-03 1990-08-15 Toshiba Corp 画像表示装置
US5367552A (en) * 1991-10-03 1994-11-22 In Vision Technologies, Inc. Automatic concealed object detection system having a pre-scan stage
JP3313397B2 (ja) * 1992-05-21 2002-08-12 株式会社東芝 Ｘ線ｃｔ装置
EP0825457A3 (en) 1996-08-19 2002-02-13 Analogic Corporation Multiple angle pre-screening tomographic systems and methods
JP2001520376A (ja) * 1997-10-10 2001-10-30 アナロジック コーポレーション 法線曲面を用いたｃｔ目標物の検出
US5907593A (en) * 1997-11-26 1999-05-25 General Electric Company Image reconstruction in a CT fluoroscopy system
JP2000107173A (ja) 1998-10-08 2000-04-18 Fuji Photo Film Co Ltd ３次元用放射線画像形成装置
JP2000271110A (ja) * 1999-03-26 2000-10-03 Hitachi Medical Corp 医用ｘ線装置
DE19954664B4 (de) * 1999-11-13 2006-06-08 Smiths Heimann Gmbh Vorrichtung zur Bestimmung von kristallinen und polykristallinen Materialien eines Gegenstandes
DE19954662B4 (de) * 1999-11-13 2004-06-03 Smiths Heimann Gmbh Vorrichtung und Verfahren zum Detektieren von unzulässigen Reisegepäckgegenständen
JP2001149368A (ja) * 1999-11-29 2001-06-05 Hitachi Medical Corp 医用ｘ線装置
JP3950782B2 (ja) * 2002-05-22 2007-08-01 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー ３次元逆投影方法およびｘ線ｃｔ装置
JP2004012407A (ja) * 2002-06-11 2004-01-15 Hitachi Ltd 透過画像提供システムおよびｘ線ｃｔ・ｄｒ撮影サービスシステム
JP2004177138A (ja) * 2002-11-25 2004-06-24 Hitachi Ltd 危険物探知装置および危険物探知方法
WO2005001457A1 (en) * 2003-06-27 2005-01-06 Tsinghua University Gamma radiation imaging system for non-destructive inspection of the luggage
US7215734B2 (en) * 2004-06-30 2007-05-08 General Electric Company Method and system for three-dimensional reconstruction of images
CN1841391A (zh) 2004-11-23 2006-10-04 通用电气公司 使用用户映射的感兴趣区域使图像同步的方法和设备
CN2774421Y (zh) 2005-04-12 2006-04-26 车延军 多屏幕显示的手术导航装置
US7496222B2 (en) * 2005-06-23 2009-02-24 General Electric Company Method to define the 3D oblique cross-section of anatomy at a specific angle and be able to easily modify multiple angles of display simultaneously
JP4711759B2 (ja) * 2005-07-08 2011-06-29 パナソニック株式会社 Ｘ線検査装置
CN100495439C (zh) * 2005-11-21 2009-06-03 清华大学 采用直线轨迹扫描的图像重建系统和方法
JP2007280354A (ja) * 2006-03-16 2007-10-25 Ricoh Co Ltd 三次元形状処理装置、三次元形状処理方法、三次元形状処理プログラム、記録媒体、パーツカタログ・システム、パーツカタログ作成方法、及びプログラム
CN101071111B (zh) * 2006-05-08 2011-05-11 清华大学 一种多视角航空集装箱安全检查系统及方法
CN101112316A (zh) * 2006-07-28 2008-01-30 Ge医疗系统环球技术有限公司 X射线混合诊断系统
CN200956017Y (zh) * 2006-09-08 2007-10-03 清华大学 集装箱多重dr/ct检测装置
US7894649B2 (en) * 2006-11-02 2011-02-22 Accuray Incorporated Target tracking using direct target registration
CN101271075A (zh) * 2007-03-22 2008-09-24 清华大学 一种ct扫描安全检查方法及其装置
WO2008120136A1 (en) * 2007-03-30 2008-10-09 Koninklijke Philips Electronics N.V. 2d/3d image registration
US8355595B2 (en) * 2007-05-15 2013-01-15 Xerox Corporation Contrast enhancement methods and apparatuses
CN101382507A (zh) * 2007-09-05 2009-03-11 同方威视技术股份有限公司 一种检查航空货运集装箱中违禁物品的装置
CN101403710B (zh) * 2007-10-05 2013-06-19 清华大学 液态物品检查方法和设备
CN102608135B (zh) * 2007-10-05 2015-02-18 清华大学 在危险品检查系统中确定ct扫描位置的方法和设备
JP5171215B2 (ja) * 2007-11-08 2013-03-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置
CN101470082B (zh) * 2007-12-27 2011-03-30 同方威视技术股份有限公司 物品检测装置及其检测方法
US8180139B2 (en) * 2009-03-26 2012-05-15 Morpho Detection, Inc. Method and system for inspection of containers
CN101936924B (zh) 2009-06-30 2012-07-04 同方威视技术股份有限公司 物品检查系统
US9401047B2 (en) * 2010-04-15 2016-07-26 Siemens Medical Solutions, Usa, Inc. Enhanced visualization of medical image data
US8615118B2 (en) * 2010-05-28 2013-12-24 The University Of Maryland, Baltimore Techniques for tomographic image by background subtraction
JP5536607B2 (ja) * 2010-10-13 2014-07-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置
JP2012170736A (ja) * 2011-02-23 2012-09-10 Toshiba Corp Ｘ線コンピュータ断層撮影装置
US20140010437A1 (en) * 2011-03-22 2014-01-09 Ram C. Naidu Compound object separation
US11284846B2 (en) * 2011-05-12 2022-03-29 The John Hopkins University Method for localization and identification of structures in projection images
CN102651145B (zh) * 2012-04-06 2014-11-05 哈尔滨工业大学 股骨三维模型可视化方法
JP2015520846A (ja) * 2012-05-01 2015-07-23 アナロジック コーポレイション Ｃｔ密度画像とスパースマルチエネルギーデータを用いてボクセルセットのｚ−有効値を決定する方法およびシステム｛ｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｚ−ｅｆｆｅｃｔｉｖｅｖａｌｕｅｆｏｒｓｅｔｏｆｖｏｘｅｌｓｕｓｉｎｇｃｔｄｅｎｓｉｔｙｉｍａｇｅａｎｄｓｐａｒｓｅｍｕｌｔｉ−ｅｎｅｒｇｙｄａｔａ｝
JP6153105B2 (ja) * 2013-01-10 2017-06-28 東芝Ｉｔコントロールシステム株式会社 Ｃｔ装置

Also Published As

Publication number Publication date
JP2018128458A (ja) 2018-08-16
CN105094725A (zh) 2015-11-25
JP2019144250A (ja) 2019-08-29
EP3021283B1 (en) 2020-12-09
EP3021283A2 (en) 2016-05-18
HK1217796A1 (zh) 2017-01-20
WO2015172726A1 (zh) 2015-11-19
JP2015227872A (ja) 2015-12-17
CN105094725B (zh) 2019-02-19
US9805504B2 (en) 2017-10-31
US20150332498A1 (en) 2015-11-19
RU2599596C1 (ru) 2016-10-10
EP3021283A3 (en) 2016-09-14

Similar Documents

Publication Publication Date Title
JP6770121B2 (ja) 画像表示方法
CN105784731B (zh) 一种定位三维ct图像中的目标的方法和安检系统
JP5737858B2 (ja) 画像処理装置、画像処理方法、及びプログラム
US8705690B2 (en) Imaging method with improved display of a tissue region, imaging device, and computer program product
JP6968826B2 (ja) 投影画像を生成するためのシステムおよび方法
CN105451657A (zh) 用于导航断层合成堆叠的包括自动聚焦的系统和方法
US10032296B2 (en) Volumertric image data visualization
JP2014508954A (ja) 合成オブジェクト分割方法およびシステム｛ｃｏｍｐｏｕｎｄｏｂｊｅｃｔｓｅｐａｒａｔｉｏｎ｝
US10672517B2 (en) Method and device for monitoring a breast examination
JP2019522775A5 (ja)
US9361711B2 (en) Lesion-type specific reconstruction and display of digital breast tomosynthesis volumes
CN103608843A (zh) 图像显示装置
CN102132322B (zh) 用于确定对象的尺寸改变的设备
WO2016095776A1 (zh) 一种定位三维ct图像中的目标的方法和安检ct系统
JP6073562B2 (ja) 医用画像処理装置
WO2016095798A1 (zh) 一种定位三维ct图像中的目标的方法和安检系统
CN110313928A (zh) 用于确定放射性结果数据的分析方法和分析单元
US9460534B2 (en) Labeling a rib cage by placing a label based on a center point and position of a rib
US20110074781A1 (en) Intermediate image generation method, apparatus, and program
WO2005104953A1 (ja) 画像診断支援装置及びその方法
JP2010240141A (ja) 情報処理装置、画像診断システム、および診断用画像の表示制御方法
US10395395B2 (en) Virtual projection images for tomosynthesis artifact reduction
US10929990B2 (en) Registration apparatus, method, and program
JP2016041227A (ja) 画像表示装置、方法、及びプログラム
US9536346B2 (en) Medical image display apparatus, medical image display method, and medical image display program

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190402

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190402

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200129

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200501

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200924

R150 Certificate of patent or registration of utility model

Ref document number: 6770121

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250