GMT pstext

GMT pstext

     pstext常常用来注释图。pstext  的命令形式如下

     pstext textfile  −Jparameters−Rwest/east/south/north[r][ −A ] [−B[p|s]parameters ] [−Cdx/dy ] [−D[j]dx[/dy][v[pen]] [−Eazim/elev[+wlon/lat[/z]][+vx0/y0]] [ −Gcolor ] [−H[i][nrec]        ]     [−Jz|Zparameters ] [−K ] [ −L ] [ −N ] [−O ] [ −P ] [−Spen ] [−U[just/dx/dy/][c|label]] [ −V ] [−W[color,][o|O|c|C[pen]]] [−X[a|c|r][x-shift[u]]] [−Y[a|c|r][y-shift[u]]] [ −Z[zlevel|+] ] [−:[i|o] ] [−ccopies ] [−f[i|o]colinfo ] [−m[flag] ]
     textfile的每一行为一个字符串信息,基本格式如下

     x y size angle fontno justify text

     x,y 字体位置

     size 大小

     angle 与水平方向夹角

     fontno 字符类型

     justify 对齐方式

     text  文本

     应用时需要注意的是x,y的与你设置 -R 参数的范围,我刚开始就是没有注意,怎么都看不到字。

这里主要说下特殊字符的实现需要例子的可以去我资源里面下载:

      @~ :转换为希腊字符

      @%no%    设置字符类型为no

      @%%:        设置字符类型为起始字符类型

      @-               下标开关

      @+              上标开关

      @#               字符大写

      @;color;      改变填充颜色

      @;;               颜色设置复原

      @:size:       改变字体大小

      @::              字体大小设置复原

      @_ text@_ 字符串加下划线               

      @@              输出@字符

      @!<char1><char2>  char2与char1重叠                

      @e, @o, @a, @E, @O, @A   与字符类型配合可以有很多特殊字符       

GMT(Generic Mapping Tools)是一款开源的地图制图软件,提供了丰富的地图制图功能和API接口。Fortran可以通过调用GMT的API接口来实现地图制图。 以下是一个简单的Fortran程序,通过调用GMT的API接口绘制了一张简单的地图: ``` program gmt_example implicit none integer :: i, j integer :: iret real :: lon(3), lat(3), z(3) character(len=80) :: command ! 创建GMT工作区 command = "gmt begin gmt_example" call system_trim(command, iret) ! 设置地图区域和投影方式 command = "gmt basemap -R-180/180/-90/90 -Jm1" call system_trim(command, iret) ! 绘制三角形 lon = (/ -100.0, -80.0, -90.0 /) lat = (/ 30.0, 20.0, 10.0 /) z = (/ 1.0, 2.0, 3.0 /) command = "gmt psxy -R -J -W1p -Gred -K -O <<EOF" do i = 1, 3 write(command, "(3f10.5)") lon(i), lat(i), z(i) end do command = trim(command) // "EOF" call system_trim(command, iret) ! 添加标题和注释 command = "gmt pstext -R -J -N -F+f16p+jTL -O <<EOF" command = trim(command) // "-100.0 30.0 Example" call system_trim(command, iret) ! 结束GMT工作区 command = "gmt end" call system_trim(command, iret) end program gmt_example ``` 在程序中,我们首先使用`gmt begin`命令创建GMT工作区,然后使用`gmt basemap`命令设置地图区域和投影方式。接着,我们使用`gmt psxy`命令绘制一个三角形,并使用`gmt pstext`命令添加标题和注释。最后,我们使用`gmt end`命令结束GMT工作区。 需要注意的是,在Fortran中调用GMT API需要使用`call system_trim(command, iret)`来执行GMT命令,并检查返回值`iret`是否为0,以确保命令执行成功。 编译程序时需要链接GMT库,可以使用以下命令: ``` $ gfortran gmt_example.f90 -o gmt_example -lgmt ``` 运行程序,将会生成一个名为`gmt_example.ps`的PostScript格式地图文件。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值