gnuplot 绘制热图(heatmap) 先处理数据符合格式要求

gnuplot 画热图之前数据一定需要满足以下条件：按照某一列（第一列）排序，而且第一列中相邻的两行不相同则加入一段空行。

如果你的数据格式是错误的，就会报错！！！

Warning: Single isoline (scan) is not enough for a pm3d plot. Hint: Missing blank lines in the data file? See 'help pm3d' and FAQ.

是因为你缺少空行！！！！！

类似问题见：

gnuplot - 空图 gnuplot - IT工具网

heatmap - Empty plot gnuplot - Stack Overflow

解决办法：按照第一列（某一列排序），然后在第一列相邻行不等时插入空格，代码如下

############project
sort -k2 -n 1.txt > tem  
#将1.txt的数据按照第列排序并输出到新的文件tem
awk '{printf("%.5f %.5f %.5f %.5f\n",$1,$2,$3,$4)}' tem > temp.txt 
# tem包含四列内容，保留五位有效数字到tem.text
awk 'NR==1{print $0} NR>1{if($2!=prev){print ""} print $0; prev=$2}' temp.txt > program  #temp.txt的按照第二列排序，并在第二列相邻行但值不同的行之间加入空行，导入program中

gnuplot -p fig.gnu #绘图脚本
rm tem*   #删除临时文件

program 就是可用的数据，绘图如下

fig.gnu 内容如下：

#sort -k1 1.txt | awk 'NR>1 && $1!=prev {print ""} {prev=$1; print}' > 1 # 按照第一列从小到大排序并在第一列相同的地方插入空格
#awk '{printf "%.5f %.5f %.5f %.5f\n", $1, $2, $3 $4}' 1 > 4
set ylabel 'Frequency (Hz)'
set xlabel '{/symbol q} (deg)'
set xrange[0:32]
set yrange[0:33]
set xlabel font "Times New Roman,25" offset 0,-1
set ylabel font "Times New Roman,25" offset -2.5,0
set xtics font "Times New Roman,18"
set ytics font "Times New Roman,18"
set xtics nomirror  ###右上无刻度线
set ytics nomirror
set cbtics nomirror

set pm3d
set dgrid3d 200,200
set palette  rgb 22,13,-31 
#set cbtics ("0.5" -0.5,"-0.25" -0.25, "0" 0,"0.25" 0.25, "0.5" 0.5)
set cbrange [0:1]
#set cbtics 0.1
set cbtics font "Times New Roman, 18"
#
set view map
splot 'program' using 2:3:4 w pm3d t ''

颜色带palette可参考科学网—谈谈gnuplot（二十九）：色板（palette）设置 - 马欢的博文