問題タブ [fourier-descriptors]

閉じた平面曲線の楕円フーリエ係数 (EFC) を見つけるプログラムを書いていますが、係数の正規化に問題があります。

閉じた平面ポリライン p は、一連の点 m_points によって記述されます。m_points[i][0] は xi 座標を保持し、m_points iは yi 座標を保持します。フォーム 0 から m_num_points-1 まで開始します。

ポリラインの EFC は、そのアルゴリズムによって計算されます (係数は EFD 配列にあります)。

高調波のセットからポリラインを復元するには、そのアルゴリズムを使用します。

ここで、EFD を正規化し、それらを新しい配列 NEFD に配置する必要があります。私はそうします：

半長軸を X に平行にして coefficeint を回転させようとすると、見苦しい結果が得られます。復元された形状が z 軸を中心に回転しているように見えます。(左が元の形状、右が元の形状)

コードの何が問題になっていますか?

アップデート

@MBo の応答が成功した後、次のコードの変更が必要です。

R初心者です。

楕円フーリエ係数 (高調波) の事前計算されたシーケンスを含むテキスト ファイルがあります。最初の列はオブジェクト名、その他 - スペースで区切られた高調波のセット: A0 B0 C0 D0 ... An Bn Cn Dn を保持します。それはファイルの一部です：

read.delim2コマンドを使用してファイルを R にロードします。

次に、 calibration_harmonicpower、plot、PCAなどのMomocs パッケージの関数を使用して、取得したデータを分析したいと思います。

そのためには読み込んだデータを Coe(?) などのオブジェクトに変換する必要があります。（正確にはどちらかわかりません。）

Momocs R パッケージで分析用にロードされたデータを準備するにはどうすればよいですか?

グレースケール画像の直接および逆DFTを取得するためにmatlabでプログラムを実行しようとしていますが、逆を適用した後に元の画像を復元できません。逆出力として複素数を取得しています。私は情報を失っているようなものです。これに関するアイデアはありますか？これが私のコードです：

pmdarima パッケージの autoARIMA を使用して、700 種類の製品を予測しようとしています。季節性については、すべての製品のパターンが異なるため、フーリエ級数が簡単になります。

しかし、ループ内の製品に基づいて、「k」の異なる値を選択するにはどうすればよいですか。そのためのテストまたは最適化機能はありますか?

提案してください。ありがとうございました。