Gnuplot pm3d地图:不同的x2tics,y2tics

时间:2018-03-31 15:01:09

标签: gnuplot

我可以用:

绘制一个“ploe figure”
degtorad(x)=x/180*pi
radtodeg(x)=x/pi*180

set pm3d map
set pm3d interpolate 0,0
set logscale zcb

set palette model RGB rgbformulae 33,13,10

set size ratio 1

set xtics("{/Symbol b}=180{/Symbol \260}" 0)

set ytics("0" -90, "10" -80, "20" -70, "30" -60, "40" -50, "50" -40, "60" -30, "70" -20, "80" -10, "90" 0, "80" 10, "70" 20, "60" 30, "50" 40, "40" 50, "30" 60, "20" 70, "10" 80, "0" 90)

splot[-90:90][-90:90] 'TSC3.txt' u (radtodeg(cos(degtorad($1+90))*degtorad(90-$2))):(radtodeg(degtorad(90-$2)*sin(degtorad($1+90)))):3 notitle 

第1列包含“beta”,第2列是“alpha”,但无关紧要。

问题是

  1. 如何删除正确的标记(unset y2tics不起作用)?
  2. 如何在上方设置“beta = 0°”(set x2tics("{/Symbol b}=0{/Symbol \260}" 0)不起作用)?
  3. Here is I am now

1 个答案:

答案 0 :(得分:2)

右边界和上边界上未标记的抽搐与x2和y2轴无关。它们是x轴和y轴的镜像抽动。要关闭正确的抽搐,请使用set ytics ... nomirror

如果我没记错的话,早期版本splot根本不适用于x2和y2轴。至少使用5.2,您可以将x2轴链接到与x轴相同的范围set link x2,然后您可以将xics放在x2轴上:

degtorad(x)=x/180*pi
radtodeg(x)=x/pi*180

set pm3d map
set pm3d interpolate 0,0
set logscale zcb

set palette model RGB rgbformulae 33,13,10

set size ratio 1

set link x2
set x2tics ("{/Symbol b}=180{/Symbol \260}" 0)

set ytics nomirror ("0" -90, "10" -80, "20" -70, "30" -60, "40" -50, "50" -40, "60" -30, "70" -20, "80" -10, "90" 0, "80" 10, "70" 20, "60" 30, "50" 40, "40" 50, "30" 60, "20" 70, "10" 80, "0" 90)

splot[-90:90][-90:90] 'TSC3.txt' u (radtodeg(cos(degtorad($1+90))*degtorad(90-$2))):(radtodeg(degtorad(90-$2)*sin(degtorad($1+90)))):3 notitle