問題タブ [polar-coordinates]

プロジェクトのディスク上のポイントのランダムな位置をマッピングするプログラムを Python で作成しました。ポイントの数が異なると、ランダムな位置が png として保存され、後で使用できるようになります。コードは次のとおりです。

問題は、このコードを実行すると、最初の数字は 3 ポイントで保存されますが、2 番目の数字は 10 ではなく 13 で保存されることです! これは、連続する N について続きます。

画像を投稿することはできないので、私の言葉を信じるか、コードを実行する必要があります。

これを克服する方法を知っている人はいますか？

実行したいシミュレーションを容易にするために必要なデカルト座標 (緯度、経度) を極座標に変換するソリューションを探していましたが、これを行うための質問や回答が見つかりませんでしたR. Matlab の組み込み関数 cart2pol を含む多くのオプションがありますが、私のデータはすべて R にあり、このフレームワークで快適に作業を続けたいと考えています。

質問：

タグ付けデータからの緯度/経度座標があり、これらを極座標 (つまり、ジャンプのサイズと角度: http://en.wikipedia.org/wiki/Polar_coordinate_system ) に変換して、それらをシャッフルまたはブートストラップできるようにしたい (どちらにするかは決めていません）約1,000回、シミュレートされた各トラックの開始点からの直線距離を計算します。私は実際のトラックを持っています。この動物がサイト アフィニティを示しているかどうかを判断することに興味があります。ジャンプ サイズと回転角度は同じですが、順序と組み合わせがまったく異なる 1,000 個のランダムなトラックをシミュレートします。したがって、距離の分布を作成し、これを実際のデータ セットの直線距離と比較するには、原点から 1,000 の直線距離が必要です。

私はブートストラップを快適に行うことができますが、デカルトの緯度/経度座標を極座標 (ジャンプ サイズと回転角度) に変換する最初のステップで立ち往生しています。Matlab などの他のプログラムでこれを行うための関数が組み込まれていることは知っていますが、R でそれを行う方法が見つかりません。for ループで手動で行うこともできますが、パッケージが出ている場合そこに、またはそれを行う簡単な方法があれば、私はそれを望んでいます。

理想的には、データを極座標に変換してシミュレーションを実行し、各ランダム トラックの終点をデカルト座標 (緯度/経度) として出力して、移動した直線距離を計算できるようにしたいと考えています。

緯度と経度の座標の 2 列のデータ フレームになるだけなので、サンプル データは投稿しませんでした。

ご協力いただきありがとうございます。このサイトや私が見逃した他のサイトのどこかに簡単な説明がある場合は、その方向を教えてください! 何も見つかりませんでした。

乾杯

ラジアル軸の開始をオフセットしたり、グラフの外に移動したりすることは可能ですか?

これは私が達成したいと思っていることです:

そして、これが私が今持っているものです。

SOに関するドキュメントとさまざまなトピックを読みましたが、役立つものは何も見つかりませんでした。それは、どこにも言及されていなければ、それも不可能だということですか。

前もって感謝します。

編集（プロットの作成に使用されるコードのスニペットを追加）：

現在、次の計算を使用して、想定される斜辺 150 と既知のヨー ピッチとロールで x、y、z 座標を取得しています。

視聴者が 0,0,0 の位置に立って P を見上げていると仮定します。Q は 150 単位離れており、ビューのヨーとピッチを知っています)。

私の数学は、ピッチがまっすぐ上下に近づくまでうまくいくようです。その時点で、x と y は何らかの方法で z を考慮する必要があることに気付きます..助けてください

1. やりたいこと:
(i)入力 n を使用して n*n デカルト グリッドを生成する

(ii)極座標の生成

(iii)円筒座標で関数を評価する

（iv）（たとえば）pcolor（またはsurfなど）を使用して結果をプロットします

2. できること... プロットを機能させる唯一の方法は、極パラメータから始めて、そこからデカルトに戻ることです。
(i)パラメータを定義する

(ii)両方のグリッドを作成し、f を評価する

(iii)プロット

問題は、現在私のグリッドが円形であり、四角いグリッドのどこでも関数を評価したいということです (私の分析は正方形のピクセル配列に依存しているため)。繰り返しますが、上記の方法 2 を使用すると、正方形に囲まれた円形のプロットが得られます。エッジに沿って白の黒い円を想像してください...しかし、この「白い」領域で関数を評価したいと思います。ただし、方法1を使用しても機能しません-プロットすると関数がすべて台無しになります(プロットがどのように見えるかを確認するには、ラゲールガウスモードをグーグルで検索してください)。

極座標から始めてすべてのデカルト点を割り当てるのではなく、長方形グリッドから始めてすべての点に極座標を割り当てられるようにしたいと考えています。

私はこれを長い間いじっていましたが、この一見単純な問題を回避する方法がわかりません。

編集 1
問題は、座標行列の生成方法にあるようです。以下に、アプローチ 1 とアプローチ 2 が異なる数値を生成する方法を示す単純な 3by3 の例のスクリーンショットを投稿しました。

これらの数値に互換性を持たせるにはどうすればよいですか?

評判ポイントがないので、画像を直接アップロードすることはできません...以下のリンクは3by3の例を示しています...作成しようとしているラゲールガウスプロットの実際の画像へのリンクについては、コメントを参照してください...

cart2pol
を適用する pol2cart を適用する

編集 2

現在、アプローチの結果 ( 1. ) は、次に示すように間違った結果をもたらします。

望ましいアプローチ、間違った結果

2 番目のアプローチでは適切な画像が得られますが、残念ながらこれは円であり、正方形全体ではありません。ここに示されています：

実施されたアプローチ、限られた結果

両方のアプローチの 3D プロットをここに示します。上の図のカラフルな部分だけが正しいです。

編集 3

これは、上記で使用されている関数fのスクリーンショットです。だけでなく、より多くの入力パラメーターを要求することに注意してくださいr,theta。一般的な値は次のとおりです。

この関数Cは正規化を行いL、ラゲール多項式です。これらの機能は両方とも徹底的にテストされており、期待どおりの結果が得られています。

編集 4ここに私の例を明示的
に実行するのに十分なコードがあります。z=U(0,5,r,phi,w0)+U(0,-5,r,phi,w0);プロット自体は で与えられpcolor(x,y,abs(z).^2)ます。

Lpl()関数はコメントとして挿入されることに注意してください。U 関数を適切に実行するには、これを独自の m ファイルとして保存する必要があります。

R の画像 (行列として表される) を、原点が 0,0 (左上隅) の極座標空間に変換しようとしています。次のような215x215行列が与えられます。x

ただし、 xn は次のようになります。

私がこれを期待するとき：

私が間違っていることは何ですか？

matplotlib次のような極座標プロットを使用してパイをプロットしようとしています:

math.pi/3.0との間のグラフの部分だけを見たいです2*math.pi/3.0。これは可能ですか？