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

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

Info

Publication number
JP7257924B2
JP7257924B2 JP2019163979A JP2019163979A JP7257924B2 JP 7257924 B2 JP7257924 B2 JP 7257924B2 JP 2019163979 A JP2019163979 A JP 2019163979A JP 2019163979 A JP2019163979 A JP 2019163979A JP 7257924 B2 JP7257924 B2 JP 7257924B2
Authority
JP
Japan
Prior art keywords
ray
calculating
calculated
locations
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.)
Active
Application number
JP2019163979A
Other languages
English (en)
Other versions
JP2021043014A (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.)
Mitutoyo Corp
Original Assignee
Mitutoyo Corp
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 filed Critical Mitutoyo Corp
Priority to JP2019163979A priority Critical patent/JP7257924B2/ja
Priority to US17/014,308 priority patent/US11346660B2/en
Priority to DE102020005518.4A priority patent/DE102020005518A1/de
Priority to CN202010940511.7A priority patent/CN112461165A/zh
Publication of JP2021043014A publication Critical patent/JP2021043014A/ja
Application granted granted Critical
Publication of JP7257924B2 publication Critical patent/JP7257924B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G01B15/045—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 by measuring absorption
    • 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
    • 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
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/03—Investigating materials by wave or particle radiation by transmission
    • 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/60—Specific applications or type of materials
    • G01N2223/646—Specific applications or type of materials flaws, defects

Landscapes

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

Description

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

Claims (13)

  1. Ｘ線を用いて被測定物を三次元形状計測するＸ線計測装置の校正方法であって、
    前記Ｘ線計測装置は、前記Ｘ線を発生させるＸ線源と、前記被測定物を回転可能に載置する回転テーブルと、該被測定物を透過した該Ｘ線を検出するＸ線画像検出器と、を備え、
    該Ｘ線画像検出器への投影像によって特定可能な形状の基準物体を特定の相対位置間隔でＮ箇所（Ｎ≧４）に配置可能な校正治具を、前記回転テーブルに載置する載置工程と、
    前記基準物体の特定の相対位置間隔を変化させずにＮ箇所に配置された該基準物体の平行移動を複数回行い、該平行移動の前後に前記Ｘ線を前記校正治具に照射し前記Ｘ線画像検出器の出力から、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置を特定する前段特徴位置算出工程と、
    前記平行移動の前後の位置と対応するＮ箇所にある該基準物体それぞれの投影像の特徴点の位置とから、該基準物体の前記Ｘ線画像検出器の検出面への射影変換を行う個別変換行列をＮ箇所にある該基準物体それぞれに対して算出する個別行列算出工程と、
    前記個別変換行列に基づきＮ箇所にある前記基準物体それぞれの移動位置を算出する個別位置算出工程と、
    前記平行移動の前後の位置それぞれに、Ｎ箇所にある該基準物体それぞれの移動位置を付加し、Ｎ箇所にある該基準物体の相対位置を算出する座標同一工程と、
    前記Ｘ線を前記校正治具に照射し前記Ｘ線画像検出器の出力から、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置を特定する後段特徴位置算出工程と、
    Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置と前記相対位置から、該基準物体の前記Ｘ線画像検出器の検出面への射影変換を行う第１変換行列を算出する変換行列算出工程と、
    前記回転テーブルを所定角度で２回以上回転させ、前記後段特徴位置算出工程から前記変換行列算出工程までを繰り返す回転検出工程と、
    前記第１変換行列に基づき該所定角度の回転毎の前記基準物体の絶対位置を算出する位置算出工程と、
    前記回転テーブルの回転で生じる該基準物体の絶対位置の変化から、前記回転テーブルの回転中心位置を算出する中心位置算出工程と、
    を含むことを特徴とするＸ線計測装置の校正方法。
  2. 請求項１において、
    前記個別行列算出工程で、Ｎ箇所に配置された前記基準物体の代わりに前記Ｘ線源と前記Ｘ線画像検出器とが平行移動したと想定し、前記個別変換行列に基づきＮ箇所毎に該Ｘ線源の移動位置を算出して、該Ｘ線源の移動位置に基づきＮ箇所にある前記基準物体それぞれの移動位置を算出することを特徴とするＸ線計測装置の校正方法。
  3. 請求項１または２において、
    前記座標同一工程で、前記平行移動の前後の位置それぞれに、Ｎ箇所にある前記基準物体それぞれの移動位置を付加した結果を補正後移動位置とし、Ｎ箇所にある該基準物体それぞれの投影像の特徴点の位置と該補正後移動位置から、該基準物体の前記Ｘ線画像検出器の検出面への射影変換を行う第２変換行列を算出することで、Ｎ箇所にある該基準物体の相対位置を算出することを特徴とするＸ線計測装置の校正方法。
  4. 請求項３において、
    前記座標同一工程で、算出された前記補正後移動位置及び前記第２変換行列を、変数とした該補正後移動位置及び該第２変換行列の初期値とし、該補正後移動位置と該第２変換行列との関係に基づいて算出されたＮ箇所にある前記基準物体それぞれの投影像の特徴点の位置と、実際に検出されたＮ箇所にある前記基準物体それぞれの投影像の特徴点の位置との位置誤差を評価することで、前記補正後移動位置と前記第２変換行列とを算出することを特徴とするＸ線計測装置の校正方法。
  5. 請求項１乃至４のいずれかにおいて、
    前記回転テーブルを特定の角度で複数回回転させ、前記前段特徴位置算出工程から前記座標同一工程までを繰り返し、
    繰り返しで得られる複数回の前記相対位置の平均を算出する、あるいは次に繰り返される際の前記平行移動の前後の位置を直前で算出された該相対位置と対応させることを特徴とするＸ線計測装置の校正方法。
  6. 請求項１乃至５のいずれかにおいて、
    前記校正治具において前記基準物体の全てが１つの平面上にのみ載置されている場合には、前記第１変換行列は射影変換行列とされ、該基準物体が三次元的に載置されている場合には、該第１変換行列は射影行列とされることを特徴とするＸ線計測装置の校正方法。
  7. 請求項１乃至６のいずれかにおいて、
    前記中心位置算出工程で、更に、前記回転テーブルの回転軸を算出することを特徴とするＸ線計測装置の校正方法。
  8. 請求項１乃至７のいずれかにおいて、
    前記位置算出工程で、前記回転テーブルの代わりに前記Ｘ線源と前記Ｘ線画像検出器とが回転したと想定し、前記第１変換行列に基づき前記所定角度の回転毎の該Ｘ線源の絶対位置を算出して、該Ｘ線源の絶対位置を座標変換することで、前記基準物体の絶対位置を算出することを特徴とするＸ線計測装置の校正方法。
  9. 請求項８において、
    前記回転テーブルを前記所定角度で３回以上回転させ、前記Ｘ線源の絶対位置を算出する際に、該Ｘ線源と前記Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを変数化し、前記第１変換行列に基づき算出される該Ｘ線源の絶対位置を真円にフィッティングした仮真円の軌跡上の位置と該Ｘ線源の絶対位置との距離誤差を評価することで、該Ｘ線源と該Ｘ線画像検出器との距離と、該Ｘ線源からの該Ｘ線画像検出器への垂線の足の位置とを算出することを特徴とするＸ線計測装置の校正方法。
  10. 請求項１乃至９のいずれかにおいて、
    前記中心位置算出工程で、前記基準物体の絶対位置の変化から、真円でフィッティングされた軌跡の中心位置を算出し、該中心位置を前記回転テーブルの回転中心位置とすることを特徴とするＸ線計測装置の校正方法。
  11. 請求項１０において、
    更に、前記回転テーブルの回転軸を算出する際には、前記軌跡の水平面からの傾斜角度を算出し、該傾斜角度と前記回転中心位置とから該回転軸を算出することを特徴とするＸ線計測装置の校正方法。
  12. 請求項１乃至１１のいずれかにおいて、
    前記基準物体は、球とされていることを特徴とするＸ線計測装置の校正方法。
  13. 請求項１乃至１２のいずれかにおいて、
    前記基準物体の投影像の特徴点の位置は、該投影像の重心位置とされていることを特徴とするＸ線計測装置の校正方法。
JP2019163979A 2019-09-09 2019-09-09 Ｘ線計測装置の校正方法 Active JP7257924B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2019163979A JP7257924B2 (ja) 2019-09-09 2019-09-09 Ｘ線計測装置の校正方法
US17/014,308 US11346660B2 (en) 2019-09-09 2020-09-08 Calibration method of x-ray measuring device
DE102020005518.4A DE102020005518A1 (de) 2019-09-09 2020-09-09 Kalibrierverfahren einer röntgen-messvorrichtung
CN202010940511.7A CN112461165A (zh) 2019-09-09 2020-09-09 X射线测量装置的校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019163979A JP7257924B2 (ja) 2019-09-09 2019-09-09 Ｘ線計測装置の校正方法

Publications (2)

Publication Number Publication Date
JP2021043014A JP2021043014A (ja) 2021-03-18
JP7257924B2 true JP7257924B2 (ja) 2023-04-14

Family

ID=74644307

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019163979A Active JP7257924B2 (ja) 2019-09-09 2019-09-09 Ｘ線計測装置の校正方法

Country Status (4)

Country Link
US (1) US11346660B2 (ja)
JP (1) JP7257924B2 (ja)
CN (1) CN112461165A (ja)
DE (1) DE102020005518A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7257925B2 (ja) * 2019-09-10 2023-04-14 株式会社ミツトヨ Ｘ線計測装置の校正方法
CN113791095B (zh) * 2021-11-18 2022-02-11 煤炭科学研究总院 一种ct扫描的试样精准位置调节方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005127881A (ja) 2003-10-24 2005-05-19 Shimadzu Corp Ｘ線ｃｔ装置における回転中心軸の較正方法およびｘ線ｃｔ装置
JP2006003200A (ja) 2004-06-17 2006-01-05 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2007078469A (ja) 2005-09-13 2007-03-29 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2007333596A (ja) 2006-06-15 2007-12-27 Shimadzu Corp Ｘ線ｃｔ装置
JP2012233919A (ja) 2004-05-26 2012-11-29 Werth Messtechnik Gmbh 座標測定装置及び測定物の測定方法
JP2015525665A (ja) 2012-08-20 2015-09-07 オレンジデンタル・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング・ウント・コムパニー・コマンディットゲゼルシャフトｏｒａｎｇｅｄｅｎｔａｌ ＧｍｂＨ ＆ Ｃｏ． ＫＧ コーンビームコンピュータートモグラフィー装置の幾何学的特徴付け及び較正
JP2016538552A (ja) 2013-11-28 2016-12-08 ニコン・メトロロジー・エヌヴェ コンピュータ断層撮影の較正装置および方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3792432B2 (ja) 1999-04-14 2006-07-05 東芝Ｉｔコントロールシステム株式会社 コンピュータ断層撮影装置
JP5572521B2 (ja) * 2010-10-26 2014-08-13 株式会社日立メディコ Ｘ線ｃｔ装置、ｘ線ｃｔ装置の画像再構成方法
WO2017052431A1 (en) * 2015-09-23 2017-03-30 Prismatic Sensors Ab Determining the orientation of an edge-on x-ray detector with respect to the direction of incoming x-rays
US10835199B2 (en) * 2016-02-01 2020-11-17 The University Of North Carolina At Chapel Hill Optical geometry calibration devices, systems, and related methods for three dimensional x-ray imaging
JP6693533B2 (ja) * 2016-02-16 2020-05-13 株式会社ニコン Ｘ線装置、ｘ線計測方法および構造物の製造方法
CN106526686B (zh) * 2016-12-07 2019-05-07 同方威视技术股份有限公司 螺旋ct设备和三维图像重建方法
US11402524B2 (en) * 2017-11-14 2022-08-02 Imatrex, Inc. Geometric calibration of X-ray imaging systems
JP2019128163A (ja) * 2018-01-19 2019-08-01 株式会社ミツトヨ 計測用ｘ線ｃｔ装置、及び、その校正方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005127881A (ja) 2003-10-24 2005-05-19 Shimadzu Corp Ｘ線ｃｔ装置における回転中心軸の較正方法およびｘ線ｃｔ装置
JP2012233919A (ja) 2004-05-26 2012-11-29 Werth Messtechnik Gmbh 座標測定装置及び測定物の測定方法
JP2006003200A (ja) 2004-06-17 2006-01-05 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2007078469A (ja) 2005-09-13 2007-03-29 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2007333596A (ja) 2006-06-15 2007-12-27 Shimadzu Corp Ｘ線ｃｔ装置
JP2015525665A (ja) 2012-08-20 2015-09-07 オレンジデンタル・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング・ウント・コムパニー・コマンディットゲゼルシャフトｏｒａｎｇｅｄｅｎｔａｌ ＧｍｂＨ ＆ Ｃｏ． ＫＧ コーンビームコンピュータートモグラフィー装置の幾何学的特徴付け及び較正
JP2016538552A (ja) 2013-11-28 2016-12-08 ニコン・メトロロジー・エヌヴェ コンピュータ断層撮影の較正装置および方法

Also Published As

Publication number Publication date
CN112461165A (zh) 2021-03-09
JP2021043014A (ja) 2021-03-18
US11346660B2 (en) 2022-05-31
DE102020005518A1 (de) 2021-03-11
US20210072022A1 (en) 2021-03-11

Similar Documents

Publication Publication Date Title
US9625257B2 (en) Coordinate measuring apparatus and method for measuring an object
US20150369757A1 (en) Method and apparatus for generating a three-dimensional model of a region of interest using an imaging system
JP7257924B2 (ja) Ｘ線計測装置の校正方法
CN103784160B (zh) 一种锥束ct***几何位置的校正装置及其校正方法
EP3306309B1 (en) Image acquisition device, image acquisition method, and image correction program
US9606072B2 (en) Radiation inspecting apparatus
JP7257925B2 (ja) Ｘ線計測装置の校正方法
JP2019158534A (ja) 計測用ｘ線ｃｔ装置、及び、その断層画像生成方法
JP6711410B2 (ja) 放射線断層撮影装置の撮像倍率校正方法
TW561241B (en) Method and apparatus for calibrating laser three-dimensional digitizing sensor
CN110686624B (zh) 利用由x射线ct装置获得的投影像的尺寸测定方法
CN114596222A (zh) 适用于一般轨迹锥束ct***几何校正的模体与标定方法
CN209032406U (zh) 一种锥束ct***几何校准装置
CN112603346B (zh) 一种基于标记物成像的探测器偏转校正方法
JP2022010983A (ja) Ｘ線計測装置の校正方法
CN113812971A (zh) 一种多自由度四维双能锥束ct成像***及方法
Blumensath et al. Calibration of robotic manipulator systems for cone-beam tomography imaging
JP2005233760A (ja) トモシンセシス装置
JP2023095304A (ja) Ｘ線計測装置の校正方法
JP2019158541A (ja) 計測用ｘ線ｃｔ装置、及び、その量産ワーク測定方法
JP6015404B2 (ja) 放射線検査装置
CN114486955A (zh) 确定至少一个评估测量数据所需的几何参数的计算机实现方法
Meli et al. A cone-beam CT geometry correction method based on intentional misalignments to render the projection images correctable
Wicklein et al. An object-independent measure for improving misalignment correction in c-arm ct
JP2004301860A (ja) コンピュータ断層撮影装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220818

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230322

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230328

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230404

R150 Certificate of patent or registration of utility model

Ref document number: 7257924

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150