我试图从插值点数据中编写R中的轮廓极坐标图.换句话说,我在极坐标中有数据,我想绘制一个幅度值并显示插值.我想批量生产类似于以下的图(在OriginPro中生成):

我在R中最接近这一点的尝试基本上是:
### Convert polar -> cart
# ToDo #
### Dummy data
x = rnorm(20)
y = rnorm(20)
z = rnorm(20)
### Interpolate
library(akima)
tmp = interp(x,y,z)
### Plot interpolation
library(fields)
image.plot(tmp)
### ToDo ###
#Turn off all axis
#Plot polar axis ontop
Run Code Online (Sandbox Code Playgroud)
产生的东西如下:

虽然这显然不是最终产品,但这是在R中创建等高极坐标图的最佳方法吗?
除了2008年的档案邮件列表,我找不到任何关于这个主题的内容.我想我并没有完全致力于将R用于图(尽管这是我拥有数据的地方),但我反对手动创建.所以,如果有另一种语言具有此功能,请建议它(我确实看到了Python示例).
关于使用ggplot2的建议 - 我似乎无法使用geom_tile例程在polar_coordinates中绘制插值数据.我在下面列出了代码,说明了我的位置.我可以用直角坐标和极坐标绘制原始图形,但我只能得到插值数据以笛卡尔坐标绘制.我可以使用geom_point绘制极坐标中的插值点,但我无法将该方法扩展到geom_tile.我唯一的猜测是这与数据顺序有关 - 即geom_tile期望排序/有序数据 - 我已经尝试过每次迭代我都可以想到将数据排序为升序/降序方位角和天顶而没有变化.
## Libs
library(akima)
library(ggplot2)
## Sample data in az/el(zenith)
tmp = seq(5,355,by=10)
geoms <- data.frame(az = tmp,
zen = runif(length(tmp)),
value …Run Code Online (Sandbox Code Playgroud) 如何将matplotlib图的原点翻转到左上角 - 而不是默认的左下角?我正在使用matplotlib.pylab.plot来生成绘图(尽管如果有另一个绘图例程更灵活,请告诉我).
我正在寻找相当于matlab的命令:axis ij;
此外,我花了几个小时冲浪matplotlib帮助和谷歌,但没有得到答案.关于我可以在哪里找到答案的一些信息也会有所帮助.