将pandas数据帧转换为geoTIFF (python)

import pandas as pd
from osgeo import gdal, osr

# 创建一个示例的pandas数据帧
data = {'Latitude': [40.7128, 34.0522, 37.7749],
        'Longitude': [-74.0060, -118.2437, -122.4194],
        'Value': [10, 20, 30]}
df = pd.DataFrame(data)

# 创建一个空的栅格数据集
driver = gdal.GetDriverByName('GTiff')
dataset = driver.Create('output.tif', 3, 3, 1, gdal.GDT_Float32)

# 设置栅格数据集的地理参考信息
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)  # 使用WGS84坐标系
dataset.SetProjection(srs.ExportToWkt())
geotransform = (-180, 360/3, 0, 90, 0, -180/3)
dataset.SetGeoTransform(geotransform)

# 将pandas数据帧中的值写入栅格数据集
band = dataset.GetRasterBand(1)
band.WriteArray(df['Value'].values.reshape(3, 3))

# 关闭栅格数据集
dataset = None