(57)【要約】【課題】受精卵の撮影画像を経時的に解析し、前核の状態を解析する手法を提供する。【解決手段】画像解析装置は、受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出し、次に前記候補領域内に存在する前核候補の数である前核候補数を検出する。そして、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する。【選択図】図３
2018097796.tif 000002
【特許請求の範囲】
【請求項１】
受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出手段と、
前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出手段と、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定手段と、
を備えることを特徴とする画像解析装置。
【請求項２】
前記判定手段は、
前記複数の撮影画像毎に前記候補領域のうちの１つを前核領域と判定する第１判定手段と、
当該前核領域内に存在する前核候補数の経時変化に基づいて前記前核の数を判定する第２判定手段と、
を備えることを特徴とする請求項１に記載の画像解析装置。
【請求項３】
前記第１判定手段は、前記候補領域のうち、経時変化において前核候補数が２つになった期間を有する候補領域を前記前核領域と判定することを特徴とする請求項２に記載の画像解析装置。
【請求項４】
前記判定手段は、
各候補領域内に存在する前核候補数の経時変化に基づいて、前記候補領域のうちの１つを前核領域と判定する第１判定手段と、
当該前核領域内における前核候補数の経時変化に基づいて前記前核の数を判定する第２判定手段と、
を備えることを特徴とする請求項１に記載の画像解析装置。
【請求項５】
前記第１判定手段は、経時変化において前記前核候補領域内における前核候補数が２つになった期間を有する場合、当該前核候補領域を前記前核領域と判定することを特徴とする請求項４に記載の画像解析装置。
【請求項６】
前記第１判定手段は、各候補領域において、異なる前核候補数毎の尤度を算出し、尤度が最も高い前核候補数を当該候補領域の前核候補数とすることを特徴とする請求項４に記載の画像解析装置。
【請求項７】
画像解析装置により実行される画像解析方法であって、
受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出工程と、
前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出工程と、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程と、
を備えることを特徴とする受精卵の画像解析方法。
【請求項８】
コンピュータを備える画像解析装置により実行されるプログラムであって、
受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出工程、
前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出工程、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程、
を前記コンピュータに実行させることを特徴とするプログラム。
【請求項９】
受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得手段と、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定手段と、
を備えることを特徴とする画像解析装置。
【請求項１０】
画像解析装置により実行される画像解析方法であって、
受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得工程と、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程と、
を備えることを特徴とする受精卵の画像解析方法。
【請求項１１】
コンピュータを備える画像解析装置により実行されるプログラムであって、
受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得工程、
前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程、
を前記コンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、受精卵における前核の状態を判定する技術に関する。
【背景技術】
【０００２】
不妊治療分野では、治療率や患者のＱＯＬ（Ｑｕａｌｉｔｙ  Ｏｆ  Ｌｉｆｅ）の観点から、胚の着床率を予測し、移植対象の胚を選別することが行われている。従来、胚を非侵襲的に評価する手法として、顕微鏡による形態観察と、その際の形態から発育状況を判定することが行われてきた。
【０００３】
近年、タイムラプス装置の普及により、上記に加えて経時的な変化をも捉えることが可能となり、ある工程中・ある工程間の形態の時間的長さを指標とする試みもなされてきた。このような発育予測指標としては、これまで、胚の形態、卵割球の数や断片化の程度、これらの時間的関係性などが報告されている。また、形態以外の指標としては、胚の呼吸量（酸素消費量）や代謝物濃度などが知られている。この種の手法は、例えば特許文献１、２などに記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第５４８１６９６号公報
【特許文献２】特許第５８０７２８８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来は、上記のような様々な指標が提案されてきたが、胚の発育の最も初期のイベントである受精の確認については、確立された方法がない状況であった。本発明は、受精卵の撮影画像を経時的に解析し、前核の状態を解析する手法を提供することを主な目的とする。
【課題を解決するための手段】
【０００６】
本発明の１つの観点では、画像解析装置は、受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出手段と、前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出手段と、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定手段と、を備える。
【０００７】
上記の画像解析装置は、受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出し、次に候補領域内に存在する前核候補の数である前核候補数を検出する。そして、前核候補数の経時変化に基づいて、受精卵における前核の数を判定する。こうして、時系列に撮影した複数の撮影画像の画像解析により、受精卵における前核の数を判定することができる。
【０００８】
上記の画像解析装置の一態様では、前記判定手段は、前記複数の撮影画像毎に前記候補領域のうちの１つを前核領域と判定する第１判定手段と、当該前核領域内に存在する前核候補数の経時変化に基づいて前記前核の数を判定する第２判定手段と、を備える。この態様では、撮影画像毎に前核領域を決定し、その前核領域について前核数の判定を行う。
【０００９】
好適には、前記第１判定手段は、前記候補領域のうち、経時変化において前核候補数が２つになった期間を有する候補領域を前記前核領域と判定する。
【００１０】
上記の画像解析装置の他の一態様では、前記判定手段は、各候補領域内に存在する前核候補数の経時変化に基づいて、前記候補領域のうちの１つを前核領域と判定する第１判定手段と、当該前核領域内における前核候補数の経時変化に基づいて前記前核の数を判定する第２判定手段と、を備える。この態様では、候補領域の１つを前核領域と決定し、その前核領域内における前核候補数の経時変化を見て前核の数を判定する。
【００１１】
好適には、前記第１判定手段は、経時変化において前記前核候補領域内における前核候補数が２つになった期間を有する場合、当該前核候補領域を前記前核領域と判定する。また、好適には、前記第１判定手段は、各候補領域において、異なる前核候補数毎の尤度を算出し、尤度が最も高い前核候補数を当該候補領域の前核候補数とする。
【００１２】
本発明の他の観点では、画像解析装置により実行される画像解析方法は、受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出工程と、前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出工程と、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程と、を備える。この方法によっても、時系列に撮影した複数の撮影画像の画像解析により、受精卵における前核の数を判定することができる。
【００１３】
本発明のさらに他の観点では、コンピュータを備える画像解析装置により実行されるプログラムは、受精卵を時系列に撮影した複数の撮影画像から、前核の候補領域を検出する第１検出工程、前記候補領域内に存在する前核候補の数である前核候補数を検出する第２検出工程、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程、を前記コンピュータに実行させる。このプログラムをコンピュータで実行することにより、上記の画像解析装置を実現することができる。
【００１４】
本発明のさらに他の観点では、画像解析装置は、受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得手段と、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定手段と、を備える。この画像解析装置によれば、時系列に撮影した複数の撮影画像の画像解析により、受精卵における前核の数を判定することができる。
【００１５】
本発明のさらに他の観点では、画像解析装置により実行される画像解析方法は、受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得工程と、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程と、を備える。この画像解析方法によれば、時系列に撮影した複数の撮影画像の画像解析により、受精卵における前核の数を判定することができる。
【００１６】
本発明のさらに他の観点では、コンピュータを備える画像解析装置により実行されるプログラムは、受精卵を時系列に撮影した複数の撮影画像から検出された前核の候補領域内に存在する前核候補の数である前核候補数を取得する取得工程、前記前核候補数の経時変化に基づいて、前記受精卵における前核の数を判定する判定工程、を前記コンピュータに実行させる。このプログラムをコンピュータで実行することにより、上記の画像解析装置を実現することができる。
【図面の簡単な説明】
【００１７】
【図１】実施形態に係る受精卵の画像解析システムの構成を示す。
【図２】受精卵の撮影画像の例を示す。
【図３】第１実施形態による前核状態解析処理のフローチャートである。
【図４】第１実施形態による前核の解析方法を説明する図である。
【図５】第２実施形態による前核状態解析処理のフローチャートである。
【図６】第２実施形態による前核の解析方法を説明する図である。
【発明を実施するための形態】
【００１８】
以下、本発明の好適な実施形態について説明する。
［背景］
これまでは、受精後一定の時刻に前核の数を観察し、その数によって受精結果を判定してきていた。通常、前核の数は２個（２ＰＮ）であり、前核の数が０，１，３，４以上のような数の場合は異常と見なされてきたが、２個以外の数的異常が観察手法によるアーティファクトなのか、実際上の異常なのかの切り分けが難しかった。
【００１９】
しかし、最近、タイムラプス装置を用いて前核の出現パターンを詳細に調べた結果が報告されている。その報告によると、前核が０個や１個の場合には異常である場合と正常だが観察漏れしている場合の２種類が混在していること、前核が３個以上の場合には正常胚はなく、異常と見なして良いこと、などがわかってきた。これにより、受精のイベントを正確に診断へ活かすための方法として、前核の数を判定すれば良いことがわかってきた。そこで、本実施形態では、受精卵の時系列的な撮影画像を解析することにより、前核の数を判定する。
【００２０】
［システム構成］
図１は、実施形態に係る受精卵の画像解析システムの構成を示す。図示のように、画像解析システムは、大別して撮影装置２と、画像解析装置３とを備える。撮影装置２は、例えば電子顕微鏡などであり、培養容器１０内の受精卵の画像を撮影し、撮影画像を画像解析装置３へ出力する。具体的には、撮影装置２は、受精後一定の期間において、同一の受精卵について所定時間間隔で時系列的に複数の撮影画像を生成し、画像解析装置３に出力する。
【００２１】
画像解析装置３は、ＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）などにより構成され、同一の受精卵についての異なる時刻の撮影画像を解析することにより、受精卵に存在する前核の数を判定する。培養容器１０内には複数の受精卵が載置されており、撮影装置２は必要に応じて、受精卵を識別するための識別情報、例えば、受精卵が収容されているウェルの識別情報などと対応付けて撮影画像を画像解析装置３へ出力する。
【００２２】
［撮影画像］
図２は、培養容器１０及び受精卵の撮影画像の例を示す。培養容器１０は円形の容器であり、中央部に収容部１１が形成されている。収容部１１内には、複数（図２の例では５×５＝２５個）のウェル１２が形成されている。ウェル１２は受精卵を収容するための窪みであり、１つのウェル１２に１つの受精卵が収容されている。図１に示すように、収容部１１は培養液１３で満たされている。
【００２３】
図２には、１つのウェル１２の部分の撮影画像２０が示されている、撮影画像２０は、ウェル１２と、ウェル１２内にある受精卵３０とを含んでいる。撮影装置２は、同一の受精卵（即ち、同一のウェル１２）についての撮影画像を所定時間間隔で複数生成し、画像解析装置３へ出力する。なお、撮影装置２は、ウェル１２毎に１枚の画像を撮影してもよく、多数のウェル１２を含む画像を撮影した後、それをウェル１２毎に分割して個々のウェル１２に対応する撮影画像としても良い。
【００２４】
［前核状態解析処理］
（第１実施形態）
図３は、第１実施形態に係る前核状態解析処理のフローチャートである。この処理は、撮影装置２から同一の受精卵についての複数の撮影画像を受信した画像解析装置３により実行される。まず、画像解析装置３は、受精卵の撮影画像の画像解析により、ウェル１２の領域（以下、「ウェル領域」と呼ぶ。）を検出する（ステップＳ１１）。画像解析装置３は、例えばエッジ抽出処理によりウェル領域を検出する。なお、ウェル１２のサイズは既知であるので、画像解析装置３はウェル領域を検出する際に、ウェル１２のサイズも考慮することができる。
【００２５】
次に、画像解析装置３は、撮影画像に基づいて、受精卵３０の領域（以下、「受精卵領域」と呼ぶ。）を検出する（ステップＳ１２）。画像解析装置３は、例えばウェル領域の内側の領域における色の違いに着目して受精卵領域を検出する。通常、ウェル１２の内側の領域において、受精卵以外の領域は白であり、受精卵３０の領域には多少の色がある。よって、画像解析装置３は、色の違いにより受精卵領域を検出することができる。
【００２６】
次に、画像解析装置３は、撮影画像に基づいて、受精卵領域内にある前核候補領域３２を検出する（ステップＳ１３）。ここで、「前核候補領域」とは、前核の可能性があると考えられる領域であり、受精卵領域内において周囲と輝度や色が異なる領域である。画像解析装置３は、例えば受精卵領域における輝度勾配を算出したり、輪郭を抽出したりすることにより前核候補領域を検出する。図２に例示する撮影画像２０では、受精卵３０の内側に４つの前核候補領域３２ａ〜３２ｄが検出されている。
【００２７】
次に、画像解析装置３は、各前核候補領域３２について、前核候補の数を解析する（ステップＳ１４）。ここで、「前核候補」とは、その撮影画像における形状、色などの特徴から、前核らしきもの、前核の可能性があると認められる部分を指し、「前核候補の数」（以下、「前核候補数」とも呼ぶ。）は、前核候補領域３２内に存在する前核候補の数を指す。
【００２８】
図２に示す撮影画像２０における前核候補数を図４（Ａ）に示す。図２に示されるように、前核候補領域３２ａには前核候補を２つ見つけることができ、前核候補領域３２ｂには前核候補を１つ見つけることができ、前核候補領域３２ｃには前核候補を１つ見つけることができ、前核候補領域３２ｄには前核候補を３つ見つけることができる。このように、画像解析装置３は、各前核候補領域３２ａ〜３２ｄについて、そこに含まれる前核候補の数を検出する。
【００２９】
次に、画像解析装置３は、撮影画像毎に前核領域を決定する（ステップS１５）。即ち、画像解析装置３は、撮影画像に含まれる１又は複数の前核候補領域３２から、そのうちの１つを前核領域と決定する。前核候補領域３２内で検出される前核候補は、本当に前核である場合もあるが、実際には前核ではなく、撮影画像上に表れるノイズその他のアーティファクトであることもある。そこで、画像解析装置３は、各前核候補領域３２内の前核候補を解析し、最も確からしい（尤度が高い）前核候補を含む前核候補領域を前核領域と決定する。
【００３０】
具体的には、前述のように正常な受精卵の前核数は２個であるので、基本的には２個の前核候補が検出された前核候補領域３２が前核領域として決定されることになる。一方、いずれの前核候補領域３２においても２個の前核候補が検出されなかった場合には、検出された前核候補数が０個、１個などであっても、その撮影画像において検出された複数の前核候補領域３２のうち最も尤度が高い前核候補領域３２が前核領域として決定される。図４（Ａ）は、図２の例における前核候補領域３２ａ〜３２ｄについての前核領域の判定結果を示す。この例では、４つの前核候補領域３２ａ〜３２ｄのうち、２個の前核候補が検出された前核候補領域３２ａが前核領域と判定されている。こうして、画像解析装置３は、１つ以上の前核候補領域３２が検出された撮影画像については、１つの前核領域を決定する。また、画像解析装置３は、決定された前核領域に含まれる前核候補を前核とみなす。
【００３１】
次に、画像解析装置３は、１つの受精卵について生成された全ての撮影画像の処理が終了したか否かを判定する（ステップＳ１６）。全ての撮影画像の処理が終了していない場合には（ステップＳ１６：Ｎｏ）、画像解析装置３は、ステップＳ１１へ戻って次の撮影画像について上記のステップＳ１１〜Ｓ１５の処理を行う。そして、全ての撮影画像の処理が終了した場合（ステップＳ１６：Ｙｅｓ）、画像解析装置３はステップＳ１７へ進む。
【００３２】
なお、複数の撮影画像についてステップＳ１１〜Ｓ１５の処理を行う際、ある撮影画像についてステップＳ１３で前核候補領域が１つも検出できなかった場合には、ステップＳ１４、Ｓ１５の処理は行わない。即ち、前核候補領域が検出できなかった撮影画像については、前核領域は決定されない。また、ステップＳ１１でウェル領域を検出できない撮影画像やステップＳ１２で受精卵領域を検出できない撮影画像は、処理の対象から除外する。
【００３３】
ステップＳ１７では、画像解析装置３は、各撮影画像について、ステップＳ１５で決定された前核領域内の前核数の経時変化に基づいて、解析の対象となっている受精卵（以下、「対象受精卵」とも呼ぶ。）の最終的な前核数を判定する（ステップＳ１７）。図４（Ｂ）は、各撮影画像における前核領域内の前核数の経時変化を示す。横軸の時刻ｔ１〜ｔ８は、各撮影画像の撮影時刻を示す。即ち、この例では、時刻ｔ１に撮影された撮影画像から時刻ｔ８に撮影された撮影画像までの８枚の撮影画像について、それぞれの前核領域内の前核数が示されている。前核数の経時変化を考えた場合、正常な受精卵であれば、前核数は時間の経過とともに、０個→１個→２個→０個と推移することが知られている。よって、画像解析装置３は、各撮影画像の前核領域内の前核数がこのように推移した場合には、前核数は２個であると判定する。前核数が２個である受精卵は、正常な受精卵であると判定される。
【００３４】
なお、前核数の推移は理論上は上記のようになるが、画像解析装置３は所定時間間隔で撮影した画像に基づいて解析を行うので、例えば前核数が１個である期間において撮影画像が存在しないような場合には、前核数の推移は、０個→２個→０個となることも考えられる。よって、このような場合でも前核数は２個であると判定してよい。但し、いずれの撮影画像においても前核数が２個に達しないような場合には、前核数を２個と判定することはない。
【００３５】
こうして、最終的に対象受精卵の前核数が決定されると、画像解析装置３は解析結果を出力する（ステップＳ１８）。例えば、画像解析装置３は、図示しない表示部などに解析結果を表示する。解析結果は、少なくともステップＳ１７で判定された、対象受精卵の前核数を含む。また、必要に応じて、各撮影画像における前核候補領域の数、前核候補の撮影画像内における位置、前核の撮影画像内における位置、どの前核候補領域が前核領域と判定されたか、どの前核候補が前核と判定されたかなどの情報を解析結果に含めても良い。
【００３６】
（第２実施形態）
第１実施形態では、基本的に撮影画像毎に１つの前核領域を決定し、その前核領域内の前核数の経時変化を解析して対象受精卵の前核数を決定している。これに対し、第２実施形態では、前核候補領域毎に、複数の撮影画像にわたる前核候補数の経時変化を解析して前核数を決定し、その結果に基づいていずれかの前核候補領域を前核領域と決定し、決定された前核領域における前核数を対象受精卵の前核数と決定する。
【００３７】
図５は、第２実施形態に係る前核状態解析処理のフローチャートである。この処理は、撮影装置２から同一の受精卵について複数の撮影画像を受信した画像解析装置３により実行される。ここで、ステップＳ２１〜Ｓ２４は、図３に示す第１実施形態のステップＳ１１〜Ｓ１４と同一であるので、説明を省略する。
【００３８】
次に、画像解析装置３は、１つの受精卵について生成された全ての撮影画像の処理が終了したか否かを判定する（ステップＳ２５）。全ての撮影画像の処理が終了していない場合には（ステップＳ２５：Ｎｏ）、画像解析装置３は、ステップＳ２１へ戻って次の撮影画像について上記のステップＳ２１〜Ｓ２４の処理を行う。そして、全ての撮影画像の処理が終了した場合（ステップＳ２５：Ｙｅｓ）、画像解析装置３はステップＳ２６へ進む。
【００３９】
ステップＳ２６では、画像解析装置３は、複数の撮影画像にわたって、前核候補領域毎に前核候補数の経時変化を解析して前核数を決定する。具体的には、画像解析装置３は、各前核候補領域内で検出された前核候補数、及び、少なくともその前後の前核候補数についての尤度を算出し、最も高い尤度を有する前核候補数をその前核候補領域における前核数と決定する。例えば、ある前核候補領域で１個の前核候補領域が検出された場合、画像解析装置３は、その前核候補領域の画像の特徴などに基づいて、その前核候補領域における前核数が１個である場合の尤度に加えて、前核数が０個である場合の尤度、及び、前核数が２個である場合の尤度を算出する。そして、それらのうち最も高い尤度を有する前核数を、その前核候補領域における前核数と決定する。
【００４０】
このとき、画像解析装置３は、必要に応じて、前核候補数の経時変化に基づいて前核数を決定する。具体的には、画像解析装置３は、時間的に前後する撮影画像における前核候補数を考慮して前核数を決定する。例えば、ある連続する５枚の撮影画像において同一の前核候補領域で検出された前核数が、
０個→１個→１個→２個→１個→２個→２個
と変化した場合、前後の前核数を考慮すると、５番目の撮影画像の前核数は１個ではなく、２個であると判断することもできる。なお、その場合には、５番目の撮影画像における前核数が１個である場合の尤度と、２個である場合の尤度とを比較し、それらの差が所定値より小さい場合に５番目の画像における前核数を１個から２個へと修正するようにしてもよい。
【００４１】
図６は、上記のようにして、複数の撮影画像にわたり前核候補領域毎に決定された前核数の例を示す。なお、この例は、図２に示す撮影画像２０とは異なる撮影画像の例であるものであり、図２に示す撮影画像２０に図示した前核候補数とは一致しない。
【００４２】
図６の例において、前核候補領域３２ａについては、時刻ｔ１３の撮影画像から前核候補領域３２ａが現れ、時刻ｔ２１の画像以降は前核候補領域３２ａが消えている。また、時刻ｔ１３〜ｔ１５の撮影画像では前核候補領域３２ａ内に前核は検出されず、時刻ｔ１６の撮影画像では前核候補領域３２ａ内に１個の前核が検出され、時刻ｔ１７〜ｔ１９の撮影画像では前核候補領域３２ａ内に２個の前核が検出され、時刻ｔ２０の撮影画像では前核候補領域３２ａ内に前核は検出されなくなっている。
【００４３】
前核候補領域３２ｂについては、時刻ｔ１２の撮影画像から前核候補領域３２ｂが現れ、時刻ｔ２３以降は前核候補領域３２ｂが消えている。また、時刻ｔ１２〜ｔ２２の撮影画像において前核候補領域３２ｂ内に１個の前核が検出されている。
【００４４】
前核候補領域３２ｃについては、時刻ｔ１１の撮影画像から前核候補領域３２ｃが現れ、時刻ｔ１７以降は前核候補領域３２ｃが消えている。また、時刻ｔ１１〜ｔ１２の撮影画像においては前核候補領域３２ｃ内に前核は検出されず、時刻ｔ１３〜ｔ１４の撮影画像においては前核候補領域３２ｃ内に１個の前核が検出され、時刻ｔ１５〜ｔ１６の撮影画像においては前核候補領域３２ｃ内に前核が検出されなくなっている。
【００４５】
前核候補領域３２ｄについては、時刻ｔ１９の撮影画像から前核候補領域３２ｄが現れ、時刻ｔ２４以降は前核候補領域３２ｄが消えている。また、時刻ｔ１９〜２１の撮影画像においては前核候補領域３２ｃ内に前核は検出されず、時刻ｔ２２〜ｔ２３の撮影画像においては前核候補領域３２ｄ内に１つの前核が検出されている。
【００４６】
次に、画像解析装置３は、各前核候補領域における前核数の経時変化から、前核領域を判定して前核数を決定する（ステップＳ２７）。具体的には、画像解析装置３は、図６に例示した各前核候補領域における前核数の推移に基づいて、複数の前核候補領域のうちの１つの前核候補領域を前核領域と決定する。前述のように、正常な受精卵であれば、前核数は時間の経過とともに、０個→１個→２個→０個と推移することが知られている。よって、画像解析装置３は、複数の前核候補領域のうち、このような正常な推移に近い経時変化を示している前核候補領域を前核領域であると判定する。図６の例においては、画像解析装置３は、前核候補領域３２ａを前核領域と判定し、対象受精卵の前核数を２個と決定する。
【００４７】
こうして前核領域が判定され、対象受精卵の前核数が決定されると、画像解析装置３は解析結果を出力する（ステップＳ２８）。なお、解析結果の出力方法は、基本的に第１実施形態のステップＳ１８と同様であるので、説明を省略する。
【００４８】
以上のように、本実施形態によれば、受精卵の時系列的な撮影画像を解析することにより前核の数を判定することができ、これにより正常な受精卵を効率的に確認することができる。
【００４９】
［変形例］
上記の実施形態では、受精卵を時系列に撮影した複数の撮影画像から前核の候補領域を検出して前核候補数を検出する処理と、前核候補数の経時変化に基づいて前核の数を判定する処理を同一の画像解析システムにより実施しているが、これらの処理を別々の装置を用いて実施するようにしてもよい。即ち、第１の装置により、受精卵を時系列に撮影した複数の撮影画像から前核の候補領域を検出して前核候補数を検出する処理を実施する。そして、その結果を第２の装置に供給し、第２の装置により、前核候補数の経時変化に基づいて前核の数を判定する処理を実施すればよい。
【符号の説明】
【００５０】
２  撮影装置
３  画像解析装置
１０  培養容器
１２  ウェル
２０  撮影画像
３０  受精卵
３２  前核候補領域
【図１】
2018097796.tif 000003
【図２】
2018097796.tif 000004
【図３】
2018097796.tif 000005
【図４】
2018097796.tif 000006
【図５】
2018097796.tif 000007
【図６】
2018097796.tif 000008