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

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条回答 默认 最新

  • 爱晚乏客游 2021-03-14 15:17
    关注

    应该是你报错的这个.nc文件不完整,格式有问题才会这样的

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

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

悬赏问题

  • ¥20 一个python博客项目的相关图例
  • ¥15 轮廓提取也提取不到,有没有别的方法,如何解决?
  • ¥50 Js和c++如何将含有图片的excel文件上传到后台服务器
  • ¥15 光电神经网络,FPGA
  • ¥20 通过防火墙出入站阻止游戏程序联网失效
  • ¥15 鼠标是可以在QT界面上移动的,但是热拔插鼠标无法移动了同时板子上是没问题的,如何解决?
  • ¥15 iframe嵌套显示问题
  • ¥20 【UE4】别人打包好的ue4游戏我该如何在自己的ue4引擎上运行
  • ¥15 power bi发布的链接地址打不开
  • ¥15 pip list列表中有库,但是编译时就显示缺少库