import commands
mean=commands.getstatusoutput('gmt grdinfo -L1 -L2 -C topo1_tra.nc')[1].split('\t')[13]
std=commands.getstatusoutput('gmt grdinfo -L1 -L2 -C topo1_tra.nc')[1].split('\t')[14]
print mean,std
!gmt grdinfo topo1_tra.nc
commands.getstatusoutput('gmt grdinfo topo1_tra.nc')[1].split('\n')[7].split(':')[-1].lstrip().rstrip()
注意:getstatusoutput容易发送阻塞,当取不到值时,一直处于等待状态。