JP5911296B2 - 画像処理装置、撮像装置、顕微鏡システム、画像処理方法、及び画像処理プログラム - Google Patents

Description

本発明は、被写体を撮像した画像に対して画像処理を施す画像処理装置、撮像装置、顕微鏡システム、画像処理方法、及び画像処理プログラムに関する。

カメラや顕微鏡等の光学機器においては、通常、レンズ等の結像光学系の特性により、光軸と直交する面内の中心に比べて周辺領域の光量が低下する現象が生じる。この現象は、一般に、シェーディングと呼ばれる。このため、従来、撮像後の画像に対して経験的に得られた補正値や実測値等に基づく画像処理を施すことにより、画像の劣化を抑制していた。このような画像処理は、シェーディング補正と呼ばれる（例えば、特許文献１及び２を参照）。

ところで、顕微鏡は、高倍率且つ高分解能の被写体（標本）観察が可能である反面、倍率を高くするほど、一度に観察可能な視野が狭くなるという問題がある。このため、１つの標本に対して視野をずらしながら撮像を行って複数の画像を取得し、これらの画像を繋ぎ合わせることにより、標本全体に対応するサイズに視野が拡大された画像を合成する画像処理が行われることがある。このような視野拡大処理がなされた連結画像は、バーチャルスライド画像と呼ばれ、バーチャルスライド画像の取得が可能な顕微鏡システムは、バーチャルスライドシステム又はバーチャル顕微鏡システムと呼ばれる（例えば、特許文献３及び４を参照）。

特開２００９−１５９０９３号公報 特開２００４−２７２０７７号公報 特開２００８−１９１４２７号公報 特開２０１１−１４１３９１号公報

ところが、複数の画像を繋ぎ合わせた連結画像においては、個別の画像に生じたシェーディングの影響が規則的に現れる。このため、シェーディングが連結画像の画質に与える影響は特に大きい。従って、個別の画像に対して精度の良いシェーディング補正を行うことが重要となる。

一般的に、顕微鏡画像のシェーディング補正は、ステージに被写体を配置しない状態で視野を撮像して得られたキャリブレーション画像を用いて行われる。（特許文献１参照）。しかしながら、同じ顕微鏡であっても、キャリブレーション画像の撮像時と、被写体を撮像する時とでは、シェーディングの状態が変化している場合がある。そこで、予め用意されたキャリブレーション画像ではなく、被写体が写った補正対象画像からシェーディング補正用のデータを作成し、このシェーディング補正用データを用いて当該補正対象画像をシェーディング補正することも提案されている（特許文献２参照）。

しかしながら、補正対象画像には、シェーディングに由来する輝度成分の他に、被写体に由来する輝度成分等の様々な成分が含まれており、シェーディングに由来する輝度成分と、その他の成分とを区別することは困難である。そのため、補正対象画像から作成されたシェーディング補正用データを用いると、被写体に由来する輝度成分のように、本来補正の影響が及ぶべきでない成分に対しても補正が行われてしまう。そして、このようにシェーディング補正の精度が低下した結果、補正された画像にアーティファクトが生じるなどして画質がかえって低下してしまうおそれがある。

本発明は、上記に鑑みてなされたものであって、補正対象画像の画像データを用いて、精度の高いシェーディング補正を行うことができる画像処理装置、撮像装置、顕微鏡システム、画像処理方法、及び画像処理プログラムを提供することを目的とする。

上述した課題を解決し、目的を達成するために、本発明に係る画像処理装置は、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力手段と、前記複数の画像を補正する際に用いられる補正用画像を作成する補正画像作成手段と、を備え、前記補正画像作成手段は、前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得手段と、前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出手段と、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定するシェーディング推定手段と、を有することを特徴とする。

上記画像処理装置において、前記補正画像作成手段は、前記任意画像の位置調整を行い、前記共通部分における前記被写体を一致させる画像調整手段をさらに有し、前記比率取得手段は、前記画像調整手段により位置調整された状態で、前記比を取得することを特徴とする。

上記画像処理装置において、前記補正画像作成手段は、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定するシェーディング推定手段をさらに有することを特徴とする。

上記画像処理装置において、前記複数の画像は、撮像視野を平行移動させて撮像された画像であることを特徴とする。
上記画像処理装置において、前記複数の画像は、撮像倍率を変化させて撮像された画像であることを特徴とする。
上記画像処理装置において、前記複数の画像は、撮像視野を回転させて撮像された画像であることを特徴とする。

本発明に係る撮像装置は、上記画像処理装置と、前記被写体を撮像する撮像手段とを備えることを特徴とする。

本発明に係る顕微鏡システムは、上記画像処理装置と、前記被写体としての標本を載置可能なステージと、前記ステージに対向して設けられる光学系と、前記光学系を介して前記標本上に設定される視野を撮像して画像を取得する画像取得部と、前記ステージと前記光学系との内の少なくとも一方を、前記光学系の光軸と直交する方向に移動させることにより、前記撮像視野を変化させるステージ位置変更部と、を有する顕微鏡装置とを備えることを特徴とする。

本発明に係る画像処理方法は、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得ステップと、前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定して前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、を含むことを特徴とする。

本発明に係る画像処理プログラムは、被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得ステップと、前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定して前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、をコンピュータに実行させることを特徴とする。

本発明によれば、複数の画像間で被写体が共通する共通部分を持つ２つの任意画像から画素の輝度の比を取得し、該比に基づいて検出されたシェーディング成分に基づいて、上記複数の画像の補正に用いられる補正用画像を作成するので、精度の高いシェーディング補正を行うことが可能となる。

図１は、本発明の実施の形態１に係る顕微鏡システムの構成例を示す模式図である。 図２は、図１に示す顕微鏡装置の構成を概略的に示す模式図である。 図３は、図１に示す顕微鏡システムの動作を示すフローチャートである。 図４は、実施の形態１における画像の撮像方法を説明する模式図である。 図５は、実施の形態１における画像の撮像方法を説明する模式図である。 図６は、一部が共通する複数の画像を示す模式図である。 図７は、位置調整された複数の画像を示す模式図である。 図８は、図１に示すシェーディング検出部の動作を示すフローチャートである。 図９は、本発明の実施の形態２における画像の撮像方法を説明する模式図である。 図１０は、複数の画像間で対応する画素を示す模式図である。 図１１は、本発明の実施の形態３における画像の撮像方法を説明する模式図である。 図１２は、一部が共通する複数の画像を示す模式図である。 図１３は、位置調整された複数の画像を示す模式図である。 図１４は、本発明の実施の形態４に係る顕微鏡システムの構成例を示すブロック図である。 図１５は、図１４に示す顕微鏡システムの動作を示すフローチャートである。 図１６は、図１４に示すＶＳ画像作成部が作成したバーチャルスライド画像を示す模式図である。

以下、本発明に係る実施の形態について、図面を参照しながら詳細に説明する。なお、これらの実施の形態により本発明が限定されるものではない。また、以下の説明においては、本発明に係る画像処理装置を顕微鏡システムに適用した例を説明するが、本発明に係る画像処理装置は、ディジタルカメラ等、撮像機能を有する各種機器に適用することが可能である。

（実施の形態１）
図１は、本発明の実施の形態１に係る顕微鏡システムの構成を示す模式図である。図１に示すように、実施の形態１に係る顕微鏡システム１は、顕微鏡装置１０と、該顕微鏡装置１０の動作を制御すると共に、該顕微鏡装置１０によって取得された画像を処理する画像処理装置１１とを備える。

図２は、顕微鏡装置１０の構成を概略的に示す模式図である。図２に示すように、顕微鏡装置１０は、略Ｃ字形のアーム１００と、該アーム１００に取り付けられ、標本ＳＰが載置される標本ステージ１０１と、鏡筒１０３の一端側に三眼鏡筒ユニット１０６を介して標本ステージ１０１と対向するように設けられた対物レンズ１０２と、鏡筒１０３の他端側に設けられた画像取得部１０４と、標本ステージ１０１を移動させるステージ位置変更部１０５とを有する。三眼鏡筒ユニット１０６は、対物レンズ１０２から入射した標本の観察光を、画像取得部１０４と後述する接眼レンズユニット１０７に分岐する。接眼レンズユニット１０７は、ユーザが標本を直接観察するためのものである。なお、以下において、対物レンズ１０２の光軸Ｌ方向をＺ軸方向とし、このＺ軸方向と直交する平面をＸＹ平面とする。図２において、顕微鏡装置１０は、標本ステージ１０１の主面が概ねＸＹ平面と一致するように設置されている。

対物レンズ１０２は、倍率が互いに異なる複数の対物レンズ（例えば、対物レンズ１０２’）を保持可能なレボルバ１０８に取り付けられている。このレボルバ１０８を回転させて、標本ステージ１０１と対向する対物レンズ１０２、１０２’を変更することにより、画像取得部１０４が撮像する画像の倍率を変化させることができる。

鏡筒１０３の内部には、複数のズームレンズと、これらのズームレンズの位置を変化させる駆動部（いずれも図示せず）とを含むズーム部が設けられている。ズーム部は、各ズームレンズの位置を調整することにより、撮像視野内の被写体を拡大又は縮小させる。なお、鏡筒１０３内の駆動部にエンコーダをさらに設けても良い。この場合、エンコーダの出力値を画像処理装置１１に出力し、画像処理装置１１において、エンコーダの出力値からズームレンズの位置を検出して撮像視野の倍率を自動算出するようにしても良い。

画像取得部１０４は、例えばＣＣＤやＣＭＯＳ等の撮像素子を含み、該撮像素子が備える各画素においてＲ（赤）、Ｇ（緑）、Ｂ（青）の各バンドにおける画素レベル（画素値）を持つカラー画像を撮像可能なカメラである。画像取得部１０４は、対物レンズ１０２から鏡筒１０３内の光学系を介して入射した光（観察光）を受光し、観察光に対応する画像データを生成して画像処理装置１１に出力する。

ステージ位置変更部１０５は、例えば、モータ１０５ａを含み、標本ステージ１０１の位置をＸＹ平面内で移動させることにより、撮像視野を変化させる。また、ステージ位置変更部１０５には、標本ステージ１０１をＺ軸に沿って移動させることにより、対物レンズ１０２の焦点を標本ＳＰに合わせる。

また、ステージ位置変更部１０５には、標本ステージ１０１の位置を検出して検出信号を画像処理装置１１に出力する位置検出部１０５ｂが設けられている。位置検出部１０５ｂは、例えばモータ１０５ａの回転量を検出するエンコーダによって構成される。
或いは、ステージ位置変更部１０５を、画像処理装置１１の制御部１６０（後述）の制御に従ってパルスを発生するパルス発生部とステッピングモータとによって構成しても良い。

画像処理装置１１は、当該画像処理装置１１に対する指示や情報の入力を受け付ける入力部１１０と、画像取得部１０４から出力された画像の入力を受け付けるインタフェースである画像入力部１２０と、顕微鏡画像やその他の情報を表示する表示部１３０と、記憶部１４０と、顕微鏡装置１０により取得された画像に対して所定の画像処理を施す演算部１５０と、これらの各部の動作を制御する制御部１６０とを備える。

入力部１１０は、キーボード、各種ボタン、各種スイッチ等の入力デバイスや、マウスやタッチパネル等のポインティングデバイスを含み、これらのデバイスを介して入力された信号を受け付けて制御部１６０に入力する。

表示部１３０は、例えば、ＬＣＤ（Liquid Crystal Display）やＥＬ（Electro Luminescence）ディスプレイ、ＣＲＴ（Cathode Ray Tube）ディスプレイ等の表示装置によって構成され、制御部１６０から出力された制御信号に従って、各種画面を表示する。

記憶部１４０は、更新記録可能なフラッシュメモリ、ＲＡＭ、ＲＯＭ等の半導体メモリや、内蔵若しくはデータ通信端子で接続されたハードディスク、ＭＯ、ＣＤ−Ｒ、ＤＶＤ−Ｒ等の記録媒体及び該記録媒体に記録された情報を読み取る読取装置等によって構成される。記憶部１４０は、画像取得部１０４から出力された画像データや、演算部１５０及び制御部１６０がそれぞれ実行する各種プログラムや各種設定情報を記憶する。具体的には、記憶部１４０は、画像取得部１０４が取得した画像に対してシェーディング補正を施す画像処理プログラム１４１を記憶する。

演算部１５０は、例えばＣＰＵ等のハードウェアによって構成され、記憶部１４０に記憶された画像処理プログラム１４１を読み込むことにより、記憶部１４０に記憶された画像データに対応する画像に対してシェーディング補正を施す画像処理を実行する。

より詳細には、演算部１５０は、画像をシェーディング補正するための補正用画像を作成する補正用画像作成部１５１と、補正用画像を用いて画像を補正する画像補正部１５６とを有する。このうち、補正用画像作成部１５１は、被写体の少なくとも一部が共通する複数の画像間で共通部分が一致するように画像の位置調整を行う画像調整部１５２と、位置調整がなされた複数の画像間で対応する画素の輝度の比を取得する比率取得部１５３と、輝度の比に基づいて当該画像内に生じているシェーディング成分を検出するシェーディング検出部１５４と、該シェーディング成分に基づき、上記共通部分以外の領域におけるシェーディングの影響を推定して補正用画像を作成するシェーディング推定部１５５とを有する。

制御部１６０は、例えばＣＰＵ等のハードウェアによって構成され、記憶部１４０に記憶された各種プログラムを読み込むことにより、記憶部１４０に記憶された各種データや入力部１１０から入力される各種情報に基づき、画像処理装置１１及び顕微鏡装置１０の各部に指示やデータの転送を行い、顕微鏡システム１全体の動作を統括的に制御する。

次に、本実施の形態１におけるシェーディング補正の原理について説明する。
顕微鏡装置１０によって撮像された画像には、被写体に対応する輝度成分（以下、被写体成分という）Ｔ（ｘ，ｙ）とシェーディングに対応する輝度成分（以下、シェーディング成分という）Ｓ（ｘ，ｙ）とが含まれる。なお、座標（ｘ，ｙ）は、画像内における各画素の位置座標を示す。そこで、各画素の輝度をＩ（ｘ，ｙ）とすると、輝度Ｉ（ｘ，ｙ）は次式（１−１）によって表すことができる。
Ｉ（ｘ，ｙ）＝Ｔ（ｘ，ｙ）×Ｓ（ｘ，ｙ） …（１−１）

また、撮像条件の一部を変更して撮像した複数の画像間では、共通の被写体成分Ｔ（ｘ，ｙ）に対するシェーディング成分Ｓ（ｘ，ｙ）の影響が、画像内の位置座標に応じて変化する。そこで、複数の画像間で被写体成分Ｔ（ｘ，ｙ）が共通する画素の輝度から、被写体成分Ｔ（ｘ，ｙ）を打ち消すことにより、画像内の位置座標に応じたシェーディング成分Ｓ（ｘ，ｙ）の変化を算出することができる。そして、このシェーディング成分の変化から、画像内におけるシェーディング成分Ｓ（ｘ，ｙ）の分布を抽出することができる。

さらに、次式（１−２）に示すように、抽出されたシェーディング成分Ｓ（ｘ，ｙ）を輝度Ｉ（ｘ，ｙ）から除去することにより、被写体成分Ｔ（ｘ，ｙ）のみの画像を取得することができる。
Ｔ（ｘ，ｙ）＝Ｉ（ｘ，ｙ）／Ｓ（ｘ，ｙ） …（１−２）

なお、撮像条件を変更する方法としては、例えば、標本ＳＰに対して撮像視野を平行移動させる、倍率を変更する、標本ＳＰに対して撮像視野を回転させるといった方法がある。本実施の形態１においては、標本ＳＰに対して撮像視野を平行移動させる方法により、撮像条件を変更する。

次に、顕微鏡システム１の動作について説明する。図３は、顕微鏡システム１の動作を示すフローチャートである。
まず、ステップＳ１０において、顕微鏡装置１０は、制御部１６０の制御の下で、標本ＳＰに対して撮像視野を変化させながら、撮像視野の一部が重複する複数の画像を撮像する。より詳細には、顕微鏡装置１０は、図４に示すように、標本ＳＰが載置された標本ステージ１０１を平行移動させつつ、画像取得部１０４に撮像視野を撮像させる。

その際、本実施の形態１において、顕微鏡装置１０は、図５に示すように、撮像を行う毎に標本ステージ１０１をＸ方向に移動量ｄｘ、及び、Ｙ方向に移動量ｄｙだけ移動させ、隣接する撮像視野Ｖ_j、Ｖ_j+1間で一部の領域Ｃ１が重なるようにして撮像を行う（ｊは撮像順序を示す番号であり、ｊ＝１、２、…）。なお、撮像視野は、Ｘ方向又はＹ方向のいずれか一方にのみ移動させても良い。画像取得部１０４は、このように撮像を行って取得した画像データを画像処理装置１１に出力する。画像処理装置１１は、画像取得部１０４から出力された画像データを、記憶部１４０に一旦記憶させる。

ステップＳ１１において、制御部１６０は、記憶部１４０から画像データを読み出し、該画像データに対応する複数の画像を演算部１５０に入力する。具体的には、図６（ａ）及び（ｂ）に示すように、撮像視野Ｖ_jを写した画像Ｍ_jと、撮像視野Ｖ_j+1を写した画像Ｍ_j+1とが演算部１５０に入力される。画像Ｍ_j、Ｍ_j+1間においては、画像Ｍ_j内の座標（ｘ_i，ｙ_i）の画素と、画像Ｍ_j+1内の座標（ｘ’_i，ｙ’_i）（＝（ｘ_i−ｄｘ，ｙ_i−ｄｙ））の画素とが、互いに画素となる。互いに対応する画素には、同じ被写体が写っている。以下、画像Ｍ_j、Ｍ_j+1の水平方向のサイズをｗとし、垂直方向のサイズをｈとする。

続くステップＳ１２において、画像調整部１５２は、位置検出部１０５ｂからの出力値に基づいて、複数の画像間で共通部分が一致するように画像の位置調整を行う。例えば、画像Ｍ_j、Ｍ_j+1の場合、図７に示すように、画像Ｍ_jに対して画像Ｍ_j+1を撮像視野の移動量ｄｘ及びｄｙの分だけずらし、領域Ｃ１同士が重なるように位置調整がなされる。

ステップＳ１３において、差分取得部１５３は、各画像の共通部分に含まれる画素の画素値から輝度を算出し、複数の画像間で互いに対応する画素の輝度の比を算出する。例えば、画像Ｍ_j、Ｍ_j+1の場合、画像Ｍ_j内の座標（ｘ_i，ｙ_i）における画素の輝度Ｉ_j、画像Ｍ_j+1内の座標（ｘ’_i，ｙ’_i）における画素の輝度Ｉ_j+1との比が算出される。

ここで、輝度Ｉ_j、Ｉ_j+1はそれぞれ、次式（１−３）、（１−４）により与えられる。
Ｉ_j（ｘ_i，ｙ_i）＝Ｔ（ｘ_i，ｙ_i）×Ｓ（ｘ_i，ｙ_i） …（１−３）
Ｉ_j+1（ｘ’_i，ｙ’_i）＝Ｔ（ｘ’_i，ｙ’_i）×Ｓ（ｘ’_i，ｙ’_i）
…（１−４）

また、対応する画素間において被写体成分は等しいから、次式（１−５）が成り立つ。
Ｔ（ｘ_i，ｙ_i）＝Ｔ（ｘ’_i，ｙ’_i） …（１−５）

さらに、座標（ｘ’_i，ｙ’_i）＝（ｘ_i−ｄｘ，ｙ_i−ｄｙ）を用いると、式（１−４）は次式（１−６）のように書き換えることができる。
Ｉ_j+1（ｘ’_i，ｙ’_i）＝Ｔ（ｘ_i，ｙ_i）×Ｓ（ｘ_i−ｄｘ，ｙ_i−ｄｙ）
…（１−６）

式（１−３）及び（１−６）より、次式（１−７）に示す関係が成り立つ。
Ｉ_j+1（ｘ’_i，ｙ’_i）／Ｉ_j（ｘ_i，ｙ_i）
＝Ｓ（ｘ_i−ｄｘ，ｙ_i−ｄｙ）／Ｓ（ｘ_i，ｙ_i） …（１−７）
即ち、輝度の比Ｉ_j+1（ｘ’_i，ｙ’_i）／Ｉ_j（ｘ_i，ｙ_i）は、シェーディング成分の変化に相当する。

ステップＳ１３においては、このような比Ｉ_j+1（ｘ’_i，ｙ’_i）／Ｉ_j（ｘ_i，ｙ_i）が、共通部分である領域Ｃ１内の全ての座標に関して算出される。

ステップＳ１４において、シェーディング検出部１５４は、画像Ｍ_j、Ｍ_j+1内のシェーディングの影響を近似するシェーディングモデルを作成し、領域Ｃ１内の各座標（ｘ_i，ｙ_i）における比Ｉ_j+1（ｘ’_i，ｙ’_i）／Ｉ_j（ｘ_i，ｙ_i）に基づいて、シェーディングモデルを修正する。

各画像に生じるシェーディングの影響は、原理的には式（１−７）に従うはずであるが、実際には、対応する画素同士においても微妙に輝度が異なったり、シェーディングにばらつきがあったりするので、領域Ｃ１内の全ての座標において式（１−７）が成り立つわけではない。そこで、シェーディング成分Ｓを表すモデルを設定し、次式（１−８）に示す誤差評価関数により、評価値Ｋ１が最小となるときのシェーディング成分Ｓを求める。

具体的なシェーディングモデルの一例として、本実施の形態１においては、次式（１−９）に示すように、画像の中心座標（ｗ／２，ｈ／２）を通る２次曲面を表す関数を用いる。ここで、２次曲面を用いるのは、シェーディング成分は一般に、画像の中心付近で少なく、画像の中心から離れるに従って増加するからである。
Ｓ（ｘ，ｙ；ａ）＝１−ａ｛（ｘ−ｗ／２）²＋（ｙ−ｈ／２）²｝ …（１−９）

従って、式（１−９）を用いた誤差評価関数（１−１０）によって与えられる評価値Ｋ１’が最小となるときの２次の係数（パラメータ）ａを求めることにより（式（１−１０’））、シェーディング成分Ｓを取得することができる。

シェーディング検出部１５４は、このパラメータａを、図８に示す演算処理により求める。図８は、シェーディング検出部１５４の動作を示すフローチャートである。
まず、ステップＳ１５１において、シェーディング検出部１５４は、パラメータａを初期化する。

続くステップＳ１５２において、シェーディング検出部１５４は、パラメータａ及びステップＳ１３において算出された領域Ｃ１内の全ての座標（ｘ_i，ｙ_i）における輝度の比Ｉ_j+1（ｘ’_i，ｙ’_i）／Ｉ_j（ｘ_i，ｙ_i）を誤差評価関数（１−１０）に代入して、評価値Ｋ１’を算出する。

ステップＳ１５３において、シェーディング検出部１５４は、評価値Ｋ１’が所定の閾値未満であるか否かを判定する。閾値の大きさは、後に続く繰り返し処理の各回において代入されるパラメータａに基づいて補正画像を作成した場合に、パラメータａが互いに異なる補正画像間で差異が明確に認識されない程度に、充分小さい値を経験的に設定するものとする。この閾値は、繰り返し処理を終了させるために設定するものであるので、繰り返し処理を終了する条件として、繰り返しにおける評価値の変化が充分に小さくなったときに終了するなど、他の終了条件を設定しても良い。評価値Ｋ１’が所定の閾値以上である場合（ステップＳ１５３：Ｎｏ）、シェーディング検出部１５４は、パラメータａを修正する（ステップＳ１５４）。その後、動作はステップＳ１５２に戻る。一方、評価値Ｋ１’が所定の閾値未満である場合（ステップＳ１５３：Ｙｅｓ）、シェーディング検出部１５４は、そのときのパラメータａを、式（１−９）におけるパラメータとして決定する（ステップＳ１５５）。この決定されたパラメータａを含む式（１−９）が、領域Ｃ１における修正済みのシェーディングモデルを表す式である。
その後、動作はメインルーチンに戻る。

ステップＳ１５において、シェーディング推定部１５５は、修正されたシェーディングモデルの適用範囲を、領域Ｃ１以外の画像Ｍ_j、Ｍ_j+1内の領域に拡張し、画像内の全領域におけるシェーディングを補正するための補正用画像を作成する。補正用画像は、シェーディングモデルＳ（ｘ，ｙ；ａ）を各画素の輝度とする画像である。

ステップＳ１６において、画像補正部１５６は、ステップＳ１５において作成された補正用画像を用いて、画像の補正を行う。即ち、補正対象である画像内の各画素の輝度Ｉ（ｘ，ｙ）を取得し、次式（１−１１）により、被写体成分Ｔ（ｘ，ｙ）を算出する。
Ｔ（ｘ，ｙ）＝Ｉ（ｘ，ｙ）／Ｓ（ｘ，ｙ；ａ） …（１−１１）
それにより、シェーディングの影響が除去された補正済みの画像が取得される。なお、補正対象の画像は、シェーディング補正用画像の作成に用いた画像Ｍ_j、Ｍ_j+1に限らず、顕微鏡装置１０において撮像した他の画像であっても良い。

さらに、ステップＳ１７において演算部１５０は、補正済みの画像を出力する。それに応じて、制御部１６０は、補正済みの画像を表示部１３０に表示させると共に、補正済みの画像に対応する画像データを記憶部１４０に記憶させる。

以上説明したように、実施の形態１によれば、補正対象である画像自体の輝度に基づいてシェーディング補正用画像を作成するので、顕微鏡装置１０において経時的なシェーディングの変化が生じている場合であっても、当該画像に対して精度の高いシェーディング補正を施すことができる。また、実施の形態１によれば、画像内の輝度を被写体成分とシェーディング成分とを分離し、シェーディング成分のみを抽出してシェーディング補正用画像を作成するので、シェーディング補正の精度を向上させることが可能となる。

（実施の形態２）
次に、本発明の実施の形態２について説明する。
実施の形態２に係る顕微鏡システムの全体的な構成及び動作は、実施の形態１と共通であり、実施の形態２においては、撮像視野の倍率を変化させて複数の画像を撮像し、これらの画像を用いてシェーディング補正用画像を作成することを特徴とする。従って、以下においては、互いに倍率が異なる複数の画像を用いてシェーディング補正用画像を作成する処理についてのみ説明する。

まず、顕微鏡システム１は、図９に示すように、標本ＳＰの視野領域Ｖ_jを撮像して画像Ｍ_jを取得する。次いで、ズームを調整して倍率を変更し、同じ視野領域Ｖ_jを撮像して画像Ｍ_j+1を取得する。この場合、画像Ｍ_j、Ｍ_j+1間では、画像全体が、同じ被写体が写った共通部分となる。

演算部１５０において、画像調整部１５２は、画像Ｍ_jの中心Ｏと、Ｍ_j+1の中心Ｏ’とが一致するように、画像Ｍ_j、Ｍ_j+1の位置調整を行う。図１０に示すように、画像Ｍ_j、Ｍ_j+1間では、画像Ｍ_jの中心Ｏから距離ｒ_iだけ離れ、所定の軸からの回転角φ_iに位置する画素Ｐ_iと、画像Ｍ_j+2の中心Ｏ’から距離ｒ’_iだけ離れ、所定の軸からの回転角φ_iに位置する画素Ｐ’_iとが互いに対応する画素となる。

続いて、比率取得部１５３は、画像Ｍ_j、Ｍ_j+1間で対応する画素の輝度の比を算出する。ここで、画像Ｍ_jの倍率をｍ_j、画像Ｍ_j+1の倍率をｍ_j+1とすると、距離ｒ’_iは次式（２−１）によって表すことができる。
ｒ’_i＝（ｍ_j+1／ｍ_j）×ｒ_i …（２−１）

上述したとおり、各画素の輝度Ｉは、被写体成分Ｔとシェーディング成分Ｓとで構成されるから、画素Ｐ_iの輝度Ｉ_j（ｒ_i，φ_i）及び画素Ｐ’_iの輝度Ｉ_j+1（ｒ’_i，φ_i）は、それぞれ、次式（２−２）及び（２−３）によって表すことができる。
Ｉ_j（ｒ_i，φ_i）＝Ｔ（ｒ_i，φ_i）×Ｓ（ｒ_i，φ_i） …（２−２）
Ｉ_j+1（ｒ’_i，φ_i）＝Ｔ_j+1（ｒ’_i，φ_i）×Ｓ_j+1（ｒ’_i，φ_i）
…（２−３）

ここで、一般には、撮像倍率を変更するとシェーディングも変化する。しかしながら、倍率の変化が小さい場合には、シェーディングの変化は極めて小さく、無視することができる。それに対して、画像内の被写体像の変化は倍率の変化と厳密に一致する。即ち、次式（２−４）のとおりである。
Ｔ_j+1（ｒ’_i，φ_i）＝Ｔ（ｒ_i，φ_i） …（２−４）

従って、式（２−２）〜（２−４）より、輝度の比Ｉ_j+1（ｒ’_i）／Ｉ_j（ｒ_i）は次式（２−５）によって与えられる。
Ｉ_j+1（ｒ’_i，φ_i）／Ｉ_j（ｒ_i，φ_i）＝Ｓ_j+1（ｒ’_i，φ_i）／Ｓ（ｒ_i，φ_i）
…（２−５）

式（２−５）に示すとおり、輝度の比Ｉ_j+1（ｒ’_i，φ_i）／Ｉ_j（ｒ_i，φ_i）は、シェーディング成分の変化に相当する。ここで、一般にシェーディング成分は、光学系の光軸Ｌからの距離に応じて軸対称に変化するから、式（２−５）は次式（２−６）に書き換えることができる。
Ｉ_j+1（ｒ’_i，φ_i）／Ｉ_j（ｒ_i，φ_i）＝Ｓ_j+1（ｒ’_i）／Ｓ（ｒ_i）
…（２−６）
比率取得部１５３は、このような比Ｉ_j+1（ｒ’_i，φ_i）／Ｉ_j（ｒ_i，φ_i）を、画像Ｍ_j内の全ての画素に対して算出する。

続いて、シェーディング検出部１５４は、画像Ｍ_j、Ｍ_j+1内のシェーディングの影響を近似するシェーディングモデルを作成し、各座標における比Ｉ_j+1（ｒ’_i，φ_i）／Ｉ_j（ｒ_i，φ_i）に基づいて、シェーディングモデルを修正する。本実施の形態２の場合、シェーディングモデルを修正するための誤差評価関数は次式（２−７）によって与えられ、評価値Ｋ２を最小にするシェーディングモデルＳ（ｒ_i）を求める演算が行われる。

また、具体的なシェーディングモデルの一例として、本実施の形態２においては、次式（２−８）に示すように、画像Ｍ_jの中心からの距離ｒに依存して変化する２次曲面を表す関数を用いる。
Ｓ（ｒ；ｂ）＝１−ｂ×ｒ^２ …（２−８）

従って、式（２−７）を用いた誤差評価関数（２−９）によって与えられる評価値Ｋ２’が最小となるときの２次の係数であるパラメータｂを求めることにより（式（２−９’））、シェーディング成分Ｓを取得することができる。

シェーディング検出部１５４は、このパラメータｂを、図８に示す演算処理により求める。ただし、図８において、パラメータａはパラメータｂに読み替えられ、評価値Ｋ１’は評価値Ｋ２’に読み替えられる。それによって決定されたパラメータｂを含む式（２−８）が、修正済みのシェーディングモデルを表す式である。

続いて、シェーディング推定部１５５は、修正されたシェーディングモデルの適用範囲を画像Ｍ_j+1全体（即ち、画像Ｍ_jよりも若干広い範囲）に拡張し、画像Ｍ_j及び画像Ｍ_j+1内の全領域におけるシェーディングを補正するための補正用画像を作成する。
その後の処理については、実施の形態１と同様である。

以上説明したように、実施の形態２によれば、画像Ｍ_j全体の画素の輝度に基づいて補正用画像を作成するので、精度の高いシェーディング補正を行うことができる。

（実施の形態３）
次に、本発明の実施の形態３について説明する。
実施の形態３に係る顕微鏡システムの全体的な構成及び動作は、実施の形態１と共通であり、実施の形態３においては、撮像視野を回転させて複数の画像を撮像し、これらの画像を用いてシェーディング補正用画像を作成することを特徴とする。従って、以下においては、ＸＹ平面における座標軸が互いに交差する複数の画像を用いてシェーディング補正用画像を作成する処理についてのみ説明する。

まず、顕微鏡システム１は、図１１及び図１２に示すように、標本ＳＰの視野領域Ｖ_jを撮像して画像Ｍ_jを取得する。次いで、標本ステージ１０１を所定の回転中心点に対して角度ｄθだけＸＹ平面内で回転させ、視野領域Ｖ_j+1を撮像して画像Ｍ_j+1を取得する。なお、実施の形態３においては、回転中心点を撮像領域Ｖ_jの中心Ｏに設定する。この場合、画像Ｍ_j、Ｍ_j+1間では、図１１に示す領域Ｃ３が、同じ被写体が写った共通部分となる。また、図１２（ａ）及び（ｂ）に示すように、画像Ｍ_jの回転中心点（中心Ｏ）から距離ｒ_iだけ離れ、所定の軸からの角度θ_iに位置する画素Ｑ_iと、画像Ｍ_j+1の回転中心点（中心Ｏ）から距離ｒ_iだけ離れ、所定の軸からの角度θ’_i（＝θ_i−ｄθ）に位置する画素Ｑ’_iとが、互いに対応する画素となる。

演算部１５０において、画像調整部１５２は、図１３に示すように、画像Ｍ_jに対して画像Ｍ_j+1を角度ｄθだけ回転させ、画像Ｍ_j、Ｍ_j+1間で領域Ｃ３が一致するよう、位置調整を行う。

続いて、比率取得部１５３は、画像Ｍ_j、Ｍ_j+1間で対応する画素の輝度の差分を算出する。具体的には、画素Ｑ_iの輝度Ｉ_j（ｒ_i，θ_i）と、画素Ｑ’_iの輝度Ｉ_j+1（ｒ_i，θ’_i）との比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）を算出する。

上述したとおり、各画素の輝度Ｉは、被写体成分Ｔとシェーディング成分Ｓとで構成されるから、画素Ｑ_iの輝度Ｉ_j（ｒ_i，θ_i）及び画素Ｑ’_iの輝度Ｉ_j+1（ｒ_i，θ’_i）は、それぞれ、次式（３−１）及び（３−２）によって表すことができる。
Ｉ_j（ｒ_i，θ_i）＝Ｔ（ｒ_i，θ_i）×Ｓ（ｒ_i，θ_i） …（３−１）
Ｉ_j+1（ｒ_i，θ’_i）＝Ｔ_j+1（ｒ_i，θ’_i）×Ｓ_j+1（ｒ_i，θ’_i）
…（３−２）

また、対応する画素間において被写体成分は等しいから、次式（３−３）が成り立つ。
Ｔ（ｒ_i，θ_i）＝Ｔ_j+1（ｒ_i，θ’_i） …（３−３）
従って、式（３−１）〜（３−３）より、比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）は次式（３−４）によって与えられる。
Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）＝Ｓ_j+1（ｒ_i，θ’_i）／Ｓ_j（ｒ_i，θ_i）
…（３−４）
即ち、輝度の比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）は、シェーディング成分の変化に相当する。

ここで、シェーディング成分Ｓには、レンズに起因して生じるシェーディング成分（以下、シェーディング成分Ｓｌとする）と、照明等、レンズ以外の要因によって生じるシェーディング成分（以下、シェーディング成分Ｓｍとする）とが含まれる。即ち、シェーディング成分Ｓは、次式（３−５）によって与えられる。
Ｓ（ｒ_i，θ_i）＝Ｓｍ（ｒ_i，θ_i）×Ｓｌ（ｒ_i，θ_i） …（３−５）

この内、シェーディング成分Ｓｌは光軸Ｌを中心とする軸対称に生じるため、次式（３−６）に示すように、画像の中心から画素までの距離ｒ_iのみによってモデル化することができる。
Ｓ（ｒ_i，θ_i）＝Ｓｍ（ｒ_i，θ_i）×Ｓｌ（ｒ_i） …（３−６）

式（３−５）及び（３−６）より、式（３−４）を次式（３−７）のように書き換えることができる。

式（３−７）より、対応する画素間の比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）は、照明等に起因するシェーディング成分Ｓｍの変化を表すと言える。比率取得部１５３は、このような比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）を、領域Ｃ３内の全ての画素に対して算出する。

続いて、シェーディング検出部１５４は、画像Ｍ_j、Ｍ_j+1内のシェーディングの影響を近似するシェーディングモデルを作成し、各座標における比Ｉ_j+1（ｒ_i，θ’_i）／Ｉ_j（ｒ_i，θ_i）に基づいて、シェーディングモデルを修正する。本実施の形態３の場合、シェーディングモデルを修正するための誤差評価関数は次式（３−８）によって与えられ、評価値Ｋ３を最小にするシェーディングモデルＳｍ（ｒ_i，θ_i）を求める演算が行われる。

具体的なシェーディングモデルの一例として、本実施の形態３においては、次式（３−９）に示すように、画像Ｍ_jの中心からの距離ｒ及び角度θに依存して変化する２次曲面を表す関数を用いる。

式（３−９）において、ｒ₀及びθ₀は所定の定数である。

従って、式（３−９）を用いた誤差評価関数（３−１０）によって与えられる評価値Ｋ３’が最小となるときの２次の係数であるパラメータｃを求めることにより（式（３−１０’））、シェーディング成分Ｓを取得することができる。

シェーディング検出部１５４は、このパラメータｃを、図８に示す演算処理により求める。ただし、図８において、パラメータａはパラメータｃに読み替えられ、評価値Ｋ１’は評価値Ｋ３’に読み替えられる。それによって決定されたパラメータｃを含む式（３−９）が、照明等に起因するシェーディング成分Ｓｍの修正済みモデルを表す式である。

続いて、シェーディング推定部１５５は、画像Ｍ_j及び画像Ｍ_j+1内の全領域におけるシェーディングを補正するための補正用画像を作成する。ここで、レンズに起因するシェーディング成分Ｓｌは経時的な変化が少ないため、コサイン四乗則により正確に予測することができる。コサイン四乗則は、レンズに照度Ｉ₀の光が入射するとき、入射光の光軸Ｌに対する角度Θと、入射後の光の照度Ｉ’との関係を表す。
Ｉ’＝Ｉ₀ｃｏｓ⁴Θ …（３−１１）

シェーディング推定部１５５は、式（３−１１）に基づいて、レンズに起因するシェーディング成分Ｓｌを表すシェーディングモデルＳｌ（ｒ）を作成する。
シェーディングモデルＳｌ(ｒ)としては、例えば、式（３−１１）に基づいて、光軸中心からの距離ｒの値に応じたシェーディング量をルックアップテーブル化したものを用いることができる。記憶部１４０は、レンズごとに作成された複数種類のルックアップテーブルを格納しており、シェーディング推定部１５５は、選択されたレンズに対応するルックアップテーブルを記憶部１４０から読み出し、このルックアップテーブルを使用して、シェーディング成分Ｓｌを求める。

さらに、シェーディング推定部１５５は、式（３−９）によって与えられるシェーディングモデルを画像全体（領域Ｃ１以外の画像Ｍ_j及び画像Ｍ_j+1の領域）に拡張すると共に、上記ルックアップテーブルに基づくシェーディングモデルＳｌ（ｒ）を用いて、次式（３−１２）に示すように、画像Ｍ_j及び画像Ｍ_j+1内の全領域におけるトータルのシェーディングを補正するための補正用画像を作成する。
Ｓ_TOTAL＝Ｓｌ（ｒ）×Ｓｍ（ｒ，θ） …（３−１２）
その後の処理については、実施の形態１と同様である。

以上説明したように、実施の形態３によれば、照明等に起因する経時的な変化が比較的大きいシェーディングについても精度良く補正を行うことができる。従って、多量の画像を撮像した場合や、撮像時間が長引いたりした場合においても、各画像に対して正確なシェーディング補正を施すことが可能となる。

（実施の形態４）
次に、本発明の実施の形態４について説明する。
図１４は、実施の形態４に係る顕微鏡システムの構成を示すブロック図である。図１４に示すように、実施の形態４に係る顕微鏡システム２は、図１に示す画像処理装置１１の代わりに画像処理装置２０を備える。画像処理装置２０は、図１に示す演算部１５０に対してバーチャルスライド（ＶＳ）画像作成部２０１をさらに有する演算部２００を備える。ＶＳ画像作成部２０１以外の画像処理装置２０及び顕微鏡システム２の構成については、実施の形態１において説明したものと同様である。

ＶＳ画像作成部２０１は、顕微鏡装置１０が標本ＳＰに対する撮像視野を平行移動させながら撮像した複数の画像を繋ぎ合わせて、標本ＳＰ全体に対応する画像（ＶＳ画像）を作成する。

次に、顕微鏡システム２の動作について説明する。図１５は、顕微鏡システム２の動作を示すフローチャートである。また、図１６は、ＶＳ画像作成部が作成したバーチャルスライド画像を示す模式図である。

まず、ステップＳ１０において、顕微鏡装置１０は、制御部１６０の制御の下で、標本ＳＰに対する撮像視野を平行移動させながら、撮像視野の一部が重複する複数の画像を撮像する。
続くステップＳ１１及びＳ１２の動作は、実施の形態１と同様である。

ステップＳ１２に続くステップＳ２０において、ＶＳ画像作成部２０１は、図１６に示すように、画像調整部１５２による位置調整済みの複数の画像を繋ぎ合わせて、バーチャルスライド画像ＶＳを作成する。このバーチャルスライド画像ＶＳにおいて、Ｘ方向で隣接する画像Ｍ_(k,l)とＭ_(K+1,l)間では、領域Ｃ_kが互いに重複する（即ち、被写体が一致する）共通部分であり、Ｙ方向で隣接する画像Ｍ_(k,l)とＭ_(k,l+1)間では、領域Ｃ_lが互いに重複する共通部分である（ｋ，ｌ＝１、２、…）。

続くステップＳ１３〜Ｓ１５の動作は、実施の形態１と同様である。ただし、この際の演算対象は、バーチャルスライド画像ＶＳ内の全ての共通部分（領域Ｃ_k及び領域Ｃ_l）であっても良いし、一部の共通部分であっても良い。後者の場合、具体的には、全ての領域Ｃ_kのみ、全ての領域Ｃ_lのみ、最初と最後の領域Ｃ_k若しくは領域Ｃ_lのみ、バーチャルスライド画像ＶＳの中央近傍の領域Ｃ_k若しくは領域Ｃ_lのみ、所定間隔で抽出した領域Ｃ_k若しくは領域Ｃ_l、又は、ランダムに選択した領域Ｃ_k若しくは領域Ｃ_lのいずれであっても良い。複数の共通部分が演算対象となる場合には、それらの共通部分に対応する複数の補正用画像が作成される。

ステップＳ１５に続くステップＳ３０において、画像補正部１５６は、ステップＳ１５において作成された補正用画像を用いて、バーチャルスライド画像ＶＳを構成する各画像Ｍ_(k,l)を補正する。この際、ステップＳ１５において複数の補正用画像が作成された場合には、補正用画像の基となった共通部分の位置に応じて補正用画像を決定すると良い。例えば、全ての領域Ｃ_kから補正用画像が作成された場合には、Ｘ方向において隣接する画像Ｍ_(k,l)、Ｍ_(K+1,l)を、両画像の共通部分である領域Ｃ_kに基づく補正用画像を用いて補正する。また、全ての領域Ｃ_lから補正用画像が作成された場合には、Ｙ方向において隣接する画像Ｍ_(k,l)、Ｍ_(K,l+1)を、両画像の共通部分である領域Ｃ_lに基づく補正用画像を用いて補正する。或いは、所定間隔で抽出された領域Ｃ_k若しくは領域Ｃ_l、又は、ランダムに選択された領域Ｃ_k若しくは領域Ｃ_lから補正用画像が作成された場合には、領域Ｃ_k若しくは領域Ｃ_lから所定範囲内の画像Ｍ_(k,l)を、当該領域Ｃ_k若しくは領域Ｃ_lに基づく補正用画像を用いて補正することとしても良い。

続くステップＳ３１において、演算部２００は、補正済みの画像を繋ぎ合わせたバーチャルスライド画像ＶＳを出力する。それに応じて、制御部１６０は、バーチャルスライド画像ＶＳを表示部１３０に表示させると共に、バーチャルスライド画像ＶＳに対応する画像データを記憶部１４０に記憶させる。

以上説明したように、実施の形態４によれば、個別の撮像視野に生じたシェーディングが補正された画像によってバーチャルスライド画像を構成するので、画質の良いバーチャルスライド画像を取得することが可能となる。

（変形例）
上記実施の形態１〜３においては、顕微鏡装置１０によって取得された画像に対する処理を説明したが、画像処理装置１１は、顕微鏡装置１０以外の各種撮像機器によって取得された画像を処理することも可能である。例えば、画像処理装置１１を、パノラマ画像の撮像が可能なディジタルカメラに適用しても良い。この場合、視野の端部が重複する複数の画像を撮像し、視野の重複部分に基づいて作成した補正用画像を用いて各画像を補正して、それらの画像を繋ぎ合わせることにより、画質の良いパノラマ画像を生成することが可能である。

本発明は、上述した各実施の形態１〜４及び変形例そのままに限定されるものではなく、各実施の形態１〜４に開示されている複数の構成要素を適宜組み合わせることによって、種々の発明を形成することができる。例えば、実施の形態１〜４に示される全構成要素からいくつかの構成要素を除外して形成してもよい。あるいは、異なる実施の形態に示した構成要素を適宜組み合わせて形成してもよい。

１、２ 顕微鏡システム
１０ 顕微鏡装置
１１、２０ 画像処理装置
１００ アーム
１０１ 標本ステージ
１０２ 対物レンズ
１０３ 鏡筒
１０４ 画像取得部
１０５ ステージ位置変更部
１０５ａ モータ
１０５ｂ 位置検出部
１０６ 三眼鏡筒ユニット
１０７ 接眼レンズユニット
１０８ レボルバ
１１０ 入力部
１２０ 画像入力部
１３０ 表示部
１４０ 記憶部
１４１ 画像処理プログラム
１５０、２００ 演算部
１５１ 補正用画像作成部
１５２ 画像調整部
１５３ 比率取得部
１５４ シェーディング検出部
１５５ シェーディング推定部
１５６ 画像補正部
１６０ 制御部
２０１ ＶＳ画像作成部

Claims (9)

1. 被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力手段と、
   前記複数の画像を補正する際に用いられる補正用画像を作成する補正画像作成手段と、
   を備え、
   前記補正画像作成手段は、
   前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得手段と、
   前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出手段と、
   前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定するシェーディング推定手段と、
   を有することを特徴とする画像処理装置。
2. 前記補正画像作成手段は、前記任意画像の位置調整を行い、前記共通部分における前記被写体を一致させる画像調整手段をさらに有し、
   前記比率取得手段は、前記画像調整手段により位置調整された状態で、前記比を取得することを特徴とする請求項１に記載の画像処理装置。
3. 前記複数の画像は、撮像視野を平行移動させて撮像された画像であることを特徴とする請求項１又は２に記載の画像処理装置。
4. 前記複数の画像は、撮像倍率を変化させて撮像された画像であることを特徴とする請求項１又は２に記載の画像処理装置。
5. 前記複数の画像は、撮像視野を回転させて撮像された画像であることを特徴とする請求項１又は２に記載の画像処理装置。
6. 請求項１〜５のいずれか１項に記載の画像処理装置と、
   前記被写体を撮像する撮像手段と、
   を備えることを特徴とする撮像装置。
7. 請求項３又は５に記載の画像処理装置と、
   前記被写体としての標本を載置可能なステージと、前記ステージに対向して設けられる光学系と、前記光学系を介して前記標本上に設定される視野を撮像して画像を取得する画像取得部と、前記ステージと前記光学系との内の少なくとも一方を、前記光学系の光軸と直交する方向に移動させることにより、前記撮像視野を変化させるステージ位置変更部と、を有する顕微鏡装置と、
   を備えることを特徴とする顕微鏡システム。
8. 被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、
   前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得ステップと、
   前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、
   前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定して前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、
   を含むことを特徴とする画像処理方法。
9. 被写体の少なくとも一部が共通する部分を持つ複数の画像を入力する画像入力ステップと、
   前記複数の画像のうち、前記被写体が共通する共通部分を持つ２つの任意画像から、画素の輝度の比を取得する比率取得ステップと、
   前記比に基づいて前記任意画像におけるシェーディング成分を検出するシェーディング検出ステップと、
   前記シェーディング成分に基づき、前記共通部分以外の各画像内の領域におけるシェーディングの影響を推定して前記複数の画像を補正する際に用いられる補正用画像を作成する補正用画像作成ステップと、
   をコンピュータに実行させることを特徴とする画像処理プログラム。

Images