JP2022010983A - Ｘ線計測装置の校正方法 - Google Patents

Ｘ線計測装置の校正方法 Download PDF

Info

Publication number
JP2022010983A
JP2022010983A JP2020111821A JP2020111821A JP2022010983A JP 2022010983 A JP2022010983 A JP 2022010983A JP 2020111821 A JP2020111821 A JP 2020111821A JP 2020111821 A JP2020111821 A JP 2020111821A JP 2022010983 A JP2022010983 A JP 2022010983A
Authority
JP
Japan
Prior art keywords
ray
rotation
measuring device
rotary table
ray source
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.)
Pending
Application number
JP2020111821A
Other languages
English (en)
Inventor
正人 今
Masato Kon
洋武 真家
Hiromu Maie
誠治 佐々木
Seiji Sasaki
常太 宮倉
Jota Miyakura
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.)
Mitutoyo Corp
Mitsutoyo Kiko Co Ltd
Original Assignee
Mitutoyo Corp
Mitsutoyo Kiko Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitutoyo Corp, Mitsutoyo Kiko Co Ltd filed Critical Mitutoyo Corp
Priority to JP2020111821A priority Critical patent/JP2022010983A/ja
Priority to US17/360,304 priority patent/US11573190B2/en
Priority to DE102021003348.5A priority patent/DE102021003348A1/de
Priority to CN202110724806.5A priority patent/CN113932741A/zh
Publication of JP2022010983A publication Critical patent/JP2022010983A/ja
Pending legal-status Critical Current

Links

Images

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/06—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 measuring the absorption
    • G01N23/083—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 measuring the absorption the radiation being X-rays
    • 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
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • G01B21/042—Calibration or calibration artifacts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • G01B15/04—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring contours or curvatures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • G01B21/047—Accessories, e.g. for positioning, for tool-setting, for measuring probes
    • 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/06—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 measuring the absorption
    • G01N23/18—Investigating the presence of flaws defects or foreign matter
    • 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/30—Accessories, mechanical or electrical features
    • G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Toxicology (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pulmonology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Electromagnetism (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Length-Measuring Devices Using Wave Or Particle Radiation (AREA)

Abstract

【課題】被測定物を回転可能に載置する回転テーブルの回転中心位置を単純な工程で容易に算出可能とする。【解決手段】校正治具１０２を回転テーブル１２０に載置する載置工程と、Ｘ線１１８を校正治具１０２に照射しＸ線画像検出器１２４の出力から、重心位置ＩｍＰｏｓ（１～Ｑ）＿Ｓｐｈｒ＿(１～Ｎ)を特定する特徴位置算出工程と、重心位置ＩｍＰｏｓ（１～Ｑ）＿Ｓｐｈｒ＿（１～Ｎ）と既知の相対位置間隔Ｐｕ、Ｐｖから、射影変換行列Ｈｋ（ｋ＝１～Ｑ）を算出する変換行列算出工程と、回転テーブル１２０を所定角度αで２回以上回転させ、特徴位置算出工程と変換行列算出工程とを繰り返す回転検出工程と、射影変換行列Ｈｋ（ｋ＝１～Ｑ）に基づいて回転テーブル１２０の回転中心位置Ｃｐを算出する中心位置算出工程と、を含む。【選択図】図４

Description

本発明は、Ｘ線計測装置の校正方法に係り、特に、被測定物を回転可能に載置する回転テーブルの回転中心位置を単純な工程で容易に算出可能なＸ線計測装置の校正方法に関する。
従来、Ｘ線計測装置（計測用Ｘ線ＣＴ装置）は、Ｘ線を用いて被測定物を三次元形状計測することができ、外観からでは確認困難な鋳物部品の鬆、溶接部品の溶接不良、および電子回路部品の回路パターンの欠陥など、主に観察・検査に用いられてきた。しかし、近年、３Ｄプリンタの普及も手伝い、加工品内部の３Ｄ寸法計測とその高精度化の需要が増大しつつある。このような需要に対して、Ｘ線計測装置にさらなる寸法計測の高精度化が望まれている。
Ｘ線計測装置における寸法計測をより高精度に実施するためには、特許文献１に記載されているように、測定開始前に装置固有の各種校正を行うことが重要となっている。
特開２０００－２９８１０５号公報
しかしながら、特許文献１に記載の校正方法では、例えば、回転テーブルの回転中心位置が、校正治具を１８０度回転させてＸ線検出器の検出面上の回転中心に対応する中心チャンネルを算出したうえで、それとは別の工程によりＸ線源と回転中心位置との距離を算出しなければならなかった。即ち、従来、回転テーブルの回転中心位置の算出は複雑で困難であった。
本発明は、前記従来の問題点を解決するべくなされたもので、被測定物を回転可能に載置する回転テーブルの回転中心位置を単純な工程で容易に算出可能とするＸ線計測装置の校正方法を提供することを課題とする。
本願の請求項１に係る発明は、Ｘ線を用いて被測定物を三次元形状計測するＸ線計測装置の校正方法であって、前記Ｘ線計測装置が、前記Ｘ線を発生させるＸ線源と、前記被測定物を回転可能に載置する回転テーブルと、該被測定物を透過した該Ｘ線を検出するＸ線画像検出器と、を備え、該Ｘ線画像検出器への投影像によって特定可能な形状の基準物体を既知の相対位置間隔でＮ箇所（Ｎ≧４）以上に配置可能な校正治具を、前記回転テーブルに載置する載置工程と、前記Ｘ線を該校正治具に照射し該Ｘ線画像検出器の出力から、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置を特定する特徴位置算出工程と、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置と前記既知の相対位置間隔から、該基準物体の前記Ｘ線画像検出器の検出面への射影変換を行う変換行列を算出する変換行列算出工程と、前記回転テーブルを所定角度で２回以上回転させ、前記特徴位置算出工程と前記変換行列算出工程とを繰り返す回転検出工程と、前記変換行列に基づいて前記回転テーブルの回転中心位置を算出する中心位置算出工程と、を含むことにより、前記課題を解決したものである。
本願の請求項２に係る発明は、前記校正治具において前記基準物体の全てが１つの平面上にのみ載置されている場合には、前記変換行列を射影変換行列とし、該基準物体が三次元的に載置されている場合には、Ｎ＝６とし且つ該変換行列を射影行列としたものである。
本願の請求項３に係る発明は、前記中心位置算出工程で、更に、前記回転テーブルの回転軸を算出するようにしたものである。
本願の請求項４に係る発明は、前記中心位置算出工程で、前記回転テーブルの代わりに前記Ｘ線源と前記Ｘ線画像検出器とが回転したと想定し、前記変換行列に基づいて前記所定角度の回転毎の該Ｘ線源の絶対位置を算出することで、前記回転テーブルの回転中心位置を算出するようにしたものである。
本願の請求項５に係る発明は、前記回転テーブルを前記所定角度で４回以上回転させ、前記Ｘ線源の絶対位置を算出する際に、該Ｘ線源と前記Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを変数とし、前記変換行列に基づいて算出される該Ｘ線源の絶対位置を真円にフィッティングした仮真円の軌跡上の位置と該Ｘ線源の絶対位置との距離誤差を評価することで、該Ｘ線源と該Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを算出するようにしたものである。
本願の請求項６に係る発明は、前記中心位置算出工程で、前記Ｘ線源の絶対位置の変化から、真円でフィッティングされた軌跡の中心位置を算出し、該中心位置を前記回転テーブルの回転中心位置とするようにしたものである。
本願の請求項７に係る発明は、更に、前記回転テーブルの回転軸を算出する際には、前記軌跡の水平面からの傾斜角度を算出し、該傾斜角度と前記回転中心位置とから該回転軸を算出するようにしたものである。
本願の請求項８に係る発明は、前記基準物体を、球としたものである。
本願の請求項９に係る発明は、前記基準物体の投影像の特徴点の位置を、該投影像の重心位置としたものある。
本発明によれば、被測定物を回転可能に載置する回転テーブルの回転中心位置を単純な工程で容易に算出することが可能となる。
本発明の実施形態に係るＸ線計測装置の基本的な構成を示す概略側面図 図１のＸ線計測装置を主要部のみを示す概略上面図 図１の校正治具を示す図（正面図（Ａ）、上面図（Ｂ）） 本発明の実施形態に係るＸ線計測装置の校正手順を示すフロー図 図４における球の絶対位置を算出する工程の詳細フロー図 図５におけるＸ線源の絶対位置の算出をした後に、Ｘ線源とＸ線画像検出器との距離と、Ｘ線源からのＸ線画像検出器への垂線の足の位置とを算出するフロー図 球の絶対位置とＸ線源の絶対位置との関係を示す図
以下、図面を参照して、本発明の実施の形態について詳細に説明する。なお、本発明は以下の実施形態及び実施例に記載した内容により限定されるものではない。又、以下に記載した実施形態及び実施例における構成要件には、当業者が容易に想定できるもの、実質的に同一のもの、いわゆる均等の範囲のものが含まれる。更に、以下に記載した実施形態及び実施例で開示した構成要素は適宜組み合わせてもよいし、適宜選択して用いてもよい。
図１に本発明の第１実施形態を示す。なお、図１では、紙面に対して左右方向をｚ軸方向、紙面に対して上下方向をｙ軸方向とし、紙面に対して垂直な方向をｘ軸方向として説明する。
Ｘ線計測装置１００は、Ｘ線を用いて被測定物を三次元形状計測する装置であり、図１に示す如く、本体部１０８、ホストコンピュータ１２８およびモーションコントローラ１３０を備えている。
なお、図１、図２では校正治具１０２が被測定物の代わりに回転テーブル１２０上に載置されている。校正治具１０２は、図３（Ａ）、（Ｂ）に示す如く、Ｘ線１１８の透過可能な材質（例えばアルミニウム等）でできており、板状部材１０４上に直径Ｄの球（基準物体）１０６を既知の相対位置間隔Ｐｕ、Ｐｖで複数（例えば４個＊３個でＮ＝１２個）備えている（即ち、球１０６は既知の相対位置間隔Ｐｕ、Ｐｖで１２箇所に配置されている）。つまり、校正治具１０２において球１０６の全てが１つの平面上にのみ載置されている。同時に、１２個の球１０６、言い換えれば、１２箇所にある球１０６の相対位置Ｘ（１～１２）は既知であるともいえる（Ｘ（１～１２）とＸ１～Ｘ１２とは同じことをいう、以下同様の表記をする）。なお、球１０６は、形状がシンプルであり、Ｘ線画像検出器１２４への投影像によって容易に特定可能な形状である。なお、図３（Ａ）では、紙面に対して左右方向をｕ軸方向、紙面に対して上下方向をｖ軸方向とし、紙面に対して垂直な方向をｗ軸方向として説明する。
前記本体部１０８は、図１に示す如く、ベース１１２上に、Ｘ線１１８の漏れを防ぐＸ線遮蔽カバー１１０と、Ｘ線１１８を発生させるＸ線源１１６と、被測定物（図示せず）を回転可能に載置する回転テーブル１２０と、被測定物を透過したＸ線１１８を検出するＸ線画像検出器１２４と、を備えている。Ｘ線源１１６は、ベース１１２上の線源支持台１１４に設けられている。線源支持台１１４は、Ｘ線源１１６をｘｙｚの３軸方向へ移動可能とする直動機構を備えることができる。回転テーブル１２０は、ベース１１２上のテーブル支持台１２２に設けられている。なお、テーブル支持台１２２は、被測定物をｘｙｚの３軸方向へ移動可能とする直動機構を備える。更には、テーブル支持台１２２に、回転テーブル１２０の回転軸Ａｘを傾斜調整可能とする傾斜機構が設けられていてもよい。Ｘ線画像検出器１２４は、Ｘ線１１８に感度がある二次元の検出面１２４Ａを有する。Ｘ線画像検出器１２４は、ベース１１２上の検出器支持台１２６に支持されている。検出器支持台１２６も、Ｘ線画像検出器１２４をｘｙｚの３軸方向へ移動可能とする直動機構を備えてもよい。Ｘ線源１１６からのＸ線１１８の放射ビームは、ｚ軸方向に円錐状に広がり、その中心線が回転テーブル１２０の回転軸Ａｘと交差し、Ｘ線画像検出器１２４の検出面１２４Ａの垂線となるように調整される。
図１に示す前記ホストコンピュータ１２８は、本体部１０８の線源支持台１１４、Ｘ線源１１６、回転テーブル１２０、テーブル支持台１２２、Ｘ線画像検出器１２４及び検出器支持台１２６を制御する。また、ホストコンピュータ１２８は、図示せぬ記憶部に格納されたプログラムを読み出して実行することで、Ｘ線計測装置１００の計測動作および校正を自動あるいは半自動で行うこともできる。つまり、ホストコンピュータ１２８は、Ｘ線計測装置１００の計測動作において、例えば、Ｘ線画像検出器１２４で得られた投影像のデータを再構成して、被測定物の三次元ボリュームデータを作成する。
そして、ホストコンピュータ１２８は、Ｘ線計測装置１００の校正において、例えば、以下に示す一連の式で示される関係を用いて、ｋ（１≦ｋ≦Ｑ；Ｑ≧３）番目の回転テーブル１２０の回転位置Ｐｏｓｋにおける３行＊３列の射影変換行列（変換行列）Ｈｋから、回転テーブル１２０の代わりにＸ線源１１６が回転したと想定して、ｋ番目の想定回転位置におけるＸ線源１１６の絶対位置Ｘｓを算出することができる（図７（Ａ）、（Ｂ））。そして、ホストコンピュータ１２８は、回転テーブル１２０の回転中心位置Ｃｐを算出する（即ち、ホストコンピュータ１２８は、射影変換行列Ｈｋに基づいて回転テーブル１２０の回転中心位置Ｃｐを算出する）。
具体的に説明すると、まずは、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とから、式（１）で示す内部パラメータ行列Ａが定義される。なお、内部パラメータ行列Ａの１行１列目の距離ｆと２行２列目の距離ｆとは、Ｘ線画像検出器１２４の画素の縦横比が異なる場合には若干値が異なることとなる。また、内部パラメータ行列Ａの１行２列目には画像の歪みに関わるスキューＳが用いられることもあるが、本実施形態では、スキューＳを０としている。
Figure 2022010983000002
このとき、ｋ番目の想定回転位置の回転行列Ｒｋは、３つの列ベクトルｒｋ１、ｒｋ２、ｒｋ３からなる３行＊３列であり、式（２）の如く表すことができる。
Ｒｋ＝［ｒｋ１ ｒｋ２ ｒｋ３］ （２）
ここで、ｋ番目の想定回転位置の並進行列Ｔｋ（１列の並進ベクトル）、式（１）、（２）を用いることで、射影変換行列Ｈｋは式（３）のように分解することができる。
Ｈｋ＝Ａ［ｒｋ１ ｒｋ２ Ｔｋ］ （３）
式（３）より、式（４）に示すように、ｋ番目の想定回転位置におけるＸ線源１１６の絶対位置Ｘｓを算出することができる。なお、符号－ｉｎｖ（ ）は逆行列を示す。
Ｘｓ＝－ｉｎｖ（Ｒｋ）＊Ｔｋ （４）
更に、ホストコンピュータ１２８は、内部パラメータ行列Ａや射影変換行列Ｈｋを使って、球１０６の絶対位置Ｘａ（１～Ｎ）を算出することもできる。イメージ的には、球１０６の絶対位置Ｘａ（１～Ｎ）は、球１０６の既知の相対位置間隔Ｐｕ、Ｐｖや、距離ｆや、位置Ｃｃ（ｃｘ，ｃｙ）の情報を使って空間上にそれぞれの球１０６を配置し、射影変換行列Ｈｋに含まれる回転・並進情報を使って実際の位置へ座標変換することで算出することができる。
なお、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とが不明である際には、ホストコンピュータ１２８は、距離ｆと位置Ｃｃ（ｃｘ，ｃｙ）とを変数として、それらに適当な値を代入し式（４）を用いて、即ち射影変換行列Ｈｋ（ｋ＝１～Ｑ）に基づいてＱ箇所のＸ線源１１６の絶対位置Ｘｓを算出する。ホストコンピュータ１２８は、算出したＱ箇所のＸ線源１１６の絶対位置Ｘｓから例えば最小二乗法で真円（仮真円）をフィッテイングさせる。そして、ホストコンピュータ１２８は、仮真円の軌跡上の位置とＱ箇所のＸ線源１１６の絶対位置Ｘｓとの距離誤差を評価し、最も距離誤差が少ない場合の距離ｆと位置Ｃｃ（ｃｘ，ｃｙ）とを算出する。なお、このときには、回転位置の総数Ｑは４以上、即ち、回転テーブル１２０を所定角度αで４回以上回転させることとなる。
図１に示す前記モーションコントローラ１３０は、ホストコンピュータ１２８に接続され、本体部１０８のＸ線源１１６や回転テーブル１２０の回転・移動、及び各種機構を制御する。
次に、Ｘ線計測装置１００の測定動作を簡略に説明する。
測定に当たっては、Ｘ線１１８を発生させた状態で回転テーブル１２０上の被測定物を回転させ、複数の角度方向（例えば角度分割数１０００～６０００程度）から投影像を収集する。収集された投影像は、被測定物を水平に横断するスライス面を基準面として再構成処理がなされ、被測定物の三次元ボリュームデータ（三次元像）を作るようになされている。
次に、Ｘ線計測装置１００の校正手順を図４から図７を用いて説明する。ここでは、ホストコンピュータ１２８で全ての演算がなされている。なお、例えばｋ＝１のときにｋ番目の回転位置Ｐｏｓｋが回転位置Ｐｏｓ１を示す。また、ｋ番目の回転位置ＰｏｓｋにおいてＮ＝１２で球１０６の数Ｎのとき、重心位置ＩｍＰｏｓｋ＿Ｓｐｈｒ＿（１～Ｎ）が１２個の球１０６の重心位置ＩｍＰｏｓ１＿Ｓｐｈｒ＿１～ＩｍＰｏｓ１＿Ｓｐｈｒ＿１２それぞれを示すものとする。
最初に、球１０６を既知の相対位置間隔Ｐｕ、Ｐｖで複数備えた校正治具１０２を、回転テーブル１２０に載置する（図４のステップＳ２；載置工程）。そして、回転テーブル１２０のまだ回転していない状態をｋ＝１とする（図４のステップＳ４）。
次に、Ｘ線１１８を校正治具１０２に照射する（図４のステップＳ６）、そして、Ｘ線画像検出器１２４の出力から、Ｎ個（Ｎ＝１２）の球１０６それぞれの投影像の重心位置（特徴点の位置）ＩｍＰｏｓｋ＿Ｓｐｈｒ＿（１～１２）を特定する（図４のステップＳ８；なお、ステップＳ６～ステップＳ８が特徴位置算出工程）。
次に、１２個の球１０６それぞれの投影像の重心位置ＩｍＰｏｓｋ＿Ｓｐｈｒ＿（１～１２）と球１０６の相対位置Ｘ（１～１２）から、球１０６のＸ線画像検出器１２４の検出面１２４Ａへの射影変換を行う射影変換行列Ｈｋを算出する（図４のステップＳ１０；変換行列算出工程）。
次に、回転位置Ｐｏｓｋの数ｋがＱ以上（本実施形態では３回以上であればよい）であるかを判断する（図４のステップＳ１２）。回転位置Ｐｏｓｋの数ｋがＱ（Ｑ≧３）以上でなければ（図４のステップＳ１２でＮｏ）、回転テーブル１２０を所定角度αで回転させる（図４のステップＳ１４）。そして、回転位置Ｐｏｓｋの数ｋを１つ増加させ（図４のステップＳ１６）、ステップＳ６からステップＳ１２までを繰り返す（ステップＳ６～ステップＳ１６；回転検出工程）。回転位置Ｐｏｓｋの数ｋがＱ（Ｑ≧３）以上となった際には（図４のステップＳ１２でＹｅｓ）、ステップＳ１８に進む。つまり、回転検出工程では、回転テーブル１２０を所定角度αで２回以上回転させ、特徴位置算出工程と変換行列算出工程とを繰り返す。なお、本実施形態では、所定角度αは、例えば一定の３０度としているが、特に限定されず、より小さな角度でもよいし、所定角度αが毎回変化してもよい。
次に、射影変換行列Ｈｋ（ｋ＝１～Ｑ）に基づいて回転テーブル１２０の回転中心位置Ｃｐ及び回転軸Ａｘを算出する（中心位置算出工程）。この中心位置算出工程の詳細を詳細に説明する。
まず、図５に示す如く、回転テーブル１２０の代わりにＸ線源１１６とＸ線画像検出器１２４とが回転したと想定する（図５のステップＳ３０）。ちなみに、図７（Ａ）には、回転テーブル１２０が回転した際の、所定角度αと球１０６の軌跡Ｆｂとが示されている。そして、図７（Ｂ）には、Ｘ線源１１６とＸ線画像検出器１２４とが回転したと想定した際のＸ線源１１６の絶対位置Ｘｓの軌跡Ｆｓが記載されている。
次に、射影変換行列Ｈｋ（ｋ＝１～Ｑ）に基づいて所定角度αの回転毎、即ちＱ箇所のＸ線源１１６の絶対位置Ｘｓを算出する（図５のステップＳ３２）。
なお、上記中心位置算出工程で、Ｑ箇所のＸ線源１１６の絶対位置Ｘｓが算出される際に、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃとが不明である場合を、図６を用いて、以下に説明する。
まず、ｋ番目の想定回転位置におけるＸ線源１１６の絶対位置Ｘｓを算出する際に、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とを変数とする（図６のステップＳ４０）。そして、射影変換行列Ｈｋに基づいて算出されるｋ番目の想定回転位置におけるＸ線源１１６の絶対位置Ｘｓを真円にフィッティングした仮真円の軌跡Ｆｓ上の位置とＸ線源１１６の絶対位置Ｘｓとの距離誤差を評価する（図６のステップＳ４２）。そして、この距離誤差が最小の誤差となる距離ｆ及び位置Ｃｃとを算出する（図６のステップＳ４４）。
具体的には、例えば、距離ｆを適当な値で仮決めし、位置Ｃｃを変化させ、その際に最小の距離誤差となる位置Ｃｃを算出する。次に、その最小の距離誤差となる位置Ｃｃで仮決めし、今度は距離ｆを変化させ、その際に最小の距離誤差となる距離ｆを算出する。また、その最小の距離誤差となる距離ｆに仮決めし、位置Ｃｃを変化させ、その際に最小の距離誤差となる位置Ｃｃを算出する。また、その最小の距離誤差となる位置Ｃｃで仮決めし、再度距離ｆを変化させ、その際に最小の距離誤差となる距離ｆを算出する。これを何度か繰り返すことで、最小の距離誤差となる距離ｆ及び位置Ｃｃとを算出することができ、距離ｆ及び位置Ｃｃの最適化を図ることができる。
次に、図４に戻り、Ｘ線源１１６の絶対位置Ｘｓ（１～Ｎ）の変化から、真円（＝仮真円）でフィッティングされた軌跡Ｆｓの中心位置Ｃｐを算出し、その中心位置Ｃｐを回転テーブル１２０の回転中心位置Ｃｐとする。より詳しく説明するならば、Ｑ箇所のＸ線源１１６の絶対位置Ｘｓを真円にフィッティングさせる（図４のステップＳ２０）。このとき、Ｑ＞３以上であれば、例えば、最小二乗法にて、真円の中心位置Ｃｐを算出する。Ｑ＝３であれば、例えば、連立方程式にて、真円の中心位置Ｃｐを算出する。
そして、例えば、真円でフィッティングされた軌跡Ｆｓの水平面（ｘｚ平面）からの傾斜角度を算出する。そして、回転テーブル１２０の回転中心位置Ｃｐ及び回転軸Ａｘを算出する（図４のステップＳ２２）。この時点で、ホストコンピュータ１２８により、Ｎ＝１２個の球１０６それぞれで、真円の中心位置Ｃｐとその軌跡Ｆｂを算出することもできる。このため、１２個の球１０６の真円の中心位置Ｃｐを平均化して、回転中心位置Ｃｐを算出するとともに、それらの軌跡Ｆｂの水平面からの傾斜を平均化することで、回転軸Ａｘの傾斜角度を算出でき、回転軸Ａｘを算出することもできる。
このように、本実施形態では、Ｘ線画像検出器１２４への投影像によって特定可能な形状の球１０６を既知の相対位置間隔Ｐｕ、Ｐｖで１２個備えた校正治具１０２を回転テーブル１２０に載置して、回転テーブル１２０を３箇所の回転角度に合わせて、校正治具１０２の投影像を取るだけの極めて簡単な一連の工程で、回転テーブル１２０の回転中心位置Ｃｐを算出することができる。つまり、本実施形態では、回転中心位置Ｃｐを算出するのに、三次元ボリュームデータを作成する必要がない。
また、本実施形態では、校正治具１０２において球１０６の全てが１つの平面上にのみ載置されているので、ｋ番目の回転位置Ｐｏｓｋにおける球１０６のＸ線画像検出器１２４の検出面１２４Ａへの射影変換を行う変換行列が射影変換行列Ｈｋとされていた。このため、１２個の球１０６のうち４個の球１０６だけを各工程における算出対象として、回転テーブル１２０の回転中心位置Ｃｐを算出することもでき、更なる校正時間の短縮が可能である。なお、本実施形態では、４個の球１０６だけでなく、１２個全ての球１０６を各工程における演算対象としたことで、回転テーブル１２０の回転中心位置Ｃｐを極めて正確に算出することができる。
また、本実施形態では、中心位置算出工程で、更に、回転テーブル１２０の回転軸Ａｘを算出する。このため、仮に回転テーブル１２０の回転軸Ａｘに対して当初校正不要と想定していても、実際に回転テーブル１２０の回転軸Ａｘを算出した結果と比較して、校正の必要性を正当に評価することができる。
また、本実施形態では、中心位置算出工程で、回転テーブル１２０の代わりにＸ線源１１６とＸ線画像検出器１２４とが回転したと想定し、射影変換行列Ｈｋから所定角度αの回転毎のＸ線源１１６の絶対位置Ｘｓを算出することで、回転テーブル１２０の回転中心位置Ｃｐを算出する。即ち、球１０６の絶対位置Ｘａを算出するのではなく、Ｘ線源１１６の絶対位置Ｘｓを算出している。このため、射影変換行列Ｈｋを直接的に使用することにより、結果的に、演算量を低減し、迅速に校正を実現することができる。なお、これに限定されず、球１０６の絶対位置Ｘａを算出することで、回転テーブル１２０の回転中心位置Ｃｐを算出するようにしてもよい。
また、本実施形態では、回転テーブル１２０を所定角度αで４回以上回転させ、Ｑ箇所のＸ線源１１６の絶対位置Ｘｓを算出する際に、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とを変数とする。そして、射影変換行列Ｈｋ（ｋ＝１～Ｑ）に基づいて算出されるＱ箇所のＸ線源１１６の絶対位置Ｘｓを真円にフィッティングした仮真円の軌跡Ｆｓ上の位置とＱ箇所のＸ線源１１６の絶対位置Ｘｓとの距離誤差を評価する。これにより、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とを算出する。このため、Ｘ線源１１６とＸ線画像検出器１２４との距離ｆと、Ｘ線源１１６からのＸ線画像検出器１２４への垂線の足の位置Ｃｃ（ｃｘ，ｃｙ）とを校正しようとした際には、これらの値を算出することができ、より正確な校正を行うことができる。
また、本実施形態では、中心位置算出工程で、Ｘ線源１１６の絶対位置Ｘｓの変化から、真円（＝仮真円）でフィッティングされた軌跡Ｆｓの中心位置Ｃｐを算出し、この中心位置Ｃｐを回転テーブル１２０の回転中心位置Ｃｐとする。つまり、真円でフィッティングすることで、回転位置の総数Ｑを低減でき、且つ中心位置Ｃｐを一義的に算出することができる。なお、これに限定されず、他の手法にて、回転テーブル１２０の回転中心位置Ｃｐを算出するようにしてもよい。
また、本実施形態では、更に、回転テーブル１２０の回転軸Ａｘを算出する際には、真円でフィッティングされた軌跡Ｆｓの水平面からの傾斜角度を算出し、その傾斜角度と回転中心位置Ｃｐとから回転軸Ａｘを算出する。このため、回転軸Ａｘを算出するのに、球１０６は１個あればよいので、回転軸Ａｘを算出する工程を簡略化でき且つ短時間で行うことができる。なお、これに限定されず、例えば、各球１０６において真円でフィッティングされた軌跡Ｆｓを算出して、その中心位置のずれから、回転軸Ａｘを算出するようにしてもよい。
また、本実施形態では、校正治具１０２上の基準物体は、球１０６とされている。このため、球１０６は、いずれの方向から投影されても輪郭が円となる。即ち、球１０６は、基準物体として、Ｘ線画像検出器１２４への投影像によって最も容易に特定可能な形状である。なお、これに限定されず、基準物体が、例えば、正多面体や変形した菱形状体を含む多面体であってもよいし、楕円体や円錐体などの曲面を含む形状とされていてもよい。
また、本実施形態では、基準物体である球１０６の投影像の特徴点の位置が、投影像の重心位置とされている。球１０６の投影像は円であることから、重心位置を算出することが容易であり、少ない位置誤差で算出することができる。なお、これに限らず、基準物体である球１０６の投影像の特徴点の位置が中心位置であってもよい。そもそも基準物体が球ではなくてもよい。例えば、校正治具１０２に設けられた局所的に特徴的な凹部や凸部を基準物体として、その特徴的な凹部や凸部を投影像の特徴点に関連付けるようにしてもよい。
即ち、本実施形態では、被測定物を回転可能に載置する回転テーブル１２０の回転中心位置Ｃｐを単純な工程で容易に算出することが可能である。
なお、上記実施形態では、校正治具１０２において球１０６の全てが１つの平面上にのみ載置されていたが、本発明はこれに限定されない。例えば、校正治具１０２において球１０６の全てが１つの平面上に載置された状態とならず、三次元的に載置されていてもよい。その際には、球１０６が少なくとも６個以上とされ、射影変換行列Ｈｋの代わりに射影行列Ｐｋが用いられる。ホストコンピュータ１２８は、式（３）の代わりに、以下の射影行列Ｐｋに係る式（５）を用いて、ｋ番目の想定回転位置の３行＊４列の射影行列Ｐｋから、Ｘ線源１１６の絶対位置Ｘｓを算出することができる。
Ｐｋ＝Ａ［ｒｋ１ ｒｋ２ ｒｋ３ Ｔｋ］ （５）
この場合には、射影行列Ｐｋを用いることで、仮に校正治具１０２の平面精度がよくなくても、正確な校正を行うことができる。
なお、上記実施形態では、球１０６が少なくとも４個（あるいは６個）とされていたが、本発明はこれに限定されない。例えば、１個の球１０６が少なくとも４箇所（あるいは６箇所）に移動して配置されるような校正治具１０２の構成であってもよい。
なお、上記実施形態では、図３（Ａ）、（Ｂ）に示す如く、ｕ、ｖ軸方向にそれぞれ、一定の相対位置間隔Ｐｕ、Ｐｖで球１０６が配置されていたが、本発明はこれに限定されない。例えば、相対位置間隔Ｐｕ、Ｐｖは各球１０６間で異なっていてもよく、ｕ、ｖ軸方向で並んでいる必要はない。これは、平面上に二次元的に球１０６が配置された場合だけでなく、三次元的に球１０６が配置された場合でも同様である。
本発明は、Ｘ線計測装置の校正に広く適用することができる。
１００…Ｘ線計測装置
１０２…校正治具
１０４…板状部材
１０６…球
１０８…本体部
１１０…Ｘ線遮蔽カバー
１１２…ベース
１１４…線源支持台
１１６…Ｘ線源
１１８…Ｘ線
１２０…回転テーブル
１２２…テーブル支持台
１２４…Ｘ線画像検出器
１２４Ａ…検出面
１２６…検出器支持台
１２８…ホストコンピュータ
１３０…モーションコントローラ
Ａ…内部パラメータ行列
Ａｘ…回転軸
Ｃｃ…Ｘ線源からのＸ線画像検出器への垂線の足の位置
Ｃｐ…回転中心位置
Ｄ…球の直径
Ｆｂ、Ｆｓ…軌跡
ＦＣＤ…Ｘ線源と回転中心位置との距離
ｆ…Ｘ線源とＸ線画像検出器との距離
Ｈｋ…射影変換行列
ＩｍＰｏｓｋ＿Ｓｐｈｒ＿（１～Ｎ）…ｋ番目の回転位置におけるＮ個の球それぞれの重心位置
ｋ、Ｎ、Ｑ…数
Ｐｋ…射影行列
Ｐｏｓｋ…ｋ番目の回転位置
Ｐｕ、Ｐｖ…球の既知の相対位置間隔
Ｒｋ…ｋ番目の回転行列
Ｘ…球の相対位置
Ｘａ…球の絶対位置
Ｘｓ…Ｘ線源の絶対位置
α…所定角度

Claims (9)

  1. Ｘ線を用いて被測定物を三次元形状計測するＸ線計測装置の校正方法であって、
    前記Ｘ線計測装置は、前記Ｘ線を発生させるＸ線源と、前記被測定物を回転可能に載置する回転テーブルと、該被測定物を透過した該Ｘ線を検出するＸ線画像検出器と、を備え、
    該Ｘ線画像検出器への投影像によって特定可能な形状の基準物体を既知の相対位置間隔でＮ箇所（Ｎ≧４）以上に配置可能な校正治具を、前記回転テーブルに載置する載置工程と、
    前記Ｘ線を該校正治具に照射し該Ｘ線画像検出器の出力から、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置を特定する特徴位置算出工程と、
    Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置と前記既知の相対位置間隔から、該基準物体の前記Ｘ線画像検出器の検出面への射影変換を行う変換行列を算出する変換行列算出工程と、
    前記回転テーブルを所定角度で２回以上回転させ、前記特徴位置算出工程と前記変換行列算出工程とを繰り返す回転検出工程と、
    前記変換行列に基づいて前記回転テーブルの回転中心位置を算出する中心位置算出工程と、
    を含むことを特徴とするＸ線計測装置の校正方法。
  2. 請求項１において、
    前記校正治具において前記基準物体の全てが１つの平面上にのみ載置されている場合には、前記変換行列は射影変換行列とされ、該基準物体が三次元的に載置されている場合には、Ｎ＝６とされ且つ該変換行列は射影行列とされることを特徴とするＸ線計測装置の校正方法。
  3. 請求項１または２において、
    前記中心位置算出工程で、更に、前記回転テーブルの回転軸を算出することを特徴とするＸ線計測装置の校正方法。
  4. 請求項１乃至３のいずれかにおいて、
    前記中心位置算出工程で、前記回転テーブルの代わりに前記Ｘ線源と前記Ｘ線画像検出器とが回転したと想定し、前記変換行列に基づいて前記所定角度の回転毎の該Ｘ線源の絶対位置を算出することで、前記回転テーブルの回転中心位置を算出することを特徴とするＸ線計測装置の校正方法。
  5. 請求項４において、
    前記回転テーブルを前記所定角度で４回以上回転させ、前記Ｘ線源の絶対位置を算出する際に、該Ｘ線源と前記Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを変数とし、前記変換行列に基づいて算出される該Ｘ線源の絶対位置を真円にフィッティングした仮真円の軌跡上の位置と該Ｘ線源の絶対位置との距離誤差を評価することで、該Ｘ線源と該Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを算出することを特徴とするＸ線計測装置の校正方法。
  6. 請求項４又は５において、
    前記中心位置算出工程で、前記Ｘ線源の絶対位置の変化から、真円でフィッティングされた軌跡の中心位置を算出し、該中心位置を前記回転テーブルの回転中心位置とすることを特徴とするＸ線計測装置の校正方法。
  7. 請求項６において、
    更に、前記回転テーブルの回転軸を算出する際には、前記軌跡の水平面からの傾斜角度を算出し、該傾斜角度と前記回転中心位置とから該回転軸を算出することを特徴とするＸ線計測装置の校正方法。
  8. 請求項１乃至７のいずれかにおいて、
    前記基準物体は、球とされていることを特徴とするＸ線計測装置の校正方法。
  9. 請求項１乃至８のいずれかにおいて、
    前記基準物体の投影像の特徴点の位置は、該投影像の重心位置とされていることを特徴とするＸ線計測装置の校正方法。
JP2020111821A 2020-06-29 2020-06-29 Ｘ線計測装置の校正方法 Pending JP2022010983A (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2020111821A JP2022010983A (ja) 2020-06-29 2020-06-29 Ｘ線計測装置の校正方法
US17/360,304 US11573190B2 (en) 2020-06-29 2021-06-28 Calibration method for X-ray measuring device
DE102021003348.5A DE102021003348A1 (de) 2020-06-29 2021-06-29 Kalibrierverfahren für eine röntgenstrahl-messvorrichtung
CN202110724806.5A CN113932741A (zh) 2020-06-29 2021-06-29 X射线测量装置的校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020111821A JP2022010983A (ja) 2020-06-29 2020-06-29 Ｘ線計測装置の校正方法

Publications (1)

Publication Number Publication Date
JP2022010983A true JP2022010983A (ja) 2022-01-17

Family

ID=78827261

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020111821A Pending JP2022010983A (ja) 2020-06-29 2020-06-29 Ｘ線計測装置の校正方法

Country Status (4)

Country Link
US (1) US11573190B2 (ja)
JP (1) JP2022010983A (ja)
CN (1) CN113932741A (ja)
DE (1) DE102021003348A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114460110B (zh) * 2022-03-08 2023-06-06 中国电子科技集团公司第三十八研究所 一种伺服系统误差补偿方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2779556B1 (fr) * 1998-06-05 2000-08-11 Ge Medical Syst Sa Dispositif et procede de simulation des vaisseaux sanguins d'un patient
JP3792432B2 (ja) 1999-04-14 2006-07-05 東芝Ｉｔコントロールシステム株式会社 コンピュータ断層撮影装置
TWI398160B (zh) * 2009-12-01 2013-06-01 Ind Tech Res Inst 攝影機校正系統與座標資料產生系統及其方法
US8764290B2 (en) * 2012-01-30 2014-07-01 Hexagon Metrology, Inc. X-ray computed tomography device calibration and verification apparatus
EP3421086B1 (en) * 2017-06-28 2020-01-15 OptiNav Sp. z o.o. Determination of geometrical information about a medical treatment arrangement comprising a rotatable treatment radiation source unit
JP2019158534A (ja) * 2018-03-12 2019-09-19 株式会社ミツトヨ 計測用ｘ線ｃｔ装置、及び、その断層画像生成方法
KR20200087549A (ko) 2019-01-11 2020-07-21 캐논 톡키 가부시키가이샤 성막장치, 전자 디바이스 제조장치, 성막방법, 및 전자 디바이스 제조방법
JP7257925B2 (ja) * 2019-09-10 2023-04-14 株式会社ミツトヨ Ｘ線計測装置の校正方法

Also Published As

Publication number Publication date
US20210404976A1 (en) 2021-12-30
US11573190B2 (en) 2023-02-07
DE102021003348A1 (de) 2021-12-30
CN113932741A (zh) 2022-01-14

Similar Documents

Publication Publication Date Title
US9625257B2 (en) Coordinate measuring apparatus and method for measuring an object
JP6735667B2 (ja) コンピュータ断層撮影の較正装置および方法
JP6767045B2 (ja) 計測用ｘ線ｃｔ装置と座標測定機の座標合せ治具
JP4415762B2 (ja) 断層撮影装置
US20040234025A1 (en) Processes and a device for determining the actual position of a structure of an object to be examined
JP5615260B2 (ja) 機械的ワークを断層撮影法によって測定するための方法
CN107684435A (zh) 锥束ct系统几何校准方法及其校准装置
JP2019158534A (ja) 計測用ｘ線ｃｔ装置、及び、その断層画像生成方法
GB2533659A (en) Method and device for correcting computed tomography measurements, comprising a coordinate measuring machine
US9606072B2 (en) Radiation inspecting apparatus
JP7257925B2 (ja) Ｘ線計測装置の校正方法
JP7257924B2 (ja) Ｘ線計測装置の校正方法
CN110686624B (zh) 利用由x射线ct装置获得的投影像的尺寸测定方法
CN114596222A (zh) 适用于一般轨迹锥束ct系统几何校正的模体与标定方法
JP2022010983A (ja) Ｘ線計測装置の校正方法
JP2023088704A (ja) 放射線撮影装置、画像生成方法、放射線撮影システム、及びプログラム
CN209032406U (zh) 一种锥束ct系统几何校准装置
Blumensath et al. Calibration of robotic manipulator systems for cone-beam tomography imaging
JP2019158541A (ja) 計測用ｘ線ｃｔ装置、及び、その量産ワーク測定方法
US7430271B2 (en) Ray tracing kernel
JP2023095304A (ja) Ｘ線計測装置の校正方法
CN110501360B (zh) 一种用于显微ct系统位姿校正的标准器及实现方法
JP6015404B2 (ja) 放射線検査装置
JP2022150461A (ja) 公差誤差推定装置、方法、プログラム、再構成装置および制御装置
CN112748135B (zh) 一种外部ct扫描模式下转台旋转中心标定方法及装置