JP2019513225A - Ｘ線検出方法及びｘ線検出器 - Google Patents

Ｘ線検出方法及びｘ線検出器 Download PDF

Info

Publication number
JP2019513225A
JP2019513225A JP2018545666A JP2018545666A JP2019513225A JP 2019513225 A JP2019513225 A JP 2019513225A JP 2018545666 A JP2018545666 A JP 2018545666A JP 2018545666 A JP2018545666 A JP 2018545666A JP 2019513225 A JP2019513225 A JP 2019513225A
Authority
JP
Japan
Prior art keywords
energy
ray
energy window
substance
linear attenuation
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
Application number
JP2018545666A
Other languages
English (en)
Other versions
JP6630841B2 (ja
Inventor
シン，ユーシャン
グオ，シャオユエ
ダァン，ヂー
シェン，ラ
チャン，リー
チェン，ジーチアン
リ，リアン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Publication of JP2019513225A publication Critical patent/JP2019513225A/ja
Application granted granted Critical
Publication of JP6630841B2 publication Critical patent/JP6630841B2/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
    • 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
    • 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • 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/48—Diagnostic techniques
    • A61B6/482—Diagnostic techniques involving multiple energy imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/10—Different kinds of radiation or particles
    • G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
    • G01N2223/1016—X-ray
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/419—Imaging computed tomograph
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/423—Imaging multispectral imaging-multiple energy imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/64—Specific applications or type of materials multiple-sample chamber, multiplicity of materials
    • 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]

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measurement Of Radiation (AREA)
  • Spectroscopy & Molecular Physics (AREA)

Abstract

新規なＸ線検出方法及びＸ線検出器が提供される。本発明の実施例に係るＸ線検出方法は、Ｘ線源によって放射された光子のエネルギー範囲をＮ（Ｎは０より大きい整数である）個のエネルギーウィンドウに分割するステップと、イメージング対象の注目物質と背景物質との線減弱係数に基づいて、Ｎ個のエネルギーウインドウのそれぞれの重み係数を取得するステップと、前記Ｎ個のエネルギーウィンドウのそれぞれの重み係数に基づいて、Ｘ線検出器のＭ（Ｍは０より大きい整数である）個の出力チャネルの重み係数行列を取得するステップと、エネルギー範囲が前記Ｎ個のエネルギーウィンドウのそれぞれに入る光子の数と前記重み係数行列に基づいて、前記Ｍ個の出力チャネルの出力結果を取得するステップとを備える。【選択図】図1

Description

本発明は、Ｘ線イメージング分野に関し、特に、Ｘ線検出方法及びＸ線検出器に関する。
レントゲンがＸ線を発見して以来、Ｘ線によるＣＴイメージング技術は急速に発展し、医用イメージング及び安全検査などの分野に広く使用されている。現在、Ｘ線イメージングに適用される検出器は、主に、エネルギー積分と光子計数という両種の信号取得方法を使用している。近年、材料組成の区別におけるマルチエネルギースペクトルイメージングの利点により、光子計数検出器が関連分野の研究ホットスポットとなっている。世界には、複数の研究機関と企業がＸ線イメージングにおける光子計数検出器の研究と製造を鋭意に行っている。光子計数検出器は、対象物を透過する光子のエネルギーの差異に応じて光子に対してエネルギーごとに計数する。このようにして、閾値を設定することにより、低エネルギー光子をフィルタリングして、ノイズ及び放射線量を低減することができる。同時に、複数の閾値を設定することにより、より広いエネルギースペクトル分布を有するＸ線をエネルギー領域ごとにカウントすることができ、異なるエネルギー領域のイメージング結果を直接に得ることができる。光子計数検出器は、エネルギー積分検出器よりも、イメージングへのノイズの影響を排除して、画質を向上することができる一方、エネルギー領域を細く分けて、Ｘ線イメージングための取得データの関連性を減少させることができる。従って、光子計数検出器を用いたＸ線イメージングは、物質の識別精度、成分の定量精度などを改善することができる。
Ｘ線イメージングにおいて、光子計数検出器がますます重要な役割を果たす。しかし、そのコストは高い一方、計数率及びエネルギーウィンドウチャネルの数に大きく制限されるため、データ収集の効率を最適化することは困難である。理想的な状況を分析すると、異なるエネルギーのＸ線光子が持つ情報は異なり、異なる光子を十分なエネルギー区間によって細分化し、且つエネルギーによって異なる処理を行うことができれば、画質を最大限に最適化することができる。Ｘ線イメージング検出器に適用される上記両種の信号取得方法は、それぞれ、二つの種類の「エネルギー重み付け」信号処理方法に対応する。すなわち、エネルギー積分検出器はエネルギーを重みとして使用してすべての光子信号を加重合計することで出力を得て、光子計数検出器は、設定されたエネルギー区間の範囲において定数を重みとして光子を一つずつ累算して出力を得る。この二つの種類の検出器のそれぞれの出力について、信号の最適化の方法を使用することはできない。これらによるＣＴイメージングシステムの画質は限られている。
本発明は、上記課題の少なくとも一つの課題を解決するために、新規なＸ線検出方法及びＸ線検出装置を提供する。
本発明の実施例のＸ線検出方法は、Ｘ線源によって放射された光子のエネルギー範囲をＮ（Ｎは０より大きい整数である）個のエネルギーウィンドウに分割するステップと、イメージング対象の注目物質と背景物質との線減弱係数に基づいて、Ｎ個のエネルギーウインドウのそれぞれの重み係数を取得するステップと、Ｎ個のエネルギーウィンドウのそれぞれの重み係数に基づいて、Ｘ線検出器のＭ（Ｍは０より大きい整数である）個の出力チャネルの重み係数行列を取得するステップと、エネルギー範囲が前記Ｎ個のエネルギーウィンドウのそれぞれに入る光子の数と前記重み係数行列に基づいて、Ｍ個の出力チャネルの出力結果を取得するステップと、を備える。
本発明の実施例のＸ線検出器は、それぞれ、Ｘ線源によって放射されたエネルギー範囲を選別範囲に入る光子の数をカウントするＮ個のエネルギーウィンドウ選別器と、ｉ組目の乗法器におけるｋ番目の乗算器が、前記Ｎ個のエネルギーウィンドウ選別器におけるｉ番目のエネルギーウィンドウ選別器でカウントして得た光子の数とＭ行×Ｎ列の重み係数行列におけるｋ行目ｉ列目の重み係数とを乗算する、それぞれＭ個の乗算器を備えるＮ（Ｍ、Ｎ、ｋ、ｉのいずれもは０より大きい整数である）組の乗算器と、ｋ番目の加算器が、各組の乗算器におけるｋ番目の乗算器で得た積を累算して、Ｘ線検出器のＭ個の出力チャネルにおけるｋ番目の出力チャネルの出力結果を得るＭ個の加算器と、を備え、重み係数行列は、イメージング対象の注目物質の線減弱係数と背景物質の線減弱係数とに基づいて取得される。
本発明の実施例のＸ線検出方法とＸ線検出器は、複種のエネルギーウィンドウ信号を柔軟に取得でき、Ｘ線検出器の出力チャネルに対するニーズを低減し、出力データのビットレートを低減することができるので、効率がより良いのエネルギースペクトルＣＴ信号収集に用いられることができる。
本発明は、以下の図面を参照して行う本発明の実施形態の説明からよりよく理解される。
図１は本発明の一実施形態に係るＸ線検出方法のフローチャートである。 図２は本発明の一実施形態に係るＸ線検出器の概略構成図である。 図３Ａは、光子計数検出器を用いたＣＴ画像を示す図である。 図３Ｂは、本発明による実施例のＸ線検出器を用いたＣＴ画像を示す図である。 図４Ａは、光子計数検出器を用いたＣＴ画像を示す図である。 図４Ｂは、本発明による実施例のＸ線検出器を用いたＣＴ画像を示す図である。
以下、本発明の各方面の特徴及び例示の実施例をさらに詳細に説明する。以下の詳細の説明において、本発明を全面に理解するように、多くの特定の詳細が記載されている。しかし、本発明は、これらの詳細における一部がなくても実施できることは、当業者には明らかである。以下、実施例に対する説明は、単に本発明の例を示すことによって本発明をより明確に理解させるために提供される。本発明は、以下に説明される特定の構成及びアルゴリズムに限定されず、本発明の思想及び範囲から逸脱しないかぎり、素子、部材及びアルゴリズムのいかなる修正、交換、改進を含む。図面及び以下の説明では、本発明を不明瞭にすることを避けるために、公知の構成及び技術が示されていない。
本発明は、上記課題の少なくとも一つの課題を解決するために、新規なＸ線検出方法及びＸ線検出装置を提供する。以下、添付図面を参照して、本発明による実施例のＸ線検出方法及びＸ線検出器を詳細に説明する。
図１は、本発明による実施例のＸ線検出方法のフローチャートである。図１に示すように、当該Ｘ線検出方法は、以下のステップを備える。
Ｓ１０２では、Ｘ線源によって放射された光子のエネルギー範囲をＮ個のエネルギーウィンドウに分割し、ただし、Ｎは０より大きい整数である。ここで、Ｘ線源によって放射される光子のエネルギー範囲はEmin＜E＜Emaxであり、当該Ｎ個のエネルギーウィンドウの境界は[E0,En]であり、ここで、Emin≦E0＜En≦Emaxである。これらのエネルギーウィンドウは、等間隔で分けられてもよく、又は不等間隔で分けられてもよい。
Ｓ１０４では、イメージング対象の注目物質と背景物質との線減弱係数に基づいて、Ｎ個のエネルギーウィンドウのそれぞれの重み係数を取得する。ここで、イメージングのニーズに応じてイメージング対象の注目物質と背景物質が設定されてもよく、例えば、イメージング対象がヨード造影剤の濃度であり、注目物質としてヨウ素が選択され、背景物質として水が選択されてもよい。各エネルギーウィンドウの重み係数は、以下の三種の方法のいずれかにより得ることができる。
ａ）注目物質の線減弱係数をμ1(E)とし、背景物質の線減弱係数をμ2(E)とする。Ｎ個のエネルギーウインドウにおけるｉ（ｉは０より大きい任意の整数である）番目のエネルギーウインドウに対して、ｉ番目のエネルギーウインドウの重み係数ωiを以下のように算出することができる。
ここで、αは、物理的意味を持たず、[０，１]の範囲内の任意の数であってもよい。
つまり、注目物質の線減弱係数と背景物質の線減弱係数とをｉ番目のエネルギーウインドウのエネルギー範囲でそれぞれ積分してから、注目物質の線減弱係数の積分結果と背景物質の線減弱係数の積分結果とを利用し、加算処理によりｉ番目のエネルギーウインドウの重み係数を算出する。
ｂ）注目物質の線減弱係数をμ1(E)とし、背景物質の線減弱係数をμ2(E)とし、Ｘ線源によって放射されたＸ光が透過する注目物質の平均又は典型的な厚さをd1とし、Ｘ線源によって放射されたＸ光が透過する背景物質の平均又は典型的な厚さをd2とする。Ｎ個のエネルギーウインドウのｉ番目（ｉは０より大きい任意の整数である）のエネルギーウインドウ[Ei,Ei+1]に対して、ｉ番目のエネルギーウインドウωiの重み係数を以下のように算出することができる。
Φi(E)は、Ｘ線源によって放射されたＸ光のエネルギースペクトルである。
ｃ）注目物質の線減弱係数をμ1(E)とし、背景物質の線減弱係数をμ2(E)とし、Ｘ線源によって放射されたＸ光が透過する注目物質の平均又は典型的な厚さをd1とし、Ｘ線源によって放射されたＸ光が透過する背景物質の平均又は典型的な厚さをd2とする。Ｎ個のエネルギーウインドウのｉ番目（ｉは０より大きい任意の整数である）のエネルギーウインドウ[Ei,Ei+1]に対して、Ei=ei、Ei+1=ei+pとし、即ち、ｉ番目のエネルギーウインドウ[Ei,Ei+1]がＰ＋１個のエネルギーポイントを有し、これらのエネルギーポイントがそれぞれエネルギースペクトルΦi(ej)に対応し、ただし、ｊ＝０，１，...，ｐであるとるすと、ｉ番目のエネルギーウィンドウの重み係数ωiを次のように算出することができる。
ここで、αは任意の定数であってもよい。
つまり、ｂ）及びｃ）というエネルギーウィンドウの重み係数を取得する両種の態様では、注目物質の線減弱係数、背景物質の線減弱係数、Ｘ線源によって放射されたＸ光がｉ番目のエネルギーウインドウにあるエネルギースペクトル、Ｘ線源によって放射されたＸ光が透過する注目物質の平均厚さ、及びＸ線源によって放射されたＸ光が透過する背景物質の平均厚さに基づいて、ｉ番目のエネルギーウインドウの重み係数を算出する。
Ｓ１０６では、Ｎ個のエネルギーウインドウのそれぞれに対する重み係数に基づいて、Ｘ線検出器のＭ個の出力チャネルの重み係数行列（便宜のため、以下は重み係数行列Ａと呼ぶ）を取得する。ただし、Ｍは、０よりも大きい整数である。
いくつかの実施例では、Ｘ線検出器のｋ番目の出力チャネル（ｋは０より大きく且つＭ未満の任意の整数である）に対して、Ｎ個のエネルギーウインドウにおけるｉ番目のエネルギーウィンドウがｋ番目の出力チャネルに属すれば、即ち、ｉ番目のエネルギーウインドウがｋ番目の出力チャネルに対応するエネルギーウインドウの集合Ωｋに属すれば、以下の式により、ｉ番目のエネルギーウィンドウの重み係数ωiを正規化することで、重み係数行列Ａにおける、ｉ番目のエネルギーウィンドウおよびｋ番目の出力チャネルに対応する重み係数ak,iを算出することができる。
ここで、ωi maxは、ｋ番目の出力チャネルに対応するエネルギーウィンドウの集合Ωｋにおける最大の重み係数であり、round(g)は、ＩＮＴ関数を示し、Ｂは、ωiに対する定量精度の桁を示す。
ｉ番目のエネルギーウインドウがｋ番目の出力チャンネルに属さない場合、即ち、ｉ番目のエネルギーウインドウがｋ番目の出力チャンネルに対応するエネルギーウインドウの集合Ωｋに属さない場合、重み係数行列Ａにおける、ｉ番目のエネルギーウインドウおよびｋ番目の出力チャネルに対応する重み係数ak,iをゼロに設定することができる。
Ｓ１０８では、エネルギー範囲がＮ個のエネルギーウィンドウにおける各エネルギーウィンドウに入る光子の数及び重み係数行列Ａに基づいて、Ｘ線検出器のＭ個の出力チャネルの出力結果を取得する。
いくつかの実施例では、エネルギー範囲がＮ個のエネルギーウィンドウににおける各エネルギーウィンドウに入る光子の数で構成された列ベクトルｘと重み係数行列Ａとを乗算することで、Ｍ個の出力チャネルの出力結果を取得することができる。以下、列ペクトルｘと重み係数行列Ａを示して、Ｍ個の出力チャネルの出力結果をｙ＝ＡＸとする。
図２は、本発明の実施例に係るＸ線検出器の概略構成図であり、当該Ｘ線検出器は、上記の方法に基づいて設計されている。図２に示すように、当該Ｘ線検出器は、Ｘ線源によって放射されたエネルギー範囲が選別範囲に入る光子の数をカウントする（エネルギーウィンドウの選別器のそれぞれが上記のＮ個のエネルギーウィンドウの一つに対応する）Ｎ個のエネルギーウィンドウ選別器と、ｉ組目の乗算器におけるｋ個目の乗算器がＮ個のエネルギーウインドウ選別器におけるｉ番目のエネルギーウィンドウ選別器でカウントして得た光子の数とＭ行×Ｎ列の重み係数行列（即ち、上記重み係数行列Ａ）おけるｋ行目ｉ列目の重み係数とを乗算するＭ個の乗算器を一組として備えるＮ組の乗算器と、k個目の加算器が各組の乗算器におけるｋ個目の乗算器から得られる積を累算することで、Ｘ線検出器のＭ個の出力チャネルにおけるｋ番目の出力チャネルの出力結果を得るＭ個の加算器と、を備える。
いくつかの実施例では、図２に示すＸ線検出器は、更に、図１に記載されている方法を利用してＮ個のエネルギーウィンドウにおける何れかの一つのエネルギーウィンドウを算出する行列取得手段を備える。例えば、ｉ番目のエネルギーウィンドウωi、および重み係数行列Ａにおけるｉ番目のエネルギーウィンドウとｋ番目の出力チャネルに対応する重み係数ak,iがある。
現在の光子計数検出器は、各出力チャネルの出力結果をすべて出力する。出力チャネルが１２８個である場合、１２８個の結果が出力され、その後、１２８個の結果を利用して重み付けやその他の演算を行う。本発明の実施例に係るＸ線検出方法及びＸ線検出器は、実際のニーズに応じて出力チャンネルの数を柔軟に設定することができ、データ伝送帯域を節約することができる。ここで、重み係数はプログラマブルであるため、柔軟に使うことができる。
以下、上述したＸ線検出方法及びＸ線検出器を適用してＸ線イメージングを行う例を示す。
Ｘ線源から放射された光子のエネルギー範囲を２１ｋV〜５１ｋV（即ち、[２１，５１]）とし、対応するエネルギースペクトルをΦ(E)（即ち、各エネルギースに対応する光子の数がΦ(E)である）とし、７つのエネルギーウィンドウに分け、即ちＮ＝７とする。
Ｘ線検出器には、７つの出力チャネルが設けられ、即ちＭ=7であり、各出力チャネルに対応するエネルギーウィンドウの集合は、以下である。
注目物質（即ち、材料１）の線減弱係数がμ1(E)であり、背景物質（即ち、材料２）の線減弱係数がμ2(E)であり、Ｘ線源によって放射されるＸ線が透過する材料１の平均又は典型的な厚さはd1センチメートルであり、放射源から放射されるＸ線が透過する材料２の平均又は典型的な厚さはd2センチメートルである。
材料２は充填半径が２ｃｍであり、厚さが２.５ｃｍである円柱体であり、材料１は充填半径が０．５ｃｍであり、厚さがそれぞれ２ｃｍ、１．５ｃｍ、１ｃｍである３つの円柱体である。
上記７つのエネルギーウインドウにおけるｉ番目のエネルギーウインドウは、[Ei,Ei+1]、Ei=eikV、Ei+1=ei+1kVであり、即ち、当該エネルギーウインドウは、二つのエネルギーポイントを有し、対応するエネルギースペクトルがΦi(ej)であり、ただし、j＝０，１であり、ｉ番目のエネルギーウインドウの重み付け係数ωiは、以下の通りである。
ωi（ｉ＝１，２，...，７）のそれぞれに対して正規化すると、下記になる。
フィルタード・バックプロジェクション(ｆｉｌｔｅｒｅｄ ｂａｃｋ ｐｒｏｊｅｃｔｉｏｎ)法で画像を再構成する。本発明の実施例によるＸ線検出器を利用したＣＴ再構成画像（重み付けられた再構成画像）と重みが定数、即ち１である光子計数検出器を利用したＣＴ再構成画像（重み付けられていない再構成画像）に対して、信号対雑音比で定量に比較する。
例１
まず、注目物質（材料１）として乳腺組織を選択し、背景物質（材料２）として水を選択する。次に、選択されたエネルギースペクトルの範囲である２１ｋＶ〜５１ｋＶにおいて、各出力チャネルに対応するエネルギウィンドウを固定すると、それぞれ、チャネル１に対応するエネルギーウインドウ[２１，２２]、チャンネル２に対応するエネルギーウインドウ[２２，２７]、チャネル３に対応するエネルギーウインドウ[２７，３２]、チャネル４に対応するエネルギーウインドウ[３２，３７]、チャネル５に対応するエネルギーウインドウ[３７，４２]、チャネル６に対応するエネルギーウインドウ[４２，４７]、チャネル７に対応するエネルギーウインドウ[４７，５１]である。
このように、各チャネルに対応する重みは、それぞれ、チャネル１に対応する重みが１であり、チャネル２に対応する重みが０．７５であり、チャネル３に対応する重みが０．６８７５であり、チャネル４に対応する重みが０．５、チャネル５に対応する重みが０．５であり、チャネル６に対応する重みが０．４３７５であり、チャネル７に対応する重みが０．４３７５である。
光子計数検出器による再構成画像及び本発明の実施例のＸ線検出器による再構成画像における、注目物質及び背景物質の信号対雑音比をそれぞれ計算すると、それぞれ、光子計数検出器による再構成画像の信号対雑音比が２０２．５７であり、本発明の実施例のＸ線検出器による再構成画像の信号対雑音比が２０７．７５である。図３Ａ及び図３Ｂは、光子計数検出器を用いたＣＴ画像と、本発明による実施例のＸ線検出器を用いたＣＴ画像をそれぞれ示す。
例２
まず、注目物質（材料１）として乳腺組織を選択し、背景物質（材料２）として脂肪を選択する。次に、選択されたエネルギースペクトルの範囲である２１ｋＶ〜５１ｋＶにおいて、各出力チャネルに対応するエネルギウィンドウを固定すると、それぞれ、チャネル１に対応するエネルギーウインドウ[２１，２２]、チャンネル２に対応するエネルギーウインドウ[２２，２７]、チャネル３に対応するエネルギーウインドウ[２７，３２]、チャネル４に対応するエネルギーウインドウ[３２，３７]、チャネル５に対応するエネルギーウインドウ[３７，４２]、チャネル６に対応するエネルギーウインドウ[４２，４７]、チャネル７に対応するエネルギーウインドウ[４７，５１]である。
このように、各チャネルに対応する重みは、それぞれ、チャネル１に対応する重みが１であり、チャネル２に対応する重みが０．７５であり、チャネル３に対応する重みが０．６２５であり、チャネル４に対応する重みが０．５、チャネル５に対応する重みが０．４３７５であり、チャネル６に対応する重みが０．４３７５であり、チャネル７に対応する重みが０．３７５である。
光子計数検出器による再構成画像及び本発明の実施例のＸ線検出器による再構成画像における、注目物質及び背景物質の信号対雑音比をそれぞれ計算すると、それぞれ、光子計数検出器による再構成画像の信号対雑音比が２１９．８５であり、本発明の実施例のＸ線検出器による再構成画像の信号対雑音比が２２７．４５である。図４Ａ及び図４Ｂは、光子計数検出器を用いたＣＴ画像と、本発明による実施例のＸ線検出器を用いたＣＴ画像をそれぞれ示す。
以上の説明から分かるように、本発明の実施例に係るＸ線検出方法及びＸ線検出装置を使用すると、複数の種類のエネルギーウインドウ信号を柔軟に得ることができ、出力信号の信号対雑音比を向上させることができる一方、出力チャネルの数を減らすことで、出力データのビットレートを低減することができ、効率の良いエネルギースペクトルＣＴ信号の収集に適合する。また、イメージングのニーズに応じて、重みを便利に調整して、ジョブによる信号収集の最適化を実現でき、例えば、エネルギーウィンドウ選別器によって散乱光子を判別し、比例によってその重みを低減させることで、散乱による影響を低減又は排除することができる。
上述の行列取得手段は、ハードウェア、ソフトウェア、ファームウェア、又はそれらの組み合わせで実現されてもよい。ハードウェアで実現される場合、例えば、電子回路、特定用途向け集積回路（ＡＳＩＣ）、適切なファームウェア、プラグイン、機能カードなどである。ソフトウェアで実現される場合、本発明の要素は、必要なタスクを実行するためのプログラム又はコードセグメントである。プログラム又はコードセグメントは、機械で読み取り可能な媒体に格納されてもよく、又はキャリアで伝送されるデータ信号を利用して伝送媒体又は通信リンクを介して伝送されてもよい。「機械で読み取り可能な媒体」は、情報を記憶又は転送可能な任意の媒体を含むことができる。機械で読み取り可能な媒体の例には、電子回路、半導体メモリデバイス、ＲＯＭ、フラッシュメモリ、消去可能ＲＯＭ（ＥＲＯＭ）、フレキシブルディスク、ＣＤ−ＲＯＭ、光ディスク、ハードディスク、光ファイバメディア、高周波（ＲＦ）リンクなどが含まれる。コードセグメントは、インターネット、イントラネットなどのコンピュータネットワークを介してダウンロードすることができる。
本発明は、その思想及び本質的な特徴から逸脱することなく、他の特定の形態で実施することができる。例えば、特定の実施例に記載されたアルゴリズムが修正されても、システムアーキテクチャは本発明の基本的な精神から逸脱しない。従って、目前の実施例はあくまでも例示的なものであり、限定的なものではなく、本発明の範囲は、上記に記載された内容に限定されるのではなく、添付した特許請求の範囲により限定される。また、特許請求の範囲及び均等物の範囲内のすべての変更は、本発明の範囲に含まれる。

Claims (10)

  1. Ｘ線源によって放射された光子のエネルギー範囲をＮ（Ｎは０より大きい整数である）個のエネルギーウィンドウに分割するステップと、
    イメージング対象の注目物質と背景物質との線減弱係数に基づいて、前記Ｎ個のエネルギーウインドウのそれぞれの重み係数を取得するステップと、
    前記Ｎ個のエネルギーウィンドウのそれぞれの重み係数に基づいて、Ｘ線検出器のＭ（Ｍは０より大きい整数である）個の出力チャネルの重み係数行列を取得するステップと、
    エネルギー範囲が前記Ｎ個のエネルギーウィンドウのそれぞれに入る光子の数および前記重み係数行列に基づいて、前記Ｍ個の出力チャネルの出力結果を取得するステップと
    を備えるＸ線検出方法。
  2. 前記Ｍ個の出力チャネルにおけるｋ番目の出力チャネルについて、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウインドウが前記ｋ番目の出力チャンネルに属する場合、
    前記ｉ番目のエネルギーウィンドウの重み係数を正規化することによって、前記重み係数行列における、前記ｉ番目のエネルギーウィンドウ及び前記ｋ番目の出力チャネルに対応する重み係数を取得し、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウインドウが前記ｋ番目の出力チャンネルに属さない場合、
    前記重み係数行列における、前記ｉ番目のエネルギーウインドウ及び前記ｋ番目の出力チャネルに対応する重み係数を０に設定し、ただし、ｋ及びｉは両方とも０より大きい整数である、
    請求項１に記載のＸ線検出方法。
  3. エネルギー範囲が前記Ｎ個のエネルギーウィンドウのそれぞれに入る光子の数で構成される列ペクトルと前記重み係数行列とを乗算して、前記Ｍ個の出力チャンネルの出力結果を取得する
    請求項１に記載のＸ線検出方法。
  4. 前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウィンドウについて、
    前記注目物質の線減弱係数および前記背景物質の線減弱係数に対して、前記ｉ番目のエネルギーウインドウのエネルギー範囲でそれぞれ積分し、
    前記注目物質の線減弱係数の積分結果と前記背景物質の線減弱係数の積分結果を利用して、加算により前記ｉ（ｉが０より大きい整数である）番目のエネルギーウインドウの重み係数を算出する
    請求項１に記載のＸ線検出方法。
  5. 前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウィンドウについて、
    前記注目物質の線減弱係数と、前記背景物質の線減弱係数と、前記Ｘ線源によって放射されたＸ線の前記ｉ番目のエネルギーウィンドウにおけるエネルギースペクトルと、及び前記Ｘ線源によって放射されたＸ線が透過する前記注目物質の平均厚さと、前記Ｘ線源によって放射されたＸ線が透過する前記背景物質の平均厚さに基づいて、前記ｉ（ｉは、０より大きい整数である）番目のエネルギーウィンドウの重み係数を算出する
    請求項１に記載のＸ線検出方法。
  6. それぞれ、Ｘ線源によって放射されたエネルギー範囲が選別範囲に入る光子の数をカウントするＮ個のエネルギーウィンドウ選別器と、
    ｉ組目の乗法器におけるｋ番目の乗算器が、前記Ｎ個のエネルギーウィンドウ選別器におけるｉ番目のエネルギーウィンドウ選別器でカウントして得た光子の数とＭ行×Ｎ列の重み係数行列におけるｋ行目ｉ列目の重み係数とを乗算する、それぞれＭ個の乗算器を備えるＮ（Ｍ、Ｎ、ｋ、ｉのいずれもは０より大きい整数である）組の乗算器と、
    ｋ番目の加算器が、各組の乗算器におけるｋ番目の乗算器で得た積を累算して、前記Ｘ線検出器のＭ個の出力チャネルにおけるｋ番目の出力チャネルの出力結果を得るＭ個の加算器と、を備え、
    前記重み係数行列は、イメージング対象の注目物質の線減弱係数と背景物質の線減弱係数とに基づいて取得される
    ことを特徴とするＸ線検出器。
  7. 前記注目物質の線減弱係数と前記背景物質の線減弱係数に基づき、前記Ｎ個エネルギーウィンドウにおける各エネルギーウィンドウの重み係数を取得し、且つ、前記Ｎ個のエネルギーウィンドウにおける各エネルギーウィンドウの重み係数に基づいて前記重み係数行列を取得する行列取得手段を更に備える
    請求項６に記載のＸ線検出器。
  8. 前記行列取得手段は、
    前記Ｍ個の出力チャネルにおけるｋ番目の出力チャネルについて、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウインドウが前記ｋ番目の出力チャンネルに属する場合、
    前記ｉ番目のエネルギーウィンドウの重み係数を正規化することによって、前記重み係数行列における、前記ｉ番目のエネルギーウィンドウ及び前記ｋ番目の出力チャネルに対応する重み係数を取得し、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウインドウが前記ｋ番目の出力チャンネルに属さない場合、
    前記重み係数行列における、前記ｉ番目のエネルギーウインドウ及び前記ｋ番目の出力チャネルに対応する重み係数が０に設定される、
    ように構成されている
    請求項７に記載のＸ線検出器。
  9. 前記行列取得手段は、さらに、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウィンドウについて、
    前記注目物質の線減弱係数および前記背景物質の線減弱係数に対して、前記ｉ番目のエネルギーウインドウのエネルギー範囲でそれぞれ積分し、
    前記注目物質の線減弱係数の積分結果と前記背景物質の線減弱係数の積分結果を利用して、加算により前記ｉ番目のエネルギーウインドウの重み係数を算出する
    ように構成されている
    請求項７に記載のＸ線検出器。
  10. 前記行列取得手段は、さらに、
    前記Ｎ個のエネルギーウィンドウにおけるｉ番目のエネルギーウィンドウについて、
    前記注目物質の線減弱係数と、前記背景物質の線減弱係数と、前記Ｘ線源によって放射されたＸ線の前記ｉ番目のエネルギーウィンドウにおけるエネルギースペクトルと、及び前記Ｘ線源によって放射されたＸ線が透過する前記注目物質の平均厚さと、前記Ｘ線源によって放射されたＸ線が透過する前記背景物質の平均厚さに基づいて、前記ｉ（ｉは、０より大きい整数である）番目のエネルギーウィンドウの重み係数を算出する
    ように構成されている
    請求項７に記載のＸ線検出器。
JP2018545666A 2016-05-05 2017-02-03 Ｘ線検出方法及びｘ線検出器 Active JP6630841B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201610294662.3A CN107345923B (zh) 2016-05-05 2016-05-05 X射线探测方法、和x射线探测器
CN201610294662.3 2016-05-05
PCT/CN2017/072832 WO2017190538A1 (zh) 2016-05-05 2017-02-03 X射线探测方法、和x射线探测器

Publications (2)

Publication Number Publication Date
JP2019513225A true JP2019513225A (ja) 2019-05-23
JP6630841B2 JP6630841B2 (ja) 2020-01-15

Family

ID=60202699

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018545666A Active JP6630841B2 (ja) 2016-05-05 2017-02-03 Ｘ線検出方法及びｘ線検出器

Country Status (6)

Country Link
US (1) US10724969B2 (ja)
EP (1) EP3454048B1 (ja)
JP (1) JP6630841B2 (ja)
CN (1) CN107345923B (ja)
AU (1) AU2017259380B2 (ja)
WO (1) WO2017190538A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022508691A (ja) * 2018-10-12 2022-01-19 ケアストリーム デンタル エルエルシー 再構成された画像における歯をセグメント化する方法
EP3667370A1 (en) 2018-12-11 2020-06-17 Koninklijke Philips N.V. Energy weighting of photon counts for conventional imaging
CN109801243B (zh) * 2019-01-28 2020-11-27 中国科学院高能物理研究所 材料成像权重的确定方法、装置、介质及电子设备
CN109875596A (zh) * 2019-03-27 2019-06-14 毕节市七星关区人民医院 一种针对乳腺专用ct设备的乳腺模体
DE102019111567A1 (de) * 2019-05-03 2020-11-05 Wipotec Gmbh Verfahren und Vorrichtung zur Röntgeninspektion von Produkten, insbesondere von Lebensmitteln
JP7440271B2 (ja) * 2020-01-10 2024-02-28 富士フイルムヘルスケア株式会社 放射線撮像装置および光子計数型検出器の較正方法
CN115097513B (zh) * 2022-06-20 2023-03-21 上海市计量测试技术研究院 能谱计数窗口的定位方法、装置、电子设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0385482A (ja) * 1989-08-30 1991-04-10 Shimadzu Corp シンチレーションカメラ
JP2002524126A (ja) * 1998-09-08 2002-08-06 ヴェリタス・ファーマソイティカルズ 組織の放射線撮影画像形成システムおよび方法
JP2004512512A (ja) * 2000-10-25 2004-04-22 フィリップス　メディカル　システムズ　（クリーヴランド），インコーポレイテッド 特徴が強調された透過撮像を用いた核撮像のシステム及び方法
JP2006101926A (ja) * 2004-09-30 2006-04-20 M & C:Kk 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法
JP2011503570A (ja) * 2007-11-09 2011-01-27 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 減衰マップを形成するための装置及び方法
JP2014158713A (ja) * 2013-02-19 2014-09-04 Toshiba Corp コンピュータ断層撮影装置、校正プログラム、及び光子数校正装置
JP2015200649A (ja) * 2014-04-07 2015-11-12 コミッサリア ア レネルジー アトミーク エ オ ゼネルジ ザルタナテイヴ 分光センサによって取得された減衰量測定値を線形化する方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6373919B1 (en) * 1999-08-10 2002-04-16 Ge Yokogawa Medical Systems, Limited Projection data correction method and apparatus, and radiation tomographic imaging method and apparatus
DE102004029010A1 (de) * 2004-06-16 2006-01-19 Siemens Ag Vorrichtung und Verfahren für die Streustrahlungskorrektur in der Projektionsradiographie, insbesondere der Mammographie
US7583779B2 (en) * 2004-11-24 2009-09-01 General Electric Company System and method for acquisition and reconstruction of contrast-enhanced, artifact-reduced CT images
DE102005027436B4 (de) * 2005-06-14 2008-09-04 Siemens Ag Verfahren zur Berechnung von absorberspezifischen Gewichtungskoeffizienten und Verfahren zur Verbesserung eines von einem Absorber abhängigen Kontrast-zu-Rausch-Verhältnisses in einem von einer Röntgeneinrichtung erzeugten Röntgenbild eines zu untersuchenden Objektes
US7473900B2 (en) * 2005-06-16 2009-01-06 Siemens Medical Solutions Usa, Inc. Acquisition window compensation for nuclear medical image reconstruction attenuation coefficient maps
WO2007029129A2 (en) * 2005-09-06 2007-03-15 Koninklijke Philips Electronics, N.V. Data handling and analysis in computed tomography with multiple energy windows
JP5049937B2 (ja) * 2008-09-30 2012-10-17 株式会社日立製作所 エネルギー情報を用いたｘ線ｃｔ画像再構成法
US9128584B2 (en) * 2013-02-15 2015-09-08 Carl Zeiss X-ray Microscopy, Inc. Multi energy X-ray microscope data acquisition and image reconstruction system and method
WO2014176328A1 (en) * 2013-04-23 2014-10-30 Virginia Tech Intellectual Properties, Inc. Hybrid detector modules and dynamic thresholding for spectral ct
CN103472074B (zh) * 2013-06-19 2016-01-20 清华大学 Ct成像系统和方法
CN104422704B (zh) * 2013-08-21 2017-07-25 同方威视技术股份有限公司 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法
CN104215997B (zh) * 2014-09-03 2016-09-21 北京市射线应用研究中心 一种α和γ双通道能谱连续监测方法及监测系统
CN105092617A (zh) * 2015-09-18 2015-11-25 重庆大学 一种基于x射线能谱ct和x射线荧光ct技术的双模态分子成像系统
CN105510364B (zh) * 2015-12-30 2019-02-19 苏州科耐视智能科技有限公司 基于x射线的工业部件缺陷的无损检测系统及其检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0385482A (ja) * 1989-08-30 1991-04-10 Shimadzu Corp シンチレーションカメラ
JP2002524126A (ja) * 1998-09-08 2002-08-06 ヴェリタス・ファーマソイティカルズ 組織の放射線撮影画像形成システムおよび方法
JP2004512512A (ja) * 2000-10-25 2004-04-22 フィリップス　メディカル　システムズ　（クリーヴランド），インコーポレイテッド 特徴が強調された透過撮像を用いた核撮像のシステム及び方法
JP2006101926A (ja) * 2004-09-30 2006-04-20 M & C:Kk 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法
JP2011503570A (ja) * 2007-11-09 2011-01-27 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 減衰マップを形成するための装置及び方法
JP2014158713A (ja) * 2013-02-19 2014-09-04 Toshiba Corp コンピュータ断層撮影装置、校正プログラム、及び光子数校正装置
JP2015200649A (ja) * 2014-04-07 2015-11-12 コミッサリア ア レネルジー アトミーク エ オ ゼネルジ ザルタナテイヴ 分光センサによって取得された減衰量測定値を線形化する方法

Also Published As

Publication number Publication date
CN107345923A (zh) 2017-11-14
EP3454048B1 (en) 2020-12-09
US20190086562A1 (en) 2019-03-21
AU2017259380A1 (en) 2018-11-22
CN107345923B (zh) 2020-05-19
US10724969B2 (en) 2020-07-28
JP6630841B2 (ja) 2020-01-15
EP3454048A1 (en) 2019-03-13
EP3454048A4 (en) 2020-02-19
WO2017190538A1 (zh) 2017-11-09
AU2017259380B2 (en) 2019-07-18

Similar Documents

Publication Publication Date Title
JP2019513225A (ja) Ｘ線検出方法及びｘ線検出器
US9861324B2 (en) Hybrid detector modules and dynamic thresholding for spectral CT
EP2537048B1 (en) Improving image quality in photon counting-mode detector systems
Tapiovaara et al. SNR and noise measurements for medical imaging: I. A practical approach based on statistical decision theory
JP7532464B2 (ja) 光子計数に基づくｘ線検出システム
Xu et al. Cascaded systems analysis of photon counting detectors
US10117628B2 (en) Photon counting apparatus
US7480362B2 (en) Method and apparatus for spectral computed tomography
Brombal et al. Large-area single-photon-counting CdTe detector for synchrotron radiation computed tomography: a dedicated pre-processing procedure
US8965095B2 (en) Noise balance pre-reconstruction data decomposition in spectral CT
WO2020162296A1 (ja) 画像処理装置および画像処理方法
EP3631763B1 (en) Method and devices for image reconstruction
US20160084968A1 (en) Direct conversion radiation detector digital signal processing electronics
CN110175957A (zh) 多能量ct基材料物质分解方法
Liu et al. Spectral response model for a multibin photon-counting spectral computed tomography detector and its applications
Wang Correction for the flux measurement bias in X-Ray source detection
Feng et al. Neural-networks-based photon-counting data correction: pulse pileup effect
Taguchi Assessment of Multienergy Interpixel Coincidence Counters (MEICC) for charge sharing correction or compensation for photon counting detectors with boxcar signals
JP6479772B2 (ja) Ｘ線撮像装置及び方法
CN106056541B (zh) 在多谱ct成像中的适应于结构的噪声抑制
Jumanazarov et al. Significance of the spectral correction of photon counting detector response in material classification from spectral x-ray CT
Atharifard et al. Per-pixel energy calibration of photon counting detectors
Avila et al. Contrast cancellation technique applied to digital x‐ray imaging using silicon strip detectors
Odinaka et al. Spectrally grouped total variation reconstruction for scatter imaging using ADMM
US11506487B2 (en) X-ray imaging method and system thereof

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180830

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190730

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190829

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191209

R150 Certificate of patent or registration of utility model

Ref document number: 6630841

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