应用背景
最近在做蒸散发的相关数据处理,上一篇用ArcPy进行处理,发现太慢了,于是继续寻找简便方法。此篇文章为可运行代码。
参考文献:《使用Python语言将.nc格式影像转换成TIF格式,.nc文件涉及到多个参数变量、12个月均的数据(多图层提取))_Jack_33学习之路的博客-程序员宅基地》
使用数据:E_2006_GLEAM_v3.6a.nc(数据放在百度网盘了)
链接:https://pan.baidu.com/s/1jb1f2lYyoP-gZlY07LPQvA?pwd=1111
提取码:1111
1. 原始代码
# -*- coding: utf-8 -*-
# 模块导入
import numpy as np
import netCDF4 as nc
from osgeo import gdal,osr,ogr
import os
import glob
import datetime
def NC_to_tiffs(data,Output_folder):
nc_data_obj = nc.Dataset(data)
# print(nc_data_obj,type(nc_data_obj)) # 了解NC_DS的数据类型,<class 'netCDF4._netCDF4.Dataset'>
# print(nc_data_obj.variables) # 了解变量的基本信息
# print(nc_data_obj)
Lon = nc_data_obj.variables['longitude'][:]
Lat = nc_data_obj.variables['latitude'][:]
u_arr = np.asarray(nc_data_obj.variables['v10']) # 这里根据需求输入想要转换的波段名称
#影像的左上角和右下角坐标
LonMin,LatMax,LonMax,LatMin = [Lon.min(),Lat.max(),Lon.max(),Lat.min()]
#分辨率计算
N_Lat = len(Lat)
N_Lon = len(Lon)
Lon_Res = (LonMax - LonMin) /(float(N_Lon)-1)
Lat_Res = (LatMax - LatMin) / (float(N_Lat)-1)
for i in range(len(u_arr[