"光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理装置であって、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得手段と、前記結像特性取得手段で取得された結像特性に基づいて、焦点位置から光軸方向に離れた物体が画像に対して与える影響が所定の条件よりも大きくなる影響範囲を求め、前記複数の第1の画像データの中から前記影響範囲内に焦点位置が含まれる第1の画像データを選択する選択手段と、前記選択手段により選択された第1の画像データを用いて、前記第2の画像データを生成する生成手段と、を有することを特徴とする画像処理装置。前記指定された被写界深度に基づいて、前記第2の画像データの焦点位置に相当する合成焦点面を決定する決定手段をさらに有し、前記決定手段は、前記合成焦点面を中心とする前記被写界深度内に物体が最も多く含まれるように、前記合成焦点面の位置を決定し、前記選択手段は、前記決定手段で決定された合成焦点面を中心として前記影響範囲を設定することを特徴とする請求項1に記載の画像処理装置。光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理装置であって、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得手段と、前記指定された被写界深度に基づいて、前記第2の画像データの焦点位置に相当する合成焦点面を決定する決定手段と、前記結像特性取得手段で取得された結像特性と前記決定手段で決定された合成焦点面とに基づいて、前記前記第2の画像データを生成する生成手段と、を有することを特徴とする画像処理装置。前記決定手段は、前記合成焦点面を中心とする前記被写界深度内に物体が最も多く含まれるように、前記合成焦点面の位置を決定することを特徴とする請求項3に記載の画像処理装置。前記決定手段は、前記複数の第1の画像データの各々について、画像内に焦点が合った領域がどの程度存在するかを表す画像評価値を計算し、複数の焦点位置のそれぞれについて、前記焦点位置を中心とする前記被写界深度内で第1の画像データの画像評価値を積分した画像評価積分値を計算し、前記画像評価積分値が最も良い値となる焦点位置を前記合成焦点面の位置に決定することを特徴とする請求項3又は4に記載の画像処理装置。前記決定手段で決定された合成焦点面を基準として複数の合成焦点面を設定し、前記複数の合成焦点面のそれぞれに対応する第2の画像データが生成されることを特徴とする請求項3~5のうちいずれか1項に記載の画像処理装置。前記複数の合成焦点面は、その間隔が前記指定された被写界深度に等しくなるように設定されることを特徴とする請求項6に記載の画像処理装置。前記被写界深度を指定する情報は、対物レンズ倍率、開口絞り、又はその両方を指定する情報であることを特徴とする請求項1~7のうちいずれか1項に記載の画像処理装置。前記被写界深度を指定する情報は、開口数を指定する情報であることを特徴とする請求項1~7のうちいずれか1項に記載の画像処理装置。前記生成手段は、前記仮想的な光学系の結像特性に基づくコンボリューション処理を適用することによって、前記第1の画像データの各々に対し焦点位置からの距離に応じたぼけを付与した後、ぼけを付与した第1の画像データを合成して前記第2の画像データを生成することを特徴とする請求項1~9のうちいずれか1項に記載の画像処理装置。前記生成手段は、前記仮想的な光学系の結像特性に基づいて決定される重みを用いて、周波数空間上で前記第1の画像データを重み付け加算することにより、前記第2の画像データを生成することを特徴とする請求項1~9のうちいずれか1項に記載の画像処理装置。前記生成手段は、前記複数の第1の画像データを撮像した光学系の結像特性によるぼけをキャンセルし、且つ、前記仮想的な光学系の結像特性によるぼけを付与するフィルタを用いて、周波数空間上で前記第1の画像データに対し前記フィルタを適用することにより、前記第2の画像データを生成することを特徴とする請求項1~9のうちいずれか1項に記載の画像処理装置。前記結像特性は、PSFまたはOTFであることを特徴とする請求項1~12のうちいずれか1項に記載の画像処理装置。前記被写体を複数の領域に分割することで得られる複数の分割領域毎に前記複数の第1の画像データが取得され、前記分割領域毎に第2の画像データが生成されることを特徴とする請求項1~13のうちいずれか1項に記載の画像処理装置。第2の画像データの被写界深度として、互いに異なる複数の被写界深度を指定可能であることを特徴とする請求項1~14のうちいずれか1項に記載の画像処理装置。前記指定された複数の被写界深度のそれぞれに対応する複数の第2の画像データが生成されることを特徴とする請求項15に記載の画像処理装置。光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理方法であって、コンピュータが、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得ステップと、コンピュータが、前記結像特性取得ステップで取得された結像特性に基づいて、焦点位置から光軸方向に離れた物体が画像に対して与える影響が所定の条件よりも大きくなる影響範囲を求め、前記複数の第1の画像データの中から前記影響範囲内に焦点位置が含まれる第1の画像データを選択する選択ステップと、コンピュータが、前記選択ステップにより選択された第1の画像データを用いて、前記第2の画像データを生成する生成ステップと、を有することを特徴とする画像処理方法。光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理方法であって、コンピュータが、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得ステップと、コンピュータが、前記指定された被写界深度に基づいて、前記第2の画像データの焦点位置に相当する合成焦点面を決定する決定ステップと、コンピュータが、前記結像特性取得ステップにより取得された結像特性と前記決定ステップにより決定された合成焦点面とに基づいて、前記第2の画像データを生成する生成ステップと、を有することを特徴とする画像処理方法。請求項17又は18に記載の画像処理方法の各ステップをコンピュータに実行させることを特徴とするプログラム。しかしながら上述した従来の技術においては、以下のような問題があった。特許文献1の方法では、被写界深度が浅い画像か、もしくは全焦点画像かの2通りしか選定できない。いずれの画像も光学顕微鏡の目視観察像とは被写界深度或いは像のぼけが大きく異なる。病理医は、画像に含まれるぼけ情報から組織や細胞の3次元的な立体構造の把握をしており、倍率に応じたぼけは診断上有益な情報である。しかし特許文献1の方法では、被写界深度の浅い高倍画像か、全焦点画像の低倍画像のいずれかしか生成できないため、その中間倍率例えば対物レンズ20倍においてぼけ情報を元にした立体構造の推測が難しいという課題があった。特許文献2の方法によれば、Zスタック画像から任意の被写界深度の画像を生成することができる。しかしながら、特許文献2には、Zスタック画像の全てのレイヤー画像を使用する方法しか開示されていない。被写界深度の拡大処理に全てのレイヤー画像を使用すると、演算コストの増加、応答速度の低下が発生する課題がある。特にレイヤー画像の高解像度化や大サイズ化が進むほど、またレイヤー画像の枚数が増えるほど、この問題は深刻なものとなる。本発明は、このような問題点に鑑みてなされたものであり、バーチャルスライドシステムでのユーザ設定に応じた被写界深度の拡大処理において演算コストを抑制し、装置のスループット単位時間当たりの処理枚数や応答速度の向上を目的とする。本発明に係る画像処理装置は、光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理装置であって、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得手段と、前記結像特性取得手段で取得された結像特性に基づいて、焦点位置から光軸方向に離れた物体が画像に対して与える影響が所定の条件よりも大きくなる影響範囲を求め、前記複数の第1の画像データの中から前記影響範囲内に焦点位置が含まれる第1の画像データを選択する選択手段と、前記選択手段により選択された第1の画像データを用いて、前記第2の画像データを生成する生成手段と、を有することを特徴とする画像処理装置である。本発明に係る別の画像処理装置は、光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理装置であって、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得手段と、前記指定された被写界深度に基づいて、前記第2の画像データの焦点位置に相当する合成焦点面を決定する決定手段と、前記結像特性取得手段で取得された結像特性と前記決定手段で決定された合成焦点面とに基づいて、前記前記第2の画像データを生成する生成手段と、を有することを特徴とする画像処理装置である。本発明に係る画像処理方法は、光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理方法であって、コンピュータが、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得ステップと、コンピュータが、前記結像特性取得ステップで取得された結像特性に基づいて、焦点位置から光軸方向に離れた物体が画像に対して与える影響が所定の条件よりも大きくなる影響範囲を求め、前記複数の第1の画像データの中から前記影響範囲内に焦点位置が含まれる第1の画像データを選択する選択ステップと、コンピュータが、前記選択ステップにより選択された第1の画像データを用いて、前記第2の画像データを生成する生成ステップと、を有することを特徴とする画像処理方法である。本発明に係る別の画像処理方法は、光軸方向に焦点位置を変えながら被写体を撮像することにより得られた複数の第1の画像データから、前記第1の画像データよりも被写界深度が深い第2の画像データを生成する画像処理方法であって、コンピュータが、生成する第2の画像データの被写界深度を指定する情報に基づいて、前記指定された被写界深度をもつ仮想的な光学系の結像特性を取得する結像特性取得ステップと、コンピュータが、前記指定された被写界深度に基づいて、前記第2の画像データの焦点位置に相当する合成焦点面を決定する決定ステップと、コンピュータが、前記結像特性取得ステップにより取得された結像特性と前記決定ステップにより決定された合成焦点面とに基づいて、前記第2の画像データを生成する生成ステップと、を有することを特徴とする画像処理方法である。本発明に係るプログラムは、本発明に係る画像処理方法の各ステップをコンピュータに実行させることを特徴とするプログラムである。本発明では、バーチャルスライドシステムでのユーザ設定に応じた被写界深度の拡大処理において、演算コストを低減し、装置のスループットや応答速度を向上することができる。[実施例1]システムの全体構成図1は、本発明の撮像システムの一実施例である、バーチャルスライドシステムの構成を示している。バーチャルスライドシステムは、検体の撮像データを取得する撮像装置バーチャルスライドスキャナとも呼ぶ120とそのデータ処理制御を行う画像処理装置ホストコンピュータとも呼ぶ110およびその周辺機器から構成されている。画像処理装置110にはユーザからの入力を受け付けるキーボードやマウスなどの操作入力デバイス111、処理画像を表示するディスプレイ112が接続される。また画像処理装置110には、記憶装置113、他のコンピュータシステム114が接続されている。多数の検体プレパラートの撮像をバッチ処理で行う場合、画像処理装置110の制御の下で撮像装置120が各検体を順に撮像し、画像処理装置110が各検体の画像データに対し必要な処理を施す。そして得られた各検体の画像データは、大容量のデータストレージである記憶装置113又は他のコンピュータシステム114へ伝送され、蓄積される。撮像装置120での撮像プレ計測および本計測は、ユーザの入力を受けて画像処理装置110がコントローラ108へ指示を送り、次にコントローラ108が本計測ユニット101とプレ計測ユニット102を制御することで実現される。本計測ユニット101はプレパラート内の検体診断のための高精細な画像を取得する撮像ユニットである。プレ計測ユニット102は本計測に先立って撮像を行う撮像ユニットで、本計測で精度の良い画像取得をするための撮像制御情報取得を目的とした画像取得を行う。コントローラ108には変位計103が接続され、本計測ユニット101またはプレ計測ユニット102内のステージに設置されるプレパラートの位置や距離が測定できる構成となっている。変位計103は、本計測およびプレ計測を行うにあたり、プレパラート内の検体の厚みを計測するために用いられる。またコントローラ108には本計測ユニット101およびプレ計測ユニット102の撮像条件を制御するための、開口絞り制御104、ステージ制御105、照明制御106、センサ制御107が接続されている。そして、それぞれはコントローラ108からの制御信号に従って、開口絞り、ステージ、照明、イメージセンサの動作を制御する構成となっている。ステージにはプレパラートを光軸に垂直な方向に移動するXYステージと光軸に沿った方向に移動するZステージがある。それぞれ、XYステージは光軸に垂直な方向に分布する検体画像を、Zステージは奥行き方向に焦点位置を変えた画像を撮像するために用いる。図示しないが、撮像装置120は、複数のプレパラートをセット可能なラックと、ラックからステージ上の撮像位置へとプレパラートを送り出す搬送機構とが設けられている。バッチ処理の場合は、コントローラ108が搬送機構を制御することで、ラックからプレパラートを1枚ずつ、プレ計測ユニット102のステージ、本計測ユニット101のステージの順に送り出す。本計測ユニット101、プレ計測ユニット102には撮像した画像を用いてオートフォーカスを実現するAFユニット109が接続されている。AFユニット109は、コントローラ108を介して、本計測ユニット101、プレ計測ユニット102のステージの位置を制御することで合焦位置を探し出すことが出来る。オートフォーカスの方式は画像を用いるパッシブ型であり、公知の位相差検出方式やコントラスト検出方式が用いられる。本計測ユニット図2aは実施例1における本計測ユニット101の内部構成を示す図である。光源201aの光は照明光学系202aを通じ、光量ムラが無いように均一化され、ステージ203aの上に設置されたプレパラート204を照射する。プレパラート204は観察対象となる組織の切片や塗抹した細胞をスライドグラス上に貼り付け、封入剤とともにカバーグラスの下に固定したものであり、検体被写体を観察可能な状態に準備したものである。結像光学系205aは検体の像を拡大して撮像手段である撮像部207aに導くものである。プレパラート204を通った光は結像光学系205aを介して、撮像部207a上の撮像面で結像する。結像光学系205aの中には開口絞り206aが存在し、開口絞り206aを調整することで被写界深度が制御できる。撮像にあたっては、光源201aを点灯させ、プレパラート204に光を照射する。そして、照明光学系202a、プレパラート204、結像光学系205aを通って撮像面に結像した像を撮像部207aのイメージセンサで受光する。モノクログレースケール撮影時には光源201aを白色で露光し、1回撮像を行う。カラー撮影時には、RGBの3つの光源201aで順番に露光し、3回撮像を行うことで、カラー画像を取得する。撮像面で結像した検体の像は、複数のイメージセンサから構成される撮像部207aで光電変換され、A/D変換を実行した後、像は電気信号として画像処理装置110に送られる。本実施例ではA/D変換を実行した後のノイズ除去や色変換処理、鮮鋭化処理に代表される現像処理は画像処理装置110内部で行うとして説明する。しかし、現像処理は撮像部207aに接続された専用の画像処理ユニット不図示で行い、その後画像処理装置110にデータを送信することも可能であり、そのような形態での実施も本発明の範疇となる。プレ計測ユニット図2bは実施例1におけるプレ計測ユニット102の内部構成を示す図である。光源201bの光は照明光学系202bを通じ、光量ムラが無いように均一化され、ステージ203bの上に設置されたプレパラート204を照射する。プレパラート204を通った光は結像光学系205bにより、撮像部207b上の撮像面で結像する。結像光学系205bの中には開口絞り206bが存在し、開口絞り206bを調整することで被写界深度が制御できる。撮像にあたっては、光源201bを点灯させ、プレパラート204に光を照射する。そして、照明光学系202b、プレパラート204、結像光学系205bを通って撮像面に結像した像を撮像部207bのイメージセンサで受光する。モノクログレースケール撮影時には光源201bを白色で露光し、1回撮像を行う。カラー撮影時には、RGBの3つの光源201bで順番に露光し、3回撮像を行うことで、カラー画像を取得する。撮像面で結像した検体の像は、イメージセンサを備えた撮像部207bで光電変換され、A/D変換を実行した後、像は電気信号として画像処理装置110に送られる。本実施例ではA/D変換を実行した後のノイズ除去や色変換処理、鮮鋭化処理に代表される現像処理は画像処理装置110内部で行うとして説明する。しかし、現像処理は撮像部207bに接続された専用の画像処理ユニット不図示で行い、その後画像処理装置110にデータを送信することも可能であり、そのような形態での実施も本発明の範疇とする。画像処理装置図3は本発明における画像処理装置ホストコンピュータ110の内部構成を示す図である。CPU301はRAM302やROM303に格納されているプログラムやデータを用いて画像処理装置全体の制御を行う。またCPU301は、以降の実施例で説明する各種演算処理、データ処理、例えば、後述する被写界深度拡大処理、現像補正処理、合成処理、圧縮処理等を行う。RAM302は記憶装置113からロードされたプログラムやデータ、他のコンピュータシステム114からネットワークI/Fインターフェース304を介してダウンロードしたプログラムやデータを一時的に記憶するエリアを備える。またRAM302は、CPU301が各種の処理を行うために必要とするワークエリアを備える。ROM303はコンピュータの機能プログラムや設定データなどを記憶する。ディスプレイ制御装置306は画像や文字等をディスプレイ112に表示させるための制御処理を行う。ディスプレイ112はユーザに入力を求めるための画面表示を行うとともに、撮像装置120から取得しCPU301で処理した画像データを画像表示する。操作入力デバイス111はキーボードやマウスなどCPU301に各種の指示を入力することのできるデバイスにより構成される。ユーザは撮像装置120の動作を制御する情報を操作入力デバイス111により入力する。308は操作入力デバイス111を介して入力された各種の指示等をCPU301に通知するためのI/Oである。記憶装置113はハードディスクなどの大容量情報記憶装置であり、OSオペレーティングシステムや以降の実施例で説明する処理をCPU301に実行させるためのプログラムやバッチ処理によりスキャンした画像データなどを記憶する。記憶装置113への情報の書き込みや記憶装置113からの情報の読み出しはI/O310を介して行われる。制御I/F312は、撮像装置120を制御するためのコントローラ108と制御コマンド信号をやりとりするためのI/Fである。コントローラ108は本計測ユニット101およびプレ計測ユニット102を制御する機能を持つ。画像I/Fインターフェース313には上述以外のインターフェース、例えば、CMOSイメージセンサやCCDイメージセンサの出力データを取り込むための外部インターフェースが接続されている。なおインターフェースとしてはUSB、IEEE1394などのシリアルインターフェースやカメラリンクなどのインターフェースを使うことが出来る。この画像I/F313を通じて本計測ユニット101やプレ計測ユニット102が接続される。深度拡大処理フロー本実施例のバーチャルスライドシステムは、図4に示す深度拡大処理フローを実行し、ユーザが指定する被写界深度を持つ深度拡大画像データを生成する。以下、深度拡大処理について図4を用いて説明する。まず、深度拡大設定S401では、画像処理装置110に対し、操作入力デバイス111を通じてユーザが所望する被写界深度の拡大量つまり深度拡大画像データの被写界深度を指定する情報を、深度拡大原情報として設定する。深度拡大原情報には様々な設定を割り当てることが出来る。一例としては、対物レンズ倍率および/または開口絞りがある。その他には実在する光学顕微鏡の対物レンズの型番や仮想的な対物レンズの種類、被写界深度拡大後のNA開口数などがある。なお、NAを設定する場合、結像光学系205aのNAよりも小さくする必要がある。もちろん被写界深度そのものを指定することもできる。深度拡大設定S401では、画像処理装置110上の撮像装置120を制御するアプリケーションを通じて、画像処理装置110に接続されたディスプレイ112に設定画面を表示し、ユーザに前述の深度拡大原情報を設定させる。例えば、対物レンズ倍率を設定させる場合には、光学顕微鏡で良く使用される倍率、例えば40倍、20倍、10倍等を選択させる構成を取ると、ユーザにとって被写界深度が感覚的に把握しやすくなる。また開口絞りを設定させる場合、NAの変更率を設定させると良い。次に、プレ計測処理S402では、後段の本計測処理S403でのデータ取得に必要なXY分割撮像位置情報およびZスタック撮像位置情報を決定する。ここでXY分割撮像位置情報は標本を光軸Z軸に垂直な面XY平面に複数回に渡って分割撮像する際に必要な情報で、具体的には各分割撮像領域の中心を結像光学系205aと重ねるために必要なXYステージの位置や移動量を表す。Zスタック撮像位置情報は分割撮像領域毎に光軸Z軸に沿って移動してZスタック画像データを取得する際のZステージ移動に必要な位置情報を表す。なお、Zスタック撮像位置は分割撮像領域毎に異なってもよい。次に本計測処理S403では、前記プレ計測処理S402で決定した情報を用いて、コントローラ108を制御し、分割撮像領域毎にZスタック画像データを取得し、画像処理装置110に送信する。次に深度拡大処理S404では、前記本計測処理S403で撮影した分割撮像領域毎のZスタック画像データと深度拡大設定S401で設定した深度拡大原情報被写界深度の拡大量を用いて、被写界深度を拡大した画像データを生成する。上記のプレ計測処理S402、本計測処理S403、深度拡大処理S404に関しては、後に詳しく説明する。次に現像補正処理S405では、前述の被写界深度を拡大した画像データに対し、画像処理装置110が色変換処理、鮮鋭化処理、ノイズ低減処理などを実施する。これらの処理により、ディスプレイ112に表示される画像の色を標本の実際の色に近づけ、観察に好ましい鮮鋭化処理を施し、画像のノイズを抑制する。次に合成処理S406では、画像処理装置110が分割撮像領域毎に深度拡大処理および現像補正処理を行った画像データをつなぎ合わせ、プレパラート上の観察対象領域全体例えば20mm角の領域の像を作る。次に圧縮処理S407では、画像処理装置110が合成処理を行った画像データを圧縮し、データ容量を削減する。圧縮方法としては静止画圧縮符号化方式のJPEGやJPEG2000などを用いることが可能である。続いて、伝送処理S408では、画像処理装置110が圧縮処理した画像データを記憶装置113に送り、データを保存する。画像処理装置110はネットワークI/Fを介してネットワーク上のコンピュータシステム114や画像サーバに画像データを送信してもよい。プレ計測及び本計測における装置の動作図5a~図5d、図6a、図6bを用いてプレ計測処理S402及び本計測処理S403における装置の動作について説明する。図5aはプレパラートの模式図である。スライドグラス510上にはカバーグラス511によって検体500を封入した領域とラベルエリア512がある。本計測処理S403では、プレ計測で求めた検体500が存在する外接矩形領域の内部を撮像対象とする。図5bは本計測処理S403で検体500が存在する外接矩形領域を複数領域に分割して撮像する様子を表している。図5cは撮像面を表しており、有効視野502は本計測ユニット101の結像光学系205aを通して像が映る範囲を、センサ有効領域503は撮像部207aのイメージセンサで撮像可能な領域を示している。物体面における撮像領域501斜線領域は本計測ユニット101の結像光学系205aを通じて結像し、撮像面における撮像領域504に対応する。図5cに示すように、センサ有効領域503には、撮像領域504に対して若干広い領域が割り当てられている。これは、結像光学系205aが有する光学的な収差やイメージセンサの取付位置のズレを許容するためのマージンである。すなわち、光学的な収差やセンサの取付位置のズレがあったとしても、物体面の撮像領域501がセンサ有効領域503内に収まるようにしてある。合成処理S406では、センサ有効領域503の画像に対して収差補正や位置ズレ補正を行い、補正した画像から撮像領域504に対応する部分を切り出し、その画像を合成に用いる。図5dはプレ計測におけるプレパラート204の撮像領域505を示している。プレ計測は本計測を精度よく撮像するための撮像制御情報取得を目的としており、画像の特徴が大まかに把握できれば良く、本計測ほどの解像度倍率は必要ない。本実施例のプレ計測では低倍率でプレパラート204の全体を撮像する。また、本計測と異なり複数の領域に分割することなく1つのイメージセンサでプレパラート204の全体を一括撮像する。これによりプレ計測ユニット102の構成を簡易化することができるとともに、プレ計測に要する時間、ひいてはプレ計測と本計測を含む撮像処理全体の時間を短縮することができる。図6aは図5bで示す分割撮像領域をラスタースキャン順に撮像する場合のステージ203aのXY方向の移動方向と順序を表している。プレパラートの左上から右下の領域を撮像するには、プレパラートを搭載するステージ203aは逆方向の右下から左上に移動する。このように撮像領域を複数領域に分割し、ステージ203aを移動しながら繰り返し撮像することで比較的小さなイメージセンサでも広い領域を撮像することができる。図6bは本計測処理S403で、焦点位置観察深さ、フォーカス位置、或いはピント位置ともいうの異なる複数の画像を取得する場合のステージ203aのZ方向奥行き方向の移動方向を示している。図に示すように、プレパラート204内の検体の上側カバーグラス裏面側に焦点位置を移動するには、結像光学系205aの光軸方向に沿って、ステージ203aをZ方向下向きに移動させる。逆に、検体の下側スライドグラス表面側に焦点位置を移動するには、ステージ203aをZ方向上向きに移動させる。撮像部207a、或いは、撮像部207aとステージ203aの両方を結像光学系205aの光軸方向に沿って移動することでも焦点位置を変更することができる。また、結像光学系205aのレンズを制御することにより光学的に焦点距離を変更することでも、焦点位置を変更可能である。なお、プレ計測ユニット102のステージ機構は本計測ユニット101とほぼ同じため、説明は省略する。プレ計測処理S402続いて、図7を用いてプレ計測処理S402の内部処理フローについて説明する。プレ計測処理ではまずプレ計測データ取得S701で標本の画像プレ計測データと呼ぶを取得する。取得時の装置の動作については既に述べたため、省略する。続いてXY分割撮像位置決定処理S702では、標本のXY方向の撮像領域を決定し、XY分割撮像位置情報を計算する。具体的にはプレ計測データの画像から検体500を検出し、検体500の外接矩形を計算し、標本の撮像範囲を決定する。この矩形範囲を本計測ユニット101で1回で撮像出来ない場合には、図5bで示すように、矩形範囲を分割し本計測の分割撮像領域501を決定する。その後、求めた分割撮像領域毎に撮影時のXYステージ移動に必要な情報、即ちXY分割撮像位置情報を計算する。なお、外接矩形の計算は画像輝度の分析等、公知の技術で実現できるため、詳細説明は省く。またユーザから撮像領域の指定がある場合には、その領域に対して分割撮像領域の決定とXY分割撮像位置情報の計算を行う。続いて、AF合焦位置取得処理S703では、AFユニット109を用いて分割撮像領域毎にオートフォーカスを実施し、フォーカスが合うZ位置AF合焦位置と呼ぶを取得する。続いて、Zスタック撮像位置決定処理S704では本計測ユニット101を用いて分割撮像領域毎にZスタック撮像位置を決定する。この処理では、各分割撮像領域に対し、AF合焦位置から光軸方向に沿って前後にステージ203aを動かし、画像内に一定以上の鮮鋭度コントラストを持つZ方向の範囲を求め、その範囲をZスタック撮像領域に決定する。なおS704では本計測ユニット101の代わりにプレ計測ユニット102を用いても良いし、変位計103を用いて、プレパラート内の検体の表面と裏面位置を調べ、その間をZスタック撮像領域に定めても良い。その後、AF合焦位置を中心として所定のZスタック間隔で動かすことにより、求めたZスタック撮像領域がカバー出来るよう、レイヤー画像の枚数と撮影時のZステージの移動に必要な情報、即ち、Zスタック撮像位置情報を求める。なお、Zスタック間隔は結像光学系205aの被写界深度が好ましいが、ユーザが事前に設定した値を用いることもできる。例えば、Zスタック撮像領域はAF合焦位置を中心に±5μmの領域で、結像光学系205aの被写界深度は0.5μmとする。その場合、Zスタック画像を構成するレイヤー画像の枚数は21枚となり、Zスタック撮像位置情報は、AF合焦位置±0.5n[μm]n=0~10となる。このように分割撮像領域毎にZスタック撮像位置情報を決めるのは、標本の凹凸やカバーグラスのうねり等の影響により、観察対象となる検体細胞や組織のZ方向の位置がばらついているからである。上記のXY分割撮像位置情報およびZスタック撮像位置情報は、画像処理装置110で撮影時にステージ203aの移動量が少なく済むように、昇順または降順に位置を並べ替え、画像処理装置110内のRAM302や記憶装置113に格納する。本計測処理S403本計測処理S403では後段の深度拡大処理S404での被写界深度拡大処理に必要となる分割領域毎のZスタック画像の取得を行う。Zスタック画像とは、光軸方向に焦点位置を変えながら被写体プレパラートを撮像することにより得られる画像群データであり、1つのZスタック画像は各焦点位置で得られた複数枚のレイヤー画像データで構成される。本計測処理S403の内部処理について図8を用いて説明する。まずS801では、画像処理装置110からプレ計測処理S402で求めたXY分割撮像位置情報の初期位置を読み出し、コントローラ108が、ステージ制御105を介して、初期の分割撮像領域が撮影できる位置にステージ203aを移動させる。次にS802では、画像処理装置110から初期の分割撮像領域におけるZスタック撮像位置情報を読み出し、同様にコントローラ108が、ステージ制御105を介してステージ203aを初期のZ位置へ移動する。例えば前述の例であれば、AF合焦位置0.5×10[μm]の位置に焦点位置を合わせるようにステージ203aが移動する。続いてS803で画像データの取得を行う。取得した画像データは本計測ユニット101から画像処理装置110に転送され、画像処理装置110内のRAM302あるいは記憶装置113に格納される。続いてS804では、現在の分割撮像領域においてZスタック画像の全てのレイヤー画像のデータを取得したかを判定し、未取得のデータが存在すれば、S805に進み、ステージ203aを次のZ位置に移動する。例えば前述の例であれば、焦点位置を0.5μmずつ小さくしていき、全部で21枚のレイヤー画像を取得する。S804で全て取得済みであれば、次のS806に進む。続いてS806では、全ての分割撮像領域のZスタック画像を取得したかを判定し、未取得のデータが存在すれば、S807に進み、ステージ203aを次の分割撮像領域に移動させる。全て取得済みであれば、処理を終了する。以上の処理を繰り返し、本計測処理S403では、各分割撮像領域のZスタック画像のデータを取得する。深度拡大処理S404深度拡大処理S404では、分割撮像領域毎に、本計測処理S403で取得したZスタック画像を構成する複数のレイヤー画像データから、このレイヤー画像データよりも被写界深度が深い深度拡大画像データを生成する処理を行う。以下、撮像により得られたレイヤー画像データを第1の画像データ、深度拡大処理により生成される画像データを第2の画像データとも呼ぶ。本実施例では説明の簡略化のため、分割撮像領域が1つで検体500が、本計測撮像領域501よりも小さい場合、本計測処理S403でモノクログレースケール撮像する場合について説明する。ただし、複数の分割撮像領域が存在する場合でも、以降で説明する深度拡大処理S404を分割撮像領域毎に繰り返すことで、検体500全体の被写界深度の拡大が可能である。またカラー撮像においてもRGB毎に被写界深度の拡大が可能である。図9aを用いて深度拡大処理S404について説明する。深度拡大処理S404では、画像処理装置110は、まず、深度拡大特性情報取得処理S901において深度拡大設定S401で設定された深度拡大原情報から、被写界深度と結像特性情報を取得する。本処理については後に詳しく説明する。次に合成焦点面対象面決定処理S902では、画像処理装置110は、後段の深度合成データ生成処理S903で用いる、画像合成時の焦点面、即ち合成焦点面と、画像合成時の処理対象面、即ち合成対象面を決定する。合成焦点面対象面決定処理S902は本実施例の深度拡大アルゴリズムにおける特徴的な処理の一つであり、後段の深度合成データ生成処理903での被写界深度の拡大処理の計算負荷を軽減する。本処理については後に詳しく説明する。次に深度合成データ生成処理S903では、前段の合成焦点面対象面決定処理S902で決定した合成焦点面および合成対象面を用いて、被写界深度を拡大したデータを生成する。本処理の内部で用いる被写界深度の拡大処理としては様々な方式が考えられるが、本実施例では後に3つの手法に絞って説明する。しかし説明した手法に限定されないため、それ以外の手法を用いる場合でも本発明の範疇となる。深度拡大特性情報取得処理S901図9bは深度拡大特性情報取得処理S901の内部処理フローである。被写界深度取得処理S1001では、画像処理装置110が、深度拡大設定S401で設定した深度拡大原情報から被写界深度を取得する。例えば、深度拡大原情報として対物レンズ倍率を設定する場合、画像処理装置110は、画像処理装置110内にある対物レンズ倍率と被写界深度の関係を記録する対応表を用いて、対物レンズ倍率から被写界深度を求めると良い。対応表には事前の実験または計算により求めた対物レンズ倍率と被写界深度の対応関係を記載しておく。深度拡大原情報として実在する光学顕微鏡の対物レンズの型番や仮想的な対物レンズの種類を設定する場合も同様に対応表を用いて被写界深度を求められる。深度拡大原情報として、対物レンズ倍率および開口絞りを設定する場合には、画像処理装置110は、対応表を用いて対物レンズ倍率に対応するNAを取得し、開口絞りによるNAの変更率を掛け算してトータルでのNAを求める。そして、下記の<1>式により、被写界深度Dを求めると良い。±D=nλ/{2NA2}...<1>式n:媒質の屈折率、λ:光の波長[μm]、NA:対物レンズのNA開口数屈折率と光の波長はそれぞれ一般的な標本の屈折率、基準波長0.55μmを用いることで大まかに計算できる。続いて結像特性取得処理S1002では、画像処理装置110は、深度拡大原情報を用いて結像特性情報を求める。結像特性情報とは、深度拡大原情報で指定された被写界深度拡大後の被写界深度をもつ仮想的な光学系の結像特性を表す情報である。具体的には、焦点位置からの距離に従って変化する、複数の2次元PSFPointSpreadFunctionの情報や、3次元PSFの情報を結像特性情報として用いることができる。PSFは周波数変換してOTFOpticalTransferFunctionの形で与えても構わない。これらの情報は深度拡大原情報と対応付けて、前述の対応表に格納しておくことで実現できる。また深度拡大原情報から求めた値例えばNAから数式で被写界深度拡大後の結像特性である3次元PSFや3次元OTFを計算しても良い。合成焦点面対象面決定処理S902図9cを用いて合成焦点面対象面決定処理S902の内部処理フローを説明する。合成焦点面対象面決定処理S902では、拡大後の被写界深度に基づいて深度拡大画像の焦点位置に相当する合成焦点面を決定する処理と、Zスタック画像を構成する全レイヤー画像の中から深度拡大画像の生成に使う画像を合成対象面として選択する処理を行う。まず、合成焦点面の決定処理について説明する。ここでは、合成焦点面を中心とする被写界深度拡大後の被写界深度の内に観察対象としての物体が最も多く含まれるように、合成焦点面のZ方向の位置を決定する。画像処理装置110は、まず画像評価値計算処理S1101でZスタック画像を構成する各レイヤー画像の画像評価値を計算する。画像評価値は、レイヤー画像内に焦点があった領域がどの程度存在するかを分析する評価値であり、画像に対する鮮鋭度評価値の面内積分値が好適である。鮮鋭度評価値としては、フーリエ変換やウェーブレット変換等の画像の周波数解析を含むコントラスト評価値、ラプラシアンフィルタに代表されるエッジ検出フィルタによる演算値、局所領域内の画素値の分散を画像全体に渡って求めた分散値等が挙げられる。図10a及び図10bは画像評価値の計算の様子を示す図である。図10aはZ方向に焦点位置を変えて撮像装置で撮像した標本のZスタック画像を示す。このZスタック画像は16枚のレイヤー画像で構成されている。Z軸に記載の番号はレイヤー画像の位置を識別するためのインデックスであり、単位は持たない。図10bは各レイヤー画像に対して計算した画像評価値を示す。横軸は画像評価値、縦軸は前述のレイヤー画像の位置を表すインデックスを示す。図10aのレイヤー画像1201aの画像評価値は図10bの1201bの値に対応している。図10bのグラフは、標本内でインデックス2および0のZ位置に、焦点が合う物体が最も多く在在することを示している。次の画像評価積分値計算処理S1102では、画像処理装置110は、深度拡大処理S404において被写界深度を拡大した後の画像における画像評価値を予測する処理を行う。図11a~図11cを用いて画像評価積分値計算処理S1102の処理を説明する。図11a、b、cはそれぞれ被写界深度が前後1枚、2枚、3枚それぞれ合計3枚、5枚、7枚の場合における画像評価積分値を示すグラフである。これらのグラフは、図10bの各インデックスにおいて、各インデックスを中心として被写界深度の枚数分だけ隣接するインデックスの画像評価値を積分した結果であり、被写界深度拡大後の各Z位置における画像評価値の推定値となっている。上記の処理により、撮像時の対物レンズ倍率や開口絞りが被写界深度拡大後のそれらと異なる場合でも、被写界深度拡大後に最もピントが合うZ位置を予測できる。図11a~図11cから分かるように、被写界深度が変わると、ピントが合う物体が最も多く含まれるZ位置画像評価値が最も大きくなるZ位置も変わり得る。なお、被写界深度拡大後の被写界深度が、撮像時のZスタック間隔レイヤー画像の間隔の整数倍で表現できない場合には、図10bの各Z位置の画像評価値を補間して積分区間内の画像評価値を求め、それらを積分することにより画像評価積分値を求めると良い。なお、補間方式としては公知の線形補間やスプライン補間等が適用可能である。次に合成焦点面決定処理S1103の処理を説明する。本処理では前段の画像評価積分値計算処理S1102で求めた画像評価積分値が最大となるZ位置を求め、そのZ位置を被写界深度拡大後の焦点面に設定する。図11aではインデックス1で表されるZ位置で画像評価積分値が最大になっており、そのZ位置が被写界深度拡大後に最も焦点が合う領域と推測される。そのため、前後1枚合計3枚の範囲まで被写界深度を拡大する場合には、合成焦点面はインデックス1で表されるZ位置に決定する。同様に、前後2枚合計5枚の範囲まで被写界深度を拡大する場合にはインデックス2のZ位置、前後3枚合計7枚の範囲まで被写界深度を拡大する場合にはインデックス3のZ位置がそれぞれ合成焦点面に選ばれる。次に合成対象面決定処理S1104の処理を説明する。合成対象面決定処理S1104では、画像処理装置110が、前述の結像特性情報に基づいて、焦点位置合成焦点面から光軸方向に離れた物体が画像に対して与える影響が所定の条件よりも大きくなる範囲これを影響範囲と呼ぶを求める。そして、Zスタック画像を構成するレイヤー画像の中から、影響範囲内に焦点位置が含まれるレイヤー画像を合成対象面として選択する。選択された合成対象面は、後段の深度合成データ生成処理S903での被写界深度拡大処理に用いられる。合成対象面の選択基準となる、所定の条件及び影響範囲の決め方の一例について、図12a及び図12bを用いて説明する。図12a及び図12bは結像特性取得処理S1002で求めた、被写界深度拡大後つまり仮想的な光学系の結像特性PSFを示す図である。図12aの円錐は、焦点位置図の原点から離れる程、PSFが拡がる様子を示している。図12bは図12aで示した3次元PSFの、原点からの距離が異なる2つの平面Z=aおよびZ=ba<bにおける2次元PSF横軸:r、縦軸:強度を示している。Z=aは被写界深度内の平面の2次元PSFを示し、Z=bは被写界深度外の平面の2次元PSFを示している。Z=bにおける原点中心画素の像強度はZ=aにおけるそれと比べてかなり小さくなっている。このようなPSF形状から予測できるように、被写界深度拡大後の焦点面における画像において、|Z|>bの位置にある物体は大きくぼけて重畳されるのみで、焦点面での画質に大きな影響を及ぼさないことが推定される。ここで焦点位置から離れた位置の物体の、焦点位置に対する画質影響評価指標の例として、焦点位置および焦点位置から離れた位置のPSFの中心イメージセンサの中心画素の強度比を定義する。そして、中心の強度比が1/N以下なら焦点面の画質への影響が無視できるとする。Nには任意の値正の実数を設定することができるが、被写界深度拡大後の画像の画質評価実験やシミュレーションなどから適切な値を求めると良い。本実施例では、N=10に設定する。図12aにおいてPSFの中心の強度比が1/Nとなる、焦点位置からの距離をZ=bとすると、少なくとも|Z|≦bの範囲内に存在する画像を用いれば、合成焦点面の画像が精度よく再現できる。従って、|Z|≦bの範囲内に存在するレイヤー画像を合成対象面と定義し、深度合成データ生成処理S903に入力すれば、画質を維持した上で、合成枚数の削減による演算コストの低減が実現できる。この例では、PSFの中心の強度比が1/Nが前述の所定の条件に対応し、|Z|≦bの範囲が影響範囲に対応する。上記のように、本実施例では、合成対象面決定処理S1104において、深度拡大像の画質に与える寄与を考慮し、深度拡大処理に用いるレイヤー画像の枚数を削減する。これにより、従来の方法Zスタック画像を構成する全てのレイヤー画像を用いる方法に比べて、後段の深度合成データ生成処理S903における処理の高速化が実現できる。なお、焦点位置における画質影響評価指標は前記の方法に限らない。例えば、焦点位置から離れた位置の、指定周波数におけるMTFModulationTransferFunctionを評価指標と定め、10以上となる範囲を合成対象面と定めても良い。またMTFと人間の視覚特性の積の積分値を評価指標とし、人間の知覚特性に基づく閾値から合成対象面を定めても良い。その他にも様々な変形例が考えられる。従って、前記以外の評価指標を適用し、合成焦点面を定める場合も本発明の範疇となる。また、画質と速度のバランスに応じて、前記で定めた合成対象面から一部の画像を間引くことでさらなる処理枚数の低減と高速化を実現することも可能である。なお、前述の画質影響評価指標を用いた合成対象面決定処理は必ずしも必須ではない。すなわち、合成焦点面を定めた後、合成対象面決定処理の代わりに、合成対象面として分割領域毎のZスタック画像を構成するレイヤー画像の全てを設定したり、合成焦点面を中心として予め定めた一定範囲のレイヤー画像を一律に設定することでも実現できる。この方法は前述の画質影響評価指標を用いて合成対象面を選択する方法に比べ、処理枚数が増加する可能性があるが、条件判定による分岐処理を減らし、回路構成の簡素化や安定動作を実現出来るメリットがある。いずれの方法においても合成焦点面を定めることにより、後段の深度合成データ生成性処理S903の計算負荷を低減できるメリットが実現できる。図13は合成焦点面対象面決定処理S902で決定した情報を説明する模式図である。1501は合成焦点面、1502は合成対象面、1503は拡大後の被写界深度を表している。前述の数値と対応づけると、1502はb=5合計11枚、1503はa=1合計3枚となっている。深度合成データ生成処理S903図14は深度合成データ生成処理S903の内部処理を示すブロック図である。画像処理装置110は、画像取得処理1601において、合成焦点面の情報と合成対象面のレイヤー画像のデータを、画像処理装置110内のRAM301や記憶装置113から取得する。次に画像処理装置110は、画像合成処理1602において、合成対象面のレイヤー画像データに対し処理を施し、被写界深度を拡大した画像データを出力する。画像合成処理1602の内部処理として様々な方式を用いることが出来る。例えば、1評価値に基づく選択統合手法、2重み付け合成手法、33次元フィルタリング手法などが存在する。本実施例では前記3つの手法を用いる場合に絞って説明するが、それ以外の被写界深度拡大方法を用いた場合も本発明の範疇に含まれるとする。1評価値に基づく選択統合手法図15aを用いて1評価値に基づく選択統合手法の例について説明する。この手法1は、拡大後の被写界深度に対応する仮想的な光学系の結像特性に基づくコンボリューション処理を適用することによって合成対象面の各々に対し焦点位置合成焦点面からの距離に応じたぼけを付与する手法である。まずコンボリューション処理S1701では、画像処理装置110は、結像特性取得処理S1002で求めた結像特性情報、つまり合成焦点面を原点としたときの当該Z位置における被写界深度拡大後の2次元PSFであるwnx,yを取得する。続いて、画像処理装置110は、合成対象面のレイヤー画像gnx,yに対し、前記のwnx,yをコンボリューションし、各レイヤー画像に合成焦点面からの距離に応じたぼけを与える。数式で表現すると以下の<2>式で表現できる。なお、<2>式ではコンボリューション演算を表し、nはレイヤー画像の位置を識別するためのインデックスを表す。これにより、各々のZ位置に焦点を合わせて撮像された各Z位置のレイヤー画像は、合成焦点面からの距離に応じて、ユーザの設定する被写界深度拡大量に対応するぼけを持つ画像に補正される。次に評価値計算処理S1702では、画像処理装置110は、コンボリューション処理を適用したレイヤー画像のそれぞれに対し、画像評価値を計算する。画像評価値は画像内の鮮鋭度の高い位置、即ち、焦点の合った位置の抽出に適した値であり、合成焦点面対象面決定処理S902で述べた鮮鋭度評価値が活用できる。画像評価値は、レイヤー画像内の各画素の水平XY位置について計算される。次に選択処理S1703では、画像処理装置110は、合成対象面のレイヤー画像の間で画像評価値を比較し、画像評価値が最も高くなるZ位置を特定する。画像評価値の比較は画像内の同一のXY位置について行われ、XY位置毎に画像評価値が最も高くなるZ位置が選択される。その結果は、XY位置と選択されたZ位置との対応を表すZ位置マップとして出力される。このとき、各XY位置において最も鮮鋭度評価値が高いZ位置を抽出すると、ノイズ等の影響により画像面内でZ位置が急激に変化し、合成画像の画質を低下させる可能性がある。そこで例えば、Z位置マップに対して空間フィルタノイズ抑制フィルタを適用し、Z位置の急激な変化を抑えることが好ましい。ノイズ抑制フィルタとしては、メディアンフィルタやバイラテラルフィルタ等を用いることができる。次に統合処理S1704では、画像処理装置110が、XY位置毎に、合成対象面のレイヤー画像群からZ位置マップに記載されたZ位置の画素値を抽出し、それらを統合した合成画像を生成する。なお、ここでは画素単位で処理を行ったが、所定サイズの画素ブロック単位で同様の処理を行うことで計算量を削減することも可能である。2重み付け合成手法次に図15bを用いて2重み付け合成手法の例について説明する。手法2は、仮想的な光学系の結像特性に基づいて決定される重みを用いて、周波数空間上で合成対象面のレイヤー画像を重み付け加算する手法である。具体的には、合成対象面となる複数のレイヤー画像の各々を2次元離散フーリエ変換し、それらを被写界深度拡大後の焦点距離に応じて重み付け加算した後、2次元離散逆フーリエ変換することで合成画像深度拡大画像を生成する。なお、以降の説明では、離散フーリエ変換をDFT、離散逆フーリエ変換をIDFTと記載する。まず、周波数変換処理S1801では、画像処理装置110が、各Z位置nのレイヤー画像gnx,yに対し、2次元DFTを行う。その結果、合成対象面の枚数分の2次元DFT画像Gnu,vが生成される。次に、重み付け合成処理S1802では、画像処理装置110が、各2次元DFT画像Gnu,vに対し、空間周波数成分毎に係数を持つ重み付け関数Wnu,vを掛け算した後、全画像を合計することで合成周波数データAu,vを生成する。数式で表現すると、下記の<3>式のようになる。Ns:合成対象面の画像の枚数被写界深度の変更はこの空間周波数成分毎に係数を持つ重み付け関数Wnu,vを変更することで実現できる。一般に焦点位置から離れるほど画像のMTFは低下するため、重み付け関数Wnu,vは合成焦点面からの距離に応じて小さくなるように設定するとよい。重み付け関数Wnu,vには様々な関数が適用できるが、例えば、結像特性取得処理S1002で取得した焦点位置からの距離に従って変化する2次元OTFを用いて重み付け関数Wnu,vを決定することが可能である。最後に逆周波数変換処理S1803では、画像処理装置110は、合成周波数データAu,vを2次元IDFTし、被写界深度を拡大した像ax,yを得る。33次元フィルタリング手法次に図15cを用いて33次元フィルタリング手法の例について説明する。なお、ここで3次元フィルタリング手法と呼ぶ手法は、特許文献2に記載される手法である。まず本手法では、画像処理装置110が、合成対象面の各レイヤー画像gnx,yに座標変換処理S1901を施し、3次元フィルタ処理で被写界深度拡大が可能な形式gX,Y,Zに変換する。なお、本処理はZ位置により像倍率が変わらない場合には不要で、その場合にはS1901の処理は省略する。次に3次元フィルタ計算S1902では、画像処理装置110が、被写界深度拡大に必要な3次元フィルタを計算する。結像光学系205aの撮像時の3次元OTFをHu,v,w,被写界深度を拡大した撮影像の3次元OTFをHau,v,wとすると、3次元フィルタはHau,v,wHu,v,w-1で計算できる。このフィルタは、レイヤー画像を撮像した結像光学系205aの結像特性によるぼけをキャンセルした上で、仮想的な光学系の結像特性によるぼけを付与する機能をもつフィルタである。結像光学系205aの撮像時の3次元OTFであるHu,v,wに関しては、予め求めた値を画像処理装置110内のROM303や記憶装置113に格納しておき、画像合成処理1602でその値を読み出して用いる。被写界深度を拡大した撮影像の3次元OTFであるHau,v,wには結像特性取得処理S1002で取得した値を設定する。次に、3次元フィルタ処理S1903では、画像処理装置110が、gX,Y,Zに対し、3次元DFTを実施し、Gu,v,wを計算する。その後、3次元フィルタ計算S1902で求めた3次元フィルタHau,v,wHu,v,w-1と積算し、被写界深度を拡大した3次元DFT画像Au,v,wを求める。次に、焦点面データ算出処理S1904では、画像処理装置110が、原理的には、Au,v,wに対し3次元IDFTを実施することで、座標変換空間における被写界深度拡大像aX,Y,Zを求めることが出来る。F-1{}は3次元IDFTを示す。しかし、上記のように、3次元IDFTを変数u,v,wそれぞれX,Y,Zの周波数変換後の変数に適用し、その後、合成焦点面座標変換前z=p、座標変換後Z=Pの画像aX,Y,Pを抽出すると演算コストが増加し、計算時間がかかる。3次元DFTIDFTは、それぞれの変数u,v,wに対し、順次1次元DFTIDFTを行うことで実現できる。そのため、Au,v,wの変数wの全要素に1次元IDFTを行うのではなく、合成焦点面Z=Pのみの1次元IDFT値を求め、変数u,vに対しそれぞれ1次元IDFTを実施してaX,Y,Pを求めることで演算コストを抑えることができる。最後に座標変換処理S1905おいて、画像処理装置110が合成焦点面の画像データに対して座標変換処理を行い、ax,y,pを求める。なお、S1905においてもS1901同様、空間変換処理が必要ない場合には本処理は省略できる。以上のように深度合成データ生成処理S903で適用可能な3つの手法について述べた。いずれの手法においても、前段の合成焦点面対象面決定処理S903において合成焦点面および合成対象面を定めることで、演算コストを抑え、計算時間を短縮できることが分かる。予め合成焦点面を定めない場合、想定される複数の焦点面に対し、被写界深度拡大処理を行う必要があり、演算コストがかかることが課題である。本実施例の方法では比較的良く標本の内部構造を表している被写界深度の浅い画像データを用いて合成焦点面を予測することにより、被写界深度拡大処理の演算コストを大幅に低減し、撮像装置のスループットや応答速度を高めることが出来る。なお、深度合成データ生成処理S903において、2および3の手法で用いるDFTは計算負荷がかかるため、本計測処理S403と深度拡大処理S404をパイプライン的に処理しても良い。例えば、本計測処理S403において、次のZ位置の画像データを取得している間に、取得済みの画像データに対して2次元DFTを実行すると良い。さらに次の分割撮像領域のZスタック画像取得に移った後、取得済みの分割撮像領域のZスタック画像に対し、深度拡大処理S404を実行すると良い。上記の並列処理により、深度拡大処理S404の深度合成データ生成処理S903における計算時間を短縮することが出来、応答速度を高めることが出来る。本実施例では本計測処理S403でモノクロ撮像する場合について説明したが、カラー撮像する場合について補足する。カラー撮像の場合、合成焦点面および合成対象面はRGB共通に設定するのが好ましい。例えば、いずれか1つのチャネルGの画像を用いて合成焦点面および合成対象面を決定し、それらをRGBで共通に用いると良い。結像特性情報PSFはRGB共通でもRGB毎に変えても良い。共通に設定する場合には、例えば、いずれか1つのチャネルGのPSFをRGBで共通に用いると良い。結像特性情報PSFをRGB毎に変える場合には、深度拡大特性情報取得処理S901でRGBそれぞれの結像特性情報を求め、それらを用いてRGB毎に深度合成データ生成処理S903を実施する。またカラー撮像はRGB3色に限らない。光源201aが複数波長の光源またはカラーフィルタを備える場合、複数波長での撮像マルチバンド撮像が可能となり、RGBと同様に複数波長毎の深度合成データ生成処理が可能となる。本実施例の構成によれば、高い解像度を持つ結像光学系205aを用いて撮像した像をユーザが所望する設定値に応じて被写界深度拡大した像に高速に補正出力することが可能となり、ユーザ病理医の標本の診断効率を高められる効果が実現できる。[実施例2]本実施例では、本計測処理S403における分割撮像領域間で合成焦点面のZ位置が大きく異なる場合の被写界深度の拡大方法について述べる。本実施例で述べる構成により、プレパラートのうねりや凹凸、標本内の組織や細胞の分布位置のばらつきが大きい場合であっても、合焦位置では焦点ぼけの無い明瞭な像を生成してユーザ病理医の診断効率を高められる効果が実現できる。図16aは本実施例における深度拡大処理S404の内部処理フローを示す。深度拡大特性情報取得処理S2001では、画像処理装置110は、S901と同様の方法で深度拡大処理後の被写界深度および結像特性情報を取得する。これらの処理は既に述べた為、詳細は省略する。続いて領域内サブ分割処理S2002では、画像処理装置110は、本計測処理S403で定めた分割撮像領域に対し、さらに領域内のサブ分割処理を行う。続いてS2003では、画像処理装置110が、S2002で作成された複数のサブ分割領域の中から最初の処理対象となるサブ分割領域を選ぶ。画像処理装置110は、処理対象のサブ分割領域に対し、合成焦点面と合成対象面を決定しS2004、被写界深度を拡大した深度合成データを生成するS2005。合成焦点面対象面決定処理S2004および深度合成データ生成処理S2005の処理は実施例1のS902およびS903と同様のため、説明は省略する。ループ処理S2006,S2007により、領域内サブ分割処理S2002で定めた全てのサブ分割領域に対し、S2004とS2005の処理が実行される。続いて図16b、図17aおよび図17bを用いて領域内サブ分割処理S2002の詳細を説明する。初期分割処理S2101では、画像処理装置110は予め定めた所定の領域サイズで領域内のサブ分割を行う。例えば、本計測撮像領域504の画素数を2048×2048画素とし、サブ分割領域サイズは512×512画素とすると、この処理により16のサブ分割領域に分割できる。続いて、画像評価値計算処理S2102では、画像処理装置110は、各レイヤー画像の各サブ分割領域16個に対する画像評価値を計算する。続いて、画像評価積分値計算処理S2103では、画像処理装置110は、S2102で求めた画像評価値から、各サブ分割領域のZ位置毎に画像評価積分値を計算する。続いて、合成焦点面決定処理S2104では、画像処理装置110は、サブ分割領域毎に画像評価積分値が最大となるZ位置を求め、これを被写界深度拡大後の焦点面である合成焦点面に決定する。S2102、S2103、S2104の処理は、それぞれ、実施例1で述べた図9cのS1101、S1102、S1103の処理に対応している。続いて、隣接領域差計算処理S2105では、画像処理装置110は、隣接するサブ分割領域の合成焦点面のZ位置の差を計算する。続いて、閾値判定S2106では、画像処理装置110は、全てのサブ分割領域において隣接する合成焦点面のZ位置の差が閾値未満ならば、領域分割処理を終了する。一つでもZ位置の差が閾値以上であれば再分割処理S2107に進み、さらに細かな分割領域サイズ例えば、512×512画素の半分の256×256画素で分割し、画像評価値計算処理S2102からの処理を続行する。なお、閾値判定S2106で用いる閾値の一例としては、深度拡大処理後の被写界深度がある。被写界深度を拡大した後でもサブ分割領域間で合成焦点面の差が被写界深度内に収まっていれば、サブ領域間の繋ぎ目を目立たなくできる。または閾値として事前に設定した一定の距離値を用いることも可能である。再分割処理S2107において新たに適用するサブ分割領域サイズとして様々な値が設定可能であるが、一例として現在の分割領域サイズの半分とする。例えば、512×512画素の分割で閾値以上と判定された場合、256×256画素で再分割する。その後、128×128画素、64×64画素、32×32画素、のように、隣接するサブ分割領域間での合成焦点面のZ位置の差が閾値未満となるまで再分割が繰り返される。図17aおよび図17bは領域内サブ分割処理S2002でのサブ分割処理を説明する1次元の模式図である。図17aの領域1、2はそれぞれ初期分割処理S2101で本計測撮像領域504を分割したサブ分割領域を示している。また点線22は、標本内の物体組織や細胞の分布を模式的に示している。図17aにおいて領域1、2の隣接する合成焦点面の差はaとなり、差aは被写界深度拡大後の被写界深度tよりも大きい。このtを閾値とした場合、領域1、2はそれぞれ図17bに示すように半分のサブ分割領域に再分割される。領域1、2を再分割した領域1-1、1-2、2-1、2-2において隣接する合成焦点面の差はそれぞれa1、a2、a3となり、いずれも閾値である被写界深度t未満に収まっている。本実施例では1つの分割撮像領域の中を再分割する方法について説明したが、隣接する分割撮像領域の間で合成焦点面の差が閾値t以上となる場合も想定される。そこで、領域内サブ分割処理S2002内で隣接する分割撮像領域の合成焦点面の差を計算し、閾値t以上ならば、それらの分割撮像領域を再分割する処理を追加すると良い。本実施例によれば、領域内サブ分割処理S2002を導入することにより、プレパラートのうねりや凹凸、標本内の組織や細胞の分布位置のばらつきにかかわらず、標本の全域で焦点が合い、かつ、繋ぎ目が滑らかな被写界深度拡大像を生成することが出来る。[実施例3]本実施例では深度拡大設定S401において、複数の被写界深度に対応する深度拡大原情報を指定可能とし、複数の被写界深度拡大像を生成する処理について述べる。被写界深度が異なる複数枚の被写界深度拡大像を生成することにより、光学顕微鏡を操作して得られる画像を一度に生成でき、バーチャルスライドを用いたユーザの診断効率を高める効果が実現できる。図9aを用いて処理の流れを説明する。まず、深度拡大特性情報取得処理S901において、画像処理装置110が、深度拡大設定S401において設定した深度拡大原情報から複数の被写界深度と結像特性情報を取得する。次に合成焦点面対象面決定処理S902では、画像処理装置110が、深度拡大原情報毎に合成焦点面と合成対象面を設定する。次に深度合成データ生成処理S903では、画像処理装置110が、S901で求めた複数の被写界深度と結像特性情報を用いて、複数の被写界深度拡大画像を生成する。S901~903のいずれの処理も複数の深度拡大原情報毎に計算する以外には実施例1、2との違いはないため、説明は省略する。次に上記で生成した画像データについて説明する。図18aは、対物レンズ倍率40倍で撮影したZスタック画像から、対物レンズ倍率40倍、20倍、10倍、5倍に対応する被写界深度の画像を同一の画像サイズで計算した例を示している。2301、2302、2303、2304はそれぞれ対物レンズ40倍、20倍、10倍、5倍に対応する被写界深度を持つ画像を表す。40倍の画像データ2301については、被写界深度拡大処理は必要なく、合成焦点面の位置のレイヤー画像がそのまま利用される。このように画像サイズを同一にした、被写界深度を変えた複数枚の画像は、光学顕微鏡の操作において開口絞りを切替えて得られる画像に近く、ぼけの情報によってユーザに標本の奥行き情報を提示する効果がある。図18bは対物レンズ倍率40倍で撮影したZスタック画像から、対物レンズ倍率40倍、20倍、10倍、5倍に対応する被写界深度の画像を、対物レンズ倍率に応じて縮小して計算した場合の例を示している。2311、2312、2313、2314はそれぞれ対物レンズ40倍、20倍、10倍、5倍に対応する被写界深度を持つ画像を表す。図に示すように2312、2313、2314の画像の一辺はそれぞれ2311の画像の一辺の1/2、1/4、1/8となっている。図18bのように画像サイズを撮影時の対物レンズ倍率との比に応じて縮小することでディスプレイ112で画像を画面全面に表示したときに広い視野を表示できるようになる。従って、対物レンズ倍率に応じて画像サイズを縮小したデータをファイル内に階層的に格納することで、データ容量を抑えて、光学顕微鏡の低倍観察像被写界深度が深く視野が広いと高倍観察像被写界深度が浅く視野が狭いを高速に表示できるようになる。図18bで示した画像サイズの縮小は、深度合成データ生成処理S903で実現する。このとき、画像サイズの縮小処理のタイミングは複数存在する。単純な実装では、図14の画像合成処理S1602の後に対物レンズ倍率に応じた縮小処理不図示を挿入することで実現できる。しかし、画像合成処理S1602において、計算負荷の大きい処理の前に縮小することで計算速度を高めることも可能である。例えば、深度合成データ生成処理S903の内部処理として1評価値に基づく選択統合手法を用いる場合、対物レンズ倍率に応じて、レイヤー画像や2次元PSFを縮小してから処理するとコンボリューションの演算コストを低減できる。また2重み付け合成手法を用いる場合や33次元フィルタリング手法を用いる場合には、周波数領域上で原点を中心とした縮小した画像サイズに相当する低周波成分領域のみを処理対象とし、その処理対象領域のみ逆変換を行うと良い。これにより、逆フーリエ変換の計算負荷を低減して縮小した合成画像が得られる。2を例に説明すると、撮像時の対物レンズ倍率と被写界深度拡大後の対物レンズ倍率の比が1/Nとすると、合成周波数データAu,vの内、原点を中心として画像サイズが1/Nの範囲の周波数データを取得し、その成分のみを逆フーリエ変換すると良い。3においても3次元フーリエ変換像のうち、逆フーリエ変換するUV平面の範囲は原点を中心として画像サイズが1/Nの範囲の周波数データで良い。また、本実施例では、複数の対物レンズ倍率や開口絞りにおいて、共通の合成焦点面になるように設定することも可能である。合成焦点面を共通にすることにより、複数の被写界深度拡大像を画面上で切り替えて表示したときに、光学顕微鏡の対物レンズ倍率や開口絞りを調整した場合と同じような視覚効果が得られるという利点がある。簡単な方法としては、最も被写界深度が浅い場合における合成焦点面を共通の焦点面として用いるとよい。別の方法としては、対物レンズ倍率毎に求めた画像評価積分値に対し、ウェイトを掛けて合計した値を新たな画像評価積分値とし、そこから合成焦点面を定めるとよい。対物レンズ倍率が40倍、20倍、10倍、5倍における画像評価積分値をそれぞれEaz、Ebz、Ecz、Edzとする。このとき、下記の<6>式に示す新たな画像評価積分値fzを用いて、倍率にかかわらず合成焦点面を一つに定めることが出来る。Wa、Wb、Wc、Wdは重み係数<6>式では、重視する倍率の重み係数を大きくすることで、合成焦点面の位置を制御出来る。以上、本実施例の構成によれば、Zスタック画像から被写界深度を変えた複数の画像を生成し、ファイルに格納することにより、画像処理装置110上で光学顕微鏡に近い操作性を提供でき、ユーザ病理医の診断効率を高める効果が実現できる。[実施例4]実施例1~3では一つの被写界深度に対して一つの合成焦点面を設定する場合について説明したが、本実施例では、一つの被写界深度に対し、複数の合成焦点面を設定し、被写界深度を拡大したZスタック画像を生成する場合について説明する。図9aを用いて処理の流れを説明する。まず、深度拡大特性情報取得処理S901において、画像処理装置110は、深度拡大設定S401において設定した深度拡大原情報から被写界深度と結像特性情報を取得する。次に合成焦点面対象面決定処理S902では、画像処理装置110は、実施例1と同様の手法で合成焦点面と合成対象面を決定する。本実施例では、事前にユーザが指定するZスタック枚数から、複数の合成焦点面と合成対象面を決定していく。なお、Zスタック枚数は深度拡大設定S401で設定する。図9cを用いて本実施例における合成焦点面対象面決定処理S902の内部処理を説明する。画像評価値計算処理S1101、画像評価積分値計算処理S1102については実施例1と同様のため、説明を省略する。合成焦点面決定処理S1103では、画像処理装置110が、実施例1と同じ方法で合成焦点面を決定し、これを第1の合成焦点面とする。次に、画像処理装置110は、深度拡大設定S401で設定したZスタック枚数に従い、複数の合成焦点面を決定していく。例えば、ユーザからm[μm]間隔でZスタック枚数3枚と指定された場合、画像処理装置110は、第1の合成焦点面を中心として、Z方向の上下にm[μm]離れた第2、第3の合成焦点面を設定する。なお、指定された位置にレイヤー画像が存在しない場合には、指定された位置に最も近いレイヤー画像のZ位置に合成焦点面を設定する。また深度拡大設定S401での設定により、Zスタック間隔には被写界深度拡大後の被写界深度距離を設定することも可能である。このとき被写界深度の値として深度拡大特性情報取得処理S901で求めた値を用いる。次に合成対象面決定処理S1104では、画像処理装置110が、前段で定めた複数の合成焦点面に対して、それぞれ合成対象面を設定する。合成対象面はそれぞれの合成焦点面に対し、個別に求める。なお、合成対象面の求め方は実施例1と同様のため、説明は省略する。上記のように、合成焦点面対象面決定処理S902で定めた複数の合成焦点面の情報および合成対象面のデータを用いて、深度合成データ生成処理S903では、被写界深度を拡大した複数の画像データを生成し、これらをZスタック画像として出力する。ここでは各々の合成焦点面および合成対象面からの被写界深度拡大処理を複数実行する。個々の処理については実施例1と同様のため、説明は省略する。以上、一つの被写界深度に対して、被写界深度を拡大した複数の画像データZスタックを生成する方法について述べた。本実施例においても、実施例1と同様に、事前に合成焦点面、合成対象面を定めることにより、深度合成データ生成処理S903における演算コストを抑制することで撮像装置のスループットや応答速度を高めることが出来る。また、本実施例を実施例2、3と組合せることで、従来のバーチャルスライドの機能を強化できる。実施例2との組合せにより、検体の注目箇所でピンぼけが無く、さらにZスタック画像により奥行きまで観察できる。また、実施例3との組合せにより、被写界深度を変えながらの観察が可能となる。[実施例5]実施例1~4では、深度拡大処理S404を現像補正処理S405の前段で行う場合について説明したが、深度拡大処理S404は本計測処理S403の直後ではなく伝送処理S408の後で実施しても良い。例えば、本計測ユニット101で撮像し、現像補正処理や合成処理、圧縮処理を施した画像データを画像処理装置110上の画像表示アプリケーションで表示しているとする。実施例1~4で述べた被写界深度拡大処理を、前記の画像表示アプリケーション上で高速に実行することで、被写界深度拡大像の対話的な表示が可能となり、光学顕微鏡の観察に近い操作性が実現でき、ユーザ病理医の診断効率を高める効果が期待できる。"