問題タブ [spectral]

私は現在、スペクトルクラスタリングを学んでいます。

L = D - W で計算されたラプラシアン行列を分解します。

W は隣接行列です。

ただし、スペクトル クラスタリングのような多くのコードをオンラインで見つけました。

diag(sum(W)) によって D を直接計算します。

D は次数行列である必要があることを知っています。これは、対角線上の各値が各点の次数であることを意味します。

しかし、W が加重グラフの場合、diag(sum(W)) は実際の「次数行列」と等しくありません...

なぜ彼らはまだこれを行うのですか。

私は次のモデルを持っています

通常、このモデルは 4 つのピークを与えるはずですが、ピリオドグラムを使用した場合

しかし、私は次の写真を取得しています

私が得ることができない問題は何ですか、私を助けてください

3000x50 の特徴ベクトル マトリックスがあります。sklearn.metrics.pairwise_distances「Similarity_Matrix」として使用して、これの類似度マトリックスを取得しました。ここnetworkxで、前のステップで生成された類似度マトリックスを として使用してグラフを作成していましたG=nx.from_numpy_matrix(Similarity_Matrix)。今、このグラフでスペクトル クラスタリングを実行したいのですGが、いくつかの Google 検索では、このグラフで scikit Learn スペクトル クラスタリングの適切な例を提供できませんでした :(少なくとも私のような初心者には。

これ、またはnetworkx、scikit Learnなどを使用したグラフカットまたはグラフパーティショニングのコードサンプルを誰かに教えてもらえますか.

どうもありがとう！

現在、.bil 形式 (ESRI BIL だと思います) で提供されているPRISM Climate Groupからの降水データがあり、これらのデータセットを Python で読み取れるようにしたいと考えています。

スペクトルパッケージをインストールしましたが、open_image()メソッドはエラーを返します。

IOError: Unable to determine file type or type not supported.

スペクトルのドキュメントには、BIL ファイルをサポートしていると明確に記載されています。同様/同等のESRI EHdr形式をサポートしていると思われるGDALの使用にもオープンですが、開始するための適切なコードスニペットが見つかりません。

球面調和関数を立方デカルト グリッドに補間しようとしています。

球状の疑似スペクトル シミュレーションの出力データには、とNrの間の動径レベルがrMinありrMax、それぞれに経度と緯度の一連の有限次球面調和関数が含まれています。球面調和関数は、三角形の切り捨てによってNi緯度と経度を含む物理的な球面グリッドにマッピングされます。Nj

ドメインは次のとおりです。

• ラジアル レベル: rMin <= r(k) <= rMax、インデックス付き1 <= k <= Nr
• 球面調和関数 (三角形の切り捨て、変換によるエイリアシングなし):
  • Nm = (Nj-1)/3
  • 0 <= m <= Nm
  • m <= l <= Nm
  • nlm == (nm+1)*(nm+2)/2l(のm組み合わせの総数)

データ配列:

• スペクトル形式:complex*16, dimension( 1:nlm, 1:Nr ) :: foo_spectral
• デカルト形式:real*8, dimension( 1:Nx, 1:Ny, 1:Nz ) :: foo_cartesian

2*rMax球状ドメインが完全に内側に収まるように、データをスペクトル表現から edge-length の立方デカルト グリッドに補間する正確かつ効率的な方法を探しています。ただし、球内で補間したいだけです。r<rMinまたはに対応する点についてrMax<rは、立方体グリッドにOUTSIDE_DOMAIN値が必要です。

現在、データをスペクトル表現 (球面ハーモニクス: foo(Nr,nlm)) から物理表現 (球面グリッド: ) に変換し、IDL のQHULLfoo(Nr,Ni,Nj)ルーチンを使用して、物理的な球面グリッドから物理的な立方体グリッド ( ) (立方体グリッドの場合に注意してください)。foo(Nx,Ny,Nz)Nx==Ny==Nz

私のデータのサイズは既存のコード (IDL で記述) が処理できるサイズよりも大きく、私の目的には球状空間への変換は不要です。たとえば、IDL に依存しない、スタンドアロンのより直接的な方法が必要です。

これをどのように行うことができるかについて何か考えはありますか? 私はオープンソースのライブラリを使用したいと思っていますが、そうしなくてもよかったと思っています。

前もって感謝します！

WAV ファイルのスペクトル重心 (SC) を検査したいと考えています。

そのために、次の MATLAB コードを使用しています。

「SpectralCentroid2(366383, 1024, 128, 44100)」と入力すると、MATLAB は次のように表示します。

なぜこの問題が発生したのかわかりません。