極地圖像情節，如何用R做？

Question

我找不到一種直接的方法來在R中制作漂亮的圖像圖，但是在極坐標中。 我基本上試圖在MATLAB中找到'polarplot3d'函數的R等價物。 我一直在玩ggplot2包但沒有太多運氣。 我錯過了一個包含我正在嘗試的功能的包嗎？ 提前感謝任何指針。

好的，我正在努力讓自己更清楚自己要做什么。 假設我想定義一個極坐標網格，徑向增量為50m和θ為2.5度。 這應該看起來像一個飛鏢。

我的數據（r和下面代碼中的角度）對應於徑向距離測量和角度。 我期望的z值是在定義網格的上述增量內的r和角度之間的二元直方圖的計數。

我的數據如下：

# synthetic data for angle and distance #
angle <- rnorm(500,mean=90,sd=15)
r <- rnorm(500,mean=700,sd=200)

# bivariate histogram #
observations <- table(cut(angle,breaks=c(seq(0,360,by=2.5))),cut(r,breaks=c(seq(0,1400,by=50))))


# the 'z' data are in observations for each bin of bivariate histogram #
# hot to plot a polar coord image? #

Answer 1

在我的系統上渲染速度非常慢，但是

library(reshape2)
library(ggplot2)
mm <- melt(counts)
ggplot(mm,aes(Var1,Var2,fill=value))+geom_tile()+coord_polar()
ggsave("polar1.png")

似乎工作。

Answer 2

我認為以下可行。 使用mapproject()從maproj庫改造我的xy坐標acording到極投影（或其他），然后使用as.image()從fields包）函數從我的新coordiantes和我的Z值建立一個圖像對象。 最終使用image.plot() 。

library("mapproj")
xyProj <- mapproject(x, y, projection="conic", parameters=-90)
library("fields")
im <- as.image(z, x=xyProj)
image.plot(im)