(57)【要約】【課題】主観的判断をすることなく、受精卵の画像データだけで正常児獲得の可否を予測可能な診断結果が得られる受精卵の画像診断システムを提供する。【解決手段】  受精卵の画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、解析手段６は、画像データを年齢別の複数のグループに区分し、複数のグループから一つを選んだ解析対象グループについて、複数の識別器を作成して、テストデータの正診率を演算し、正診率に基づいて複数の識別器から年齢別の最適識別器を決定する。【選択図】図１
2020067680.tif 000002
【特許請求の範囲】
【請求項１】
機械学習の手法を用いた受精卵の画像診断システムであって、
妊娠予後が判明している受精卵の画像データを格納する記憶手段と、
受精卵の画像データを記憶手段に取り込むための受信手段と、
前記画像データを解析する解析手段と、
前記解析手段による解析結果を出力する出力手段とを備えており、
前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、
前記解析手段は、
前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、
前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、
前記学習用データから複数の識別器を作成し、
前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、
前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、
前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする受精卵の画像診断システム。
【請求項２】
前記解析手段は、前記最適識別器を決定する前に、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データの個数を変えて繰り返し、個数の異なる前記学習用データについて演算された前記正診率に基づいて、前記学習用データの最適個数を決定し、前記学習用データの数を前記最適個数として、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データを変えて再度繰り返し、前記再度の繰り返しにより得られた複数の前記正診率に基づいて、前記再度の繰り返しにより得られた複数の前記識別器から最適識別器を決定する請求項１に記載の受精卵の画像診断システム。
【請求項３】
コンピュータに、機械学習の手法を用いて受精卵の画像診断を実行させるための受精卵の画像診断プログラムであって、
受精卵の画像データを記憶手段に取り込む取り込みステップと、
前記画像データを解析する解析ステップと、
前記解析手段による解析結果を出力する出力ステップとを前記コンピュータに実行させ、
前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、
前記解析ステップにおいて、
前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、
前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、
前記学習用データから複数の識別器を作成し、
前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、
前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、
前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする受精卵の画像診断プログラム。
【請求項４】
前記解析ステップにおいて、前記解析手段は、前記最適識別器を決定する前に、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データの個数を変えて繰り返し、個数の異なる前記学習用データについて演算された前記正診率に基づいて、前記学習用データの最適個数を決定し、前記学習用データの数を前記最適個数として、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データを変えて再度繰り返し、前記再度の繰り返しにより得られた複数の前記正診率に基づいて、前記再度の繰り返しにより得られた複数の前記識別器から最適識別器を決定する請求項３に記載の受精卵の画像診断プログラム。
【請求項５】
機械学習の手法を用いた受精卵の画像診断用識別器の作成方法であって、
受精卵の画像データを記憶手段に取り込む取り込みステップと、
前記画像データを解析する解析ステップと、
前記解析手段による解析結果を出力する出力ステップとを備え、
前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、
前記解析ステップにおいて、
前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、
前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、
前記学習用データから複数の識別器を作成し、
前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、
前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、
前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする受精卵の画像診断用識別器の作成方法。
【請求項６】
前記解析ステップにおいて、前記解析手段は、前記最適識別器を決定する前に、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データの個数を変えて繰り返し、個数の異なる前記学習用データについて演算された前記正診率に基づいて、前記学習用データの最適個数を決定し、前記学習用データの数を前記最適個数として、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データを変えて再度繰り返し、前記再度の繰り返しにより得られた複数の前記正診率に基づいて、前記再度の繰り返しにより得られた複数の前記識別器から最適識別器を決定する請求項５に記載の受精卵の画像診断用識別器の作成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、受精卵の画像診断システム、受精卵の画像診断プログラム及び受精卵の画像診断用識別器の作成方法に関する。
【背景技術】
【０００２】
我が国は少子化の危機に瀕している一方で、不妊に悩む患者は多い。不妊患者の妊娠成立を図る治療として、体外受精又は顕微授精が挙げられる。体外受精は、卵巣から卵細胞を取り出し（採卵）、体外で精子を与えて受精卵を作る治療であり、顕微授精は、体外で顕微鏡を見ながら精子を直接卵細胞に注入して受精させ受精卵を作る治療である。いずれの治療においても、受精卵を培養して成長させ、一旦凍結して保存し、患者の条件が良い時に解凍して、子宮に戻して妊娠成立を図るというものである。
【０００３】
子宮に戻すべき受精卵の選択は、解凍前と解凍後に顕微鏡で見て医師、検査技師又は培養士が顕微鏡観察により主観的に判断しており、正常と判断した受精卵であっても、現在、妊娠に至るのは高々３０％に留まっている。すなわち、主観的な判断は精度が低く、判断手法を工夫したとしても精度を上げるには一定の限界があった。
【０００４】
他方、１回当たりの治療費は高額であり、かつ自費負担となる。このため、受精卵選択の判断の精度が低いと、治療回数も増え時間が浪費され、心理的ストレスに加え、経済的ストレスも多大となる。このような背景の下、不妊の治療においては、できるだけ早く、かつ経済的負担が少なく妊娠に至ることが患者にも社会的にも求められている。
【０００５】
この点、医療の分野においては、コンピュータを用いて診断画像の分析、評価等を行う各種技術が提案されている。例えば、特許文献１には、診断画像の読影に際して、画像の特徴量を用いた類似症例の検索をコンピュータで実現させ、診断精度を向上させる診断支援技術が提案されている。特許文献２には、撮像された画像のみを用いて、受精卵内の細胞を識別することなく受精卵の品質評価に有益な情報の提示をコンピュータで実現させる受精卵品質評技術が提案されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００２−２３０５１８号公報
【特許文献２】特開２０１０−１８１４０２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１に記載の技術は、受精卵の画像を対象とするものではなく、対象画像が臓器を対象としたＣＴ画像であり、類似症例の検索を可能にして診断精度を向上させるというものであった。すなわち、特許文献１の記載からは、受精卵の画像に対して妊娠成立の可否を判断する手法を導き出すのは困難であった。他方、特許文献２に記載の技術は、受精卵の画像を対象とするものであるが、受精卵を品質評価する技術に留まり、妊娠成立の可否を直接的に判断する技術ではなかった。
【０００８】
本発明は、前記のような従来の問題を解決するものであり、主観的判断をすることなく、受精卵の画像データだけで正常児獲得の可否又は妊娠成立を予測可能な診断結果が得られる受精卵の画像診断システム、受精卵の画像診断プログラム及び受精卵の画像診断用識別器の作成方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
前記目的を達成するために、本発明の画像診断システムは、機械学習の手法を用いた受精卵の画像診断システムであって、妊娠予後が判明している受精卵の画像データを格納する記憶手段と、受精卵の画像データを記憶手段に取り込むための受信手段と、前記画像データを解析する解析手段と、前記解析手段による解析結果を出力する出力手段とを備えており、前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、前記解析手段は、前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、前記学習用データから複数の識別器を作成し、前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする。
【００１０】
本発明の画像診断プロブラムは、コンピュータに、機械学習の手法を用いて受精卵の画像診断を実行させるための受精卵の画像診断プログラムであって、受精卵の画像データを記憶手段に取り込む取り込みステップと、前記画像データを解析する解析ステップと、前記解析手段による解析結果を出力する出力ステップとを前記コンピュータに実行させ、前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、前記解析ステップにおいて、前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、前記学習用データから複数の識別器を作成し、前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする。
【００１１】
本発明の画像診断用識別器の作成方法は、機械学習の手法を用いた受精卵の画像診断用識別器の作成方法であって、受精卵の画像データを記憶手段に取り込む取り込みステップと、前記画像データを解析する解析ステップと、前記解析手段による解析結果を出力する出力ステップとを備え、前記画像データは、生児獲得に至ったことが判明している正常画像データと、生児獲得に至らなかったことが判明している異常画像データとで構成される画像データ、又は妊娠成立に至ったことが判明している正常画像データと、妊娠成立に至らなかったことが判明している異常画像データとで構成される画像データであり、かつ患者の年齢情報を含んでおり、前記解析ステップにおいて、前記画像データを、前記患者の年齢情報に基づいて、年齢別の複数のグループに区分し、前記複数のグループから一つを選んだ解析対象グループについて、テストデータと学習用データとに分類し、前記学習用データから複数の識別器を作成し、前記各識別器で前記テストデータを診断して、前記テストデータの正診率を演算し、前記正診率に基づいて前記複数の前記識別器から最適識別器を決定し、前記分類から前記最適識別器の決定までの処理を、他の前記グループを解析対象グループに選んで実行し、これを繰り返して年齢別の前記最適識別器を決定することを特徴とする。
【００１２】
本発明によれば、受精卵の画像データに対し、撮影や測定等の機材を用いることなく、コンピュータのプログラムを実行させるだけで、主観的判断をすることなく診断が可能になる。そして、データを診断する識別器は、妊娠予後が判明している妊娠予後判明画像データに基づいて機械学習の手法により統計的に作成することに加え、年齢別に作成するので、診断結果は正常児獲得又は妊娠成立の可否が高い精度で予測可能なものとなる。
【００１３】
本発明の画像診断システムにおいては、前記解析手段は、前記最適識別器を決定する前に、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データの個数を変えて繰り返し、個数の異なる前記学習用データについて演算された前記正診率に基づいて、前記学習用データの最適個数を決定し、前記学習用データの数を前記最適個数として、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データを変えて再度繰り返し、前記再度の繰り返しにより得られた複数の前記正診率に基づいて、前記再度の繰り返しにより得られた複数の前記識別器から最適識別器を決定することが好ましい。
【００１４】
本発明の画像診断プロブラム及び画像診断用識別器の作成方法においては、前記解析ステップにおいて、前記解析手段は、前記最適識別器を決定する前に、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データの個数を変えて繰り返し、個数の異なる前記学習用データについて演算された前記正診率に基づいて、前記学習用データの最適個数を決定し、前記学習用データの数を前記最適個数として、前記識別器の作成及び前記正診率の演算を、選択する前記学習用データを変えて再度繰り返し、前記再度の繰り返しにより得られた複数の前記正診率に基づいて、前記再度の繰り返しにより得られた複数の前記識別器から最適識別器を決定することが好ましい。
【００１５】
前記本発明の好ましい構成によれば、正診率に基づいて、学習用データの最適個数を決定するので、機械学習の過学習による診断精度の低下を防止することができる。
【発明の効果】
【００１６】
本発明の効果は前記のとおりであり、本発明によれば、受精卵の画像データに対し、撮影や測定等の機材を用いることなく、コンピュータのプログラムを実行させるだけで、主観的判断をすることなく診断が可能になる。そして、データを診断する識別器は、妊娠予後が判明している妊娠予後判明画像データに基づいて機械学習の手法により統計的に作成することに加え、年齢別に作成するので、診断結果は正常児獲得又は妊娠成立の可否が高い精度で予測可能なものとなる。
【図面の簡単な説明】
【００１７】
【図１】本発明の一実施形態に係る受精卵の画像診断システムの構成図。
【図２】本発明の一実施形態に係る画像診断システムのフローチャート。
【図３】本発明の一実施形態に係る年齢区分の一例を示した図。
【図４】受精卵の正常画像データの一例の模式図。
【図５】受精卵の異常画像データの一例の模式図。
【図６】図３に示した正常画像データに加え、これを増幅した画像データを示した模式図。
【図７】本発明の一実施形態において、学習用データの最適個数の決定のプロセスを示すフローチャート。
【図８】本発明の一実施形態において、グループ１の学習用データの枚数と正診率との関係を示す図。
【図９】本発明の一実施形態において、最適識別器の決定のためのプロセスを示すフローチャート。
【図１０】本発明の一実施形態において、グループ２の学習用データの枚数と正診率との関係を示す図。
【図１１】本発明の一実施形態において、グループ３の学習用データの枚数と正診率との関係を示す図。
【図１２】本発明の一実施形態において、グループ４の学習用データの枚数と正診率との関係を示す図。
【図１３】本発明の一実施形態において、グループ５の学習用データの枚数と正診率との関係を示す図。
【図１４】本発明の一実施形態において、最適識別器と診断対象画像データとの関係を示す図。
【図１５】本発明の一実施形態に係る最適識別器と比較用識別器との正診率の比較結果を示した図。
【発明を実施するための形態】
【００１８】
以下、本発明の一実施形態について図面を参照しながら説明する。図１は、本発明の一実施形態に係る受精卵の画像診断システム（以下、単に「画像診断システム」という。）の構成図を示している。画像診断システム１は、受精卵の画像データに基づいて、妊娠成立の可否に関する判断情報を出力するシステムである。画像診断システム１は、インターネット８を介して、複数の医療機関１、２等に対して、受精卵の画像データ等の各種データの受け渡しを行う。解析手段６（ＣＰＵ：中央演算処理装置）により各種演算処理が行われる。キーボード等の入力手段４により各種データの入力、設定を行う。ディスプレイ、プリンタ等の出力手段５により、出力結果が出力される。受信手段２で受信した医療機関からの画像データは、医療機関毎に記憶手段７に格納される。解析手段６による解析結果は、送信手段３により複数の医療機関１、２等に向けて送信される。
【００１９】
画像診断システム１が取り扱う画像データは、体外受精又は顕微授精を行ってその妊娠予後が判明している画像データ（以下、「妊娠予後判明画像データ」という。）と、診断対象の画像データ（以下、「診断対象画像データ」である。診断対象画像データは、妊娠成立の可否が不明なデータである。図１において、受信手段２は、診断対象画像データ１１を複数の医療機関１、２等からインターネット８を介して個別に受信する。これとは別に、受信手段２は妊娠予後判明画像データ１０を受信する。両画像データは、例えばjpegフォーマットであり、いずれも記憶手段７に格納される。
【００２０】
前記の画像データはいずれも解析する前に予め画像処理を行なう。画像処理は、例えば背景のアーチファクトや格納容器の影などの消去、数値化された画像の標準化処理である。よって、解析手段６で解析される画像データは、予め画像処理がされたものである。画像処理は、解析前に実施されていればよく、同一プログラムで実施してもよく、画像処理用の別のプログラムで実施してもよい。通常はいずれの場合も、画像処理後の画像データを解析手段６が解析することになるが、同一プログラムで画像処理を実施する場合は、解析手段６に解析前の画像データに対し画像処理を行う機能を持たせてもよい。
【００２１】
解析手段６は診断対象画像データについて、妊娠成立の可否に関する判断情報を導き出す。このような判断情報を導き出す前提として、画像診断システム１には、妊娠予後判明画像データ１０に基づいて、後述する画像診断用の最適識別器が作成されている。すなわち、画像診断システム１による情報処理は、妊娠予後判明画像データ１０に基づいて最適識別器を作成する情報処理と、最適識別器を用いて診断対象画像データ診断する情報処理とに大別される。本実施形態では、この２つの情報処理を１つのプログラムで実行する例で説明するが、各情報処理を別のプログラムで実行してもよい。
【００２２】
画像診断システム１による情報処理の手順については、最初に図２を参照しながらその概要を説明する。図２は、本発明の一実施形態に係る画像診断システム１のフローチャートを示している。後に説明する図７は学習用データの最適個数の決定のプロセスを示すフローチャーを示しており、図９は最適識別器の決定のためのプロセスを示すフローチャートを示している。画像診断システム１は、コンピュータに画像診断プログラムを実行させて使用される。本実施形態では、本発明を画像診断システム１として説明するが、図２、図７及び図９に示した画像診断システム１の情報処理の流れは、画像診断プログラムがコンピュータに実行させるステップを示したものでもある。したがって、本実施形態で説明する画像診断システム１による情報処理の内容は、画像診断プログラムがコンピュータに実行させる内容でもあり、画像診断用識別器の作成方法のプロセスでもある。
【００２３】
図２において、初期設定（ステップ１００）の後、画像データを読み込む（ステップ１０１）。前記のとおり、画像診断システム１が取扱う画像データは、妊娠予後が判明している妊娠予後判明画像データ１０（図１参照）と、診断対象の診断対象画像データ１１（図１参照）の２種類である。両画像データを連続して読み込んでもよいが、少なくとも妊娠予後判明画像データ１０を読み込んでおけば、最適識別器の作成が可能になり（ステップ１０３）、この後に診断対象画像データ１１を読み込んでもよい。
【００２４】
ステップ１０２及びステップ１０３は、診断対象画像データ１１を診断する前の前段階のステップであり、本実施形態では、診断の精度を高めるために、妊娠予後判明画像データ１０（図１参照）を年齢区分し（ステップ１０２）、年齢区分毎に最適識別器を決定するようにしている（ステップ１０３）。このため、診断対象画像データ１１は患者の年齢情報を含んでいる。図３に年齢区分の一例を示している。図３の例では、図３に示した年齢区分で診断対象画像データ１１を５つのグループに分類する。
【００２５】
年齢区分毎に最適識別器が決定されると、診断対象画像データ１１（図１参照）の年齢に応じた最適識別器で同画像データを診断し（ステップ１０４）、解析結果を出力する（ステップ１０５）。例えば、図３のグループ１（４２歳以上）に属する診断対象画像データ１１に対しては、グループ１に対応した最適識別器で、同画像データを診断する。
【００２６】
最適識別器の決定のための解析処理及び最適識別器による画像データの診断のための解析処理はいずれも、人工知能による教師付機械学習のパターン認識を活用したものである。解析処理のプログラム言語としては例えばWolfram言語が挙げられる。機械学習のパターン認識についてはlogistic regression、naive Bayes、random forest、nearest neighbors、neural network、deep learningという認識方法が挙げられる。Wolfram言語にはこれらの認識方法を扱う関数がすでに用意されており、各認識方法から最適なものを選択すればよい。
【００２７】
次に画像診断システム１が取り扱う画像について説明する。図１において、妊娠予後判明画像データ１０は、正常画像データと異常画像データとに分類される。正常画像データは、生児獲得に至ったことが判明している画像データであり、異常画像データは、生児獲得に至らなかったことが判明している画像データである。生児獲得に至らなかった理由としては、臨床的には染色体異常、子宮因子などが考えられるが、異常画像データの分類に際し、生児獲得に至らなかった理由は問わない。
【００２８】
正常画像データ及び異常画像データのいずれについても、解析手段６で解析する際には形式を統一する。例えば、画像解像度を例えば７２pixel／inchで１１３×１１３pixelsの正方形になるように変換する。図４の符合１２は正常画像データの一例の模式図を示しており、図５の符合１３は異常画像データの一例の模式図を示している。各画像データは正方形に変換されているため、回転及び反転により、同じ正方形の画像データを増幅することができる。
【００２９】
図６は図４に示した正常画像データに加え、これを増幅した画像データを示している。１個の画像データは、図６に示したように、９０度回転、１８０度回転、２７０度回転により３個の画像が追加され、反転、反転９０度回転、反転１８０度回転、反転２７０度回転により４枚の画像が追加され、元の画像データを含めると８枚の画像データが得られる。このことは、図５に示した異常画像データ１３についても同様である。このように、画像データを増幅することにより、１個の画像データを８倍の枚数の画像データに増幅することができるので、画像データの特徴を多面的に表現でき、後に説明する識別器の精度を高めることができる。
【００３０】
画像データの形式の統一は、前記のものに限られるのではなく、適宜変更してもよい。例えば、各画像データを円形に変換し、回転及び反転により、画像データを増幅してもよい。円形画像データの場合、画像データを小さな角度で小刻みに回転させることにより、増幅枚数を増やすことができ、画像データの特徴をより多面的に表現でき、識別器の精度向上に有利になる。
【００３１】
画像診断システム１は、正常画像データ及び異常画像データを含む妊娠予後判明画像データから選んだ複数の学習用データを用いて識別器を作成する。本実施形態では、機械学習の過学習による診断精度の低下を防止するために、予め学習用データの最適個数を決定する。以下、図７を参照しながら、学習用データの最適個数の決定について説明する。
【００３２】
図７は、学習用データの最適個数の決定のプロセスを示すフローチャートである。便宜のため数値例を挙げながら説明するが、一例でありこれに限定されるものではない。本実施形態では、画像診断システム１は、図１において妊娠予後判明画像データ１０として、媒精後５日目の受精卵（胚盤胞）の画像を５６９１枚受信した。これらの画像は図３の例のように、５グループに年齢区分した（図２のステップ１０３）。前記のとおり、画像データの増幅により、画像データは任意の枚数に増幅される。
【００３３】
以下、５グループのそれぞれについて、学習用データの最適枚数を決定する。図７は、グループ１（４２歳以上）についての学習用データの最適枚数の決定のプロセスを示している。最初に、グループ１の年齢区分の画像データ６６９枚をテストデータ及び学習用データに分け、学習用データの中から一定枚数を検証用データとする（ステップ１１０）。この分け方は無作為であり、全画像データの２０％をテストデータ、残り８０％を学習用データとし、学習用データの２０％を検証用データとした。この学習用データからｎ１枚を選択し、このｎ１枚を増幅して学習用データｎ枚を作成し（ステップ１１１）、検証用データを含むｎ枚の学習用データを用いて識別器を作成する（ステップ１１２）。このステップでは、検証用データで調整しながら識別器を学習させて、ｎ枚の学習用データから最良の識別器を作成するようにしている。ｎ枚の学習用データの正常画像データと異常画像データの枚数の比率は限定せず様々な値をとる。
【００３４】
識別器とは入力されたデータをいくつかのクラスに分類する識別関数プログラムのことであり、ここではクラスは正常と異常の２分類である。本実施形態では識別器は画像診断用識別器であり、識別器の作成に検証用データを含む学習用データを用い、識別器の診断能の評価にテストデータを用いる（ステップ１１３）。学習用データから画像特徴を抽出し、正常と異常とをできるだけ区別できるようなパラメータを設定することで識別器を求める。識別器に診断対象の画像データを入力すると、正常か異常かのいずれかの結果が出力される。識別器の作成に際しては、様々なパラメータ（Ｌ２正則化値、学習画像サイズ、解像度等）を調整する。これらの各識別器を用いて、テストデータを診断し、正診率が高くかつ標準偏差が小さい識別器を暫定的な識別器とする。
【００３５】
診断は、正常又は異常の２種類であり、診断した画像データの枚数分の診断結果が得られる。検証用データは正常又は異常が判明している画像データであるため、診断結果の正否の判定が可能である。すなわち、診断結果が正常であり、当該テストデータも正常である場合と、診断結果が異常であり、当該テストデータも異常である場合は、正判定となる。他方、診断結果が正常であり、当該テストデータが異常である場合と、診断結果が異常であり、当該テストデータが正常である場合は否判定となる。本実施形態では、正診率を次の式（１）で定義する。
式（１）正診率＝正判定の画像データ枚数／画像データ総枚数
【００３６】
暫定的な識別器を作成した後は、この識別器でテストデータを診断し（ステップ１１３）、正診率を演算する（ステップ１１４）。テストデータについても、正常又は異常が判明している画像データであるため、診断結果の正否の判定が可能である。正診率の演算は、式（１）のとおりである。正診率の演算の後は、学習用データｍ枚を増加させ（ステップ１１５）、改めて学習用データｎ枚を選択し（ステップ１１１）、ステップ１１２〜ステップ１１１４を繰り返す。すなわち、新たな暫定的な識別器を繰り返し作成し、新たな識別器毎にテストデータを診断して正診率を演算する。
【００３７】
ステップ１１１〜ステップ１１５を繰り返すことにより、学習用データの枚数毎に暫定的な識別器が複数作成され、各識別器による正診率の演算結果も得られることになる。すなわち、学習用データの枚数毎に正診率が得られる。解析手段６は、正診率が高くかつ標準偏差が小さくなる観点から最適学習用データ枚数Ｎを決定する。
【００３８】
図８はグループ１（４２歳以上）の学習用データの枚数と正診率との関係を示す図である。学習用データが１７５０枚近傍の部分を例にとると、下端１６は平均値から標準偏差値を引いた値、上端１７は平均値に標準偏差値を加えた値を示しており、点１５は平均値を示している。下端１６と上端１７との間が離れているほど、ばらつきが大きく標準偏差も大きくなる。図８の結果では、学習用データの枚数が１２０００枚近傍のときに、正診率は高い値を維持し、かつばらつきも小さくなっており（Ａ部参照）、高い正診率を安定して確保できる。解析手段６は、正診率の平均値及び標準偏差を尺度として、最適学習用データ枚数Ｎとし１２１４４枚を決定した（ステップ１１６）。
【００３９】
グループ１〜５の学習用データの最適枚数Ｎを決定した後は、解析手段６は最適識別器の決定のための演算処理を行う。図９は、最適識別器の決定のためのプロセスを示すフローチャートである。最適識別器の決定はグループ１〜５の各グループについて行う。以下、グループ１の最適識別器の決定について説明する。最初に、グループ１の年齢区分の画像データ６６９枚をテストデータ及び学習用データに分け、学習用データの中から一定枚数を検証用データとする（ステップ１２０）。この分け方は無作為であり、全画像データの２０％をテストデータ、残り８０％を学習用データ（５３５枚）とし、学習用データの２０％を検証用データとした。
【００４０】
次に、５３５枚の学習用データを増幅してグループ１の最適学習用データ枚数である１２１４４枚の画像データとする（ステップ１２１）。１２１４４枚の学習用データの正常画像データと異常画像データの比率は限定せず様々な値をとる。この１２１４４枚の学習用データ（検証用データ含む）から識別器を作成する（ステップ１２２）。このステップでは、図７のステップ１１２と同様に、検証用データで調整しながら識別器を学習させて、Ｎ枚の学習用データから最良の識別器を作成する。
【００４１】
次に、識別器で診断を行うテストデータＮ２枚を抽出する（ステップ１２３）。テストデータは、グループ１の全データのうち学習用データに用いなかった残りの２０％に相当する１３５枚の画像データである。１３５枚のテストデータの内訳は、グループ１の全画像データの正常画像データと異常画像データの比率に応じた枚数とした。この１３５枚のテストデータについて、ステップ１２２で作成した識別器で診断を行う（ステップ１２４）。
【００４２】
診断要領は、図７のステップ１１３と同様であり、正常又は異常の２種類であり、１３５枚のテストデータに対し１３５個の診断結果が得られ、前記式（１）で定義した正診率の計算式により、１３５枚のテストデータの診断結果について正診率を演算する（ステップ１２５）。
【００４３】
正診率演算後に、パラメータ（Ｌ２正則化値、学習用データ等）を変化させて新たな識別器を作成し、改めて正診率を演算し、正診率の平均値及び標準偏差を尺度として、ステップ１２１で作成した識別器をより結果の良い識別器に置き換えてもよい。
【００４４】
ステップ１２４の正診率の演算が完了すると、新たにランダムに選んだ５３５枚の画像データから１２１４４枚の画像データを新たに作成し、ステップ１２０〜１２５を繰り返し、改めて１３５個の正診率が得られる。以後、同じ手順を繰り返すことにより、新たに選択した学習用データについてそれぞれ１３５個の正診率が得られる。解析手段６は、最も高い正診率をもたらした識別器を最適識別器として採用する（ステップ１２６）。このステップは、図２のステップ１０３に相当する。
【００４５】
以上の解析プロセスを経て、グループ１（４２歳以上）の年齢区分の最適識別器が作成されたことになる。以下、図７、図９に示した解析プロセスを順次繰り返して、グループ２〜５の年齢区分の最適識別器を作成する。すなわち、各グループの画像データを解析して、各グループ毎に最適学習用データ枚数Ｎを決定した後（図７のステップ１１６）、各グループ毎に最適識別器を決定する（図９のステップ１２６）。
【００４６】
各グループについて、図７のステップ１１１〜ステップ１１６を繰り返すことにより、各グループ毎に学習用データの枚数と正診率との関係が得られる。図１０〜図１３に、グループ２〜５の学習用データの枚数と正診率との関係を示している。図１０は、グループ２（４０歳以上４２歳未満）の年齢区分における学習用データの枚数と正診率との関係を示している。この結果を用いて、最適学習用データ枚数Ｎを決定する。決定要領はグループ１（４２歳以上）の場合と同様であり、正診率の平均値及び標準偏差を尺度として、最適学習用データ枚数Ｎを決定する。図１０の結果では、学習用データの枚数が８０００枚近傍のときに、正診率は高い値を維持し、かつばらつきも小さくなっており（Ｂ部参照）、解析手段６は、最適学習用データ枚数Ｎとして８０８５枚を決定した（ステップ１１６）。
【００４７】
最適学習用データ枚数Ｎの決定は、他のグループ３〜５についても、以下の通り同様である。図１１は、グループ３（３８歳以上４０歳未満）の年齢区分における学習用データの枚数と正診率との関係を示しており、図１１の結果では、学習用データの枚数が８０００枚近傍のときに、正診率は高い値を維持し、かつばらつきも小さくなっており（Ｃ部参照）、解析手段６は、最適学習用データ枚数Ｎとして７８４８枚を決定した。
【００４８】
図１２は、グループ４（３５歳以上３８歳未満）の年齢区分における学習用データの枚数と正診率との関係を示しており、図１２の結果では、学習用データの枚数が１７０００枚近傍のときに、正診率は高い値を維持し、かつばらつきも小さくなっており（Ｄ部参照）、解析手段６は、最適学習用データ枚数Ｎとして１７１１２枚を決定した。
【００４９】
図１３は、グループ５（３５歳未満）の年齢区分における学習用データの枚数と正診率との関係を示しており、図１３の結果では、学習用データの枚数が１６０００枚近傍のときに、正診率は高い値を維持し、かつばらつきも小さくなっており（Ｅ部参照）、解析手段６は、最適学習用データ枚数Ｎとして１６６６５枚を決定した。
【００５０】
各グループについて、最適学習用データ枚数Ｎが決定されると、図９のプロセスを経て、最適識別器が決定され、グループ１〜５の５つの年齢区分のそれぞれに対応した５つの最適識別器が決定されることになる。以後は診断対象画像データについて解析が行われる。すなわち、図２において、年齢区分毎に最適識別器を決定した後は（ステップ１０３）、この最適識別器で診断対象画像データを診断する（ステップ１０４）。
【００５１】
最適識別器による診断は、診断対象画像データの年齢に応じた最適識別器で診断する。図１４は、最適識別器と診断対象画像データとの関係を示している。最適識別器１〜５は、グループ１〜５（図３参照）に対応している。最適識別器１は、診断対象画像データのうち４２歳以上のデータを診断し、同様に最適識別器２〜５は、図１４に示した年齢区分の診断対象画像データを診断する。
【００５２】
最適識別器１〜５の効果を検証するために、比較解析を行った。この比較解析のために、前記の解析に用いた妊娠予後判明画像データ（媒精後５日目の受精卵（胚盤胞）の画像５６９１枚）を用いて、比較用識別器を新たに作成した。比較用識別器の作成手順は、最適識別器１〜５の作成手順と同様であるが、最適識別器１〜５の作成が年齢区分された画像データを用いたのに対して、比較用識別器の作成は年齢区分していない画像データを用いた点が異なっている。
【００５３】
図１５は、最適識別器１〜５と比較用識別器との正診率の比較結果を示している。図９を参照しながら説明したとおり、最適識別器１〜５はその導出過程において最も高い正診率が得られた識別器であり、同様に比較用識別器についてもその導出過程において最も高い正診率が得られた識別器である。図１５に示した最適識別器１〜５の正診率は最適識別器１〜５の導出過程で得られたものであり、比較用識別器の正診率は比較用識別器の導出過程で得られたものである。図１５によれば、最適識別器１〜４においては、正診率は比較用識別器の正診率を上回っており、年齢別に最適識別器を作成することの効果が確認できた。最適識別器５の正診率０．６３９は比較用識別器の正診率０．７２１を下回っていたが、正診率０．６３９は６割を超える正診率であり、十分高い正診率であるということができる。このため、本実施形態に係る年齢別の最適識別器による診断結果は正常児獲得の可否を高い精度で予測したものであるということができる。
【００５４】
前記実施形態では、妊娠予後判明画像データ１０は、正常画像データと異常画像データとに分類され、正常画像データは、生児獲得に至ったことが判明している画像データであり、異常画像データは、生児獲得に至らなかったことが判明している画像データであった。これに代えて、正常画像データを妊娠成立に至ったことが判明している画像データとし、異常画像データを妊娠成立に至らなかったことが判明している画像データとしてもよい。この場合、年齢別の最適識別器の作成手順は前記実施形態と同じであり、作成された年齢別の最適識別器による診断結果は妊娠成立の可否を高い精度で予測したものであるということができる。
【００５５】
以上のとおり、本発明は機械学習の手法を用いたものであり、受精卵の画像データに対し、撮影や測定等の機材を用いることなく、コンピュータのプログラムを実行させるだけで、主観的判断をすることなく診断が可能になる。そして、データを診断する識別器は、妊娠予後が判明している妊娠予後判明画像データに基づいて機械学習の手法により統計的に作成することに加え、年齢別に作成するので、診断結果は正常児獲得又は妊娠不成立の可否が高い精度で予測可能なものとなる。
【００５６】
以上、本発明の一実施形態について説明したが、最適識別器の決定のためのプロセスは適宜変更したものであってもよい。例えば前記実施形態では、識別器を作成する際に（図７のステップ１１２）、検証用データで調整しながら識別器を学習させて作成するという例で説明したがこれを省き、診断はテストデータでの診断（図７のステップ１１３）だけとしてもよい。また、図３では５グループの年齢区分の例を示したが、これに限るものではない。
【００５７】
前記実施形態においては、妊娠予後判明画像データ１０をインターネット８を介して受信する例で説明したが（図１参照）、これに限るものではなく、インターネット８を介在さずに、受信手段２で直接受信してもよい。また、妊娠予後判明画像データ１０は、画像診断システム１の外部にある記憶手段に一旦保存したものを受信手段２で受信してもよいが、顕微鏡に画像診断システム１を接続することで画像データを保存することなくリアルタイムで受信し、結果についてもリアルタイムで画像診断システム１の出力手段５で表示するようにしてもよい。
【００５８】
また、図７の学習用データｎ１枚（図７のステップ１１１）、図９の学習用データＮ１枚（図９のステップ１２１）における正常画像データと異常画像データの枚数の比率は限定せず様々な値をとる旨説明したが、両画像データの比率は、識別器の正診率がより高まる最適な比率を決定すればよく、例えば１：１、２：１又は１：２等でもよい。具体的には、最終的な比率を決定する前に、試験的に画像診断システム１を稼働させて、適宜比率を変更して、識別器の正診率を確認しながら、最適な比率を決定してもよい。
【符号の説明】
【００５９】
１  受精卵の画像診断システム
２  受信手段
３  送信手段
４  入力手段
５  出力手段
６  解析手段
７  記憶手段
【図１】
2020067680.tif 000003
【図２】
2020067680.tif 000004
【図３】
2020067680.tif 000005
【図４】
2020067680.tif 000006
【図５】
2020067680.tif 000007
【図６】
2020067680.tif 000008
【図７】
2020067680.tif 000009
【図８】
2020067680.tif 000010
【図９】
2020067680.tif 000011
【図１０】
2020067680.tif 000012
【図１１】
2020067680.tif 000013
【図１２】
2020067680.tif 000014
【図１３】
2020067680.tif 000015
【図１４】
2020067680.tif 000016
【図１５】
2020067680.tif 000017