Matlab に画像があり、その中のすべての重心を見つけています。画像には複数の重心があり、「imellipse」を使用してそのうちの 1 つを隔離しています。重心の1つのピクセル位置を見つけようとしています。Matlab のコマンド ウィンドウに場所のポップアップを表示するだけで十分です。これが私のコードです：

これは、重心マーカーのピクセル位置を推定しようとしている写真です。 http://imgur.com/DnSQTUi

Matlab で画像のセクションの周りにボックスを描画し、その画像の重心を見つける方法を見つけようとしています。「新しい」セクションを元の画像の隣に表示することが理想的です。その後、その「新しい」画像の重心に移動できます。

これが私のコードです：

これは、マスクを重ねてそのセクションの重心を見つけようとしている画像です。

私は生態学者で、主にビーガン R パッケージを使用しています。

私は2つの行列（サンプルx存在量）を持っています（以下のデータを参照）：

マトリックス 1/ nrow= 6replicates*24sites、ncol=15 種の存在量 (魚) マトリックス 2/ nrow= 3replicates*24sites、ncol=10 種の存在量 (無脊椎動物)

サイトは両方のマトリックスで同じです。サイトのペア間で全体的なブレイ・カーティスの非類似度 (両方の行列を考慮) を取得したいと考えています。2 つのオプションが表示されます。

オプション 1、レプリケート (サイト スケールで) の魚類と大型無脊椎動物の存在量を平均し、2 つの平均存在量行列 (nrow=24sites、ncol=15+10 平均存在量) をバインドし、bray-curtis を計算します。

オプション 2、集合体ごとに、サイトのペア間のブレイカーティス非類似度を計算し、サイト重心間の距離を計算します。次に、2 つの距離行列を合計します。

よくわからない場合は、以下の R コードでこれら 2 つの操作を行いました。

選択肢 2 が正しく、選択肢 1 よりも適切かどうか教えてください。

前もって感謝します。

ピエール

ここにRコードの例があります

データの生成

オプション 1、存在量と cbind の平均化

オプション 2、重心間の集合距離ごとに計算し、2 つの距離行列を合計する

Gavin Simpson の fuse() を使用して 2 つの距離行列を合計する

2 つのユークリッド距離行列の合計 (Jari Oksanen 補正のおかげ)

さらに距離ベースの分析を行うには、以下の「coord.centroid」を使用します (正しいですか?)

オプション 1 と 2 の比較

私はプログラミングに少し慣れていないので、すみません。基本的に、私は「分析」し、ユーザーが 2 辺の長さとそれらの間の角度のサイズを指定する三角形の画像を生成するように割り当てられ、プログラムが実行され、3 番目の辺の長さと 2 辺の長さを見つけます他の角度 (余弦の法則を使用)。次に、すべての辺の長さ、すべての角度測定値 (オプション) のテキスト出力が必要です。また、画像に「ここにあなたの三角形があります\n x 平方ピクセルの領域があります」と出力するタートルの形式で領域を出力する必要があります。同じように。また、三角形の重心は (0,0) でなければなりません。これは私がこれまでに持っているものです:

また、これにもタートル グラフィックを使用しています。

私はこれに本当に苦労しています。いくつかの助けをいただければ幸いです。助けてください！ありがとう！

問題は、バイナリ イメージ内のいくつかの白い輪郭の重心を計算できる単純化されたバージョンのアルゴリズムが必要なことです。たとえば、白い輪郭が 1 つしかない場合、輪郭の中心の座標 Xc と Yc は次の式を使用して計算されます。

ここで、M は強度 m_i の合計、m_i はピクセル強度値、x_i と y_i は画像上のピクセル位置、n はピクセルの総数です。

誰かがいくつかの輪郭で作業するのに似たものを提案できますか、または同じ式を使用してそれらの1つを計算するときに他のものを無視する方法はありますか?

PathIterator を使用して任意の Shape オブジェクトの中心を計算し、曲線のパスを考慮できるようにしようとしていますが、標準の 1x1 長方形の中心を見つけると、getCenter() メソッドがポイントを返します。

私の getCenter() メソッド:

points.toString() を印刷すると、コンソールに次のように表示されることがわかりました。

入力 Shape オブジェクトが Rectangle2D.Double(0, 0, 1, 1) であることを考えると、points 配列には 4 つではなく、6 つのエントリがあることに気付きました。明らかに、それはポイント (0, 0) を私が望むよりも 2 回多く説明していて、それがなぜなのか混乱しています。PathIterator.isDone() メソッドの結果ですか? 私はそれを間違って使用していますか？PathIterator ができない場合、どうすれば問題を解決できますか?

このような点のプロットがあります。

これらの点が形成するトラックは、円または楕円になります。上の 2 つの画像の円形トラックの中心が明らかに異なっています。

これらのトラック (円形/楕円形) の中心点を見つけるにはどうすればよいですか? 中心である (x,y) 座標を見つけたいのですが、プロットされたデータ セット内の点である必要はありません。つまり、medoid は必要ありません。

編集:また、これらの点の大部分を包む円/楕円の方程式を見つけることができる方法はありますか? 楕円トラックでは、トラック上のポイントを囲む楕円を追加しました。値は試行錯誤によって計算されました。中心も、プロットを目で見て計算されました。プログラムでこれを行うにはどうすればよいですか?

kmeans アルゴリズムの使用例を次に示します: http://mnemstudio.org/clustering-k-means-example-1.htm

この例では、著者は最初の重心として「(ユークリッド距離測定を使用して) 最も離れた 2 人の個人」を使用しました。

2 つのクラスターではなく、10 のクラスターが必要な場合はどうすればよいでしょうか。最初の 10 個の重心を選択するにはどうすればよいですか? 最も離れた 10 人の個人を選択する方法はありますか? または、別の方法でそれらを選択する必要があります。

PS:私の場合、ランダムな選択を使用するのは良くないと思います。また、最初の 10 人を重心として使用しようとしてきましたが、それらを選択するより良い方法を探しています。