(57)【要約】  【課題】２次元放射線画像を検出するＸ線撮像システムを評価するために撮像された画像の縮率を自動かつ高い精度で測定可能な画像の縮率計算方法を提供する。  【解決手段】長さが既知の矩形片を撮像して、その画像を得る工程と、矩形片の画像のエッジを表わす４本の直線の式を求め、その式からエッジ間の距離を求め、エッジ間の距離と矩形片の長さとの比較によって画像の縮率を算出することにより、上記課題を解決する。【選択図】図１
2005058315.tif 000002
【特許請求の範囲】
【請求項１】
２次元放射線画像を撮像するシステムを評価するために検出された画像の縮率を計算するに際し、
長さが既知の矩形片を撮像して、その画像を得る工程と、
前記矩形片の画像のエッジを表わす４本の直線の式を求め、その式からエッジ間の距離を求める工程と、
前記エッジ間の距離と前記矩形片の長さとの比較によって前記画像の縮率を算出する工程と、
を含むことを特徴とする画像の縮率計算方法。
【請求項２】
前記矩形片内の領域の画像データの平均値と前記矩形片外の領域の画像データの平均値とを閾値として画像を２値化処理した後、主副方向にずらした画像を元の画像から引き、その結果が０以外のときは１とすることでエッジのみを抽出することを特徴とする請求項１に記載の画像の縮率計算方法。
【請求項３】
２次元放射線画像を撮像するシステムを評価するために検出された画像の縮率を計算するに際し、
長さが既知の矩形片を撮像して、その画像を得る工程と、
前記矩形片の画像内の縮率方向と直交する方向の第１のエッジ点を求める工程と、
前記第１のエッジ点間の中点を求める工程と、
前記矩形片の画像において前記中点から縮率方向の第２のエッジ点を求める工程と、
前記矩形片の画像の画像データ軸に対する傾きを求める工程と、
前記中点から前記第２のエッジ点までの距離および前記矩形片の画像の傾きから、前記矩形片の画像の軸方向の長さを求める工程と、
前記矩形片の画像の長さを前記矩形片の実際の長さと比較し、画像の縮率を算出する工程と、
を含むことを特徴とする画像の縮率計算方法。
【請求項４】
前記第１のエッジ点および前記第２のエッジ点を求める工程において、前記矩形片内の領域の画像データの平均値と前記矩形片外の領域の画像データの平均値との半値を求め、この半値を横切る前後のデータ値を直線近似して、前記第１のエッジ点および前記第２のエッジ点の位置を決定することを特徴とする請求項３に記載の画像の縮率計算方法。
【請求項５】
前記矩形片端と直交方向に、メディアンフィルター処理または／および移動平均処理を行なうことを特徴とする請求項３または４に記載の画像の縮率計算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像の縮率計算方法に関し、詳しくは、２次元Ｘ線（放射線）画像を検出するＸ線撮像システムを評価するためにＸ線撮像装置等によって撮影された２次元Ｘ線画像の縮率を計算する画像の縮率計算方法に関するものである。
【背景技術】
【０００２】
例えば、医療診断や製造業分野の非破壊試験を目的とするＸ線（放射線）画像撮影において、Ｘ線フィルムやＸ線量を記録する蓄積性蛍光体シート（イメージングプレート：ＩＰ）およびＸ線を検出して電気信号に変換するフラットパネルデバイス（ＦＰＤ）等を使用したＸ線撮像システムがある。
【０００３】
このようなＸ線撮像システムの性能を評価するために、大きさおよび材質が既知で、感度、分解能、表示精度およびダイナミックレンジ等の確認を行なうことを目的とするＱＣファントムが一般に用いられている。
一方、被写体の大きさに対するその画像の縮率は、副走査速度／光（主）走査速度と電気信号の検出周波数との関係、装置天板と蓄積性蛍光体シートとの距離による幾何学的な拡大等によって変化するため、Ｘ線撮像システムによって得られる２次元Ｘ線画像の信頼性を保証するためには、Ｘ線撮像システムによって得られたＸ線画像の縮率を測定することが必要となる。
【０００４】
このため、Ｘ線画像の縮率を測定サンプルとして、上述のＱＣファントムのファントム画像に加え、このファントム画像とは別にＸ線不透過目盛りを持つ定規を撮像し、その定規の目盛りを計算して、ファントム画像の指定された直線の幾何学的歪、例えば読取光学系の故障や異常などのトラブルにより走査光が波状に走査されることによって生じる幾何学的歪や、スキャン幅を計測し、画像の縮率を測定するものが以下に示す特許文献１に開示されている。
【特許文献１】特開平８−３０７６０５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、上記の従来の画像の縮率の計算方法では、オペレータが目視で画像の測定したい直線を指定し、この直線の長さを定規を用いて目視で計測し縮率を求めているため、画素単位でしか物体の直線位置を指定することができず、寸法誤差が大きいという事情がある。
特に、Ｘ線による画像は、Ｘ線の散乱等により輪郭がぼやけてしまうため、測定したい直線の境界部分を明確に決定することができず、この境界部分を拡大して決定しようとすると、個人差が出て、計測にばらつきが生じるという事情がある。
【０００６】
また、上記公報に記載の方法では、ＱＣファントムのファントム画像のとは別にＸ線不透過目盛りを持つ定規を撮像する必要があり、撮像された定規の目盛りを計算する必要があるという問題がある。
【０００７】
本発明の技術的課題は、前記のような事情に鑑みてなされたものであり、その目的は、２次元放射線画像を検出するＸ線撮像システムを評価するために撮像された画像の縮率を自動かつ高い精度で測定可能な画像の縮率計算方法を提供することにある。
なお、本発明の前記ならびにその他の目的と新規な特徴は、本明細書の記述および添付図面から明らかになるであろう。
【課題を解決するための手段】
【０００８】
前記課題を達成するため、本発明は、２次元放射線画像を撮像するシステムを評価するために検出された画像の縮率を計算するに際し、長さが既知の矩形片を撮像して、その画像を得る工程と、前記矩形片の画像のエッジを表わす４本の直線の式を求め、その式からエッジ間の距離を求める工程と、前記エッジ間の距離と前記矩形片の長さとの比較によって前記画像の縮率を算出する工程とを含むことを特徴とする画像の縮率計算方法を提供するものである。
ここで、前記矩形片内の領域の画像データの平均値と、前記矩形片外の領域の画像データの平均値とを閾値として画像を２値化処理した後、主副方向にずらした画像を元の画像から引き、その結果が０以外のときは１とすることでエッジのみを抽出することが好ましい。
【０００９】
また、本発明は、２次元放射線画像を撮像するシステムを評価するために検出された画像の縮率を計算するに際し、長さが既知の矩形片を撮像して、その画像を得る工程と、前記矩形片の画像内の縮率方向と直交する方向の第１のエッジ点を求める工程と、前記第１のエッジ点間の中点を求める工程と、前記矩形片の画像において、前記中点から縮率方向の第２のエッジ点を求める工程と、前記矩形片の画像の画像データ軸に対する傾きを求める工程と、前記中点から前記第２のエッジ点までの距離および前記矩形片の画像の傾きから、前記矩形片の画像の軸方向の長さを求める工程と、前記矩形片の画像の長さを前記矩形片の実際の長さと比較し、前記画像の縮率を算出する工程とを含むことを特徴とする画像の縮率計算方法を提供するものである。
【００１０】
ここで、前記第１のエッジ点および前記第２のエッジ点を求める工程において、前記矩形片内の領域の画像データの平均値と前記矩形片外の領域の画像データの平均値との半値を求め、この半値を横切る前後のデータ値を直線近似して、前記第１のエッジ点および前記第２のエッジ点の位置を決定することが好ましい。
また、前記矩形片端と直交方向に、メディアンフィルター処理または／および移動平均処理を行なうことが好ましい。
【発明の効果】
【００１１】
本発明によれば、試験片のエッジを表す直線の式から試験片の長さを求め、画像の縮率が算出されるので、画像の縮率が自動的に求められる。その際、試験片内領域の画像データの平均値と試験片外領域の画像データの平均値との半値を横切る前後の値を直線近似してエッジ位置が決定され、エッジ位置の推定誤差が小さくなるため、縮率の寸法誤差が小さく抑えられ、個人差によるばらつきもなく、縮率を高精度で正確に測定することができる。さらに、メディアンフィルター処理および移動平均処理することで、ノイズが除去され、計算精度が向上する。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施の形態を図面に基づいて詳細に説明する。
図１（ａ）および（ｂ）は、それぞれ本発明に係る画像の縮率計算方法を示す一実施例に係るフローチャートを例示したものであり、図２および図３は、それぞれ図１（ａ）および（ｂ）の工程説明図を示す。
【００１３】
図１（ａ）に示す本発明の第１実施形態の画像の縮率計算方法では、まず、ステップＳ１において、長さが分かっている長方形の金属板を試験片として撮影し、図２（ａ）に示すように、金属板の画像Ｇ1 を含む処理領域の画像Ａ1 （のＱＬ（クォンタムレベル））を得る。
なお、本発明法を実施する際には、この長方形の金属板部分（領域）を含み、この金属板部分より少し広い領域を、本発明法を適用する対象となる処理領域とするのが良い。
また、図２（ａ）中、Ｇ1 は、Ｘ線が透過しない金属板の画像であるので低濃度（透過または白）であり、Ｂ1 は、Ｘ線に曝露されて露光された部分（バックグラウンド）の画像であるので高濃度（黒）であり、Ａ1 は、金属板の画像Ｇ1 およびバックグラウンドの画像Ｂ1 の両方を含む本発明法の適用対象の処理領域全体の画像Ａ1 である。
【００１４】
ここで、例えば、ステップＳ１で用いる長さが分かっている長方形の金属板を試験片として、ＣＲ（コンピューテッドラジオグラフィー）において、感度、分解能、表示精度およびダイナミックレンジ等の機能の評価を行なうためのＱＣファントム内の長方形の金属板を用いることができる。
【００１５】
この場合には、ステップＳ１において、Ｘ線を照射してＱＣファントムを蓄積性蛍光体からなる蓄積性蛍光体シート（イメージングプレート）に撮影して、ＱＣファントムが撮影されたイメージングプレートをポリゴン走査機構を用いてレーザ光等の励起光によって光走査し、光走査によって発生する輝尽発光をフォトマルチプライヤ等の光検出器で検出することにより、ＱＣファントムの各ファントム画像のＱＬ（クォンタムレベル）を読み取る。こうして、ＱＣファントムの長方形の金属板部分の画像Ｇ1 を含む処理領域の画像Ａ1 のＱＬをデジタル画像データとして取得することができる。
なお、本実施形態では、処理領域は、ＱＣファントムの長方形の金属板部分（領域）を含み、この金属板部分より少し広い領域とするのが良い。
【００１６】
ここで、本実施形態においては、処理領域の画像Ａ1 のＱＬ（画像データ）がログデータの場合には、リニアデータに変換する。ところで、本発明法においては、処理領域の画像Ａ1 のＱＬとしてログデータを用いて縮率の計算を行っても良いが、縮率の計算精度を高めるためには、リニアデータを用いて縮率の計算をするのが好ましい。なお、ＣＲにおいて得られるデジタル画像データは、一般に濃度データ（ログ（Ｌｏｇ）データ）であるので、リニアデータに変換するのが良い。
なお、この時、処理領域の画像Ａ1 のＱＬ（画像データ）に対して、ノイズを除去するために、画像の主／副両方向に所定のメディアンフィルター長でメディアンフィルター処理を行なっても良いし、計算精度を向上させ、ノイズの影響を少なくするため、主／副両方向に所定のマスクサイズ（画素数）の移動平均処理を行なっても良い。
【００１７】
その後、ステップ２において、処理領域内の金属板領域の画像Ｇ1 のＱＬ（画像データ）の平均値Ｄｍと、処理領域内の金属板領域以外のバックグランド領域の画像Ｂ1 のＱＬ（以下、画像データとする）の平均値Ｄｂとを求める（図２（ａ）参照）。
これらの平均値の半値ｍｍ（＝（Ｄｍ＋Ｄｂ）／２）を閾値として、処理領域内において、この閾値以上の画像データ値を１、この閾値未満の画像データ値を０として、２値画像Ａ2 に変換する（図２（ｂ）参照）。この時、上述したように、ＣＲにおいては、金属板領域の画像Ｇ1 は低濃度、バックグランド領域の画像Ｂ1 は高濃度となるので、一般にＤｂ＞ｍｍ＞Ｄｍとなる。このため、バックグランド領域の２値画像Ｂ2 の画像データは１（黒）、金属板領域の２値画像Ｇ2 の画像データは、０（白）となる。
【００１８】
続いて、ステップＳ３において、処理領域の２値画像Ａ2 の主および副両方向（例えば主走査方向および副走査方向）にそれぞれ１画素ずらした２値画像Ａ3 を作成する。この時、処理領域の２値画像Ａ3 は、金属板領域の２値画像Ｇ2 （画像データ＝０）およびバックグラウンド領域の２値画像Ｂ2 （画像データ＝１）を主／副両方向にそれぞれ１画素ずらした金属板領域の２値画像Ｇ3 （画像データ＝０）およびバックグラウンド領域の２値画像Ｂ3 （画像データ＝１）を含む（図２（ｃ）参照）。
【００１９】
こうして得られた処理領域の２値画像Ａ3 を元の２値画像Ａ2 より引き算し、その結果が、１または−１のときは１とし、０のときは０とすることで、金属板領域の矩形状の４辺からなるエッジ画像Ｇ4 を作成する。すなわち、金属板領域の２値画像Ｇ2 とバックグラウンド領域の２値画像Ｂ2 との境界となる金属板領域の矩形状の周辺の４辺を構成する画素のみが１または−１となり、金属板領域の２値画像Ｇ2 の内部領域の画素およびバックグラウンド領域の２値画像Ｂ2 の内部領域の画素では、０となるので、金属板領域の矩形状の４辺を構成する画素の画像データのみが１となるエッジ画像Ｇ4 が作成される（図２（ｃ）参照）。
このように、主／副両方向にずらした画像を元の画像から引くのは、それぞれの方向へ差分をとる（微分する）ことに等しく、データが変化している部分であるエッジを抽出することになる。
【００２０】
次に、ステップＳ４において、こうして作成されたエッジ画像Ｇ4 に対して、例えば、ハフ変換処理を施し、エッジ画像Ｇ4 の４辺を構成する４本の直線の式を求める（図２（ｄ）参照）。なお、これらの４本の直線の式を画像データから求める方法は、ハフ変換処理に限定されず、画像データから直線の式を求めるための従来公知の画像データ処理を用いれば良い。また、この時に、ノイズ等の影響をなくすため、エッジ画像の主／副両方向にメディアンフィルター処理を行なっても良い。
【００２１】
続いて、ステップＳ５において、４本の直線の式から４つの交点Ｐ1 〜Ｐ4 を求める。その後、ステップＳ６では、交点Ｐ1 とＰ2 との間の中点Ｏ1 と、交点Ｐ3 とＰ4 との間の中点Ｏ2 を求める。そして、ステップＳ７において、交点Ｐ1 とＰ2 間の中点Ｏ1 と、交点Ｐ3 とＰ4 間の中点Ｏ2 との距離Ｌを求める（図２（ｅ）参照）。これによって、金属板（試験片）の画像Ｇ1 の中点Ｏ1 Ｏ2 間の長さＬが求められる。
【００２２】
最後に、ステップＳ８において、金属板の画像Ｇ1 の中点Ｏ1 Ｏ2 間の長さＬと金属板の実際の長さＬ0 から画像Ｇ1 の縮率ε（＝（Ｌ−Ｌ0 ）／Ｌ0 あるいは１００×（Ｌ−Ｌ0 ）／Ｌ0 ％）が求められ、画像の縮率計算の全工程が終了する。
なお、上記実施形態では、金属板試験片の画像Ｇ1 の縮率εを画像Ｇ1 の長手方向（または主方向）についてのみ計算しているが、本発明はこれに限定されず、短手方向（または副方向）についても縮率を計算しても良いのはもちろんである。すなわち、縮率計算方向を、長手方向および／または短手方向としても良いし、あるいは主方向および／または副方向としても良い。
【００２３】
以上のように、本発明の第１実施形態では、長さが既知の金属板等の試験片の画像Ｇ1 のエッジを表す４本の直線の式を求め、求められた４本の直線の式からエッジ間の距離を求めるので、画像Ｇ1 の縮率が自動計算によって容易に得られる。
【００２４】
また、試験片の画像Ｇ1 およびバックグラウンドの画像Ｂ1 をそれぞれ２値画像Ｇ2 およびＢ2 に変換した後、２値画像Ｇ2 およびＢ2 をそれぞれ主／副両方向に１画素ずらした２値画像Ｇ3 およびＢ3 を元の２値画像Ｇ2 およびＢ2 から引き算して、その結果が０でないのときは１とすることで、エッジ画像Ｇ4 が容易に抽出される。
さらに、エッジ画像Ｇ4 から４本の直線の式を求めるのにハフ変換を用いたことで、直線の式が容易に求められる。
本発明の第１実施形態の縮率計算方法は、基本的に以上のように構成される。
【００２５】
次に、本発明の第２実施形態の画像の縮率計算方法について説明する。
図１（ｂ）に示す本発明の第２実施形態の画像の縮率計算方法では、まず、図１（ａ）に示す本発明の第１実施形態と同様に、ステップＳ１において、長さが分かっている長方形の金属板を試験片として撮影し、図３（ａ）に示すように、金属板の画像Ｇ1 のＱＬ（画像データ）を得る。
この時、上記第１実施形態と同様に、ＱＬ（画像データ）が、ログデータである場合には、リニアデータとするのが好ましい。
なお、金属板の画像Ｇ1 のＱＬ（画像データ）の取得およびログデータのリニアデータへの変換は、上記第１実施形態と同様に行えば良い。
【００２６】
その後、こうして得られた画像データ（リニアデータ）に対して、ノイズを除去するために、主（走査）方向または副（走査）方向、好ましくは主／副両方向に各ラインに沿って所定のメディアンフィルター長でメディアンフィルター処理を行なうのが好ましい。ここで、メディアンフィルター長は、特に制限的ではなく、画像の縮率の計算精度やノイズに応じて適宜決定すれば良いが、例えば、３〜１１とするのが好ましい。
【００２７】
次に、こうして得られたメディアンフィルター処理済みの画像データに対して、計算精度を向上させ、ノイズの影響を少なくするため、主方向または副方向、好ましくは、主／副両方向に各ラインに沿って所定の画素数（マスクサイズ）の移動平均処理を行なう。ここで、移動平均の画素サイズは、特に制限的ではなく、画像の縮率の計算精度やノイズに応じて適宜決定すれば良いが、例えば、３〜１１とするのが好ましい。この移動平均処理は、画像データの高周波成分を除去することができ、画像データの変化を滑らかにすることができる。その結果、画像の縮率の計算精度を向上させ、ノイズの影響を少なくすることができる。
【００２８】
次に、ステップＳ９において、画像Ｇ1 の縮率の計算方向と直交方向の４つのエッジ点Ｅ1 〜Ｅ4 を求める（図３（ａ）参照）。この後、ステップＳ１０において、２組のエッジ点Ｅ1 とＥ2 との間の中点Ｏ3 およびエッジ点Ｅ3 とＥ4 との間の中点Ｏ4 をそれぞれ求める。
【００２９】
続いて、ステップＳ１１において、中点Ｏ3 およびＯ4 からそれぞれ求めたい縮率計算方向へのエッジ点Ｅ5 およびＥ6 を探す（図３（ｂ）参照）。なお、本実施形態においても、ステップＳ１において、第１実施形態における図２（ａ）に示すように、金属板領域の画像Ｇ1 およびバックグラウンド領域の画像Ｂ1 を含む処理領域全体の画像Ａ1 を含む処理領域の画像Ａ1 のＱＬ（画像データ）を求めておいて、本実施形態における処理エッジ点Ｅ1 〜Ｅ6 を求める各工程において、処理領域内の金属板領域の画像Ｇ1 のＱＬ（画像データ）の平均値Ｄｍと、処理領域内の金属板領域以外のバックグランド領域の画像Ｂ1 のＱＬ（画像データ）の平均値Ｄｂとを求め、次いで、これらの平均値の半値ｍｍ（＝（Ｄｍ＋Ｄｂ）／２）を求め、この半値ｍｍを横切る前後のデータ値を直線近似して、エッジ点Ｅ1 〜Ｅ6 の位置を決定するのが好ましい。
【００３０】
次に、ステップＳ１２において、金属板領域の画像Ｇ1 の軸方向の傾きθを求める。そして、ステップＳ１３において、金属板領域の画像Ｇ1 の軸方向の長さＬを、次式（１）より算出する（図３（ｃ）参照）。
Ｌ＝（Ｌ1 ＋Ｌ3 ）cos θ＋Ｌ2           （１）
ここで、上記式（１）において、Ｌ1 は、エッジ点Ｅ1 とＥ2 との間の中点Ｏ3 とエッジ点Ｅ5 との距離、Ｌ3 は、エッジ点Ｅ3 とＥ4 との間の中点Ｏ4 とエッジ点Ｅ6 との距離であり、Ｌ2 は、中点Ｏ3 と中点Ｏ4 との間の距離である。
【００３１】
こうして、最後に、ステップＳ１４において、金属板領域の画像Ｇ1 の軸方向の長さＬと金属板の実際の長さＬ0 から画像Ｇ1 の縮率ε（＝（Ｌ−Ｌ0 ）／Ｌ0 あるいは１００×（Ｌ−Ｌ0 ）／Ｌ0 ％）が求められ、画像の縮率計算の全工程が終了する。
なお、上記実施形態では、金属板試験片の画像Ｇ1 の縮率εを画像Ｇ1 の長手方向（または主方向）についてのみ計算しているが、本発明はこれに限定されず、短手方向（または副方向）についても縮率を計算しても良いのはもちろんである。すなわち、縮率計算方向を、長手方向および／または短手方向としても良いし、あるいは主方向および／または副方向としても良い。
【００３２】
また、本発明の第２実施形態では、同様に、長さが既知の金属板等の試験片の画像Ｇ1 において、求めたい縮率計算方向と直交する方向のエッジ点Ｅ1 とＥ2 との間の中点Ｏ3 およびエッジ点Ｅ3 とＥ4 との間の中点Ｏ4 をそれぞれ求め、中点Ｏ3 からエッジ点Ｅ5 までの距離Ｌ1 、中点Ｏ4 からエッジ点Ｅ6 までの距離Ｌ3 、中点Ｏ3 とＯ4 との間の距離Ｌ2 および試験片の画像Ｇ1 の軸方向の傾きθから、試験片の画像Ｇ1 の長さＬが求められ、画像Ｇ1 の縮率が自動計算によって得られる。
【００３３】
さらに、試験片領域内の画像Ｇ1 の画像データの平均値Ｄm と試験片領域外のバックグラウンド領域の画像Ｂ1 の画像データの平均値Ｄb とを求め、これらの平均値Ｄm およびＤb の半値ｍｍを横切る前後の値を直線近似してエッジ位置を決定するので、エッジ位置の推定誤差が小さくなる。この場合、画像データに対して、試験片端と直交する方向に、メディアンフィルター処理を行えば、ノイズが除去でき、さらに、移動平均処理を行なえば、ノイズがより除去され、計算精度がさらに向上する。
本発明の第２実施形態の縮率計算方法は、基本的に以上のように構成される。
【００３４】
以上、本発明の第１および第２実施形態の縮率計算方法について詳細に説明したが、本発明はこれらの実施形態に限定されるものではなく、本願の特許請求の範囲に記載されている発明の要旨を逸脱しない範囲で、設計において種々の変更ができるものである。
例えば、上記各実施形態では、試験片として長さが既知の長方形の金属板を用いたが、本発明はこれに限定されず、長さが既知の矩形の金属板であればよく、特に、正方形の金属板を用いれば、主副の長さが同時に求められ、高効率の縮率計算が可能となる。
【実施例】
【００３５】
本発明法を実施例に基づいて具体的に説明する。
本発明の第２実施形態の画像の縮率計算法を図１（ｂ）に示すフローチャートに従って行った。
ＦＣＲ（富士写真フィルム株式会社製）において、図３（ａ）に示すように、長さが既知の長方形の金属板を含む試験片（ＱＣファントムの一部）が撮影された蓄積性蛍光体シートを光走査して金属板領域の画像Ｇ1 、さらには金属板領域の画像Ｇ1 および金属板領域外のバックグラウンド領域の画像Ｂ1 を含む処理領域の画像Ａ1 （図２（ａ）参照）のＱＬ（デジタル画像データ）を取得し（ステップＳ１参照）、得られたデジタル画像データを用いて、本発明法のアルゴリズムを行い、撮影した試験片の画像Ｇ1 の長さを求め、主および副（走査）方向の縮率を確認した。
【００３６】
ここで、ＥＤＲ条件  Ｓｅｍｉ、Ｌ＝１とし、メディアンフィルター処理の際のメディアンフィルター長Ｆを、Ｆ＝３とした。
また、主方向および副方向について、処理領域は、主方向および副方向開始位置（主Ｍｓ、副Ｓｓ）から主／副方向終了位置（主Ｍｅ、副Ｓｅ）とし、バックグラウンド領域１は、主／副方向バックグラウンド１開始位置（主Ｍｂ１ｓ、副Ｓｂ１ｓ）から主／副方向バックグラウンド１終了位置（主Ｍｂ１ｅ、副Ｓｂ１ｅ）とし、バックグラウンド領域２は、主／副方向バックグラウンド２開始位置（主Ｍｂ２ｓ、副Ｓｂ２ｓ）から主／副方向バックグラウンド２終了位置（主Ｍｂ２ｅ、副Ｓｂ２ｅ）とし、計算領域１は、主／副方向計算１開始位置（主Ｍｍ１ｓ、副Ｓｍ１ｓ）から主／副方向計算１終了位置（主Ｍｍ１ｅ、副Ｓｍ１ｅ）とし、計算領域２は、主／副方向計算２開始位置（主Ｍｍ２ｓ、副Ｓｍ２ｓ）から主／副方向計算２終了位置（主Ｍｍ２ｅ、副Ｓｍ２ｅ）とした。
また、出力データは、主方向画素数（ｉｎｔ）および副方向画素数（ｉｎｔ）とした。
【００３７】
ここで、本発明法のアルゴリズムを、最初に主方向について、次に副方向について主方向と全く同様に、図１（ｂ）に示すフローチャートに従って、以下のように実施した。
まず、ＣＲにおいて得られるデジタル画像データは、濃度データ（Ｌｏｇデータ）であるので、処理領域のＬｏｇデータＤｏ「」「」を下記式に従ってリニアデータｆｌｏａｔＤ「」「」に変換した。ここで、「」「」は、主走査方向（ｘ）のデータの配列および副走査方向（ｙ）のデータの配列を示す。
Ｄ「」「」＝１０＾（Ｄｏ「」「」／（（２＾画像ｂｉｔ）／Ｌ））
【００３８】
次に、ノイズ等を除去し、ノイズ等の影響を無くすために、こうして得られたリニアデータＤ「」「」に対して主／副両方向にメディアンフィルター処理を行った。
続いて、計算精度を向上させ、ノイズ等の影響を無くすために、メディアンフィルター処理されたリニアデータＤ「」「」に対して主／副両方向に移動平均処理を行った。
【００３９】
その後、ステップＳ９において、バックグラウンド領域１の平均値（Ａｂ）と計算領域１の平均値（Ａｏ）を求め、半値Ｈ＝（Ａｂ＋Ａｏ）／２を求めた。
次に、計算領域１の（（Ｍｍ１ｓ＋Ｍｍ１ｅ）／２，（Ｓｍ１ｓ＋Ｓｍ１ｅ）／２）をスタートポイントとして、副走査方向の両方向へ探索していき、初めて半値Ｈを下回る画素（（Ｍｍ１ｓ＋Ｍｍ１ｅ）／２，Ｓ１ｔ），（Ｍｍ１ｓ＋Ｍｍ１ｅ）／２，Ｓ１ｂ）を求めた。
続いて、ステップＳ１０において、計算領域１の中心画素（（Ｍｍ１ｓ＋Ｍｍ１ｅ）／２，（Ｓ１ｔ＋Ｓ１ｂ）／２）を求め、ステップＳ１１において、この中心画素をスタートポイントとして、求めたい縮率計算方向である主走査方向へ検索していき初めて半値Ｈを下回る画素（Ｍ１，（Ｓ１ｔ＋Ｓ１ｂ）／２）を求めた。
【００４０】
同様にして、ステップＳ９〜１１を繰り返して、計算領域２の中心画索（（Ｍｍ２Ｓ＋Ｍｍ２ｅ）／２，（Ｓ２ｔ＋Ｓ２ｂ）／２）を求めるとともに、この点をスタートポイントとして、主走査方向へ探索していき、初めて半値Ｈを下回る画素（Ｍ２，（Ｓ２ｔ＋Ｓ２ｂ）／２）を決定した。
続いて、ステップＳ１２において、傾き角θ＝ａｒｃｔａｎｇｅｎｔ（（（Ｓ２ｔ＋Ｓ２ｂ）／２−（Ｓ１ｔ＋Ｓ１ｂ）／２）／（（Ｍｍ２ｓ＋Ｍｍ２ｅ）／２−（Ｍｍｌｓ＋Ｍｍｌｅ）／２））を求めた。
【００４１】
この後、ステップＳ１３において、出力データとして、金属板試験片の画素数（ｉｎｔ）（小数点以下は四捨五入）を下記式によって求めた。
（（Ｍｍｌｓ＋Ｍｍｌｅ）／２‐Ｍｌ）×ｃｏｓθ＋√（（（Ｓ２ｔ＋Ｓ２ｂ）
／２−（Ｓ１ｔ＋Ｓ１ｂ）／２）2 ＋（（Ｍｍ２ｓ＋Ｍｍ２ｅ）／２−
（Ｍｍｌｓ＋Ｍｍｌｅ）／２）2 ）＋（Ｍ２−（Ｍｍ２ｓ＋Ｍｍ２ｅ）／２）
×ｃｏｓθ
最後に、得られた主／副画素数金属板試験片
【図面の簡単な説明】
【００４２】
【図１】（ａ）および（ｂ）は、それぞれ本発明の第１および第２実施形態に係る画像の縮率計算方法を示すフローチャートである。
【図２】（ａ）〜（ｅ）は、それぞれ図１（ａ）に示す本発明の第１実施形態に係る画像の縮率計算方法を示す工程説明図である。
【図３】（ａ）〜（ｃ）は、それぞれ図１（ｂ）に示す本発明の第２実施形態に係る画像の縮率計算方法を示す工程説明図である。
【符号の説明】
【００４３】
Ｄb ，Ｄm   平均値
Ｅ1 ，Ｅ2 ，Ｅ3 ，Ｅ4 ，Ｅ5 ，Ｅ6   エッジ点
Ａ1 ，Ｂ1 ，Ｇ1   画像
Ａ2 ，Ａ3 ，Ｂ2 ，Ｂ3 ，Ｇ2 ，Ｇ3   ２値画像
Ｇ4   エッジ画像
Ｌ、Ｌ1 、Ｌ2   距離
Ｏ1 ，Ｏ2 ，Ｏ3 ，Ｏ4   中点
Ｐ1 ，Ｐ2 ，Ｐ3 ，Ｐ4   交点
【図１】
2005058315.tif 000003
【図２】
2005058315.tif 000004
【図３】
2005058315.tif 000005