FY-4A的nc数据批量转tif并投影

前段时间处理了一些fy4a的沙尘检测数据，综合了一些网上的帖子搞了一个批量处理的方法，主要实现方法是用python读取fy4a数据，然后用arcgis的模型构建器实现批量投影。

FY4A数据读取代码来自(7条消息) FY4A数据读取与插值_爱吃肉的悟空的博客-CSDN博客

模型构建器的设计参考(7条消息) arcgis_nc批量转tif（建模）(亲测适用)_红叶的舞台的博客-CSDN博客_nc转tif

 

1、稍微修改代码，加入循环

rootdir='D:\\fydata\\0603input'
savedir='D:\\fydata\\0603output\\'
list=os.listdir(rootdir)
for i in range(0,len(list)):
    file_path=os.path.join(rootdir,list[i])
    if os.path.isfile(file_path):
        # 存储nc的路径
        sv_path=savedir+os.path.basename(file_path)
        # 需要插值的经纬度
        new_lon, new_lat = utils.gen_lat_lon(70, 140, 60, 0, 0.04)
        # 得到数据块，并存储为nc
        fy4a.get_data_from_data_name(file_path, 'DST', new_lon, new_lat, sv_path=sv_path)
        # 数据对象
        logging.debug(fy4a.read_FY4A(file_path))

 

2、nc数据批量转tif

 

 

3、批量投影

 

 

4、结果