Polar contour plot

[Matteo Ottaviani]

Hallo,

does anyone know a direct way to obtain filled contour plots in polar coordinates, like here:

http://www.idlcoyote.com/cg_tips/polarwind.php

WITHOUT exploiting the Coyote library?

The POLAR_CONTOUR function produces results on a cartesian grid....

Thank you!

[Matteo Ottaviani]

I forgot to mention that I do not need to plot over a map projection, just simply produce a filled contour plot.

[markus.sc...@gmail.com]

Have you tried 

m = MAP('PolarStereographic')
and otherwise following the code on 

https://www.nv5geospatialsoftware.com/docs/CONTOUR.html

You may need to manually add data to get an overlap at longitude -180 and 180, and manually add the poles at all longitudes.

Hopefully, the map sets up all the rest you need.

[Matteo Ottaviani]

Hi Markus! Thanks much for answering. As a matter of fact I did try what you suggested. The problem is that I can't find a way to get it with 0deg at the center (Polar stereographic has 90deg at the center). 'Center latitude' and "limits" do not seem to work correctly in my case.

Basically I'd like a direct translation of something like this, with the classical coordinate system (r, theta), with 0<r<90 and 0<theta<360:

https://stackoverflow.com/questions/70315757/customise-polar-contour-plot-axes-style