绘图脚本#

1
#!/bin/sh
2
## Copyright (c) UPC lgl. All rights reserved.
3
# *************************************************#
4
# ** Filename:  imagesc.sh                       **#
5
# ** Author:Guolong Li(lgl)                      **#
6
# ** Date:2023.02.28                             **#
7
# ** E-mail:1712343961@qq.com                    **#
8
# ** Aim: use psimage to plot image              **#
9
# ** Environment: shell (Linux), Seismic Unix    **#
10
# *************************************************#
11

12
psfile=image.dat    # binaryfile
13
outname=image.eps
14
nx=481                # number of samples in 2nd (slow) dimension nx->n2
15
nz=234                # number of samples in 1st (fast) dimension nz->n1
16
dx=10                 # sampling interval X
17
dx=$dz
18
#dz=10                 # sampling interval Z
19

20

21
title=''     # title of plot
22
#xlabel='Distance(m)'  # X轴标签 'Distance(m)' 'Distance(km)'
23
#zlabel='Depth(m)'     # Z轴标签 'Time(s)' 'Time(ms)' 'Depth(km)' 'Depth(m)'
24
#legend=1              #  1 显示colorbar 0 不显示
25
#cchar='wrgb=1.0,0,0 grgb=1.0,1.0,1.0 brgb=0,0,1.0'
26
cchar='wrgb=0.6,0,0 grgb=0,0,0.8 brgb=0.376,0.376,0.376'
27
#units='m/s'        # unit label for legend
28

29
width=6.0           # width in inches of axes box  6.0
30
height=6.0           # height in inches of axes box  8.0  图像大小
31
labelfont=Times-Roman    # font name for axes label Times-Roman Helvetica
32
labelsize=18            # font size for axes labels  18
33
titlefont=Helvetica-Bold #font name for title
34
titlesize=24      #font size for title      24
35
xgrid2=none     # grid lines on axis 2 - none, dot, dash, or solid
36
zgrid1=none     # grid lines on axis 1 - none, dot, dash, or solid    图 网格
37
lstyle=vertright  # vertleft vertright   horibottom             legend 位置
38
legendfont=times_roman10  # times_roman10    font name for title
39
bckgr=0.5                 # background gray value   0.5
40
perc=99                  # percentile used to determine clip 100
41
lwidth=1.2     # colorscale (legend) width in inches 1.2
42
lheight=$height/3        # colorscale (legend) height in inches height/3   legend 大小
43
legendfont=times_roman10    # font name for legend title
44
lgrid=none   #  grid lines on legend axis - none, dot, dash, or solid  legend 网格
45
lx=     # colorscale (legend) x-position in inches           legend 位置
46
nxf=0                 # first sample in 2nd dimension
47
nzf=0                 # first sample in 1st dimension  坐标轴起始绘制
48
nxs=1.0        # factor by which to scale d2 before imaging  X
49
nzs=1.0        # factor by which to scale d1 before imaging  Z  缩放因子
50

51
## laplain 滤波
52

53
#suaddhead < $psfile n1=151>xx.su
54
#sukfrac <xx.su power=2 phasefac=-1 d1=$dz d2=$dx>xx.bin
55
#sustrip <xx.bin > lap.dat
56
#psfile=lap.dat
57
psimage  <  $psfile n1=$nz n2=$nx d2=$dx d1=$dz  width=$width height=$height \
58
    label1=$zlabel label2=$xlabel labelfont=$labelfont labelsize=$labelsize \
59
    title=$title titlefont=$titlefont titlesize=$titlesize \legend=$legend\
60
    lstyle=$lstyle legendfont=$legendfont $cchar  bckgr=$bckgr perc=$perc \
61
    lwidth=$lwith lheight=$lheight units=$units legendfont=$legendfont \
62
    grid2=$xgrid2 grid1=$zgrid1 lgrid=$lgrid lx=$lx f2=$nxf  f1=$nzf d1s=$nzs d2s=$nxs  > $outname
63

64
#Some example colormap settings:
65
# red white blue: wrgb=1.0,0,0 grgb=1.0,1.0,1.0 brgb=0,0,1.0
66
# white red blue: wrgb=1.0,1.0,1.0 grgb=1.0,0.0,0.0 brgb=0,0,1.0
67
# blue red white: wrgb=0.0,0.0,1.0 grgb=1.0,0.0,0.0 brgb=1.0,1.0,1.0
68
# red green blue: wrgb=1.0,0,0 grgb=0,1.0,0 brgb=0,0,1.0
69
# orange light-blue green: wrgb=1.0,.5,0 grgb=0,.7,1.0 brgb=0,1.0,0
70
# red light-blue dark blue: wrgb=0.0,0,1.0 grgb=0,1.0,1.0 brgb=0,0,1.0
71

72

73

74
#######
75
# brgb=0.0,0.0,0.0   red, green, blue values corresponding to black
76
# grgb=1.0,1.0,1.0   red, green, blue values corresponding to grey
77
# wrgb=1.0,1.0,1.0   red, green, blue values corresponding to white
78
# bhls=0.0,0.0,0.0   hue, lightness, saturation corresponding to black
79
# ghls=0.0,1.0,0.0   hue, lightness, saturation corresponding to grey
80
# whls=0.0,1.0,0.0   hue, lightness, saturation corresponding to white
81

82

83
#Legal font names are:
84
# AvantGarde-Book AvantGarde-BookOblique AvantGarde-Demi AvantGarde-DemiOblique
85
# Bookman-Demi Bookman-DemiItalic Bookman-Light Bookman-LightItalic
86
# Courier Courier-Bold Courier-BoldOblique Courier-Oblique
87
# Helvetica Helvetica-Bold Helvetica-BoldOblique Helvetica-Oblique
88
# Helvetica-Narrow Helvetica-Narrow-Bold Helvetica-Narrow-BoldOblique
89
# Helvetica-Narrow-Oblique NewCentrySchlbk-Bold NewCenturySchlbk-BoldItalic
90
# NewCenturySchlbk-Roman Palatino-Bold
91
# Palatino-BoldItalic Palatino-Italics Palatino-Roman
92
# SanSerif-Bold SanSerif-BoldItalic SanSerif-Roman
93
# Symbol Times-Bold Times-BoldItalic
94
# Times-Roman Times-Italic ZapfChancery-MediumItalic