(57)【特許請求の範囲】
【請求項１】
培養中の受精卵を撮影した画像を解析対象とする画像処理装置であって、
同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像群の中から任意の画像を基準画像とし、前記基準画像における正確な受精卵の輪郭を基準輪郭とし、前記基準輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を基準プロファイルとして保存する基準プロファイル作成手段と、
前記入力画像群の中から任意の画像を処理画像とし、前記基準輪郭に基づいて前記処理画像における受精卵の輪郭の候補の一つである第１の候補輪郭を算出し、前記第１の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第１の候補プロファイル作成手段と、
前記第１の候補輪郭の位置および／または大きさを変化させて前記処理画像における受精卵の輪郭の候補の一つである第２の候補輪郭を算出し、前記第２の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第２の候補プロファイル作成手段と、
前記基準プロファイルと前記候補プロファイルとの類似性を示す類似度を算出し、前記基準プロファイルと最も類似する候補プロファイルに基づいて、前記処理画像における受精卵の輪郭を決定する輪郭決定手段と、
を具備することを特徴とする画像処理装置。
【請求項２】
前記基準プロファイル作成手段、前記第１の候補プロファイル作成手段、および前記第２の候補プロファイル作成手段における分割領域は、輪郭の内側と外側のそれぞれに輪郭に沿う方向に複数定義され、かつ輪郭と垂直の方向にそれぞれ一定の幅を有するものであることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記輪郭決定手段における類似度は、輪郭を介して対向する内側の分割領域内と外側の分割領域内の画素値の統計量の差に基づいて算出することを特徴とする請求項２に記載の画像処理装置。
【請求項４】
前記一定の幅は、予め推定された受精卵を覆う透明帯の厚さの略半分とすることを特徴とする請求項２に記載の画像処理装置。
【請求項５】
前記一定の幅は、前記基準輪郭に基づく受精卵の大きさと、前記基準画像における脂肪小滴を閾値処理によって特定することで推定する受精卵中心部の大きさとに基づいて算出することを特徴とする請求項４に記載の画像処理装置。
【請求項６】
前記第２の候補プロファイル作成手段は、前記処理画像の撮像時刻と培養開始時刻との差が大きくなるにつれて、候補プロファイルのパターン数を増やすことを特徴とする請求項１に記載の画像処理装置。
【請求項７】
前記第１の候補プロファイル作成手段は、前記第１の候補輪郭の大きさを前記基準輪郭の大きさと同一と仮定し、前記第１の候補輪郭を移動させながら、内部の画素値に基づいて前記基準輪郭との比較処理を行うことで前記第１の候補輪郭の最適な位置を算出することを特徴とする請求項１に記載の画像処理装置。
【請求項８】
前記第１の候補プロファイル作成手段は、前記第１の候補輪郭の大きさを前記基準輪郭の大きさと同一と仮定し、前記第１の候補輪郭を移動させながら、前記第１の候補輪郭の内側と外側を含み、かつ輪郭と垂直の方向にそれぞれ一定の幅を有する領域の画素値に基づいて前記基準輪郭との比較処理を行うことで前記第１の候補輪郭の最適な位置を算出することを特徴とする請求項１に記載の画像処理装置。
【請求項９】
前記基準プロファイル作成手段は、前記基準画像を表示し、ユーザに対して前記基準画像における受精卵の輪郭を基準輪郭として指定させることを特徴とする請求項１に記載の画像処理装置。
【請求項１０】
前記基準プロファイル作成手段は、予め保持された輪郭検出済のテンプレート画像と前記基準画像とのテンプレートマッチング処理を行うことで推定する受精卵の輪郭を基準輪郭として指定することを特徴とする請求項１に記載の画像処理装置。
【請求項１１】
前記プロファイル比較手段は、更に、決定した受精卵の輪郭に基づいて出力画像を生成することを特徴とする請求項１に記載の画像処理装置。
【請求項１２】
前記プロファイル比較手段は、同一の受精卵に対する全ての時刻の入力画像に対して輪郭を決定した後、全ての時刻の輪郭を囲むことができる矩形を算出し、各入力画像において前記矩形の中心と受精卵の中心を一致させて、前記矩形に沿って入力画像を切り抜くことで出力画像を生成することを特徴とする請求項１１に記載の画像処理装置。
【請求項１３】
培養中の受精卵を撮影した画像を解析対象とする画像処理方法であって、
同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像群の中から任意の画像を基準画像とし、前記基準画像における正確な受精卵の輪郭を基準輪郭とし、前記基準輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を基準プロファイルとして保存する基準プロファイル作成ステップと、
前記入力画像群の中から任意の画像を処理画像とし、前記基準輪郭に基づいて前記処理画像における受精卵の輪郭の候補の一つである第１の候補輪郭を算出し、前記第１の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第１の候補プロファイル作成ステップと、
前記第１の候補輪郭の位置および／または大きさを変化させて前記処理画像における受精卵の輪郭の候補の一つである第２の候補輪郭を算出し、前記第２の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第２の候補プロファイル作成ステップと、
前記基準プロファイルと前記候補プロファイルとの類似性を示す類似度を算出し、前記基準プロファイルと最も類似する候補プロファイルに基づいて、前記処理画像における受精卵の輪郭を決定する輪郭決定ステップと、
を含むことを特徴とする画像処理方法。
【請求項１４】
コンピュータを請求項１に記載の画像処理装置として機能させるためのプログラム。
【請求項１５】
コンピュータを請求項１に記載の画像処理装置として機能させるためのプログラムを記憶した記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、培養中の受精卵を撮影した画像を解析対象とする画像処理装置等に関するものである。
【背景技術】
【０００２】
現在、家畜などの受精卵をインキュベータ内で培養し、移植可能な状態になるまで発育させることが行われている。品質が低い受精卵は移植しても受胎しにくい為、何らかの仕組みで受精卵の品質を評価する必要がある。牛の体外受精卵の場合、従来は、熟練者が培養７日目に受精卵をインキュベータから取り出し、顕微鏡で観察することで、受精卵の最終品質を評価していた。しかしながら、インキュベータから取り出して観察することは、受精卵にとって大きなストレスとなる。
【０００３】
特許文献１では、インキュベータから受精卵を取り出すことなく撮像された画像のみを用いて、受精卵の品質評価に有益な情報を提示する受精卵品質評価支援システムが開示されている。特許文献１に記載の技術では、受精卵が配置された培養容器を定期的に撮像し、撮像した画像から受精卵を含む最小近似円領域を切り出す。そして、切り出した画像の経時変化を解析することで、受精卵の品質評価に有益な情報を提示している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特願２０１０−００１０１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
一般に、解析対象が自律的に動く場合、画像から経時変化を解析する為には、複数の画像間における解析対象の位置合わせが必要不可欠である。
位置合わせを行う従来技術としては、解析対象の輪郭全体の形状を比較する方法、解析対象の模様の変わり目や輪郭線の一部の形状など特徴的な部分同士を比較する方法、テンプレートと呼ばれる見本画像と各画像とを比較し、一致度が最も高いテンプレートに基づいて位置合わせを行う方法などが存在する。
これらの従来技術は、解析対象の輪郭や模様の変わり目など位置合わせ処理の基準に用いられる境界が識別可能な画像に対しては十分な性能を発揮する。しかしながら、受精卵に対しては、いずれの従来技術も十分な性能を発揮することができない。
【０００６】
受精卵は略球体であり、受精卵の輪郭線が現れるべき球の表面部分は、ある程度の厚みを持った半透明、ゼラチン質の殻である透明帯によって構成されている。透明帯部分は、光が内部で複雑に反射や屈折を繰り返すので、背景よりも暗くなる部位、背景よりも明るくなる部位、背景と略同じ明るさになる部位が混在する。これらの部位には、撮像時の光の加減によって、背景と同じ明るさとなる不明りょうな部分が存在し、公知の画像処理技術によって透明帯と背景の境界（＝受精卵表面の輪郭）を検出することは困難である。
一方、受精卵の内部には、脂肪小滴と呼ばれる光を通しにくい部分（画像においては黒い点）が無数に含まれている。脂肪小滴の画素値は、他の部位の画素値と明確な差がある為、公知の画像処理技術によって検出し易い。しかし、脂肪小滴の分布は、受精卵の内部で均一ではない。また、各脂肪小滴は、それぞれランダムに移動し続けている。従って、受精卵画像における脂肪小滴は、位置合わせ処理の基準とすることはできない。そればかりか、前述の従来技術によって受精卵表面の輪郭を検出しようとすると、脂肪小滴に由来する受精卵の内部の輪郭が強すぎるために処理の妨げとなる。
【０００７】
本発明は、前述した問題点に鑑みてなされたもので、その目的とすることは、受精卵の経時変化を解析する際、複数の画像間における受精卵の位置を正確に合わせることができる画像処理装置等を提供することである。
【課題を解決するための手段】
【０００８】
前述した目的を達成するために第１の発明は、培養中の受精卵を撮影した画像を解析対象とする画像処理装置であって、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像群の中から任意の画像を基準画像とし、前記基準画像における正確な受精卵の輪郭を基準輪郭とし、前記基準輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を基準プロファイルとして保存する基準プロファイル作成手段と、前記入力画像群の中から任意の画像を処理画像とし、前記基準輪郭に基づいて前記処理画像における受精卵の輪郭の候補の一つである第１の候補輪郭を算出し、前記第１の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第１の候補プロファイル作成手段と、前記第１の候補輪郭の位置および／または大きさを変化させて前記処理画像における受精卵の輪郭の候補の一つである第２の候補輪郭を算出し、前記第２の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第２の候補プロファイル作成手段と、前記基準プロファイルと前記候補プロファイルとの類似性を示す類似度を算出し、前記基準プロファイルと最も類似する候補プロファイルに基づいて、前記処理画像における受精卵の輪郭を決定する輪郭決定手段と、を具備することを特徴とする画像処理装置である。
第１の発明によって、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像間における受精卵の位置を正確に合わせることができる。
【０００９】
第１の発明における前記基準プロファイル作成手段、前記第１の候補プロファイル作成手段、および前記第２の候補プロファイル作成手段における分割領域は、輪郭の内側と外側のそれぞれに輪郭に沿う方向に複数定義され、かつ輪郭と垂直の方向にそれぞれ一定の幅を有するものであることが望ましい。
これによって、輪郭に沿う方向に対する輪郭周辺の画素値の変化に基づいて輪郭を決定することができる。
【００１０】
また、第１の発明における前記輪郭決定手段における類似度は、輪郭を介して対向する内側の分割領域内と外側の分割領域内の画素値の統計量の差に基づいて算出することが望ましい。
これによって、輪郭の内側と外側のコントラストの強弱を用いて輪郭を特徴付けることができる。
【００１１】
また、第１の発明における前記一定の幅は、例えば、予め推定された受精卵を覆う透明帯の厚さの略半分とする。
これによって、内側の分割領域が受精卵中心部を含まないので、輪郭決定手段において、受精卵中心部の脂肪小滴が自律的に移動することによる影響をなくすことができる。
【００１２】
また、第１の発明における前記一定の幅は、例えば、前記基準輪郭に基づく受精卵の大きさと、前記基準画像における脂肪小滴を閾値処理によって特定することで推定する受精卵中心部の大きさとに基づいて算出する。
これによって、画像に応じた適切な値を一定の幅として設定することができる。
【００１３】
また、第１の発明における前記第２の候補プロファイル作成手段は、前記処理画像の撮像時刻と培養開始時刻との差が大きくなるにつれて、候補プロファイルのパターン数を増やすことが望ましい。
これによって、培養期間の後半における受精卵の大きさの変化の増大に応じて、候補プロファイルのパターン数を適切な数に設定することができる。
【００１４】
また、第１の発明における前記第１の候補プロファイル作成手段は、例えば、前記第１の候補輪郭の大きさを前記基準輪郭の大きさと同一と仮定し、前記第１の候補輪郭を移動させながら、内部の画素値に基づいて前記基準輪郭との比較処理を行うことで前記第１の候補輪郭の最適な位置を算出する。
これによって、受精卵の輪郭の候補の初期値である第１の候補輪郭を自動的に算出することができる。
【００１５】
また、第１の発明における前記第１の候補プロファイル作成手段は、例えば、前記第１の候補輪郭の大きさを前記基準輪郭の大きさと同一と仮定し、前記第１の候補輪郭を移動させながら、前記第１の候補輪郭の内側と外側を含み、かつ輪郭と垂直の方向にそれぞれ一定の幅を有する領域の画素値に基づいて前記基準輪郭との比較処理を行うことで前記第１の候補輪郭の最適な位置を算出する。
これによって、受精卵中心部の脂肪小滴が自律的に移動することによる影響を受けずに、受精卵の輪郭の候補の初期値である第１の候補輪郭を自動的に算出することができる。
【００１６】
また、第１の発明における前記基準プロファイル作成手段は、例えば、前記基準画像を表示し、ユーザに対して前記基準画像における受精卵の輪郭を基準輪郭として指定させる。
これによって、基準画像における基準輪郭を正確に指定することができる。
【００１７】
また、第１の発明における前記基準プロファイル作成手段は、予め保持された輪郭検出済のテンプレート画像と前記基準画像とのテンプレートマッチング処理を行うことで推定する受精卵の輪郭を基準輪郭として指定する。
これによって、基準画像における基準輪郭を自動的に指定することができる。
【００１８】
また、第１の発明における前記プロファイル比較手段は、更に、決定した受精卵の輪郭に基づいて出力画像を生成することが望ましい。これによって、解析内容に応じた出力画像（受精卵観測画像）を得ることができる。
【００１９】
第１の発明における前記プロファイル比較手段は、例えば、同一の受精卵に対する全ての時刻の入力画像に対して輪郭を決定した後、全ての時刻の輪郭を囲むことができる矩形を算出し、各入力画像において前記矩形の中心と受精卵の中心を一致させて、前記矩形に沿って入力画像を切り抜くことで出力画像を生成する。
これによって、受精卵周辺の部位の画像情報が損失せず、管理し易い出力画像を得ることができる。
【００２０】
第２の発明は、培養中の受精卵を撮影した画像を解析対象とする画像処理方法であって、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像群の中から任意の画像を基準画像とし、前記基準画像における正確な受精卵の輪郭を基準輪郭とし、前記基準輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を基準プロファイルとして保存する基準プロファイル作成ステップと、前記入力画像群の中から任意の画像を処理画像とし、前記基準輪郭に基づいて前記処理画像における受精卵の輪郭の候補の一つである第１の候補輪郭を算出し、前記第１の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第１の候補プロファイル作成ステップと、前記第１の候補輪郭の位置および／または大きさを変化させて前記処理画像における受精卵の輪郭の候補の一つである第２の候補輪郭を算出し、前記第２の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイルとして保存する第２の候補プロファイル作成ステップと、前記基準プロファイルと前記候補プロファイルとの類似性を示す類似度を算出し、前記基準プロファイルと最も類似する候補プロファイルに基づいて、前記処理画像における受精卵の輪郭を決定する輪郭決定ステップと、を含むことを特徴とする画像処理方法である。
第２の発明によって、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像間における受精卵の位置を正確に合わせることができる。
【００２１】
第３の発明は、コンピュータを第１の発明の画像処理装置として機能させるためのプログラムである。
ネットワークを介して第３の発明を配布し、汎用的なコンピュータにインストールすることで、第１の発明の画像処理装置を提供することができる。
【００２２】
第４の発明は、コンピュータを第１の発明の画像処理装置として機能させるためのプログラムを記憶した記憶媒体である。
汎用的なコンピュータが備える記憶媒体読取装置に第４の発明を装着し、汎用的なコンピュータに第４の発明が記憶するプログラムをインストールすることで、第１の発明の画像処理装置を提供することができる。
【発明の効果】
【００２３】
本発明により、受精卵の経時変化を解析する際、複数の画像間における受精卵の位置を正確に合わせることができる画像処理装置等を提供することができる。
【図面の簡単な説明】
【００２４】
【図１】１細胞卵の受精卵を撮像した画像を示す模式図
【図２】２細胞卵の受精卵を撮像した画像を示す模式図
【図３】１細胞卵の受精卵を撮像した画像に輪郭が重畳された状態を示す模式図
【図４】２細胞卵の受精卵を撮像した画像に輪郭が重畳された状態を示す模式図
【図５】画像処理装置１のハードウエア構成図
【図６】画像処理装置１のソフトウエア構成図
【図７】基準プロファイル作成処理を示すフローチャート
【図８】基準輪郭６２を指定する処理の第１の例を説明する図
【図９】基準輪郭６２を指定する処理の第２の例を説明する図
【図１０】分割領域６３を説明する図
【図１１】一定の幅を決定する処理を説明する図
【図１２】基準プロファイル６の例を示す図
【図１３】候補プロファイル作成処理を示すフローチャート
【図１４】候補輪郭７２の初期値を決定する処理の第１の例を説明する図
【図１５】候補輪郭７２の初期値を決定する処理の第２の例を説明する図
【図１６】候補プロファイル７の第１の例を示す図
【図１７】候補プロファイル７の第２の例を示す図
【図１８】候補プロファイル７の第３の例を示す図
【図１９】輪郭決定処理を示すフローチャート
【発明を実施するための形態】
【００２５】
以下図面に基づいて、本発明の実施形態を詳細に説明する。本発明の実施形態において参照する図面では、説明を分かりやすくする為に、解析対象である培養中の受精卵を撮像した画像を模式的に示す。しかしながら、実際の画像に対して本発明を適用できることは言うまでもない。
また、本発明の実施形態において模式的に示す画像は、白から黒までの明暗だけで表現するグレースケール画像であるが、本発明はいわゆるカラー画像に対しても同様に適用可能であることは言うまでもない。
【００２６】
最初に、図１から図４を参照しながら、受精卵を撮像した画像の特徴について説明する。
図１は、１細胞卵の受精卵を撮像した画像を示す模式図である。図１に示すように、受精卵の中心（受精卵中心部）には、主に脂肪小滴に由来する強い輪郭を伴った部分（図１においては黒に近い画素値を持つ領域）が存在する。一方、受精卵は、ある程度の厚みを持った透明帯（図１においては灰色の画素値を持つドーナッツ状の領域）が覆っている。透明帯には、撮像時の光の加減によって、受精卵の輪郭が背景と同じ明るさとなる不明りょうな部分がある。
【００２７】
ここで、本発明の実施形態における用語の定義について説明する。
「受精卵の輪郭」とは、画像（二次元）におけるドーナツ状の透明帯の外縁（二重円のうち外円）を意味するものとする。
「受精卵中心部」とは、画像（二次元）におけるドーナツ状の透明帯の内縁（二重円のうち内円）よりも内側の領域を意味するものとする。
「受精卵の大きさ」とは、透明帯と受精卵中心部の両方を含む大きさであり、画像（二次元）における透明帯の外縁を真円と仮定する場合には、半径、直径、または面積を意味するものとし、楕円と仮定する場合には、長径および短径、または面積を意味するものとする。
「受精卵中心部の大きさ」とは、画像（二次元）における受精卵中心部を真円と仮定する場合には、半径、直径、または面積を意味するものとし、楕円と仮定する場合には、長径および短径、または面積を意味するものとする。
【００２８】
図２は、２細胞卵の受精卵を撮像した画像を示す模式図である。図２に示す画像の撮像時刻は、図１に示す画像の撮像時刻よりも後である。時間が経過したことに伴い、図２の画像における受精卵中心部は、図１の画像と異なっており、受精卵の大きさも変化している。
特に、脂肪小滴に相当する領域は、図１の画像と大きく異なり、受精卵中心部の中で偏在している。従って、従来技術によって算出可能な領域の重心は、受精卵の中心点と一致しない。
一方、透明帯については、光源の方向が変化していない為、概ね図１の画像と同様の特徴（明るさの分布など）を保っている。
本発明の実施の形態では、図１と図２の画像のように、撮像時の光源の方向が変化していない画像同士を対象とするものとする。
【００２９】
図３は、１細胞卵の受精卵を撮像した画像に輪郭が重畳された状態を示す模式図、図４は、２細胞卵の受精卵を撮像した画像に輪郭が重畳された状態を示す模式図である。
図３、図４ともに、図１、図２の画像から検出が期待される受精卵の輪郭を示している。一般に、受精卵は球体をしている為、図３、図４ともに、円形の輪郭を重畳している。図３、図４に示す円の中心位置は、異なる撮像時刻において撮像された画像間において受精卵の位置を示す共通基準として使うことができる。受精卵の位置を正確に同定することができれば、受精卵の大きさの変化、受精卵内部の状態の変化などを精緻に解析することが可能となる。
尚、本発明の実施の形態では、受精卵を楕円体と仮定し、楕円形の輪郭を重畳することも可能である。
【００３０】
次に、図５から図１９を参照しながら、本発明の実施形態における画像処理装置について説明する。
以下では、混乱を避けるために、受精卵および受精卵中心部は、球体と仮定して説明する。
【００３１】
最初に、図５、図６を参照しながら、画像処理装置１の構成について説明する。
図５は、画像処理装置１のハードウエア構成図である。尚、図５のハードウエア構成は一例であり、用途、目的に応じて様々な構成を採ることが可能である。
【００３２】
画像処理装置１は、制御部１１、記憶部１２、メディア入出力部１３、通信制御部１４、入力部１５、表示部１６、周辺機器Ｉ／Ｆ部１７等が、バス１８を介して接続される。
【００３３】
制御部１１は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）等で構成される。
【００３４】
ＣＰＵは、記憶部１２、ＲＯＭ、記録媒体等に格納されるプログラムをＲＡＭ上のワークメモリ領域に呼び出して実行し、バス１８を介して接続された各装置を駆動制御し、画像処理装置１が行う後述する処理を実現する。
ＲＯＭは、不揮発性メモリであり、コンピュータのブートプログラムやＢＩＯＳ等のプログラム、データ等を恒久的に保持している。
ＲＡＭは、揮発性メモリであり、記憶部１２、ＲＯＭ、記録媒体等からロードしたプログラム、データ等を一時的に保持するとともに、制御部１１が各種処理を行う為に使用するワークエリアを備える。
【００３５】
記憶部１２は、ＨＤＤ（ハードディスクドライブ）であり、制御部１１が実行するプログラム、プログラム実行に必要なデータ、ＯＳ（オペレーティングシステム）等が格納される。プログラムに関しては、ＯＳ（オペレーティングシステム）に相当する制御プログラムや、後述する処理をコンピュータに実行させるためのアプリケーションプログラムが格納されている。
これらの各プログラムコードは、制御部１１により必要に応じて読み出されてＲＡＭに移され、ＣＰＵに読み出されて各種の手段として実行される。
【００３６】
メディア入出力部１３（ドライブ装置）は、データの入出力を行い、例えば、ＣＤドライブ（−ＲＯＭ、−Ｒ、−ＲＷ等）、ＤＶＤドライブ（−ＲＯＭ、−Ｒ、−ＲＷ等）等のメディア入出力装置を有する。
通信制御部１４は、通信制御装置、通信ポート等を有し、コンピュータとネットワーク間の通信を媒介する通信インタフェースであり、ネットワークを介して、他のコンピュータ間との通信制御を行う。ネットワークは、有線、無線を問わない。
【００３７】
入力部１５は、データの入力を行い、例えば、キーボード、マウス等のポインティングデバイス、テンキー等の入力装置を有する。
入力部１５を介して、コンピュータに対して、操作指示、動作指示、データ入力等を行うことができる。
表示部１６は、ＣＲＴモニタ、液晶パネル等のディスプレイ装置、ディスプレイ装置と連携してコンピュータのビデオ機能を実現するための論理回路等（ビデオアダプタ等）を有する。
【００３８】
周辺機器Ｉ／Ｆ（インタフェース）部１７は、コンピュータに周辺機器を接続させるためのポートであり、周辺機器Ｉ／Ｆ部１７を介してコンピュータは周辺機器とのデータの送受信を行う。周辺機器としては、例えば、受精卵を撮像するための撮像装置等がある。周辺機器Ｉ／Ｆ部１７は、ＵＳＢやＩＥＥＥ１３９４やＲＳ−２３２Ｃ等で構成されており、通常複数の周辺機器Ｉ／Ｆを有する。周辺機器との接続形態は有線、無線を問わない。
バス１８は、各装置間の制御信号、データ信号等の授受を媒介する経路である。
【００３９】
図６は、画像処理装置１のソフトウエア構成図である。画像処理装置１は、図６に示す各プログラムを記憶部１２に記憶し、各プログラムコードを制御部１１により必要に応じて読み出してＲＡＭに移し、ＣＰＵに読み出すことで、各種の手段として機能する。
【００４０】
画像処理装置１は、基準プロファイル作成プログラム２、候補プロファイル作成プログラム３、輪郭決定プログラム４等を実行可能なプログラムコードとして記憶部１２に記憶する。また、画像処理装置１は、処理の進行に応じて、入力画像５、基準プロファイル６、候補プロファイル７、出力画像８等を記憶部１２またはＲＡＭに記憶する。
【００４１】
入力画像５は、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた画像群である。
基準プロファイル作成プログラム２は、基準プロファイル作成処理を実現する為のプログラムである。基準プロファイル作成処理では、入力画像５の中から任意の画像を基準画像とし、基準画像における正確な受精卵の輪郭を基準輪郭とし、基準輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を基準プロファイル６として保存する。基準プロファイル作成処理については、図７を参照しながら後述する。
【００４２】
候補プロファイル作成プログラム３は、候補プロファイル作成処理を実現する為のプログラムである。候補プロファイル作成処理では、入力画像５の中から任意の画像を処理画像とし、基準輪郭に基づいて処理画像における受精卵の輪郭の候補の一つである初期候補輪郭（第１の候補輪郭）を算出し、初期候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイル７として保存する。更に、候補プロファイル作成処理では、初期候補輪郭の位置および／または大きさを変化させて処理画像における受精卵の輪郭の候補の一つであるその他の候補輪郭（第２の候補輪郭）を算出し、その他の候補輪郭に沿って複数の分割領域を定義し、分割領域内の画素値の統計量を算出し、算出結果を候補プロファイル７として保存する。候補プロファイル作成処理については、図１３を参照しながら後述する。
【００４３】
輪郭決定プログラム４は、輪郭決定処理を実現する為のプログラムである。輪郭決定処理では、基準プロファイル６と候補プロファイル７との類似性を示す類似度を算出し、基準プロファイル６と最も類似する候補プロファイル７に基づいて、処理画像における受精卵の輪郭を決定する。更に、輪郭決定処理では、決定した受精卵の輪郭に基づいて出力画像８（受精卵観測画像）を生成する。輪郭決定処理については、図１９を参照しながら後述する。
【００４４】
次に、図７から図１２を参照しながら、基準プロファイル作成処理について説明する。
図７は、基準プロファイル作成処理を示すフローチャートである。
図７に示すように、画像処理装置１の制御部１１は、入力画像５の中から基準画像を選択し（Ｓ１０１）、基準輪郭の指定を行う（Ｓ１０２）。
【００４５】
図８は、基準輪郭６２を指定する処理の第１の例を説明する図である。図８に示す例では、制御部１１は、表示部１６に基準画像６１を表示し、ユーザに対して基準画像６１における受精卵の輪郭を基準輪郭６２として選択させる。具体的には、ユーザは、入力部１５の一つであるマウスを用いて、マウスポインタ１９が指し示す円を拡大、縮小、移動させて、受精卵の輪郭（大きさおよび位置）と、マウスポインタ１９が指し示す円の大きさおよび位置とを合わせる。制御部１１は、マウスポインタ１９が指し示す円の中心点の座標および半径を算出し、基準輪郭６２を示す円の中心点の座標（ｘｂ、ｙｂ）と半径「ｒｂ」としてＲＡＭ等に記憶する。尚、添え字の「ｂ」は、「ｂａｓｅ」の頭文字である。
【００４６】
図９は、基準輪郭６２を指定する処理の第２の例を説明する図である。図９に示す例では、予め保持された輪郭検出済のテンプレート画像９１と基準画像６１とのテンプレートマッチング処理を行うことで推定する受精卵の輪郭を基準輪郭６２として指定する。
テンプレート画像９１は、記憶部１２に予め記憶しておいても良いし、必要に応じてメディア入出力部１３、通信制御部１４、周辺機器Ｉ／Ｆ部１７等を介して入力するようにしても良い。また、処理に使用するテンプレート画像９１は、予め決められていても良いし、制御部１１が所定の条件（受精卵の種類、培養環境などの情報に基づく条件）に基づいて複数の画像の中から決定しても良いし、ユーザが選択するようにしても良い。
【００４７】
図９に示すように、テンプレート画像９１は、正確な輪郭であるテンプレート輪郭９２が重畳されている。テンプレート輪郭９２を示す円は、中心点が（ｘｔ、ｙｔ）、半径がｒｔである。尚、添え字の「ｔ」は、「ｔｅｍｐｌａｔｅ」の頭文字である。
テンプレート画像９１と基準画像６１とのテンプレートマッチング処理は、公知技術の中から適切な手法を用いれば良い。制御部１１は、例えば、テンプレート輪郭９２の内部の全画素、テンプレート輪郭９２の内部の一部の画素、テンプレート輪郭９２の周辺の画素などを比較対象とすることができる。
【００４８】
図９に示す処理によって基準輪郭６２を指定する場合、基準画像６１としては、培養初期に撮像された画像を選択することが望ましい。培養初期であれば、受精卵の大きさ、受精卵中心部の状態などに関する個体差が小さく、受精卵の種類や培養環境が同じテンプレート画像９１を選択することで、基準輪郭６２を精度良く検出することができる。
【００４９】
図７の説明に戻る。
次に、制御部１１は、Ｓ１０２において指定した基準輪郭６２の内側と外側に複数の分割領域を定義する（Ｓ１０３）。分割領域は、基準輪郭６２の内側と外側のそれぞれに基準輪郭６２に沿う方向に複数定義され、かつ基準輪郭６２と垂直の方向にそれぞれ一定の幅を有するものである。
【００５０】
図１０は、分割領域６３を説明する図である。
図１０に例示するように、制御部１１は、基準輪郭６２の内側に、基準輪郭６２と垂直の方向に一定の幅（図１０では６４ａにて図示）を有する分割領域６３を複数定義する。言い換えると、まず、中心座標が基準輪郭６２と同一、半径が基準輪郭６２の半径「ｒｂ」から一定の値「ｗ」を引いた値を有する円６５ａ（図１０では点線にて図示）を定義する。そして、基準輪郭６２を外円、円６５ａを内円とするドーナツ形状の二重円領域を複数に分割し、分割領域６３とする。
また、図１０に例示するように、制御部１１は、基準輪郭６２の外側に、基準輪郭６２と垂直の方向に一定の幅（図１０では６４ｂにて図示）を有する分割領域６３を複数定義する。言い換えると、まず、中心座標が基準輪郭６２と同一、半径が基準輪郭６２の半径「ｒｂ」に一定の値「ｗ」を足した値を有する円６５ｂ（図１０では点線にて図示）を定義する。そして、基準輪郭６２を内円、円６５ｂを外円とするドーナツ形状の二重円領域を複数に分割し、分割領域６３とする。
尚、一定の値「ｗ」は、内側と外側ともに同じ値としても良いし、異なる値としても良い。
【００５１】
図１０に示す例では、制御部１１は、基準輪郭６２の中心座標（ｘｂ、ｙｂ）を通る８本の分割線を定義し、内側と外側ともに、８等分に分割している。図１０に示す「Ｒｉ０」〜「Ｒｉ７」が、内側の分割領域６３ａである。また、図１０に示す「Ｒｏ０」〜「Ｒｏ７」が、内側の分割領域６３ｂである。尚、添え字の「ｉ」は、「ｉｎ」の頭文字、「ｏ」は「ｏｕｔ」である。
【００５２】
尚、８等分だけでなく、ｎ等分（ｎは２以上の自然数）に一般化できることは言うまでもない。
また、均等に分割するだけでなく、何らかの基準に基づいて、不均等に分割するようにしても良い。図１、図２を参照すると分かるように、透明帯の領域は、円周方向に沿って一定の範囲で同程度の明るさが続き、ある境界で異なる明るさとなる。このことに着目し、制御部１１は、円周方向に画素を走査していき、一定の大きさのトーンジャンプが発生した箇所を分割領域６３の境界とするようにしても良い。
【００５３】
一定の値「ｗ」は、例えば、受精卵中心部を覆う透明帯の厚さの半分程度の値とすることが望ましい。透明帯の厚さの半分程度とする、すなわち、内側の分割領域６３ａが受精卵中心部を含まないようにすることで、後述する輪郭決定処理において、受精卵中心部の脂肪小滴が自律的に移動することによる影響をなくすことができる。透明帯の厚さは、例えば、受精卵の種類ごとに予め実測しておき、ユーザが入力部１５を介して入力しても良い。
また、一定の値「ｗ」（分割領域６３の一定の幅）は、制御部１１が、基準画像６１に基づいて自動的に決定しても良い。
【００５４】
図１１は、一定の幅を決定する処理を説明する図である。
基準画像６１に対しては、Ｓ１０２において基準輪郭６２が決定されている。図１１に示す「ｒｂ」が、基準輪郭６２の半径である。
制御部１１は、閾値処理（２値化処理）によって脂肪小滴の画素を特定する。更に、制御部１１は、例えば、脂肪小滴の画素を全て含む円の中で最も小さい半径を有する円を算出し、算出した円の内部（図１１では縦線のハッチングによって図示）を推定受精卵中心部６７とする。図１１に示す「ｒｐ」が、推定受精卵中心部６７の半径である。尚、添え字の「ｐ」は、「ｐｒｅｓｕｍｐｔｉｏｎ」の頭文字である。
そして、制御部１１は、基準輪郭６２の半径「ｒｂ」（基準輪郭６２に基づく受精卵の大きさ）から、推定受精卵中心部６７の半径「ｒｐ」（推定した受精卵中心部の大きさ）を引いた値を「ｔ」とする。
「ｔ」は、透明帯の厚さ（ｔｈｉｃｋｎｅｓｓ）と推定できる為、制御部１１は、「ｔ／２」を一定の値「ｗ」として決定する。
【００５５】
図７の説明に戻る。
次に、制御部１１は、分割領域６３内の画素値の統計量を算出し（Ｓ１０４）、算出結果を基準プロファイル６として記憶部１２またはＲＡＭに記憶する（Ｓ１０５）。
分割領域６３内の画素値の統計量は、例えば、平均値とすることができる。図１０に示す例では、分割領域Ｒｉ０〜Ｒｉ７、Ｒｏ０〜Ｒｏ７内の画素値に対して、それぞれ平均値Ａｉ０〜Ａｉ７、Ａｏ０〜Ａｏ７を算出し、ベクトルＶ＝（Ａｉ０、・・・、Ａｉ７、Ａｏ０、・・・、Ａｏ７）とする。基準プロファイル６には、ベクトルＶ、基準輪郭６２の中心座標（ｘｂ、ｙｂ）、半径「ｒｂ」などが含まれる。
【００５６】
分割領域６３内の画素値の統計量として、平均値の他に有用な例について説明する。
図１０に示す例であれば、分割領域Ｒｉ０〜Ｒｉ７、Ｒｏ０〜Ｒｏ７内の画素値に対して、それぞれ平均値Ａｉ０〜Ａｉ７、Ａｏ０〜Ａｏ７を算出する。次に、Ａａｌｌ＝（Ａｉ０＋・・・＋Ａｉ７＋Ａｏ０＋・・・＋Ａｏ７）／１６（平均値の平均）を算出する。そして、Ｂｉ０＝Ａｉ０−Ａａｌｌ、・・・、Ｂｉ７＝Ａｉ７−Ａａｌｌ、Ｂｏ０＝Ａｏ０−Ａａｌｌ、・・・、Ｂｏ７＝Ａｏ７−Ａａｌｌを算出し、ベクトルＶ＝（Ｂｉ０、・・・、Ｂｉ７、Ｂｏ０、・・・、Ｂｏ７）とする。これによって、撮像環境の変動によって基準画像６１とその他の入力画像５との間で画像全体の明るさが変動している場合にも、図１９にて詳述する輪郭決定処理において、適切にプロファイルの比較をすることができる。
【００５７】
また、ベクトルＶｂの要素は、１つの統計量だけでなく、複数の統計量を含めても良い。
図１０に示す例であれば、分割領域Ｒｉ０〜Ｒｉ７、Ｒｏ０〜Ｒｏ７内の画素値に対して、それぞれ平均値Ａｉ０〜Ａｉ７、Ａｏ０〜Ａｏ７、および分散Ｖｉ０〜Ｖｉ７、Ｖｏ０〜Ｖｏ７を算出し、ベクトルＶ＝（Ａｉ０、・・・、Ａｉ７、Ａｏ０、・・・、Ａｏ７、Ｖｉ０、・・・、Ｖｉ７、Ｖｏ０、・・・、Ｖｏ７）とする。これによって、図１９にて詳述する輪郭決定処理において、より詳細にプロファイルの比較をすることができる。
【００５８】
図１２は、基準プロファイル６の例を示す図である。図１２には、視覚的に判別し易いように、図１に示す画像に対して基準輪郭６２が重畳された画像、および各分割領域Ｒｉ０〜Ｒｉ７、Ｒｏ０〜Ｒｏ７内の画素値の平均値に対応する明るさを有する矩形を示している。
【００５９】
次に、図１３から図１８を参照しながら、候補プロファイル作成処理について説明する。
図１３は、候補プロファイル作成処理を示すフローチャートである。
図１３に示すように、画像処理装置１の制御部１１は、入力画像５の中から任意の画像（受精卵の輪郭が不明な画像）を処理画像とし、基準輪郭６２に基づいて処理画像における受精卵の輪郭の候補の一つである初期候補輪郭（第１の候補輪郭）を算出する（Ｓ２０１）。
【００６０】
図１４は、初期候補輪郭７２を決定する処理の第１の例を説明する図である。図１４に示す例では、初期候補輪郭７２の大きさを基準輪郭６２の大きさと同一と仮定し、初期候補輪郭７２を移動させながら、輪郭の内部の画素値に基づいて基準輪郭６２との比較処理を行うことで初期候補輪郭７２の最適な位置を算出する。
【００６１】
基準画像６１に対しては、基準輪郭６２の中心座標（ｘｂ、ｙｂ）および半径「ｒｂ」が決定されている。制御部１１は、基準輪郭６２の内部を比較領域６８ａとする。また、制御部１１は、初期候補輪郭７２の半径「ｒ０」を基準輪郭６２の半径「ｒｂ」とし、仮の中心座標（ｘ０、ｙ０）を定義する。制御部１１は、初期候補輪郭７２の内部を比較領域７８ａとする。
そして、制御部１１は、初期候補輪郭７２の中心座標（ｘ０、ｙ０）をずらしながら、比較領域６８ａと比較領域７８ａとの画素ごとの二乗誤差の総和を算出していき、二乗誤差の総和が最小となる中心座標（ｘ０、ｙ０）を決定する。
【００６２】
図１５は、初期候補輪郭７２を決定する処理の第２の例を説明する図である。図１５に示す例では、初期候補輪郭７２の大きさを基準輪郭６２の大きさと同一と仮定し、初期候補輪郭７２を移動させながら、輪郭の内側と外側を含み、かつ輪郭と垂直の方向にそれぞれ一定の幅を有する領域の画素値に基づいて基準輪郭６２との比較処理を行うことで初期候補輪郭７２の最適な位置を算出する。
【００６３】
基準画像６１は、基準輪郭６２の中心座標（ｘｂ、ｙｂ）および半径「ｒｂ」を決定済である。制御部１１は、基準輪郭６２の内側と外側を含み、かつ輪郭と垂直の方向にそれぞれ一定の幅を有する領域を比較領域６８ｂとする。また、制御部１１は、初期候補輪郭７２の半径「ｒ０」を基準輪郭６２の半径「ｒｂ」とし、仮の中心座標（ｘ０、ｙ０）を定義する。制御部１１は、初期候補輪郭７２の内側と外側を含み、かつ輪郭と垂直の方向にそれぞれ一定の幅を有する領域を比較領域７８ｂとする。
そして、制御部１１は、初期候補輪郭７２の中心座標（ｘ０、ｙ０）をずらしながら、比較領域６８ｂと比較領域７８ｂとの画素ごとの二乗誤差の総和を算出していき、二乗誤差の総和が最小となる中心座標（ｘ０、ｙ０）を決定する。
一定の幅は、比較領域６８ｂ、７８ｂが、受精卵中心部の領域を含まない程度（例えば、図７のＳ１０３にて前述した一定の値「ｗ」）とする。これによって、受精卵中心部の脂肪小滴が自律的に移動することによる影響をなくすことができる。
【００６４】
図１３の説明に戻る。
次に、制御部１１は、初期候補輪郭７２の内側と外側に複数の分割領域を定義する（Ｓ２０２）。分割領域は、初期候補輪郭７２の内側と外側のそれぞれに初期候補輪郭７２に沿う方向に複数定義され、かつ初期候補輪郭７２と垂直の方向にそれぞれ一定の幅を有するものである。
処理画像における分割領域の定義（Ｓ２０２）は、前述した基準画像６１における分割領域の定義（図７のＳ１０３）と同様に行う。
【００６５】
次に、制御部１１は、分割領域内の画素値の統計量を算出し（Ｓ２０３）、算出結果を候補プロファイル７として記憶部１２またはＲＡＭに記憶する（Ｓ２０４）。処理画像における分割領域内の画素値の統計量の算出（Ｓ２０３）は、前述した基準画像６１における分割領域内の画素値の統計量の算出（図７のＳ１０４）と同様に行う。
【００６６】
次に、制御部１１は、候補プロファイルの全パターンを算出済かどうか確認する（Ｓ２０５）。候補プロファイルのパターンについては、Ｓ２０６にて後述する。
全パターンを算出済でなければ（Ｓ２０５のＮｏ）、制御部１１はＳ２０６に進む。
全パターンを算出済であれば（Ｓ２０５のＹｅｓ）、制御部１１は処理を終了する。
【００６７】
Ｓ２０６において、制御部１１は、候補輪郭の変更を行う。候補輪郭の変更とは、２０１にて算出した初期候補輪郭７２の位置および／または大きさを変化させて初期候補輪郭７２とは異なる候補輪郭を算出することである。
具体的には、制御部１１は、初期候補輪郭７２の中心座標のｘ座標「ｘ０」、ｙ座標「ｙ０」、半径「ｒ０」の３つについて、例えば、±１だけ変化させ、全ての組み合わせを候補プロファイルのパターンとして算出する。±１だけ変化させる場合には、３×３×３＝２７通りの候補プロファイルのパターンが算出される。尚、「２７通り」は、初期候補領域７２を含むパターン数である。
【００６８】
図１６〜図１８は、それぞれ候補プロファイル７の第１の例、第２の例、第３の例を示す図である。
図１６と図１７に示す画像に重畳された候補輪郭を比較すると、中心座標および半径が異なる。また、図１６と図１８に示す画像に重畳された候補輪郭を比較すると、少なくとも中心座標が異なる。
そして、中心座標や半径が異なることから、図１６〜図１８に示す各分割領域Ｒｉ０〜Ｒｉ７、Ｒｏ０〜Ｒｏ内の画素値の平均値も、それぞれ異なる明るさを示していることが分かる。
【００６９】
前述した例では、初期候補領域７２の中心座標のｘ座標「ｘ０」、ｙ座標「ｙ０」、半径「ｒ０」の３つについて、±１だけ変化させるとしたが、変化させる幅を広げて、例えば、±２などとしても良い。但し、変化させる幅を広げ過ぎると候補プロファイルのパターンが膨大になってしまい、処理時間が遅くなる。
そこで、適切な基準に基づいて、候補プロファイルのパターン数を増やすことが望ましい。例えば、培養期間の後半になると、受精卵の大きさの変化が激しくなる為、処理画像の撮像時刻と培養開始時刻との差が大きくなるにつれて、候補プロファイルのパターン数を増やすことが望ましい。
【００７０】
図１３の説明に戻る。
Ｓ２０６の後は、Ｓ２０２から処理を繰り返す。制御部１１は、変更した候補輪郭の内側と外側に複数の分割領域を定義する（Ｓ２０２）。分割領域は、変更した候補輪郭の内側と外側のそれぞれに変更した候補輪郭に沿う方向に複数定義され、かつ変更した候補輪郭と垂直の方向にそれぞれ一定の幅を有するものである。次に、制御部１１は、分割領域内の画素値の統計量を算出し（Ｓ２０３）、算出結果を候補プロファイル７として記憶部１２またはＲＡＭに記憶する（Ｓ２０４）。
以下、Ｓ２０５において候補プロファイルの全パターンを算出済と判定するまで、Ｓ２０６、Ｓ２０２〜Ｓ２０４を繰り返す。
【００７１】
次に、図１９を参照しながら、輪郭決定処理について説明する。
図１９は、輪郭決定処理を示すフローチャートである。
図１９に示すように、画像処理装置１の制御部１１は、候補プロファイル７を１つ選択し（Ｓ３０１）、基準プロファイル６と候補プロファイル７との類似度を算出する（Ｓ３０２）。
【００７２】
類似度は、例えば、基準プロファイル６と候補プロファイル７の両方に含まれるベクトルの各要素（分割領域６３ごとの統計量）同士の二乗誤差の総和Σ１に基づいて定義する。二乗誤差を求めるペアは、同じ位置の分割領域６３同士とする。
この場合、例えば、基準プロファイル６のベクトルがＶｂ＝（Ｂi０、・・・Ｂiｎ、Ｂｏ０、・・・Ｂｏｎ）、候補プロファイル７のベクトルがＶｃ＝（Ｃi０、・・・、Ｃiｎ、Ｃｏ０、・・・、Ｃｏｎ）（ｎは２以上の自然数）とすると、Σ１＝（Ｂi０−Ｃi０）＾２＋・・・＋（Ｂiｎ−Ｃiｎ）＾２＋（Ｂｏ０−Ｃｏ０）＾２＋・・・＋（Ｂｏｎ−Ｃｏｎ）＾２である。
【００７３】
また、その他に有用な類似度としては、輪郭を介して対向する内側の分割領域６３ａの統計量と外側の分割領域６３ｂの統計量の差に関する二乗誤差の総和Σ２に基づいて定義するものがある。
この場合、例えば、基準プロファイル６のベクトルがＶｂ＝（Ｂi０、・・・Ｂiｎ、Ｂｏ０、・・・Ｂｏｎ）、候補プロファイル７のベクトルがＶｃ＝（Ｃi０、・・・、Ｃiｎ、Ｃｏ０、・・・、Ｃｏｎ）（ｎは２以上の自然数）とすると、Σ２＝（｜Ｂi０−Ｂｏ０｜−｜Ｃi０−Ｃｏ０｜））＾２＋・・・＋（｜Ｂiｎ−Ｂｏｎ｜−｜Ｃiｎ−Ｃｏｎ｜））＾２である（絶対値の記号は外しても良い）。
類似度としてΣ２を用いることは、輪郭の内側と外側のコントラストの強弱を用いて輪郭を特徴付けることを意味する。
【００７４】
次に、制御部１１は、全ての候補プロファイル７を処理済かどうか確認する（Ｓ３０３）。
全ての候補プロファイル７を処理済ではない場合（Ｓ３０３のＮｏ）、制御部１１は、Ｓ３０１から処理を繰り返す。
全ての候補プロファイル７を処理済の場合（Ｓ３０３のＹｅｓ）、制御部１１は、Ｓ３０４に進む。
【００７５】
次に、制御部１１は、Ｓ３０２の算出結果から、基準プロファイル６と最も類似する候補プロファイル７、すなわち最も類似度が高い（前述のΣ１またはΣ２の値が最も小さい）候補プロファイル７に基づいて、図１３のＳ２０１にて選択した処理画像の輪郭を決定する（Ｓ３０４）。処理画像は、入力画像５の中から選択した受精卵の輪郭が不明な画像である。
【００７６】
そして、制御部１１は、決定した輪郭に基づいて出力画像８を生成する（Ｓ３０５）。
例えば、制御部１１は、同一の受精卵に対する全ての撮影時刻の入力画像５に対して輪郭を決定した後、どの撮影時刻の輪郭（すなわち、受精卵を覆うドーナツ状の透明帯の外縁）についても多少の余裕を持って囲むことができる全撮影時刻共通の大きさの矩形を算出する。そして、制御部１１は、各撮影時刻の入力画像５において矩形の中心と受精卵の中心を一致させて、矩形に沿って入力画像５を切り抜くことで出力画像８を生成し、記憶部１２に記憶する。また、制御部１１は、決定した輪郭の情報（半径、中心座標などの情報）については、出力画像８と組にして記憶部１２に記憶する。これによって、受精卵周辺の部位の画像情報が損失せず、管理し易い出力画像８（受精卵観測画像）を得ることができる。
【００７７】
また、例えば、制御部１１は、処理画像に決定した輪郭を重畳し、出力画像８を生成しても良い。
また、例えば、制御部１１は、処理画像から決定した受精卵の輪郭の内部を切り抜くことで出力画像８を生成しても良い。
その他、制御部１１は、解析内容に応じて、決定した輪郭に基づいて様々な方法で出力画像８を生成することができる。
【００７８】
画像処理装置１は、受精卵の輪郭が不明な全ての画像（基準画像６を除く入力画像５）を処理画像とし、以上の処理を行うことで、同一の受精卵を異なる時刻で光源方向を一定にして撮像することによって得られた入力画像５の全ての画像に対して、受精卵の位置を正確に合わせた出力画像８を生成することができる。
尚、基準画像６は、固定する必要はなく、培養開始時刻からの時間の経過につれて、新たな画像を基準画像６として指定するようにしても良い。この場合、図１９に示す処理によって決定済の輪郭を有する出力画像８を基準画像６とすることもできる。
【００７９】
前述の図５から図１９の説明では、受精卵および受精卵中心部は、球体と仮定して説明したが、楕円体と仮定して、本発明を適用することも可能である。
受精卵の回転を無視する場合には、前述した円における半径のパラメータを、楕円の長径と短径の２つのパラメータに置き換えることで、円に対する処理とほぼ同様に、前述した処理を行うことができる。
受精卵を楕円体と仮定する場合、候補プロファイル８のパターンを生成するパラメータは、中心座標のｘ座標、ｙ座標、長径、短径の４つとなる。
【００８０】
また、受精卵の回転を考慮する場合についても、本発明の技術的思想を拡張することが可能である。受精卵が回転角度θだけ回転している場合、分割領域６３を定義する際の角度の計測開始点（図１０の例であれば、分割領域「Ｒｉ０」、「Ｒｏ０」の円周方向の端部）を回転角度θに従って変更する。
受精卵を楕円体と仮定し、かつ受精卵の回転を考慮する場合、候補プロファイル８のパターンを生成するパラメータは、中心座標のｘ座標、ｙ座標、長径、短径、回転角度の５つとなる。
【００８１】
尚、回転に起因して輪郭周辺の統計量に変化がある場合には、回転による統計量の変化をモデル化する必要がある。すなわち、基準プロファイル６に対して回転角度θを適用して、候補プロファイル８と比較する必要がある。
モデル化の方法の一例としては、回転角度θ＝０、ｘ、２ｘ、３ｘ、・・・（ｘは、例えば３６０の約数）のように所定の間隔を設けて、回転角度ごとに楕円の中心座標、長径、短径、分割領域６３の統計量を算出し、複数の基準プロファイル６を生成し、それらを線形補間する方法が挙げられる。
これによって、回転に起因する輪郭周辺の統計量の変化の影響を抑えることができる。
【００８２】
更に、円や楕円における中心座標に代えて輪郭内部の領域の重心座標とし、円や楕円における半径、長径、短径に代えて輪郭内部の領域の面積とすることで、一般形状に対して本発明の技術的思想を拡張することも可能である。回転を考慮する場合については、前述の楕円の場合と同様である。
【００８３】
以上、添付図面を参照しながら、本発明に係る画像処理装置等の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【００８４】
１………画像処理装置
２………基準プロファイル作成プログラム
３………候補プロファイル作成プログラム
４………輪郭決定プログラム
５………入力画像
６………基準プロファイル
７………候補プロファイル
８………出力画像
６１………基準画像
６２………基準輪郭
６３………分割領域
７１………処理画像
７２………初期候補輪郭
９１………テンプレート画像
９２………テンプレート輪郭
【図５】
0005418324.tif 000002
【図６】
0005418324.tif 000003
【図７】
0005418324.tif 000004
【図８】
0005418324.tif 000005
【図９】
0005418324.tif 000006
【図１０】
0005418324.tif 000007
【図１１】
0005418324.tif 000008
【図１３】
0005418324.tif 000009
【図１４】
0005418324.tif 000010
【図１５】
0005418324.tif 000011
【図１９】
0005418324.tif 000012
【図１】
0005418324.tif 000013
【図２】
0005418324.tif 000014
【図３】
0005418324.tif 000015
【図４】
0005418324.tif 000016
【図１２】
0005418324.tif 000017
【図１６】
0005418324.tif 000018
【図１７】
0005418324.tif 000019
【図１８】
0005418324.tif 000020