本文整理汇总了Python中nco.Nco.ncap2方法的典型用法代码示例。如果您正苦于以下问题：Python Nco.ncap2方法的具体用法？Python Nco.ncap2怎么用？Python Nco.ncap2使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类nco.Nco的用法示例。
在下文中一共展示了Nco.ncap2方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: str
# 需要导入模块: from nco import Nco [as 别名]
# 或者: from nco.Nco import ncap2 [as 别名]
dirs.append(dir_ini)
dir_hs = 'processed_hillshade'
dirs.append(dir_hs)
for dir_processed in dirs:
if not os.path.isdir(os.path.join(idir, dir_processed)):
os.mkdir(os.path.join(idir, dir_processed))
pvars = ('thk', 'usurf', 'velsurf_mag', 'velbase_mag')
fill_value = -2e9
v_str = ' '.join('='.join([x, str(fill_value) + ';']) for x in pvars)
m_str = 'sftgif=mask*0; where(thk>10) {sftgif=1;}; where(usurf>300) {sftgif=3;};'
ncap2_str = 'where(thk<10) {{ {} }}; {}'.format(v_str, m_str)
exp_files = glob(os.path.join(idir, 'state', '*.nc'))
for exp_file in exp_files:
exp_basename = os.path.split(exp_file)[-1].split('.nc')[0]
exp_nc_wd = os.path.join(idir, dir_nc, exp_basename + '.nc')
exp_gtiff_wd = os.path.join(idir, dir_gtiff, exp_basename + '.tif')
logger.info('masking variables where ice thickness < 10m')
nco.ncap2(input='-s "{}" {}'.format(ncap2_str, exp_file), output=exp_nc_wd, overwrite=True)
opt = [c.Atted(mode="o", att_name="_FillValue", var_name=myvar, value=fill_value) for myvar in pvars]
nco.ncatted(input=exp_nc_wd, options=opt)
logger.info('extracting ice-noice transition')
exp_ini_wd = os.path.join(idir, dir_ini, exp_basename + '.shp')
cmd = ['extract_interface.py', '-m', 'sftgif', '-t', 'ice_noice', '--epsg', '26710', '-o', exp_ini_wd, exp_nc_wd]
sub.call(cmd)
for mvar in pvars:
m_exp_nc_wd = 'NETCDF:{}:{}'.format(exp_nc_wd, mvar)
m_exp_gtiff_wd = os.path.join(idir, dir_gtiff, mvar + '_' + exp_basename + '.tif')
logger.info('Converting variable {} to GTiff and save as {}'.format(mvar, m_exp_gtiff_wd))
gdal.Translate(m_exp_gtiff_wd, m_exp_nc_wd, options=gdal_gtiff_options)