极地图像情节，如何用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)