riverByName('GTiff')
out_tif_name = Output_folder + '\\'+ data.split('\\')[-1].split('.')[0] + '_' + str(i+1) + '.tif'
out_tif = driver.Create(out_tif_name,N_Lon,N_Lat,1,gdal.GDT_Float32)
# 设置影像的显示范围
#-Lat_Res一定要是-的
geotransform = (LonMin,Lon_Res, 0, LatMax, 0, -Lat_Res)
out_tif.SetGeoTransform(geotransform)
#获取地理坐标系统信息,用于选取需要的地理坐标系统
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326) # 定义输出的坐标系为"WGS 84",AUTHORITY["EPSG","4326"]
out_tif.SetProjection(srs.ExportToWkt()) # 给新建图层赋予投影信息
#数据写出
out_tif.GetRasterBand(1).WriteArray(ndvi_arr_float[i]) # 将数据写入内存,此时没有写入硬盘
out_tif.FlushCache() # 将数据写入硬盘
out_tif = None # 注意必须关闭tif文件
def main():
Input_folder = 'F:\\data___python_test\\nc_to_tif\\nc'
Output_folder = 'F:\\data___python_test\\nc_to_tif\\tif_result'
# 读取所有nc数据
data_list = glob.glob(Input_folder + '\\*.nc4')
for i in range(len(data_list)):
data = data_list[i]
NC_to_tiffs(data,Output_folder)
print data + '-----转tif成功'
print'----转换结束----'
main()
本文作者:DQTDQT 限于作者水平有限,如文中存在任何错误,欢迎不吝指正、交流。
联系方式: QQ:1426097423 e-mail:duanquntaoyx@163.com
本文版权归作者和博客园共有,欢迎转载、交流,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文链接,如果觉得本文对您有益,欢迎点赞、探讨。
|