"試料を互いに異なる高さから撮像して、焦点部位の異なる多数の2次元試料画像を得るステップと、前記試料画像に対して、離散ウェーブレット変換を行うステップと、それぞれの試料画像に対して離散ウェーブレット変換を行うことにより得られた詳細サブーバンドにおいて、各ピクセルの詳細サブーバンド係数値を比較して、最大係数値を示す画像の撮像高さで、各ピクセルに対する初期の高さ地図を作成するステップと、それぞれの入力画像において、離散ウェーブレット変換の行われた画像の近似サブーバンド領域に対して、ピントの度合い焦点整合度を計算するステップと、前記焦点整合度にフィルターをかけて、合焦点のピクセル境界点の高さ情報は残留させ、非合焦点のピクセル非境界点の高さ情報は除去するステップと、前記フィルターを通過できずに除去されたピクセル非境界点の高さを、フィルターを通過したピクセル境界点の高さ値から内挿して算出するステップと、前記高さ情報の除去されたピクセルに対して前記内挿によって算出された高さを代入して高さ地図を作成するステップと、を含む3次元プロファイルマップの作成方法。前記初期の高さ地図を作成するステップ後に、各ピクセルを中心とする検証ウィンドウ内において、同じ高さとして計算されたピクセルの詳細サブーバンドウェーブレット係数を加算し、高さに応じたウェーブレット係数の和を比較して、最も大きな値を与える高さを当該ピクセルの実際の高さで補正するステップをさらに含む、請求項1に記載の3次元プロファイルマップの作成方法。前記検証ウィンドウは、各ピクセルを中心として3x3のピクセルサイズを有するものである、請求項2に記載の3次元プロファイルマップの作成方法。前記詳細サブーバンドウェーブレット係数が加算される、すなわち、グループ化されるピクセルを選定する方法は、検証ウィンドウ内において同じ高さに相当するピクセルを同じグループとして選定することである、請求項2に記載の3次元プロファイルマップの作成方法。前記近似サブーバンド領域の焦点整合度は、次式1に示すように、当該ピクセル座標:x,yの明るさIx,yと、周辺4点座標:x,y1、x,y-1、x1,y、x-1,yの明るさとの差の絶対値の和として計算される、請求項1に記載の3次元プロファイルマップの作成方法。前記焦点整合度が、隣り合う両ピクセルの境界が現れる明るさの差の3倍未満であれば、非合焦点のピクセルと見なして高さ情報を除去する、請求項5に記載の3次元プロファイルマップの作成方法。前記内挿される点の高さは、境界点の高さの重み付け平均から求められ、前記重み付け平均の重み付け値関数は、距離逆数1/r、ここで、rは、内挿される点と境界点との間の距離を表わす。の関数であるか、あるいは、距離のガウス関数である、請求項1に記載の3次元プロファイルマップの作成方法。前記重み付け平均の重み付け値関数は、1/rまたはexp-r2/1000である、請求項7に記載の3次元プロファイルマップの作成方法。前記非境界点の高さを内挿するための境界点の選択は、内挿される点境界点を中心として同心円または同心四角形を描き、中心から外方に向けて検索するが、境界点が見出されれば、見出された境界点から±30°の範囲境界点が見出された方位角内においてはそれ以上の境界点を検索しないような方式により行われる、請求項1に記載の3次元プロファイルマップの作成方法。本発明の目的は、試料画像において、合焦点の領域の正確な位置情報を把握して、位置に応じた高さ変化を正しく反映する、3次元プロファイルマップの作成方法を提供することである。本発明の他の目的は、合焦点の画像部分を有さない非合焦点の領域が関心領域の情報把握を妨げず、周辺の合焦点区間とは高さの連続性を有する、3次元プロファイルマップの作成方法を提供することである。本発明のさらに他の目的は、画像の焦点整合度を正確に反映することのできる、3次元プロファイルマップの作成方法を提供することである。前記目的を達成するために、本発明は、試料を互いに異なる高さから撮像して、焦点部位の異なる多数の2次元試料画像を得るステップと、多数の前記試料画像に対して、離散ウェーブレット変換DiscreteWaveletTransform:DWTを行うステップと、それぞれの試料画像に対して離散ウェーブレット変換DWTを行うことにより得られた詳細サブーバンドsubbandにおいて、各ピクセルの詳細サブーバンドsubband係数値を比較して、最大係数値を示す画像の撮像高さで、各ピクセルに対する初期の高さ地図を作成するステップと、それぞれの入力画像において、離散ウェーブレット変換DWTの行われた画像の近似サブーバンドsubband領域に対して、ピントの度合い焦点整合度を計算するステップと、前記焦点整合度にフィルターをかけて、合焦点のピクセル境界点の高さ情報は残留させ、非合焦点のピクセル非境界点の高さ情報は除去するステップと、前記フィルターを通過できずに除去されたピクセル非境界点の高さを、フィルターを通過したピクセル境界点の高さ値から内挿して算出するステップと、前記高さ情報の除去されたピクセルに対して前記内挿によって算出された高さを代入して高さ地図を作成するステップと、を含む3次元プロファイルマップの作成方法を提供する。本発明に係る3次元プロファイルマップの作成方法によれば、互いに異なる高さ情報とともに入力された多数の入力画像を用いて、合焦点の全領域から、試料の当該部分の高さ情報を得、非合焦点の部分に対して不連続的な高さを入力することに伴う、不正確な高さプロファイルの問題がない、整合性の良い高さ地図を提供することができる。本発明に従い製作された3次元プロファイルマップを用いて、試料の正確な測定および観察が可能になり、高さに応じて焦点部位が異なってくる試料の画像を、整合性よく1枚の合成画像として合成することができる。以下、添付図面に基づき、本発明を詳述する。図1は、本発明に係る3次元プロファイルマップの作成方法が実現可能な光学顕微鏡の構成ブロック図である。同図に示すように、本発明において採用可能な光学顕微鏡は、試料に照明光を照射する照明と、試料の画像を拡大する光学鏡筒および試料画像を結像するカメラを備え、試料の拡大された画像を得る顕微鏡部10と、前記顕微鏡部10の高さ、すなわち、焦点距離を調節する顕微鏡運転制御部12と、前記顕微鏡部10から試料の拡大された画像を受け取ってデジタル画像信号に変換する画像入力部22と、前記顕微鏡運転制御部12から試料の画像が得られる高さ顕微鏡の高さが入力される高さ入力部24、および前記画像入力部22からデジタル画像信号が入力され、前記高さ入力部24から試料画像が得られた高さ情報が入力されて、入力された画像の焦点整合度を評価し、試料の3次元プロファイルマップ高さ地図を生成する画像処理部26と、を備える。なお、前記画像処理部26は、前記3次元プロファイルマップを用いて、合焦点の画像からなる、深度の拡張された画像を生成することができる。図2は、本発明の一実施形態による3次元プロファイルマップの作成方法アルゴリズムを説明するためのフローチャートである。図1および図2に示すように、本発明に係る3次元プロファイルマップの作成方法は、試料を互いに異なる高さから撮像して、焦点部位の異なる多数の2次元試料画像を得るステップS10と、前記試料画像に対して、離散ウェーブレット変換DWTを行うステップS20と、それぞれの試料画像に対して離散ウェーブレット変換DWTを行うことにより得られた詳細サブーバンドsubbandにおいて、各ピクセルの詳細サブーバンド係数値を比較して、最大係数値を示す画像の撮像高さで、各ピクセルに対する初期の高さ地図を作成するステップと、それぞれの入力画像において、離散ウェーブレット変換DWTの行われた画像の近似サブーバンドsubband領域に対して、焦点整合度を計算するステップと、前記焦点整合度にフィルターをかけて、合焦点のピクセルの高さ情報は残留させ、非合焦点のピクセルの高さ情報は除去するステップと、前記フィルターを通過できずに除去されたピクセルの高さを、フィルターを通過したピクセルの高さ値から内挿して算出するステップと、前記高さ情報の除去されたピクセルに対して前記内挿によって算出された高さを代入して高さ地図を作成するステップと、を含む。ここで、前記ステップとステップとの間に、各ピクセルを中心とする検証ウィンドウ内において、同じ高さとして計算されたピクセルの詳細サブーバンドウェーブレット係数を加算し、高さに応じたウェーブレット係数の和を比較して、最も大きな値を与える高さを当該ピクセルの実際の高さで補正するステップおよび/または前記補正済みの高さ地図に中央値フィルターをかけて、インパルスノイズの除去された高さ地図を作成するステップをさらに含むことが好ましい。ここで、前記試料画像を得るステップS10は、顕微鏡部10によって行われ、得られた試料画像は、画像入力部22においてデジタル画像信号に変換された後、画像処理部26において離散ウェーブレット変換DWTなどの残りの過程が行われても良い。必要に応じて、前記カラーデジタル画像信号は、画像処理部26において、通常の方法によって黒白のデジタル画像信号に変換された後S12、離散ウェーブレット変換DWTなどの残りの過程が行われてもよい。このような黒白画像の変換は、今後の画像信号データの数を減らして、画像処理部26にかかる計算の負担を低減することができる。以下、本発明に係る3次元プロファイルマップの作成方法をより具体的に説明する。本発明に係る3次元プロファイルマップの作成方法は、a試料画像の全体に対して高さ情報を得るステップ初期の高さ地図の作成と、b試料の境界線境界点に相当する領域の情報を抽出し、前記ステップaにおいて求められた境界線領域の高さ情報を用いて境界線以外の領域の高さを内挿法により算出するステップ、といった2つのステップに大きく分けられる。前記ステップaは、合焦点の領域に関する正確な高さ情報を得ることが目的である。本発明は、前記ステップaにおいてより一層正確な高さ情報を得、ステップbにおいて高さ情報の整合性を判断し、整合性の低いピクセルの高さ情報を整合性よく補正して、非合焦点の領域を処理する。前記ステップaにおいて、画像の高さ情報を得るために、まず、試料を互いに異なる高さ距離から撮像して、焦点部位の異なる多数例えば、2~8枚、好ましくは、2~5枚の2次元画像を得S10、前記2次元画像から、合焦点部分の位置情報と高さ情報を取り合わせて一つの高さ地図例えば、xy平面上に当該位置の高さ値入力画像番号が表示された形のデータである2次元地図。区別し易くするために、カラーごとに高さを表示する場合もある。を作成する。このために、得られた2次元画像の各ピクセルに対して、i焦点整合度の評価、ii焦点整合度焦点情報からの各ピクセルの高さ情報の算出、およびiii取り合わせられた高さ情報を用いた高さ地図の作成といった3ステップを行う。まず、i試料画像の各部分に対する焦点整合度の評価方法は、下記の通りである。一般に、合焦点の画像は境界がはっきりしており、最もピントの合った合焦点が最高の位置において、光の強さと明暗対比が最大となるため、ピクセル間の明暗対比を用いる境界検出方法により焦点整合度を評価することができる。本発明において、境界検出方法として、通常の信号または画像の分析方法である離散ウェーブレット変換DWT法を採用する。離散ウェーブレット変換DWT法は、局所化したウェーブレットを用いて信号を分析する方法であり、多層構造分析によって、様々な解像度の信号分析を行うことができ、局所的な信号の分析に役立つ。特に、離散ウェーブレット変換DWTを繰り返し行うことにより多重離散ウェーブレット変換、上位レベルに上がるにつれて、平均値平滑化効果が高まるため、同型ノイズなどの非インパルス型ノイズが除去されるというメリットがあるが、インパルス型ノイズは依然として残留するため、高さ地図が整合性を有するように後加工してインパルス型ノイズを除去する。試料の2次元画像に対して離散ウェーブレット変換DWTを行うとS20、図3AのAに示すように、LL、HL、LH、HHからなる4個のサブーバンドに画像が分けられる。ここで、LLは解像度が1/4に減少された原画像の縮小画像であり、近似サブーバンドApproximationSub-bandと呼ばれ、HLサブーバンドは、垂直方向の画像情報の変化の度合い例えば、原画像に対して離散ウェーブレット変換DWTを1回行った場合、HLサブーバンドの係数は垂直方向への画像の明るさの変化の度合いを示し、離散ウェーブレット変換DWTを2回以上行った場合、HLサブーバンドの係数は縮小された画像の情報変化を表わす。、LHは、水平方向の画像情報の変化の度合い、HHは、対角線方向の画像情報の変化の度合いを表わす。例えば、LHが大きいということは、水平方向への画像情報の変化が大きいことを意味するため、垂直方向の境界が存在することを意味する。前記LH、HL、HHを詳細サブーバンドDetailSub-bandと称する例えば、米国特許6151415号などを参照されたい。かかる離散ウェーブレット変換DWTは、原画像の解像度に応じて、1枚の原画像に対して2回図3AのBおよび図3B、3回図3AのCなど、多数回に亘って行うことができる。既に離散ウェーブレット変換DWTの行われた画像に重ねて離散ウェーブレット変換DWTを行う場合、現在得られたLLサブーバンドに対して離散ウェーブレット変換DWTを行うことにより、図3AのBまたはCなどの多層マルチレベル構造の画像を得る。このように、離散ウェーブレット変換DWTを行うことにより得られた詳細サブーバンドのウェーブレット係数例えば、図3Bには、16個のピクセルに対するウェーブレット係数が示してある。は、画像の明暗対比の特性を反映するため、詳細サブーバンドのウェーブレット係数のうち、絶対値が最大となる値を取り合わせると、境界検出が可能となる。なお、近似サブーバンド値は、ノイズの除去された低解像度の画像であるため、これを用いても、境界検出などの画像処理を行うことが可能である。すなわち、本発明は、従来より画像の境界検出や信号分析に用いられてきたフーリエ分析の代わりに、様々な解像度で局所的なデータアクセスが行える離散ウェーブレット変換DWT法を応用して、焦点整合度を評価計算する。このような離散ウェーブレット変換DWTを用いた焦点整合度の評価方法の一例は、本出願人による大韓民国特許出願第10-2010-0054861号出願日:2010年06月10日にも開示されており、前記特許出願の内容は、参照としてこの明細書中に取り込まれる。例えば、入力解像度が1600x1200の画像から、解像度が約100x100の画像を得るために、約3回に亘って離散ウェーブレット変換DWTを行い、離散ウェーブレット変換DWTを行うことにより得られた詳細サブーバンドLH、HL、HHのピクセル係数の最大値すなわち、LH、HL、HHサブーバンド係数値のうちの最大値を、当該画像における各ピクセルの焦点整合度であると判断する。同様に、入力されたそれぞれの画像に対して、ピクセル係数の最大値を得る。本発明において、離散ウェーブレット変換DWTは、1回以上の多重に行われ、変換回数は、原画像の解像度、合成されるピクセルの大きさなどに応じて、1~10回、好ましくは、2~5回行われる。次いで、各ピクセルに対して、前記ピクセルの係数値が最大となる画像の高さ撮像または焦点距離を、当該ピクセルに対する合焦点の位置高さであると判断し、当該ピクセルにこの高さを指定する。これは、前記ii焦点整合度焦点情報から各ピクセルの高さ情報を算出する過程である。つまり、画像の各ピクセルに対する焦点整合度焦点情報から、各ピクセルの高さ情報を得るのである。次いで、iii取り合わせられた高さ情報を用いて高さ地図を作成するステップにおいて、ピクセル単位で画像の高さ地図を作成し、このようにして作成された初期の高さ地図の一例が図4に示してある。図4は、4枚の金属板と円形の金属構造物がほとんど同じ高さに構成されたサンプルの初期の高さ地図であり、前記円形の金属構造物の中央部図4における真ん中の黒色部が部分的に凹んでいる窪んでいるような構造を有する。このようにして得られた初期の高さ地図においては、上述したように、インパルス型ノイズが除去されていない状態であり、ピクセル単位の明暗対比値は外部の環境要素に敏感であるため、得られた高さの整合性を確認することが好ましい。ピクセル単位の敏感度を緩和しつつも、ウェーブレット係数の境界値情報を保持するには、ピクセルの周辺にウィンドウを設定し、このウィンドウ内においてウェーブレット係数を比較することが有効であるAwavelet-basedimagefusiontutorial,G.PajaresandJ.M.delaCruz,PatternRecognition3720041855~1872参照。かようなウィンドウ基盤の合成方法は、高さ地図を作成するに際して焦点整合度を判断するために用いられるか、あるいは、整合性を確認するに際して単にフィルターを適用するために用いられる。しかしながら、このようなウィンドウ基盤の方式を用いると、全ての入力画像の全てのピクセルに対してウィンドウ演算が必要となるため、速度が遅くなってしまう。このため、本発明においては、ウェーブレット係数の意味焦点整合度を保持しつつ、計算が簡単に行われるように、ウィンドウ例えば、当該ピクセルを中心とする3x3ピクセルサイズの検証ウィンドウ内のウェーブレット係数を入力高さに応じてグループ化する。一般に、1つ若しくは2つのピクセルのみ特異的な高さを有する試料は多くないため、ウィンドウ中に境界が存在すれば、ほとんどの場合は境界点ではない境界線の形で存在し、同じ高さを有する境界線のウェーブレット係数の和は、異なる高さを有するウェーブレット係数の和よりも大きい。たとえウィンドウ内のいずれかのピクセルが周辺と高さが大きく異なるインパルスノイズの形を帯びて大きいウェーブレット係数を有するとしても、ウィンドウ中に境界が存在すれば、一般に境界点の数がインパルス型ノイズよりも多く、境界点のウェーブレット係数は大きな値を有するため、同じ高さにある境界点のウェーブレット係数の和は、普通、同じ高さにあるインパルスノイズのウェーブレット係数よりも大きい。このため、本発明においては、ウィンドウ内において同じ高さとして計算されたピクセルをグループ化し、グループ化されたピクセルのウェーブレット係数を加算し、高さに応じたウェーブレット係数の和を比較して、最も大きな値を与える高さを当該ピクセルの実際の高さで補正するS24。以下、この方式をグループ化した係数和検証GroupedCoefficientsSummationVerification:GCSV方式と呼ぶ。検証ウィンドウ内において、同じ高さのウェーブレット係数のグループ化は、前記ステップiiにおいて得られた当該ピクセルの入力画像の番号高さを基準としてグループ化される。すなわち、同じ高さのウェーブレット係数は、同じグループにくくる。例えば、互いに異なる3つの高さ1、2および3から試料を撮像して、焦点部位の異なる3つの2次元画像を得た場合、高さ1に相当するウェーブレット係数を同じグループにくくり、高さ2に相当するウェーブレット係数を同じグループにくくり、高さ3に相当するウェーブレット係数を同じグループにくくる。なお、上述した方法だけでは、インパルスノイズが完全に除去されないため、中央値フィルターリングメディアンフィルターリングを1回以上行うことにより、インパルスノイズの除去された1次補正済み高さ地図を作成することが好ましいS26。前記中央値フィルターリングを行うためには、フィルターをかけようとするピクセルを中心として適当な領域のウィンドウ例えば、5x5ピクセルを設定し、前記ウィンドウ内に存在する全ての高さ情報を昇順若しくは降順に並べ替え、その中央値で最初に選択したピクセルの高さ情報を補正取り替えする。このように、境界領域に対する正確な高さが反映されて、ピクセルの高さが周辺ピクセルの高さと整合されるように1次補正された画像地図の一例を図5に示す。図4は、離散ウェーブレット変換DWT後に、単に焦点整合度を用いてプロットした高さ地図であり、図5は、図4の高さ地図に対して、前記グループ化された係数和検証方式GCSVと中央値フィルターリングを行うことにより得られた高さ地図である。図5は、図4に比べて、隣り合うピクセル間の高さのばらつきが低減されていることが分かる。このようにして作成された1次補正済み高さ地図は、図5に示すように、非合焦点の領域において高さのばらつきが激しく現れ、高さ情報が不連続的に入力されるため、高さの変化を正しく表示することができないという問題がある。高さ情報が不連続的に離散化して入力され、実際の試料が図6に示すように連続的な高さ分布を有すると、図7横軸は試料の位置、縦軸は試料の高さを表すに示すように、図6の横軸に平行な点線よりも高い個所はH値、それよりも低い個所はL値を有する高さ分布を示す。図6および図7中、グラフは実際のサンプルのプロファイルを示すものであり、ボックス内の数は、デジタルの形で入力されたサンプルの高さ情報であり、例えば、0はL、1はHに相当し、中央値はHに設定してもよく、Lに設定してもよいため、0.5という値を使用した。ここで、一旦中央値が決定された限り、これを変えることはできない。このため、この方式では、試料が階段状の勾配を有するか、あるいは、連続的な勾配を有するかを区別することはできない。ところが、勾配が不連続的であるピクセルの位置を比較すれば、連続的な高さの変化と不連続的な高さの変化とを区別することができる。図7中、不連続的な勾配の位置は、階段が始まる1点であるか、あるいは、階段の勾配がやや緩やかな場合、階段が上がり始まる点と、階段が上がり終わる点といった2点である。連続的に高さが変化する試料の場合にも、高さが上がり始まる点と、それ以上上がらない点といった2点である。2つの場合を比較すれば、不連続的な勾配が現れる点が近づいた場合には勾配が急峻であり、そうではない場合には勾配が緩やかである。このため、勾配の不連続点が勾配を判断する上で重要な根拠となる。高さ分布の判断に際して重要な根拠となる高さ勾配の不連続点は、焦点整合度により把握することができる。このために、それぞれの入力画像において、離散ウェーブレット変換DWTの行われた、好ましくは、最高レベルHighestlevelまで離散ウェーブレット変換DWTが行われ、必要に応じて、非インパルス型ノイズが除去された画像の近似サブーバンド領域に対して、焦点整合度を計算するS28。図8は、L、Hといった2種類の高さを測定するとき、焦点が高さLに結ばれる3つの光と、高さが連続的に変わる試料を示す図である。Lの位置に達した光1、2は焦点が正確に合ったため、焦点整合度が高く評価されるのに対し、ピントが合っていない非合焦点3の光は焦点整合度が低く評価される。前記近似サブーバンド領域において、特定のピクセルの焦点整合度は、次式1に示すように、当該ピクセル座標:x,yの明るさIx,yと、周辺4点座標:x,y1、x,y-1、x1,y、x-1,yの明るさとの差の絶対値の和として計算され得る。このとき、前記焦点整合度にフィルターをかけて、例えば、閾値を設定して、高域通過フィルターをかけるとS30、合焦点の点1、2の点は残留し、非合焦点の点3の点は消える。例えば、隣り合う両ピクセルの明るさintensityの差が40以上であれば、ピクセル間の境界がはっきりと現れるが、このような明るさの差を示すピクセルが上下左右の4方向に少なくとも2つ存在するとき、すなわち、フィルターをかけようとするピクセルが孤立されていない線としての意味を有するときを閾の条件とする。実験の結果、この場合、閾値が80~120の場合に、輪郭線の形態のみが全く残る。このため、この場合の閾値は80、好ましくは、100に設定することができる。よって、本発明においては、前記焦点整合度が、隣り合う両ピクセルの境界が現れるすなわち、明らかになる明るさの差の3倍未満、好ましくは、2.5倍未満、さらに好ましくは、2倍未満であれば、非合焦点のピクセルと見なして、高さ情報を除去する。ここで、隣り合う両ピクセルの境界が現れる明るさの差は、画像の輪郭線が現れる明るさの差を実験的に決定することができる。このように、焦点整合度にフィルターをかけると、残留した点は境界がはっきりしている境界点であり、ここには勾配が不連続的な点図8における2番の位置が含まれる。このとき、前記閾値をあまりにも低く設定すれば例えば、50、非合焦点の部分が残留し図9参照、前記閾値を適度に設定すれば例えば、100、境界点は残留し、非合焦点の領域の点を消去することができる図10参照。図9および図10において、円形の試料画像内の白抜き部は境界線として認識されて消去されずに残留した点であり、黒色部は境界線以外の部分として認識されて消去された部分である。図9に示すように、未消去点白抜き部が多ければ、画像は全体的に明るくなるが、境界としての意味は薄くなり、図10においては、境界以外の点が十分に除去された状態を示す。前記ピクセルの高さ情報の消去において、境界に相当するピクセルの高さ情報のみ残留して、画像の境界線がいずれも鮮やかに見えることが好ましく、もし、一部の境界情報が失われるとしても、全体的な境界情報は残留し、且つ、境界線以外のピクセルはいずれも消去されることが好ましい。観察される試料の画像、高さの変化などによるが、焦点が合わなくて高さ情報が消去されるピクセルの数の非限定的な例は、全体のピクセルに対して10~70、具体的に20~50、さらに具体的には30~40である。ここで、閾値を通過した境界点がいずれも勾配が連続していない点ではないものの、図11および図12に示すように、同じ高さに存在し、且つ、隣り合う境界点の集まりの周縁にある点は、勾配が連続していない点である。具体的に、図11および図12中、1、2、3、4は境界点であり、1′、2′、3′、4′は高域通過フィルターによって消去された点である。このように、前記フィルターを通過できずに消去された点の高さを、フィルターを通過した点の高さ値例えば、前記ステップS26において求められた1次補正済みの高さ値を内挿して算出しS32、これを用いて最終的に高さ地図を完成する。高さの不連続点は高さ分布を分析する上で重要な役割を果たすが、実際に高さ地図を作成するに際しては、別途に抽出を行う必要がない。境界点は高さ値をそのまま使用し、消去された点非境界点は、周辺において相対する最も近い2つの境界点の高さを内挿して、高さを求めることにより、勾配の変化に関する情報が十分に得られるためである。例えば、図11および図12において、1′のように、境界点ではないものの、境界点として残留した点と同じ高さにある点は、周辺にある境界点の内挿を通じて、周辺境界点と同じ高さを有する。これに対し、2′のように、高さの異なる境界点の間にある点は、周辺の境界点の高さから内挿して、中間の高さを有する。このとき、内挿される点の高さは、境界点の高さの重み付け平均を求めるものであり、重み付け値は近い境界値に一層大きく影響されるため、距離逆数1/rの関数または距離のガウス関数として与えられる。例えば、ガウス関数と重み付け値関数を1/rで展開したときここで、rは、内挿される点と境界点との距離を表わす。、最も低い次数である1/rを重み付け値として適用するか、あるいは、ガウス関数の場合、exp-r2/1000を重み付け値として用いて、好適な結果を得ることができた。一方、図6及び図7などに示す1次元試料とは異なり、2次元の実際の試料においては、境界値が閉曲線の形だけではなく、開曲線、点の集まりなどの形で存在してもよい。このため、本発明においては、高さが算出されるべき消去された点から最も近い個所に位置している3~7個、例えば、5個の境界点から内挿を行うことが好ましい。内挿に用いられる境界点の検出には、内挿される点消去された点を中心とする同心円の検索と方位角の検索とを組み合わせて用いる。同心円の検索と方位角の検索が組み合わせられた方式は、図13に示すように、内挿される点を中心として同心円または同心四角形を描き、中心から外方に向けて検索を行うが、境界点が見出されれば、見出された境界点から、例えば、±30°、好ましくは、±15°の範囲すなわち、境界点が見出された方位角内においてはそれ以上の境界点を検索しない方式である。図13中、中空の四角形は、内挿される点の位置、中実の四角形は境界点、点線の四角形は、検索する同心円、薄色の点線で描かれた三角形の部分は、境界点の見出後に検索の対象から除外された方位角の領域を示す。前記方法は、内挿される点を含む単一閉曲線を検索する方法であり、1次元である図11において2′の点を内挿するに際して、左側に検索する際に境界点2を検索したならば、左側にはそれ以上の検索を行わないことと同じ意味である。このように、境界点が検索されれば、1/r若しくはexp-r2/1000を重み付け値に設定し、内挿を行うことにより、消去された点の高さを算出し、前記高さ情報が除去されたピクセルに対して前記内挿によって算出された高さを代入して最終的に高さ地図を完成するS34。図14は、本発明に従い最終的に完成された3次元プロファイル高さマップの一例を示す。このように、試料画像の最終的な高さ地図が完成されれば、当該高さの離散ウェーブレット変換DWTのサブーバンドで最高レベルのウェーブレットサブーバンド値を再構成し、R、G、Bそれぞれのカラーに対して、最高レベルのウェーブレットの近似サブーバンドを重み付け平均として、柔らかく連結した後、多重逆離散ウェーブレット変換InverseDWTを行うことにより、入力画像の解像度レベルで最終的な合成画像が得られる。生成された最終的な合成画像は、画像の全領域に亘ってピントの合った合焦点、深度の拡張された2次元画像となる。"