(57)【特許請求の範囲】
【請求項１】
検体を撮像して画像データを取得する撮像システムであって、
検体の像を拡大する結像光学系、及び、前記結像光学系により拡大された検体の像を撮像する撮像手段を有する撮像ユニットと、
前記撮像ユニットから得られる検体の画像データを分析することにより、当該検体の染色方法を推定する染色方法推定手段と、
前記染色方法推定手段により推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するデータ削減処理手段と、を有する
ことを特徴とする撮像システム。
【請求項２】
前記撮像手段はイメージセンサを有し、
前記データ削減処理手段は、前記イメージセンサの複数の画素の電荷を統合し１つの画素の情報として読み出すビニング処理によって、前記撮像手段から読み出す画像データの解像度を低下させる手段を含む
ことを特徴とする請求項１に記載の撮像システム。
【請求項３】
前記データ削減処理手段は、前記撮像手段から読み出された画像データに対し解像度変換処理を適用することによって、前記画像データの解像度を低下させる手段を含む
ことを特徴とする請求項１または２に記載の撮像システム。
【請求項４】
前記データ削減処理手段は、前記撮像ユニットから得られる画像データを圧縮する圧縮手段を含む
ことを特徴とする請求項１〜３のうちいずれか１項に記載の撮像システム。
【請求項５】
前記圧縮手段は、前記第１の染色方法と推定された検体の画像データの圧縮率を、前記第２の染色方法と推定された検体の画像データの圧縮率よりも高める
ことを特徴とする請求項４に記載の撮像システム。
【請求項６】
前記圧縮手段は、圧縮率を高めるために画像データの色空間を変換する色変換手段を含み、
前記色変換手段は、前記染色方法推定手段により推定された染色方法に応じて色変換式を変更する
ことを特徴とする請求項４または５に記載の撮像システム。
【請求項７】
前記撮像ユニットは、第１の撮像ユニットと前記第１の撮像ユニットよりも高精細の撮像を行う第２の撮像ユニットとを有し、
前記染色方法推定手段は、前記第１の撮像ユニットから得られる画像データを用いて染色方法を推定し、
前記データ削減処理手段は、前記第２の撮像ユニットから得られる画像データのデータ量を削減する
ことを特徴とする請求項１〜６のうちいずれか１項に記載の撮像システム。
【請求項８】
前記第２の撮像ユニットは、倍率を変更可能な結像光学系を有しており、前記第１の染色方法と推定された検体を撮像する場合は、前記第２の染色方法と推定された検体を撮像する場合よりも、前記結像光学系の倍率を低くする
ことを特徴とする請求項７に記載の撮像システム。
【請求項９】
検体の画像データを処理する画像処理装置であって、
検体の画像データを取得する画像取得手段と、
前記画像取得手段で取得した検体の画像データを分析することにより、当該検体の染色方法を推定する染色方法推定手段と、
前記染色方法推定手段により推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するデータ削減処理手段と、を有し、
前記データ削減処理手段は、前記画像取得手段で取得した画像データに対し解像度変換処理を適用することによって、前記画像データの解像度を低下させる手段を含む
ことを特徴とする画像処理装置。
【請求項１０】
前記染色方法推定手段は、前記画像取得手段で取得した画像データの縮小画像のデータを用いて染色方法を推定する
ことを特徴とする請求項９に記載の画像処理装置。
【請求項１１】
前記データ削減処理手段は、前記画像取得手段で取得した画像データのファイルフォーマットに応じてデータ量の削減方法を決定する
ことを特徴とする請求項９または１０に記載の画像処理装置。
【請求項１２】
検体の画像データを処理するプログラムであって、
検体の画像データを取得するステップと、
取得した検体の画像データを分析することにより、当該検体の染色方法を推定するステップと、
推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するステップと、
をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、多数の検体を撮像して画像データを取得する撮像システム、および、多数の検体の画像データを処理する画像処理装置に関する。
【背景技術】
【０００２】
近年、病理分野において、病理診断のツールである光学顕微鏡の代替として、プレパラートに載置された被検試料の撮像と画像のデジタル化によってディスプレイ上での病理診断を可能とするバーチャル・スライド・システムが注目を集めている。バーチャル・スライド・システムを用いた病理診断画像のデジタル化により、従来の被検試料の光学顕微鏡像をデジタルデータとして取り扱うことが可能となる。その結果、遠隔診断の迅速化、デジタル画像を用いた患者への説明、希少症例の共有化、教育・実習の効率化、などのメリットが得られると期待されている。
【０００３】
光学顕微鏡と同等程度の操作をバーチャル・スライド・システムで実現するためには、プレパラート上の被検試料全体をデジタル化する必要がある。被検試料全体のデジタル化により、バーチャル・スライド・システムで作成したデジタルデータをＰＣやワークステーション上で動作するビューワソフトで観察することができる。被検試料全体をデジタル化した場合の画素数は、通常、数億画素から数十億画素と非常に大きなデータ量となる。
【０００４】
バーチャル・スライド・システムで作成したデータ量は膨大であるが、それゆえ、ビューワで拡大・縮小処理を行うことでミクロ（細部拡大像）からマクロ（全体俯瞰像）まで観察することが可能となり、種々の利便性を提供する。必要な情報を予めすべて取得しておくことで、低倍画像から高倍画像までユーザが求める解像度・倍率による即時の表示が可能となる。
【０００５】
一般に病理診断では医師は多数のスライドを観察するため、医療の現場では、大量のプレパラートを短時間（例えば一晩）でデジタル化できるバッチ処理機能を備えたバーチャル・スライド・システムが必要とされている。すなわち、バーチャル・スライド・システムでは、スループット（単位時間あたりの処理枚数）の向上が望まれている。また、バッチ処理の実現のためには、人の判断を必要としない、撮像条件の自動判断の仕組みも必要となる。
【０００６】
特許文献１では、バッチ処理の効率化のために、まず検体を低解像度のハイスループットモードで撮像し、撮像した画像内に特定の信号を検出したときのみ高精細モードに移行する、という構成が開示されている。また撮像条件の自動判断に関して、特許文献２には、撮影した画像から標本の染色方法を推定し、それぞれの染色方法に適した色変換処理を施すことにより、標本毎の高画質化を図る方法が開示されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１０−２１６９１９号公報
【特許文献２】特開２００９−２９４４６３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１の方法は、高解像度で撮像すべき検体の割合が少ない場合には、処理時間の短縮を期待できる。しかしながら、前述のように病理診断用の検体をバッチ処理する場合
には、原則として全ての検体（プレパラート）を高解像度で撮像することとなるため、特許文献１の方法を適用したとしてもスループットの向上を図ることはできない。
【０００９】
また、病理診断には大別して組織診と細胞診とがあり、それぞれ標本の作り方や観察（診断）の目的が相違する。しかしながら従来のシステムでは、組織診用と細胞診用の検体が混在していたとしても、検体の種類によらず同じ条件で処理を行うことしかできなかった。そのため、検体によっては、必要以上のスペック（例えば解像度、画質、圧縮率など）の画像データが生成、記録されることとなり、データ量の増大による処理効率の低下を招いている。特に画像の高解像度化や大サイズ化が進むにつれ、この問題は深刻となる。とはいえ、データ量を削減するために、全ての検体について、一律に解像度や画質を落としたり、圧縮率を高めたりしてしまうと、観察や診断に必要な情報が欠落してしまうおそれもある。
【００１０】
本発明は、このような問題点に鑑みてなされたものであり、多数の検体の画像データをバッチ処理する際のスループットの向上を図るための技術を提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明の第１態様は、検体を撮像して画像データを取得する撮像システムであって、検体の像を拡大する結像光学系、及び、前記結像光学系により拡大された検体の像を撮像する撮像手段を有する撮像ユニットと、前記撮像ユニットから得られる検体の画像データを分析することにより、当該検体の染色方法を推定する染色方法推定手段と、前記染色方法推定手段により推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するデータ削減処理手段と、を有する撮像システムを提供する。
【００１２】
本発明の第２態様は、検体の画像データを処理する画像処理装置であって、検体の画像データを取得する画像取得手段と、前記画像取得手段で取得した検体の画像データを分析することにより、当該検体の染色方法を推定する染色方法推定手段と、前記染色方法推定手段により推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するデータ削減処理手段と、を有し、前記データ削減処理手段は、前記画像取得手段で取得した画像データに対し解像度変換処理を適用することによって、前記画像データの解像度を低下させる手段を含む画像処理装置を提供する。
【００１３】
本発明の第３態様は、検体の画像データを処理するプログラムであって、検体の画像データを取得するステップと、取得した検体の画像データを分析することにより、当該検体の染色方法を推定するステップと、推定された染色方法が組織診用の検体の染色に用いられる第１の染色方法であった場合に、当該検体の画像データのデータ量を、前記画像データの解像度を当該推定された染色方法が細胞診用の検体の染色に用いられる第２の染色方法であった場合よりも低下させることで、削減するステップと、をコンピュータに実行させるプログラムを提供する。
【発明の効果】
【００１４】
本発明によれば、多数の検体の画像データをバッチ処理する際のスループットの向上を図ることができる。
【図面の簡単な説明】
【００１５】
【図１】バーチャル・スライド・システムの構成図。
【図２】本計測ユニットとプレ計測ユニットの構成図。
【図３】画像処理装置の内部構成図。
【図４】実施例１の本計測処理とプレ計測推定制御処理のフローチャート。
【図５】本計測とプレ計測の撮像領域を説明する図。
【図６】ＸＹステージおよびＺステージの移動方向を説明する図。
【図７】実施例１のプレ計測データ取得処理と染色方法推定処理のフローチャート。
【図８】実施例１の染色方法推定に用いる色ヒストグラムを説明する図。
【図９】実施例１の染色方法推定方法を説明する図。
【図１０】実施例１の撮像条件設定処理のフローチャート。
【図１１】実施例１のＣＣＤイメージセンサにおけるビニング処理を説明する図。
【図１２】実施例２のＣＭＯＳイメージセンサにおけるビニング処理を説明する図。
【図１３】実施例３のＪＰＥＧ圧縮符号化方式の符号化・復号化装置のブロック図。
【図１４】実施例４の撮像条件設定処理のフローチャート。
【図１５】実施例４の本計測ユニットの結像光学系の構成を示す図。
【図１６】実施例５における染色推定制御処理の処理手順を示すフローチャート。
【発明を実施するための形態】
【００１６】
［実施例１］
（システムの全体構成）
図１は、本発明の撮像システムの一実施例である、バーチャル・スライド・システムの構成を示している。
【００１７】
バーチャル・スライド・システムは、検体の撮像データを取得する撮像装置（バーチャル・スライド・スキャナとも呼ぶ）１２０とそのデータ処理・制御を行う画像処理装置（ホストコンピュータとも呼ぶ）１１０およびその周辺機器から構成されている。
画像処理装置１１０にはユーザからの入力をキーボードやマウスなどの操作デバイスを通じて受け付ける操作部１１１、処理画像を表示する表示部１１２が接続される。また画像処理装置１１０には、記憶装置１１３、他のコンピュータシステム１１４が接続されている。
【００１８】
多数の検体（プレパラート）の撮像をバッチ処理で行う場合、画像処理装置１１０の制御の下で撮像装置１２０が各検体を順に撮像し、画像処理装置１１０が各検体の画像データに対し必要な処理を施す。そして得られた各検体の画像データは、大容量のデータストレージである記憶装置１１３又は他のコンピュータシステム１１４へ伝送され、蓄積される。
【００１９】
撮像装置１２０での撮像（プレ計測および本計測）は、ユーザの入力を受けて画像処理装置１１０がコントローラ１０８へ指示を送り、次にコントローラ１０８が本計測ユニット１０１とプレ計測ユニット１０２を制御することで実現される。
本計測ユニット１０１はプレパラート内の検体診断のための高精細な画像を取得する撮像ユニットである。プレ計測ユニット１０２は本計測に先立って撮像を行う撮像ユニットで、本計測で精度の良い画像取得をするための撮像制御情報取得を目的とした画像取得を行う。詳しくは後述するが、本実施例において特徴的な処理は、第１の撮像ユニットであるプレ計測ユニット１０２で撮像したデータを用いて、第２の撮像ユニットである本計測ユニット１０１を制御することで、検体に応じたデータ量の削減を行う点にある。
【００２０】
コントローラ１０８には変位計１０３が接続され、本計測ユニット１０１またはプレ計測ユニット１０２内のステージに設置されるプレパラートの位置や距離が測定できる構成となっている。変位計１０３は、本計測およびプレ計測を行うにあたり、プレパラート内の検体の厚みを計測するために用いられる。
【００２１】
またコントローラ１０８には本計測ユニット１０１およびプレ計測ユニット１０２の撮像条件を制御するための、開口絞り制御１０４、ステージ制御１０５、照明制御１０６、センサ制御１０７が接続されている。そして、それぞれはコントローラ１０８からの制御信号に従って、開口絞り、ステージ、照明、イメージセンサの動作を制御する構成となっている。
【００２２】
ステージにはプレパラートを光軸に垂直な方向に移動するＸＹステージと光軸に沿った方向に移動するＺステージがある。それぞれ、ＸＹステージは光軸に垂直な方向に分布する検体画像を、Ｚステージは奥行き方向に焦点位置を変えた画像を撮像するために用いる。図示しないが、撮像装置１２０は、複数のプレパラートをセット可能なラックと、ラックからステージ上の撮像位置へとプレパラートを送り出す搬送機構とが設けられている。バッチ処理の場合は、コントローラ１０８が搬送機構を制御することで、ラックからプレパラートを１枚ずつプレ計測ユニット１０２のステージ、本計測ユニット１０１のステージの順に送り出す。
【００２３】
本計測ユニット１０１、プレ計測ユニット１０２には撮像した画像を用いてオートフォーカスを実現するＡＦユニット１０９が接続されている。ＡＦユニット１０９は、コントローラ１０８を介して、本計測ユニット１０１、プレ計測ユニット１０２のステージの位置を制御することで合焦位置を探し出すことが出来る。オートフォーカスの方式は画像を用いるパッシブ型であり、公知の位相差検出方式やコントラスト検出方式が用いられる。
【００２４】
（本計測ユニット）
図２（ａ）は実施例１における本計測ユニット１０１の内部構成を示す図である。
光源２０１の光は照明光学系２０２を通じ、光量ムラが無いように均一化され、ステージ２０３の上に設置されたプレパラート２０４を照射する。プレパラート２０４は観察対象となる組織の切片や塗抹した細胞をスライドグラス上に貼り付け、封入剤とともにカバーグラスの下に固定したものであり、検体（被写体）を観察可能な状態に準備したものである。
【００２５】
結像光学系２０５は検体の像を拡大して撮像手段である撮像部２０７に導くものである。プレパラート２０４を通った光は結像光学系２０５を介して、撮像部２０７上の撮像面で結像する。結像光学系２０５の中には開口絞り２０６が存在し、開口絞り２０６を調整することで被写界深度が制御できる。
【００２６】
撮像面で結像した検体の像は、複数のイメージセンサから構成される撮像部２０７で光電変換され、Ａ／Ｄ変換を実行した後、像は電気信号として画像処理装置１１０に送られる。本実施例ではＡ／Ｄ変換を実行した後のノイズ除去や色変換処理、鮮鋭化処理に代表される現像処理は画像処理装置１１０内部で行うとして説明する。しかし、現像処理は撮像部２０７に接続された専用の画像処理ユニット（不図示）で行い、その後画像処理装置１１０にデータを送信することも可能であり、そのような形態での実施も本発明の範疇とする。
【００２７】
（プレ計測ユニット）
図２（ｂ）は実施例１におけるプレ計測ユニット１０２の内部構成を示す図である。
光源３０１の光は照明光学系３０２を通じ、光量ムラが無いように均一化され、ステージ３０３の上に設置されたプレパラート２０４を照射する。プレパラート２０４を通った光は結像光学系３０５により、撮像部３０７上の撮像面で結像する。結像光学系３０５の中には開口絞り３０６が存在し、開口絞り３０６を調整することで被写界深度が制御できる。
【００２８】
撮像面で結像した検体の像は、イメージセンサを備えた撮像部３０７で光電変換され、Ａ／Ｄ変換を実行した後、像は電気信号として画像処理装置１１０に送られる。本実施例ではＡ／Ｄ変換を実行した後のノイズ除去や色変換処理、鮮鋭化処理に代表される現像処理は画像処理装置１１０内部で行うとして説明する。しかし、現像処理は撮像部３０７に接続された専用の画像処理ユニット（不図示）で行い、その後画像処理装置１１０にデータを送信することも可能であり、そのような形態での実施も本発明の範疇とする。
【００２９】
（画像処理装置）
図３は本発明における画像処理装置（ホストコンピュータ）１１０の内部構成を示す図である。
ＣＰＵ４０１はＲＡＭ４０２やＲＯＭ４０３に格納されているプログラムやデータを用いて画像処理装置全体の制御を行う。またＣＰＵ４０１は、以降の実施例で説明する各種演算処理、データ処理、例えば、後述する染色方法推定処理や撮像条件設定処理を行う。
【００３０】
ＲＡＭ４０２は外部記憶装置４１１からロードされたプログラムやデータ、他のコンピュータシステム４０５からＩ／Ｆ（インターフェース）４０４を介してダウンロードしたプログラムやデータを一時的に記憶するエリアを備える。またＲＡＭ４０２は、ＣＰＵ４０１が各種の処理を行うために必要とするワークエリアを備える。ＲＯＭ４０３はコンピュータの機能プログラムや設定データなどを記憶する。ディスプレイ制御装置４０６は画像や文字等をディスプレイ４０７に表示させるための制御処理を行う。ディスプレイ４０７はユーザに入力を求めるための画面表示を行うとともに、撮像装置１２０から取得しＣＰＵ４０１で処理した画像データを画像表示する。
【００３１】
操作入力デバイス４０９はキーボードやマウスなどＣＰＵ４０１に各種の指示を入力することのできるデバイスにより構成される。ユーザは撮像装置１２０の動作を制御する情報を操作入力デバイス４０９により入力する。４０８は操作入力デバイス４０９を介して入力された各種の指示等をＣＰＵ４０１に通知するためのＩ／Ｏである。
【００３２】
外部記憶装置４１１はハードディスクなどの大容量情報記憶装置であり、ＯＳ（オペレーティングシステム）や以降の実施例で説明する処理をＣＰＵ４０１に実行させるためのプログラムやバッチ処理によりスキャンした画像データなどを記憶する。
外部記憶装置４１１への情報の書き込みや外部記憶装置４１１からの情報の読み出しはＩ／Ｏ４１０を介して行われる。コントローラ４１３は撮像装置１２０を制御するためのユニットであり、Ｉ／Ｆ（インターフェース）４１２を介してＣＰＵ４０１と制御信号および応答信号のやり取りをする。
【００３３】
コントローラ４１３は本計測ユニット１０１およびプレ計測ユニット１０２を制御する機能を持つ。Ｉ／Ｆ（インターフェース）４１４には上述以外のインターフェース、例えば、ＣＭＯＳイメージセンサやＣＣＤイメージセンサの出力データを取り込むための外部インターフェースが接続されている。なおインターフェースとしてはＵＳＢ、ＩＥＥＥ１３９４などのシリアルインターフェースやカメラリンクなどのインターフェースを使うことが出来る。このＩ／Ｆ４１４を通じて様々な周辺機器が接続できる。
【００３４】
（本計測処理）
本実施例のバーチャル・スライド・システムは、検体（標本）の染色方法を推定するための「プレ計測」と、推定した染色方法に応じた解像度で検体を撮像する「本計測」とを実行する。ここでは、プレ計測で得られた低解像度の画像データを分析して染色方法を推定し、その結果に従って本計測ユニット１０１の動作を制御する処理を「プレ計測推定制御処理」とよぶ。以下、実際の処理の順番とは逆になるが、最初に本計測処理の概要について説明し、その後で本実施例の特徴であるプレ計測推定制御処理について詳しく説明する。
【００３５】
図４（ａ）は本計測処理の処理フローを示す図である。
本計測データ取得処理Ｓ５０１では、コントローラ１０８の制御により、本計測ユニット１０１がプレパラートを撮像し、画像処理装置１１０に画像データを送信する。
次に現像・補正処理Ｓ５０２において、画像処理装置１１０が画像データに対し色変換
処理、鮮鋭化処理、ノイズ低減処理などを実施する。これらの処理により、モニタ表示される画像の色を標本の実際の色に近づけることができるとともに、画像のノイズが抑えられる。
次に合成処理Ｓ５０３では、画像処理装置１１０が、分割して撮像した物体面の画像をつなぎ合わせ、プレパラート上の対象領域（例えば２０ｍｍ角の領域）を映した像を作る。
次に圧縮処理Ｓ５０４では、画像処理装置１１０が、合成したデータを圧縮し、データ容量を削減する。圧縮方法としては静止画圧縮符号化方式のＪＰＥＧやＪＰＥＧ２０００などを用いることが可能である。続いて、伝送処理Ｓ５０５では、画像処理装置１１０が記憶装置１１３に画像データを送り、保存する。或いは画像処理装置１１０は、ネットワークＩ／Ｆを介してネットワーク上のコンピュータシステム１１４や画像サーバに画像データを送信してもよい。
【００３６】
（本計測処理：本計測データ取得処理Ｓ５０１）
図５、図６を用いて本計測データ取得処理Ｓ５０１について説明する。
図５（ａ）はプレパラートの模式図である。スライドグラス６１０上にはカバーグラス６１１によって検体６００を封入した領域とラベルエリア６１２がある。本実施例の本計測データ取得処理Ｓ５０１では、カバーグラス６１１があると想定される領域を撮像対象とする。なお、プレ計測で検体６００が存在する外接矩形領域を求め、本計測ではその領域のみを撮像することでデータ量を削減することも好ましい。
【００３７】
図５（ｂ）は本計測データ取得処理Ｓ５０１でカバーグラス６１１が存在する領域を複数領域に分割して撮像する様子を表している。図５（ｃ）は撮像面を表しており、有効視野６０２は本計測ユニット１０１の結像光学系２０５を通して像が映る範囲を、センサ有効領域６０３は撮像部２０７のイメージセンサで撮像可能な領域を示している。
物体面における撮像領域６０１（斜線領域）は本計測ユニット１０１の結像光学系２０５を通じて結像し、撮像面における撮像領域６０４に対応する。
【００３８】
図５（ｃ）に示すように、センサ有効領域６０３には、撮像領域６０４に対して若干広い領域が割り当てられている。これは、結像光学系２０５が有する光学的な収差やイメージセンサの取付位置のズレを許容するためのマージンである。すなわち、光学的な収差やセンサの取付位置のズレがあったとしても、物体面の撮像領域６０１がセンサ有効領域６０３内に収まるようにしてある。合成処理Ｓ６０３では、センサ有効領域６０３の画像に対して収差補正や位置ズレ補正を行い、補正した画像から撮像領域６０４に対応する部分を切り出し、その画像を合成に用いる。
【００３９】
図６（ａ）は図５（ｂ）で示す分割領域をラスタースキャン順に撮像する場合のステージ２０３のＸＹ方向の移動方向と順序を表している。プレパラートの左上から右下の領域を撮像するには、プレパラートを搭載するステージ２０３は逆方向の右下から左上に移動する。
このように撮像領域を複数領域に分割し、ステージ２０３を移動しながら繰り返し撮像することで比較的小さなイメージセンサでも広い領域を撮像することができる。
【００４０】
図６（ｂ）は本計測データ取得処理Ｓ５０１で、焦点位置（観察深さ、フォーカス位置、或いはピント位置ともいう）の異なる複数の画像を取得する場合のステージ２０３のＺ方向（奥行き方向）の移動方向を示している。図に示すように、プレパラート２０４内の検体の上側（カバーグラス裏面側）に焦点位置を移動するには、結像光学系２０５の光軸方向に沿って、ステージ２０３をＺ方向下向きに移動させる。逆に、検体の下側（スライドグラス表面側）に焦点位置を移動するには、ステージ２０３をＺ方向上向きに移動させる。焦点位置を変更しながら検体を複数回撮像することで、焦点位置の異なる複数枚の画
像データを取得する処理は一般に「Ｚスタック」とよばれている。
撮像部２０７、或いは、撮像部２０７とステージ２０３の両方を結像光学系２０５の光軸方向に沿って移動することでも焦点位置を変更することができる。また、結像光学系２０５のレンズを制御することにより光学的に焦点距離を変更することでも、焦点位置を変更可能である。なお、プレ計測ユニット１０２のステージ機構は本計測ユニット１０１とほぼ同じため、説明は省略する。
【００４１】
（プレ計測推定制御処理）
図４（ｂ）はプレ計測推定制御処理の処理フローを示す図である。
プレ計測データ取得処理Ｓ９０１では、コントローラ１０８の制御により、プレ計測ユニット１０２がプレパラートを撮像し、画像処理装置１１０に画像データを送信する。
次に、染色方法推定処理Ｓ９０２では、画像処理装置１１０が、撮像部３０７で撮像した画像を解析し、検体（標本）に施された染色方法を推定する。
次に、撮像条件設定処理Ｓ９０３では、画像処理装置１１０が、染色方法推定処理Ｓ９０２で推定した情報を元に、本計測処理Ｓ９０４での撮像条件を本計測ユニット１０１に設定する。
最後に本計測処理Ｓ９０４では、撮像条件設定処理Ｓ９０３で設定された撮像条件を使って、スライドグラスとカバーグラスで挟まれた検体内の所望の位置にフォーカスを合わせ、本計測ユニット１０１で画像取得を行う。そして図４（ａ）で説明した処理により、高解像度の合成画像を生成する。
【００４２】
図５（ｄ）はプレ計測におけるプレパラート２０４の撮像領域６０５を示している。プレ計測は本計測を精度よく撮像するための撮像制御情報取得を目的としており、画像の特徴が大まかに把握できれば良く、本計測ほどの倍率は必要ない。また、プレ計測では被写界深度を深くするとよい。これにより、検体にピントを合わせやすくできる。
【００４３】
本実施例のプレ計測では低倍率でプレパラート２０４の全体を撮像する。また、本計測と異なり複数の領域に分割することなく１つのイメージセンサでプレパラート２０４の全体を一括撮像する。これによりプレ計測ユニット１０２の構成を簡易化することができるとともに、プレ計測に要する時間、ひいてはプレ計測と本計測を含む撮像処理全体の時間を短縮することができる。ただし、プレ計測でも本計測相当の解像度が必要な場合、本計測同様に倍率を上げて、物体面での撮像対象領域を複数に分割して撮像しても良い。
【００４４】
（プレ計測推定制御処理：プレ計測データ取得処理Ｓ９０１）
図７（ａ）は、本実施例のプレ計測データ取得処理Ｓ９０１の詳細を示している。
ステージ設定処理Ｓ１００１では、コントローラ１０８が搬送機構を制御して、プレ計測ユニット１０２のステージ３０３にプレパラート２０４をセットする。
【００４５】
照明照射処理Ｓ１００２では、光源３０１が点灯し、プレパラート２０４に光を照射する。撮像処理Ｓ１００３では、照明光学系３０２、プレパラート２０４、結像光学系３０５を通って撮像面に結像した像を撮像部３０７のイメージセンサで画像化する。なお、本実施例では、ＲＧＢの３つの光源３０１で順番に露光し、３回撮像を行うことで、カラー画像を取得する。すなわち、Ｓ１００２とＳ１００３が３回繰り返される。
【００４６】
現像・合成処理Ｓ１００４では、撮像処理Ｓ１００３で取得した生データ（ＲＡＷデータ）に対し、画像処理装置１１０が現像・合成処理を行う。現像・合成処理Ｓ１００４では、色変換、ノイズ除去処理などを行う。ｓＲＧＢやＡｄｏｂｅＲＧＢなどの様々な色空間規格が存在し、そのいずれを用いても良いが、本実施例では代表的なｓＲＧＢ色空間に色変換する。
【００４７】
（プレ計測推定制御処理：染色方法推定処理Ｓ９０２）
図７（ｂ）は、本実施例の特徴である染色方法推定処理Ｓ９０２の処理内容を示している。本実施例の染色方法推定処理は、プレ計測で得られた画像データの色に基づいて検体の染色方法を推定する処理である。
【００４８】
まず色空間変換処理Ｓ１１０１では、プレ計測で得られた画像データの色空間変換が行われる。色空間としては、ｘｙＹ表色系（ｘｙ色度図）、輝度色差信号ＹＵＶ、均等色空間ＣＩＥ  Ｌ＊ａ＊ｂ＊、ＨＳＶ色空間、ＨＬＳ色空間などがある。本実施例ではＣＩＥ
Ｌ＊ａ＊ｂ＊色空間に変換される。なお、ｓＲＧＢのまま以降の処理を行う場合には、Ｓ１１０１の処理を省略することもできる。
【００４９】
ヒストグラム生成処理Ｓ１１０２では、画像処理装置１１０が、色空間変換した画像データから色ヒストグラム（色出現分布情報）を作成する。
図８（ａ）、図８（ｂ）はヒストグラム生成処理Ｓ１１０２の一例を示す。図８（ａ）に示すように、例えば、Ｌ＊ａ＊ｂ＊色空間をＬ＊軸を中心軸として３０度ずつ１２等分し、それぞれの領域Ａ１〜Ａ１２に出現する画素の数を数える。そして図８（ｂ）に示すように、プレ計測で得られた画像データに対する一次元ヒストグラムを描く。図８（ｂ）の横軸は領域Ａ１〜Ａ１２であり、縦軸は画素の出現頻度（画素の数）である。
【００５０】
図５（ａ）と図５（ｄ）から分かるように、プレ計測撮像領域６０５の中には検体６００が存在しない部分もある。検体６００が存在しない部分は、照明の色、即ち、無彩色の色となる。検体以外の画素を取り除いた方が精度良く検体の染色方法を推定できるため、Ｌ＊軸から所定の距離にある画素（ほぼ無彩色の画素）をヒストグラムから取り除くと良い。
【００５１】
マッチング度算出処理Ｓ１１０３では、画像処理装置１１０が、データベース１３０４から染色方法毎の色ヒストグラム（色出現分布情報）を取得する。なおデータベース１３０４には、各染色方法の検体サンプルを用いて作成された、各染色方法の典型的な色出現分布を示す色ヒストグラムがあらかじめ格納されている。図９の１３０２はＡ染色方法の色ヒストグラム、１３０３はＢ染色方法の色ヒストグラムを示している。そして、画像処理装置１１０は、ヒストグラム生成処理Ｓ１１０２で計算したヒストグラム１３０１と、各染色方法の典型的なヒストグラム１３０２、１３０３とを比較し、染色方法毎のマッチング度を計算する。
【００５２】
ヒストグラム同士のマッチング度（類似度）は、例えば、ヒストグラムの内積や、ヒストグラム・インターセクションなどで評価できる。正規化相互相関（内積を計算する前にそれぞれの１次元ヒストグラムの総和が０になるように正規化してから内積を計算）を用いれば、内積の最大値は１に抑えられるため、次のステップで判定する際の閾値が導入しやすく、確度の低い推定が排除できる。
【００５３】
次に染色方法推定処理Ｓ１１０４では、マッチング度の最も大きい染色方法を、プレ計測した検体の染色方法として選定する。例えば、図９の例では、プレ計測した検体の色ヒストグラム１３０１は、Ａ染色方法の色ヒストグラム１３０２との内積（相関）が最も大きくなる。
【００５４】
染色方法推定処理Ｓ１１０４において、最大のマッチング度が所定の閾値より小さい場合には、「染色方法不明」と判定すると良い。誤判定により検体の染色方法と異なる染色方法を選定してしまうと、適切な撮像条件が設定できない可能性があるため、閾値の導入により誤判定の確率を下げることが望ましい。
【００５５】
染色方法推定処理Ｓ１１０４で決定された染色方法の情報は、ＲＡＭ４０２や外部記憶装置４１１など画像処理装置１１０がアクセスできる適切な場所に記憶される。なお、染色方法毎の色ヒストグラムを格納するデータベース１３０４は、外部記憶装置４１１に存在しても良いし、他のコンピュータシステム４０５に存在しても良い。
【００５６】
（プレ計測推定制御処理：撮像条件設定処理Ｓ９０３）
図１０は、撮像条件設定処理Ｓ９０３の詳細を示している。
まずＳ１４０１で、画像処理装置１１０は、有意な染色推定が出来たか否かを判定する。もし確度の高い推定ができていればＳ１４０２に進み、画像処理装置１１０は、データベース１４００にアクセスし、格納された染色方法毎の撮像制御情報を取得する。そして、Ｓ１４０３にて、画像処理装置１１０は、Ｓ１４０２で取得した撮像制御情報に基づき、本計測において撮像データ量が削減できる撮像条件を設定する。本実施例においては、撮像データ削減設定処理Ｓ１４０３は、撮像部２０７における電荷読み出し方法を変更する処理である。後で詳しく説明する。
【００５７】
Ｓ１４０１で有意な染色推定が出来なかった（すなわち染色方法不明）と判定された場合、Ｓ１４０４において、画像処理装置１１０は、予め決められている事前条件を設定する処理を行う。本実施例では、重要情報の見落としを防止するため、事前条件は「データ削減を行わない」という設定としている。ただし、図４（ａ）のステップＳ５０４の圧縮処理においてデフォルトの圧縮パラメータが設定されている場合には、その圧縮パラメータを用いて画像データは圧縮される。
【００５８】
（撮像データ量の削減）
図１１はＣＣＤイメージセンサの受光素子から画素信号を読み出す様子を模式的に示した図である。はじめに、図１１を用いてＣＣＤイメージセンサの画素信号の読み出し方法について説明する。
【００５９】
１５０１はＣＣＤイメージセンサにおける１画素の受光素子を表している。光電変換によって各画素１５０１に蓄えられた電荷を読み出すには、まず１行分の電荷を垂直方向にバケツリレー方式で伝送し、水平レジスタに移動する。次に、水平レジスタ１５０２に保持された電荷を１画素ずつ水平方向に移動させる。この操作により、各画素に蓄えられた電荷を出力ノード１５０３から順に読みだすことができる。
水平レジスタ１５０２の各画素および出力ノード１５０３は、スーパーピクセルと呼ばれる画素であり、複数の受光素子が蓄えた電荷を統合して蓄積できる大きな蓄積電子数（ウェル深度）を備えている。
【００６０】
続いてＣＣＤイメージセンサのビニング処理について説明する。ビニング処理とはイメージセンサの複数画素の受光素子の出力（電荷）を統合して、１つの画素の信号として出力する処理である。ビニング処理を行うと、解像度は低下するものの、ノイズが少なく階調性の豊かな画像を得ることができる。本明細書では、最終的に１画素に統合される受光素子の縦横の画素数を「ビニング設定」と表現することとする。
【００６１】
図１１（ａ）〜（ｄ）は２×２のビニング処理でＣＣＤイメージセンサの受光素子から信号値を取り出すまでの流れを示している。
図１１（ａ）に示すように、イメージセンサの１５１０ａで示す部分に光が照射され、４つの画素それぞれにａ個の電子が蓄積されたとする。まず図１１（ｂ）に示すように、２行分の電荷を水平レジスタ１５０２に移動することで、水平レジスタ１５０２の左２画素１５１０ｂにそれぞれ２×ａ個の電子が蓄積される。図１１（ｃ）の１５１０ｃは、水平レジスタ１５０２内の電荷を２画素分右に移動させた状態を示している。さらに水平レジスタ１５０２内の電荷を２画素分右に移動させると、図１１（ｄ）の１５１０ｄに示す
ように、出力ノード１５０３に４×ａ個の電子が集められる。これにより、４画素に蓄積された電荷の合計値を１画素の情報として読み出すことが可能となる。ビニング処理によって複数画素分の電荷を加算することにより、ランダムノイズが低減されるためＳＮ比が向上する。また、階調性に優れた画像信号を得ることができる。
図１１では、ビニング設定が２×２の例を示したが、３×３、４×４、あるいは２×４など、統合する縦横の画素数は任意に設定することが可能である。
【００６２】
（染色方法毎の撮像制御情報取得処理Ｓ１４０２）
次に、図１０の染色方法毎の撮像制御情報取得処理Ｓ１４０２について詳細に説明する。
本実施例における染色方法毎の撮像制御情報とは、ＣＣＤイメージセンサの受光素子から画素信号を読み出す際のビニング設定であり、この点が本実施例の特徴である。
【００６３】
病理診断においては、組織診か細胞診かで標本の作り方が異なる。また染色方法毎に染色される部位、色も異なる。そのため、染色方法毎に見たい個所や観察するポイントが違っている。例えば、組織診においては、立体構造を持つ検体をパラフィンで固定し、４μｍ程度の薄い切片を作って染色する。染色方法としては、ＨＥ染色（ヘマトキシリン・エオシン染色）が一般に用いられる。組織診では、細胞一つ一つの特性よりも複数の細胞が集まった組織の状態が重視されるために比較的解像度が必要とされない傾向にあり、組織間の差を見分ける為に微妙な色や階調を識別できることが重視される。そのため、データ量削減の観点でも、画質の観点でもビニング処理が適している。
【００６４】
一方、細胞診においては、剥離または穿刺吸引により取得した細胞をスライドに塗布し、固定、染色する。染色方法としては、パパニコロウ染色やギムザ染色が一般に用いられる。細胞診では、組織診に比べて、細胞一つ一つの情報が重視される傾向にあり、画像の解像度は高い方が望ましい。
【００６５】
まとめると、診断に十分な解像力を持ったバーチャル・スライド・システムにおいては、組織診で用いる染色方法（例えば、ＨＥ染色）ではＣＣＤイメージセンサのビニング処理の適用が適している。一方、細胞診で用いる染色方法（例えば、パパニコロウ染色、ギムザ染色）ではＣＣＤイメージセンサのビニング処理の必要性はさほど高くないと言える。（ただし、ノイズを抑える用途ではビニング処理は有効である）
【００６６】
次に、例を挙げて、染色方法毎の撮像制御情報取得処理Ｓ１４０２の処理を説明する。
データベース１４００に格納されている染色方法毎の撮像制御情報の一例を以下に示す。このように本実施例では、染色方法毎の撮像制御情報として、染色方法とビニング設定が関連付けられたデータを用いる。
【表１】
0005372068.tif 000002
【００６７】
例えば、図７（ｂ）のＳ１１０４において染色方法が「パパニコロウ染色」であると推定された場合、図１０のＳ１４０２では、上記テーブルを用いてビニング設定「１×１」が得られる。続いて、Ｓ１４０３において、画像処理装置１１０はコントローラ１０８を通じて、センサ制御１０７にビニング設定「１×１」の設定を行う指示を出す。以上で撮像条件設定処理Ｓ９０３は終了し、本計測処理Ｓ９０４が実行される。このとき、ビニング設定「１×１」に従って、ＣＣＤイメージセンサから画像信号が出力される。
【００６８】
一方、Ｓ１１０４において染色方法が「ＨＥ染色」であると推定された場合、Ｓ１４０２では、上記テーブルを用いてビニング設定「２×２」が得られる。この場合は、本計測処理Ｓ９０４において、ビニング設定「２×２」に従って、ＣＣＤイメージセンサから画像信号が出力される。したがって、ビニング設定「１×１」（つまりビニング処理無し）の場合に比べて、取得される画像の縦・横の画素数は半分となり、全体のサイズは４分の１となる。これにより、検体に応じて効率的なデータ量削減が実現できることが分かる。
【００６９】
以上述べた実施例１においては、複数の染色方法を施したプレパラートを混在させてバッチ処理を行った場合においても、プレ計測において染色方法を推定し、それに応じてビニング処理の設定を選ぶことにより、本計測で取得するデータ量を削減できる。そのため、過剰な解像度でのデータ取得を回避し、装置のスループットを高めることができる。また画像データは適切な解像度に抑えられるため、表示部１１２への表示や画像処理装置１１０上での編集、記憶装置１１３への保存において応答速度を高め、ユーザに使い勝手の良い環境を提示できる。
【００７０】
なお、本実施例では、染色方法毎の撮像制御情報として、染色方法とビニング設定とが対応付けられたテーブルを用いたが、撮像制御情報はこのテーブルの例に限られない。例えば、各染色方法のビニング設定（統合される縦、横それぞれの画素数）は任意に設定できる。また、撮像制御情報のテーブルに、ビニング設定自体ではなく、相当倍率を設定することもできる。相当倍率とは、光学系の撮像倍率にビニング処理による一辺の画素数の縮小率を掛けた値である。相当倍率のテーブルの一例を以下に示す。
【表２】
0005372068.tif 000003
【００７１】
例えば、結像光学系２０５の倍率が４０倍であるときに、Ｓ１１０４において染色方法が「ＨＥ染色」であると推定された場合、一辺の画素数の縮小率は２０／４０＝１／２となり、ビニング数は２×２と求まる。このように相当倍率から求めたビニング数に従って、ＣＣＤイメージセンサの出力に対しビニング処理を施すことで、上記と同様の結果を得ることができる。
【００７２】
撮像制御情報として相当倍率を用いる方法は、結像光学系２０５の倍率が変更可能な構成の場合に好ましく利用できる。例えば、撮像倍率が４０倍、相当倍率が２０倍の場合は、ビニング数は２×２となり、撮像倍率が６０倍、相当倍率が２０倍の場合は、ビニング数は３×３となる。このように、用いられる撮像倍率に合わせてビニング数が適応的に決定されるため、実際の撮像倍率にかかわらず、常に一定の解像度（データ量）の撮像画像が得られるという利点がある。
【００７３】
［実施例２］
実施例２では、ＣＭＯＳイメージセンサを用いたビニング処理によりデータ量を削減する例を説明する。なお、本実施例におけるＣＭＯＳイメージセンサは画素毎に増幅した電圧信号が読み出し可能なＡＰＳ（Ａｃｔｉｖｅ  Ｐｉｘｅｌ  Ｓｅｎｓｏｒ）型であるとし、以降の説明を行う。
【００７４】
ＣＣＤイメージセンサではバケツリレー式に電荷を移動し、最後にまとめて増幅するため、電荷移動時に電荷の統合が可能である。しかし、ＣＭＯＳイメージセンサの場合、各
画素に蓄積された電荷を、電圧に変換・増幅して画素毎に読み出している。そのため、電荷移動時にビニングを行うことはできない。そこで本実施例では、ＣＭＯＳイメージセンサの後段に解像度変換処理（画素統合処理）を行う演算回路を実装し、ＣＭＯＳイメージセンサの出力信号に対してデータ演算を行うことで、ビニング処理を実現する。
【００７５】
図１２はＣＭＯＳイメージセンサ１６００に接続されたデータ演算回路１６０３を用いてビニング処理を実現する例を示している。
ＣＭＯＳイメージセンサ１６００内には受光素子１６０１と読み出し回路１６０２が存在する。受光素子１６０１の各画素はフォトダイオードからなる受光素子で、それぞれに読み出し回路１６０２が接続されており、蓄積電荷を電圧で読みだす構成となっている。
【００７６】
読み出し回路１６０２は、ＣＤＳ（相関二重サンプリング）によるノイズ低減処理回路およびＡ／Ｄ変換回路を有しており、読み出し回路１６０２の出力はデジタル信号である。
データ演算回路１６０３は、指定されたビニング設定に応じて、統合する画素内のデジタル信号の加算を行うことでビニング処理を実行する。ビニング処理の結果はデータ演算回路１６０３に接続するメモリ１６０４に書きだされる。
【００７７】
例えば、図１２のＣＭＯＳセンサ１６００に一様に光があたり、受光素子１６０１の各画素から電圧値ｄ１［Ｖ］が出力され、読み出し回路１６０２においてデジタル信号値ｄ２に変換されたとする。ビニング設定が２×２の場合、データ演算回路１６０３はハッチングで示した４つの画素の信号値を加算し、加算結果の値４×ｄ２を画素統合後の画素の位置に対応するメモリ位置に結果が書き込む。
【００７８】
上記の処理により、ＣＣＤイメージセンサ以外のセンサにおいてもビニング処理が可能であり、実施例１同様に、推測した染色方法を元にビニング設定を行うことで、データ量の削減が可能である。
以上述べた実施例２の構成によっても、実施例１と同様の効果が実現できる。すなわち、複数の染色方法を施したプレパラートを混在させてバッチ処理を行った場合においても、プレ計測において染色方法を推定し、それによって本計測で取得するデータ量を削減できる。そのため、過剰な解像度でデータを撮像する必要がなく、高いスループットが実現できる。また画像データは適切な解像度に抑えられているため、表示や編集、保存において応答速度を高め、ユーザに使い勝手の良い環境を提示できる。
【００７９】
なお、本実施例では画像データの解像度変換をデータ演算回路１６０３で実現したが、同様の処理を画像処理装置１１０内にあるＣＰＵ４０１および外部記憶装置に格納したプログラムを用いて実現することも可能であり、そうした実施形も本発明の範疇である。
【００８０】
［実施例３］
実施例３では、実施例１、２とは異なり、図４（ａ）に示した本計測処理の内部フローの圧縮処理Ｓ５０４において染色方法に応じて画像圧縮パラメータを変更することでデータ量を削減する例を説明する。
【００８１】
本実施例では圧縮処理Ｓ５０４で用いられる画像圧縮方式の一例として、ＪＰＥＧ圧縮符号化方式について説明する。なお、ＪＰＥＧ２０００やＪＰＥＧ−ＸＲなどの他の圧縮符号化方式でも画像変換処理、量子化処理を行う点ではＪＰＥＧと類似の構成を取るため、本実施例を適用できる。（例えば、ＪＰＥＧ２０００では離散コサイン変換ではなくウェーブレット変換を用いる。）
【００８２】
まず初めにＪＰＥＧ圧縮方式の概要について簡単に説明を行う。
図１３はＪＰＥＧ圧縮符号化方式の符号化および復号化装置のブロック図を表している。図１３の上段はＪＰＥＧ画像符号化装置のブロック図である。
【００８３】
画像入力部１７０１に対して符号化対象となる画像を構成する画素信号が、最小符号化単位ごとに入力され、その出力は色変換部１７０２に入力される。色変換部１７０２は、入力した画像信号に対し、カラーの場合にはＹＣｒＣｂ（輝度色差）などの圧縮に適した色空間へ、色変換処理を行う。色変換に関しては、規定および任意の色変換式（色変換アルゴリズム）を用いることができる。このとき、圧縮時に用いた色変換式の情報を圧縮符号化データに付与することで、復号化装置で復元（逆色変換）が可能となる。色変換部１７０２の出力は離散コサイン変換部１７０３に入力される。
離散コサイン変換部１７０３は、入力した画像信号に対して、２次元離散コサイン変換処理を行い、離散コサイン変換係数を計算して出力する。離散コサイン変換の詳細は公知であるのでここでは説明を省略する。
【００８４】
次に量子化部１７０４は、入力した係数を量子化テーブル１７１３に基づいて量子化し、その量子化値に対するインデックスを出力する。量子化テーブル１７１３は、画質と圧縮率を決める要素である。量子化部１７０４に入力された変換係数は、ＤＣ成分は差分を計算した上で、ジグザグ・スキャンされ、後続のエントロピ符号化部１７０５に出力される。
エントロピ符号化部１７０５は、符号化テーブルに従い、入力された量子化後の離散コサイン変換係数を符号化し、ビットストリームを出力する。符号出力部１７０６は、そのビットストリームを含む符号列を伝送路を介して符号入力部１７０７へ出力可能である。
【００８５】
次に、上記のビットストリームを復号化する復号化装置について説明する。
図１３の下段はＪＰＥＧ画像復号化装置の構成を表すブロック図であり、１７０７が符号入力部、１７０８はエントロピ復号化部、１７０９は逆量子化部、１７１０は逆離散コサイン変換部、１７１１は逆色変換部、１７１２は画像出力部である。
【００８６】
符号入力部１７０７は、符号列を入力し、それに含まれるヘッダを解析して後続の処理に必要なパラメータを抽出し、必要な場合は処理の流れを制御し、あるいは後続の処理ユニットに対して該当するパラメータを送出するものである。
また、符号列に含まれるビットストリームはエントロピ復号化部１７０８に出力される。エントロピ復号化部１７０８はビットストリームを復号化し、出力する。復元された量子化インデックスは逆量子化部１７０９に出力される。
【００８７】
逆量子化部１７０９は入力した量子化インデックスおよびヘッダから読み込まれた量子化テーブル１７１３から、離散コサイン変換係数を復元する。変換係数は後続の逆離散コサイン変換部１７１０に出力される。
逆離散コサイン変換部１７１０は、２次元の逆離散コサイン変換を実行し、変換係数から最小符号化単位の色変換された原画像の信号を計算する。逆色変換部１７１１は入力した最小符号化単位の色変換された原画像の信号に対し、逆の色変換を行い、最小符号化単位の原画像の信号を復元する。
原画像を構成する全ての最小符号化単位について上記の処理を行い、原画像全体を復元する。画像出力部１７１２は、復元された画像信号を出力する。
【００８８】
（撮像データ量の削減）
次に、本実施例の特徴である、染色方法に応じた圧縮パラメータの制御について述べる。
前述したように組織診と細胞診を比較すると、組織診の方が細胞診に比べて解像度が重視されない傾向がある。そのため、ＨＥ染色などの組織診で用いる染色方法で染色された
標本では、パパニコロウ染色やギムザ染色などの細胞診で用いる染色方法に比べて、量子化テーブル内の高周波成分の量子化ステップを粗く設定すると良い。それにより本計測で撮像する画像データのデータ量（圧縮後の画像データのサイズ）を削減することができる。
【００８９】
また、染色された画像では偏った色の分布が起きる為、染色方法毎に色変換式を最適化することで、偏った色成分に多くのビット数を割り当て、その他の成分に少ないビット数を割り当てることで画像全体を効率的に圧縮することができる。
【００９０】
以下に、公知の技術を用いてデータ圧縮率を高める色変換式を決定する方法の一例を述べる。
ＪＰＥＧ圧縮符号化で良く用いられるＹＣｂＣｒ変換のＣｂＣｒ平面において、黄色（Ｙｅ）および青（Ｂ）はＣｂ軸に比較的近い位置に出現する。そのため、黄色（Ｙｅ）および青（Ｂ）の出現頻度が多い画像では、離散コサイン変換後のＣｒ軸に対応する空間周波数の量子化テーブルを粗くすることで画質を維持して画像圧縮を行える。一方、マゼンタ（Ｍｇ）および緑（Ｇ）はそれぞれＣｂＣｒ平面の第１象現および第３象現の中心位置に出現する。そのため、マゼンタ（Ｍｇ）および緑（Ｇ）の出現頻度の多い画像では、画質を維持しつつ、離散コサイン後のＣｂ軸、Ｃｒ軸に対応する空間周波数の量子化テーブルを粗くするには限界がある。
【００９１】
これに対応するには、出現頻度が多い成分が軸上にのるよう、Ｃｂ軸及びＣｒ軸を変換した新たなＣｂ´軸及びＣｒ´軸を求め、色の座標をＣｂ´Ｃｒ´平面で表現すると良い。これにより、出現座標をＣｂ´またはＣｒ´の軸のいずれか一方に偏らせることができる。そして離散コサイン変換後にＣｂ´またはＣｒ´軸に対応する空間周波数の量子化テーブルを粗くすることで画質を維持しつつ、画像圧縮の効率を上げられる。
【００９２】
色変換部１７０２におけるＲＧＢ色空間からＣ１Ｃ２Ｃ３色空間に対する色変換は以下のアフィン変換（一次変換と平行移動からなる）で表されるとする。（lt1、lt2、lt3は
平行移動成分）
【数１】
0005372068.tif 000004
【００９３】
Ｃ１Ｃ２Ｃ３色空間からＣ１´Ｃ２´Ｃ３´色空間への変換式はアフィン変換で表されるとすると、以下の式で表現できる。（mt1、mt2、mt3は平行移動成分）
【数２】
0005372068.tif 000005
【００９４】
従って、ＲＧＢ色空間からＣ１´Ｃ２´Ｃ３´色空間への変換式は以下で表される。
【数３】
0005372068.tif 000006
【００９５】
逆色変換部１７１１でＣ１´Ｃ２´Ｃ３´色空間からＲＧＢ色空間での値を求めるには、色変換部１７０２で用いた行列の逆行列を計算し、符号出力部１７０６が出力する符号化データに逆行列の情報を記載し、逆色変換部１７１１でその情報を用いればよい。
【００９６】
例えば、前述のマゼンタ（Ｍｇ）および緑（Ｇ）の出現頻度が高い染色方法において、ＹＣｂＣｒ色空間から、よりデータ量を圧縮できる色変換式を求めるには、まず、その染色が施された多数の検体の画像データから、２つの出現頻度の高いＣｂＣｒ座標を抽出する。そして、その座標を通る直線を求め、軸と直線の回転角θを求める。そして回転角θからアフィン変換行列を計算することで、ＲＧＢ色空間からＹＣｂ´Ｃｒ´への変換式が求められる。
【００９７】
なお、ここではＹＣｂＣｒ色空間を例に、データ圧縮率を高める色空間変換の求め方を説明したが、本実施例はその他の色空間にも適用できる。そのため、ｓＲＧＢ、ＡｄｏｂｅＲＧＢ、ＸＹＺ、Ｌ＊ａ＊ｂ＊、ＨＳＶなどの色空間を対象として色空間変換を行う方法も本発明の範疇とする。
【００９８】
以上説明したように、染色方法毎の撮像制御情報として、染色方法毎に対応付けて色変換式および量子化テーブルを保持し、推定した染色方毎に対応する値を用いて画像データの圧縮符号化を行うことで、本計測のデータ量を削減できる。
【００９９】
以下、本実施例において実施例１，２と異なる部分を説明する。本実施例では、染色方法毎の撮像制御情報取得処理Ｓ１４０２において、上記の染色方法毎の撮像制御情報から、推定した染色方法に対応する色変換式および量子化テーブルを取得する。そして、撮像データ削減設定処理Ｓ１４０３において、本計測処理における圧縮処理Ｓ５０４の画像圧縮パラメータとして、取得した色変換式、量子化テーブルをそれぞれ色変換部１７０２および量子化部１７０４で用いるパラメータとして設定する。そして本計測処理Ｓ９０４において設定した色変換式および量子化テーブルを用いて撮像した画像データを圧縮する。これにより伝送処理Ｓ５０５で画像処理装置に送付される画像データのサイズが小さくなるため、装置のスループットが向上する。また、本実施例で述べた方法は、実施例１，２で述べたビニングと異なり、量子化テーブルの設定しだいで、より細かなデータ削減ができるというメリットがある。
【０１００】
以上述べたように、実施例３においては、複数の染色方法を施したプレパラートを混在させてバッチ処理を行った場合においても、プレ計測において染色方法を推定し、それに応じて画像圧縮処理の設定を選ぶことにより、本計測で取得するデータ量を削減できる。そのため、過剰な解像度でのデータ取得を回避し、装置のスループットを高めることができる。また画像データのサイズが抑えられるため、記憶装置１１３への保存や他のコンピュータシステム１１４への送信レスポンスを高め、ユーザに使い勝手の良い環境を提示できる。
【０１０１】
なお、本実施例では、染色方法に応じて圧縮率（画質）を異ならせるという構成を採用
したが、本発明の構成はこれに限られない。例えば、データ量削減のために組織診用の染色の画像データのみを圧縮し、それ以外の画像データについては画質の劣化を防ぐために非圧縮（若しくは可逆圧縮）としてもよい。
【０１０２】
［実施例４］
実施例４においては、実施例１〜３で述べた撮像データ削減設定処理の前段に撮像倍率変更処理を加えることで、後段のデータ削減の効果を高める効果を実現する。
【０１０３】
図１４は実施例４における撮像条件設定処理Ｓ９０３の詳細を示すフローチャートである。実施例４の特徴は、染色方法毎の撮像制御情報が、実施例１〜３で述べた情報（ビニング設定、圧縮設定など）に加え、結像光学系２０５の対物レンズの倍率の情報を含む点、ならびに、撮像倍率変更処理Ｓ１８００を有する点である。下記のテーブルは、撮像制御情報に含まれる対物レンズ倍率情報の例を示している。前述のように細胞診よりも組織診の方が解像度が重視されない傾向にあるため、組織診で用いられるＨＥ染色の方が細胞診で用いられる染色方法に比べて撮像倍率が低い値に設定されている。それ以外の処理Ｓ１４０１〜Ｓ１４０４については、図１０で説明した実施例１〜３の処理と実質的に同じであるため、詳しい説明を省略する。
【表３】
0005372068.tif 000007
【０１０４】
図１５は結像光学系２０５の物体側の構造を示した模式図である。結像光学系２０５の物体側には異なる倍率を持つ複数の対物レンズ（図１５では４０倍のレンズ１９０１と２０倍のレンズ１９０２）が取り付けられたリボルバ１９００が装着されている。本実施例のバーチャル・スライド・システムにおいては画像処理装置１１０の指示をコントローラ１０８を介して不図示のリボルバ制御部に伝えることで、リボルバを回転させ、対物レンズの倍率を変更する。
【０１０５】
なお、対物レンズを交換する際、焦点位置が移動する場合があるが、焦点調整は本計測処理Ｓ９０４においてＡＦユニット１０９を用いて行い、コントローラ１０８はステージ制御１０５を通じて観察に適切な合焦位置にステージを移動するとする。
【０１０６】
以上述べた実施例４においても、実施例１〜３と同様、染色方法に応じて本計測の撮像データのデータ量を適切に削減できるという効果が得られる。さらに本実施例では、撮像倍率変更処理Ｓ１８０３を導入したことにより、次のような効果が得られる。
【０１０７】
１つは、デジタル画像の見えを、ユーザである病理医が見慣れている光学顕微鏡の観察像の見えにできる限り近づけることができるという点である。像の見えに影響するパラメータの典型例は、被写界深度やぼけ方である。光学顕微鏡を用いた診断とバーチャル・スライド・システムを用いた診断を併用するユーザの場合、被写界深度やぼけ方は光学顕微鏡の光学性能と類似している方が両者で見た画像を比較しやすい。そのため、光学顕微鏡で診断を行う際の倍率に合わせて、撮像制御情報における染色方法毎の撮像倍率を設定しておくことで、バーチャル・スライド・システムでの診断の利便性を向上することができる。
【０１０８】
また、ＨＥ染色の場合に低倍で撮像することにより、データ削減性能を高めることができる。例えば、レンズ倍率を低くするほど画像の高周波成分が低下するため、画像データ
圧縮処理におけるデータ削減性能（圧縮率）を高めることができる。なお、画像の高周波成分を低下させることで、実施例１、２で述べたビニング処理において画素間隔のサンプリングによって発生するエイリアシング（折り返し歪）を抑制する効果も得られる。
【０１０９】
［実施例５］
実施例１〜４では、プレ計測で得られた画像を用いて染色方法を推定したが、他の情報を用いて染色方法を推定することもできる。実施例５では、プレ計測で撮像した画像以外の画像から染色方法を推定し、データ量を削減する方法について述べる。なお本実施例の方法を用いる場合には、プレ計測ユニットは不要となる。
【０１１０】
図１６は本実施例における染色方法推定制御処理の内部処理フローを示すフローチャートである。以下、図２０の各ステップについて説明する。
画像取得処理Ｓ２００１では、画像処理装置１１０が、染色方法の推定に用いる画像を取得する。取得方法は、本計測ユニット１０１から撮像データを取得する方法と、既に他の撮像装置で撮像され所定の画像フォーマット形式で格納されている画像データファイルを取得する方法の２つがあり、いずれの方法を用いてもよい。後者の場合、画像処理装置１１０は、例えば図１に示す他のコンピュータシステム１１４からネットワークを通じて取得することもできるし、記憶装置１１３からファイルを読み込むこともできる。
【０１１１】
次に、染色方法推定処理Ｓ２００２では、画像処理装置１１０が、Ｓ２００１で取得した画像データを解析し、画像データ中の検体に施された染色方法を推定する。具体的な推定処理は実施例１〜４で既述のものと同様である。
【０１１２】
Ｓ２００１で取得した画像データが、本計測ユニット１０１で撮像された画像データのように高精細でデータサイズの大きい画像データの場合には、オリジナルサイズではなく、低解像度の縮小画像のデータを用いて染色方法の推定処理を行うとよい。染色方法を推定する目的であればそれほど高精細の画像データは必要ないため、データ量の削減により処理速度の向上を図るほうが望ましいからである。
【０１１３】
具体的には、解像度変換により低解像度の縮小画像を作成すればよい。解像度変換の方法としては、ｍ×ｍの画素ブロック単位で平均化することで１／ｍのサイズの低解像度画像を作成する方法や、バイキュービック、バイリニア、ニアレストネイバー等の補間を利用して低解像度画像を作成する方法など、いかなる方法を用いることもできる。また、一部の画像ファイルフォーマットでは、高速表示のため、縮小画像（サムネール画像）や元の画像の半分や１／４のサイズなど複数のサイズの画像データを格納する構造を備えている。そのようなファイルフォーマットの画像の場合、画像処理装置１１０はファイル内に存在する低解像度の縮小画像を抽出し、染色方法の推定に用いるとよい。
【０１１４】
なお、染色方法推定処理Ｓ２００２の処理はネットワーク上のコンピュータで行っても良い。その場合、画像処理装置１１０で低解像度の画像データを計算あるいは抽出し、その画像をネットワーク上のコンピュータに送る。続いて、ネットワーク上のコンピュータは染色方法を推定し、推定結果を画像処理装置１１０に送る。
また、画像処理装置１１０がヒストグラム生成処理Ｓ１１０２（図７（ｂ））まで実行してその結果をネットワーク上のコンピュータに送り、ネットワーク上のコンピュータがマッチング度算出処理Ｓ１１０３以降の処理を実行してもよい。この方法では、画像データそのものを伝送する必要がないため、通信時間を短縮することができる。
【０１１５】
次に、データ削減処理Ｓ２００３では、画像処理装置１１０は、染色方法推定処理Ｓ２００２で推定した染色方法に従って、画像データのデータ量を削減する。
【０１１６】
画像取得処理Ｓ２００１において本計測ユニット１０１を用いて撮像した画像データを入力する場合、実施例１〜３と類似の構成を取ることが出来る。すなわち、画像処理装置１１０は、染色方法毎に設定されたビニング設定や色変換式と量子化テーブルの組み合わせを参照し、その設定に従って（必要に応じて）Ｓ２００１で取得した画像データの解像度変換や圧縮処理を実行する。例えば、組織診用には、光学倍率２０倍相当の画像データで十分であったとする。Ｓ２００１で取得した画像データの撮像倍率が４０倍であり、且つ、Ｓ２００２で染色方法がＨＥ染色と推定された場合は、画像処理装置１１０がデジタル画像処理で２×２画素を１画素に統合（実施例１、２のビニング処理に相当）することで、データ量を削減すればよい。
【０１１７】
一方、画像取得処理Ｓ２００１において、記憶装置１１３又は他のシステム１１４から画像データを読み込む場合には、画像処理装置１１０は画像データのファイルフォーマットに応じてデータ量の削減方法を決定するとよい。
【０１１８】
例えば、複数のサイズの画像データを保有しているファイルフォーマットの画像の場合は、そのファイル内の一部のデータを削除することで、同じようなデータ量削減効果を得ることができる。例えば、画像データファイルの中に、等倍画像（倍率４０倍相当）、１／２倍画像（倍率２０倍相当）、１／４倍画像（倍率１０倍相当）、１／１０倍画像（倍率４倍相当）の４種類の解像度の画像が含まれているとする。染色方法がＨＥ染色と推定された場合は、不必要な解像度の画像データ（この例では等倍画像データ）を削除することで、大幅なデータ量削減が可能となる。この方法は、解像度変換が不要なため、計算負荷がかからないという利点がある。
【０１１９】
また他の一例として、画像データファイルがＪＰＥＧ２０００圧縮符号化方式の場合を考える。ＪＰＥＧ２０００では、圧縮符号化時に画像データに対して離散ウェーブレット変換を施し、画像を４つの帯域（ＬＬ、ＬＨ、ＨＬ、ＨＨ）に分解し、最低域のＬＬの画像を同様の操作で再帰的に分解している。反対に、復号時は最低域のＬＬから再帰的に合成することで徐々に低解像度画像から高解像度画像が復元できる仕組みとなっている。そのため、元の画像データが４０倍相当の解像度であり、２０倍相当の解像度に削減する場合にはＬＨ、ＨＬ，ＨＨの帯域の符号化データを削除する。これにより最低域ＬＬの画質劣化を伴うことなく、効率的にデータ量を削減できる。一部のデータを削除した後、ファイルの整合性が取れるように、ヘッダ情報等のファイル構成を表す情報を更新する。
【０１２０】
このように所定の画像フォーマット形式で格納された画像データファイルであっても、染色方法と画像フォーマット形式に応じて、効率的にデータ削減ができる。
【０１２１】
以上、実施例１〜実施例５で述べたように、染色方法の推定結果に応じて画像データのデータ量を適宜削減することにより、過剰な解像度でのデータ取得を回避し、装置のスループット並びにデータ伝送・保存・編集・閲覧の応答を高め、ユーザの利便性向上を図ることができる。
【０１２２】
なお、実施例１〜実施例５の構成は本発明の具体例にすぎず、本発明はこれらの実施例の構成以外にも様々な構成を採用し得る。例えば、撮像制御情報を格納するデータベースは、画像処理装置１１０の外部記憶装置ではなく、ネットワークサーバ上にあってもよい。また、撮像制御情報は拠点（病院、企業、学術機関、分析センター、その他の施設など）ごとに同じ値に設定することもできるし、ユーザである医師がバッチ処理毎に撮像制御情報を設定してもよい。
また実施例１では、ＣＣＤイメージセンサの各画素に蓄積された電荷を読み出す際に、電荷を加算するビニング処理を行い、データ量を削減する例を説明したが、画素単位の増幅機能を持たないパッシブ型のＣＭＯＳセンサにおいても同様の構成が可能である。その
ため、ＣＣＤ、ＣＭＯＳセンサの種類に関わらず、染色方法に応じてビニング設定を変えることは本発明の範疇とする。
また各実施例で説明した構成を互いに組み合わせることで、データ削減効果をさらに高めることもできる。例えば、実施例１、２のビニング処理と実施例３の圧縮処理を組み合わせたり、実施例３の圧縮処理と実施例４の撮像倍率変更を組み合わせたりするとよい。
【符号の説明】
【０１２３】
１０１：本計測ユニット、１０２：プレ計測ユニット、１０８：コントローラ、１１０：画像処理装置
【図１】
0005372068.tif 000008
【図２】
0005372068.tif 000009
【図３】
0005372068.tif 000010
【図４】
0005372068.tif 000011
【図５】
0005372068.tif 000012
【図６】
0005372068.tif 000013
【図７】
0005372068.tif 000014
【図８】
0005372068.tif 000015
【図９】
0005372068.tif 000016
【図１０】
0005372068.tif 000017
【図１１】
0005372068.tif 000018
【図１２】
0005372068.tif 000019
【図１３】
0005372068.tif 000020
【図１４】
0005372068.tif 000021
【図１５】
0005372068.tif 000022
【図１６】
0005372068.tif 000023