基本周波数とそのすべての高調波を自動的に識別して削除するにはどうすればよいですか?
削除したいハートビートアーティファクトで汚染された顕微鏡データがあります。データは、約60Hzでキャプチャされた大量の時系列画像で構成されています。 GIF形式の小さなサンプルクリップを次に示します。 時間の経過に伴う平均ピクセル強度を取得し、ウェルチ法を使用してピリオドグラムを計算しました。 ご覧のように、〜1.8Hzに鋭いピークがあり、心拍数(〜108ビート/分)に対応する可能性があります。また、1.8Hzの整数倍に多数の調和ピークがあります。正確な心拍数はデータセットごとに異なる可能性がありますが、ピリオドグラムの影付きの領域で示されているように、生物学的にもっともらしい範囲を指定できます。 私ができるようにしたいのは: 心拍に対応する基本周波数とそのすべての高調波を自動的に検出します 基本波とすべての高調波を除去するようにデータをフィルタリングします。 現時点では、ピリオドグラムで最大のピークを見つけ、それをで乗算することにより、ポイント1を非常に大まかに解くことができは、調和ピークの推定数ですが、このハックよりも良い方法である。1 、2 、。。。、N1,2,...,N1, 2, ..., NNNN ポイント2に関しては、コムフィルターを使用して基本波とそのすべての高調波を削除することについて言及しているこの質問に出くわしました。これは使用するのに最適な方法ですか?重要な考慮事項の1つは、大きな配列の各ピクセルの時系列にフィルターを適用する必要があるため、計算効率の高い方法が非常に望ましいことです。 データの例 ここに、.csvフォーマットの平均ピクセル値の例の1D時系列をアップロードしました。 ここでは、圧縮.matファイルに1000フレームのサンプル(空間的に2:1に間引かれた)も追加しました