from osgeo import gdal
import pandas
import numpy as np
csv_data = pandas.read_csv(r'D:\MOD\breakpoint.csv')#CSV所在位置
np_data = np.array(csv_data)
print(np_data)
np_data[np_data == -9999] = np.nan
b4dataset = gdal.Open(r'D:\MOD\breakpoint\breakpointtimes\vegetation2001zong.tif')#投影信息——将现有投影信息赋给将要生成的tif
im_geotrans = b4dataset.GetGeoTransform() # 获取仿射矩阵信息
im_proj = b4dataset.GetProjection() # 获取投影信息
datatype = gdal.GDT_Float64
driver = gdal.GetDriverByName("GTiff")
dataset = driver.Create(r'D:\MOD\biggestyear.tif', np_data.shape[1], np_data.shape[0], 1, datatype)
#tif影像存放位置
dataset.SetGeoTransform(im_geotrans) # 写入仿射变换参数
dataset.SetProjection(im_proj) # 写入投影
dataset.GetRasterBand(1).WriteArray(np_data)
CSV格式转换为TIFF格式(GEOTIFF)代码
最新推荐文章于 2024-08-13 14:46:55 发布