(57)【要約】【課題】所定の染色色素で染色された染色標本の色素量を精度良く推定すること。【解決手段】本発明のある実施の形態において、記憶部１６は、基準スペクトル情報１６３として、観察対象の染色標本を染色している染色色素について異なる染色状態での複数の基準スペクトルを記憶する。スペクトル取得部１４１は、染色標本を撮像した染色標本画像の画素毎に分光特性値を取得する。標本作成条件推定処理部１４２は、例えば染色標本を作成した際の作成条件を推定して染色標本の作成条件を取得する。基準スペクトル決定部１４６は、基準スペクトル情報１６３から染色標本の作成条件に応じた基準スペクトルを選択し、染色色素の最適基準スペクトルを決定する。色素量推定部１４７は、スペクトル取得部１４１が取得した分光特性値をもとに、染色色素の最適基準スペクトルを用いて染色標本を染色している染色色素の色素量を推定する。【選択図】図１
2011095225.tif 000002
【特許請求の範囲】
【請求項１】
所定の染色色素で染色された観察対象の染色標本を撮像した染色標本画像を用いて前記染色標本の色素量を推定する画像処理装置であって、
前記染色色素についての異なる染色状態での複数の色素分光特性値を記憶する色素分光特性記憶手段と、
前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得手段と、
前記染色標本の作成条件を取得する作成条件取得手段と、
前記色素分光特性記憶手段に記憶されている前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定手段と、
前記分光特性取得手段によって取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定手段と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記作成条件取得手段は、前記分光特性取得手段によって取得された分光特性値をもとに前記染色標本を作成した際の作成条件を推定する作成条件推定手段を有し、該作成条件推定手段によって推定された作成条件を前記染色標本の作成条件として取得することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記作成条件取得手段は、前記染色標本を作成した際の作成条件の入力を依頼する作成条件入力依頼手段を有し、前記作成条件入力依頼手段による依頼に対して入力された作成条件を前記染色標本の作成条件として取得することを特徴とする請求項１に記載の画像処理装置。
【請求項４】
前記作成条件推定手段は、
前記染色標本画像中に所定の解析領域を設定する解析領域設定手段と、
前記解析領域を構成する画素について前記分光特性取得手段が取得した分光特性値をもとに、前記解析領域の特徴量を取得する特徴量取得手段と、
を有し、
前記特徴量をもとに前記染色標本を作成した際の作成条件を推定することを特徴とする請求項２に記載の画像処理装置。
【請求項５】
前記解析領域は、前記染色標本画像中に映る前記染色標本内の核の領域であることを特徴とする請求項４に記載の画像処理装置。
【請求項６】
前記分光特性取得手段は、所定の波長毎または所定の波長帯域毎に前記分光特性値を取得し、
前記特徴量取得手段は、前記解析領域の画素の前記所定の波長毎または前記所定の波長帯域毎の前記分光特性値をグラフ化して分光特性グラフを作成し、該分光特性グラフのグラフ形状を解析して前記特徴量を算出することを特徴とする請求項４に記載の画像処理装置。
【請求項７】
前記特徴量取得手段は、前記分光特性グラフのグラフ形状が平坦な平坦波長区間および／または前記平坦波長区間における分光特性値平均と前記分光特性グラフのピーク波長における分光特性値との差分を前記特徴量として算出することを特徴とする請求項６に記載の画像処理装置。
【請求項８】
前記色素分光特性記憶手段は、前記染色色素についての異なる染色状態での複数の色素分光特性値として、前記染色色素で個別に染色された単一染色標本であって、作成条件の異なる複数の単一染色標本から取得した前記作成条件毎の色素分光特性値を記憶しており、
前記色素分光特性決定手段は、前記染色標本の作成条件をもとに、前記色素分光特性記憶手段に記憶されている作成条件毎の色素分光特性値の中から該当する色素分光特性値を選択することを特徴とする請求項１に記載の画像処理装置。
【請求項９】
前記作成条件推定手段は、２以上の作成条件によって定義される所定の作成条件空間内に前記作成条件毎の色素分光特性値それぞれについての作成条件判定パラメータを分布させた作成条件判定パラメータ分布を参照し、前記解析領域についての作成条件判定パラメータに該当する作成条件判定パラメータを選出することによって前記染色作成条件を作成した際の作成条件を推定することを特徴とする請求項８に記載の画像処理装置。
【請求項１０】
前記色素分光特性決定手段は、前記染色標本の作成条件が前記色素分光特性記憶手段に記憶されているいずれかの色素分光特性値の作成条件と一致している場合には、該一致している作成条件に対応する色素分光特性値を前記染色色素の最適色素分光特性値として決定することを特徴とする請求項１に記載の画像処理装置。
【請求項１１】
前記色素分光特性決定手段は、前記染色標本の作成条件が前記色素分光特性記憶手段に記憶されているいずれの色素分光特性値の作成条件とも一致していない場合には、該作成条件が前記染色標本の作成条件と最も類似する色素分光特性値を選択し、選択した色素分光特性値を該色素分光特性値の作成条件と前記染色標本の作成条件との差分に従って補正して前記染色色素の最適色素分光特性値を決定することを特徴とする請求項１に記載の画像処理装置。
【請求項１２】
表示対象色素を指定する色素選択手段と、
前記色素量推定手段によって推定された前記染色標本上の標本点の前記染色色素の色素量をもとに、前記表示対象色素を含む前記染色標本画像中の画素位置を他の画素位置と識別可能に表した表示画像を生成する表示画像生成手段と、
前記表示画像を表示部に表示処理する表示処理手段と、
を備えることを特徴とする請求項１に記載の画像処理装置。
【請求項１３】
前記色素量推定手段によって推定された前記染色標本上の標本点の前記染色色素の色素量を補正する色素量補正手段と、
前記色素量補正手段による補正後の前記染色色素の色素量をもとに分光特性値を合成する分光特性合成手段と、
前記分光特性合成手段によって合成された分光特性値をもとにＲＧＢ画像を合成し、前記染色色素による染色状態を補正後の前記染色色素の色素量で表した表示画像を生成する表示画像生成手段と、
前記表示画像を表示部に表示処理する表示処理手段と、
を備えることを特徴とする請求項１に記載の画像処理装置。
【請求項１４】
表示対象色素を指定する色素選択手段を備え、
前記色素量補正手段は、前記表示対象色素の色素量を補正することを特徴とする請求項１３に記載の画像処理装置。
【請求項１５】
表示対象色素を指定する色素選択手段を備え、
前記色素量補正手段は、前記表示対象色素以外の染色色素の色素量をゼロに補正し、
前記表示画像生成手段は、前記表示対象色素以外の染色色素を非表示として前記表示画像を生成することを特徴とする請求項１３に記載の画像処理装置。
【請求項１６】
所定の染色色素についての異なる染色状態での複数の色素分光特性値が記憶された画像処理装置において、前記染色色素で染色された観察対象の染色標本を撮像した染色標本画像を用いて前記染色標本の色素量を推定する画像処理方法であって、
前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得工程と、
前記染色標本の作成条件を取得する作成条件取得工程と、
前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定工程と、
前記分光特性取得工程で取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定工程と、
を含むことを特徴とする画像処理方法。
【請求項１７】
顕微鏡を用い、所定の染色色素で染色された観察対象の染色標本を撮像して染色標本画像を取得する画像取得手段と、
少なくとも前記染色色素についての異なる染色状態での複数の色素分光特性値を記憶する分光特性記憶手段と、
前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得手段と、
前記染色標本の作成条件を取得する作成条件取得手段と、
前記分光特性記憶手段に記憶されている前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定手段と、
前記分光特性取得手段によって取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定手段と、
を備えることを特徴とする顕微鏡システム。
【請求項１８】
前記分光特性記憶手段には、染色標本の属性を示す属性値毎に前記複数の色素分光特性値が記憶されており、
前記観察対象の染色標本の属性値についての前記複数の色素分光特性値をもとに、前記観察対象の染色標本を撮像する際の前記画像取得手段の動作環境を設定するためのシステムパラメータを設定するシステム環境設定部を備えることを特徴とする請求項１７に記載の顕微鏡システム。
【請求項１９】
前記分光特性記憶手段には、染色標本の属性を示す属性値毎に予め測定されたシステム分光特性値がさらに記憶されており、
前記観察対象の染色標本の属性値についての前記システム分光特性値をもとに、前記観察対象の染色標本を撮像する際の前記画像取得手段の動作環境を設定するためのシステムパラメータを設定するシステム環境設定部を備えることを特徴とする請求項１７に記載の顕微鏡システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、所定の染色色素で染色された染色標本を撮像した染色標本画像を処理する画像処理装置、画像処理方法および顕微鏡システムに関するものである。
【背景技術】
【０００２】
例えば病理診断では、臓器摘出や針生検によって得た組織検体を厚さ数ミクロン程度に薄切して標本を作成し、様々な所見を得るために顕微鏡を用いて拡大観察することが広く行われている。中でも光学顕微鏡を用いた透過観察は、機材が比較的安価で取り扱いが容易である上、歴史的に古くから行われてきたこともあって、最も普及している観察方法の一つである。ここで、生体から採取した検体は光をほとんど吸収および散乱せず、無色透明に近い。このため、標本の作成に際し、色素による染色を施すのが一般的である。
【０００３】
染色手法としては種々のものが提案されており、その総数は１００種類以上にも達するが、特に病理標本に関しては、青紫色のヘマトキシリン（以下、「色素Ｈ」と表記する。）および赤色のエオジン（以下、「色素Ｅ」と表記する。）の２つの色素を染色色素として用いるヘマトキシリン−エオジン染色（以下、「Ｈ＆Ｅ染色」と呼ぶ。）が標準的に用いられている。このＨ＆Ｅ染色後の標本（染色標本）では、細胞核や骨組織等が青紫色に、細胞質や結合織、赤血球等が赤色に染色され、これらが容易に視認できる。この結果、観察者は、細胞核等の組織を構成する要素の大きさや位置関係等を把握でき、標本の状態を形態学的に判断することが可能となる。
【０００４】
ところで、染色標本の観察は、目視によって行われる他、染色標本を撮像して表示装置に画面表示することによっても行なわれており、染色標本画像を画像処理して解析し、医師等による観察・診断を支援しようとする試みが提案されている。その一例として、従来から、染色標本をマルチバンド撮像して得た染色標本画像をもとに染色標本上の点（標本点）を染色している染色色素の色素量を定量推定する手法が知られており、様々な目的に応用されている。例えば、非特許文献１には、色素量を推定し、推定した色素量をもとに染色標本画像の色情報を補正する手法が開示されている。また、非特許文献２には、推定した色素量をもとに標本の染色状態を定量評価する手法が開示されている。また、特許文献１には、推定した色素量をもとに標本内の組織を分類し、画像を組織毎の領域に領域分割する手法が開示されている。
【０００５】
ここで、染色標本画像（マルチバンド画像）から色素量を定量推定する方法について、Ｈ＆Ｅ染色された染色標本を例に挙げて説明する。色素量の推定に先立ち、例えば照明光を照射した状態で標本なしの背景を撮像して背景（照明光）のマルチバンド画像を取得しておく。先ず、背景のマルチバンド画像をＩ0とし、観察対象の染色標本のマルチバンド画像をＩとして、次式（１）に従って各画素位置の分光透過率ｔ（ｘ，λ）を算出する。ｘはマルチバンド画像の画素を表す位置ベクトルであり、λは波長である。また、Ｉ（ｘ，λ）はマルチバンド画像Ｉの波長λにおける画素位置（ｘ）の画素値を表し、Ｉ0（ｘ，λ）はマルチバンド画像Ｉ0の波長λにおける画素位置（ｘ）の画素値を表す。
【数１】
2011095225.tif 000003
【０００６】
分光透過率ｔ（ｘ，λ）に関しては、ランベルト・ベール（Lambert-Beer）の法則が成り立つ。例えば、染色標本が色素Ｈおよび色素Ｅの２種類の染色色素で染色されている場合、ランベルト・ベールの法則により、各波長λにおいて次式（２）が成り立つ。
【数２】
2011095225.tif 000004
【０００７】
式（２）において、ｋH（λ）およびｋE（λ）は、波長λに依存して決まる物質固有の係数であり、ｋH（λ）は色素Ｈに対応する係数、ｋE（λ）は色素Ｅに対応する係数である。例えば、ｋH（λ），ｋE（λ）の各値は、染色標本を染色している色素Ｈおよび色素Ｅの色素分光特性値（以下、染色標本を染色している染色色素の色素分光特性値を「基準スペクトル」と呼ぶ。）である。また、ｄH（ｘ），ｄE（ｘ）は、マルチバンド画像の各画素位置（ｘ）に対応する染色標本の各標本点における色素Ｈおよび色素Ｅの色素量に相当する。より詳細には、ｄH（ｘ）は、色素Ｈのみで染色された染色標本における色素Ｈの色素量を「１」としたときのこの色素量に対する相対的な値として求められる。同様に、ｄE（ｘ）は、色素Ｅのみで染色された染色標本における色素Ｅの色素量を「１」としたときのこの色素量に対する相対的な値として求められる。なお、色素量は、濃度とも呼ばれる。
【０００８】
ここで、上記式（２）は、波長λ毎に独立して成り立つ。また、式（２）はｄH（ｘ），ｄE（ｘ）の線形式であり、これを解く手法は一般には重回帰分析として知られている。例えば、２つ以上の異なる波長について式（２）を連立させれば、これらを解くことができる。
【０００９】
例えば、Ｍ個（Ｍ≧２）の波長λ1，λ2，・・・，λMについて式を連立させると、次式（３）のように表記できる。［］tは転置行列、［］-tは逆行列を表す。
【数３】
2011095225.tif 000005
【００１０】
そして、最小二乗推定を用いて上記式（３）を解くと次式（４）が得られ、色素Ｈの色素量の推定値ｄ＾H（ｘ）および色素Ｅの色素量の推定値ｄ＾E（ｘ）がそれぞれ求まる。なお、ｄ＾は、ｄの上に推定値を表す記号「＾（ハット）」が付いていることを示す。
【数４】
2011095225.tif 000006
【００１１】
そして、この式（４）によって、染色標本上の任意の標本点における色素Ｈおよび色素Ｅの色素量の推定値が得られる。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００５−３３１３９４号公報
【非特許文献】
【００１３】
【非特許文献１】“Color  Correction  of  Pathological  Images  Based  on  Dye  Amount  Quantification”，OPTICAL  REVIEW  Vol.12，No.4(2005)，p.293-300
【非特許文献２】“Development  of  support  systems  for  pathology  using  spectral  transmittance  -  The  quantification  method  of  stain  conditions”，Proceedings  of  SPIE  -  Image  Processing  Vol.4684(2002)，p.1516-1523
【発明の概要】
【発明が解決しようとする課題】
【００１４】
ところで、上記したｋH（λ）やｋE（λ）といった色素量の推定に用いる染色色素の基準スペクトルは、例えば、その染色色素を用いてそれぞれ個別に染色した標本（以下、「単一染色標本」と呼ぶ。）を予め用意し、測定機器等を用いてその分光特性値（スペクトル）を測定することによって得られる。しかしながら、実際に染色標本を染色している染色色素のスペクトルは、観察対象とする染色標本によって変動する場合がある。このため、観察対象の染色標本を染色している染色色素のスペクトルが事前に取得された基準スペクトルと異なる場合が発生してしまい、色素量の推定精度が低下するという問題があった。
【００１５】
ここで、染色色素のスペクトルは、染色標本の染色状態（濃い／薄いといった染色の程度）によって変動することが知られている。この染色状態のばらつきは、例えば染色に要する時間（染色時間）や標本の厚みといった染色標本の作成条件によって生じるが、標本を作成する臨床検査技師の経験等の個人差が影響するため、常に均一に保つのは難しい。一方で、染色標本は医師によって観察・診断されるため、医師の好みに応じた染色状態が求められる場合もある。また、染色色素のスペクトルは、薬剤（染色液）の水素イオン指数（ｐＨ）によって変動することも知られている。このｐＨ等の薬剤の状態は、標本を染色する施設によって変化し得る。同じ施設であっても、使用する薬剤を開封してからの時間に応じて変化し得る。
【００１６】
本発明は、上記した従来の問題点に鑑みて為されたものであり、所定の染色色素で染色された染色標本の色素量を精度良く推定することができる画像処理装置、画像処理方法および顕微鏡システムを提供することを目的とする。
【課題を解決するための手段】
【００１７】
上記した課題を解決し、目的を達成するための、本発明のある態様にかかる画像処理装置は、所定の染色色素で染色された観察対象の染色標本を撮像した染色標本画像を用いて前記染色標本の色素量を推定する画像処理装置であって、前記染色色素についての異なる染色状態での複数の色素分光特性値を記憶する色素分光特性記憶手段と、前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得手段と、前記染色標本の作成条件を取得する作成条件取得手段と、前記色素分光特性記憶手段に記憶されている前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定手段と、前記分光特性取得手段によって取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定手段と、を備えることを特徴とする。
【００１８】
この態様にかかる画像処理装置によれば、観察対象の染色標本を染色している染色色素について異なる染色状態での複数の色素分光特性値を記憶しておき、これら複数の色素分光特性値の中から染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、染色色素の最適色素分光特性値を決定することができる。そして、染色標本画像を構成する各画素について取得した分光特性値をもとに、決定した染色色素の最適色素分光特性値を用いて染色標本上の標本点における色素量を推定することができる。したがって、観察対象の染色標本の色素量を精度良く推定することができる。
【００１９】
また、本発明の別の態様にかかる画像処理方法は、所定の染色色素についての異なる染色状態での複数の色素分光特性値が記憶された画像処理装置において、前記染色色素で染色された観察対象の染色標本を撮像した染色標本画像を用いて前記染色標本の色素量を推定する画像処理方法であって、前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得工程と、前記染色標本の作成条件を取得する作成条件取得工程と、前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定工程と、前記分光特性取得工程で取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定工程と、を含むことを特徴とする。
【００２０】
また、本発明の別の態様にかかる顕微鏡システムは、顕微鏡を用い、所定の染色色素で染色された観察対象の染色標本を撮像して染色標本画像を取得する画像取得手段と、少なくとも前記染色色素についての異なる染色状態での複数の色素分光特性値を記憶する分光特性記憶手段と、前記染色標本画像を構成する画素の画素値をもとに、該画素に対応する前記染色標本上の標本点における分光特性値を取得する分光特性取得手段と、前記染色標本の作成条件を取得する作成条件取得手段と、前記分光特性記憶手段に記憶されている前記複数の色素分光特性値の中から前記染色標本の作成条件に応じた染色状態での色素分光特性値を選択し、前記染色色素の最適色素分光特性値を決定する色素分光特性決定手段と、前記分光特性取得手段によって取得された分光特性値をもとに、前記染色色素の最適色素分光特性値を用いて前記染色標本上の標本点における前記染色色素の色素量を推定する色素量推定手段と、を備えることを特徴とする。
【発明の効果】
【００２１】
本発明によれば、所定の染色色素で染色された染色標本の色素量を精度良く推定することができる。
【図面の簡単な説明】
【００２２】
【図１】図１は、実施の形態１における画像処理装置の機能構成を示すブロック図である。
【図２】図２は、作成条件が同じ１組の色素Ｈおよび色素Ｅの基準スペクトルグラフを示す図である。
【図３】図３は、作成条件判定パラメータの決定方法を説明する説明図である。
【図４】図４は、部分基準スペクトルグラフを示す図である。
【図５】図５は、作成条件判定パラメータ分布の一例を示す図である。
【図６】図６は、実施の形態１の画像処理装置が行う処理手順を示すフローチャートである。
【図７】図７は、標本作成条件推定処理の詳細な処理手順を示すフローチャートである。
【図８】図８は、解析領域選択画面の一例を示す図である。
【図９】図９は、解析領域確認画面の一例を示す図である。
【図１０】図１０は、吸光度グラフの一例を示す図である。
【図１１】図１１は、図１０に示した吸光度グラフをもとに作成した平均化グラフの一例を示す図である。
【図１２】図１２は、二次微分平均二乗値をグラフ化した二次微分二乗グラフを示す図である。
【図１３】図１３は、平坦波長区間選択画面の一例を示す図である。
【図１４】図１４は、作成条件修正画面の一例を示す図である。
【図１５】図１５は、実施の形態１における画像表示処理の詳細な処理手順を示すフローチャートである。
【図１６】図１６は、実施の形態１における表示画像の観察画面の一例を示す図である。
【図１７】図１７は、実施の形態２における画像処理装置の機能構成を示すブロック図である。
【図１８】図１８は、実施の形態２における画像表示処理の詳細な処理手順を示すフローチャートである。
【図１９】図１９は、実施の形態２における観察画面の一例を示す図である。
【図２０】図２０は、実施の形態３の顕微鏡システムの全体構成を説明する模式図である。
【図２１】図２１は、実施の形態３の顕微鏡システムの機能構成を示すブロック図である。
【図２２】図２２は、特性データのデータ構成例を説明する図である。
【図２３】図２３は、特性データのデータ構成例を説明する他の図である。
【図２４】図２４は、実施の形態３の観察システム制御部が行う処理手順を示すフローチャートである。
【図２５】図２５は、実施の形態３における染色標本属性指定画面の一例を示す図である。
【発明を実施するための形態】
【００２３】
以下、図面を参照し、本発明の好適な実施の形態について詳細に説明する。本実施の形態では、観察対象とする染色標本としてＨ＆Ｅ染色された病理標本（生体組織標本）を例示する。そして、この染色標本を被写体としてマルチバンド撮像し、得られたマルチバンド画像をもとに染色標本の各点（標本点）の色素量を推定する画像処理装置について説明する。なお、この実施の形態によって本発明が限定されるものではない。また、図面の記載において、同一部分には同一の符号を付して示している。
【００２４】
ここで、本実施の形態では、上記したようにＨ＆Ｅ染色された染色標本を観察対象とする。このため、染色標本を染色している色素は、色素Ｈおよび色素Ｅであるが、実際の染色標本内には、これら染色色素の吸収成分の他に、無染色時において吸収成分を持つ例えば赤血球等の組織が存在し得る。すなわち、赤血球は、染色を施さない状態であってもそれ自身が特有の色を有しており、Ｈ＆Ｅ染色後において赤血球自身の色として観察される。そこで、以下では、染色色素を色素Ｈ、色素Ｅおよび赤血球自身の色（以下、「色素Ｒ」と表記する。）の３種類として説明する。
【００２５】
（実施の形態１）
図１は、実施の形態１における画像処理装置１の機能構成を示すブロック図である。実施の形態１の画像処理装置１は、染色標本画像を撮像する染色標本画像撮像部１１と、操作部１２と、表示部１３と、画像処理部１４と、記憶部１６と、装置各部を制御する制御部１７とを備える。ここで、染色標本画像撮像部１１を除く構成は、例えばワークステーションやパソコン等の汎用コンピュータを用いて実現できる。
【００２６】
染色標本画像撮像部１１は、撮像対象の染色標本の観察像をマルチバンド撮像するマルチバンドカメラで構成され、例えば、チューナブルフィルタや二次元ＣＣＤカメラ、チューナブルフィルタを透過する光の波長を調整するフィルタ制御器、二次元ＣＣＤカメラを制御するカメラ制御器等で構成される。実施の形態１では、観察対象の染色標本（以下、「観察染色標本」と呼ぶ。）を撮像対象とする。
【００２７】
ここで、染色標本画像撮像部１１は、染色標本を透過観察可能な光学顕微鏡と接続されている。この光学顕微鏡は、照明光を射出する光源や対物レンズ、染色標本を載置して対物レンズの光軸方向およびこの光軸方向と垂直な面内を移動する電動ステージ、電動ステージ上の染色標本を透過照明するための照明光学系、対物レンズ、染色標本の観察像を結像させるための観察光学系等を備え、照明光学系によって光源からの照明光を染色標本に照射するとともに、対物レンズと協働し、観察光学系によって染色標本の観察像を結像させる。
【００２８】
染色標本画像撮像部１１は、この光学顕微鏡によって観察される染色標本の観察像をチューナブルフィルタを介して二次元ＣＣＤカメラの撮像素子上に投影し、マルチバンド撮像して染色標本画像を得る。チューナブルフィルタは、透過光の波長を電気的に調整可能なフィルタであって、実施の形態１では、１〔ｎｍ〕以上の任意の幅（以下、「選択波長幅」と呼ぶ。）の波長帯域を選択可能なものを用いる。例えば、ケンブリッジリサーチアンドインストルメンテーション社製の液晶チューナブルフィルタ「ＶａｒｉＳｐｅｃ（バリスペック）」等、市販のものを適宜用いることができる。例えば、染色標本画像撮像部１１は、このチューナブルフィルタによって所定の選択波長幅毎に順次波長帯域を選択しながら染色標本の観察像を撮像することで、染色標本画像をマルチバンド画像として得る。
【００２９】
この染色標本画像撮像部１１によって得られる染色標本画像の画素値は、チューナブルフィルタが選択した波長帯域における光の強度に相当し、染色標本の各標本点について選択した波長帯域の画素値が得られる。ここで、染色標本の各標本点とは、投影された撮像素子の各画素に対応する染色標本上の各位置のことであり、以下では、染色標本上の各標本点が染色標本画像の各画素位置に対応しているものとする。
【００３０】
なお、染色標本画像撮像部１１の構成としてチューナブルフィルタを用いた構成を例示したが、これに限定されるものではなく、撮像対象とする染色標本の各標本点における光の強度情報が取得できればよい。例えば、特開平７−１２０３２４号公報に開示されている撮像方式を用い、所定枚数（例えば１６枚）のバンドパスフィルタをフィルタホイールで回転させて切り替えながら、面順次方式で染色標本をマルチバンド撮像する構成としてもよい。
【００３１】
操作部１２は、例えば、キーボードやマウス、タッチパネル、各種スイッチ等によって実現されるものであり、操作入力に応じた操作信号を制御部１７に出力する。表示部１３は、ＬＣＤやＥＬディスプレイ等のフラットパネルディスプレイ、あるいはＣＲＴディスプレイ等の表示装置によって実現されるものであり、制御部１７から入力される表示信号に従って各種画面を表示する。
【００３２】
画像処理部１４は、ＣＰＵ等のハードウェアによって実現される。この画像処理部１４は、分光特性取得手段としてのスペクトル取得部１４１と、作成条件取得手段としての標本作成条件推定処理部１４２と、色素量推定手段としての色素量推定部１４７と、表示画像生成手段としての表示画像生成部１４８とを含む。
【００３３】
スペクトル取得部１４１は、染色標本画像撮像部１１が観察染色標本をマルチバンド撮像して得た染色標本画像（以下、「観察染色標本画像」と呼ぶ。）を構成する各画素位置のスペクトルを取得する。
【００３４】
標本作成条件推定処理部１４２は、観察染色標本の作成条件を推定するための処理を行う。この標本作成条件推定処理部１４２は、解析領域設定手段としての解析領域設定部１４３と、特徴量取得手段としての特徴量取得部１４４と、作成条件推定手段としての作成条件推定部１４５と、色素分光特性決定手段としての基準スペクトル決定部１４６とを備える。
【００３５】
解析領域設定部１４３は、解析領域選択入力依頼部１７１による選択入力依頼に応答して操作部１２から入力されたユーザ操作に従って、観察染色標本画像中の解析領域を設定する。特徴量取得部１４４は、解析領域設定部１４３が設定した解析領域の特徴量を取得する。作成条件推定部１４５は、解析領域の特徴量をもとに、観察染色標本を作成した際の作成条件を推定する。そして、基準スペクトル決定部１４６は、作成条件推定部１４５が推定した作成条件をもとに、基準スペクトル情報１６３に記憶されている基準スペクトルの中から各染色色素（色素Ｈ、色素Ｅおよび色素Ｒ）についてそれぞれ１つの基準スペクトルを選択し、染色色素毎の最適色素分光特性値（以下、最適色素分光特性値を「最適基準スペクトル」と呼ぶ。）を決定する。
【００３６】
色素量推定部１４７は、スペクトル取得部１４１が観察染色標本の各画素位置について取得したスペクトルをもとに、基準スペクトル決定部１４６が決定した色素Ｈ、色素Ｅおよび色素Ｒの最適基準スペクトルを用いて観察染色標本の色素量を推定する。表示画像生成部１４８は、観察染色標本の表示用の画像（表示画像）を生成する。
【００３７】
記憶部１６は、更新記憶可能なフラッシュメモリ等のＲＯＭやＲＡＭといった各種ＩＣメモリ、内蔵或いはデータ通信端子で接続されたハードディスク、ＣＤ−ＲＯＭ等の情報記憶媒体およびその読取装置等によって実現されるものである。この記憶部１６には、画像処理装置１を動作させ、この画像処理装置１が備える種々の機能を実現するためのプログラムや、このプログラムの実行中に使用されるデータ等が一時的または永続的に記憶される。例えば、記憶部１６には、観察染色標本上の各標本位置における色素量を推定するための画像処理プログラム１６１が記憶される。そして、記憶部１６は、色素分光特性記憶手段として、基準スペクトル情報１６３を記憶する。
【００３８】
基準スペクトル情報１６３は、予め取得される各染色色素（色素Ｈ、色素Ｅおよび色素Ｒ）の基準スペクトルのデータを記憶する。この基準スペクトル情報１６３は、異なる複数の染色状態での色素Ｈおよび色素Ｅの基準スペクトルの組み合わせを記憶する。これら異なる複数の染色状態での色素Ｈおよび色素Ｅの染色色素毎の基準スペクトルは、例えば作成条件を異ならせて作成した色素Ｈおよび色素Ｅそれぞれの単一染色標本から取得する。作成条件としては、染色状態を変化させる要因となる例えば観察染色標本の染色に要した染色時間や標本の厚み、染色に使用した薬剤のｐＨ等が挙げられる。また、各基準スペクトルそれぞれについて事前に分散σ2baseが算出され、併せて記憶される。
【００３９】
ここで、色素Ｈおよび色素Ｅの基準スペクトルの取得方法について説明する。例えば、事前に基準スペクトルの組み合わせを取得する複数通りの作成条件（実施の形態１では、染色時間、使用する薬剤のｐＨおよび標本の厚みの組み合わせの３つ）を定義しておく。そして、定義された複数通りの作成条件毎に、その作成条件に従って色素Ｈのみで染色した単一染色標本（以下、「Ｈ単一染色標本」と呼ぶ。）を作成するとともに、その作成条件に従って色素Ｅのみで染色した単一染色標本（以下、「Ｅ単一染色標本」と呼ぶ。）を作成する。
【００４０】
そして、得られた作成条件毎のＨ単一染色標本を例えば染色標本画像撮像部１１を用いる等してマルチバンド撮像する。その後、各Ｈ単一染色標本を順次処理対象としてそのスペクトルを推定する。例えば、後述する図６のステップａ３でスペクトル取得部１４１が行う処理手順と同様の手順で、処理対象のＨ単一染色標本上の複数の標本点における分光透過率ｔ（ｘ，λ）をそのマルチバンド画像をもとに算出し、この分光透過率ｔ（ｘ，λ）を分光吸光度ａ（ｘ，λ）に変換する処理を行う。なお、サンプリングする標本点は、処理対象のＨ単一染色標本上の任意の位置とすることができるが、色素Ｈとして代表的な色分布を持つ画素をサンプリングするのが望ましい。そして、得られた複数の標本点における分光吸光度ａ（ｘ，λ）の平均を算出し、処理したＨ単一染色標本の作成条件についての色素Ｈの基準スペクトルとする。色素Ｅについても同様の処理を行って作成条件毎の基準スペクトルを取得する。そして、作成条件が同じ色素Ｈおよび色素Ｅの基準スペクトルの各組み合わせをそれぞれその作成条件と対応付けて基準スペクトル情報１６３に設定する。
【００４１】
また、基準スペクトル情報１６３は、色素Ｒの基準スペクトルについては、作成条件に関わらず例えば次のようにして取得したものを記憶する。すなわち、染色が施されていない無染色標本を用意し、マルチバンド撮像する。そして、このマルチバンド画像をもとに、無染色標本上の複数の標本点における分光透過率ｔ（ｘ，λ）を算出し、これを分光吸光度ａ（ｘ，λ）に変換する処理を行う。このとき、サンプリングする標本点は、赤血球の領域を選ぶこととする。そして、得られた複数の標本点における分光吸光度ａ（ｘ，λ）の平均を算出し、色素Ｒの基準スペクトルとする。
【００４２】
なお、各染色色素の基準スペクトルの取得方法はこれに限定されるものではなく、分光器等の測定機器を用い、予め定義される各作成条件に該当する染色状態での各染色色素（色素Ｈ、色素Ｅおよび色素Ｒ）のスペクトルを測定することによって取得する構成としてもよい。
【００４３】
また、記憶部１６には、作成条件推定部１４５が観察染色標本を作成した際の作成条件を推定する処理で用いる作成条件判定パラメータＡｄｊの作成条件判定パラメータ分布（不図示）が記憶される。ここで、作成条件判定パラメータＡｄｊは、事前に基準スペクトル情報１６３に記憶されている作成条件が同じ色素Ｈおよび色素Ｅの基準スペクトルの組み合わせ毎に決定され、作成条件判定パラメータ分布は、これら決定された作成条件毎の作成条件判定パラメータＡｄｊをもとに作成される。
【００４４】
後述するように、実施の形態１では、解析部位を核としている。ここで、色素Ｈは、標本内の組織のうち、特に核を染色する。このため、染色標本内の核の領域のスペクトルは、主に色素Ｈの基準スペクトルによって特徴付けられる。そこで、本例では、作成条件が同じ色素Ｈおよび色素Ｅの基準スペクトルの組み合わせを順次決定対象とする。そして、決定対象の組み合わせの色素Ｈおよび色素Ｅの基準スペクトルをグラフ化した基準スペクトルグラフの形状を解析して色素Ｅに対する色素Ｈの基準スペクトル特徴（以下、「Ｈ基準特徴」と呼ぶ。）を抽出し、抽出したＨ基準特徴をもとに、決定対象の組み合わせについての（決定対象の組み合わせに該当する作成条件についての）作成条件判定パラメータＡｄｊを決定する。
【００４５】
図２は、作成条件が同じ１組の色素Ｈおよび色素Ｅの基準スペクトルグラフを示す図である。この基準スペクトルグラフは、横軸を波長とし、縦軸を吸光度値として波長毎の基準スペクトル値をプロットしたものであり、図２では、色素Ｈの基準スペクトルグラフを一点鎖線で示し、色素Ｅの基準スペクトルグラフを二点鎖線で示している。
【００４６】
ここで、作成条件判定パラメータＡｄｊの決定方法について、図２に示した色素Ｈおよび色素Ｅの基準スペクトルの組み合わせに着目し、図３を参照しながら説明する。先ず、図３に示すように、色素Ｅの基準スペクトルグラフから、基準スペクトル値が最大となるピーク波長ＰEを検出する。続いて、検出したピーク波長ＰEよりも長波長側で色素Ｈの基準スペクトルグラフと色素Ｅの基準スペクトルグラフとか交差する波長ＨSを取得する。この波長ＨSは、色素Ｈの波長毎の基準スペクトル値を近似した近似式（以下、「Ｈスペクトル近似式」と呼ぶ。）と、色素Ｅの波長毎の基準スペクトル値を近似した近似式（以下、「Ｅスペクトル近似式」と呼ぶ。）とを求め、ピーク波長ＰEよりも長波長側の波長帯域におけるこれら近似式の交点を求めることで取得する。
【００４７】
図４は、決定対象の色素Ｈおよび色素Ｅの基準スペクトルグラフのピーク波長ＰEよりも長波長側の所定の波長帯域における部分基準スペクトルグラフを示す図である。また、図４中において、Ｈスペクトル近似式のグラフＧ１およびＥスペクトル近似式のグラフＧ３を併せて示している。図４に示すように、色素Ｅの基準スペクトルのピーク波長ＰEよりも長波長側におけるグラフＧ１とグラフＧ３との交点Ｐ１を算出する。そして、算出した交点Ｐ１の波長を、図３に示すように、色素Ｈおよび色素Ｅの基準スペクトルグラフが交差する波長ＨSとして取得する。
【００４８】
続いて、図３に示すように、取得した波長ＨSよりも長波長側で色素Ｈの基準スペクトル値が最大となるピーク波長ＰHを検出する。
【００４９】
また、決定対象の組み合わせを構成する色素Ｈの基準スペクトルの分散σ2base_Hによる誤差を考慮し、その標準偏差σbase_Hをもとに、次式（５）に従ってピーク波長ＰHに対して長波長側に所定幅の分散波長幅Ｗbase_Hを設定する。
【数５】
2011095225.tif 000007
【００５０】
そして、図３に示すように、この分散波長幅Ｗbase_Hと上記したように取得・検出した波長ＨSおよびピーク波長ＰHとをもとに、１つ目のＨ基準特徴ＲWを算出する。このＨ基準特徴ＲWである波長区間は、次式（６）で表される。
ＲW＝（ＰH＋Ｗbase_H）−ＨS  ・・・（６）
【００５１】
また、図３に示すように、波長ＨSにおける基準スペクトル値ａ（ＨS）と、ピーク波長ＰEにおける基準スペクトル値ａ（ＰE）との波長間変化（差分）を２つ目のＨ基準特徴ＲΔPとして算出する。このＨ基準特徴ＲΔPは、次式（７）で表される。
【数６】
2011095225.tif 000008
【００５２】
その後、次式（８）に従い、Ｈ基準特徴ＲWおよびＲΔPの値を用いて決定対象の組み合わせについての作成条件判定パラメータＡｄｊを算出する。ｋは、任意の値が設定される係数である。
【数７】
2011095225.tif 000009
【００５３】
なお、係数ｋは、波長ＨSを算出する際に求めたＨスペクトル近似式およびＥスペクトル近似式をもとに定義する構成としてもよい。例えば、これらの近似式をもとに波長ＨSおよびピーク波長ＰEの波長間の変化率ηを算出し、この変化率ηを用いて次式（９）に従って定義するようにしてもよい。
ｋ＝１−η  ・・・（９）
【００５４】
以上の処理手順で、基準スペクトル情報１６３に記憶されている作成条件毎の色素Ｈおよび色素Ｅの基準スペクトルの各組み合わせそれぞれについて作成条件判定パラメータＡｄｊを決定したならば、作成条件判定パラメータ分布を作成する。図５は、作成条件判定パラメータ分布の一例を示す図である。図５に示すように、作成条件判定パラメータ分布は、作成条件である染色時間、標本の厚みおよび水素イオン濃度指数（ｐＨ）を各軸とした作成条件空間における作成条件判定パラメータＡｄｊの分布を表したものである。
【００５５】
図１に戻り、制御部１７は、ＣＰＵ等のハードウェアによって実現される。この制御部１７は、操作部１２から入力される操作信号や染色標本画像撮像部１１から入力される画像データ、記憶部１６に記憶されるプログラムやデータ等に基づいて画像処理装置１を構成する各部への指示やデータの転送等を行い、画像処理装置１全体の動作を統括的に制御する。
【００５６】
また、制御部１７は、解析領域選択入力依頼部１７１と、作成条件入力依頼部１７２と、色素選択手段としての色素選択入力依頼部１７３と、表示処理手段としての画像表示処理部１７５とを含む。解析領域選択入力依頼部１７１は、解析領域の候補領域（解析領域候補）の選択入力を依頼する処理を行い、操作部１２を介して入力される例えば病理医や臨床検査技師等のユーザ操作に従って解析領域候補を選択する。作成条件入力依頼部１７２は、作成条件推定部１４５によって推定された作成条件の修正入力を依頼する処理を行い、操作部１２を介してユーザ操作を受け付ける。色素選択入力依頼部１７３は、表示対象色素の選択入力を依頼する処理を行い、操作部１２を介してユーザ操作を受け付ける。画像表示処理部１７５は、例えば観察染色標本の表示画像等を表示部１３に表示する処理を行う。
【００５７】
図６は、実施の形態１の画像処理装置１が行う処理手順を示すフローチャートである。なお、ここで説明する処理は、記憶部１６に記憶された画像処理プログラム１６１に従って画像処理装置１の各部が動作することによって実現される。
【００５８】
実施の形態１では、図６に示すように先ず、制御部１７が染色標本画像撮像部１１の動作を制御し、観察染色標本を順次マルチバンド撮像する（ステップａ１）。得られた観察染色標本の染色標本画像の画像データは、記憶部１６に記憶される。
【００５９】
続いて、スペクトル取得部１４１が、観察染色標本画像の各画素位置のスペクトルを取得する（ステップａ３）。例えば、スペクトル取得部１４１は、観察染色標本画像を構成する画素毎に対応する観察染色標本上の標本点におけるスペクトルを推定することによって、各画素位置のスペクトルを取得する。
【００６０】
ここで、詳細なスペクトル推定の手順について説明する。染色標本上の標本点における分光透過率ｔ（ｘ，λ）は、背景技術で次式（１）に示したように、マルチバンド撮像した染色標本画像の位置ベクトルｘで表される任意の画素位置（ｘ）の画素値Ｉ（ｘ，λ）を背景（照明光）のマルチバンド画像の対応する画素位置（ｘ）の画素値Ｉ0（ｘ，λ）で除算することによって得られる。
【数８】
2011095225.tif 000010
【００６１】
実際には、波長λは離散的にしか観測できない。このため、波長方向のサンプル点数をＭとすると、分光透過率ｔ（ｘ，λ）は、次式（１０）に示すようにＭ次元のベクトルとして表される。［］ｔは、転置行列を表す。
【数９】
2011095225.tif 000011
【００６２】
得られた分光透過率ｔ（ｘ，λ）は、次式（１１）に従って分光吸光度ａ（ｘ，λ）に変換できる。以下、分光吸光度を単に「吸光度」と呼ぶ。
ａ（ｘ，λ）＝−ｌｏｇ（ｔ（ｘ，λ））  ・・・（１１）
【００６３】
実施の形態１では、スペクトル取得部１４１は、式（１０）に従って分光透過率ｔ（ｘ，λ）を算出し、式（１１）に従って分光透過率ｔ（ｘ，λ）を吸光度ａ（ｘ，λ）に変換する処理を観察染色標本画像の全ての画素について行い、各画素位置（ｘ）のスペクトルとして吸光度ａ（ｘ，λ）を取得する。取得した観察染色標本画像の各画素位置（ｘ）のスペクトル（吸光度ａ（ｘ，λ））のデータは、取得過程で算出した各画素位置（ｘ）の分光透過率ｔ（ｘ，λ）のデータとともに記憶部１６に記憶される。
【００６４】
その後、図６に示すように、スペクトル取得部１４１は、取得した観察標本画像の各画素位置のスペクトルをもとに、観察染色標本のＲＧＢ画像（以下、「観察染色ＲＧＢ画像」と呼ぶ。）を合成する（ステップａ５）。合成した観察染色ＲＧＢ画像の画像データは、記憶部１６に記憶され、適宜表示部１３に表示処理されてユーザに提示される。
【００６５】
具体的には、スペクトル取得部１４１は、観察染色標本画像の各画素位置についてスペクトルを取得する過程で算出した分光透過率をＲＧＢ値に変換し、観察染色ＲＧＢ画像を合成する。観察染色標本画像上の任意の画素位置（ｘ）における分光透過率をＴ（ｘ）とすると、ＲＧＢ値ＧRGB（ｘ）は、次式（１２）で表される。
ＧRGB（ｘ）＝ＨＴ（ｘ）  ・・・（１２）
【００６６】
ここで、式（１２）のＨは、次式（１３）で定義される行列である。この行列Ｈはシステム行列とも呼ばれており、Ｆはチューナブルフィルタの分光透過率、Ｓはカメラの分光感度特性、Ｅは照明の分光放射特性をそれぞれ表す。
Ｈ＝ＦＳＥ  ・・・（１３）
【００６７】
続いて、図６に示すように、標本作成条件推定処理に移る（ステップａ７）。図７は、標本作成条件推定処理の詳細な処理手順を示すフローチャートである。
【００６８】
図７に示すように、標本作成条件推定処理では先ず、解析領域選択入力依頼部１７１が、ユーザによる解析領域候補の選択操作を受け付けて解析領域候補を選択する（ステップｂ１）。例えば、解析領域選択入力依頼部１７１は、解析領域選択画面を表示部１３に表示し、解析領域候補の選択入力依頼をユーザに通知する処理を行う。そして、解析領域選択入力依頼部１７１は、この選択入力依頼の通知に応答してユーザが操作入力した解析領域候補の選択情報を解析領域設定部１４３に通知する。
【００６９】
図８は、解析領域選択画面の一例を示す図である。図８に示すように、解析領域選択画面は、観察染色画像表示部Ｗ１１を備える。この観察染色画像表示部Ｗ１１には、図６のステップａ５で合成された観察染色ＲＧＢ画像が表示される。また、解析領域選択画面は、解析部位メニューＭ１１と、選択モードメニューＭ１３と、ｍａｎｕａｌ設定メニューＭ１５と、グラフモードメニューＭ１７とを備える。この他、解析領域選択画面には、操作を確定するＯＫボタンＢ１１が配置されている。
【００７０】
解析部位メニューＭ１１は、解析領域候補として選択する領域に映る部位（組織）を指定するためのものであり、「核」「細胞質」「線維」「その他」のいずれか１つを択一的に選択可能なラジオボタンが配置されている。
【００７１】
選択モードメニューＭ１３には、解析領域候補の選択モードとして「ｍａｎｕａｌ」または「ｔｉｓｓｕｅ」を択一的に選択可能なラジオボタンＲＢ１３１，ＲＢ１３３が配置されている。「ｍａｎｕａｌ」は、ユーザ操作に従って解析領域候補を手動で選択する選択モードであり、例えばユーザが観察染色画像表示部Ｗ１１上で指定したシーズ領域を解析領域候補として選択する。「ｔｉｓｓｕｅ」は、ユーザが指定した組織の領域を画像処理によって特定し、解析領域候補として選択する選択モードである。図８の例では、ラジオボタンＲＢ１３５〜ＲＢ１３７によって「核」「細胞質」「線維」のいずれかの組織が指定できるようになっている。
【００７２】
ｍａｎｕａｌ設定メニューＭ１５では、選択モードの１つである「ｍａｎｕａｌ」に関する設定を行う。例えば、ｍａｎｕａｌ設定メニューＭ１５には、ブロックサイズを入力する入力ボックスＩＢ１５１と、ブロック数を入力する入力ボックスＩＢ１５３とが配置されており、それぞれ所望の値を設定することができる。ブロックサイズは、観察染色画像表示部Ｗ１１上で指定するシーズ領域のサイズであり、図８中に例示するように例えば入力ボックスＩＢ１５１に「２」を入力した場合には、１つのシーズ領域は２×２画素のサイズとされる。また、シーズ領域は複数個指定することができるようになっており、その数がブロック数に相当する。ユーザは、所望のブロック数を入力ボックスＩＢ１５３に入力する（図８では「１」）。
【００７３】
また、選択モードが「ｍａｎｕａｌ」の場合、観察染色画像表示部Ｗ１１上でシーズ領域を指定すると、該当する画素位置について取得されているスペクトルがグラフ表示部Ｗ１３においてグラフ表示されるようになっている。グラフモードメニューＭ１７では、このグラフ表示部Ｗ１３におけるグラフ表示に関する設定を行う。例えば、グラフモードメニューＭ１７には、グラフ表示部Ｗ１３において観察染色画像表示部Ｗ１１上で指定したシーズ領域のスペクトル平均グラフを表示させるためのチェックボックスＣＢ１７１が配置されている。チェックボックスＣＢ１７１をチェックしない場合には、シーズ領域を構成する各画素位置のスペクトルグラフがグラフ表示される。例えば、ｍａｎｕａｌ設定メニューＭ１５において入力ボックスＩＢ１５１の入力値が「２」の場合には、シーズ領域は４画素で構成されており、該当する４つの画素位置についてそれぞれ波長毎に取得されているスペクトルをプロットしたスペクトルグラフが表示される。一方、チェックボックスＣＢ１７１をチェックした場合には、これらシーズ領域を構成する各画素位置のスペクトルグラフと併せて、そのスペクトルの波長毎の平均値をプロットしたスペクトル平均グラフが表示される。例えば、図８中のグラフ表示部Ｗ１３では、破線で示す各画素位置のスペクトルグラフと、実線で示すスペクトル平均グラフとが表示されている。
【００７４】
また、グラフモードメニューＭ１７には、グラフ表示部Ｗ１３に表示させるスペクトルグラフの種類として吸光度グラフまたは分光透過率グラフを択一的に選択可能なラジオボタンＲＢ１７１，ＲＢ１７３が配置されている。ラジオボタンＲＢ１７１によって吸光度を選択すれば、シーズ領域の各画素位置について取得されているスペクトルである吸光度がグラフ表示される。一方、ラジオボタンＲＢ１７３によって分光透過率を選択した場合には、吸光度を取得する過程で算出した分光透過率がグラフ表示される。
【００７５】
例えば、選択モードとして「ｍａｎｕａｌ」を選び、解析領域候補を選択する場合の操作手順としては、ユーザは先ず、操作部１２を構成するマウスによって観察染色画像表示部Ｗ１１上の所望の位置をクリックし、シーズ領域を指定する。このとき、観察染色画像表示部Ｗ１１上には、指定した（クリックした）位置にシーズ領域を示すマーカＭＫ１１が表示されるようになっている。また、グラフ表示部Ｗ１３には、シーズ領域を構成する各画素位置のスペクトルがグラフ表示される。一旦指定したシーズ領域は、観察染色画像表示部Ｗ１１上で例えばマーカＭＫ１１をドラッグ＆ドロップ操作する等して移動させることが可能である。これによれば、ユーザは、グラフ表示部Ｗ１３にてスペクトルを確認しながらシーズ領域を指定できる。また、ブロック数として２以上の値を入力している場合には、観察染色画像表示部Ｗ１１上の別の位置を再度クリックすることで、新たなシーズ領域を指定する。操作を確定する場合には、ＯＫボタンＢ１１をクリックする。
【００７６】
以上のようにして操作を確定すると、解析領域選択入力依頼部１７１が、図７のステップｂ１の処理として、指定されたシーズ領域（観察染色画像表示部Ｗ１１上のマーカＭＫ１１の画素位置）を解析領域候補として選択する。このとき、解析領域選択入力依頼部１７１は、基準候補領域の選択情報を画像処理部１４の解析領域設定部１４３に通知する。ここで、解析領域候補の選択情報は、解析領域候補の画素位置を含む。
【００７７】
なお、図８の解析領域選択画面において、選択モードメニューＭ１３にて選択モードとして「ｔｉｓｓｕｅ」を選択した場合には、内部処理として、ユーザがラジオボタンＲＢ１３５〜ＲＢ１３７のうちのいずれかを選択することによって指定した組織が映る画素位置を教師データを用いて抽出する処理が行われ、解析領域候補が選択される。教師データは、予め各組織の代表的な分光特性パターンや色情報を測定して作成しておく。そして、観察染色標本画像あるいは観察染色ＲＧＢ画像を画像処理することによって、ユーザが指定した組織の領域を抽出する。
【００７８】
解析領域候補が選択されると、図７に示すように、続いて解析領域設定部１４３が、解析領域を設定する（ステップｂ３）。例えば、解析領域設定部１４３は、解析領域選択入力依頼部１７１から通知された解析領域候補の選択情報をもとに、観察染色ＲＧＢ画像中で解析領域候補と画素値が類似する画素を探索する処理を行い、解析領域を設定する。
【００７９】
具体的な処理手順としては先ず、解析領域設定部１４３は、観察染色ＲＧＢ画像の各画素値をＲＢ色空間に写像する。このとき、解析領域候補が複数の画素位置で構成されている場合には、例えば解析領域候補を構成する各画素の写像点の平均値（解析領域候補の各画素のＲＢ色空間における座標値の平均値）を算出し、解析領域候補代表点とする。
【００８０】
続いて、解析領域設定部１４３は、解析領域候補以外の各画素を順次処理対象とし、解析領域候補の写像点（または解析領域候補代表点）と処理対象の画素の写像点との距離Ｄｉｓｔを算出する。そして、得られた距離Ｄｉｓｔを処理対象の画素についての類似度として得る。ここで、解析領域候補の写像点（または解析領域候補代表点）をＳ（Ｒ，Ｂ）とし、処理対象の画素（ｘi，ｙi）（ｉ＝１，２，・・・，ｎ）の写像点（ＲＧ色空間における座標値）をｓ（ｒi，ｂi）とすると、ｓ（ｒi，ｂi）とＳ（Ｒ，Ｂ）との距離Ｄｉｓｔは、次式（１４）で表される。ｎは、処理対象とする解析領域候補以外の画素の数を表す。
【数１０】
2011095225.tif 000012
【００８１】
そして、解析領域設定部１４３は、得られた解析領域候補以外の各画素の類似度を閾値処理し、類似度が高い画素（例えば類似度が所定の閾値以上の画素）を抽出する。そして、解析領域設定部１４３は、ここで抽出した画素および解析領域候補の各画素で構成される領域を解析領域とする。閾値処理で用いる閾値ＳＴｈは、観察染色ＲＧＢ画像中の解析領域候補の画素値をもとに設定する。例えば、解析領域候補が複数の画素位置で構成されている場合には、各画素の写像点（ＲＢ色空間における座標値）の分散Ｖ（Ｓ）を求め、次式（１５）に従って設定する。ｋは係数であり、任意に設定できる。
【数１１】
2011095225.tif 000013
【００８２】
なお、類似度の算出方法は、上記した手法に限定されるものではなく、適宜選択して用いることができる。例えば、輝度値の類似度や色分布の類似度、あるいはスペクトルの類似度等を算出する。このとき１つの類似度を算出することとしてもよいし、これら複数の類似度を組み合わせて総合的な類似度を算出してもよい。
【００８３】
以上のようにして設定された解析領域は、表示部１３に確認表示されてユーザに提示されるようになっている。例えば解析領域選択入力依頼部１７１が、解析領域確認画面を表示部１３に表示する処理を行う。またこのとき、解析領域選択入力依頼部１７１は、解析領域の修正入力依頼をユーザに通知する処理を行う。
【００８４】
図９は、解析領域確認画面の一例を示す図である。図９に示すように、解析領域確認画面は、観察染色画像表示部Ｗ２１と、解析領域表示部Ｗ２３とを備える。観察染色画像表示部Ｗ２１には、観察染色ＲＧＢ画像が表示される。また、解析領域表示部Ｗ２３には、観察染色ＲＧＢ画像中の解析領域を識別表示した解析領域識別画像が表示される。例えば、解析領域外の画素値を所定の色（例えば白）に置き換えて解析領域外の画素を非表示とした画像が表示される。
【００８５】
この解析領域確認画面は、修正モードメニューＭ２１を備え、設定された解析領域に過不足があり不十分と判断した場合に修正できるようになっている。修正モードメニューＭ２１には、解析領域の修正モードとして「追加」または「削除」を択一的に選択可能なラジオボタンＲＢ２１１，ＲＢ２１３が配置されている。また、この他、解析領域確認画面には、操作を確定する確定ボタンＢ２１が配置されている。
【００８６】
例えばユーザは、解析領域表示部Ｗ２３の解析領域識別画像によって解析領域の設定漏れが生じていると判断した場合には、ラジオボタンＲＢ２１１を選択した上で、例えば観察染色画像表示部Ｗ２１上の解析領域として追加したい画素（追加画素）の位置をクリックする。一方、解析領域表示部Ｗ２３の解析領域識別画像によって解析領域が過剰に設定されていると判断した場合には、ラジオボタンＲＢ２１３を選択した上で、例えば観察染色画像表示部Ｗ２１上の解析領域から除外したい画素（除外画素）の位置をクリックする。操作を確定する場合には、確定ボタンＢ２１をクリックする。
【００８７】
そして、以上のようにして解析領域の修正入力依頼の通知に応答して修正操作が入力されると（図７のステップｂ５：Ｙｅｓ）、解析領域選択入力依頼部１７１が、修正情報を解析領域設定部１４３に通知するようになっている。ここで、修正情報は、追加画素として指定された画素位置や除外画素として指定された画素位置等を含む。解析領域設定部１４３は、この修正情報に従って解析領域を修正する（ステップｂ７）。
【００８８】
例えば、解析領域設定部１４３は、修正情報に追加画素が設定されている場合には、観察染色ＲＧＢ画像の画素値をもとに、追加画素と類似する画素であって、追加画素と連結する画素を抽出する。例えば、追加画素の隣接画素から順番にその輝度値を閾値処理する。閾値は、例えば追加画素の輝度値をもとに設定する。そして、追加画素と輝度値が類似する画素を、追加画素と連結している限り抽出する。解析領域設定部１４３は、このようにして抽出した画素を解析領域に追加する。一方、解析領域設定部１４３は、修正情報に除外画素が設定されている場合には、観察染色ＲＧＢ画像の画素値をもとに、除外画素と類似する画素であって、除外画素と連結する画素を抽出する。例えば、除外画素の隣接画素から順番にその輝度値を閾値処理する。閾値は、例えば除外画素の輝度値をもとに設定する。そして、除外画素と輝度値が類似する画素を、除外画素と連結している限り抽出する。解析領域設定部１４３は、このようにして抽出した画素を解析領域から除外する。
【００８９】
最終的に、解析領域設定部１４３は、以上のようにして設定・修正した解析領域を構成する各画素の画素位置に、解析部位に応じた解析領域ラベルを付与する。すなわち、図８の解析部位メニューＭ１１で指定されている部位が「核」の場合には、解析領域の画素位置に解析領域ラベルＯｂｊ＿Ｎを付与する。同様に、指定されている部位が「細胞質」であれば解析領域ラベルＯｂｊ＿Ｃを付与し、指定されている部位が「線維」であれば解析領域ラベルＯｂｊ＿Ｆを付与し、指定されている部位が「その他」であれば解析領域ラベルＯｂｊ＿Ｏを付与する。
【００９０】
ここで、解析領域は、観察染色標本内に含まれる例えば核や細胞質、線維の領域といった代表的な組織の領域であるのが望ましい。すなわち、図８の解析領域選択画面でユーザが指定するシーズ領域が、この代表的な組織の領域であるのが望ましい。以下、実施の形態１では、ユーザが、図８の観察染色画像表示部Ｗ１１上で核が映る領域をシーズ領域として指定し、解析部位メニューＭ１１にてラジオボタンＲＢ１１１によって核を指定することとして説明する。
【００９１】
続いて、特徴量取得部１４４が、設定された解析領域（観察染色標本画像中の解析領域ラベルが付与された画素位置）の特徴量を取得する。例えば、特徴量取得部１４４は、解析領域を構成する各画素についてスペクトル取得部１４１が取得した波長毎のスペクトル（吸光度）をグラフ化し、グラフ化した吸光度グラフの形状を解析して特徴量を取得する。
【００９２】
先ず、特徴量取得部１４４は、図７に示すように、解析領域内の画素について取得されているスペクトルをもとに、吸光度グラフを作成する（ステップｂ９）。ここで、解析領域内の画素をｉ（ｉ＝１，２，３，・・・，ｎ）とし、スペクトルの波長数をＤとする。そして、画素ｉのλdにおけるスペクトルをａi（λd）（ｄ＝１，２，３，・・・，Ｄ）とすると、吸光度ベクトルＡ（λ）は、次式（１６）で表される。
【数１２】
2011095225.tif 000014
【００９３】
また、平均吸光度ベクトルＡ￣（λ）は、次式（１７）で表される。なお、Ａ￣は、Ａの上に平均値を表す記号「￣」が付いていることを示す。
【数１３】
2011095225.tif 000015
【００９４】
図７のステップｂ９において、特徴量取得部１４４は、例えばこの平均吸光度ベクトルＡ￣（λ）をグラフ化し、吸光度グラフを作成する。またこのとき、特徴量取得部１４４は、その分散σ2Obj_Nを算出する。図１０は、吸光度グラフの一例を示す図である。図１０に示すように、特徴量取得部１４４は、横軸を波長（波長数）とし、縦軸を吸光度値として、算出した平均吸光度ベクトルＡ￣（λ）の値を波長毎にプロットして吸光度グラフを得る。
【００９５】
続いて、図７に示すように、特徴量取得部１４４は、吸光度グラフの形状解析を行って特徴量を取得する（ステップｂ１１）。ここで、上記したように、実施の形態１では、色素Ｈによって染色される核を解析部位としており、この核の領域のスペクトルは、主に色素Ｈの基準スペクトルによって特徴付けられる。具体的には、解析領域のスペクトルは、色素Ｈの基準スペクトル（例えば図１０に例示）の影響を受けて長波長側で吸光度値が減少する特性を有する。また、このときの吸光度値の減少量は、観察染色標本の染色状態と相関すると想定できる。
【００９６】
そこで、実施の形態１では先ず、特徴量取得部１４４は、解析領域のスペクトルをもとに作成した吸光度グラフをもとに、長波長帯域において、波長間の吸光度変化が小さく、図１０中に例示するように吸光度グラフが平坦な波長区間（以下、「平坦波長区間」と呼ぶ。）Ｚを特徴量の１つとして算出する。また、特徴量取得部１４４は、この平坦波長区間Ｚとは別の特徴量として、ピーク変化度ΔＰを算出する。例えば、特徴量取得部１４４は、ピーク波長Ｐにおける吸光度値と、平坦波長区間Ｚ内の吸光度値の平均値との差分を、ピーク変化度ΔＰとする。
【００９７】
具体的な平坦波長区間Ｚの算出手順としては先ず、吸光度グラフから、吸光度値が最大となる波長をピーク波長Ｐとして検出する。なお、解析部位として想定されている核や細胞質、線維等の組織毎に、予めピーク波長Ｐが現れる波長帯域を学習しておくようにしてもよい。そして、指定された解析部位について学習されている波長帯域の吸光度値を参照してピーク波長Ｐを検出するようにしてもよい。このようにすれば、ピーク波長Ｐを検出するために吸光度グラフ全域を参照する必要がない。
【００９８】
続いて、検出したピーク波長Ｐから最長波長である波長Ｄまでの波長帯域（Ｐ−Ｄ波長帯域）において、連続する２波長間の吸光度を平均化した平均化グラフを作成する。図１１は、図１０に示した吸光度グラフをもとに作成した平均化グラフの一例を示す図である。図１１では、実線で示す平均化グラフとともに、Ｐ−Ｄ波長帯域における吸光度グラフを破線で示している。
【００９９】
続いて、波長変化を強調するため、吸光度の二次微分Ｌ（ｉ）を次式（１８）に従って算出するとともに、二次微分Ｌ（ｉ）の波長間平均Ｌ￣（ｉ）（以下、「二次微分平均」と呼ぶ。）を次式（１９）に従って算出する。なお、Ｌ￣は、Ｌの上に記号「￣」が付いていることを示す。
【数１４】
2011095225.tif 000016
【０１００】
さらに、変化方向の影響を除外するため、二次微分平均Ｌ￣（ｉ）の二乗値を算出する。ここで、平均吸光度ベクトルＡ￣（λ）は、その分散σ2Obj_Nによる誤差を含んでいる。そこで、この誤差を考慮し、平坦波長区間Ｚを算出する際に対象とする対象波長を制限する。例えば、平均吸光度ベクトルＡ￣（λ）について算出されている分散σ2Obj_Nの標準偏差σObj_Nをもとに、次式（２０）に従って分散波長幅ＷObj_Nを算出する。ｓｔｅｐは、観察染色標本画像の波長間隔（染色標本画像撮像部１１が観察染色標本画像を撮像した際のチューナブルフィルタの選択波長幅）を表す。
【数１５】
2011095225.tif 000017
【０１０１】
そして、ピーク波長Ｐに対して算出した分散波長幅ＷObj_Nを考慮した波長Ｐ＋ＷObj_Nを対象波長とし、二次微分平均Ｌ￣（ｉ）の二乗値Ｌ￣（ｉ）2（ｉ＝Ｐ＋ＷObj_N，Ｐ＋ＷObj_N＋２，Ｐ＋ＷObj_N＋３，・・・，Ｐ＋ＷObj_N＋Ｄ）を算出する。
【０１０２】
図１２に、二次微分平均二乗値Ｌ￣（ｉ）2をグラフ化した二次微分二乗グラフを示す。特徴量取得部１４４は、この二次微分二乗グラフをもとに、Ｐ＋ＷObj_Nよりも長波長側の二次微分平均二乗値Ｌ￣（ｉ）2の値を順次参照し、最初に所定の閾値以下となった波長をＳとして検出する。また、その次に所定の閾値以下となった波長をＥとして検出する。例えば、図１２中に一点鎖線で示す閾値を用い、波長Ｓおよび波長Ｅを検出する。そして、検出した波長Ｓおよび波長Ｅをもとに、波長Ｓから波長Ｅ−１までの波長区間を平坦波長区間Ｚとする。
【０１０３】
なお、平坦波長区間Ｚの算出方法はこれに限定されるものではない。例えば、ユーザ操作に従って選択する構成としてもよい。図１３は、平坦波長区間選択画面の一例を示す図である。図１３に示すように、平坦波長区間選択画面は、観察染色ＲＧＢ画像または解析領域画像が表示される画像表示部Ｗ３１を備える。いずれの画像を表示させるのかは、画像表示メニューＭ３１で選択できるようになっている。すなわち、この画像表示メニューＭ３１で「ＨＥ」を選択すると、画像表示部Ｗ３１には観察染色ＲＧＢ画像が表示される。一方、「解析領域」を選択すると、画像表示部Ｗ３１には、観察染色ＲＧＢ画像中の解析領域を識別表示した解析領域画像が表示される。例えば、解析領域外の画素値を所定の色（例えば白）に置き換えて解析領域外の各画素を非表示とした画像が表示される。
【０１０４】
また、平坦波長区間選択画面は、グラフ表示部Ｗ３３を備え、図７のステップｂ９で作成された吸光度グラフが表示される。グラフ表示部Ｗ３３に表示させる吸光度グラフの表示波長は、表示波長メニューＭ３３で適宜選択できる。例えば、ユーザは、この表示波長メニューＭ３３において、表示させたい波長範囲の短波長側の値を入力ボックスＩＢ３３１に入力し、長波長側の値を入力ボックスＩＢ３３３に入力する。なお、吸光度グラフに限らず、吸光度グラフをもとに作成した平均化グラフ（図１１を参照）をグラフ表示部Ｗ３３に表示する構成としてもよい。
【０１０５】
そして、ユーザは、グラフ表示部Ｗ３３に表示される吸光度グラフを見ながら、区間選択メニューＭ３５の入力ボックスＩＢ３５１に所望の波長Ｓを入力し、入力ボックスＩＢ３５３に所望の波長Ｅ−１を入力する。これにより、入力ボックスＩＢ３５１の波長Ｓから入力ボックスＩＢ３５３の波長Ｅ−１までの波長区間が平坦波長区間Ｚとして選択される。また、このようにして入力ボックスＩＢ３５１，ＩＢ３５３に値を入力すると、グラフ表示部Ｗ３３において、現在選択されている平坦波長区間Ｚが例えば図１３に例示するように破線によって示され、吸光度グラフ上で識別表示されるようになっている。
【０１０６】
また、具体的なピーク変化度ΔＰの算出手順としては先ず、平坦波長区間Ｚにおける吸光度値平均ａ￣flatを算出する。なお、ａ￣は、ａの上に記号「￣」が付いていることを示す。この吸光度値平均ａ￣flatは、平坦波長区間Ｚ内の各波長の吸光度値をａ（λi）（ｉ＝Ｓ，Ｓ＋１，・・・，Ｅ）とすると、次式（２１）で表される。ｓｔｅｐは、観察染色標本画像の波長間隔を表す。
【数１６】
2011095225.tif 000018
【０１０７】
そして、算出した平坦波長区間Ｚにおける吸光度値平均ａ￣flatをもとに、次式（２２）に従ってピーク変化度ΔＰを算出する。ａ（Ｐ）は、ピーク波長Ｐにおける吸光度値である。
【数１７】
2011095225.tif 000019
【０１０８】
以上のようにして特徴量取得部１４４が算出した平坦波長区間Ｚおよびピーク変化度ΔＰの各特徴量は、記憶部１６に記憶される。
【０１０９】
そして、解析領域の特徴量として平坦波長区間Ｚおよびピーク変化度ΔＰを算出したならば、続いて、作成条件推定部１４５が、これら平坦波長区間Ｚおよびピーク変化度ΔＰをもとに、観察染色標本の作成条件を推定する。実施の形態１では、作成条件として、観察染色標本の染色に要した染色時間、染色に使用した薬剤のｐＨおよび標本の厚みを推定する。ここで行う作成条件の推定では、上記したように事前に決定されて記憶部１６に記憶される作成条件判定パラメータＡｄｊの作成条件判定パラメータ分布を用いる。
【０１１０】
上記のように、実施の形態１では解析部位が核であり、この核の領域のスペクトルは主に色素Ｈの基準スペクトルによって特徴付けられる。このため、特徴量取得部１４４が特徴量として取得した平坦波長区間Ｚおよびピーク変化度ΔＰは、観察染色標本の染色状態に応じたＨ基準特徴ＲWおよびＲΔPと相関すると想定できる。ここで、式（８）に示して上記したように、予め決定される作成条件判定パラメータＡｄｊは、Ｈ基準特徴ＲΔPをＨ基準特徴ＲWで除することによって得ている。そこで、図７に示すように、作成条件推定部１４５は先ず、作成条件パラメータ算出手段として、次式（２３）に従って解析領域の作成条件判定パラメータＡｄｊAを算出する（ステップｂ１３）。
【数１８】
2011095225.tif 000020
【０１１１】
そして、作成条件推定部１４５は、算出した解析領域の作成条件判定パラメータＡｄｊAに該当する作成条件判定パラメータＡｄｊを作成条件判定パラメータ分布から選出することによって観察染色標本の作成条件を推定する（ステップｂ１５）。解析領域の作成条件判定パラメータＡｄｊAと一致する作成条件判定パラメータＡｄｊがあれば、その作成条件判定パラメータＡｄｊを選出し、その作成条件を観察染色標本の作成条件として推定する。一致する作成条件判定パラメータＡｄｊがなければ、解析領域の作成条件判定パラメータＡｄｊAと値が類似する作成条件判定パラメータＡｄｊを選出し、その作成条件を観察染色標本の作成条件として推定する。具体的には、解析領域の作成条件判定パラメータＡｄｊAとの差の絶対値が最小となる作成条件判定パラメータＡｄｊを類似する作成条件判定パラメータＡｄｊとして選出する。ここで、差の絶対値が最小となる作成条件判定パラメータＡｄｊが複数存在する場合もあるが、このような場合には、これら複数の作成条件判定パラメータＡｄｊをもとに観察染色標本の作成条件をそれぞれ推定する（観察染色標本の作成条件を複数通り推定する）。
【０１１２】
なお、上記したステップｂ１５では、解析領域の作成条件判定パラメータＡｄｊAとの差の絶対値が最小となる作成条件判定パラメータＡｄｊを類似する作成条件判定パラメータＡｄｊとして選出することとした。これに対し、予め設定される所定の閾値を用いて差の絶対値を閾値処理し、例えば差の絶対値が閾値よりも小さい作成条件判定パラメータＡｄｊを類似する作成条件判定パラメータＡｄｊとして選出する構成としてもよい。そして、このときに作成条件判定パラメータＡｄｊを複数選出した場合には、各作成条件判定パラメータＡｄｊをもとに観察染色標本の作成条件を複数通り推定するようにしてもよい。
【０１１３】
また、異なる作成条件に対して同一の作成条件判定パラメータＡｄｊの値が定義され得るため、解析領域の作成条件判定パラメータＡｄｊAと一致する作成条件判定パラメータＡｄｊが複数存在する場合もあり得る。このような場合にも、これら複数の作成条件判定パラメータＡｄｊをもとに観察染色標本の作成条件をそれぞれ推定する。
【０１１４】
以上のようにして推定された観察染色標本の作成条件は、表示部１３に確認表示されてユーザに提示されるようになっている。例えば、作成条件入力依頼部１７２が、作成条件修正画面を表示部１３に表示する処理を行う。またこのとき、作成条件入力依頼部１７２が、作成条件の修正入力依頼をユーザに通知する処理を行う。
【０１１５】
図１４は、作成条件修正画面の一例を示す図である。図１４に示すように、作成条件修正画面は、観察染色画像表示部Ｗ４１を備える。この観察染色画像表示部Ｗ４１には、観察染色ＲＧＢ画像が表示される。また、作成条件修正画面は、作成条件表示タブＴＭ４１（ＴＭ４１−１，２，３，・・・）を備える。上記したステップｂ１５において複数の作成条件判定パラメータＡｄｊが選出されて観察染色標本の作成条件が複数通り推定された場合には、この作成条件表示タブＴＭ４１によって各作成条件を選択・表示できるようにしている。この作成条件表示タブＴＭ４１には、推定された染色時間、標本の厚みおよびｐＨの各値が修正可能に表示される。例えば、入力ボックスＩＢ４１１，ＩＢ４１３の数値を変更することによって、推定された染色時間を修正することができる。同様に、入力ボックスＩＢ４１５の数値を変更することによって推定された標本の厚みを修正でき、入力ボックスＩＢ４１７の数値を変更することによって推定されたｐＨを修正できる。この他、作成条件修正画面には、作成条件表示タブＴＭ４１での操作を確定して修正内容を登録する登録ボタンＢ４１と、修正入力を終えるＯＫボタンＢ４３とが配置されている。
【０１１６】
そして、以上のようにして作成条件の修正入力依頼の通知に応答して修正操作が入力されると（図７のステップｂ１７：Ｙｅｓ）、作成条件入力依頼部１７２が、修正情報を作成条件推定部１４５に通知するようになっている。ここで、修正情報は、修正された染色時間、標本の厚みおよびｐＨの各値を含む。作成条件推定部１４５は、この修正情報に従って推定した作成条件を修正する（ステップｂ１９）。以上のようにして推定・修正された観察染色標本の作成条件は、記憶部１６に記憶される。
【０１１７】
ここで、作成条件修正画面において登録ボタンＢ４１をクリックすると、選択している作成条件表示タブＴＭ４１の各入力ボックスＩＢ４１１，ＩＢ４１３に入力されている作成条件の各値と、解析領域の作成条件判定パラメータＡｄｊAとの組み合わせが新たな作成条件判定パラメータＡｄｊとして作成条件判定パラメータ分布に追加されるようになっている。例えば、基準スペクトル決定部１４６が行う処理として後述する最適基準スペクトルの決定手順に従って各染色色素の最適基準スペクトルを決定し、決定した最適基準スペクトルを記憶部１６に記憶するとともに、作成条件の各値と解析領域の作成条件判定パラメータＡｄｊAとの組み合わせを決定した最適基準スペクトルと関連付けた上で記憶部１６に記憶されている作成条件判定パラメータ分布に追加して更新する。
【０１１８】
続いて、基準スペクトル決定部１４６が、作成条件推定部１４５が推定・修正した観察染色標本の作成条件に従って基準スペクトル情報１６３から該当する色素Ｈ、色素Ｅおよび色素Ｒの基準スペクトルを選択し（ステップｂ２１）、観察染色標本の色素量推定に用いる最適基準スペクトルとして決定する（ステップｂ２３）。
【０１１９】
すなわち、基準スペクトル決定部１４６は、観察染色標本の作成条件が作成条件推定部１４５による推定値である場合（ユーザ操作に従って修正されていない場合）には、この作成条件と対応付けられた色素Ｈおよび色素Ｅの基準スペクトルと、色素Ｒの基準スペクトルとを基準スペクトル情報１６３から読み出して選択し、最適基準スペクトルとして決定する。
【０１２０】
一方、基準スペクトル決定部１４６は、観察染色標本の作成条件がユーザ操作によって修正されている場合には、修正後の作成条件に従って基準スペクトル情報１６３に記憶されている基準スペクトルを選択した後補正し、最適基準スペクトルを決定する。
【０１２１】
ここで、この場合の最適基準スペクトルの決定手順について説明する。先ず、基準スペクトル決定部１４６は、修正後の作成条件である染色時間をｔ、標本の厚みをｄ、水素イオン指数（ｐＨ）をｐとをもとに、作成条件判定パラメータ分布上で距離が最小となる作成条件判定パラメータＡｄｊを作成条件判定パラメータ分布から選出する。そして、基準スペクトル決定部１４６は、選出した作成条件判定パラメータＡｄｊに対する染色時間をｔ0、標本の厚みをｄ0、水素イオン指数（ｐＨ）をｐ0として取得する。
【０１２２】
続いて、基準スペクトル決定部１４６は、染色時間ｔ0とｔ、標本の厚みｄ0とｄ、水素イオン指数（ｐＨ）ｐ0とｐをそれぞれ比較し、差分が最小となる作成条件を抽出する。そして、基準スペクトル決定部１４６は、抽出した作成条件に従って補正マトリクスを作成する。例えば、標本の厚みｄ0とｄの差分が最小の場合には、次式（２４）に従って補正マトリクスを作成する。
【数１９】
2011095225.tif 000021
【０１２３】
すなわち、式（２５）に示す変換マトリクスは、標本の厚みｄ0を固定として染色時間を変化させた場合の基準スペクトルの変動を定義したものであり、例えば、予め測定しておいた標本の厚みｄ0に対する染色時間の変化に応じた各染色色素（色素Ｈ、色素Ｅおよび色素Ｒ）の基準スペクトルの変動量を近似した関数で実現できる。詳細には、複数通りの標本の厚みに対する変換マトリクスを用意しておき、標本の厚みｄ0の値に応じた変換マトリクスを選択して用いる。同様に、式（２６）に示す変換マトリクスは、標本の厚みｄ0を固定としてｐＨを変化させた場合の基準スペクトルの変動を定義したものであり、例えば、予め測定しておいた標本の厚みｄ0に対するｐＨの変化に応じた各染色色素（色素Ｈ、色素Ｅおよび色素Ｒ）の基準スペクトルの変動量を近似した関数で実現できる。
【０１２４】
なお、染色時間を固定として標本の厚みおよびｐＨをそれぞれ変化させた場合の基準スペクトルの変動を定義した変換マトリクスや、ｐＨを固定として染色時間および標本の厚みをそれぞれ変化させた場合の基準スペクトルの変動を定義した変換マトリクスについても同様に用意しておき、差分が最小となる作成条件に応じて適宜用いる。なお、変換マトリクスは、例示したものに限定されるものではない。例えば、標本の厚みを固定とした場合のこれら染色時間の変化あるいはｐＨの変化に伴う各染色色素の基準スペクトルの変動をモデル化することによって定義しておくこととしてもよい。同様に、染色時間を固定とした場合の標本の厚みの変化あるいはｐＨの変化に伴う各染色色素の基準スペクトルの変動や、ｐＨを固定とした場合の染色時間の変化あるいは標本の厚みの変化に伴う各染色色素の基準スペクトルの変動をモデル化して変換マトリクスを定義しておくこととしてもよい。
【０１２５】
そして、この補正マトリクスを用い、作成条件として染色時間ｔ0、標本の厚みｄ0および水素イオン指数（ｐＨ）ｐ0が対応付けられた各染色色素の基準スペクトルを補正したならば、基準スペクトル決定部１４６は、補正後の基準スペクトルを最適基準スペクトルとする。そして、図６のステップａ７にリターンし、その後ステップａ９に移行する。
【０１２６】
そして、ステップａ９では、色素量推定部１４７が、観察染色標本画像の各画素位置について取得されているスペクトル（吸光度ａ（ｘ，λ））をもとに、ステップａ７の標本作成条件推定処理で各染色色素についてそれぞれ決定した最適基準スペクトルを用いて観察染色標本の色素量を推定する。
【０１２７】
背景技術で式（２）に示して説明したように、分光透過率ｔ（ｘ，λ）に関しては、ランベルト・ベールの法則が成り立つ。また、分光透過率ｔ（ｘ，λ）は、上記した式（１１）を用いて吸光度ａ（ｘ，λ）に変換できる。実施の形態１においてもこれらの式を適用して色素量を推定する。すなわち、ランベルト・ベールの法則によれば、観察染色標本画像の各画素（ｘ，ｙ）に対応する観察染色標本上の各標本点における吸光度ａ（ｘ，λ）は、次式（２７）で表される。ただしこのとき、色素Ｈの基準スペクトルｋHとして、色素Ｈについて決定した最適基準スペクトルを用い、色素Ｅの基準スペクトルｋEとして、色素Ｅについて決定した最適基準スペクトルを用い、色素Ｒの基準スペクトルｋRとして、色素Ｒについて決定した最適基準スペクトルを用いる。
【数２０】
2011095225.tif 000022
【０１２８】
各画素（ｘ，ｙ）に対応する観察染色標本上の各標本点における色素Ｈ、色素Ｅおよび色素Ｒの色素量は、例えば背景技術で式（３）に示して説明した手法を適用し、重回帰分析を行うことで推定（算出）することができる。推定した色素量のデータは、記憶部１６に記憶される。
【０１２９】
以上のように色素量を推定したならば、図６に示すように、画像表示処理に移る（ステップａ１１）。この画像表示処理では、ステップａ９で推定した色素量をもとに観察染色標本の表示用の画像（表示画像）を生成し、表示部１３に表示する処理を行う。図１５は、実施の形態１における画像表示処理の詳細な処理手順を示すフローチャートである。
【０１３０】
画像表示処理では先ず、画像表示処理部１７５が、図６のステップａ５で合成された観察染色標本のＲＧＢ画像（観察染色ＲＧＢ画像）を表示画像として表示部１３に表示する処理を行う（ステップｃ１）。
【０１３１】
続いて、色素選択入力依頼部１７３が、表示対象色素の選択入力依頼の通知を表示部１３に表示する処理を行い、この選択入力依頼の通知に応答して表示対象色素の選択操作が入力されず、表示対象色素の選択がない間は（ステップｃ３：Ｎｏ）、ステップｃ９に移行する。
【０１３２】
そして、ユーザによる表示対象色素の選択操作が入力された場合には（ステップｃ３：Ｙｅｓ）、表示画像生成部１４８が、観察染色ＲＧＢ画像をもとに、表示対象色素によって染色されている観察染色標本上の領域（表示対象色素を含む画素位置）を識別可能に表した表示画像を生成する（ステップｃ５）。例えば、表示画像生成部１４８は、図６のステップａ９で観察染色標本画像の画素毎に推定された観察染色標本上の各標本点における色素量をもとに、表示対象色素を含む（表示対象色素の色素量が“０”ではない）画素位置を選出し、表示対象色素染色領域とする。例えば、表示対象色素として色素Ｈが選択された場合には、色素Ｈの色素量が“０”でなく、色素Ｈを含む画素位置を選出して表示対象色素染色領域を設定する。そして、表示画像生成部１４８は、観察染色ＲＧＢ画像をもとに、表示対象色素染色領域内の画素をその他の画素と識別可能に表した表示画像を生成する。
【０１３３】
続いて画像表示処理部１７５が、ステップｃ５で生成された表示画像を表示部１３に表示する処理を行い（ステップｃ７）、その後ステップｃ９に移行する。このとき、ステップｃ７で生成された表示画像を、既に表示処理されている表示画像にかえて表示処理することとしてもよいし、これらを並べて表示処理する構成としてもよい。
【０１３４】
そして、ステップｃ９では、画像表示の終了判定を行い、終了しない間は（ステップｃ９：Ｎｏ）、ステップｃ３に戻って表示対象色素の選択操作を受け付ける。例えばユーザによって画像表示の終了操作が入力された場合には終了すると判定する（ステップｃ９：Ｙｅｓ）。そして、図６のステップａ１１にリターンし、その後処理を終える。
【０１３５】
ここで、ユーザが表示画像を観察する際の操作例について説明する。図１６は、実施の形態１における表示画像の観察画面の一例を示す図である。図１６に示す観察画面は、２つの画像表示部Ｗ５１，Ｗ５３を備える。また、観察画面は、表示対象色素を選択するための色素選択メニューＭ５１を備え、各染色色素を表示対象色素として個別に選択できるようになっている。図１６では、チェックボックスＣＢ５１によって色素Ｈが表示対象色素として選択されている。
【０１３６】
そして、図１６中に向かって左側の画像表示部Ｗ５１には、例えば観察染色ＲＧＢ画像が表示される。一方、向かって右側の画像表示部Ｗ５３には、表示対象色素染色領域を識別表示した表示対象色素識別画像が表示される。この表示対象色素識別画像は、図１５のステップｃ５の処理で生成される表示画像の一例であり、例えば、表示対象色素染色領域を表示し、それ以外の領域を非表示とした画像として生成される。図１６では、観察染色ＲＧＢ画像中の色素Ｈについて設定した表示対象色素染色領域を表示し、それ以外の画素を非表示とした表示対象色素識別画像が表示されている。この場合の内部処理としては例えば、表示画像生成部１４８は、色素選択メニューＭ５１で選択されている表示対象色素をもとに表示対象色素領域を設定し、設定した表示対象色素染色領域外の画素値を所定の色（例えば白）に置き換えて表示対象色素識別画像を生成する。
【０１３７】
さらに、観察画面は、画像表示部Ｗ５３に表示される表示対象色素識別画像の描画モードを指定する描画メニューＭ５３を備える。例えば、図１６の例では、描画モードとして「なし」「輪郭」「色」「パターン」のいずれか１つを択一的に選択可能なラジオボタンが配置されている。図１６に例示するように、描画メニューＭ５３にて「なし」を選択すると、表示対象色素識別画像がそのまま画像表示部Ｗ５３に表示される。「輪郭」を選択すると、表示対象色素識別画像中で、表示対象色素毎の表示対象色素染色領域が輪郭線で囲まれて識別表示される。「色」を選択すると、表示対象色素識別画像中で、表示対象色素毎の表示対象色素染色領域が所定の描画色に置き換えられて識別表示される。描画色は、表示対象色素毎に予め設定しておく。「パターン」を選択すると、表示対象色素識別画像中で、表示対象色素毎の表示対象色素染色領域が所定の塗り潰しパターンで識別表示される。塗り潰しパターンは、表示対象色素毎に予め設定しておく。例えば、色素選択メニューＭ５１で表示対象色素を２つ以上選択した場合に、この描画メニューＭ５３で「色」や「パターン」を選択すれば、各表示対象色素それぞれの表示対象色素染色領域を識別可能に表示させることができる。また、描画メニューＭ５３には、ユーザ設定ボタンＢ５３が配置されており、クリックすることで表示対象色素に割り当てる色や塗り潰しパターン、あるいは描画メニューＭ５３で提示する識別表示項目の編集等が行える。
【０１３８】
なお、識別表示の方法は、上記したものに限定されるものではない。例えば、表示対象色素の色素量を含む各画素位置について、その色合いを色素量の値に応じて段階的に変化させるといったことも可能である。
【０１３９】
以上説明したように、実施の形態１によれば、観察染色標本の作成条件を推定することができる。そして、推定した作成条件に応じた基準スペクトルを基準スペクトル情報１６３に記憶されている各染色色素の基準スペクトルの組み合わせから選択し、各染色色素の最適基準スペクトルを決定することができる。そして、観察染色標本画像の各画素について取得したスペクトルをもとに、決定した各染色色素の最適基準スペクトルを用いて観察染色標本上の標本点における染色色素の色素量を推定することができる。したがって、観察染色標本の作成条件に応じた最適な各染色色素の基準スペクトルを用いて、観察対象の染色標本の色素量を精度良く推定することができる。また、ユーザは、観察染色標本の作成条件を記録しておく必要がなく、ユーザの手間を省くことができる。
【０１４０】
また、推定した各染色色素の色素量をもとにユーザが選択した表示対象色素を含む観察染色標本画像の画素位置を選出し、表示対象色素によって染色されている観察染色標本上の領域（表示対象色素を含む画素位置）を識別可能に表した表示画像を生成することができる。したがって、観察染色標本内を視認性良く表した画像をユーザに提示することができる。これによれば、ユーザによる観察効率を向上させることができる。ユーザにとっては、観察したい所望の染色色素を選択することによって、観察染色標本中の所望の染色色素の領域を個別にあるいは組み合わせて視認性良く観察できる。
【０１４１】
なお、上記した実施の形態１では、作成条件推定部１４５が観察染色標本の作成条件を推定することとした。これに対し、作成条件の推定を行わない構成としてもよい。例えば、図６のステップａ７で行う標本作成条件推定処理にかえて、観察染色標本を作成した際の作成条件をユーザ操作に従って取得する処理を行う構成としてもよい。
【０１４２】
この場合には、作成条件入力依頼部１７２が、作成条件入力依頼手段として、例えば図１４に例示した作成条件修正画面と同様の作成条件入力画面を表示部１３に表示処理し、観察染色標本を作成した際の作成条件の入力依頼をユーザに通知する処理を行う。そして、作成条件入力依頼部１７２は、この入力依頼の通知に応答してユーザが操作入力した作成条件を観察染色標本の作成条件として取得する。続いて、基準スペクトル決定部１４６が、作成条件入力依頼部１７２による処理によってユーザ操作に従って取得した観察染色標本の作成条件をもとに、作成条件判定パラメータ分布を参照して該当する作成条件判定パラメータＡｄｊを選出する。そして、基準スペクトル決定部１４６は、選出した作成条件判定パラメータＡｄｊの基準スペクトルを基準スペクトル情報１６３から読みだして色素量推定に用いる各染色色素の最適基準スペクトルを決定する。このとき、基準スペクトル決定部１４６は、取得した作成条件と一致する作成条件があれば、対応する各染色色素の基準スペクトルを最適基準スペクトルとして決定する。一方、一致する作成条件がない場合には、基準スペクトル決定部１４６は、図７のステップｂ２１，ｂ２３と同様の処理手順で各染色色素の基準スペクトルを選択し、補正して各染色色素の最適基準スペクトルを決定する。
【０１４３】
（実施の形態２）
図１７は、実施の形態２における画像処理装置１ｂの機能構成を示すブロック図である。なお、図１７において、実施の形態１で説明した画像処理装置１と同様の構成については同一の符号を付して示している。
【０１４４】
図１７に示すように、実施の形態２の画像処理装置１ｂは、染色標本画像撮像部１１と、操作部１２と、表示部１３と、画像処理部１４ｂと、記憶部１６ｂと、制御部１７ｂとを備える。
【０１４５】
画像処理部１４ｂは、スペクトル取得部１４１と、標本作成条件推定処理部１４２と、色素量推定部１４７と、色素量補正手段としての色素量補正部１４９ｂと、分光特性合成手段としてのスペクトル合成部１５０ｂと、表示画像生成手段としての表示画像生成部１４８ｂとを含む。色素量補正部１４９ｂは、色素量調整入力依頼部１７７ｂによる調整入力依頼に応答して操作部１２から入力されたユーザ操作に従って、色素量推定部１４７が推定した色素Ｈ、色素Ｅおよび色素Ｒの色素量を補正する。スペクトル合成部１５０ｂは、色素量補正部１４９ｂによる補正後の色素Ｈ、色素Ｅおよび色素Ｒの色素量をもとに分光透過率ｔ（ｘ，λ）を合成する。
【０１４６】
また、記憶部１６ｂには、観察染色標本上の各標本位置における色素量を推定し補正するための画像処理プログラム１６１ｂと、基準スペクトル情報１６３とが記憶される。
【０１４７】
そして、制御部１７ｂは、解析領域選択入力依頼部１７１と、作成条件入力依頼部１７２と、色素指定手段としての色素選択入力依頼部１７３と、色素量調整入力依頼部１７７ｂと、表示処理手段としての画像表示処理部１７５ｂとを含む。色素量調整入力依頼部１７７ｂは、色素量の調整入力を依頼する処理を行い、操作部１２を介してユーザによる色素量の調整操作を受け付ける。
【０１４８】
この実施の形態２の画像処理装置１ｂは、実施の形態１において図６に示した処理手順において、ステップａ１１の画像表示処理にかえて図１８に示す画像表示処理を行う。なお、画像処理装置１ｂが行う処理は、記憶部１６ｂに記憶された画像処理プログラム１６１ｂに従って画像処理装置１ｂの各部が動作することによって実現される。
【０１４９】
画像表示処理では先ず、図１８に示すように、画像表示処理部１７５ｂが、実施の形態１で説明した図６のステップａ５で合成されている観察染色ＲＧＢ画像を表示画像として表示部１３に表示する処理を行う（ステップｄ１）。
【０１５０】
続いて、色素選択入力依頼部１７３が、表示対象色素の選択入力依頼の通知を表示部１３に表示する処理を行う。そして、この選択入力依頼の通知に応答して表示対象色素の選択操作が入力された場合には（ステップｄ３：Ｙｅｓ）、色素量補正部１４９ｂが、表示対象色素以外の色素を非表示として色素量を補正する（ステップｄ５）。例えば、実施の形態１で説明した図６のステップａ９で観察染色標本画像の画素毎に推定された色素量のうち、表示対象色素以外の色素の色素量を例えば全て“０”に置き換えて補正する。
【０１５１】
続いて、スペクトル合成部１５０ｂが、補正後の色素Ｈ、色素Ｅおよび色素Ｒの色素量をもとに分光透過率ｔ（ｘ，λ）を合成する（ステップｄ７）。例えば、スペクトル合成部１５０ｂは、次式（２８）に従い、実施の形態１で示した図７の標本作成条件推定処理で各染色色素について決定した最適基準スペクトルを用いて画素位置（ｘ）における分光透過率ｔ（ｘ，λ）を新たに合成する。
【数２１】
2011095225.tif 000023
【０１５２】
続いて、表示画像生成部１４８ｂが、新たに合成した各画素位置（ｘ）の分光透過率ｔ（ｘ，λ）をＲＧＢ値に変換し、ＲＧＢ画像を合成することによって表示画像を生成する（ステップｄ９）。分光透過率ｔ（ｘ，λ）をＲＧＢ値に変換する処理は、図６のステップａ５の手順と同様に、実施の形態１で上記した式（１２），（１３）を用いて行う。ここで合成されるＲＧＢ画像は、表示対象色素のみの染色状態を表した（表示対象色素の色素量のみを視覚化した）画像である。
【０１５３】
そして、画像表示処理部１７５ｂが、ステップｄ９で生成された表示画像を表示部１３に表示する処理を行い（ステップｄ１１）、その後ステップｄ２３に移行する。このとき、ステップｄ９で生成された表示画像を、既に表示されている表示画像にかえて表示処理することとしてもよいし、これらを並べて表示処理する構成としてもよい。
【０１５４】
また、以上のようにして表示画像を表示している間、表示対象色素についての色素量の調整操作を受け付けるようになっている。例えば、色素量調整入力依頼部１７７ｂは、色素量の調整入力依頼の通知を表示部１３に表示する処理を行い、この調整入力依頼の通知に応答して色素量の調整操作が入力されると（ステップｄ１３：Ｙｅｓ）、入力された調整量を色素量補正部１４９ｂに通知する。
【０１５５】
そして、色素量補正部１４９ｂが、表示対象色素の色素量を、色素量調整入力依頼部１７７ｂから通知された調整量に従って補正する（ステップｄ１５）。その後、スペクトル合成部１５０ｂが、補正後の色素Ｈ、色素Ｅおよび色素Ｒの色素量をもとに、ステップｄ７と同様の処理手順で上記した式（３０）に従って新たに分光透過率ｔ（ｘ，λ）を合成する（ステップｄ１７）。そして、表示画像生成部１４８ｂが、新たに合成した各画素位置の分光透過率ｔ（ｘ，λ）を、ステップｄ９と同様の手順で上記した式（１２），（１３）を用いてＲＧＢ値に変換し、ＲＧＢ画像を合成することによって表示画像を生成する（ステップｄ１９）。
【０１５６】
続いて画像表示処理部１７５ｂが、ステップｄ１９で生成された表示画像を表示部１３に表示する処理を行い（ステップｄ２１）、その後ステップｄ２３に移行する。このとき、ステップｄ１９で生成された表示画像を、既に表示されている表示画像にかえて表示処理することとしてもよいし、これらを並べて表示処理する構成としてもよい。
【０１５７】
そして、ステップｄ２３では、画像表示の終了判定を行い、終了しない間は（ステップｄ２３：Ｎｏ）、ステップｄ３に戻って表示対象色素の選択操作を受け付ける。例えばユーザによって画像表示の終了操作が入力された場合には終了すると判定する（ステップｄ２３：Ｙｅｓ）。
【０１５８】
ここで、ユーザが表示画像を観察する際の操作例について説明する。図１９は、実施の形態２における表示画像の観察画面の一例を示す図である。図１９に示す観察画面は、３つの画像表示部Ｗ６１，Ｗ６３，Ｗ６５を備える。また、観察画面は、表示対象色素を選択するための色素選択メニューＭ６１や、画像表示部Ｗ６３，Ｗ６５に表示される表示対象色素染色画像の描画モードを指定する描画メニューＭ６３を備える。
【０１５９】
そして、図１９中に向かって左側の画像表示部Ｗ６１には、例えば観察染色ＲＧＢ画像が表示される。一方、向かって中央および右側の画像表示部Ｗ６３，Ｗ６５には、それぞれ表示対象色素の色素量のみを表した表示対象色素染色画像が表示される。この表示対象色素染色画像は、図１８のステップｄ９，ｄ１９の処理で生成される表示画像に相当し、図１９では、表示対象色素として選択されている色素Ｈのみの染色状態を表した画像が表示される。
【０１６０】
さらに、観察画面は、色素量調整メニューＭ６５を備え、表示対象色素の色素量を調整するためのスライダーバーＳＢ６５や、スライダーバーＳＢ６５での操作を確定するＯＫボタンＢ６５等が配置されている。例えば、画像表示部Ｗ６３や画像表示部Ｗ６５に表示されている表示対象色素染色画像を観察・診断中に、表示対象色素の染色状態を濃くしたい、あるいは薄くしたいといった場合に、ユーザは、この色素量調整メニューＭ６５においてスライダーバーＳＢ６５を操作し、表示対象色素の色素量の調整量を入力する。図１９中、右側の画像表示部Ｗ６５に表示されている表示対象色素染色画像は、スライダーバーＳＢ６５を用い、中央の画像表示部Ｗ６３に表示されている表示対象色素染色画像よりも色素Ｈの色素量を薄く調整することによって表示させた画像である。
【０１６１】
以上説明したように、実施の形態２によれば、実施の形態１と同様の効果を奏するとともに、ユーザの調整操作に従って推定した表示対象色素の色素量を補正することができる。そして、補正後の各染色色素の色素量をもとに各画素位置のスペクトルを合成し、ＲＧＢ画像を合成することによって表示画像を生成することができる。あるいは、表示対象色素以外の色素の色素量をゼロに補正することで、表示対象色素のみの色素量を視覚化した画像を表示画像として生成することができる。したがって、各染色色素による染色状態を調整して観察染色標本内を視認性良く表した画像をユーザに提示することができる。ユーザにとっては、所望の染色色素を選択して色素量を調整したり、観察・診断に不要な染色色素を除く等して、各染色色素による染色状態を個別に調整して視認性良く観察できる。したがって、診断精度を向上させることが可能となる。
【０１６２】
（実施の形態３）
実施の形態３は、本発明を顕微鏡システムに適用したものである。図２０は、実施の形態３の顕微鏡システム２の全体構成を説明する模式図である。また、図２１は、顕微鏡システム２の機能構成を示すブロック図である。実施の形態３の顕微鏡システム２は、画像取得手段としての観察部３と、観察システム制御部５と、分光特性記憶手段としての特性データ記憶部７とを備え、これらの間がデータの授受可能に接続されて構成されている。
【０１６３】
観察部３は、図２１に示すように、観察染色標本を観察するための染色標本観察部３１と、観察染色標本画像を撮像するための染色標本画像撮像部３３とを備える。ここで、染色標本画像撮像部３３は、実施の形態１で説明した染色標本画像撮像部１１と同様に構成される。
【０１６４】
染色標本観察部３１は、観察染色標本を透過観察可能な顕微鏡で構成され、照明光を射出する光源や対物レンズ、観察染色標本を載置して対物レンズの光軸方向およびこの光軸方向と垂直な面内を移動する電動ステージ、電動ステージ上の観察染色標本を透過照明するための照明光学系、対物レンズ、観察染色標本の観察像を結像させるための観察光学系等を備える。この染色標本観察部３１は、照明光学系によって光源からの照明光を観察染色標本に照射するとともに、対物レンズと協働し、観察光学系によって観察染色標本の観察像を結像させる。
【０１６５】
また、観察システム制御部５は、観察部３の染色標本画像撮像部３３によって撮像された観察染色標本画像を医師等が観察・診断するためのものであり、ワークステーションやパソコン等の汎用コンピュータで実現される。この観察システム制御部５は、観察部３を構成する染色標本観察部３１および染色標本画像撮像部３３に対する動作指示を行い、染色標本画像撮像部３３から入力される観察染色標本画像を処理してディスプレイに表示する。
【０１６６】
具体的には、観察システム制御部５は、操作部５２と、表示部５３と、画像処理部５４と、記憶部５６と、制御部５７とを備える。ここで、実施の形態３の顕微鏡システム２は、実施の形態１の画像処理装置１の構成をもとにしたものであり、画像処理部５４は、スペクトル取得部５４１と、標本作成条件推定処理部５４２と、色素量推定部５４７と、表示画像生成部５４８とを備え、標本作成条件推定処理部５４２は、解析領域設定部５４３と、特徴量取得部５４４と、作成条件推定部５４５と、基準スペクトル決定部５４６とを含む。
【０１６７】
また、制御部５７は、解析領域選択入力依頼部５７１と、作成条件入力依頼部５７２と、色素選択入力依頼部５７３と、画像表示処理部５７５と、染色標本属性入力依頼部５７６と、特性データ選択部５７７と、システム環境設定部５７８とを含む。
【０１６８】
なお、画像処理部５４の各部５４１〜５４８および制御部５７を構成する解析領域選択入力依頼部５７１、作成条件入力依頼部５７２、色素選択入力依頼部５７３および画像表示処理部５７５は、実施の形態１で説明した同一名称の機能部と同様の処理を行う。なお、ここでは、実施の形態１を適用する場合を例示したが、実施の形態１の変形例や実施の形態２の構成をもとにして画像処理部５４や制御部５７を構成することとしてもよい。
【０１６９】
この制御部５７において、染色標本属性入力依頼部５７６は、観察染色標本の属性を示す属性値をユーザ操作に従って指定する。ここでは、観察染色標本の属性（染色標本属性）が染色種類、臓器、注目対象組織および施設の４つの属性項目で構成されることとし、染色標本属性入力依頼部５７６は、観察染色標本に関するこれら４つの属性項目の属性値をユーザ操作に従って指定する。また、実施の形態３では、ユーザは、この観察染色標本の染色標本属性の指定と併せて、その観察染色標本を観察する際の顕微鏡（染色標本観察部３１）の観察倍率を指定するようになっている。
【０１７０】
また、特性データ選択部５７７は、染色標本属性入力依頼部５７６によって指定された染色標本属性をもとに、特性データ記憶部７に記憶されている特性データの中から１つ以上の特性データを選択する。
【０１７１】
システム環境設定部５７８は、観察部３の動作環境（システム環境）を設定するためのシステムパラメータを設定する。例えば、システム環境設定部５７８は、染色標本観察部３１の動作環境を設定するための観察パラメータや、染色標本画像撮像部３３の動作環境を設定するための撮像パラメータをシステムパラメータとして設定する。
【０１７２】
記憶部５６には、観察システム制御部５７を動作させ、この観察システム制御部５７が備える種々の機能を実現するためのプログラムや、このプログラムの実行中に使用されるデータ等の他、観察染色標本上の各標本位置における色素量を推定するための画像処理プログラム５６１が記憶される。
【０１７３】
特性データ記憶部７は、染色標本属性の各属性項目の属性値に応じた特性データを記憶する。この特性データ記憶部７は、例えばネットワークを介して観察システム制御部５と接続されるデータベース装置で実現され、観察システム制御部５とは離れた別の場所に設置されて特性データを記憶し、特性データを管理する。なお、この特性データは、観察システム制御部５７の記憶部５６に記憶しておく構成としてもよい。
【０１７４】
図２２および図２３は、特性データ記憶部７に記憶される特性データのデータ構成例を説明する図である。ここで、図２２は、特性データ記憶部７において、属性項目の１つである染色種類と関連付けられて記憶される特性データの一覧を示している。また、図２３は、特性データ記憶部７において、属性項目の１つである施設と関連付けられて記憶される特性データの一覧を示している。図２２および図２３に示した属性項目毎の特性データの関連付けは、例えば公知のデータベース管理ツールを用いて管理される。ただし、特性データのデータ構成はこれに限定されるものではなく、各属性項目の属性値を指定することによってその属性値に応じた特性データを取得可能な構成であればよい。
【０１７５】
具体的には、特性データ記憶部７には、図２２に示すように、染色種類についての特性データとして、その属性値毎に、属性項目である染色色素および施設と、観察パラメータである観察倍率と、作成条件である染色時間、標本の厚みおよびｐＨと、測定日付と、分光特性値とが関連付けられて記憶される。この染色種類と関連付けられる分光特性値（データセットＡ−０１〜０３，Ａ−１１〜１３，Ａ−２１，Ａ−３１等）は、対応する染色種類の染色色素について予め測定した分光特性値（スペクトルデータ）であって、対応する施設（その分光特性値の測定対象とした染色標本（単一染色標本）を採取した医療施設）で、対応する観察倍率，染色時間，標本の厚み，ｐＨを条件として、対応する測定日に測定した分光特性値を記憶する。ここで、分光特性値は、実施の形態１等で説明した基準スペクトルに相当するものであり、例えば分光吸光度として設定される。ただし、分光透過率や分光反射率等のスペクトル特徴値を用いることとしてもよい。
【０１７６】
また、特性データ記憶部７には、図２３に示すように、施設についての特性データとして、その属性値毎に、属性項目である染色種類および臓器と、観察パラメータである観察倍率と、測定日付と、白色画像信号値、照明分光特性値およびカメラ分光特性値の各システム分光特性値とが関連付けられて記憶される。この施設と関連付けられるシステム分光特性値は、対応する施設で測定されたシステムに関する分光特性値を記憶する。詳細には、白色画像信号値（データセットＣ−０１〜Ｃ−０５，・・・）は、組織が存在していない領域を撮像した画像信号値であって、対応する観察倍率で、対応する測定日に測定された画像信号値を記憶する。また、照明光分光特性値（データセットＣ−１１〜Ｃ−１５，・・・）は、対応する観察倍率で、対応する測定日に分光計等を用いて測定された染色標本観察部３１の照明の分光特性値である。また、カメラ分光特性値（データセットＣ−２１〜Ｃ−２５，・・・）は、対応する観察倍率で、対応する測定日に測定された染色標本画像撮像部３３のカメラ（二次元ＣＣＤカメラ）の分光特性値である。
【０１７７】
なお、線維領域や位相物体等について特殊撮影を行う場合に対応し、例えばＮＡ値やデフォーカス量、光量といった観察パラメータが異なる条件で分光特性値を測定し、その値を測定時の条件と対応付けて特性データ記憶部７に記憶しておくこととしてもよい。また、施設によっては、染色過程や染色に用いる試薬等が異なる場合がある。このような場合に対応し、施設毎に対応する条件で分光特性値を測定し、その値を測定時の条件と対応付けて特性データ記憶部７に記憶しておくこととしてもよい。
【０１７８】
図２４は、実施の形態３の観察システム制御部５が行う処理手順を示すフローチャートである。なお、図２４において、実施の形態１と同様の処理工程については、同一の符号を付して示している。
【０１７９】
図２４に示すように、先ず、染色標本属性入力依頼部５７６が、染色標本属性指定画面を表示部５３に表示して染色標本属性の指定依頼を通知する処理を行い、操作部５２を介してユーザによる染色標本属性および観察倍率の指定操作を受け付ける（ステップｅ１１）。
【０１８０】
図２５は、実施の形態３における染色標本属性指定画面の一例を示す図である。図２５に示す染色標本属性指定画面には、染色種類、臓器、注目対象組織および施設の４つの属性項目の属性値を指定するためのスピンボックスＢ８１〜Ｂ８４、観察倍率を指定するためのスピンボックスＢ８５、これらの各スピンボックスでの操作を確定するＯＫボタンＢ８１、操作をキャンセルするキャンセルボタンＢ８２等が配置されている。各スピンボックスＢ８１〜Ｂ８４は、該当する属性項目について指定可能な属性値の一覧を選択肢として提示し、指定操作を受け付けるものあり、特性データ記憶部７に記憶されている各属性項目の属性値が選択肢として提示される。例えば、染色種類を指定するためのスピンボックスＢ８１では、図２２または図２３に例示したＨＥ染色、ＶＢ染色、ＤＡＢ染色の他、例えばオルセイン染色、ＭＴ染色等の染色種類の属性値と、該当する属性値が含まれない場合に指定する「その他」の各選択肢が提示される。臓器、注目対象組織および施設の各スピンボックスＢ８２〜Ｂ８４も同様に、特性データ記憶部７に記憶されている属性値や、「その他」の各選択肢が提示される。また、スピンボックスＢ８５は、観察倍率の値として指定可能な値を選択肢として提示し、選択を促すものであり、特性データ記憶部７に記憶されている観察倍率の値が選択肢として提示される。
【０１８１】
この染色標本属性指定画面においてユーザは、染色種類として、観察染色標本に施されている染色の種類を指定する。また、ユーザは、臓器として、その観察染色標本が採取された臓器を指定する。また、ユーザは、その観察染色標本を観察・診断する際に注目する注目対象の組織を指定する。また、ユーザは、施設として、例えばその観察染色標本を採取した医療施設を指定する。また、ユーザは、これら４つの染色標本属性とともに、観察倍率を指定する。また、染色標本属性指定画面には、メモ欄Ｍ８１が配置されており、例えばその観察染色標本を作成した日付や、その観察染色標本を観察・診断した日付といった事項をユーザが自由に記入できるようになっている。
【０１８２】
ここで、図２５の染色標本属性指定画面においてユーザが、染色種類を「ＨＥ染色」、臓器を「腎臓」、注目対象組織を「弾性線維」、施設を「Ａ病院」、観察倍率を「２０倍」として指定した場合を例に挙げて図２４のステップｅ１２以降の処理を説明する。
【０１８３】
染色標本属性の属性値が指定されたならば、続いて図２４に示すように、特性データ選択部５７７が、特性データ記憶部７を参照し、指定された染色標本属性の属性値に応じた１つ以上の特性データを選択する（ステップｅ１２）。具体的には、例示した条件の場合、図２２に示した染色種類についての特性データの中から、染色種類が「ＨＥ染色」であり施設が「Ａ病院」であり観察倍率が「２０倍」であるレコードＲ７１〜Ｒ７４を選択し、対応する分光特性値のデータセットＡ−０１，Ａ−０２，Ａ−１１，Ａ−１２を取得する。また、図２３に示した施設についての特性データの中から、施設が「Ａ病院」であり染色種類が「ＨＥ染色」であり臓器が「腎臓」であり観察倍率が「２０倍」であるレコードＲ７５を選択し、対応するシステム分光特性値のデータセットＣ−０１，Ｃ−１１，Ｃ−２１を取得する。
【０１８４】
続いてシステム環境設定部５７８が、ステップｅ１２で選択された特性データをもとに、システムパラメータ（観察パラメータおよび撮像パラメータ）を設定する（ステップｅ１３）。ここで、撮像パラメータは、マルチバンドカメラの動作に関する値であり、システム環境設定部５７８は、設定した撮像パラメータの値を染色標本画像撮像部３３に通知して染色標本画像撮像部３３に対する動作指示を行う。染色標本画像撮像部３３は、このシステム環境設定部５７８による動作指示に応答し、通知された撮像パラメータに従って、例えばゲインの設定、露光時間の設定、チューナブルフィルタによって選択する波長帯域（選択波長幅）の設定等を行ってマルチバンドカメラを駆動する。
【０１８５】
実施の形態３では、システム環境設定部５７８は、この撮像パラメータの１つとして、チューナブルフィルタの選択波長帯域（選択波長幅）を設定する。例えば、図３を参照して説明した波長ＨSの付近（例えば前後５〔ｎｍ〕）の波長帯域における選択波長幅を、チューナブルフィルタによって選択可能な最小の波長幅である１〔ｎｍ〕に設定する。具体的には、取得した分光特性値のデータセットＡ−０１，Ａ−０２，Ａ−１１，Ａ−１２をもとに、作成条件が同じデータセットＡ−０１，Ａ−１１の分光特性値を組み合わせて波長ＨSを求める。同様に、作成条件が同じデータセットＡ−０２，Ａ−１２を組み合わせて波長ＨSを求める。そして、これら組み合わせ毎に求めた波長ＨSの付近の波長帯域における選択波長幅を、例えば１〔ｎｍ〕に設定する。一方、波長ＨSの前後５〔ｎｍ〕の波長帯域以外の波長帯域における選択波長幅を、初期値（例えば５〔ｎｍ〕）に設定する。染色標本画像撮像部３３は、ここで設定される波長帯域毎の選択波長幅に従い、チューナブルフィルタによって選択する波長帯域を順次選択し、選択した波長帯域毎に観察染色標本画像を撮像する。
【０１８６】
また、システム環境設定部５７８は、２つ目の撮像パラメータとして露光時間を設定する。例えば、システム環境設定部５７８は、ステップｅ１２で選択された白色画像信号値のデータセット（例示した条件ではデータセットＣ−０１）を用い、白色画像信号値の最大値が所定輝度値を持つように露光時間を調整し、波長ＨSの前後５〔ｎｍ〕の波長帯域以外の波長帯域における露光時間とする。一方、波長ＨSの前後５〔ｎｍ〕の波長帯域における露光時間については、システム環境設定部５７８は先ず、観察部３の染色標本観察部３１や染色標本画像撮像部３３に対する各種動作指示を行い、指定された観察倍率で白色画像信号値を取得する。そして、システム環境設定部５７８は、取得した白色画像信号値を用い、測定波長毎に露光時間を算出する。これによれば、波長ＨS付近について、観察時（観察染色標本画像の撮像時）の環境に応じて露光時間を設定できる。
【０１８７】
なお、ここでは、チューナブルフィルタの選択波長帯域（選択波長幅）および露光時間の２つを撮像パラメータとして設定することとしたが、これ以外の設定に関する値についても、必要に応じて適宜撮像パラメータとして設定することができる。
【０１８８】
一方、観察パラメータは、顕微鏡の動作に関する値であり、システム環境設定部５４４は、設定した観察パラメータの値を染色標本観察部３１に通知して染色標本観察部３１に対する動作指示を行う。染色標本観察部３１は、このシステム環境設定部５４４による動作指示に応答し、通知された観察パラメータ等に従って例えば対物レンズの観察倍率の切り換えや切り換えた観察倍率等に応じた光源の調光制御、各種光学素子の切り換え、電動ステージの移動等、観察染色標本の観察に伴う顕微鏡各部の調整を行う。
【０１８９】
実施の形態３では、システム環境設定部５４４は、この観察パラメータとして、ステップｅ１１の通知に応答して指定された観察倍率の値を設定する。なお、設定する観察パラメータは、観察倍率に限らず、例えば焦点位置や絞り等の値を適宜観察パラメータとして設定することができる。
【０１９０】
その後、システム環境設定部５７８は、順次選択波長幅とともに該当する波長帯域における露光時間を染色標本画像撮像部３３に出力するとともに、観察パラメータとして設定した観察倍率の値を染色標本観察部３１に出力する。これによれば、観察対象標本に応じて最適な観察部３の動作環境（システム環境）を自動的に設定することができる。そしてこの結果、観察部３は、システム環境設定部５７８が設定したシステムパラメータに従って動作し、選択波長幅毎に観察染色標本をマルチバンド撮像して観察染色標本画像を取得する（ステップｅ１４）。
【０１９１】
続いて、画像処理部５４のスペクトル取得部５４１が、観察染色標本画像の各画素位置のスペクトルを取得する（ステップａ３）。具体的には、スペクトル取得部５４１は、実施の形態１と同様にして観察染色標本画像を構成する画素毎に対応する観察染色標本上の標本点におけるスペクトルを推定し、各画素位置のスペクトルを取得する。
【０１９２】
そして、スペクトル取得部５４１は、取得した観察標本画像の各画素位置のスペクトルをもとに、観察染色ＲＧＢ画像を合成する（ステップａ５）。このステップａ５の処理は、実施の形態１で説明した図６のステップａ５と同様に実現できる。ただし、実施の形態３では、スペクトル取得部５４１は、上記した式（１２）を構成するシステム行列Ｈを算出する際に、ステップｅ１２で選択された特性データを用いる。具体的には、システム行列Ｈを表す上記した式（１３）において、カメラの分光感度特性Ｓとして選択したカメラ分光特性値のデータセットＣ−２１を用いる。また、照明の分光放射特性Ｅとして選択した照明光分光特性値のデータセットＣ−１１を用いる。これによれば、指定された施設に応じた適切なカメラの分光感度特性Ｓおよび照明の分光放射特性Ｅの各値を用いて観察染色ＲＧＢ画像を合成することができる。
【０１９３】
その後、標本作成条件推定処理部５４２が、ステップｅ１２で選択された特性データをもとに、標本作成条件推定処理で用いる作成条件判定パラメータ分布を作成する（ステップｅ６）。具体的には、例示した条件の場合、分光特性値のデータセットＡ−０１，Ａ−０２，Ａ−１１，Ａ−１２によって定まる作成条件判定パラメータＡｄｊをもとに作成条件判定パラメータ分布を作成する。
【０１９４】
ここで、作成条件判定パラメータＡｄｊは、実施の形態１で説明したように、作成条件が同じ色素Ｈおよび色素Ｅの基準スペクトルの組み合わせ毎に決定されるものである。例示の条件では、図２２に示すように、データセットＡ−０１とデータセットＡ−１１の作成条件が同じである。また、データセットＡ−０２とデータセットＡ−１２の作成条件が同じである。したがって、これらを組み合わせ、データセットＡ−０１とデータセットＡ−１１の各分光特性値によって定まる作成条件判定パラメータＡｄｊと、データセットＡ−０２とデータセットＡ−１２の各分光特性値によって定まる作成条件判定パラメータＡｄｊとをもとに作成条件判定パラメータ分布を作成する。
【０１９５】
なお、作成条件判定パラメータＡｄｊの決定方法および作成条件判定パラメータ分布の作成方法は、実施の形態１と同様である。また、作成条件判定パラメータＡｄｊは、事前に決定して特性データ記憶部７あるいは記憶部５６に記憶しておく構成としてよい。この場合には、事前に特性データ記憶部７においてデータセットとして記憶されている同一条件の分光特性値の組み合わせ毎に作成条件判定パラメータＡｄｊを決定しておき、ステップｅ６では、選択された特性データに応じた作成条件判定パラメータＡｄｊを作成条件空間に分布させることで作成条件判定パラメータ分布を作成すればよい。また、作成条件判定パラメータ分布についても事前に作成しておく構成としてよい。すなわち、前述のように事前に決定した全ての作成条件判定パラメータＡｄｊを作成条件空間に分布させたものを作成条件判定パラメータ分布として作成しておき、この作成条件判定パラメータ分布のうちの選択された特性データに該当する作成条件判定パラメータＡｄｊのみを参照するようにしてもよい。
【０１９６】
以上のようにして作成条件判定パラメータ分布を作成したならば、標本作成条件推定処理に移り（ステップａ７）、実施の形態１と同様の処理（図７を参照）を行って各染色色素の最適基準スペクトルを取得する。ただし、実施の形態３の標本作成条件推定処理では、図７のステップｂ１５の処理として、図２４のステップｅ６で設定した作成条件判定パラメータ分布から作成条件判定パラメータＡｄｊを選出する。
【０１９７】
そして、色素量推定部５４７が、観察染色標本画像の各画素位置についてステップａ３で取得したスペクトル（吸光度ａ（ｘ，λ））をもとに、ステップａ７の標本作成条件推定処理で各染色色素についてそれぞれ決定した最適基準スペクトルを用いて観察染色標本の色素量を推定する（ステップａ９）。その後、画像表示処理に移り（ステップａ１１）、実施の形態１と同様の処理を行う。
【０１９８】
以上説明したように、実施の形態３によれば、標本の属性を示す属性値に応じて定められた特性データを属性値毎に記憶しておき、観察染色標本の属性値等に応じた特性データを選択することができる。そして、選択した特性データを用いて作成条件判定パラメータＡｄｊを選出し、観察染色標本の作成条件を推定することによって各染色色素の最適基準スペクトルを決定することができる。したがって、観察染色標本の色素量をより精度良く推定することができる。
【０１９９】
なお、上記した各実施の形態では、Ｈ＆Ｅ染色が施された染色標本を観察対象とする場合について説明し、Ｈ＆Ｅ染色された染色標本を対象としているため、色素Ｈ、色素Ｅおよび色素Ｒの色素量を推定する場合を例示した。これに対し、本発明は、他の染色色素で染色された標本に対しても同様に適用することができ、その色素量を推定できる。また、上記した各実施の形態における色素Ｒのように、標本自体が有する固有の色についても同様に扱うことができる。
【産業上の利用可能性】
【０２００】
以上のように、本発明の画像処理装置、画像処理方法および顕微鏡システムは、所定の染色色素で染色された染色標本の色素量を精度良く推定するのに適している。
【符号の説明】
【０２０１】
１，１ｂ      画像処理装置
１１      染色標本画像撮像部
１２，５２      操作部
１３，５３      表示部
１４，１４ｂ，５４      画像処理部
１４１，５４１      スペクトル取得部
１４２，５４２      標本作成条件推定処理部
１４３，５４３      解析領域設定部
１４４，５４４      特徴量取得部
１４５，５４５      作成条件推定部
１４６，５４６      基準スペクトル決定部
１４７，５４７      色素量推定部
１４８，１４８ｂ，５４８      表示画像生成部
１４９ｂ      色素量補正部
１５０ｂ      スペクトル合成部
１６，１６ｂ，５６      記憶部
１６１，１６１ｂ，５６１      画像処理プログラム
１６３      基準スペクトル情報
１７，１７ｂ      制御部
１７１，５７１      解析領域選択入力依頼部
１７２，５７２      作成条件入力依頼部
１７３，５７３      色素選択入力依頼部
１７５，１７５ｂ，５７５      画像表示処理部
１７７ｂ      色素量調整入力依頼部
５７６      染色標本属性入力依頼部
５７７      特性データ選択部
５７８      システム環境設定部
２      顕微鏡システム
３      観察部
３１      染色標本観察部
３３      染色標本画像撮像部
５      観察システム制御部
７      特性データ記憶部
【図１】
2011095225.tif 000024
【図２】
2011095225.tif 000025
【図３】
2011095225.tif 000026
【図４】
2011095225.tif 000027
【図５】
2011095225.tif 000028
【図６】
2011095225.tif 000029
【図７】
2011095225.tif 000030
【図８】
2011095225.tif 000031
【図９】
2011095225.tif 000032
【図１０】
2011095225.tif 000033
【図１１】
2011095225.tif 000034
【図１２】
2011095225.tif 000035
【図１３】
2011095225.tif 000036
【図１４】
2011095225.tif 000037
【図１５】
2011095225.tif 000038
【図１６】
2011095225.tif 000039
【図１７】
2011095225.tif 000040
【図１８】
2011095225.tif 000041
【図１９】
2011095225.tif 000042
【図２０】
2011095225.tif 000043
【図２１】
2011095225.tif 000044
【図２２】
2011095225.tif 000045
【図２３】
2011095225.tif 000046
【図２４】
2011095225.tif 000047
【図２５】
2011095225.tif 000048