python处理nc数据转换为tif格式

嘻嘻爱学习

已于 2022-08-11 10:49:31 修改

阅读量3.8k

点赞数 14

分类专栏：遥感数据学习文章标签： python 开发语言

于 2022-08-11 10:45:56 首次发布

本文链接：https://blog.csdn.net/kkkyyyxxx/article/details/126279935

版权

应用背景

最近在做蒸散发的相关数据处理，上一篇用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[