特許 6194791 画像処理装置及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6194791

(24)【登録日】2017年8月25日

(45)【発行日】2017年9月13日

(54)【発明の名称】画像処理装置及びプログラム

(51)【国際特許分類】

   C12M 1/34 20060101AFI20170904BHJP

   G01N 33/48 20060101ALI20170904BHJP

   G01N 21/27 20060101ALI20170904BHJP

   G01N 33/49 20060101ALI20170904BHJP

   G06T 7/00 20170101ALI20170904BHJP

【ＦＩ】

   C12M1/34 B

   G01N33/48 Z

   G01N33/48 M

   G01N21/27 A

   G01N33/49 K

   G06T7/00 630

【請求項の数】15

【全頁数】20

(21)【出願番号】特願2013-271140(P2013-271140)

(22)【出願日】2013年12月27日

(65)【公開番号】特開2015-123047(P2015-123047A)

(43)【公開日】2015年7月6日

【審査請求日】2016年2月18日

(73)【特許権者】

【識別番号】000005496

【氏名又は名称】富士ゼロックス株式会社

(74)【代理人】

【識別番号】110000154

【氏名又は名称】特許業務法人はるか国際特許事務所

(72)【発明者】

【氏名】尾崎 良太

(72)【発明者】

【氏名】織田 英人

(72)【発明者】

【氏名】加藤 典司

【審査官】 小林 薫

(56)【参考文献】

【文献】 特開２００２−２１４２２８（ＪＰ，Ａ）

【文献】 特開２０１０−２１６９２０（ＪＰ，Ａ）

【文献】 特開２０１１−１７７０６３（ＪＰ，Ａ）

【文献】 特表２００３−５０４６２７（ＪＰ，Ａ）

【文献】 特開２０１３−０２７３６８（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｃ１２Ｍ １／００− ３／１０

Ｇ０１Ｎ ３３／４８−３３／９８

Ｇ０１Ｎ ２１／００−２１／０１；２１／１７−２１／６１

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

ＤＷＰＩ（Ｔｈｏｍｓｏｎ Ｉｎｎｏｖａｔｉｏｎ）

(57)【特許請求の範囲】

【請求項1】

対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段と、
前記推定された標的細胞数に基づいて、前記対象の検体から生成される複数の標本のそれぞれについて前記標的細胞が含まれているかを検出する対象とするか否かを定める情報を含む検出パラメータを設定する設定手段と、
前記検出パラメータに基づいて検出の対象とされた標本の撮像画像に対して、該検出パラメータに基づいて設定した判定対象領域に前記標的細胞が含まれるか否かを判定する判定手段と、を含む
画像処理装置。

【請求項2】

前記推定手段により推定された標的細胞数と、前記対象の検体から生成される複数の標本の数に基づき算出した１標本当たりの標的細胞数に基づいて、予め定められた閾値以上の標的細胞が含まれる標本の数を決定する決定手段と、
前記決定手段により決定した数の標本を前記複数の標本の中から選択する選択手段と、を含み、
前記設定手段は、前記選択手段により選択された標本の情報に基づいて前記検出パラメータを設定する
請求項１に記載の画像処理装置。

【請求項3】

前記設定手段は、前記決定手段により決定した標本の数と、前記算出した１標本当たりの標的細胞数との積に対する前記閾値の差が小さくなるにつれて、前記判定対象領域が多く設定されるように前記設定条件を設定する
請求項２に記載の画像処理装置。

【請求項4】

前記標本の厚さ、該標本の輝度、該標本の濃度に基づいて、該標本に含まれる標的細胞数を検出する難易度を示す検出難易度を算出する難易度算出手段を含み、
前記設定手段は、前記算出された検出難易度に基づいて前記検出パラメータを設定する
請求項１乃至３のいずれかに記載の画像処理装置。

【請求項5】

前記設定手段は、前記検出難易度が高い程、前記判定対象領域が多く設定されるように前記設定条件を設定する
請求項４に記載の画像処理装置。

【請求項6】

前記難易度算出手段は、前記標本の厚さの均一性が高い程、前記検出難易度を低く算出し、前記標本に含まれる核の輝度が低い程、前記検出難易度を低く算出し、前記標本に含まれる予め定めた濃度よりも薄い画素の割合が高い程、前記検出難易度を高く算出する
請求項４又は５に記載の画像処理装置。

【請求項7】

前記標的細胞が含まれると判定された互いに少なくとも一部が重なる複数の判定対象領域がある場合に、前記判定手段により前記標的細胞が含まれる可能性が最も高いと判定された判定対象領域に該複数の判定対象領域を統合する統合手段を含む
請求項１乃至６のいずれかに記載の画像処理装置。

【請求項8】

前記検体は母体血であり、
前記標的細胞は、有核赤血球であり、
前記検体の特徴は、母体の年齢、病歴、妊娠週数の少なくともいずれかを含み、
前記検体の特徴に基づく分類ごとに単位血液量当たりの基準の有核赤血球数が定められており、
前記推定手段は、前記対象の検体の特徴に応じた分類について定められた基準の有核赤血球数と、該対象の検体の血液量に基づいて、該対象の検体に含まれる有核赤血球数を推定する
請求項１乃至７のいずれかに記載の画像処理装置。

【請求項9】

対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段、
前記推定された標的細胞数に基づいて、前記対象の検体から生成される複数の標本のそれぞれについて前記標的細胞が含まれているかを検出する対象とするか否かを定める情報を含む検出パラメータを設定する設定手段、及び、
前記検出パラメータに基づいて検出の対象とされた標本の撮像画像に対して、該検出パラメータに基づいて設定した判定対象領域に前記標的細胞が含まれるか否かを判定する判定手段
としてコンピュータを機能させるためのプログラム。

【請求項10】

対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段と、
前記推定された標的細胞数に基づいて、前記対象の検体の撮像画像において前記標的細胞の領域が含まれているか否かを判定する判定対象領域の設定条件を定める情報を含む検出パラメータを設定する設定手段と、
前記検出パラメータに基づいて検出の対象とされた標本の撮像画像に対して、該検出パラメータに基づいて設定した判定対象領域に前記標的細胞が含まれるか否かを判定する判定手段と、を含む
画像処理装置。

【請求項11】

対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段、
前記推定された標的細胞数に基づいて、前記対象の検体の撮像画像において前記標的細胞の領域が含まれているか否かを判定する判定対象領域の設定条件を定める情報を含む検出パラメータを設定する設定手段、及び、
前記検出パラメータに基づいて検出の対象とされた標本の撮像画像に対して、該検出パラメータに基づいて設定した判定対象領域に前記標的細胞が含まれるか否かを判定する判定手段
としてコンピュータを機能させるためのプログラム。

【請求項12】

母体の年齢、病歴、妊娠週数の少なくともいずれかに基づいて、前記母体の血液に含まれる標的細胞数を推定する推定手段と、
前記推定された標的細胞数に基づいて、前記母体の血液の撮像画像から前記標的細胞数を検出する処理に関する検出パラメータを設定する設定手段と、を含む
画像処理装置。

【請求項13】

母体の年齢、病歴、妊娠週数の少なくともいずれかに基づいて、前記母体の血液に含まれる標的細胞数を推定する推定手段、及び、
前記推定された標的細胞数に基づいて、前記前記母体の血液の撮像画像から前記標的細胞を検出する処理に関する検出パラメータを設定する設定手段
としてコンピュータを機能させるためのプログラム。

【請求項14】

母体の年齢、病歴、妊娠週数の少なくともいずれかに応じた分類について定められた基準の有核赤血球数と、該母体の血液量とに基づいて、前記母体の血液に含まれる有核赤血球数を推定する推定手段と、
前記推定された有核赤血球数に基づいて、前記母体の血液の撮像画像から前記有核赤血球を検出する処理に関する検出パラメータを設定する設定手段と、を含む
画像処理装置。

【請求項15】

母体の年齢、病歴、妊娠週数の少なくともいずれかに応じた分類について定められた基準の有核赤血球数と、該母体の血液量とに基づいて、前記母体の血液に含まれる有核赤血球数を推定する推定手段、及び、
前記推定された有核赤血球数に基づいて、前記母体の血液の撮像画像から前記有核赤血球を検出する処理に関する検出パラメータを設定する設定手段
としてコンピュータを機能させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、画像処理装置及びプログラムに関する。

【背景技術】

【0002】

例えば、特許文献１では、母体血を検体、母体血中の有核赤血球を標的細胞として、検体の顕微鏡画像から抽出される核候補領域に基づいて判定対象領域を設定し、設定した判定対象領域に標的細胞が含まれるか否かを判定している。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１２−２５４０４２号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

本発明の目的は、対象の検体から標的細胞を検出する処理のパラメータを、対象の検体の特徴に応じて設定できる画像処理装置及びプログラムを提供することにある。

【課題を解決するための手段】

【0005】

請求項１に係る発明は、対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段と、前記推定された標的細胞数に基づいて、前記対象の検体の撮像画像から前記標的細胞を検出する処理に関する検出パラメータを設定する設定手段と、を含む画像処理装置である。

【0006】

請求項２に係る発明は、前記検出パラメータは、前記対象の検体から生成される複数の標本のそれぞれについて前記標的細胞が含まれているかを検出する対象とするか否かを定める情報、前記撮像画像において前記標的細胞領域が含まれているか否かを判定する判定対象領域の設定条件を定める情報を含み、前記検出パラメータに基づいて検出の対象とされた標本の撮像画像に対して、該検出パラメータに基づいて設定した判定対象領域に前記標的細胞が含まれるか否かを判定する判定手段を含む請求項１に記載の画像処理装置である。

【0007】

請求項３に係る発明は、前記推定手段により推定された標的細胞数と、前記対象の検体から生成される複数の標本の数に基づき算出した１標本当たりの標的細胞数に基づいて、予め定められた閾値以上の標的細胞が含まれる標本の数を決定する決定手段と、前記決定手段により決定した数の標本を前記複数の標本の中から選択する選択手段と、を含み、前記設定手段は、前記選択手段により選択された標本の情報に基づいて前記検出パラメータを設定する請求項２に記載の画像処理装置である。

【0008】

請求項４に係る発明は、前記設定手段は、前記決定手段により決定した標本の数と、前記算出した１標本当たりの標的細胞数との積に対する前記閾値の差が小さくなるにつれて、前記判定対象領域が多く設定されるように前記設定条件を設定する請求項３に記載の画像処理装置である。

【0009】

請求項５に係る発明は、前記標本の厚さ、該標本の輝度、該標本の濃度に基づいて、該標本に含まれる標的細胞数を検出する難易度を示す検出難易度を算出する難易度算出手段を含み、前記設定手段は、前記算出された検出難易度に基づいて前記検出パラメータを設定する請求項２乃至４のいずれかに記載の画像処理装置である。

【0010】

請求項６に係る発明は、前記設定手段は、前記検出難易度が高い程、前記判定対象領域が多く設定されるように前記設定条件を設定する請求項５に記載の画像処理装置である。

【0011】

請求項７に係る発明は、前記難易度算出手段は、前記標本の厚さの均一性が高い程、前記検出難易度を低く算出し、前記標本に含まれる核の輝度が低い程、前記検出難易度を低く算出し、前記標本に含まれる予め定めた濃度よりも薄い画素の割合が高い程、前記検出難易度を高く算出する請求項５又は６に記載の画像処理装置である。

【0012】

請求項８に係る発明は、前記標的細胞が含まれると判定された互いに少なくとも一部が重なる複数の判定対象領域がある場合に、前記判定手段により前記標的細胞が含まれる可能性が最も高いと判定された判定対象領域に該複数の判定対象領域を統合する統合手段を含む請求項２乃至７のいずれかに記載の画像処理装置である。

【0013】

請求項９に係る発明は、前記検体は母体血であり、前記標的細胞は、有核赤血球であり、前記検体の特徴は、母体の年齢、病歴、妊娠週数の少なくともいずれかを含み、前記検体の特徴に基づく分類ごとに単位血液量当たりの基準の有核赤血球数が定められており、前記推定手段は、前記対象の検体の特徴に応じた分類について定められた基準の有核赤血球数と、該対象の検体の血液量に基づいて、該対象の検体に含まれる有核赤血球数を推定する請求項１乃至８のいずれかに記載の画像処理装置である。

【0014】

請求項１０に係る発明は、対象の検体の特徴に基づいて、該対象の検体に含まれる標的細胞数を推定する推定手段と、前記推定された標的細胞数に基づいて、前記対象の検体の撮像画像から前記標的細胞を検出する処理に関する検出パラメータを設定する設定手段としてコンピュータを機能させるためのプログラムである。

【発明の効果】

【0015】

請求項１及び１０に記載の発明によれば、対象の検体から標的細胞を検出する処理のパラメータを、対象の検体の特徴に応じて設定できる。

【0016】

請求項２に記載の発明によれば、対象の検体から生成した複数の標本のうち処理対象とする標本の数と、各標本について設定する判定対象領域とを、対象の検体について推定される標的細胞数に応じて設定できる。

【0017】

請求項３に記載の発明によれば、処理対象とする標本に閾値以上の標的細胞が含まれるように選択できる。

【0018】

請求項４に記載の発明によれば、処理対象とする標本の中に含まれる標的細胞の数と閾値との差が小さい場合には判定対象とする領域数を増やして、検出逃れの発生を抑制できる。

【0019】

請求項５に記載の発明によれば、対象の検体から標的細胞を検出する処理のパラメータを、処理対象とする検体から標的細胞を検出することの難易度に応じて設定できる。

【0020】

請求項６に記載の発明によれば、処理対象とする検体から標的細胞を検出することの難易度が高い場合には判定対象とする領域数を増やして、検出逃れの発生を抑制できる。

【0021】

請求項７に記載の発明によれば、本構成を有さない場合と比較して、処理対象とする検体から標的細胞を検出することの難易度を精度良く算出できる。

【0022】

請求項８に記載の発明によれば、複数の判定対象領域を統合しない場合と比較して、同一対象についての複数の判定対象領域の判定結果に基づく精度の良い判定結果を得ることができる。

【0023】

請求項９に記載の発明によれば、本構成を有さない場合と比較して、対象の母体血に含まれる有核赤血球の数を精度良く推定できる。

【図面の簡単な説明】

【0024】

【図1】本実施形態に係る画像処理システムのシステム構成例を示す図である。

【図2】画像処理装置が備える機能の一例を示す機能ブロック図である。

【図3】検出パラメータ管理テーブルの一例を示す図である。

【図4】標的細胞候補領域の統合例を説明する図である。

【図5】第１の実施例に係る処理のフロー図である。

【図6】検出難易度の算出処理のフロー図である。

【図7】標的細胞の検出処理のフロー図である。

【図8】標的細胞領域の統合処理のフロー図である。

【図9】第２の実施例に係る処理のフロー図である。

【図10】第３の実施例に係る処理のフロー図である。

【発明を実施するための形態】

【0025】

以下、本発明を実施するための実施の形態（以下、実施形態という）を、図面に従って説明する。

【0026】

［１．システム構成の説明］
図１には、本実施形態に係る画像処理システム１のシステム構成例を示した。図１に示されるように画像処理システム１は、光学顕微鏡１０、画像処理装置２０、入力装置５０、表示装置６０を含み、画像処理装置２０と光学顕微鏡１０、入力装置５０及び表示装置６０とはデータ通信可能に接続されている。

【0027】

光学顕微鏡１０は、試料台に配置されたスライドグラス１１上の試料を、対物レンズ１２等の光学系を介してＣＣＤカメラ１４で撮像する。光学顕微鏡１０には、スライドグラス１１と対物レンズ１２との距離を変化させる焦準機構１３を備えており、スライドグラス１１上の試料を複数の異なる焦点距離で撮像するようになっている。本実施形態では、試料には、母体血をスライドグラス１１に塗布し、メイ・ギムザ染色を施したものを用いる。これにより、母体血中の胎児由来有核赤血球（ＮＲＢＣｓ）が青紫色に染色される。以下、ＮＲＢＣｓを標的細胞と称する。

【0028】

画像処理装置２０は、光学顕微鏡１０により撮像された撮像画像を取得するとともに、当該取得した撮像画像の中から標的細胞を検出する。例えば、画像処理装置２０は、光学顕微鏡１０により撮像された撮像画像の中に設定した判定対象領域に、標的細胞が含まれている可能性を示すスコア（例えば確率）を、標的細胞の識別条件を学習した識別器に基づき判定することとしてよい。なお、画像処理装置２０において行われる標的細胞の検出処理の詳細については後述する。

【0029】

入力装置５０は、例えばキーボードやマウス等のデバイスであり、ユーザから受け付けた操作を画像処理装置２０に入力する。例えば、画像処理装置２０は、表示装置６０に表示された画像に対し、ユーザが入力装置５０によって指定した画像領域の情報を標的細胞の正例、負例、又は他の特定細胞の画像特徴を学習するための学習情報として取得し、学習情報に基づいて標的細胞の識別条件（識別パラメータ）を識別器に学習させることとしてよい。

【0030】

表示装置６０は、例えば液晶表示装置６０等であり、画像処理装置２０による処理の結果に基づいて画面を表示する。例えば、表示装置６０には、光学顕微鏡１０で撮像された撮像画像や、画像処理装置２０による標的細胞の検出結果等が表示される。

【0031】

［２．画像処理装置２０が備える機能の説明］
次に、本実施形態に係る画像処理装置２０に備えられる機能について説明する。

【0032】

図２には、画像処理装置２０が備える機能の一例を示す機能ブロック図を示した。図２に示されるように、画像処理装置２０は、検体情報取得部２１、検体特徴情報記憶部２２、標的細胞数推定部２３、撮像画像取得部２４、標本特徴抽出部２５、検出難易度判定部２６、検出対象標本設定部２７、検出パラメータ設定部２８、核候補領域抽出部２９、判定対象領域設定部３０、画像特徴生成部３１、標的細胞判定部３２、標的細胞領域設定部３３、標的細胞領域統合部３４、判定結果出力部３５を備える。

【0033】

画像処理装置２０に備えられる上記の各部の機能は、ＣＰＵ等の制御手段、メモリ等の記憶手段、外部デバイスとデータを送受信する入出力手段等を備えたコンピュータである画像処理装置２０が、コンピュータ読み取り可能な情報記憶媒体に格納されたプログラムを読み込み実行することで実現されるものとしてよい。なお、プログラムは光ディスク、磁気ディスク、磁気テープ、光磁気ディスク、フラッシュメモリ等の情報記憶媒体によって画像処理装置２０に供給されることとしてもよいし、インターネット等のデータ通信網を介して画像処理装置２０に供給されることとしてもよい。以下、画像処理装置２０が備える各部の機能の詳細について説明する。

【0034】

検体情報取得部２１は、検査対象とする検体（母体）の情報（検体情報）を取得する。例えば、検体情報取得部２１は、入力装置５０を介して入力されたデータに基づいて、検体情報を取得することとしてよい。例えば、検体情報には、母体血を採取した被検査対象である人の年齢、病歴、妊娠週数等を含むこととしてよい。

【0035】

検体特徴情報記憶部２２は、検体の特徴に基づいて、検体に含まれる標的細胞を推定する際に用いる情報を記憶する。例えば、検体特徴情報記憶部２２は、検体の特徴に応じた分類に対して、例えば分類ごとの単位血液量当たりの標的細胞数を記憶することとしてよい。例えば、画像処理装置２０は、年齢、病歴、妊娠週数の各情報に基づいて検体ごとの特徴ベクトルを生成してクラスタリングし、各クラスタに属する検体について予め測定された単位血液量当たりの標的細胞（有核赤血球）数の代表値（例えば平均値）を代表有核赤血球数として検体特徴情報記憶部２２に関連付けて記憶することとしてよい。また、検体特徴情報記憶部２２には、年齢、病歴、妊娠週数と、対応する有核赤血球数との関係を定めたテーブルや数式が記憶されることとしてもよい。

【0036】

標的細胞数推定部２３は、検体情報取得部２１により取得された検体に関する検体情報（検体の特徴を表す情報）と、検体特徴情報記憶部２２に記憶される情報とに基づいて、当該検体に含まれる標的細胞数を推定する。例えば、標的細胞数推定部２３は、検体情報取得部２１により取得された検体情報に対応する分類について、検体特徴情報記憶部２２に記憶される代表標的細胞数に基づいて、検体から得られた試料（母体血）に含まれる標的細胞数を推定する。例えば、標的細胞数推定部２３は、検体情報取得部２１により取得された検体情報に対応する分類について、検体特徴情報記憶部２２に記憶される単位血液量当たりの有核赤血球数と試料の血液量とを乗じることで推定標的細胞数を得ることとしてよい。

【0037】

また、標的細胞数の推定方法は上記に限られない。例えば、検体特徴情報記憶部２２は、予め１以上の検体について、検体の特徴ベクトルと標的細胞数を関連付けて記憶しておき、標的細胞数推定部２３は、検体特徴情報記憶部２２に記憶される検体の特徴ベクトルの中から、対象の検体の特徴ベクトルと類似する（例えば最も特徴ベクトル間の距離が近い）特徴ベクトルを特定し、特定した特徴ベクトルに関連付けられている標的細胞数を対象の検体の標的細胞数としてもよい。

【0038】

撮像画像取得部２４は、検体から得られたスライドグラス上の試料（母体血）を撮像した撮像画像を取得する。例えば、撮像画像取得部２４は、光学顕微鏡１０２に備えられたＣＣＤカメラ５より試料（母体血）を撮像した撮像画像を光学顕微鏡１０２から取得する。また、撮像画像取得部２４は、後述する検出パラメータに基づいて、対象の検体に定められた検出対象フラグが真（Ｔ）である場合には、撮像画像を取得し、対象の検体に定められた検出対象フラグが偽（Ｆ）である場合には、撮像画像を取得しないこととしてよい。

【0039】

標本特徴抽出部２５は、対象の標本の特徴を抽出する。標本の特徴は、標本の作成状態を示す情報を含み、例えば、標本特徴抽出部２５は、検体から生成した標本の厚さの均一性、細胞核の染色度合い、白血球の含有率等を標本特徴として抽出することとしてよい。標本特徴抽出部２５は、対象の検体から生成した標本の全て又は一部（例えば検出対象とする標本に限定する等）の標本のそれぞれについて標本特徴を抽出することとしてよい。

【0040】

標本特徴抽出部２５は、例えば、標本の厚さの均一性（Ａ）を、標本の複数の点に関して、標本の表面からスライドグラスまでの深度を光学顕微鏡１０で測定し、測定された複数の深度の分散値の逆数として算出することとしてよい。

【0041】

標本特徴抽出部２５は、例えば、標本の細胞核の染色度合い（Ｂ）を、標本の撮像画像（光学顕微鏡１０で撮像することとしてよい）の平均輝度値で、予め定めた閾値を除した値（すなわち、閾値／平均輝度値）として算出してもよいし、標本の撮像画像において、予め定めた閾値よりも輝度値が低い画素の割合として算出してもよい。

【0042】

標本特徴抽出部２５は、例えば、標本の白血球の含有率（Ｃ）を、標本の撮像画像（光学顕微鏡１０で撮像することとしてよい）における、予め定めた濃度よりも薄い色の画素の割合として算出してもよい。

【0043】

検出難易度判定部２６は、標本特徴抽出部２５により抽出された標本特徴に基づいて、標本から標的細胞を検出する難易度（Ｄｆ）を判定する。例えば、検出難易度判定部２６は、標本についての検出難易度Ｄｆを、標本特徴抽出部２５で抽出された標本の厚さの均一性Ａ、標本の細胞核の染色度合いＢ、標本の白血球の含有率Ｃに基づいて、以下の式（１）により算出することとしてよい。なお、ｗ_１，ｗ_２，ｗ_３は、ｗ_１＋ｗ_２＋ｗ_３＝１を満たす、それぞれ０以上の予め定められた係数としてよく、Ａ_０は、標本の厚さの均一性の基準値、Ｂ_０は、標本の細胞核の染色度合いの基準値、Ｃ_０は、標本の白血球の含有率の基準値としてよい。
Ｄｆ＝ｗ_１・Ａ_０／Ａ＋ｗ_２・Ｂ_０／Ｂ＋ｗ_３・Ｃ／Ｃ_０ （１）

【0044】

上式（１）によれば、検出難易度Ｄｆは、標本の厚さの均一性Ａの値が大きくなるにつれ低く算出され、標本の細胞核の染色度合いＢの値が大きくなるにつれ低く算出され、標本の白血球の含有率Ｃの値が大きくなるにつれ高く算出される。なお、検出難易度判定部２６は、上式（１）により算出された検出難易度Ｄｆが０から１までの値を取るように補正してもよい。例えば、検出難易度判定部２６は、Ｄｆが下限閾値より低い場合には０に、Ｄｆが上限閾値（例えば１より大きい値）よりも大きい場合には１とすることとしてもよい。

【0045】

検出対象標本設定部２７は、標的細胞数推定部２３により推定された推定標的細胞数（Ｃｂ）、又は検出難易度判定部２６により判定された検出難易度（Ｄｆ）の少なくともいずれかに基づいて、対象の検体から生成する複数の標本のうち検出対象とする標本を設定する。

【0046】

検出対象標本設定部２７は、例えば、対象の検体について推定された推定標的細胞数（Ｃｂ）を試料から生成する標本数（Ｎ）で除して、１標本当たりの標的細胞数（ａ＝Ｃｂ／Ｎ）を算出する。ここで、検出対象標本設定部２７は、必要標的細胞数（Ｚ）に対して、Ｃｂ＞Ｚである場合には、ａ・Ｘ≧Ｚとなる整数Ｘを検出対象標本数に決定し、Ｃｂ≦Ｚである場合には、検出対象標本数をＮ（全標本）に決定することとしてよい。そして、検出対象標本設定部２７は、試料から生成したＮ個の標本の中から、検出対象標本数の標本を選択して（例えば、識別番号順、生成日時順等）、選択した標本の検出対象フラグを真（Ｔ）に、選択しなかった標本の検出対象フラグを偽（Ｆ）に設定することとしてよい。

【0047】

また、検出対象標本設定部２７は、検出難易度判定部２６により算出された検出難易度に基づいて検出対象とする標本を設定する場合には、例えば、対象の検体から生成された複数の標本のうち検出難易度が閾値以下の標本を検出対象に設定することとしてもよい。

【0048】

検出パラメータ設定部２８は、試料から生成した標本について検出処理時に利用する検出パラメータを設定する。例えば、検出パラメータ設定部２８は、検出対象標本設定部２７により設定された検出対象標本と、検出難易度判定部２６により判定された検出難易度との少なくとも一方に基づいて、検出パラメータを設定する。例えば、検出パラメータは、標本の検出対象の当否を示す検出対象フラグ、標本の撮像画像から核候補領域として抽出する画像領域の条件を示す核候補領域パラメータ、標本の撮像画像について設定する判定対象領域の設定条件を示す判定対象領域パラメータを含むこととしてよい。

【0049】

図３には、検出パラメータ設定部２８により設定される検出パラメータを格納する検出パラメータ管理テーブルの一例を示した。図３に示されるように、検出パラメータ管理テーブルには、検体を識別する検体ＩＤ、検体から得た複数の標本のそれぞれを識別する標本ＩＤ、各標本を検出対象とするか否かを示す検出対象フラグ、核候補領域パラメータ、判定対象領域パラメータが関連付けて記憶される。なお、核候補領域パラメータには、核候補領域として抽出する対象画素の色範囲と、連結画素数の閾値が含む。また、判定対象領域パラメータには、核候補領域において判定対象領域の基点とする画素間のシフト量を示すステップ幅、判定対象領域の基点とする画素について設定する複数サイズの判定対象領域に関し、最小の判定対象領域と最大の判定対象領域とのサイズ比を示す最大倍率、最小の判定対象領域と最大の判定対象領域間までに何段階の拡大を実行するかを示す倍率段階を含む。

【0050】

検出パラメータ設定部２８は、例えば、検出対象標本設定部２７で設定した検出対象標本数Ｘと１標本当たりの標的細胞数ａの積と、必要標的細胞数Ｚとの差である余裕度Ｙ（Ｙ＝ａ・Ｘ−Ｚ）と、検出難易度判定部２６により判定された検出難易度Ｄｆとの少なくとも一方に基づいて、検出パラメータを設定することとしてよい。より具体的には、検出パラメータ設定部２８は、Ｍ（Ｍは２以上の整数）段階のレベル（Ｌ_１〜Ｌ_Ｍ）ごとに検出パラメータを定めておき（ただし同一画像についてＬ_ｉ＋１における検出パラメータにより抽出・設定される核候補領域と判定対象領域の数が、Ｌ_ｉにおける検出パラメータにより抽出・設定される核候補領域と判定対象領域の数よりも多い）、余裕度Ｙと検出難易度Ｄｆの少なくとも一方の値に基づいてレベルを決定し、決定したレベルに基づいて検出パラメータを設定することとしてよい。例えば、検出パラメータ設定部２８は、Ｍ段階のレベルごとに余裕度の範囲を定めておき、上記計算されたＹがどのレベルの範囲に属するかに基づいてレベルを決定してもよい。また、例えば、検出パラメータ設定部２８は、Ｍ段階のレベルごとに検出難易度の範囲を定めておき、上記計算されたＤｆがどのレベルの範囲に属するかに基づいてレベルを決定してもよい。また、例えば、検出パラメータ設定部２８は、Ｍ段階のレベルごとに余裕度と検出難易度の和の範囲を定めておき、上記計算されたＹとＤｆの和がどのレベルの範囲に属するかに基づいてレベルを決定してもよい。

【0051】

なお、検出パラメータ設定部２８は、検出対象標本設定部２７で設定した標本のそれぞれについて、それぞれの標本について算出された検出難易度に基づいて標本ごとに異なる検出パラメータを設定することとしてよい。

【0052】

核候補領域抽出部２９は、処理対象の標本について撮像画像取得部２４により取得された撮像画像の中から、当該標本について検出パラメータ設定部２８により設定された核候補領域パラメータに基づいて、核候補領域を抽出する。例えば、核候補領域抽出部２９は、核候補領域パラメータに含まれる色範囲に含まれる画素を黒画素、色範囲に含まれない画素を白画素として二値化し、隣り合う黒画素を連結した連結画素群のうち、核候補領域パラメータに含まれる連結画素群の数以上の連結画素群を、核候補領域として抽出することとしてよい。

【0053】

判定対象領域設定部３０は、処理対象の標本について撮像画像取得部２４により取得された撮像画像に対して、当該標本について検出パラメータ設定部２８により設定された判定対象領域パラメータと、核候補領域抽出部２９により抽出された核候補領域とに基づいて、標的細胞の有無の判定対象とする判定対象領域を設定する。例えば、判定対象領域設定部３０は、核候補領域抽出部２９により抽出された核候補領域に含まれる１以上の画素のそれぞれに対して、当該画素を中心（又は基点）とする矩形領域を判定対象領域として設定する。ここで、判定対象領域設定部３０は、処理対象の標本について設定された判定対象領域パラメータに含まれるステップ幅（シフト量）に基づき、判定対象領域の基点となる核候補領域内の画素をステップ幅だけ順次移動させて、判定対象領域を順次設定することとしてよい。また、判定対象領域設定部３０は、判定対象領域の基点となる１つの画素に関し、処理対象の標本について設定された判定対象領域パラメータに含まれる最大倍率と倍率段階に基づき、判定対象領域のサイズを１〜最大倍率まで、倍率段階で定められた段階数に渡って変化させて複数の異なるサイズの判定対象領域を設定することとしてよい。

【0054】

画像特徴生成部３１は、判定対象領域設定部３０により設定された判定対象領域の画像特徴量を生成する。例えば、画像特徴生成部３１は、判定対象領域のＨＯＧ特徴量を算出し、ＨＯＧ特徴量を画像特徴量として得ることとしてよい。以下、二種のＨＯＧ特徴量の算出処理について具体的に説明する。

【0055】

画像特徴生成部３１は、判定対象領域内の画素ごとに輝度勾配方向と輝度勾配強度を求め、対象画像をＸ個のセルから成るＹ個のブロックに分割し、各ブロックを構成するセルごとに輝度勾配方向と輝度勾配強度から輝度勾配方向ヒストグラム（［第１勾配方向の値，第２勾配方向の値，・・・，第９勾配方向の値］）を求め、これらの２乗平均が１になるようにブロック単位で正規化を行う。そして、画像特徴生成部３１は、上記正規化した輝度勾配方向ヒストグラムをブロック内で結合して作成したＸ×９個の値をブロックの特徴量とし、さらに対象画像内のブロックをすべて結合して作成したＹ×Ｘ×９個の値を判定対象領域のＨＯＧ特徴量とする。

【0056】

また、画像特徴生成部３１は、判定対象領域内の画素ごとに輝度勾配方向と輝度勾配強度を求め、対象画像をＸ個のセルから成るＹ個のブロックに分割し、各ブロックを構成するセルごとに輝度勾配方向と輝度勾配強度から輝度勾配方向ヒストグラム（［第１勾配方向の値，第２勾配方向の値，・・・，第１８勾配方向の値］）を求める。そして、画像特徴生成部３１は、上記の輝度勾配方向ヒストグラムをブロック内で結合して作成したＸ×１８個の値をブロックの特徴量とし、さらに対象画像内のブロックをすべて結合して作成したＹ×Ｘ×１８個の値を判定対象領域のＣｅｌｌ−ＨＯＧ特徴量とする。

【0057】

標的細胞判定部３２は、判定対象領域の画像特徴量に基づいて、判定対象領域に標的細胞が含まれる可能性（信頼度）を判定する。例えば、予め標的細胞が写し出された画像領域の画像特徴量に基づいて標的細胞の識別条件（識別パラメータ）を識別器に学習させ、識別器による判定対象領域の画像特徴量の識別結果を得ることとしてよい。なお、識別器には、ＡｄａＢｏｏｓｔやＳＶＭ（サポートベクタマシン）等を利用することとしてよい。識別器は、判定対象領域の画像特徴量に基づいて、判定対象領域に標的細胞が含まれる可能性を示すスコア（信頼度）を出力する。例えば、識別器は、判定対象領域に含まれる細胞が標的細胞であれば正の値のスコアを、標的細胞でなければ負の値のスコアを出力することとしてよい。

【0058】

標的細胞領域設定部３３は、判定対象領域についての標的細胞判定部３２による判定結果に基づいて、標的細胞が含まれる候補領域（標的細胞領域）を設定する。例えば、標的細胞領域設定部３３は、判定対象領域のうち識別器により出力された信頼度が０以上であれば標的細胞の候補領域に設定することとしてよい。

【0059】

標的細胞領域統合部３４は、標的細胞領域設定部３３により設定された候補領域のうち、互いに少なくとも一部が重なる候補領域を１つに統合する。例えば、標的細胞領域統合部３４は、複数の候補領域が互いに重なっている場合は、複数の候補領域をそのうちの１つ（例えば信頼度が最大の領域）に統合することとしてよい。ここで、標的細胞領域統合部３４は、複数の候補領域が同じ核候補領域から設定した判定対象領域であれば重なっていると判断してもよいし、複数の候補領域のうち予め定めた閾値以上の面積が互いに重なっていれば複数の候補領域が重なっていると判断してもよい。また、標的細胞領域統合部３４は、統合後の候補領域は信頼度が一番高い候補領域のみを残してもよいし、重なっている複数の候補領域のすべてを含む領域を統合後の候補領域としてもよい。また、標的細胞領域統合部３４は、統合後の候補領域の信頼度には、統合された候補領域の信頼度を加算してもよいし、統合後の候補領域の信頼度に重なっている候補領域の数をかけたものを、統合後の候補領域の信頼度としてもよい。

【0060】

図４には、標的細胞候補領域の統合例を説明する図を示した。図４に示されるように、撮像画像Ｉの中に核候補領域７０が含まれており、核候補領域７０に対して互いに一部が重なる２つの標的細胞領域７１，７２が設定され、標的細胞領域７１についてスコアが３０、標的細胞領域７２についてスコアが１０であったとすると、標的細胞領域統合部３４は、標的細胞領域７１，７２のうちスコアが最大の標的細胞領域７１に標的細胞領域７２を統合し、標的細胞領域７１のスコアに標的細胞領域７２のスコアを加算して更新することとしてよい。

【0061】

判定結果出力部３５は、標的細胞領域統合部３４により統合した候補領域の情報を、試料についての標的細胞の検出結果として出力する。例えば、判定結果出力部３５は、標的細胞領域統合部３４により統合された候補領域を、その信頼度に基づいて並び替えてリスト表示する表示画面を表示装置６０に表示させることとしてよい。なお、判定結果出力部３５は、信頼度が指定された又は予め定められた閾値未満の候補領域については、リスト表示に含めないようにすることとしてよい。

【0062】

［３．処理の例示］
次に、図５乃至図１０に基づいて、本実施形態に係る画像処理システム１において行われる処理の一例について詳細に説明する。

【0063】

［３−１．第１の実施例］
まず、図５〜図８に示したフロー図に基づき、画像処理システム１により実行される第１の実施例に係る処理の詳細について説明する。

【0064】

［３−１−１（１）．メイン処理（１）］
図５に示されるように、画像処理装置２０は、入力装置５０から入力されたデータに基づいて、対象の検体から生成した標本の数（Ｎ（Ｎは１以上の整数）とする）を設定するとともに（Ｓ１０１）、対象の検体の情報（検体情報）を取得する（Ｓ１０２）。ここで、検体情報には、年齢、病歴、妊娠週数の情報の少なくともいずれかを含むこととしてよい。

【0065】

次に、画像処理装置２０は、Ｓ１０２で取得した検体情報に基づいて、対象の検体に含まれる標的細胞数（ａ）を推定する（Ｓ１０３）。例えば、画像処理装置２０は、Ｓ１０２で取得した検体情報が属する分類を特定し、特定した分類について予め定められた単位血液当たりの代表標的細胞数（基準標的細胞数）と、検体となる試料の血液量とに基づいて、検体に含有される推定標的細胞数を算出することとしてよい。

【0066】

画像処理装置２０は、Ｓ１０３で算出した標的細胞数（ａ）と、Ｓ１０１で設定した標本数（Ｎ）とに基づいて、１標本当たりの標的細胞数（ｂ）を算出する（Ｓ１０４）。そして、画像処理装置２０は、Ｓ１０４で算出した１標本当たりの標的細胞数（ｂ）に基づいて、検出に必要とされる必要細胞数（ｃ）以上となる、検出対象標本数（ｄ）を決定する（Ｓ１０５）。

【0067】

ここで、画像処理装置２０は、Ｓ１０５で決定した検出対象標本数（ｄ）の処理対象とする標本を、対象の検体から生成された複数の標本の中から選択して決定する（Ｓ１０６）。

【0068】

次に、画像処理装置２０は、Ｓ１０６で決定した処理対象とする標本から標的細胞を検出する難易度を示す検出難易度を算出する処理を実行する（Ｓ１０７）。検出難易度の算出処理の詳細については、図６に示したフロー図に基づき説明する。

【0069】

［３−１−２．検出難易度の算出処理］
図６に示されるように、画像処理装置２０は、変数ｉを１に初期化して（Ｓ２０１）、標本Ｓ_ｉの深度（Ａ）を測定する（Ｓ２０２）。例えば、標本Ｓ_ｉの複数の点に関して、標本Ｓ_ｉの表面からスライドグラスまでの深度を光学顕微鏡１０で測定することとしてよい。

【0070】

次に、画像処理装置２０は、標本Ｓｉを光学顕微鏡１０により撮像した撮像画像を取得し（Ｓ２０３）、取得した撮像画像に基づいて、標本Ｓ_ｉの染色度（Ｂ）を取得する（Ｓ２０４）。例えば、画像処理装置２０は、標本Ｓ_ｉの細胞核の染色度合い（Ｂ）を、標本Ｓ_ｉの撮像画像の平均輝度値で、予め定めた閾値を除した値（すなわち、閾値／平均輝度値）として算出してもよいし、標本Ｓ_ｉの撮像画像において、予め定めた閾値よりも輝度値が低い画素の割合として算出してもよい。

【0071】

さらに、画像処理装置２０は、Ｓ２０３で取得した撮像画像に基づいて、標本Ｓ_ｉ中の白血球の割合（Ｃ）を取得する（Ｓ２０５）。例えば、画像処理装置２０は、標本Ｓ_ｉの白血球の割合（Ｃ）を、標本Ｓ_ｉの撮像画像における、予め定めた濃度よりも薄い色の画素の割合として算出してよい。

【0072】

画像処理装置２０は、例えば、上述した式（１）に基づいて、深度（Ａ）の分散値、試料の染色度（Ｂ）、白血球割合（Ｃ）に基づいて、標本Ｓ_ｉから標的細胞の検出難易度Ｄｆ（Ｄｆ_ｉ）を算出する（Ｓ２０６）。

【0073】

ここで、画像処理装置２０は、変数ｉがｄに達していない場合には（Ｓ２０７：Ｎ）、変数ｉをインクリメント（１加算）して（Ｓ２０８）、Ｓ２０２に戻り、変数ｉがｄに達している場合には（Ｓ２０７：Ｙ）、リターンする。次に、図５のフロー図に戻り説明を続ける。

【0074】

［３−１−１（２）．メイン処理（２）］
図５に示されるように、画像処理装置２０は、処理対象の標本について標的細胞の検出難易度Ｄｆの算出を終えると、Ｓ１０７で算出された検出難易度Ｄｆと、Ｓ１０６で決定した処理対象の標本の情報との少なくとも１つに基づいて、検出パラメータを設定する（Ｓ１０８）。例えば、検出パラメータは、標本の検出対象の当否を示す検出対象フラグ、標本から核候補領域として抽出する画像領域の条件を示す核候補領域パラメータ、標本の画像について設定する判定対象領域の設定条件を示す判定対象領域パラメータを含むこととしてよい。ここで、画像処理装置２０は、例えば、検出対象フラグについては、Ｓ１０６で決定した処理対象の標本の情報に基づいて設定し、核候補領域パラメータと判定対象領域パラメータについては標本Ｓ_ｉについて算出した検出難易度Ｄｆ_ｉに基づいて設定することとしてよい。

【0075】

画像処理装置２０は、Ｓ１０８で設定した検出パラメータに基づいて、検体から標的細胞を検出する処理を実行する（Ｓ１０９）。標的細胞の検出処理の詳細については、図７に示したフロー図に基づき説明する。

【0076】

［３−１−３．標的細胞の検出処理］
図７に示されるように、画像処理装置２０は、変数ｉ，ｊ，ｋの値をそれぞれ１に初期化して（Ｓ３０１）、対象の検体についての標本Ｓ_ｉが検出対象であるか否かを判断する（Ｓ３０２）。例えば、画像処理装置２０は、標本Ｓ_ｉの検出対象フラグが真（Ｔ）であれば検出対象と判断し、偽（Ｆ）であれば検出対象でないと判断することとしてよい。

【0077】

画像処理装置２０は、標本Ｓ_ｉが検出対象である場合には（Ｓ３０２：Ｙ）、例えば光学顕微鏡１０に標本Ｓ_ｉを撮像させて、標本Ｓ_ｉの撮像画像Ｔ_ｉを取得する（Ｓ３０３）。

【0078】

次に、画像処理装置２０は、標本Ｓ_ｉについて設定された検出パラメータに含まれる核候補領域パラメータに基づいて、撮像画像Ｔ_ｉから核候補領域（Ａ_ｉ１〜Ａ_ｉＭ）を抽出する（Ｓ３０４）。なお、Ｍは撮像画像Ｔ_ｉに含まれる核候補領域の数である。なお、核候補領域の抽出は、上述した核候補領域抽出部２９により実行することとしてよい。

【0079】

画像処理装置２０は、標本Ｓ_ｉについて設定された検出パラメータに含まれる判定対象領域パラメータに基づいて、撮像画像Ｔ_ｉから抽出した核候補領域Ａ_ｉｊに対して、判定対象領域（Ｂ_ｉｊ１〜Ｂ_ｉｊＬ）を設定する（Ｓ３０５）。なお、Ｌは核候補領域Ａ_ｉｊについて設定される判定対象領域の数であり、判定対象領域の設定は、上述した判定対象領域設定部３０により実行することとしてよい。

【0080】

画像処理装置２０は、判定対象領域Ｂ_ｉｊｋの画像特徴量Ｖ_ｉｊｋを算出する（Ｓ３０６）。なお、画像特徴量の算出は、上述した画像特徴生成部３１により実行することとしてよい。

【0081】

画像処理装置２０は、Ｓ３０６で算出した画像特徴量Ｖ_ｉｊｋに基づいて、判定対象領域Ｂ_ｉｊｋに標的細胞が含まれる可能性を示すスコアＰ_ｉｊｋを算出する（Ｓ３０７）。

【0082】

ここで、画像処理装置２０は、Ｓ３０７で算出したスコアＰ_ｉｊｋが閾値Ｔｈ以上であれば（Ｓ３０８：Ｙ）、判定対象領域Ｂ_ｉｊｋを標的細胞領域に設定する（Ｓ３０９）。

【0083】

画像処理装置２０は、Ｓ３０９の後、又はＳ３０７で算出したスコアＰ_ｉｊｋが閾値Ｔｈ未満である場合には（Ｓ３０８：Ｎ）、変数ｋがＬに達しているか否かを判断し（Ｓ３１０）、変数ｋがＬに達していない場合には（Ｓ３１０：Ｎ）、ｋをインクリメント（ｋに１加算）して（Ｓ３１１）、Ｓ３０６に戻る。また、画像処理装置２０は、Ｓ３１０において、変数ｋがＬに達している場合には（Ｓ３１０：Ｙ）、変数ｊがＭに達しているか否かを判断する（Ｓ３１２）。

【0084】

画像処理装置２０は、Ｓ３１２において変数ｊがＭに達してない場合には（Ｓ３１２：Ｎ）、ｊをインクリメント（ｊに１加算）するとともに、ｋを１に初期化して（Ｓ３１３）、Ｓ３０５に戻る。また、画像処理装置２０は、Ｓ３１２において変数ｊがＭに達している場合には（Ｓ３１２：Ｙ）、標本Ｓｉについて設定された標的細胞領域の統合処理を実行する（Ｓ３１４）。ここで、標的細胞領域の統合処理の詳細については図８に示したフロー図に基づき説明する。

【0085】

［３−１−４．標的細胞領域の統合処理］
図８に示されるように、画像処理装置２０は、標本Ｓ_ｉの撮像画像Ｔ_ｉについて設定したそれぞれの標的細胞領域に対して、処理済みフラグを付与する（Ｓ４０１）。なお、処理済みフラグは、真偽値であり、処理済みが真（Ｔ）、未処理が偽（Ｆ）とし、Ｓ４０１で付与される処理済みフラグの初期値は偽（Ｆ）であるとする。

【0086】

次に、画像処理装置２０は、処理済みフラグの値が未処理の標的細胞領域を１つ選択する（Ｓ４０２）。以下、Ｓ４０２で選択された標的細胞領域をＡｓとする。

【0087】

画像処理装置２０は、撮像画像Ｔ_ｉについて設定した標的細胞領域の中から、Ｓ４０２で選択された標的細胞領域Ａｓと重なる他の標的細胞領域を抽出する（Ｓ４０３）。そして、画像処理装置２０は、Ｓ４０３で抽出した標的細胞領域のうち、最大のスコアを有する標的細胞領域を選択する（Ｓ４０４）。以下、Ｓ４０４で選択された標的細胞領域をＡｔとする。

【0088】

画像処理装置２０は、標的細胞領域Ａｓのスコアが標的細胞領域Ａｔのスコア以上である場合には（Ｓ４０５：Ｙ）、標的細胞領域ＡｓにＳ４０３で抽出された標的細胞領域を統合して（Ｓ４０６）、標的細胞領域Ａｓのスコアに、統合する他の標的細胞領域のスコアを加算する（Ｓ４０７）。例えば、標的細胞領域Ａｓに対して統合される標的細胞領域は削除することとしてもよいし、統合済みの状態に更新することとしてもよい。そして、画像処理装置２０は、標的細胞領域Ａｓと、標的細胞領域Ａｓに統合された他の標的細胞領域のそれぞれの処理済みフラグを真（Ｔ）に更新する（Ｓ４０８）。

【0089】

また、画像処理装置２０は、標的細胞領域Ａｓのスコアが標的細胞領域Ａｔのスコア未満である場合には（Ｓ４０５：Ｎ）、標的細胞領域Ａｓを、標的細胞領域Ａｔに統合して（Ｓ４０９）、標的細胞領域Ａｔのスコアに、標的細胞領域Ａｓのスコアを加算する（Ｓ４１０）。そして、画像処理装置２０は、標的細胞領域Ａｓの処理済みフラグを真（Ｔ）に更新する（Ｓ４１１）。

【0090】

画像処理装置２０は、Ｓ４０８又はＳ４１１の後に、未処理の標的細胞領域が残っている場合には（Ｓ４１２：Ｙ）、Ｓ４０２に戻り、未処理の標的細胞領域が残っていない場合には（Ｓ４１２：Ｎ）、標的細胞の統合処理を終えてリターンする。

【0091】

ここで、図７に戻り説明を続ける。画像処理装置２０は、標本Ｓ_ｉについて設定された標的細胞領域の統合処理を終えると、変数ｉがｄ（検出対象標本数又は全標本数）に達していない場合には（Ｓ３１５：Ｎ）、ｉをインクリメント（１加算）するとともに、ｊ及びｋを１に初期化して（Ｓ３１６）、Ｓ３０２に戻る。また、画像処理装置２０は、Ｓ３１５において、変数ｉがｄに達している場合には（Ｓ３１５：Ｙ）、リターンする。

【0092】

また、画像処理装置２０は、Ｓ３０２において、標本Ｓｉが検出対象でない場合であって（Ｓ３０２：Ｎ）、変数ｉがｄに達していないときには（Ｓ３１７：Ｎ）、ｉをインクリメント（１加算）するとともに、ｊ及びｋを１に初期化して（Ｓ３１６）、Ｓ３０２に戻り、変数ｉがｄに達しているときには（Ｓ３１７：Ｙ）、リターンする。

【0093】

［３−１−１（３）．メイン処理（３）］
ここで、図５のフロー図に戻り説明を続ける。画像処理装置２０は、処理対象の検体のそれぞれについて標的細胞の検出処理を終えると、各検体について統合後の標的細胞領域に基づいて、標的細胞領域の情報を出力して（Ｓ１１０）、処理を終了する。例えば、画像処理装置２０は、検体ごとに、統合後の標的細胞領域のスコアが大きい順に標的細胞領域をソートして標的細胞領域を表示装置６０に表示させることとしてよい。

【0094】

以上説明した第１の実施例に係る処理によれば、処理対象とする検体に含まれる推定の標的細胞数と、処理対象とする検体から標的細胞を検出する難易度に基づいて、検出処理に利用するパラメータを設定できる。

【0095】

［３−２．第２の実施例］
次に、図９に示したフロー図に基づき、画像処理システム１により実行される第２の実施例に係る処理の詳細について説明する。第２の実施例においては、検出難易度を用いずに検出パラメータを設定する点で第１の実施例と異なっている。

【0096】

図９に示されるように、画像処理装置２０は、例えば、入力装置５０から入力されたデータに基づいて、対象の検体から生成した標本の数（Ｎ（Ｎは１以上の整数）とする）を設定するとともに（Ｓ５０１）、対象の検体の情報（検体情報）を取得する（Ｓ５０２）。ここで、検体情報には、年齢、病歴、妊娠週数の情報の少なくともいずれかを含むこととしてよい。

【0097】

次に、画像処理装置２０は、Ｓ５０２で取得した検体情報に基づいて、対象の検体に含まれる標的細胞数（ａ）を推定する（Ｓ５０３）。例えば、画像処理装置２０は、Ｓ５０２で取得した検体情報が属する分類を特定し、特定した分類について予め定められた単位血液当たりの代表標的細胞数と、検体となる試料の血液量とに基づいて、検体に含有される推定標的細胞数を算出することとしてよい。

【0098】

画像処理装置２０は、Ｓ５０３で算出した標的細胞数（ａ）と、Ｓ５０１で設定した標本数（Ｎ）とに基づいて、１標本当たりの標的細胞数（ｂ）を算出する（Ｓ５０４）。そして、画像処理装置２０は、Ｓ５０４で算出した１標本当たりの標的細胞数（ｂ）に基づいて、検出に必要とされる必要細胞数（ｃ）以上となる、検出対象標本数（ｄ）を決定する（Ｓ５０５）。

【0099】

ここで、画像処理装置２０は、Ｓ５０５で決定した検出対象標本数（ｄ）の処理対象とする標本を、対象の検体から生成された複数の標本の中から選択して決定する（Ｓ５０６）。

【0100】

次に、画像処理装置２０は、Ｓ５０６で決定した処理対象とする標本の情報と、Ｓ５０４で算出した１標本当たりの標的細胞数（ｂ）に基づいて、検出パラメータを設定する（Ｓ５０７）。例えば、検出パラメータは、標本の検出対象の当否を示す検出対象フラグ、標本から核候補領域として抽出する画像領域の条件を示す核候補領域パラメータ、標本の画像について設定する判定対象領域の設定条件を示す判定対象領域パラメータを含むこととしてよい。ここで、画像処理装置２０は、例えば、検出対象フラグについては、Ｓ５０６で決定した処理対象の標本の情報に基づいて設定し、核候補領域パラメータと判定対象領域パラメータについては、検出対象標本数（ｄ）と１標本当たりの標的細胞数を乗じたｂ・ｄから、必要細胞数（ｃ）を減じた余裕度（ｂ・ｄ−ｃ）に基づいて設定することとしてよい。

【0101】

画像処理装置２０は、Ｓ５０７で設定した検出パラメータに基づいて、検体から標的細胞を検出する処理を実行する（Ｓ５０８）。標的細胞の検出処理の詳細については、図７に示したフロー図に示した通りであり、第１の実施例と共通するためここでの説明は省略する。

【0102】

そして、画像処理装置２０は、処理対象の検体のそれぞれについて標的細胞の検出処理を終えると、各検体について統合後の標的細胞領域に基づいて、標的細胞領域の情報を出力して（Ｓ５０９）、処理を終了する。

【0103】

以上説明した第２の実施例に係る処理によれば、処理対象とする検体に含まれる推定の標的細胞数に基づいて、検出処理に利用するパラメータを設定できる。

【0104】

［３−３．第３の実施例］
次に、図１０に示したフロー図に基づき、画像処理システム１により実行される第３の実施例に係る処理の詳細について説明する。第３の実施例においては、対象の検体について推定される標的細胞数を用いずに検出パラメータを設定する点で第１の実施例と異なっている。

【0105】

図１０に示されるように、画像処理装置２０は、例えば、入力装置５０から入力されたデータに基づいて、対象の検体から生成した標本の数（Ｎ（Ｎは１以上の整数）とする）を設定する（Ｓ６０１）。

【0106】

ここで、画像処理装置２０は、対象の検体から生成された複数の標本の中から処理対象とする標本を決定する（Ｓ６０２）。例えば、画像処理装置２０は、対象の検体から生成された複数の標本の全てを処理対象としてもよいし、対象の検体について算出された検出難易度Ｄｆに基づいて処理対象とする標本を決定することとしてもよい。例えば、画像処理装置２０は、検出難易度の範囲に応じて割合を定めておき（具体的には、検出難易度が上がるにつれて割合を高くする）、検出難易度Ｄｆに対応する割合と全標本数との積に基づいて（具体的には得られた積を切り上げにより整数にすることとしてよい）、処理対象とする標本数を決定し、決定した標本数に基づいて処理対象の標本を決定することとしてよい。ここでは、処理対象とする標本数をｄで表す。

【0107】

次に、画像処理装置２０は、処理対象とする標本から標的細胞を検出する難易度を示す検出難易度を算出する処理を実行する（Ｓ６０３）。検出難易度の算出処理の詳細については、図６に示したフロー図の通りであり、第１の実施例と共通するためここでの説明は省略する。

【0108】

次に、画像処理装置２０は、処理対象とする標本Ｓ_ｉ（ｉ＝１〜ｄ）について、Ｓ６０３で算出された検出難易度Ｄｆ_ｉに基づいて、標本Ｓｉの検出パラメータを設定する（Ｓ６０４）。例えば、検出パラメータは、標本の検出対象の当否を示す検出対象フラグ、標本から核候補領域として抽出する画像領域の条件を示す核候補領域パラメータ、標本の画像について設定する判定対象領域の設定条件を示す判定対象領域パラメータを含むこととしてよい。ここで、画像処理装置２０は、例えば、検出対象フラグについては、Ｓ６０３で決定した処理対象の標本の情報に基づいて設定し、核候補領域パラメータと判定対象領域パラメータについては標本Ｓ_ｉについて算出した検出難易度Ｄｆ_ｉに基づいて設定することとしてよい。

【0109】

画像処理装置２０は、Ｓ６０４で設定した検出パラメータに基づいて、検体から標的細胞を検出する処理を実行する（Ｓ６０５）。標的細胞の検出処理の詳細については、図７に示したフロー図の通りであり、第１の実施例と共通するためここでの説明は省略する。

【0110】

そして、画像処理装置２０は、処理対象の検体のそれぞれについて標的細胞の検出処理を終えると、各検体について統合後の標的細胞領域に基づいて、標的細胞領域の情報を出力して（Ｓ６０６）、処理を終了する。

【0111】

以上説明した第３の実施例に係る処理によれば、処理対象とする検体から標的細胞を検出する難易度に基づいて、検出処理に利用するパラメータを設定できる。

【0112】

本発明は上記の実施形態に限定されるものではない。例えば、標的細胞領域を統合せずに、標的細胞領域をスコアの順に表示させることとしても構わない。また、画像処理装置２０は、光学顕微鏡１０から検体の撮像画像を取得する場合に限らず、他のコンピュータから検体の撮像画像を取得することとしてもよい。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】