2011-12-21 15 views
0

私は、網膜の密度セルに異なる密度点を使って地形図を作ろうとしています。 volcano3dの例は本当にいいですし、私はこの例に基づいて私の言いたいことを考えています。 私は島の輪郭にx、y座標がgeom_pathのようにプロットされていればパス内にのみ等高線マップを構築できますか? 私はそれをやろうとしましたが、問題はx軸とy軸が連続していないことです。 これは私のデータの最初のポイントで、例を挙げると、xとyは座標で、freqはz軸でなければなりません。データ全体を入れないのは申し訳ありませんが、輪郭は3000ポイント以上あり、内側には約250個のサンプリングがありますので、長いデータセットです。ggplot2 geom_contour離散的オブジェクト

p <- ggplot(data = contour.xy, aes(x = x, y = y)) + 
    geom_path(data = contour.xy, aes(x = x, y = y)) + 
    geom_path(data = opticnerve.xy, aes(x = x, y = y)) + 
    geom_contour(data = counts, aes(x = "x", y = "y", z = "freq")) + 
    coord_equal() 

私は画像を入れてみましたが、アンチスパム規制は、私は写真を添付することはできませんので、私はので、新しいユーザーです。

は、私は本当にここでvolcano3d例です

おかげで、メリークリスマス、すべての

+0

のためにエイドリアン・バッデリーにすべてのあなたの助けが、特別には感謝して、私はすべてのポイントを作成することができます私たちは、あなたがこれまで持っているもの、あなたのコードを表示します。 –

+0

ggplotは+ geom_path(データ= contour.xy、AES(X = X、Yは= Y))+ geom_path(データ= opticnerve(データ= contour.xy、AES(のx = x、yはy)を=)。+ geom_contour(data = counts、aes(x = "x"、y = "y"、z = "freq"))+ coord_equal() –

+0

は次のとおりです。私は画像をプロットする必要があるコードは、輪郭を適用しようとするまでうまく動作します –

答えて

1

は、最後にいくつかの研究を行った後、私は私が探していた正確に何ができますパッケージを見つけました。それはAdrian Baddeleyによって書かれたspatstatであり、彼は非常に親切にスクリプトの一部で私を助けます。基本的には、不規則な形のウィンドウを作成することです。この場合、私の「島」と、分析する必要があるすべての点がウィンドウの内側にあります。そしてowin(ポリ=リスト(X = XP、Y = YP)) プロット(マップ)

- 地図<: は、ウィンドウを作成するために、私はYPとしてXPとYベクトルとしてXベクトルを用いて、次のコードを使用しました点パターン(PPP)は

は本当に彼のサポート

0

に任意のヘルプやアドバイスをお願いします(再現性の問題を提供するために、コードがdirectlabelsパッケージからexample(direct.labels)から取られています)。

volcano3d <- melt(volcano) 
colnames(volcano3d) <- c("x", "y", "z") 

あなたがそうのように、Maunga Whauこの場合には、データの等高線図を描くことができます:

(p <- ggplot(volcano3d, aes(x, y, z = z)) + 
    geom_contour() 
) 

aesが式を取ることで、あなたのサンプルコードに問題があります、文字列ではありません。実際aes(x, y, z = freq)で行いますので、あなたが

aes(x = x, y = y, z = freq) 

または

aes_string(x = "x", y = "y", z = "freq") 

のいずれかを持っている必要があり通常、あなたは、明示的xyの引数名を与える必要はありません。

+0

'geom_path'が有用になるのは私には分かりません。これがあなたが望むものなら私に教えてください。 –

+0

エースを訂正していただきありがとうございます。それでも、私はコマンドを実行することはできません。それは私に教えてくれます:if(empty(new))return(data.frame()): TRUE/FALSEが必要な値がありません。 –

+0

パスは私のために重要です。なぜなら、島の形を私に与えてくれるので、私は価値がないから輪郭関数を制限する必要があるからです。火山の例は本当に素晴らしく、多くのもので遊ぶことができますが、地形図には限界がありません。我々は同じ火山を見たいと思うが、形が規則的ではない島では、プロットの多くのセクションが空である。私はそれが理にかなっていることを願っています。画像は説明しやすいが、ここに投稿することはできません。 –

関連する問題