卷心没有菜 2021-03-14 13:27 采纳率: 100%
浏览 262
已结题

python nc文件转tif报错怎么破?

问题是这样的,最近在处理CLDAS-V2.0的数据

这样的nc文件,我打算按照nc文件中的要素不同分文件夹转成tif文件并保存。

搞成这样的!!!!

但是我发现,不管什么怎么样,都智能处理54个nc文件

报错如下:

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

file_list = os.listdir(r'G:\date\CLDAS')

for i in file_list:
    if i.split('-')[1] == 'GST':
        NC_DS = nc.Dataset('G:/date/CLDAS/{}'.format(i))
        name = i.split('.')[0]+'.tif'
        OutTif = r'G:\date\CLDAS处理好的TIF影像\GST\{}'.format(name)
        Lat = NC_DS.variables['LAT'][:]
        Lon = NC_DS.variables['LON'][:]
        PAIR_MEAN = NC_DS.variables['TG_MEAN'][:]
        PAIR_MEAN = np.asarray(PAIR_MEAN)
        print(PAIR_MEAN)
        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)
        print(N_Lat,N_Lon,Lon_Res,Lat_Res)

        spei_ds = gdal.GetDriverByName('Gtiff').Create(OutTif, N_Lon, N_Lat, 1, gdal.GDT_Float32)
        geotransform = (LonMin, Lon_Res, 0, LatMin, 0, Lat_Res)
        spei_ds.SetGeoTransform(geotransform)

        srs = osr.SpatialReference()  
        srs.ImportFromEPSG(4326)  
        spei_ds.SetProjection(srs.ExportToWkt())

        spei_ds.GetRasterBand(1).WriteArray(PAIR_MEAN)  
        spei_ds.FlushCache()  
        spei_ds = None  
        NC_DS.close()

这个就是我的代码!!!

编译环境是anaconda + pycharm

  • 写回答

1条回答

      报告相同问题?

      相关推荐 更多相似问题

      问题事件

      • 系统已结题 4月16日
      • 已采纳回答 4月8日

      悬赏问题

      • ¥15 “glmnet”运行出错
      • ¥30 如何用matlab表达以下公式?
      • ¥15 在arm架构芯片上基于32位linux操作系统做内存检查
      • ¥15 怎么样才能禁止VS自动调整Windows窗体布局
      • ¥15 mysql5.7.40安装到Initializing database报错,如何解决?
      • ¥30 如何降低hdfs中datanode的JVM内存用量
      • ¥15 Android URL如何转成视频/音频,可行吗?
      • ¥20 Hive SQL数据查询,子查询
      • ¥15 c++字符串分割问题
      • ¥15 vue+uniapp