"被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力手段と、前記複数の画像を補正する際に用いられる補正用画像を作成する補正画像作成手段と、を備え、前記補正画像作成手段は、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得手段と、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出手段と、を有することを特徴とする画像処理装置。前記補正画像作成手段は、前記任意画像の位置調整を行い、前記共通部分における前記被写体を一致させる画像調整手段をさらに有し、前記差分取得手段は、前記画像調整手段により位置調整された状態で、前記差分を取得することを特徴とする請求項1に記載の画像処理装置。前記補正画像作成手段は、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定するシェーディング推定手段をさらに有することを特徴とする請求項1又は2に記載の画像処理装置。前記複数の画像は、撮像視野を平行移動させて撮像された画像であることを特徴とする請求項1~3のいずれか1項に記載の画像処理装置。前記複数の画像は、撮像倍率を変化させて撮像された画像であることを特徴とする請求項1~3のいずれか1項に記載の画像処理装置。前記複数の画像は、撮像視野を回転させて撮像された画像であることを特徴とする請求項1~3のいずれか1項に記載の画像処理装置。請求項1~6のいずれか1項に記載の画像処理装置と、前記被写体を撮像する撮像手段と、を備えることを特徴とする撮像装置。請求項1~6のいずれか1項に記載の画像処理装置と、前記被写体としての標本を載置可能なステージと、前記ステージに対向して設けられる光学系と、前記光学系を介して前記標本上に設定される視野を撮像して画像を取得する画像取得部と、前記ステージと前記光学系との内の少なくとも一方を、前記光学系の光軸と直交する方向に移動させることにより、前記撮像視野を変化させるステージ位置変更部と、を有する顕微鏡装置と、を備えることを特徴とする顕微鏡システム。被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得ステップと、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づいて前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、を含むことを特徴とする画像処理方法。被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得ステップと、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づいて前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、をコンピュータに実行させることを特徴とする画像処理プログラム。ところが、複数の画像を繋ぎ合わせた連結画像においては、個別の画像に生じたシェーディングの影響が規則的に現れる。このため、シェーディングが連結画像の画質に与える影響は特に大きい。従って、個別の画像に対して精度の良いシェーディング補正を行うことが重要となる。一般的に、顕微鏡画像のシェーディング補正は、ステージに被写体を配置しない状態で視野を撮像して得られたキャリブレーション画像を用いて行われる。特許文献1参照。しかしながら、同じ顕微鏡であっても、キャリブレーション画像の撮像時と、被写体を撮像する時とでは、シェーディングの状態が変化している場合がある。そこで、予め用意されたキャリブレーション画像ではなく、被写体が写った補正対象画像からシェーディング補正用のデータを作成し、このシェーディング補正用データを用いて当該補正対象画像をシェーディング補正することも提案されている特許文献2参照。しかしながら、補正対象画像には、シェーディングに由来する輝度成分の他に、被写体に由来する輝度成分等の様々な成分が含まれており、シェーディングに由来する輝度成分と、その他の成分とを区別することは困難である。そのため、補正対象画像から作成されたシェーディング補正用データを用いると、被写体に由来する輝度成分のように、本来補正の影響が及ぶべきでない成分に対しても補正が行われてしまう。そして、このようにシェーディング補正の精度が低下した結果、補正された画像にアーティファクトが生じるなどして画質がかえって低下してしまうおそれがある。本発明は、上記に鑑みてなされたものであって、補正対象画像の画像データを用いて、精度の高いシェーディング補正を行うことができる画像処理装置、撮像装置、顕微鏡システム、画像処理方法、及び画像処理プログラムを提供することを目的とする。上述した課題を解決し、目的を達成するために、本発明に係る画像処理装置は、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力手段と、前記複数の画像を補正する際に用いられる補正用画像を作成する補正画像作成手段と備え、前記補正画像作成手段は、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得手段と、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出手段とを有することを特徴とする。上記画像処理装置において、前記補正画像作成手段は、前記任意画像の位置調整を行い、前記共通部分における前記被写体を一致させる画像調整手段をさらに有し、前記差分取得手段は、前記画像調整手段により位置調整された状態で、前記差分を取得することを特徴とする。上記画像処理装置において、前記補正画像作成手段は、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定するシェーディング推定手段をさらに有することを特徴とする。上記画像処理装置において、前記複数の画像は、撮像視野を平行移動させて撮像された画像であることを特徴とする。上記画像処理装置において、前記複数の画像は、撮像倍率を変化させて撮像された画像であることを特徴とする。上記画像処理装置において、前記複数の画像は、撮像視野を回転させて撮像された画像であることを特徴とする。本発明に係る撮像装置は、上記画像処理装置と、前記被写体を撮像する撮像手段とを備えることを特徴とする。本発明に係る顕微鏡システムは、上記画像処理装置と、前記被写体としての標本を載置可能なステージと、前記ステージに対向して設けられる光学系と、前記光学系を介して前記標本上に設定される視野を撮像して画像を取得する画像取得部と、前記ステージと前記光学系との内の少なくとも一方を、前記光学系の光軸と直交する方向に移動させることにより、前記撮像視野を変化させるステージ位置変更部と、を有する顕微鏡装置とを備えることを特徴とする。本発明に係る画像処理方法は、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得ステップと、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づいて前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップとを含むことを特徴とする。本発明に係る画像処理プログラムは、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ2つの任意画像から、画素の輝度の差分を取得する差分取得ステップと、前記差分に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づいて前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップとをコンピュータに実行させることを特徴とする。本発明によれば、複数の画像間で被写体が共通する共通部分を持つ2つの任意画像から画素の輝度の差分を取得し、該差分に基づいて検出されたシェーディング成分に基づいて、上記複数の画像の補正に用いられる補正用画像を作成するので、精度の高いシェーディング補正を行うことが可能となる。以下、本発明に係る実施の形態について、図面を参照しながら詳細に説明する。なお、これらの実施の形態により本発明が限定されるものではない。また、以下の説明においては、本発明に係る画像処理装置を顕微鏡システムに適用した例を説明するが、本発明に係る画像処理装置は、ディジタルカメラ等、撮像機能を有する各種機器に適用することが可能である。実施の形態1図1は、本発明の実施の形態1に係る顕微鏡システムの構成を示す模式図である。図1に示すように、実施の形態1に係る顕微鏡システム1は、顕微鏡装置10と、該顕微鏡装置10の動作を制御すると共に、該顕微鏡装置10によって取得された画像を処理する画像処理装置11とを備える。図2は、顕微鏡装置10の構成を概略的に示す模式図である。図2に示すように、顕微鏡装置10は、略C字形のアーム100と、該アーム100に取り付けられ、標本SPが載置される標本ステージ101と、鏡筒103の一端側に三眼鏡筒ユニット106を介して標本ステージ101と対向するように設けられた対物レンズ102と、鏡筒103の他端側に設けられた画像取得部104と、標本ステージ101を移動させるステージ位置変更部105とを有する。三眼鏡筒ユニット106は、対物レンズ102から入射した標本の観察光を、画像取得部104と後述する接眼レンズユニット107に分岐する。接眼レンズユニット107は、ユーザが標本を直接観察するためのものである。なお、以下において、対物レンズ102の光軸L方向をZ軸方向とし、このZ軸方向と直交する平面をXY平面とする。図2において、顕微鏡装置10は、標本ステージ101の主面が概ねXY平面と一致するように設置されている。対物レンズ102は、倍率が互いに異なる複数の対物レンズ例えば、対物レンズ102’を保持可能なレボルバ108に取り付けられている。このレボルバ108を回転させて、標本ステージ101と対向する対物レンズ102、102’を変更することにより、画像取得部104が撮像する画像の倍率を変化させることができる。鏡筒103の内部には、複数のズームレンズと、これらのズームレンズの位置を変化させる駆動部いずれも図示せずとを含むズーム部が設けられている。ズーム部は、各ズームレンズの位置を調整することにより、撮像視野内の被写体を拡大又は縮小させる。なお、鏡筒103内の駆動部にエンコーダをさらに設けても良い。この場合、エンコーダの出力値を画像処理装置11に出力し、画像処理装置11において、エンコーダの出力値からズームレンズの位置を検出して撮像視野の倍率を自動算出するようにしても良い。画像取得部104は、例えばCCDやCMOS等の撮像素子を含み、該撮像素子が備える各画素においてR赤、G緑、B青の各バンドにおける画素レベル画素値を持つカラー画像を撮像可能なカメラである。画像取得部104は、対物レンズ102から鏡筒103内の光学系を介して入射した光観察光を受光し、観察光に対応する画像データを生成して画像処理装置11に出力する。ステージ位置変更部105は、例えば、モータ105aを含み、標本ステージ101の位置をXY平面内で移動させることにより、撮像視野を変化させる。また、ステージ位置変更部105には、標本ステージ101をZ軸に沿って移動させることにより、対物レンズ102の焦点を標本SPに合わせる。また、ステージ位置変更部105には、標本ステージ101の位置を検出して検出信号を画像処理装置11に出力する位置検出部105bが設けられている。位置検出部105bは、例えばモータ105aの回転量を検出するエンコーダによって構成される。或いは、ステージ位置変更部105を、画像処理装置11の制御部160後述の制御に従ってパルスを発生するパルス発生部とステッピングモータとによって構成しても良い。画像処理装置11は、当該画像処理装置11に対する指示や情報の入力を受け付ける入力部110と、画像取得部104から出力された画像の入力を受け付けるインタフェースである画像入力部120と、顕微鏡画像やその他の情報を表示する表示部130と、記憶部140と、顕微鏡装置10により取得された画像に対して所定の画像処理を施す演算部150と、これらの各部の動作を制御する制御部160とを備える。入力部110は、キーボード、各種ボタン、各種スイッチ等の入力デバイスや、マウスやタッチパネル等のポインティングデバイスを含み、これらのデバイスを介して入力された信号を受け付けて制御部160に入力する。表示部130は、例えば、LCDLiquidCrystalDisplayやELElectroLuminescenceディスプレイ、CRTCathodeRayTubeディスプレイ等の表示装置によって構成され、制御部160から出力された制御信号に従って、各種画面を表示する。記憶部140は、更新記録可能なフラッシュメモリ、RAM、ROM等の半導体メモリや、内蔵若しくはデータ通信端子で接続されたハードディスク、MO、CD-R、DVD-R等の記録媒体及び該記録媒体に記録された情報を読み取る読取装置等によって構成される。記憶部140は、画像取得部104から出力された画像データや、演算部150及び制御部160がそれぞれ実行する各種プログラムや各種設定情報を記憶する。具体的には、記憶部140は、画像取得部104が取得した画像に対してシェーディング補正を施す画像処理プログラム141を記憶する。演算部150は、例えばCPU等のハードウェアによって構成され、記憶部140に記憶された画像処理プログラム141を読み込むことにより、記憶部140に記憶された画像データに対応する画像に対してシェーディング補正を施す画像処理を実行する。より詳細には、演算部150は、画像をシェーディング補正するための補正用画像を作成する補正用画像作成部151と、補正用画像を用いて画像を補正する画像補正部156とを有する。このうち、補正用画像作成部151は、被写体の少なくとも一部が共通する複数の画像間で共通部分が一致するように画像の位置調整を行う画像調整部152と、位置調整がなされた複数の画像間で対応する画素の輝度の差分を取得する差分取得部153と、輝度の差分に基づいて当該画像内に生じているシェーディング成分を検出するシェーディング検出部154と、該シェーディング成分に基づき、上記共通部分以外の領域におけるシェーディングの影響を推定して補正用画像を作成するシェーディング推定部155とを有する。制御部160は、例えばCPU等のハードウェアによって構成され、記憶部140に記憶された各種プログラムを読み込むことにより、記憶部140に記憶された各種データや入力部110から入力される各種情報に基づき、画像処理装置11及び顕微鏡装置10の各部に指示やデータの転送を行い、顕微鏡システム1全体の動作を統括的に制御する。次に、本実施の形態1におけるシェーディング補正の原理について説明する。顕微鏡装置10によって撮像された画像には、被写体に対応する輝度成分以下、被写体成分というTx,yとシェーディングに対応する輝度成分以下、シェーディング成分というSx,yとが含まれる。なお、座標x,yは、画像内における各画素の位置座標を示す。そこで、各画素の輝度をIx,yとすると、輝度Ix,yは次式1-1によって表すことができる。Ix,y=Tx,y×Sx,y...1-1また、撮像条件の一部を変更して撮像した複数の画像間では、共通の被写体成分Tx,yに対するシェーディング成分Sx,yの影響が、画像内の位置座標に応じて変化する。そこで、複数の画像間で被写体成分Tx,yが共通する画素の輝度から、被写体成分Tx,yを打ち消すことにより、画像内の位置座標に応じたシェーディング成分Sx,yの変化を算出することができる。そして、このシェーディング成分の変化から、画像内におけるシェーディング成分Sx,yの分布を抽出することができる。さらに、次式1-2に示すように、抽出されたシェーディング成分Sx,yを輝度Ix,yから除去することにより、被写体成分Tx,yのみの画像を取得することができる。Tx,y=Ix,y/Sx,y...1-2なお、撮像条件を変更する方法としては、例えば、標本SPに対して撮像視野を平行移動させる、倍率を変更する、標本SPに対して撮像視野を回転させるといった方法がある。本実施の形態1においては、標本SPに対して撮像視野を平行移動させる方法により、撮像条件を変更する。次に、顕微鏡システム1の動作について説明する。図3は、顕微鏡システム1の動作を示すフローチャートである。まず、ステップS10において、顕微鏡装置10は、制御部160の制御の下で、標本SPに対して撮像視野を変化させながら、撮像視野の一部が重複する複数の画像を撮像する。より詳細には、顕微鏡装置10は、図4に示すように、標本SPが載置された標本ステージ101を平行移動させつつ、画像取得部104に撮像視野を撮像させる。その際、本実施の形態1において、顕微鏡装置10は、図5に示すように、撮像を行う毎に標本ステージ101をX方向に移動量dx、及び、Y方向に移動量dyだけ移動させ、隣接する撮像視野Vj、Vj1間で一部の領域C1が重なるようにして撮像を行うjは撮像順序を示す番号であり、j=1、2、...。なお、撮像視野は、X方向又はY方向のいずれか一方にのみ移動させても良い。画像取得部104は、このように撮像を行って取得した画像データを画像処理装置11に出力する。画像処理装置11は、画像取得部104から出力された画像データを、記憶部140に一旦記憶させる。ステップS11において、制御部160は、記憶部140から画像データを読み出し、該画像データに対応する複数の画像を演算部150に入力する。具体的には、図6a及びbに示すように、撮像視野Vjを写した画像Mjと、撮像視野Vj1を写した画像Mj1とが演算部150に入力される。画像Mj、Mj1間においては、画像Mj内の座標xi,yiの画素と、画像Mj1内の座標x’i,y’i=xi-dx,yi-dyの画素とが、互いに画素となる。互いに対応する画素には、同じ被写体が写っている。以下、画像Mj、Mj1の水平方向のサイズをwとし、垂直方向のサイズをhとする。続くステップS12において、画像調整部152は、位置検出部105bからの出力値に基づいて、複数の画像間で共通部分が一致するように画像の位置調整を行う。例えば、画像Mj、Mj1の場合、図7に示すように、画像Mjに対して画像Mj1を撮像視野の移動量dx及びdyの分だけずらし、領域C1同士が重なるように位置調整がなされる。ステップS13において、差分取得部153は、各画像の共通部分に含まれる画素の画素値から輝度を算出し、複数の画像間で互いに対応する画素の輝度の差分を算出する。例えば、画像Mj、Mj1の場合、画像Mj内の座標xi,yiにおける画素の輝度Ij、画像Mj1内の座標x’i,y’iにおける画素の輝度Ij1との差分が算出される。ここで、輝度Ij、Ij1はそれぞれ、次式1-3、1-4により与えられる。Ijxi,yi=Txi,yi×Sxi,yi...1-3Ij1x’i,y’i=Tx’i,y’i×Sx’i,y’i...1-4また、対応する画素間において被写体成分は等しいから、次式1-5が成り立つ。Txi,yi=Tx’i,y’i...1-5さらに、座標x’i,y’i=xi-dx,yi-dyを用いると、式1-4は次式1-6のように書き換えることができる。Ij1x’i,y’i=Txi,yi×Sxi-dx,yi-dy...1-6式1-3及び1-6より、次式1-7に示す関係が成り立つ。Ij1x’i,y’i/Ijxi,yi=Sxi-dx,yi-dy/Sxi,yi...1-7即ち、輝度の差分Ij1x’i,y’i/Ijxi,yiは、シェーディング成分の変化に相当する。ステップS13においては、このような差分Ij1x’i,y’i/Ijxi,yiが、共通部分である領域C1内の全ての座標に関して算出される。ステップS14において、シェーディング検出部154は、画像Mj、Mj1内のシェーディングの影響を近似するシェーディングモデルを作成し、領域C1内の各座標xi,yiにおける差分Ij1x’i,y’i/Ijxi,yiに基づいて、シェーディングモデルを修正する。各画像に生じるシェーディングの影響は、原理的には式1-7に従うはずであるが、実際には、対応する画素同士においても微妙に輝度が異なったり、シェーディングにばらつきがあったりするので、領域C1内の全ての座標において式1-7が成り立つわけではない。そこで、シェーディング成分Sを表すモデルを設定し、次式1-8に示す誤差評価関数により、評価値K1が最小となるときのシェーディング成分Sを求める。具体的なシェーディングモデルの一例として、本実施の形態1においては、次式1-9に示すように、画像の中心座標w/2,h/2を通る2次曲面を表す関数を用いる。ここで、2次曲面を用いるのは、シェーディング成分は一般に、画像の中心付近で少なく、画像の中心から離れるに従って増加するからである。Sx,y;a=1-a{x-w/22y-h/22}...1-9従って、式1-9を用いた誤差評価関数1-10によって与えられる評価値K1’が最小となるときの2次の係数パラメータaを求めることにより式1-10’、シェーディング成分Sを取得することができる。シェーディング検出部154は、このパラメータaを、図8に示す演算処理により求める。図8は、シェーディング検出部154の動作を示すフローチャートである。まず、ステップS151において、シェーディング検出部154は、パラメータaを初期化する。続くステップS152において、シェーディング検出部154は、パラメータa及びステップS13において算出された領域C1内の全ての座標xi,yiにおける輝度の差分Ij1x’i,y’i/Ijxi,yiを誤差評価関数1-10に代入して、評価値K1’を算出する。ステップS153において、シェーディング検出部154は、評価値K1’が所定の閾値未満であるか否かを判定する。閾値の大きさは、後に続く繰り返し処理の各回において代入されるパラメータaに基づいて補正画像を作成した場合に、パラメータaが互いに異なる補正画像間で差異が明確に認識されない程度に、充分小さい値を経験的に設定するものとする。この閾値は、繰り返し処理を終了させるために設定するものであるので、繰り返し処理を終了する条件として、繰り返しにおける評価値の変化が充分に小さくなったときに終了するなど、他の終了条件を設定しても良い。評価値K1’が所定の閾値以上である場合ステップS153:No、シェーディング検出部154は、パラメータaを修正するステップS154。その後、動作はステップS152に戻る。一方、評価値K1’が所定の閾値未満である場合ステップS153:Yes、シェーディング検出部154は、そのときのパラメータaを、式1-9におけるパラメータとして決定するステップS155。この決定されたパラメータaを含む式1-9が、領域C1における修正済みのシェーディングモデルを表す式である。その後、動作はメインルーチンに戻る。ステップS15において、シェーディング推定部155は、修正されたシェーディングモデルの適用範囲を、領域C1以外の画像Mj、Mj1内の領域に拡張し、画像内の全領域におけるシェーディングを補正するための補正用画像を作成する。補正用画像は、シェーディングモデルSx,y;aを各画素の輝度とする画像である。ステップS16において、画像補正部156は、ステップS15において作成された補正用画像を用いて、画像の補正を行う。即ち、補正対象である画像内の各画素の輝度Ix,yを取得し、次式1-11により、被写体成分Tx,yを算出する。Tx,y=Ix,y/Sx,y;a...1-11それにより、シェーディングの影響が除去された補正済みの画像が取得される。なお、補正対象の画像は、シェーディング補正用画像の作成に用いた画像Mj、Mj1に限らず、顕微鏡装置10において撮像した他の画像であっても良い。さらに、ステップS17において演算部150は、補正済みの画像を出力する。それに応じて、制御部160は、補正済みの画像を表示部130に表示させると共に、補正済みの画像に対応する画像データを記憶部140に記憶させる。以上説明したように、実施の形態1によれば、補正対象である画像自体の輝度に基づいてシェーディング補正用画像を作成するので、顕微鏡装置10において経時的なシェーディングの変化が生じている場合であっても、当該画像に対して精度の高いシェーディング補正を施すことができる。また、実施の形態1によれば、画像内の輝度を被写体成分とシェーディング成分とを分離し、シェーディング成分のみを抽出してシェーディング補正用画像を作成するので、シェーディング補正の精度を向上させることが可能となる。実施の形態2次に、本発明の実施の形態2について説明する。実施の形態2に係る顕微鏡システムの全体的な構成及び動作は、実施の形態1と共通であり、実施の形態2においては、撮像視野の倍率を変化させて複数の画像を撮像し、これらの画像を用いてシェーディング補正用画像を作成することを特徴とする。従って、以下においては、互いに倍率が異なる複数の画像を用いてシェーディング補正用画像を作成する処理についてのみ説明する。まず、顕微鏡システム1は、図9に示すように、標本SPの視野領域Vjを撮像して画像Mjを取得する。次いで、ズームを調整して倍率を変更し、同じ視野領域Vjを撮像して画像Mj1を取得する。この場合、画像Mj、Mj1間では、画像全体が、同じ被写体が写った共通部分となる。演算部150において、画像調整部152は、画像Mjの中心Oと、Mj1の中心O’とが一致するように、画像Mj、Mj1の位置調整を行う。図10に示すように、画像Mj、Mj1間では、画像Mjの中心Oから距離riだけ離れ、所定の軸からの回転角φiに位置する画素Piと、画像Mj2の中心O’から距離r’iだけ離れ、所定の軸からの回転角φiに位置する画素P’iとが互いに対応する画素となる。続いて、差分取得部153は、画像Mj、Mj1間で対応する画素の輝度の差分を算出する。ここで、画像Mjの倍率をmj、画像Mj1の倍率をmj1とすると、距離r’iは次式2-1によって表すことができる。r’i=mj1/mj×ri...2-1上述したとおり、各画素の輝度Iは、被写体成分Tとシェーディング成分Sとで構成されるから、画素Piの輝度Ijri,φi及び画素P’iの輝度Ij1r’i,φiは、それぞれ、次式2-2及び2-3によって表すことができる。Ijri,φi=Tri,φi×Sri,φi...2-2Ij1r’i,φi=Tj1r’i,φi×Sj1r’i,φi...2-3ここで、一般には、撮像倍率を変更するとシェーディングも変化する。しかしながら、倍率の変化が小さい場合には、シェーディングの変化は極めて小さく、無視することができる。それに対して、画像内の被写体像の変化は倍率の変化と厳密に一致する。即ち、次式2-4のとおりである。Tj1r’i,φi=Tri,φi...2-4従って、式2-2~2-4より、輝度の差分Ij1r’i/Ijriは次式2-5によって与えられる。Ij1r’i,φi/Ijri,φi=Sj1r’i,φi/Sri,φi...2-5式2-5に示すとおり、輝度の差分Ij1r’i,φi/Ijri,φiは、シェーディング成分の変化に相当する。ここで、一般にシェーディング成分は、光学系の光軸Lからの距離に応じて軸対称に変化するから、式2-5は次式2-6に書き換えることができる。Ij1r’i,φi/Ijri,φi=Sj1r’i/Sri...2-6差分取得部153は、このような差分Ij1r’i,φi/Ijri,φiを、画像Mj内の全ての画素に対して算出する。続いて、シェーディング検出部154は、画像Mj、Mj1内のシェーディングの影響を近似するシェーディングモデルを作成し、各座標における差分Ij1r’i,φi/Ijri,φiに基づいて、シェーディングモデルを修正する。本実施の形態2の場合、シェーディングモデルを修正するための誤差評価関数は次式2-7によって与えられ、評価値K2を最小にするシェーディングモデルSriを求める演算が行われる。また、具体的なシェーディングモデルの一例として、本実施の形態2においては、次式2-8に示すように、画像Mjの中心からの距離rに依存して変化する2次曲面を表す関数を用いる。Sr;b=1-b×r2...2-8従って、式2-7を用いた誤差評価関数2-9によって与えられる評価値K2’が最小となるときの2次の係数であるパラメータbを求めることにより式2-9’、シェーディング成分Sを取得することができる。シェーディング検出部154は、このパラメータbを、図8に示す演算処理により求める。ただし、図8において、パラメータaはパラメータbに読み替えられ、評価値K1’は評価値K2’に読み替えられる。それによって決定されたパラメータbを含む式2-8が、修正済みのシェーディングモデルを表す式である。続いて、シェーディング推定部155は、修正されたシェーディングモデルの適用範囲を画像Mj1全体即ち、画像Mjよりも若干広い範囲に拡張し、画像Mj及び画像Mj1内の全領域におけるシェーディングを補正するための補正用画像を作成する。その後の処理については、実施の形態1と同様である。以上説明したように、実施の形態2によれば、画像Mj全体の画素の輝度に基づいて補正用画像を作成するので、精度の高いシェーディング補正を行うことができる。実施の形態3次に、本発明の実施の形態3について説明する。実施の形態3に係る顕微鏡システムの全体的な構成及び動作は、実施の形態1と共通であり、実施の形態3においては、撮像視野を回転させて複数の画像を撮像し、これらの画像を用いてシェーディング補正用画像を作成することを特徴とする。従って、以下においては、XY平面における座標軸が互いに交差する複数の画像を用いてシェーディング補正用画像を作成する処理についてのみ説明する。まず、顕微鏡システム1は、図11及び図12に示すように、標本SPの視野領域Vjを撮像して画像Mjを取得する。次いで、標本ステージ101を所定の回転中心点に対して角度dθだけXY平面内で回転させ、視野領域Vj1を撮像して画像Mj1を取得する。なお、実施の形態3においては、回転中心点を撮像領域Vjの中心Oに設定する。この場合、画像Mj、Mj1間では、図11に示す領域C3が、同じ被写体が写った共通部分となる。また、図12a及びbに示すように、画像Mjの回転中心点中心Oから距離riだけ離れ、所定の軸からの角度θiに位置する画素Qiと、画像Mj1の回転中心点中心Oから距離riだけ離れ、所定の軸からの角度θ’i=θi-dθに位置する画素Q’iとが、互いに対応する画素となる。演算部150において、画像調整部152は、図13に示すように、画像Mjに対して画像Mj1を角度dθだけ回転させ、画像Mj、Mj1間で領域C3が一致するよう、位置調整を行う。続いて、差分取得部153は、画像Mj、Mj1間で対応する画素の輝度の差分を算出する。具体的には、画素Qiの輝度Ijri,θiと、画素Q’iの輝度Ij1ri,θ’iとの差分Ij1ri,θ’i/Ijri,θiを算出する。上述したとおり、各画素の輝度Iは、被写体成分Tとシェーディング成分Sとで構成されるから、画素Qiの輝度Ijri,θi及び画素Q’iの輝度Ij1ri,θ’iは、それぞれ、次式3-1及び3-2によって表すことができる。Ijri,θi=Tri,θi×Sri,θi...3-1Ij1ri,θ’i=Tj1ri,θ’i×Sj1ri,θ’i...3-2また、対応する画素間において被写体成分は等しいから、次式3-3が成り立つ。Tri,θi=Tj1ri,θ’i...3-3従って、式3-1~3-3より、差分Ij1ri,θ’i/Ijri,θiは次式3-4によって与えられる。Ij1ri,θ’i/Ijri,θi=Sj1ri,θ’i/Sjri,θi...3-4即ち、輝度の差分Ij1ri,θ’i/Ijri,θiは、シェーディング成分の変化に相当する。ここで、シェーディング成分Sには、レンズに起因して生じるシェーディング成分以下、シェーディング成分Slとすると、照明等、レンズ以外の要因によって生じるシェーディング成分以下、シェーディング成分Smとするとが含まれる。即ち、シェーディング成分Sは、次式3-5によって与えられる。Sri,θi=Smri,θi×Slri,θi...3-5この内、シェーディング成分Slは光軸Lを中心とする軸対称に生じるため、次式3-6に示すように、画像の中心から画素までの距離riのみによってモデル化することができる。Sri,θi=Smri,θi×Slri...3-6式3-5及び3-6より、式3-4を次式3-7のように書き換えることができる。式3-7より、対応する画素間の差分Ij1ri,θ’i/Ijri,θiは、照明等に起因するシェーディング成分Smの変化を表すと言える。差分取得部153は、このような差分Ij1ri,θ’i/Ijri,θiを、領域C3内の全ての画素に対して算出する。続いて、シェーディング検出部154は、画像Mj、Mj1内のシェーディングの影響を近似するシェーディングモデルを作成し、各座標における差分Ij1ri,θ’i/Ijri,θiに基づいて、シェーディングモデルを修正する。本実施の形態3の場合、シェーディングモデルを修正するための誤差評価関数は次式3-8によって与えられ、評価値K3を最小にするシェーディングモデルSmri,θiを求める演算が行われる。具体的なシェーディングモデルの一例として、本実施の形態3においては、次式3-9に示すように、画像Mjの中心からの距離r及び角度θに依存して変化する2次曲面を表す関数を用いる。式3-9において、r0及びθ0は所定の定数である。従って、式3-9を用いた誤差評価関数3-10によって与えられる評価値K3’が最小となるときの2次の係数であるパラメータcを求めることにより式3-10’、シェーディング成分Sを取得することができる。シェーディング検出部154は、このパラメータcを、図8に示す演算処理により求める。ただし、図8において、パラメータaはパラメータcに読み替えられ、評価値K1’は評価値K3’に読み替えられる。それによって決定されたパラメータcを含む式3-9が、照明等に起因するシェーディング成分Smの修正済みモデルを表す式である。続いて、シェーディング推定部155は、画像Mj及び画像Mj1内の全領域におけるシェーディングを補正するための補正用画像を作成する。ここで、レンズに起因するシェーディング成分Slは経時的な変化が少ないため、コサイン四乗則により正確に予測することができる。コサイン四乗則は、レンズに照度I0の光が入射するとき、入射光の光軸Lに対する角度Θと、入射後の光の照度I’との関係を表す。I’=I0cos4Θ...3-11シェーディング推定部155は、式3-11に基づいて、レンズに起因するシェーディング成分Slを表すシェーディングモデルSlrを作成する。シェーディングモデルSlとしては、例えば、式3-11に基づいて、光軸中心からの距離rの値に応じたシェーディング量をルックアップテーブル化したものを用いることができる。記憶部140は、レンズごとに作成された複数種類のルックアップテーブルを格納しており、シェーディング推定部155は、選択されたレンズに対応するルックアップテーブルを記憶部140から読み出し、このルックアップテーブルを使用して、シェーディング成分Slを求める。さらに、シェーディング推定部155は、式3-9によって与えられるシェーディングモデルを画像全体領域C1以外の画像Mj及び画像Mj1の領域に拡張すると共に、上記ルックアップテーブルに基づくシェーディングモデルSlrを用いて、次式3-12に示すように、画像Mj及び画像Mj1内の全領域におけるトータルのシェーディングを補正するための補正用画像を作成する。STOTAL=Slr×Smr,θ...3-12その後の処理については、実施の形態1と同様である。以上説明したように、実施の形態3によれば、照明等に起因する経時的な変化が比較的大きいシェーディングについても精度良く補正を行うことができる。従って、多量の画像を撮像した場合や、撮像時間が長引いたりした場合においても、各画像に対して正確なシェーディング補正を施すことが可能となる。実施の形態4次に、本発明の実施の形態4について説明する。図14は、実施の形態4に係る顕微鏡システムの構成を示すブロック図である。図14に示すように、実施の形態4に係る顕微鏡システム2は、図1に示す画像処理装置11の代わりに画像処理装置20を備える。画像処理装置20は、図1に示す演算部150に対してバーチャルスライドVS画像作成部201をさらに有する演算部200を備える。VS画像作成部201以外の画像処理装置20及び顕微鏡システム2の構成については、実施の形態1において説明したものと同様である。VS画像作成部201は、顕微鏡装置10が標本SPに対する撮像視野を平行移動させながら撮像した複数の画像を繋ぎ合わせて、標本SP全体に対応する画像VS画像を作成する。次に、顕微鏡システム2の動作について説明する。図15は、顕微鏡システム2の動作を示すフローチャートである。また、図16は、VS画像作成部が作成したバーチャルスライド画像を示す模式図である。まず、ステップS10において、顕微鏡装置10は、制御部160の制御の下で、標本SPに対する撮像視野を平行移動させながら、撮像視野の一部が重複する複数の画像を撮像する。続くステップS11及びS12の動作は、実施の形態1と同様である。ステップS12に続くステップS20において、VS画像作成部201は、図16に示すように、画像調整部152による位置調整済みの複数の画像を繋ぎ合わせて、バーチャルスライド画像VSを作成する。このバーチャルスライド画像VSにおいて、X方向で隣接する画像MとM間では、領域Ckが互いに重複する即ち、被写体が一致する共通部分であり、Y方向で隣接する画像MとM間では、領域Clが互いに重複する共通部分であるk,l=1、2、...。続くステップS13~S15の動作は、実施の形態1と同様である。ただし、この際の演算対象は、バーチャルスライド画像VS内の全ての共通部分領域Ck及び領域Clであっても良いし、一部の共通部分であっても良い。後者の場合、具体的には、全ての領域Ckのみ、全ての領域Clのみ、最初と最後の領域Ck若しくは領域Clのみ、バーチャルスライド画像VSの中央近傍の領域Ck若しくは領域Clのみ、所定間隔で抽出した領域Ck若しくは領域Cl、又は、ランダムに選択した領域Ck若しくは領域Clのいずれであっても良い。複数の共通部分が演算対象となる場合には、それらの共通部分に対応する複数の補正用画像が作成される。ステップS15に続くステップS30において、画像補正部156は、ステップS15において作成された補正用画像を用いて、バーチャルスライド画像VSを構成する各画像Mを補正する。この際、ステップS15において複数の補正用画像が作成された場合には、補正用画像の基となった共通部分の位置に応じて補正用画像を決定すると良い。例えば、全ての領域Ckから補正用画像が作成された場合には、X方向において隣接する画像M、Mを、両画像の共通部分である領域Ckに基づく補正用画像を用いて補正する。また、全ての領域Clから補正用画像が作成された場合には、Y方向において隣接する画像M、Mを、両画像の共通部分である領域Clに基づく補正用画像を用いて補正する。或いは、所定間隔で抽出された領域Ck若しくは領域Cl、又は、ランダムに選択された領域Ck若しくは領域Clから補正用画像が作成された場合には、領域Ck若しくは領域Clから所定範囲内の画像Mを、当該領域Ck若しくは領域Clに基づく補正用画像を用いて補正することとしても良い。続くステップS31において、演算部200は、補正済みの画像を繋ぎ合わせたバーチャルスライド画像VSを出力する。それに応じて、制御部160は、バーチャルスライド画像VSを表示部130に表示させると共に、バーチャルスライド画像VSに対応する画像データを記憶部140に記憶させる。以上説明したように、実施の形態4によれば、個別の撮像視野に生じたシェーディングが補正された画像によってバーチャルスライド画像を構成するので、画質の良いバーチャルスライド画像を取得することが可能となる。変形例上記実施の形態1~3においては、顕微鏡装置10によって取得された画像に対する処理を説明したが、画像処理装置11は、顕微鏡装置10以外の各種撮像機器によって取得された画像を処理することも可能である。例えば、画像処理装置11を、パノラマ画像の撮像が可能なディジタルカメラに適用しても良い。この場合、視野の端部が重複する複数の画像を撮像し、視野の重複部分に基づいて作成した補正用画像を用いて各画像を補正して、それらの画像を繋ぎ合わせることにより、画質の良いパノラマ画像を生成することが可能である。本発明は、上述した各実施の形態1~4及び変形例そのままに限定されるものではなく、各実施の形態1~4に開示されている複数の構成要素を適宜組み合わせることによって、種々の発明を形成することができる。例えば、実施の形態1~4に示される全構成要素からいくつかの構成要素を除外して形成してもよい。あるいは、異なる実施の形態に示した構成要素を適宜組み合わせて形成してもよい。"