gnuplot 2D polar plot with heatmap from 3D dataset - possible?

前端 未结 2 1767
爱一瞬间的悲伤
爱一瞬间的悲伤 2020-12-16 04:51

I\'d like to take a plot that I have done as a surface in 3D using cartesian coordinates and view it as a heatmap in 2D in POLAR coordinates. The reason for this is I am obt

2条回答
  •  隐瞒了意图╮
    2020-12-16 05:57

    The image plotting works only for equally distributed rectangular grids, just like any bitmap image is arranged. Otherwise you must use splot with pm3d.

    The set grid polar works only for 2D, so you must use multiplot to overlay your heatmap with the polar grid. Here a, quite lengthy, example to show you how it might work:

    reset
    set terminal pngcairo size 800,800
    set output '3d-polar.png'
    
    set lmargin at screen 0.05
    set rmargin at screen 0.85
    set bmargin at screen 0.1
    set tmargin at screen 0.9
    
    set pm3d map
    unset key
    
    set multiplot
    
    # plot the heatmap
    set parametric
    set isosamples 500
    
    unset border
    unset xtics
    unset ytics
    
    set angles degree
    r = 6
    set urange[0:r] # radius
    set vrange[0:360] # angle
    set xrange[-r:r]
    set yrange[-r:r]
    set colorbox user origin 0.9,0.1 size 0.03,0.8
    splot u*cos(v), u*sin(v), (cos(v)*besj0(2*u))**2
    
    # now plot the polar grid only
    set style line 11 lc rgb 'white' lw 2
    set grid polar ls 11
    set polar
    set rrange[0:r]
    unset raxis
    set rtics format '' scale 0
    unset parametric
    set for [i=0:330:30] label at first (r+0.35)*cos(i), first (r+0.35)*sin(i)\
    center sprintf('%d', i)
    plot NaN w l
    unset multiplot
    

    The result is:

    enter image description here

    And now some details about some tricks:

    • In order to get a square size, you can't use set size ratio 1, because the margins differ for the 2D and 3D plots, even if you would specify some absolute margins. Therefore, I set a square canvas size (terminal option size 800,800), and set appropriate absolute margins.

    • You cannot unset rtics because then the grid would disappear.

    • The grid labels must be set manually.

    • The colorbox was also set manually because otherwise it would have overlapped with the 0 label.

    • Plotting NaN does only plot the grid

提交回复
热议问题