junlysky的个人博客分享 http://blog.sciencenet.cn/u/junlysky

博文

GMT-4.5.12的一个bug(grdpaste)

已有 2379 次阅读 2015-4-10 15:28 |系统分类:科研笔记|关键词:学者

     GTOPO30,即global digital elevation model (DEM),水平格网间隔为弧度30'(最大有1km)。数据块在全球的分布如下图:中国区域包含了四个DEM文件E060N90、E100N90、E060N40、E100N40。

    例如,对于地震多发的四川地区,经纬度包含在lat26-lat34,lon97-lon110,使用GTOPO30画四川地区的grd文件就需要将E60N40、E100N40两个文件制作的grd文件合并。GMT-4.5.12中可以使用grdpaste命令合并。grdpaste用法如下:

usage: grdpaste <file_a> <file_b> -G<outfile> [-V] [-f[i|o]<colinfo>]

    where file_a and file_b are to be combined into outfile.
    file_a and file_b must have same dx,dy and one edge in common.
    If in doubt, run grdinfo first and check your files.
    Use grdcut and/or grdsample to adjust files as necessary.


在运行paste时候,即使是两个相邻的有公共边的文件,有时会报错,例如:

#! /bin/tcsh -f

set rangeEN="100/140/40/60"
set rangeES="100/140/0/40"
set rangeWS="60/100/0/40"
set rangeWN="60/100/40/60"
grdraster 6 -I0.5m -R$rangeEN -GChinaEN.grd
grdraster 7 -I0.5m -R$rangeES -GChinaES.grd
grdraster 8 -I0.5m -R$rangeWS -GChinaWS.grd
grdraster 9 -I0.5m -R$rangeWN -GChinaWN.grd

grdpaste ChinaES.grd ChinaEN.grd -GChinaE.grd              
grdpaste ChinaWN.grd ChinaWS.grd -GChinaW.grd          
grdpaste ChinaW.grd ChinaE.grd -GChina.grd                    

运行该gmt脚本,会报错:              

grdpaste:  Grids do not share a common edge!


此时,就很纳闷了,明明是有公共边的却报错。

尝试在grdpaste中 a.grd与b.grd的顺序调换,则程序不会报错。

此处改为:grdpaste ChinaEN.grd ChinaES.grd -GChinaE.grd              #ok


结论:这是gmt-4.5.12的一个bug

如果语句 grdpaste <file_a> <file_b> -G<outfile> 报错 Grids do not share a common edge!

可互换<file_a>与<file_b>的位置。


























将文件E060N90、E100N90、E060N40、E100N40合并成一个大grd文件:


set rangeEN="100/140/40/60"
set rangeES="100/140/0/40"
set rangeWS="60/100/0/40"
set rangeWN="60/100/40/60"
grdraster 6 -I0.5m -R$rangeEN -GChinaEN.grd
grdraster 7 -I0.5m -R$rangeES -GChinaES.grd
grdraster 8 -I0.5m -R$rangeWS -GChinaWS.grd
grdraster 9 -I0.5m -R$rangeWN -GChinaWN.grd

grdpaste ChinaEN.grd ChinaES.grd -GChinaE.grd              #ok
grdpaste ChinaWN.grd ChinaWS.grd -GChinaW.grd           #ok
grdpaste ChinaW.grd ChinaE.grd -GChina.grd                  #ok



http://blog.sciencenet.cn/blog-1056595-881340.html

上一篇:兴趣,教师,压力教育
下一篇:mark learning python

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-5-28 21:09

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部