文章

🎨 一些简单的绘图脚本介绍

🎨 一些简单的绘图脚本介绍

在科研和日常学习中,绘图是一项重要技能。为了更高效地完成可视化任务,我整理了自己常用的GMT绘图脚本,便于日后复用和更新。

GMT中文手册,需要经产翻阅,内容非常详细 GMT中文手册

首先是我通常使用的绘图配置参数文件gmt.conf,这个文件负责定义图片的默认参数。大家根据喜好自行选择是否使用。

gmt.conf

1. ✏️ 绘制时空基线图

  • 我们在手动增删干涉对之后,还需要再绘制出对应时空基线图,可使用脚本:

    plot_baseline.csh

    1
    
      plot_baseline.csh intf.in baseline_table.dat
    

    图片说明文字


2. 🎨 绘制笛卡尔坐标系栅格数据(雷达坐标系)

  • 快速预览绘制笛卡尔坐标系的栅格数据,可使用脚本:

    plot_grd_JX.csh

    1
    
      plot_grd_JX.csh grd_file limitL limitU cpt
    

    图片说明文字


3. 🎨 绘制地理坐标系栅格数据(WGS84)

  • 快速预览绘制地理坐标系的栅格数据,可使用脚本:

    plot_grd.csh

    1
    
      plot_grd.csh grd_file limitL limitU
    

    图片说明文字


4. 📈 绘制时间序列图

  • 快速预览绘制时间序列,可使用脚本:

    plot_timeseries.csh

    1
    
      plot_timeseries.csh disp_file
    

    图片说明文字

5. 🔄 数据格式的转化

  • TIF栅格数据转GRD:

    1
    
      gdal_translate -of GMT A.tif A.grd
    
  • GRD栅格数据转TIF:

    1
    
      gdal_translate -of Gtiff A.grd A.tifff
    
  • shp要素数据转GMT:

    📢注意,在GMT6.1版本后,可直接绘制shp数据, 一个shapefile文件通常至少包含三个文件.shp.shx.dbf 。如果觉得文件太多,也可转为gmt格式,仅需一个文件绘图。

本文由作者按照 CC BY 4.0 进行授权