爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6667|回复: 4

[图形美化] (已解决)grads色标上数字出问题

[复制链接]
回帖奖励 2 金钱 回复本帖可获得 1 金钱奖励! 每人限 1 次
发表于 2016-10-19 20:16:41 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
本帖最后由 东风急流 于 2016-10-22 10:18 编辑

参照清风大神的一页多图统一色标轴,画了张图,但是色标上数字有点小问题,最后一个数字后面有个单引号。。。感觉不那么完美。
这个玩意是从哪里出来的呢
sst_lujin.png


我的脚本如下,请大家指教
'reinit'
'open d:\tu3\kz6.ctl' ;*利用此前的ARWpost输出文件读取地形数据
'open d:\tu3\cp6.ctl' ;*利用此前的ARWpost输出文件读取地形数据
*颜色设置(用清风调色板克隆)
'set rgb 20 160 0 200'
'set rgb 21 130 0 220'
'set rgb 22 30 60 255'
'set rgb 23 0 160 255'
'set rgb 24 0 200 200'
'set rgb 25 0 210 140'
'set rgb 26 0 220 0'
'set rgb 27 160 230 50'
'set rgb 28 230 220 50'
'set rgb 29 230 175 45'
'set rgb 30 240 130 40'
'set rgb 31 250 60 60'
'set rgb 32 240 0 130'
'run D:\tu3\lujin\13colors.gs';*为后续统一色标轴做准备

*画上图
'set vpage 0.5 5.5 1 8'
'set dfile 1'
'set grads off' ;*关闭grads角标 以便编辑个性的下标
'set grid off'
'set mproj latlon'
'set poli on' ;*画出国界省界线
'set mpdset hires' ;*设置高分辨率地图
'set xlint 5' ;*定义横坐标的标记间隔
'set ylint 5' ;*定义纵坐标的标记间隔
'set xlopts 1 4 0.15' ;*横坐标 颜色 线宽 大小
'set ylopts 1 4 0.15' ;*纵坐标 颜色 线宽 大小
'set gxout shaded'
'set lon 110 125'
'set lat 15 25'
'set t 7'
'd SST-273.15'

*****************************************************************
*控制试验路径
i=1
while(i<=13) ;*文件内容n行之内
read_file=read('d:/tu3/kz6.txt')
read_code=sublin(read_file,1) ;*读取文件打开代码
if(read_code>0);
break;
endif ;*读取正常为0,其他要么出错,要么结束
xbefore=x0
ybefore=y0
read_line=sublin(read_file,2) ;*读取文件行内容
lon0=subwrd(read_line,1) ;*取得第1部分 经度
lat0=subwrd(read_line,2) ;*取得第2部分 纬度
date0=subwrd(read_line,3) ;*取得第3部分 时间
'q w2xy ' lon0 ' ' lat0 ;*将经纬度坐标转换为GRADS坐标
x0=subwrd(result,3) ;*取得转换后的X坐标
y0=subwrd(result,6) ;*取得转换后的Y坐标
'draw mark 5 'x0' 'y0' 0.1' ;*作出各个坐标点 24
if(i>1) ;*第一个点作为起始点
'draw line 'xbefore' 'ybefore' 'x0' 'y0'' ;
endif;
i=i+1
Endwhile
***************************************************************

*画下图
'set vpage 5.5 10.5 1 8'
*'set parea 6 11 1 8'
'set dfile 2'
'set grads off' ;*关闭grads角标 以便编辑个性的下标
'set grid off'
'set mproj latlon'
'set poli on' ;*画出国界省界线
'set mpdset hires' ;*设置高分辨率地图
'set xlint 5' ;*定义横坐标的标记间隔
'set ylint 5' ;*定义纵坐标的标记间隔
'set xlopts 1 4 0.15' ;*横坐标 颜色 线宽 大小
'set ylopts 1 4 0.15' ;*纵坐标 颜色 线宽 大小
'set gxout shaded'
'set lon 110 125'
'set lat 15 25'
'set t 7'
'd SST.2-273.15' ;

*****************************************************************
*耦合试验路径
i=1
while(i<=13) ;*文件内容n行之内
read_file=read('d:/tu3/cp6.txt')
read_code=sublin(read_file,1) ;*读取文件打开代码
if(read_code>0);
break;
endif ;*读取正常为0,其他要么出错,要么结束
xbefore=x0
ybefore=y0
read_line=sublin(read_file,2) ;*读取文件行内容
lon0=subwrd(read_line,1) ;*取得第1部分 经度
lat0=subwrd(read_line,2) ;*取得第2部分 纬度
date0=subwrd(read_line,3) ;*取得第3部分 时间
'q w2xy ' lon0 ' ' lat0 ;*将经纬度坐标转换为GRADS坐标
x0=subwrd(result,3) ;*取得转换后的X坐标
y0=subwrd(result,6) ;*取得转换后的Y坐标
'draw mark 5 'x0' 'y0' 0.1' ;*作出各个坐标点 24
if(i>1) ;*第一个点作为起始点
'draw line 'xbefore' 'ybefore' 'x0' 'y0'' ;
endif;
i=i+1
Endwhile
***************************************************************
*'xcbar 9.5 9.7 1 1.2 -fw 0.1 -fh 0.15 -ft 3 -fs 1 -fo 0 -line off'
'run plotcbar D:\tu3\lujin\13colors.gs 1 0.9 0 5.5 1.8'
'printim d:/tu3/lujin/sst_lujin.png white' ;*图形输出
;


密码修改失败请联系微信:mofangbao
发表于 2016-10-20 11:21:21 | 显示全部楼层

回帖奖励 +1 金钱

'run D:\tu3\lujin\13colors.gs';*为后续统一色标轴做准备
32与'中间多加一个空格试试
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-10-20 15:37:49 | 显示全部楼层
zhaolanxiai 发表于 2016-10-20 11:21
'run D:\tu3\lujin\13colors.gs';*为后续统一色标轴做准备
32与'中间多加一个空格试试

谢谢,确实是这个问题,又看了清风的帖子,上面有这个的提醒,还是自己不够细心。
密码修改失败请联系微信:mofangbao
发表于 2016-10-23 15:12:18 | 显示全部楼层

回帖奖励 +1 金钱

厉害厉害,看的我都头晕了
密码修改失败请联系微信:mofangbao
发表于 2018-9-17 11:10:54 | 显示全部楼层

回帖奖励 +1 金钱

我刚好需要找这个问题,谢谢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表