【手続補正書】
【提出日】令和3年9月7日(2021.9.7)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
１つ又は複数のハードウェアプロセッサによって、ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割するステップであって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割するステップと；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉで：
１つ又は複数のハードウェアプロセッサによって、前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算するステップであって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表し、
2020533591.tif 000001
であり、
＾Ｓ（ベクトルｋ，ｔｉ）は、前記ソース波動場の前記フーリエ変換である、前記計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記ソース波動場の時間微分のフーリエ変換の正の成分である、
2020533591.tif 000002
を計算するステップであって、＾Ｓ´（ベクトルｋ，ｔｉ）は、前記ソース波動場の時間微分の前記フーリエ変換である、前記計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記ソース波動場の前記フーリエ変換の前記正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算するステップと；
１つ又は複数のハードウェアプロセッサによって、特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算するステップであって、
2020533591.tif 000003
であり、
＾Ｒ（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の前記フーリエ変換である、前記計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記レシーバ波動場の時間微分のフーリエ変換の負の成分である、
2020533591.tif 000004
を計算するステップであって、＾Ｒ´（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の時間微分の前記フーリエ変換である、前記計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記レシーバ波動場の前記フーリエ変換の前記負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の前記下降成分を合計するステップにより、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算するステップと；
１つ又は複数のハードウェアプロセッサによって、前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の前記上昇成分を合計するステップにより、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算するステップと；
１つ又は複数のハードウェアプロセッサによって、リバースタイムマイグレーション（ＲＴＭ）で偽像を除去するために、前記計算されたＳｄ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）に基づいて画像を生成するステップと；を備える、
方法。
【請求項２】
2020533591.tif 000005
であり、
Ｖは、前記特定のウィンドウにおける平均速度であり、
2020533591.tif 000006
は、前記波数ベクトルのノルムであり、
εは、ユーザが指定した数であり、
ｊは、虚数単位である、
請求項１に記載の方法。
【請求項３】
2020533591.tif 000007
である、
請求項２に記載の方法。
【請求項４】
前記複数のウィンドウは、複数の重なり合うウィンドウである、
請求項１に記載の方法。
【請求項５】
各ウィンドウにおける、前記ソース波動場の前記下降成分を合計するステップと、前記レシーバ波動場の前記上昇成分を合計するステップは、各ウィンドウｗにおけるＳＤ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）を、重なり合う領域においてテーパをつけて合計するステップ；を備える、
請求項４に記載の方法。
【請求項６】
前記フーリエ変換は、１Ｄ空間のフーリエ変換である、
請求項１に記載の方法。
【請求項７】
前記Ｓｄ（ベクトルｘ，ｔｉ）と前記Ｒｕ（ベクトルｘ，ｔｉ）とに基づき画像項Ｉｄｕ（ベクトルｘ）を計算するステップ；を更に備える、
請求項１に記載の方法。
【請求項８】
メモリと；
操作を実行するために編成された１つ又は複数のハードウェアプロセッサと；を備え、
前記操作は：
ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割する操作であって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割する操作と；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉでの：
前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算する操作であって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表し、
2020533591.tif 000008
であり、
＾Ｓ（ベクトルｋ，ｔｉ）は、前記ソース波動場の前記フーリエ変換である、前記計算する操作と；
前記ソース波動場の時間微分のフーリエ変換の正の成分である、
2020533591.tif 000009
を計算する操作であって、＾Ｓ´（ベクトルｋ，ｔｉ）は、前記ソース波動場の時間微分の前記フーリエ変換である、前記計算する操作と；
前記ソース波動場の前記フーリエ変換の前記正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算する操作と；
特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算する操作であって、
2020533591.tif 000010
であり、
＾Ｒ（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の前記フーリエ変換である、前記計算する操作と；
前記レシーバ波動場の時間微分のフーリエ変換の負の成分である、
2020533591.tif 000011
を計算する操作であって、＾Ｒ´（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の時間微分の前記フーリエ変換である、前記計算する操作と；
前記レシーバ波動場の前記フーリエ変換の前記負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の前記下降成分を合計する操作により、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の前記上昇成分を合計する操作により、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算する操作と；
リバースタイムマイグレーション（ＲＴＭ）で偽像を除去するために、前記計算されたＳｄ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）に基づいて画像を生成する操作と；を含む、
デバイス。
【請求項９】
2020533591.tif 000012
であり、
Ｖは、前記特定のウィンドウにおける平均速度であり、
2020533591.tif 000013
は、前記波数ベクトルのノルムであり、
εは、ユーザが指定した数であり、
ｊは、虚数単位である、
請求項８に記載のデバイス。
【請求項１０】
2020533591.tif 000014
である、
請求項９に記載のデバイス。
【請求項１１】
前記複数のウィンドウは、複数の重なり合うウィンドウである、
請求項８に記載のデバイス。
【請求項１２】
各ウィンドウにおける、前記ソース波動場の前記下降成分を合計する操作と、前記レシーバ波動場の前記上昇成分を合計する操作は、各ウィンドウｗにおけるＳＤ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）を、重なり合う領域においてテーパをつけて合計する操作；を備える、
請求項１１に記載のデバイス。
【請求項１３】
前記フーリエ変換は、１Ｄ空間のフーリエ変換である、
請求項８に記載のデバイス。
【請求項１４】
前記操作は、
前記Ｓｄ（ベクトルｘ，ｔｉ）と前記Ｒｕ（ベクトルｘ，ｔｉ）とに基づき画像項Ｉｄｕ（ベクトルｘ）を計算する操作；を更に備える、
請求項８に記載のデバイス。
【請求項１５】
操作を実行するために、コンピュータシステムによって実行可能な命令を格納する、非一時的なコンピュータ読取り可能媒体であって：
前記操作は：
ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割する操作であって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割する操作と；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉでの：
前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算する操作であって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表し、
2020533591.tif 000015
であり、
＾Ｓ（ベクトルｋ，ｔｉ）は、前記ソース波動場の前記フーリエ変換である、前記計算する操作と；
前記ソース波動場の時間微分のフーリエ変換の正の成分である、
2020533591.tif 000016
を計算する操作であって、＾Ｓ´（ベクトルｋ，ｔｉ）は、前記ソース波動場の時間微分の前記フーリエ変換である、前記計算する操作と；
前記ソース波動場の前記フーリエ変換の前記正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算する操作と；
特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算する操作であって、
2020533591.tif 000017
であり、
＾Ｒ（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の前記フーリエ変換である、前記計算する操作と；
前記レシーバ波動場の時間微分のフーリエ変換の負の成分である、
2020533591.tif 000018
を計算する操作であって、＾Ｒ´（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の時間微分の前記フーリエ変換である、前記計算する操作と；
前記レシーバ波動場の前記フーリエ変換の前記負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の前記下降成分を合計する操作により、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の前記上昇成分を合計する操作により、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算する操作と；
リバースタイムマイグレーション（ＲＴＭ）で偽像を除去するために、前記計算されたＳｄ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）に基づいて画像を生成する操作と；を備える、
非一時的なコンピュータ読取り可能媒体。
【請求項１６】
前記複数のウィンドウは、複数の重なり合うウィンドウである、
請求項１５に記載の媒体。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００７１
【補正方法】変更
【補正の内容】
【００７１】
Ｔｏｎｇ  Ｗ．  Ｆｅｉ，  Ｙｉ  Ｌｕｏ  ａｎｄ  Ｆｕｈａｏ  Ｑｉｎ，  ２０１４，  Ａｎ  ｅｎｄｅｍｉｃ  ｐｒｏｂｌｅｍ  ｉｎ  ｒｅｖｅｒｓｅ−ｔｉｍｅ  ｍｉｇｒａｔｉｏｎ：  ８４ｔｈ  Ａｎｎｕａｌ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｍｅｅｔｉｎｇ，  ＳＥＧ  Ｅｘｐａｎｄｅｄ  Ａｂｓｔｒａｃｔｓ．
［第１の局面］
ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割するステップであって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割するステップと；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉでの：
前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算するステップであって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表す、前記計算するステップと；
前記ソース波動場のフーリエ変換の正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算するステップと；
特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算するステップと；
前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算するステップと；
前記レシーバ波動場のフーリエ変換の負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算するステップと；
前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算するステップと；
前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の下降成分を合計することにより、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算するステップと；
前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の上昇成分を合計することにより、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算するステップと；を備える、
方法。
［第２の局面］
2020533591.tif 000019
であり、
＾Ｓ（ベクトルｋ，ｔｉ）は、前記ソース波動場のフーリエ変換であり、
＾Ｒ（ベクトルｋ，ｔｉ）は、前記レシーバ波動場のフーリエ変換である、
第１の局面に記載の方法。
［第３の局面］
複数の重なり合うウィンドウの各ウィンドウについて、前記特定の時間値ｔ＝ｔｉでの：
前記ソース波動場の時間微分のフーリエ変換の正の成分である、
2020533591.tif 000020
を計算するステップであって、＾Ｓ´（ベクトルｋ，ｔｉ）は、前記ソース波動場の時間微分のフーリエ変換である、前記計算するステップと；
前記レシーバ波動場の時間微分のフーリエ変換の負の成分である、
2020533591.tif 000021
を計算するステップであって、＾Ｒ´（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の時間微分のフーリエ変換である、前記計算するステップと；を更に備える、
第２の局面に記載の方法。
［第４の局面］
2020533591.tif 000022
であり、
2020533591.tif 000023
Ｖは、前記特定のウィンドウにおける平均速度であり、
εは、ユーザが指定した数であり、
ｊは、虚数単位である、
第３の局面に記載の方法。
［第５の局面］
2020533591.tif 000024
である、
第４の局面に記載の方法。
［第６の局面］
前記複数のウィンドウは、複数の重なり合うウィンドウである、
第１の局面に記載の方法。
［第７の局面］
各ウィンドウにおける、前記ソース波動場の下降成分を合計するステップと、前記レシーバ波動場の上昇成分を合計するステップは、各ウィンドウｗにおけるＳＤ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）を、重なり合う領域においてテーパをつけて合計するステップ；を備える、
第６の局面に記載の方法。
［第８の局面］
前記フーリエ変換は、１Ｄ空間のフーリエ変換である、
第１の局面に記載の方法。
［第９の局面］
前記Ｓｄ（ベクトルｘ，ｔｉ）と前記Ｒｕ（ベクトルｘ，ｔｉ）とに基づき画像項Ｉｄｕ（ベクトルｘ）を計算するステップ；を更に備える、
第１の局面に記載の方法。
［第１０の局面］
メモリと；
操作を実行するために編成された処理ユニットと；を備え、
前記操作は：
ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割する操作であって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割する操作と；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉでの：
前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算する操作であって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表す、前記計算する操作と；
前記ソース波動場のフーリエ変換の正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算する操作と；
特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の下降成分を合計することにより、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の上昇成分を合計することにより、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算する操作と；を含む、
デバイス。
［第１１の局面］
2020533591.tif 000025
であり、
＾Ｓ（ベクトルｋ，ｔｉ）は、前記ソース波動場のフーリエ変換であり、
＾Ｒ（ベクトルｋ，ｔｉ）は、前記レシーバ波動場のフーリエ変換である、
第１０の局面に記載のデバイス。
［第１２の局面］
前記操作は、
複数の重なり合うウィンドウの各ウィンドウについて、前記特定の時間値ｔ＝ｔｉでの：
前記ソース波動場の時間微分のフーリエ変換の正の成分である、
2020533591.tif 000026
を計算する操作であって、＾Ｓ´（ベクトルｋ，ｔｉ）は、前記ソース波動場の時間微分のフーリエ変換である、前記計算する操作と；
前記レシーバ波動場の時間微分のフーリエ変換の負の成分である、
2020533591.tif 000027
を計算する操作であって、＾Ｒ´（ベクトルｋ，ｔｉ）は、前記レシーバ波動場の時間微分のフーリエ変換である、前記計算する操作と；を更に備える、
第１１の局面に記載のデバイス。
［第１３の局面］
2020533591.tif 000028
であり、
2020533591.tif 000029
Ｖは、前記特定のウィンドウにおける平均速度であり、
εは、ユーザが指定した数であり、
ｊは、虚数単位である、
第１２の局面に記載のデバイス。
［第１４の局面］
2020533591.tif 000030
である、
第１３の局面に記載のデバイス。
［第１５の局面］
前記複数のウィンドウは、複数の重なり合うウィンドウである、
第１０の局面に記載のデバイス。
［第１６の局面］
各ウィンドウにおける、前記ソース波動場の下降成分を合計する操作と、前記レシーバ波動場の上昇成分を合計する操作は、各ウィンドウｗにおけるＳＤ（ベクトルｘ，ｔｉ）及びＲＵ（ベクトルｘ，ｔｉ）を、重なり合う領域においてテーパをつけて合計する操作；を備える、
第１５の局面に記載のデバイス。
［第１７の局面］
前記フーリエ変換は、１Ｄ空間のフーリエ変換である、
第１０の局面に記載のデバイス。
［第１８の局面］
前記操作は、
前記Ｓｄ（ベクトルｘ，ｔｉ）と前記Ｒｕ（ベクトルｘ，ｔｉ）とに基づき画像項Ｉｄｕ（ベクトルｘ）を計算する操作；を更に備える、
第１０の局面に記載のデバイス。
［第１９の局面］
操作を実行するために、コンピュータシステムによって実行可能な命令を格納する、非一時的なコンピュータ読取り可能媒体であって：
前記操作は：
ソース波動場Ｓ（ベクトルｘ，ｔ）及びレシーバ波動場Ｒ（ベクトルｘ，ｔ）のボリュームを複数のウィンドウに分割する操作であって、ベクトルｘ＝（ｘ，ｙ，ｚ）は位置ベクトルを表し、ｔは時間を表す、前記分割する操作と；
前記複数のウィンドウの各ウィンドウについて、特定の時間値ｔ＝ｔｉでの：
前記ソース波動場のフーリエ変換の正の成分である＾Ｓ＋（ベクトルｋ，ｔｉ）を計算する操作であって、ベクトルｋ＝（ｋｘ，ｋｙ，ｋｚ）は波数ベクトルを表す、前記計算する操作と；
前記ソース波動場のフーリエ変換の正の成分のヒルベルト変換である＾ＨＳ＋（ベクトルｋ，ｔｉ）を計算する操作と；
特定のウィンドウにおける前記ソース波動場の下降成分であるＳＤ（ベクトルｘ，ｔｉ）を、前記＾Ｓ＋（ベクトルｋ，ｔｉ）と前記＾ＨＳ＋（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分である＾Ｒ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記レシーバ波動場のフーリエ変換の負の成分のヒルベルト変換である＾ＨＲ−（ベクトルｋ，ｔｉ）を計算する操作と；
前記特定のウィンドウにおける前記レシーバ波動場の上昇成分であるＲＵ（ベクトルｘ，ｔｉ）を、前記＾Ｒ−（ベクトルｋ，ｔｉ）と前記＾ＨＲ−（ベクトルｋ，ｔｉ）の和の逆フーリエ変換の実部成分として計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記ソース波動場の下降成分を合計することにより、前記特定の時間値での前記ソース波動場の下降成分であるＳｄ（ベクトルｘ，ｔｉ）を計算する操作と；
前記複数のウィンドウの各ウィンドウにおける前記レシーバ波動場の上昇成分を合計することにより、前記特定の時間値での前記レシーバ波動場の上昇成分であるＲｕ（ベクトルｘ，ｔｉ）を計算する操作と；を備える、
非一時的なコンピュータ読取り可能媒体。
［第２０の局面］
前記複数のウィンドウは、複数の重なり合うウィンドウである、
第１９の局面に記載の媒体。