Register
Login
Resources
Docs Blog Datasets Glossary Case Studies Tutorials & Webinars
Product
Data Engine LLMs Platform Enterprise
Pricing Explore
Connect to our Discord channel

reduce.sh 2.4 KB

You have to be logged in to leave a comment. Sign In
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
  1. #!/usr/bin/env bash
  2. # Test grdmath stacking mode -S for available operators
  3. ps=reduce.ps
  4. gmt set PS_MEDIA letter MAP_TITLE_OFFSET 4p FONT_TITLE 12p
  5. # Create 3 small grids with integers in 0-100 range
  6. cat << EOF | gmt xyz2grd -R0/2/0/2 -I1 -r -Z -G1.grd
  7. 43
  8. 31
  9. 82
  10. 82
  11. EOF
  12. cat << EOF | gmt xyz2grd -R0/2/0/2 -I1 -r -Z -G2.grd
  13. 96
  14. 74
  15. 78
  16. 79
  17. EOF
  18. cat << EOF | gmt xyz2grd -R0/2/0/2 -I1 -r -Z -G3.grd
  19. 77
  20. 71
  21. 2
  22. 48
  23. EOF
  24. # Plot the three grids on top
  25. gmt grd2xyz 1.grd | gmt pstext -R0/2/0/2 -JX1.75i -P -B0g1 -B+t1 -K -Y8.5i -F+f18p+jCM > $ps
  26. gmt grd2xyz 2.grd | gmt pstext -R -J -O -K -B0g1 -B+t2 -X2.25i -F+f18p+jCM >> $ps
  27. gmt grd2xyz 3.grd | gmt pstext -R -J -O -K -B0g1 -B+t3 -X2.25i -F+f18p+jCM >> $ps
  28. # Plot MEAN MEDIAN MODE on next row
  29. gmt grdmath [123].grd -S MEAN = tmp.grd
  30. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -B0g1 -B+tMEAN -K -Y-2i -X-4.5i -F+f18p+jCM+z%0.1f >> $ps
  31. gmt grdmath [123].grd -S MEDIAN = tmp.grd
  32. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tMEDIAN -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  33. gmt grdmath [123].grd -S MODE = tmp.grd
  34. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tMODE -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  35. # Plot STD, MAD, LMSSCL on next row
  36. gmt grdmath [123].grd -S STD = tmp.grd
  37. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -B0g1 -B+tSTD -K -Y-2i -X-4.5i -F+f18p+jCM+z%0.1f >> $ps
  38. gmt grdmath [123].grd -S MAD = tmp.grd
  39. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tMAD -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  40. gmt grdmath [123].grd -S LMSSCL = tmp.grd
  41. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tLMSSCL -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  42. # Plot MIN, MAX, RMS on next row
  43. gmt grdmath [123].grd -S MIN = tmp.grd
  44. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -B0g1 -B+tMIN -K -Y-2i -X-4.5i -F+f18p+jCM+z%0.1f >> $ps
  45. gmt grdmath [123].grd -S MAX = tmp.grd
  46. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tMAX -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  47. gmt grdmath [123].grd -S RMS = tmp.grd
  48. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tRMS -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  49. # Plot ADD, SUB, AND on next row
  50. gmt grdmath [123].grd -S ADD = tmp.grd
  51. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -B0g1 -B+tADD -K -Y-2i -X-4.5i -F+f18p+jCM+z%0.1f >> $ps
  52. gmt grdmath [123].grd -S SUB = tmp.grd
  53. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -K -B0g1 -B+tSUB -X2.25i -F+f18p+jCM+z%0.1f >> $ps
  54. gmt grdmath [123].grd -S AND = tmp.grd
  55. gmt grd2xyz tmp.grd | gmt pstext -R -J -O -B0g1 -B+tAND -X2.25i -F+f18p+jCM+z%0.1f >> $ps
Tip!

Press p or to see the previous file or, n or to see the next file

Comments

Loading...