注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

notes

@-@

 
 
 

日志

 
 
 
 

GMT笔记(zz from crust.cn)  

2010-10-20 18:49:21|  分类: 地理 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

psbasemap GMT画底图的指令
  一定要有的参数
    -B  设定图边界尺标的间隔
      a 文字批注的间隔
      f 边框标记的间隔
      g 隔线的间隔
        如果要针对X轴或Y轴做不同的设定,就在参数间加斜线/
      NEWS  每个轴加上标记
      news  每个轴不加标记
      :.title:  在图上面加title这个标题
      :title: 在轴上加title这个标题(只有在-Jx的情况下有用)
      例如-BNEWsf2a3g2:xx:/f3a3g2:yy::.title:
     -Jx(-JX)  线性投影(Linear Projection)一般绘图用
      -Jx后面接的数字代表每单位的长度(预设长度单位为英吋,加c的话就变公分)
       -JX后面接的数字代表轴的全长
       如果要针对X轴或Y轴做不同的设定,就在数字间加斜线/
     -Jm(-JM)  麦卡托投影(Mercator Projection)地图绘图用
      -Jm后面接的数字代表每度的长度(预设长度单位为英吋,加c的话就变公分)
       -JM后面接的数字代表整张图的宽
       如果要针对X轴或Y轴做不同的设定,就在数字间加斜线/
    -R  地图的经纬度范围
      西/东/南/北  依序填入经纬度
    >out.ps  把画出来的图输出到out.ps这个档案
  不一定要有的参数
    -E  旋转地图的视野方向
    -G  把图着色,用r/g/b调
  -K  在这一行下面还有后续的指令
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -V  画面上会显示本行指令的执行结果
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式
    psbasemap -JM7 -R117/126/18/26 -BNEWsf2a3g2:xx:/f3a3g2:y::.test: -P -V -K >123.ps
pscoast GMT画陆块海洋和海岸线的程序
  一定要有的参数
    -J  投影法,和前面的一样
    -R  地图的经纬度范围,和前面的一样
  不一定要有的参数
    -E  旋转地图的视野方向
    -G  把干地着色,用r/g/b调色
    -S  把湿地着色,用r/g/b调色
    -D  图的分辨率(f,h,i,l,c)左往右为分辨率高至低
    -I  加河流上去(1-10)大河至小河(a,r,i,c)大小河都画
    -W  画海岸线
    -N  画边界(1.国界,2.州界,3.海洋,a.全部)
  -K  在这一行下面还有后续的指令
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -V  画面上会显示本行指令的执行结果
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式
    pscoast -JM7 -R117/126/18/26 -G50/50/50 -S150/150/150 -W1 -Di -V -O -K >>123.ps

psxy GMT在地图上画点线的程序
  一定要有的参数
    -J  投影法,和前面的一样
    -R  地图的经纬度范围,和前面的一样
  不一定要有的参数
    -A  画大圆线
    -B  和先前psbasemap的一样
    -C  对照颜色,后面要接文件名称
    -G  把符号着色,用r/g/b调色
  -K  在这一行下面还有后续的指令
    -N  不要跳过超出范围的点
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -M  线分段
    -Sa 在图上画星星,后面要接大小
    -Sc 在图上画圆圈,后面要接大小
    -Sd 在图上画钻石,后面要接大小
    -Se 在图上画椭圆形,在资料的第3,4,5栏要有方向长轴短轴
    -Sh 在图上画六角形,后面要接大小
    -Ss 在图上画正方型,后面要接大小
    -St 在图上画三角形,后面要接大小
    -Si 在图上画反三角,后面要接大小
    -Sv 在图上画向量,在数据的第3,4栏要有方向长度
    -Sx 在图上画叉叉,后面要接大小
    -V  画面上会显示本行指令的执行结果
    -W  笔宽
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式一  把数据直接打在档案里
    psxy -JM7 -R117/126/18/26 -W1 -O -K -V??! >>123.ps
    120 23
    124 27
    !
  表示方式二  把数据另外打在数据文件里
    psxy 123.dat -JM7 -R117/126/18/26 -W1 -M -V -O -K >>123.ps
  在123.dat的档案里是这样(这是画线 条用的档案)
    >
    120 25
    123 25
    >
    121 23
    123 27
  在123.dat的档案里是这样(这是点用的档案)
    120 25
    123 25
    121 23
    123 27
  表示方式三  把更多的数据打在数据文件里
    psxy 123.dat -R117/126/18/26 -JM7 -Sc -Ccpt >> 123.ps
    在123.dat的档案里要有四栏数据,依序为X,Y,颜色代码,圆的大小

grdimage GMT画二维等高图的指令

  一定要有的参数
    grd grd档案的档名
    -C  对照颜色,后面要接文件名称
    -J  投影方式,和前面一样
  不一定要有的参数
    -B  和先前psbasemap的一样
  -K  在这一行下面还有后续的指令
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -R  地图的经纬度范围,和前面的一样
    -V  画面上会显示本行指令的执行结果
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式
    grdimage file.grd -JOc122/21/22/69/4.8i -R117/126/18/26 -Cfile.cpt -P -K -V >> 123.ps

psscale GMT画比例尺的指令
   一定要有的参数
    -D  定义比例尺的位置和大小
        如果比例尺是水平的  -D距中央/距顶端/长度/宽度h
        如果比例尺是垂直的  -D距左端/距中央/长度/宽度
  不一定要有的参数
    -B  定义比例尺的间隔,和前面一样
    -C  对照颜色,后面要接文件名称
  -K  在这一行下面还有后续的指令
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -V  画面上会显示本行指令的执行结果
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式
    psscale -Cfile.cpt -D2.8/-1/4/0.22h -Bf1000a2000:”Elevation (m)”: -O -K >> 123.ps

grdcontour GMT画二维等高线图的指令
  一定要有的参数
    grd grd档案的档名
    -C  等高线的间距
    -J  投影方式,和前面一样
  不一定要有的参数
    -A  等高线数值标示的间隔
    -B  和先前psbasemap的一样
    -E  视角
    -G  同一条等高线数值标示的间隔
  -K  在这一行下面还有后续的指令
    -L  等高线最低或最高的限制
    -N  等高线的单位(如果有下-A的参数才有效)
    -O  把这一行指令绘制的图迭加到先前画的图上
    -P  把图旋转90度
    -R  地图的经纬度范围,和前面的一样
    -T  高地(-T+)或低地(-T-)加上标记
    -V  画面上会显示本行指令的执行结果
    -X  移动整个图的水平位置
    -Y  移动整个图的垂直位置
  表示方式
    grdcontour file.grd -JM7 -R117/126/18/26 -C0.025 -W1 -A1 -P -K -V >> 123.ps
回到目录:回到 最前面
grd数据的处理
    转换成xyz值
      grd2xyz file.grd -R117/126/18/26 > file.xyz
    块状平均(I是精确度每多少度一格)
      blockmean file.xyz -R117/126/18/26 -I0.1 > file2.xyz
    转换成grd檔
      surface file2.xyz -R117/126/18/26 -I0.1 -Gfile3.grd
    查询档案信息
      minmax file.xyz
      grdinfo file.grd
qaidam
Qdm/user/byk/SUMMER/2003
/home_gobi/byk/SECONDARY/topo

  评论这张
 
阅读(1033)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017