知財求人 - 知財ポータルサイト「IP Force」
▶ 上海美沃精密儀器股▲フン▼有限公司の特許一覧 ▶ 復旦大学附属眼耳鼻喉科医院の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2023-09-26
(54)【発明の名称】マルチモーダルデータに基づく両眼円錐角膜の検査方法
(51)【国際特許分類】
A61B 3/107 20060101AFI20230919BHJP
G06T 7/00 20170101ALI20230919BHJP
【ＦＩ】
A61B3/107
G06T7/00 612
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2022541624
(86)(22)【出願日】2021-08-05
(85)【翻訳文提出日】2022-07-04
(86)【国際出願番号】 CN2021110812
(87)【国際公開番号】W WO2023272876
(87)【国際公開日】2023-01-05
(31)【優先権主張番号】202110717940.2
(32)【優先日】2021-06-28
(33)【優先権主張国・地域又は機関】CN
(81)【指定国・地域】
(71)【出願人】
【識別番号】522267826
【氏名又は名称】上海美沃精密儀器股▲フン▼有限公司
(71)【出願人】
【識別番号】522267837
【氏名又は名称】復旦大学附属眼耳鼻喉科医院
(74)【代理人】
【識別番号】110002468
【氏名又は名称】弁理士法人後藤特許事務所
(72)【発明者】
【氏名】沈 陽
(72)【発明者】
【氏名】周 行涛
(72)【発明者】
【氏名】李 慧傑
(72)【発明者】
【氏名】王 崇陽
(72)【発明者】
【氏名】陳 文光
(72)【発明者】
【氏名】趙 ▲セイ▼
(72)【発明者】
【氏名】李 美燕
(72)【発明者】
【氏名】▲セン▼ 芸勇
(72)【発明者】
【氏名】徐 海鵬
(72)【発明者】
【氏名】牛 凌凌
(72)【発明者】
【氏名】趙 武校
(72)【発明者】
【氏名】韓 田
【テーマコード（参考）】
4C316
5L096
【Ｆターム（参考）】
4C316AA03
4C316FB27
4C316FZ01
5L096BA06
5L096BA13
5L096EA03
5L096EA15
5L096EA16
5L096FA32
5L096FA33
5L096FA64
5L096FA66
5L096FA69
5L096GA30
5L096GA51
5L096HA11
(57)【要約】
本発明は、マルチモーダルデータに基づく両眼円錐角膜の診断方法に関し、両眼角膜の４屈折マップ及び角膜絶対高さデータを利用して、両眼の間の相互関係を考慮し、深度畳み込みネットワーク方法、従来の機械学習ｓｖｍ方法及び調整可能なＢＳＦの高さマップ増幅方法に基づいた病巣の感受性及び特異性の識別を組合わせ、その感受性及び特異性のバランスを取り、多次元で患者を単位とする円錐角膜の発病率を総合的に判断する。両眼データに合わせて人工選択特徴を含むだけでなく、深度ネットワーク自体がビッグデータから学習する特徴も含むから、診断方法は、より強いロバスト性及び正確性を有する。
【選択図】図１
【選択図】図１
【特許請求の範囲】
【請求項1】
両眼マルチモーダルデータを収集するステップであって、前記データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、前記角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、前記角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、分類された症例に基づいて、前記両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の前記４屈折マップにおける各トポグラフィーデータ及び前記角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、前記４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての前記４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての前記４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、前記角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得する、ステップ４）と、
前記分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、前記ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含むことを特徴とする、
マルチモーダルデータに基づく両眼円錐角膜の診断方法。
【請求項2】
前記分岐方法Ａは、具体的には、データをスケーリングするステップであって、ステップ３）で処理された全ての前記４屈折マップデータマトリックスを線形補間の方法で２２４×２２４のサイズにスケーリングするステップＡ－１と、
データを正規化するステップであって、ステップＡ－１のデータを７：３の比率に応じてトレーニングセットと検証セットに分け、次に前記トレーニングセットでそれぞれ前記４屈折マップデータマトリックスの平均値と標準偏差を計算し、それに応じて４つの平均値と４つの標準偏差を取得して、全ての症例の４屈折マップデータマトリックスを前記平均値と前記標準偏差で正規化するステップＡ－２と、
深度畳み込みネットワークの分類ネットワーク設計に基づいて、Ｒｅｓｎｅｔ５０分類ネットワークを用いて前記４屈折マップデータマトリックスに二分類を行うことにより、単眼の正常角膜及び円錐角膜を識別するステップＡ－３と、
分類モデルをトレーニングするステップであって、前記４屈折マップデータマトリックスをチャネルに応じて接続して、４チャネルの入力を取得し、データ増幅は、回転、平行移動、ランダムファジー前処理を採用し、損失関数は、二分類の交差エントロピー関数を採用し、ＭｏｂｉｌｅＮｅｔＶ３のＩＭＡＧＥＮＥＴデータセットでのトレーニング重み値を初期重み値とし、次に微調整トレーニングを行い、最終的に前記トレーニングセットと前記検証セットのｌｏｓｓ値の差が最小のトレーニング重み値をトレーニング結果として選択するステップＡ－４と、
モデル指標を評価するステップであって、前記検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法Ａで円錐角膜の感受性及び特異性を識別し、取得するステップＡ－５と、
結果を出力するステップであって、ステップＡ－５での前記トレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ａで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ａｌ）及びｐ（Ａｒ）として記録し、分類結果Ｐ（Ａ）＝ｐ（Ａｌ）（ｐ（Ａｌ）＞ｐ（Ａｒ））、ｐ（Ａｒ）（ｐ（Ａｌ）＜ｐ（Ａｒ））を出力するステップＡ－６と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の診断方法。
【請求項3】
前記分岐方法Ｂは、具体的には、角膜前面の軸方向曲率の固有値を計算するステップであって、角膜前面の軸方向曲率データマトリックスにおいて曲率の最大最小値点及び位置座標を計算し、直径６ｍｍ位置における上下屈折力差ＩＳ値を計算し、直径が４．５ｍｍの範囲内の表面不規則性ＳＲＩ値、表面非対称性ＳＡＩ値を計算するステップＢ－１と、
角膜前面の相対高さの固有値を計算するステップであって、角膜前面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－２と、
角膜後面の相対高さの固有値を計算するステップであって、角膜後面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－３と、
角膜厚さの固有値を計算するステップであって、角膜厚さデータマトリックスにおいて厚さの最小値及び位置座標を計算し、角膜頂点での厚さを計算するステップＢ－４と、
距離の固有値を計算するステップであって、ステップＢ－２における前記角膜前面の高さの最大値位置とステップＢ－３における前記角膜後面の高さの最大値位置との間の距離を計算し、ステップＢ－２における前記角膜前面の高さの最大値位置とステップＢ－４における前記角膜厚さの最小値位置との間の距離を計算し、ステップＢ－３における前記角膜後面の高さの最大値位置とステップＢ－４における前記角膜厚さの最小値位置との間の距離を計算するステップＢ－５と、
角膜体積の固有値を計算するステップであって、前記角膜厚さデータマトリックスを半径が４．５ｍｍの範囲内で体積積分して前記角膜体積を取得するステップＢ－６と、
ステップＢ－１からステップＢ－６までの全ての前記固有値に正規化処理を行い、かつ全ての正規化された症例データの固有値を７：３の比率に応じてトレーニングセットと検証セットに分けるステップＢ－７と、
ＳＶＭサポートベクターマシン二分類方法を用いてステップＢ－７で正規化された前記トレーニングセットの特徴データに特徴トレーニングを行い、前記特徴トレーニング中にＲＢＦ ｋｅｒｎｅｌを選択し、ｃｒｏｓｓ－ｖａｌｉｄａｔｉｏｎ及びｇｒｉｄ－ｓｅａｒｃｈを用いて最適なｃ及びｇを取得してデータをトレーニングするステップＢ－８と、
モデル指標を評価するステップであって、前記検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法で円錐角膜の感受性及び特異性を識別するステップＢ－９と、
結果を出力するステップであって、ステップＢ－９での前記トレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ｂで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｂｌ）及びｐ（Ｂｒ）として記録し、分類結果Ｐ（Ｂ）＝ｐ（Ｂｌ）（ｐ（Ｂｌ）＞ｐ（Ｂｒ））、ｐ（Ｂｒ）（ｐ（Ｂｌ）＜ｐ（Ｂｒ））を出力するステップＢ－１０と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の診断方法。
【請求項4】
前記分岐方法Ｃは、具体的には、前記角膜前面の標準的相対高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜前面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、前記角膜前面データと取得された最適フィッティング球面との間の高さ差を前記角膜前面の標準的相対高さデータとするステップＣ－１と、
前記角膜前面の特徴的高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜前面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在の前記ＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得して、前記角膜前面データと前記取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を前記角膜前面の特徴的相対高さデータとするステップＣ－２と、
前記角膜前面の増幅高さデータを取得するステップであって、ステップＣ－２で取得された前記標準的相対高さデータとステップＣ－３で取得された前記１１組の特徴的相対高さデータを差分して、１１組の角膜前面増幅データを取得するステップＣ－３と、
前記角膜後面の標準的相対高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜後面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、前記角膜後面データと取得された最適フィッティング球面との間の高さ差を前記角膜後面の標準的相対高さデータとするステップＣ－４と、
前記角膜後面の特徴的高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜後面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在の前記ＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得し、前記角膜後面データと取得された異なる前記ＢＦＳ最適フィッティング球面との間の高さ差を前記角膜後面の特徴的相対高さデータとするステップＣ－５と、
前記角膜後面の増幅高さデータを取得するステップであって、ステップＣ－４で取得された前記角膜後面の標準的相対高さデータとステップＣ－５で取得された前記１１組の特徴的相対高さデータを差分して、１１組の角膜後面増幅データを取得するステップＣ－６と、
ステップＣ－３及びＣ－６で合計して取得された２２組の前記角膜前後面の増幅データマトリックスを特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するステップＣ－７と、
本分岐方法Ｃで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｃｌ）とｐ（Ｃｒ）として記録するステップであって、ｐ（Ｃｌ）とｐ（Ｃｒ）は、いずれも現在の症例に基づいて計算された各組の増幅データとステップＣ－７で取得された前記臨界閾値との間の差分値を重み比率として加算して取得され、分類結果Ｐ（Ｃ）＝ｐ（Ｃｌ）（ｐ（Ｃｌ）＞ｐ（Ｃｒ））、ｐ（Ｃｒ）（ｐ（Ｃｌ）＜ｐ（Ｃｒ））を出力するステップＣ－８と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の診断方法。
【請求項5】
前記分岐方法Ｄは、具体的には、データ方位を統一するステップであって、右眼の４屈折マップデータマトリックスを縦列方向に沿って鏡像し、前記左右眼の４屈折マップデータマトリックスの鼻側及び側頭側の方位を統一するステップＤ－１と、
４屈折マップｄｉｆｆグラフマトリックスを取得するステップであって、それぞれ前記左右眼の４屈折マップデータマトリックスに対してポイントツーポイント差分を行った後に絶対値を取得し、４屈折マップｄｉｆｆグラフデータマトリックスを取得するステップＤ－２と、
ｄｉｆｆデータ特徴を計算するステップであって、直径が６ｍｍのデータ範囲内の前記４屈折マップｄｉｆｆグラフデータマトリックスにおける全てのデータの平均値、最大値、標準偏差をそれぞれ計算して特徴量とするステップＤ－３と、
ステップＤ－３で取得された１２組の前記左右眼角膜の４屈折マップｄｉｆｆグラフの平均値、最大値、標準偏差を特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するか、又はＳＶＭ分類方法を採用し、全てのタイプのｄｉｆｆデータの特徴に正規化トレーニングとテストを行い、最適な感受性及び特異性を与えるステップＤ－４と、
本分岐方法で取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれＰ（Ｄ）として記録するステップであって、前記Ｐ（Ｄ）は、現在の症例に基づいて計算された各組のｄｉｆｆデータ固有値とステップＤ－４で取得された前記臨界閾値との間の差分値を重み比率として加算して取得されるステップＤ－５と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の診断方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、眼科データ診断技術に関し、特にマルチモーダルデータに基づく両眼円錐角膜診断方法に関する。
【背景技術】
【0002】
人間の目の円錐角膜は、角膜が拡張し、角膜の中央が薄くなって前方に突出し、円錐形を呈するという特徴がある臨床眼科疾患である。このような疾患は、屈折矯正手術の禁忌であり、単眼又は両眼に発生し、かつ一般的に視力の著しい低下をもたらす。円錐角膜は、角膜の後面で最初に発病し、さらに角膜の前面に徐々に発展することが多い。
【0003】
現在の円錐角膜の診断及び治療は、既に角膜病学、屈折矯正手術学、検眼学等が密接に協力する臨床疾患学に発展している。一般的に角膜トポグラフィーに基づいて円錐角膜診断を行う方法は、臨床統計学的方法であり、角膜トポグラフィーの形態的特徴を利用して臨床パラメータ及び病歴に合わせて円錐角膜の確定診断及び臨床病期の分類を行うことが多い。統計学のモデルで取得された結果は、サマリーパラメータであり、かつ既知の診断結果データセットに基づいてパラメータ境界線が取得され、例えば、広く応用されるＫＩＳＡ指数、ＩＳ指数、ＳＲＩ／ＳＡＩ指数等である。このような方法は、プラットフォームデータに制限されることが多く、人為的に認められた有限個の単眼特徴に過度に依存し、両眼の間の相互関係を無視し、かつ異なる指数における感受性及び特異性が異なり、早期円錐角膜と不完全型円錐角膜に良好な判断を与えることができない。
【発明の概要】
【発明が解決しようとする課題】
【0004】
本発明は、両眼円錐角膜の早期診断の問題に対して、マルチモーダルデータに基づく両眼円錐角膜の診断方法を提供する。マルチモーダル屈折トポグラフィーのマトリックスデータに基づいて、畳み込みニューラルネットワーク法、固有値ＳＶＭ法、両眼対比法及び調整可能なＢＦＳの増幅トポグラフィー法を組み合わせて、両眼の総合円錐角膜の診断結果を与え、特に早期の後円錐角膜及び不完全型円錐角膜のスクリーニング及び診断に対してより高いロバスト性及び正確性を有する。
【課題を解決するための手段】
【0005】
本発明の技術手段は、以下のとおりである。マルチモーダルデータに基づく両眼円錐角膜の診断方法は、具体的には、
両眼マルチモーダルデータを収集するステップであって、データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、
分類された症例に基づいて、両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の４屈折マップにおける各トポグラフィーデータ及び角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得するステップ４）と、
分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含む。
両眼マルチモーダルデータを収集するステップであって、データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、
分類された症例に基づいて、両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の４屈折マップにおける各トポグラフィーデータ及び角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得するステップ４）と、
分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含む。
【0006】
好ましくは、前記分岐方法Ａは、具体的には、
データをスケーリングするステップであって、ステップ３）で処理された全ての４屈折マップデータマトリックスを線形補間の方法で２２４×２２４のサイズにスケーリングするステップＡ－１と、
データを正規化するステップであって、ステップＡ－１のデータを７：３の比率に応じてトレーニングセットと検証セットに分け、次にトレーニングセットでそれぞれ４屈折マップデータマトリックスの平均値と標準偏差を計算し、それに応じて４つの平均値と４つの標準偏差を取得して、全ての症例の４屈折マップデータマトリックスを平均値と標準偏差で正規化するステップＡ－２と、
深度畳み込みネットワークの分類ネットワーク設計に基づいて、Ｒｅｓｎｅｔ５０分類ネットワークを用いて４屈折マップデータマトリックスに二分類を行うことにより、単眼の正常角膜及び円錐角膜を識別するステップＡ－３と、
分類モデルをトレーニングするステップであって、４屈折マップデータマトリックスをチャネルに応じて接続して、４チャネルの入力を取得し、データ増幅は、回転、平行移動、ランダムファジー前処理を採用し、損失関数は、二分類の交差エントロピー関数を採用し、ＭｏｂｉｌｅＮｅｔＶ３のＩＭＡＧＥＮＥＴデータセットでのトレーニング重み値を初期重み値とし、次に微調整トレーニングを行い、最終的にトレーニングセットと検証セットのｌｏｓｓ値の差が最小のトレーニング重み値をトレーニング結果として選択するステップＡ－４と、
モデル指標を評価するステップであって、検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法Ａで円錐角膜の感受性及び特異性を識別し、取得するステップＡ－５と、
結果を出力するステップであって、ステップＡ－５でのトレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ａで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ａｌ）及びｐ（Ａｒ）として記録し、分類結果Ｐ（Ａ）＝ｐ（Ａｌ）（ｐ（Ａｌ）＞ｐ（Ａｒ））、ｐ（Ａｒ）（ｐ（Ａｌ）＜ｐ（Ａｒ））を出力するステップＡ－６と、を実現する。
データをスケーリングするステップであって、ステップ３）で処理された全ての４屈折マップデータマトリックスを線形補間の方法で２２４×２２４のサイズにスケーリングするステップＡ－１と、
データを正規化するステップであって、ステップＡ－１のデータを７：３の比率に応じてトレーニングセットと検証セットに分け、次にトレーニングセットでそれぞれ４屈折マップデータマトリックスの平均値と標準偏差を計算し、それに応じて４つの平均値と４つの標準偏差を取得して、全ての症例の４屈折マップデータマトリックスを平均値と標準偏差で正規化するステップＡ－２と、
深度畳み込みネットワークの分類ネットワーク設計に基づいて、Ｒｅｓｎｅｔ５０分類ネットワークを用いて４屈折マップデータマトリックスに二分類を行うことにより、単眼の正常角膜及び円錐角膜を識別するステップＡ－３と、
分類モデルをトレーニングするステップであって、４屈折マップデータマトリックスをチャネルに応じて接続して、４チャネルの入力を取得し、データ増幅は、回転、平行移動、ランダムファジー前処理を採用し、損失関数は、二分類の交差エントロピー関数を採用し、ＭｏｂｉｌｅＮｅｔＶ３のＩＭＡＧＥＮＥＴデータセットでのトレーニング重み値を初期重み値とし、次に微調整トレーニングを行い、最終的にトレーニングセットと検証セットのｌｏｓｓ値の差が最小のトレーニング重み値をトレーニング結果として選択するステップＡ－４と、
モデル指標を評価するステップであって、検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法Ａで円錐角膜の感受性及び特異性を識別し、取得するステップＡ－５と、
結果を出力するステップであって、ステップＡ－５でのトレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ａで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ａｌ）及びｐ（Ａｒ）として記録し、分類結果Ｐ（Ａ）＝ｐ（Ａｌ）（ｐ（Ａｌ）＞ｐ（Ａｒ））、ｐ（Ａｒ）（ｐ（Ａｌ）＜ｐ（Ａｒ））を出力するステップＡ－６と、を実現する。
【0007】
好ましくは、前記分岐方法Ｂは、具体的には、
角膜前面の軸方向曲率の固有値を計算するステップであって、角膜前面の軸方向曲率データマトリックスにおいて曲率の最大最小値点及び位置座標を計算し、直径６ｍｍ位置における上下屈折力差ＩＳ値を計算し、直径が４．５ｍｍの範囲内の表面不規則性ＳＲＩ値、表面非対称性ＳＡＩ値を計算するステップＢ－１と、
角膜前面の相対高さの固有値を計算するステップであって、角膜前面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－２と、
角膜後面の相対高さの固有値を計算するステップであって、角膜後面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－３と、
角膜厚さの固有値を計算するステップであって、角膜厚さデータマトリックスにおいて厚さの最小値及び位置座標を計算し、角膜頂点での厚さを計算するステップＢ－４と、
距離の固有値を計算するステップであって、ステップＢ－２における角膜前面の高さの最大値位置とステップＢ－３における角膜後面の高さの最大値位置との間の距離を計算し、ステップＢ－２における角膜前面の高さの最大値位置とステップＢ－４における角膜厚さの最小値位置との間の距離を計算し、ステップＢ－３における角膜後面の高さの最大値位置とステップＢ－４における角膜厚さの最小値位置との間の距離を計算するステップＢ－５と、
角膜体積の固有値を計算するステップであって、角膜厚さデータマトリックスを半径が４．５ｍｍの範囲内で体積積分して角膜体積を取得するステップＢ－６と、
ステップＢ－１からステップＢ－６までの全ての固有値に正規化処理を行い、かつ全ての正規化された症例データの固有値を７：３の比率に応じてトレーニングセットと検証セットに分けるステップＢ－７と、
ＳＶＭサポートベクターマシン二分類方法を用いてステップＢ－７で正規化されたトレーニングセットの特徴データに特徴トレーニングを行い、特徴トレーニング中にＲＢＦ ｋｅｒｎｅｌを選択し、ｃｒｏｓｓ－ｖａｌｉｄａｔｉｏｎ及びｇｒｉｄ－ｓｅａｒｃｈを用いて最適なｃ及びｇを取得してデータをトレーニングするステップＢ－８と、
モデル指標を評価するステップであって、検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法で円錐角膜の感受性及び特異性を識別するステップＢ－９と、
結果を出力するステップであって、ステップＢ－９でのトレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ｂで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｂｌ）及びｐ（Ｂｒ）として記録し、分類結果Ｐ（Ｂ）＝ｐ（Ｂｌ）（ｐ（Ｂｌ）＞ｐ（Ｂｒ））、ｐ（Ｂｒ）（ｐ（Ｂｌ）＜ｐ（Ｂｒ））を出力するステップＢ－１０と、を実現する。
角膜前面の軸方向曲率の固有値を計算するステップであって、角膜前面の軸方向曲率データマトリックスにおいて曲率の最大最小値点及び位置座標を計算し、直径６ｍｍ位置における上下屈折力差ＩＳ値を計算し、直径が４．５ｍｍの範囲内の表面不規則性ＳＲＩ値、表面非対称性ＳＡＩ値を計算するステップＢ－１と、
角膜前面の相対高さの固有値を計算するステップであって、角膜前面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－２と、
角膜後面の相対高さの固有値を計算するステップであって、角膜後面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－３と、
角膜厚さの固有値を計算するステップであって、角膜厚さデータマトリックスにおいて厚さの最小値及び位置座標を計算し、角膜頂点での厚さを計算するステップＢ－４と、
距離の固有値を計算するステップであって、ステップＢ－２における角膜前面の高さの最大値位置とステップＢ－３における角膜後面の高さの最大値位置との間の距離を計算し、ステップＢ－２における角膜前面の高さの最大値位置とステップＢ－４における角膜厚さの最小値位置との間の距離を計算し、ステップＢ－３における角膜後面の高さの最大値位置とステップＢ－４における角膜厚さの最小値位置との間の距離を計算するステップＢ－５と、
角膜体積の固有値を計算するステップであって、角膜厚さデータマトリックスを半径が４．５ｍｍの範囲内で体積積分して角膜体積を取得するステップＢ－６と、
ステップＢ－１からステップＢ－６までの全ての固有値に正規化処理を行い、かつ全ての正規化された症例データの固有値を７：３の比率に応じてトレーニングセットと検証セットに分けるステップＢ－７と、
ＳＶＭサポートベクターマシン二分類方法を用いてステップＢ－７で正規化されたトレーニングセットの特徴データに特徴トレーニングを行い、特徴トレーニング中にＲＢＦ ｋｅｒｎｅｌを選択し、ｃｒｏｓｓ－ｖａｌｉｄａｔｉｏｎ及びｇｒｉｄ－ｓｅａｒｃｈを用いて最適なｃ及びｇを取得してデータをトレーニングするステップＢ－８と、
モデル指標を評価するステップであって、検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法で円錐角膜の感受性及び特異性を識別するステップＢ－９と、
結果を出力するステップであって、ステップＢ－９でのトレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ｂで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｂｌ）及びｐ（Ｂｒ）として記録し、分類結果Ｐ（Ｂ）＝ｐ（Ｂｌ）（ｐ（Ｂｌ）＞ｐ（Ｂｒ））、ｐ（Ｂｒ）（ｐ（Ｂｌ）＜ｐ（Ｂｒ））を出力するステップＢ－１０と、を実現する。
【0008】
好ましくは、前記分岐方法Ｃは、具体的には、
角膜前面の標準的相対高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜前面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、角膜前面データと取得された最適フィッティング球面との間の高さ差を角膜前面の標準的相対高さデータとするステップＣ－１と、
角膜前面の特徴的高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜前面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得して、現在のＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得して、角膜前面データと取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を角膜前面の特徴的相対高さデータとするステップＣ－２と、
角膜前面の増幅高さデータを取得するステップであって、ステップＣ－２で取得された標準的相対高さデータとステップＣ－３で取得された１１組の特徴的相対高さデータを差分して、１１組の角膜前面増幅データを取得するステップＣ－３と、
角膜後面の標準的相対高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜後面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、角膜後面データと取得された最適フィッティング球面との間の高さ差を角膜後面の標準的相対高さデータとするステップＣ－４と、
角膜後面の特徴的高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜後面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在のＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得し、角膜後面データと取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を角膜後面の特徴的相対高さデータとするステップＣ－５と、
角膜後面の増幅高さデータを取得するステップであって、ステップＣ－４で取得された角膜後面の標準的相対高さデータとステップＣ－５で取得された１１組の特徴的相対高さデータを差分して、１１組の角膜後面増幅データを取得するステップＣ－６と、
ステップＣ－３及びＣ－６で合計して取得された２２組の角膜前後面の増幅データマトリックスを特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するステップＣ－７と、
本分岐方法Ｃで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｃｌ）とｐ（Ｃｒ）として記録するステップであって、ｐ（Ｃｌ）とｐ（Ｃｒ）は、いずれも現在の症例に基づいて計算された各組の増幅データとステップＣ－７で取得された臨界閾値との間の差分値を重み比率として加算して取得され、分類結果Ｐ（Ｃ）＝ｐ（Ｃｌ）（ｐ（Ｃｌ）＞ｐ（Ｃｒ））、ｐ（Ｃｒ）（ｐ（Ｃｌ）＜ｐ（Ｃｒ））を出力するステップＣ－８と、を実現する。
角膜前面の標準的相対高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜前面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、角膜前面データと取得された最適フィッティング球面との間の高さ差を角膜前面の標準的相対高さデータとするステップＣ－１と、
角膜前面の特徴的高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜前面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得して、現在のＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得して、角膜前面データと取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を角膜前面の特徴的相対高さデータとするステップＣ－２と、
角膜前面の増幅高さデータを取得するステップであって、ステップＣ－２で取得された標準的相対高さデータとステップＣ－３で取得された１１組の特徴的相対高さデータを差分して、１１組の角膜前面増幅データを取得するステップＣ－３と、
角膜後面の標準的相対高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜後面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、角膜後面データと取得された最適フィッティング球面との間の高さ差を角膜後面の標準的相対高さデータとするステップＣ－４と、
角膜後面の特徴的高さデータを取得するステップであって、角膜前後面の絶対高さデータに対して、角膜後面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在のＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得し、角膜後面データと取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を角膜後面の特徴的相対高さデータとするステップＣ－５と、
角膜後面の増幅高さデータを取得するステップであって、ステップＣ－４で取得された角膜後面の標準的相対高さデータとステップＣ－５で取得された１１組の特徴的相対高さデータを差分して、１１組の角膜後面増幅データを取得するステップＣ－６と、
ステップＣ－３及びＣ－６で合計して取得された２２組の角膜前後面の増幅データマトリックスを特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するステップＣ－７と、
本分岐方法Ｃで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｃｌ）とｐ（Ｃｒ）として記録するステップであって、ｐ（Ｃｌ）とｐ（Ｃｒ）は、いずれも現在の症例に基づいて計算された各組の増幅データとステップＣ－７で取得された臨界閾値との間の差分値を重み比率として加算して取得され、分類結果Ｐ（Ｃ）＝ｐ（Ｃｌ）（ｐ（Ｃｌ）＞ｐ（Ｃｒ））、ｐ（Ｃｒ）（ｐ（Ｃｌ）＜ｐ（Ｃｒ））を出力するステップＣ－８と、を実現する。
【0009】
好ましくは、前記分岐方法Ｄは、具体的には、
データ方位を統一するステップであって、右眼の４屈折マップデータマトリックスを縦列方向に沿って鏡像し、左右眼の４屈折マップデータマトリクスの鼻側及び側頭側の方位を統一するステップＤ－１と、
４屈折マップｄｉｆｆグラフマトリックスを取得するステップであって、それぞれ左右眼の４屈折マップデータマトリックスに対してポイントツーポイント差分を行った後に絶対値を取得し、４屈折マップｄｉｆｆグラフデータマトリックスを取得するステップＤ－２と、
ｄｉｆｆデータ特徴を計算するステップであって、直径が６ｍｍのデータ範囲内の４屈折マップｄｉｆｆグラフデータマトリックスにおける全てのデータの平均値、最大値、標準偏差をそれぞれ計算して特徴量とするステップＤ－３と、
ステップＤ－３で取得された１２組の左右眼角膜の４屈折マップｄｉｆｆグラフの平均値、最大値、標準偏差を特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するか、又はＳＶＭ分類方法を採用し、全てのタイプのｄｉｆｆデータの特徴に正規化トレーニングとテストを行い、最適な感受性及び特異性を与えるステップＤ－４と、
本分岐方法で取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれＰ（Ｄ）として記録するステップであって、Ｐ（Ｄ）は、現在の症例に基づいて計算された各組のｄｉｆｆデータ固有値とステップＤ－４で取得された臨界閾値との間の差分値を重み比率として加算して取得されるステップＤ－５と、を実現する。
データ方位を統一するステップであって、右眼の４屈折マップデータマトリックスを縦列方向に沿って鏡像し、左右眼の４屈折マップデータマトリクスの鼻側及び側頭側の方位を統一するステップＤ－１と、
４屈折マップｄｉｆｆグラフマトリックスを取得するステップであって、それぞれ左右眼の４屈折マップデータマトリックスに対してポイントツーポイント差分を行った後に絶対値を取得し、４屈折マップｄｉｆｆグラフデータマトリックスを取得するステップＤ－２と、
ｄｉｆｆデータ特徴を計算するステップであって、直径が６ｍｍのデータ範囲内の４屈折マップｄｉｆｆグラフデータマトリックスにおける全てのデータの平均値、最大値、標準偏差をそれぞれ計算して特徴量とするステップＤ－３と、
ステップＤ－３で取得された１２組の左右眼角膜の４屈折マップｄｉｆｆグラフの平均値、最大値、標準偏差を特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するか、又はＳＶＭ分類方法を採用し、全てのタイプのｄｉｆｆデータの特徴に正規化トレーニングとテストを行い、最適な感受性及び特異性を与えるステップＤ－４と、
本分岐方法で取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれＰ（Ｄ）として記録するステップであって、Ｐ（Ｄ）は、現在の症例に基づいて計算された各組のｄｉｆｆデータ固有値とステップＤ－４で取得された臨界閾値との間の差分値を重み比率として加算して取得されるステップＤ－５と、を実現する。
【発明の効果】
【0010】
本発明の有益な効果は、以下のとおりである。本発明に係るマルチモーダルデータに基づく両眼円錐角膜の診断方法によれば、両眼の間の相互関係を考慮し、深度畳み込みネットワーク方法、従来の機械学習ｓｖｍ方法及び調整可能なＢＳＦの高さマップ増幅方法に基づいた病巣の感受性及び特異性の識別を組み合わせて、その感受性及び特異性のバランスを取り、多次元で患者を単位とする円錐角膜の発病率を総合的に判断する。両眼データに合わせて人工選択特徴を含むだけでなく、深度ネットワーク自体がビッグデータから学習する特徴も含むから、診断方法は、より強いロバスト性及び正確性を有する。
【図面の簡単な説明】
【0011】
【図1】本発明に係るマルチモーダルデータに基づく両眼円錐角膜の診断方法のフローチャートである。
【図2】角膜の４屈折マップである。
【図3】角膜前後面の絶対高さマップである。
【図4】本発明に係る方法における分類ネットワークの構造図である。
【図5】左右眼角膜の４屈折マップＤｉｆｆグラフである。
【発明を実施するための形態】
【0012】
以下、図面及び具体的な実施例を参照しながら本発明を詳細に説明する。本実施例は、本発明の技術的解決手段を前提として実施され、詳細な実施形態及び具体的な操作過程を提供するが、本発明の保護範囲は、下記実施例に限定されない。
【0013】
図１に示すマルチモーダルデータに基づく両眼円錐角膜の診断方法のフローチャートは、具体的には、以下のステップ１～５を含む。
【0014】
ステップ１では、両眼マルチモーダルデータを収集する。データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、前眼部三次元解析装置、又は前眼部ＯＣＴ測定装置により４屈折マップを取得する。図２に示すように、角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、図３に示すとおり、角膜絶対高さデータは、角膜前後面の絶対高さデータを含む。
【0015】
ステップ２では、分類された症例に基づいて、両眼マルチモーダルデータを分類された種類と関連付ける。分類数は、具体的な需要に依存し、例えば、円錐角膜及び正常角膜に二分類する。
【0016】
ステップ３では、角膜の４屈折マップ及び色対応コードに基づいて、０．０２ｍｍをデータステップサイズとしてトポグラフィー（４屈折マップにおける４枚のマップは、いずれもトポグラフィーと呼ばれる）の直径９ｍｍ範囲内の二次元フルサンプリングデータマトリクスを取得する。各トポグラフィーデータ及び角膜前後面の高さデータマトリックスのサイズは４５１×４５１である。
【0017】
ステップ４では、以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断する。４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄである。
【0018】
分岐方法Ａ：当該分岐方法は、単一の症例を基準とする。全ての４屈折マップマトリックスにおける各データは、いずれも計算に関与し、円錐角膜判定に影響を与える可能性のある全ての要因を最大限に保持する。機械の自己学習により大量データから病変を判定するという目的を達成する。これにより、人為的、かつ主観的な特徴選択に依存せず、病巣識別の特異性を強化することができる。
【0019】
ステップＡ－１では、データをスケーリングする。ステップ３で処理された全ての４屈折マップデータマトリックスを線形補間の方法で２２４×２２４のサイズにスケーリングする。
【0020】
ステップＡ－２では、データを正規化する。ステップＡ－１のデータを７：３の比率に応じてトレーニングセットと検証セットに分ける。次にトレーニングセットでそれぞれ４屈折マップデータマトリックスの平均値と標準偏差を計算し、それに応じて４つの平均値と４つの標準偏差を取得する。全ての症例の４屈折マップデータマトリックスを平均値と標準偏差で正規化する。
【0021】
ステップＡ－３では、深度畳み込みネットワークの分類ネットワーク設計に基づいて、Ｒｅｓｎｅｔ５０分類ネットワークを用いて４屈折マップデータマトリックスに二分類を行うことにより、単眼の正常角膜及び円錐角膜を識別する。構築されたＲｅｓｎｅｔ５０分類ネットワークモデルは、バックボーンネットワークが不変で、一番目の畳み込み層のチャネル入力を４に修正し、最後に全接続層を出力する出力数を２に修正し、ネットワーク構造を図４に示す。
【0022】
ステップＡ－４では、分類モデルをトレーニングする。４屈折マップデータマトリックスをチャネルに応じて接続して、４チャネルの入力を取得する。データ増幅は、回転、平行移動、ランダムファジー前処理等を採用する。損失関数は、二分類の交差エントロピー関数を採用する。ＭｏｂｉｌｅＮｅｔＶ３のＩＭＡＧＥＮＥＴデータセットでのトレーニング重み値を初期重み値とし、次に微調整トレーニングを行う。６０回の反復トレーニングを行い、初期学習率が０．０１であり、２０回目で、４０回目でそれぞれ学習率を１０倍低下させる。最終的にトレーニングセットと検証セットのｌｏｓｓ値の差が最小のトレーニング重み値をトレーニング結果として選択する。
【0023】
ステップＡ－５では、モデル指標を評価し、検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法Ａで円錐角膜の感受性及び特異性を識別する。
【0024】
ステップＡ－６では、結果を出力する。ステップＡ－５でのトレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ａで取得されたある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ａｌ）及びｐ（Ａｒ）として記録する。分類結果Ｐ（Ａ）＝ｐ（Ａｌ）、（ｐ（Ａｌ）＞ｐ（Ａｒ））；ｐ（Ａｒ）、（ｐ（Ａｌ）＜ｐ（Ａｒ））を出力する。
【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】
上述した実施例は、本発明のいくつかの実施形態を示し、その説明が具体的で詳細であるが、本発明の特許範囲を限定するものと理解すべきではない。なお、当業者にとっては、本発明の構想から逸脱しない前提で、さらにいくつかの変形及び改善を行うことができ、これらは、いずれも本発明の保護範囲に属する。したがって、本発明の特許の保護範囲は、添付された特許請求の範囲を基準とすべきである。
【図1】
【図2】
【図3】
【図4】
【図5】
【手続補正書】
【提出日】2022-07-04
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
両眼マルチモーダルデータを収集するステップであって、前記データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、前記角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、前記角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、分類された症例に基づいて、前記両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の前記４屈折マップにおける各トポグラフィーデータ及び前記角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、前記４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての前記４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての前記４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、前記角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得する、ステップ４）と、
前記分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、前記ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含むことを特徴とする、
マルチモーダルデータに基づく両眼円錐角膜の検査方法。
【請求項2】
前記分岐方法Ａは、具体的には、データをスケーリングするステップであって、ステップ３）で処理された全ての前記４屈折マップデータマトリックスを線形補間の方法で２２４×２２４のサイズにスケーリングするステップＡ－１と、
データを正規化するステップであって、ステップＡ－１のデータを７：３の比率に応じてトレーニングセットと検証セットに分け、次に前記トレーニングセットでそれぞれ前記４屈折マップデータマトリックスの平均値と標準偏差を計算し、それに応じて４つの平均値と４つの標準偏差を取得して、全ての症例の４屈折マップデータマトリックスを前記平均値と前記標準偏差で正規化するステップＡ－２と、
深度畳み込みネットワークの分類ネットワーク設計に基づいて、Ｒｅｓｎｅｔ５０分類ネットワークを用いて前記４屈折マップデータマトリックスに二分類を行うことにより、単眼の正常角膜及び円錐角膜を識別するステップＡ－３と、
分類モデルをトレーニングするステップであって、前記４屈折マップデータマトリックスをチャネルに応じて接続して、４チャネルの入力を取得し、データ増幅は、回転、平行移動、ランダムファジー前処理を採用し、損失関数は、二分類の交差エントロピー関数を採用し、ＭｏｂｉｌｅＮｅｔＶ３のＩＭＡＧＥＮＥＴデータセットでのトレーニング重み値を初期重み値とし、次に微調整トレーニングを行い、最終的に前記トレーニングセットと前記検証セットのｌｏｓｓ値の差が最小のトレーニング重み値をトレーニング結果として選択するステップＡ－４と、
モデル指標を評価するステップであって、前記検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法Ａで円錐角膜の感受性及び特異性を識別し、取得するステップＡ－５と、
結果を出力するステップであって、ステップＡ－５での前記トレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ａで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ａｌ）及びｐ（Ａｒ）として記録し、分類結果Ｐ（Ａ）＝ｐ（Ａｌ）（ｐ（Ａｌ）＞ｐ（Ａｒ））、ｐ（Ａｒ）（ｐ（Ａｌ）＜ｐ（Ａｒ））を出力するステップＡ－６と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の検査方法。
【請求項3】
前記分岐方法Ｂは、具体的には、角膜前面の軸方向曲率の固有値を計算するステップであって、角膜前面の軸方向曲率データマトリックスにおいて曲率の最大最小値点及び位置座標を計算し、直径６ｍｍ位置における上下屈折力差ＩＳ値を計算し、直径が４．５ｍｍの範囲内の表面不規則性ＳＲＩ値、表面非対称性ＳＡＩ値を計算するステップＢ－１と、
角膜前面の相対高さの固有値を計算するステップであって、角膜前面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－２と、
角膜後面の相対高さの固有値を計算するステップであって、角膜後面の相対高さデータマトリックスにおいて高さの最大値及び位置座標を計算するステップＢ－３と、
角膜厚さの固有値を計算するステップであって、角膜厚さデータマトリックスにおいて厚さの最小値及び位置座標を計算し、角膜頂点での厚さを計算するステップＢ－４と、
距離の固有値を計算するステップであって、ステップＢ－２における前記角膜前面の高さの最大値位置とステップＢ－３における前記角膜後面の高さの最大値位置との間の距離を計算し、ステップＢ－２における前記角膜前面の高さの最大値位置とステップＢ－４における前記角膜厚さの最小値位置との間の距離を計算し、ステップＢ－３における前記角膜後面の高さの最大値位置とステップＢ－４における前記角膜厚さの最小値位置との間の距離を計算するステップＢ－５と、
角膜体積の固有値を計算するステップであって、前記角膜厚さデータマトリックスを半径が４．５ｍｍの範囲内で体積積分して前記角膜体積を取得するステップＢ－６と、
ステップＢ－１からステップＢ－６までの全ての前記固有値に正規化処理を行い、かつ全ての正規化された症例データの固有値を７：３の比率に応じてトレーニングセットと検証セットに分けるステップＢ－７と、
ＳＶＭサポートベクターマシン二分類方法を用いてステップＢ－７で正規化された前記トレーニングセットの特徴データに特徴トレーニングを行い、前記特徴トレーニング中にＲＢＦ ｋｅｒｎｅｌを選択し、ｃｒｏｓｓ－ｖａｌｉｄａｔｉｏｎ及びｇｒｉｄ－ｓｅａｒｃｈを用いて最適なｃ及びｇを取得してデータをトレーニングするステップＢ－８と、
モデル指標を評価するステップであって、前記検証セットで予測を行い、次に真値と比較評価して、最後に当該分岐方法で円錐角膜の感受性及び特異性を識別するステップＢ－９と、
結果を出力するステップであって、ステップＢ－９での前記トレーニングセットに対するテストの感受性及び特異性が要求に達すると、本分岐方法Ｂで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｂｌ）及びｐ（Ｂｒ）として記録し、分類結果Ｐ（Ｂ）＝ｐ（Ｂｌ）（ｐ（Ｂｌ）＞ｐ（Ｂｒ））、ｐ（Ｂｒ）（ｐ（Ｂｌ）＜ｐ（Ｂｒ））を出力するステップＢ－１０と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の検査方法。
【請求項4】
前記分岐方法Ｃは、具体的には、前記角膜前面の標準的相対高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜前面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、前記角膜前面データと取得された最適フィッティング球面との間の高さ差を前記角膜前面の標準的相対高さデータとするステップＣ－１と、
前記角膜前面の特徴的高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜前面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在の前記ＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得して、前記角膜前面データと前記取得された異なるＢＦＳ最適フィッティング球面との間の高さ差を前記角膜前面の特徴的相対高さデータとするステップＣ－２と、
前記角膜前面の増幅高さデータを取得するステップであって、ステップＣ－２で取得された前記標準的相対高さデータとステップＣ－３で取得された前記１１組の特徴的相対高さデータを差分して、１１組の角膜前面増幅データを取得するステップＣ－３と、
前記角膜後面の標準的相対高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜後面の直径が８ｍｍの範囲内の絶対高さデータで球面フィッティングを行い、ＢＦＳ値を取得して、前記角膜後面データと取得された最適フィッティング球面との間の高さ差を前記角膜後面の標準的相対高さデータとするステップＣ－４と、
前記角膜後面の特徴的高さデータを取得するステップであって、前記角膜前後面の絶対高さデータに対して、前記角膜後面の直径が８ｍｍの範囲内の絶対高さデータを基準とし、最薄点位置の半径２ｍｍ範囲内のデータを除去して球面フィッティングを行い、ＢＦＳ値を取得し、現在の前記ＢＦＳを基準とし、０．２ｍｍをステップサイズとし、上下へそれぞれ５組オフセットして１１組の異なるＢＦＳ値を取得し、前記角膜後面データと取得された異なる前記ＢＦＳ最適フィッティング球面との間の高さ差を前記角膜後面の特徴的相対高さデータとするステップＣ－５と、
前記角膜後面の増幅高さデータを取得するステップであって、ステップＣ－４で取得された前記角膜後面の標準的相対高さデータとステップＣ－５で取得された前記１１組の特徴的相対高さデータを差分して、１１組の角膜後面増幅データを取得するステップＣ－６と、
ステップＣ－３及びＣ－６で合計して取得された２２組の前記角膜前後面の増幅データマトリックスを特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するステップＣ－７と、
本分岐方法Ｃで取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれｐ（Ｃｌ）とｐ（Ｃｒ）として記録するステップであって、ｐ（Ｃｌ）とｐ（Ｃｒ）は、いずれも現在の症例に基づいて計算された各組の増幅データとステップＣ－７で取得された前記臨界閾値との間の差分値を重み比率として加算して取得され、分類結果Ｐ（Ｃ）＝ｐ（Ｃｌ）（ｐ（Ｃｌ）＞ｐ（Ｃｒ））、ｐ（Ｃｒ）（ｐ（Ｃｌ）＜ｐ（Ｃｒ））を出力するステップＣ－８と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の検査方法。
【請求項5】
前記分岐方法Ｄは、具体的には、データ方位を統一するステップであって、右眼の４屈折マップデータマトリックスを縦列方向に沿って鏡像し、前記左右眼の４屈折マップデータマトリックスの鼻側及び側頭側の方位を統一するステップＤ－１と、
４屈折マップｄｉｆｆグラフマトリックスを取得するステップであって、それぞれ前記左右眼の４屈折マップデータマトリックスに対してポイントツーポイント差分を行った後に絶対値を取得し、４屈折マップｄｉｆｆグラフデータマトリックスを取得するステップＤ－２と、
ｄｉｆｆデータ特徴を計算するステップであって、直径が６ｍｍのデータ範囲内の前記４屈折マップｄｉｆｆグラフデータマトリックスにおける全てのデータの平均値、最大値、標準偏差をそれぞれ計算して特徴量とするステップＤ－３と、
ステップＤ－３で取得された１２組の前記左右眼角膜の４屈折マップｄｉｆｆグラフの平均値、最大値、標準偏差を特徴とし、全てのサンプルデータを組み合わせて各組のデータの円錐角膜症例と正常症例の臨界閾値を統計するか、又はＳＶＭ分類方法を採用し、全てのタイプのｄｉｆｆデータの特徴に正規化トレーニングとテストを行い、最適な感受性及び特異性を与えるステップＤ－４と、
本分岐方法で取得された前記ある症例の両眼に円錐角膜が発生する確率をそれぞれＰ（Ｄ）として記録するステップであって、前記Ｐ（Ｄ）は、現在の症例に基づいて計算された各組のｄｉｆｆデータ固有値とステップＤ－４で取得された前記臨界閾値との間の差分値を重み比率として加算して取得されるステップＤ－５と、を実現することを特徴とする、
請求項１に記載のマルチモーダルデータに基づく両眼円錐角膜の検査方法。
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】０００１
【補正方法】変更
【補正の内容】
【0001】
本発明は、眼科データ検査技術に関し、特にマルチモーダルデータに基づく両眼円錐角膜検査方法に関する。
【手続補正4】
【補正対象書類名】明細書
【補正対象項目名】０００３
【補正方法】変更
【補正の内容】
【0003】
現在の円錐角膜の検査及び治療は、既に角膜病学、屈折矯正手術学、検眼学等が密接に協力する臨床疾患学に発展している。一般的に角膜トポグラフィーに基づいて円錐角膜検査を行う方法は、臨床統計学的方法であり、角膜トポグラフィーの形態的特徴を利用して臨床パラメータ及び病歴に合わせて円錐角膜の検査及び臨床病期の分類を行うことが多い。統計学のモデルで取得された結果は、サマリーパラメータであり、かつ既知の検査結果データセットに基づいてパラメータ境界線が取得され、例えば、広く応用されるＫＩＳＡ指数、ＩＳ指数、ＳＲＩ／ＳＡＩ指数等である。このような方法は、プラットフォームデータに制限されることが多く、人為的に認められた有限個の単眼特徴に過度に依存し、両眼の間の相互関係を無視し、かつ異なる指数における感受性及び特異性が異なり、早期円錐角膜と不完全型円錐角膜に良好な判断を与えることができない。
【手続補正5】
【補正対象書類名】明細書
【補正対象項目名】０００４
【補正方法】変更
【補正の内容】
【0004】
本発明は、両眼円錐角膜の早期検査の問題に対して、マルチモーダルデータに基づく両眼円錐角膜の検査方法を提供する。マルチモーダル屈折トポグラフィーのマトリックスデータに基づいて、畳み込みニューラルネットワーク法、固有値ＳＶＭ法、両眼対比法及び調整可能なＢＦＳの増幅トポグラフィー法を組み合わせて、両眼の総合円錐角膜の検査結果を与え、特に早期の後円錐角膜及び不完全型円錐角膜のスクリーニング及び検査に対してより高いロバスト性及び正確性を有する。
【手続補正6】
【補正対象書類名】明細書
【補正対象項目名】０００５
【補正方法】変更
【補正の内容】
【0005】
本発明の技術手段は、以下のとおりである。マルチモーダルデータに基づく両眼円錐角膜の検査方法は、具体的には、
両眼マルチモーダルデータを収集するステップであって、データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、
分類された症例に基づいて、両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の４屈折マップにおける各トポグラフィーデータ及び角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得するステップ４）と、
分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含む。
両眼マルチモーダルデータを収集するステップであって、データは、両眼角膜の４屈折マップ及び角膜絶対高さデータを含み、角膜の４屈折マップは、角膜前面の軸方向曲率、角膜前面の相対高さトポグラフィー、角膜後面の相対高さトポグラフィー及び角膜厚さトポグラフィーを含み、角膜絶対高さデータは、角膜前後面の絶対高さデータを含むステップ１）と、
分類された症例に基づいて、両眼マルチモーダルデータを分類された種類と関連付け、要求に応じてデータ分類を行うステップ２）と、
ステップ２）の４屈折マップにおける各トポグラフィーデータ及び角膜前後面の高さデータを同じサイズのデータマトリックスに統一するステップ３）と、
以上のデータに基づいて４つの分岐方法でそれぞれ両眼早期円錐角膜を判断するステップであって、４つの分岐方法は、それぞれ分岐方法Ａ、分岐方法Ｂ、分岐方法Ｃ及び分岐方法Ｄであり、
分岐方法Ａは、全ての４屈折マップデータマトリックスをデータ処理した後に深度畳み込みネットワークの分類ネットワークに送信して円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ａ）を出力し、
分岐方法Ｂは、全ての４屈折マップデータマトリックスにおける各マップデータマトリックスに対して固有値計算を行い、固有値データをＳＶＭサポートベクターマシン二分類方法に送信し、円錐角膜の感受性及び特異性を識別し、ある症例を取得し、当該ある症例に対して分類結果Ｐ（Ｂ）を出力し、
分岐方法Ｃは、角膜前後面の絶対高さデータを最適フィッティング球面データと比較し、円錐角膜症例と正常症例の臨界閾値を取得することにより、ある症例を判定し、当該ある症例に対して分類結果Ｐ（Ｃ）を出力し、
分岐方法Ｄは、左右眼の４屈折マップデータマトリックスの平均値、最大値、標準偏差を特徴量とし、臨界閾値を利用するか、又はＳＶＭ分類方法を用いて、最適な感受性及び特異性、及びある症例の両眼に円錐角膜が発生する確率Ｐ（Ｄ）を取得するステップ４）と、
分岐方法Ａ、Ｂ、Ｃ、Ｄにおける最後の結果を重み付け加算することにより、ある症例の両眼に円錐角膜が発生する最終的な確率を取得するステップ５）と、を含む。
【手続補正7】
【補正対象書類名】明細書
【補正対象項目名】００１０
【補正方法】変更
【補正の内容】
【0010】
本発明の有益な効果は、以下のとおりである。本発明に係るマルチモーダルデータに基づく両眼円錐角膜の検査方法によれば、両眼の間の相互関係を考慮し、深度畳み込みネットワーク方法、従来の機械学習ｓｖｍ方法及び調整可能なＢＳＦの高さマップ増幅方法に基づいた病巣の感受性及び特異性の識別を組み合わせて、その感受性及び特異性のバランスを取り、多次元で患者を単位とする円錐角膜の発病率を総合的に判断する。両眼データに合わせて人工選択特徴を含むだけでなく、深度ネットワーク自体がビッグデータから学習する特徴も含むから、検査方法は、より強いロバスト性及び正確性を有する。
【手続補正8】
【補正対象書類名】明細書
【補正対象項目名】００１１
【補正方法】変更
【補正の内容】
【0011】
【図1】本発明に係るマルチモーダルデータに基づく両眼円錐角膜の検査方法のフローチャートである。
【図2】角膜の４屈折マップである。
【図3】角膜前後面の絶対高さマップである。
【図4】本発明に係る方法における分類ネットワークの構造図である。
【図5】左右眼角膜の４屈折マップＤｉｆｆグラフである。
【手続補正9】
【補正対象書類名】明細書
【補正対象項目名】００１３
【補正方法】変更
【補正の内容】
【0013】
図１に示すマルチモーダルデータに基づく両眼円錐角膜の検査方法のフローチャートは、具体的には、以下のステップ１～５を含む。
【国際調査報告】