前段时间处理了一些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、结果
原创文章,作者:306829225,如若转载,请注明出处:https://blog.ytso.com/267185.html