基于arcpy批量处理nc数据

欢迎客官移步到简研–简简单单科研中学习和交流

该文章主要完成：
1、定制下载ERA5数据
2、nc数据提取为各个单波段影像
3、提取后的各个波段进行栅格运算和保存

美剧词典网页版

美剧词典小程序版

1、定制下载ERA5数据
ERA在线编译器下载ERA5数据：利用数据平台的在线编译器toolbox-editor，编写脚本文件，将GRIB格式数据转为nc格式进行下载，脚本如下：

API下载ERA5数据

import cdsapi

c = cdsapi.Client()
for i in range (2019,2020):
    print (i)
    c.retrieve(
        'reanalysis-era5-single-levels-monthly-means',
        {
            'product_type': 'monthly_averaged_reanalysis',
            'variable': 'total_precipitation',
            'year': str(i),
            'month': [
                '01', '02', '03',
                '04', '05', '06',
                '07', '08', '09',
                '10', '11', '12',
            ],
            'time': '00:00',
            'format': 'netcdf',
        },
        '%sera5tp.nc' %(i))

2、nc数据处理
arcpy.MakeNetCDFRasterLayer_md工具介绍
arcpy.CopyRaster_management工具介绍
arcpy.MakeRasterLayer_management工具介绍

#coding=utf-8
import arcpy
import os
from arcpy.sa import *
#增加模块
arcpy.CheckOutExtension("spatial")

for i in range(1979,2020):
    print('start')
    #%s_ERA_ECV_air_temperature.nc 原始数据
    arcpy.MakeNetCDFRasterLayer_md(in_netCDF_file="【目录】/%s_ERA_ECV_air_temperature.nc" %str(i),
                                   variable="tas",
                                   x_dimension="lon",
                                   y_dimension="lat",
                                   out_raster_layer="tas_Layer%s" %str(i),
                                   band_dimension="time",
                                   dimension_values="",
                                   value_selection_method="BY_VALUE")

    #%stas.tif 已经将nc数据转为多波段的tif
    arcpy.CopyRaster_management(in_raster="tas_Layer%s" %str(i),
                                out_rasterdataset="【目录】/%stas.tif" %str(i),
                                config_keyword="",
                                background_value="",
                                nodata_value="1.#QNAN0e+000",
                                onebit_to_eightbit="NONE",
                                colormap_to_RGB="NONE",
                                pixel_type="",
                                scale_pixel_value="NONE",
                                RGB_to_Colormap="NONE"
                                )
    for j in range(1,13):
        print "MakeRas_tif%s_%s" %(str(i),str(j))
        arcpy.MakeRasterLayer_management("【目录】/%stas.tif" %str(i),"MakeRas_tif%s_%s" %(str(i),str(j)),"#","-180.125 -90.125 179.875 90.125","%s" %str(j))

    outPlus = Raster("MakeRas_tif%s_1" %str(i)) + Raster("MakeRas_tif%s_2" %str(i))

    outPlus.save("【目录】/%savetas.tif" %str(i))
    print('end')

这是分界线

安利个看美剧背单词，练习听说的工具|美剧词典

美剧词典是一款短情景形式的英语学习工具。首先呢，颜值是十分高的，界面简洁，配色清新；其次，有才华啊，目前共收录了近10万条字幕数据，10万条短情景和80万条词汇。短情景均为几秒的无字幕原声经典美剧视频片段，不占用多少时间就能让你get到原汁原味的单词含义和用处啦。美剧词典现包括“查询”、“单词本”、“音频"和"我的”四项功能。