ECG画像から信号を復元する方法

私のプロジェクトでは、通常のカメラ（jpeg）で撮影したECG画像をデジタル化する必要があります。たとえば、次のカメラキャプチャ画像があります。

そして、私はこのようなものを取得したい：-

そして、ECGのデジタル化に関するこのビデオのように、デジタル化されたデータ（x、yポイント）

どうすればいいかわからないので、いくつかの研究論文を検索して相談しました。アルゴリズムの一般的なアプローチは次のとおりです。

1. グレーレベル画像に変更
2. グリッド線を削除する
3. 欠落点を追加する
4. 2D画像を1D画像に変換する

2番目のポイント、つまり、グリッド線の削除に固執しています。これを行うためにいくつかの参考文献を調べたところ、ヒストグラム分析が役立つ可能性があります。

これを行う方法を教えてください（MATLAB 2010を使用しています）？任意の助けをいただければ幸いです。

申し訳ありませんが、私はMathematicaを使用していますが、このアイデアをMatlabに簡単に実装できるはずです。とにかくコードを提供するので、説明が十分に詳しくない場合は、コードから残りを取得できます。

基本的な考え方は、画像を列ごとに見ることです。ピクセルのすべての列を個別に処理します。プロットでは、グレー値を反転していることに注意してください。したがって、黒は1、白は0
です。（反転した）輝度ピクセル値をプロットすると、基本的に2つの状況しかありません。1つ目は、列が垂直グリッド線上にない場合です。プロットは次のようになります

2番目の状況は、垂直グリッド線上に直接いる場合です。次に、グリッド線は列全体の明るさに影響します

しかし、あなたが見るのは、あなたの暗いEEGが常に最大であるように見えるということです。したがって、非常に複雑なアルゴリズムは次のとおりです。すべての列を調べて、最も黒いピクセルの位置を取得します。

img = ColorConvert[
   ImagePad[
    Import["http://i.stack.imgur.com/500Kg.jpg"], {{0, -20}, {0, 0}}],
    "Grayscale"];

Image[
 Transpose[Function[With[{m = Min[#]},
     Map[Function[{v}, If[v == m, 1, 0]], #]]] /@ 
   Transpose[ImageData[img, "Real"]]
  ]
 ]

ここでは、画像が完全に白だったので、画像の右側を少しトリミングしました。結果は

これで、ポイントに参加したり、好きな方法で補間したりして、脳波を取得できます

グリッド線が赤でトレースが黒のカラー画像があります。そのため、赤のピクセルは無視してください！

画像が正確に配置されているかどうかわからない場合は、グリッド線を使用してスキューを計算できます（右に進むと、単純に画像の傾斜がピクセル/ピクセルで表されます）。

その後、黒いトレースを1d値に変えるのは簡単です。左端の最初の列から始めて、黒いピクセル（またはピクセルの小さな連結グループの重心）を見つけます-垂直位置があなたの値です。
それを画像の各列に対して行います。
欠損値がある場合、既知の値の前後で補間する必要があります。

余分なポイントについては、ランダムなドットまたはノイズスパイクを見つけることができるように、列と列の間でトレースが変化する量に制限を設定できます。