给出python代码，读取Geotiff格式的栅格数据，实现自定义四参数模型的坐标转换，生成为Geotiff栅格文件。

以下是示例代码，供参考：

import numpy as np
import gdal

# 定义自定义四参数模型的参数
dx = 10
dy = -5
sx = 0.01
sy = -0.005

# 读取原始的Geotiff栅格数据
src_file_path = 'input.tif'
src_raster = gdal.Open(src_file_path)
src_band = src_raster.GetRasterBand(1)
src_array = np.array(src_band.ReadAsArray())

# 获取原始栅格数据的地理信息
src_geotransform = src_raster.GetGeoTransform()
src_projection = src_raster.GetProjection()

# 计算新的地理信息
src_x_size = src_raster.RasterXSize
src_y_size = src_raster.RasterYSize
new_geotransform = (src_geotransform[0] + dx, src_geotransform[1] + sx, src_geotransform[2], src_geotransform[3] + dy, src_geotransform[4], src_geotransform[5] + sy)

# 创建新的Geotiff栅格文件
dst_file_path = 'output.tif'
driver = gdal.GetDriverByName('GTiff')
dst_raster = driver.Create(dst_file_path, src_x_size, src_y_size, 1, gdal.GDT_Float32)
dst_band = dst_raster.GetRasterBand(1)

# 将新的地理信息和原始的栅格数据写入新的Geotiff栅格文件
dst_raster.SetProjection(src_projection)
dst_raster.SetGeoTransform(new_geotransform)
dst_band.WriteArray(src_array)

# 释放资源
src_band = None
dst_band = None
src_raster = None
dst_raster = None

该代码使用gdal模块读取和写入Geotiff格式的栅格数据，实现了自定义四参数模型的坐标转换，并生成了新的Geotiff栅格文件。其中，dx和dy代表平移量，sx和sy代表缩放系数。经过计算后，新的地理信息通过SetGeoTransform方法设置到了新的Geotiff栅格文件中。