nc文件转tif（python）

最新推荐文章于 2024-03-12 13:02:15 发布

_LLBB

最新推荐文章于 2024-03-12 13:02:15 发布

阅读量433

点赞数

文章标签： python numpy 机器学习

本文链接：https://blog.csdn.net/weixin_46058813/article/details/130001712

版权

# -*- coding: utf-8 -*-
"""
Created on Thu Jun  9 22:31:05 2022

@author: Binb
"""

import os
import netCDF4 as nc
import numpy as np
from osgeo import gdal, osr, ogr
import glob

data = r'E:/data/GLEAM/nc/SMsurf_1980-2021_GLEAM_v3.6a_YR.nc'
nc_data_obj = nc.Dataset(data)
print(nc_data_obj.variables)



def nc2tif(data, Output_folder):
    tmp_data = nc.Dataset(data)  # 利用.Dataset()方法读取nc数据
    Lat_data = tmp_data.variables['lat'][:]
    Lon_data = tmp_data.variables['lon'][:]
    re_data = tmp_data.variables['SMsurf'][:]
    tmp_arr = np.asarray(tmp_data.variables['SMsurf'])
    #tmp_arr = np.asarray(tmp_data.variables['time'])
   
    # 影像的左上角&右下角坐标
    Lonmin, Latmax, Lonmax, Latmin = [Lon_data.min(), Lat_data.max(), Lon_data.max(), Lat_data.min()]
    print(Lonmin, Latmax, Lonmax, Latmin)
   
    # 分辨率计算
    Num_lat = len(Lat_data)  # 5146