GDAL中无法使用transform函数，提示“ERROR1:latitudeorlongitudeexceededlimits错误。

这个错误通常是由于坐标系信息不正确而导致的。在使用transform函数时，确保你的坐标系信息设置正确，尤其是经度和纬度的范围限制。以下是示例代码，可用于正确设置坐标系信息：

import ogr
import osr

# 定义转换前后的坐标系
source = osr.SpatialReference()
source.ImportFromEPSG(4326) # WGS84坐标系

target = osr.SpatialReference()
target.ImportFromEPSG(3857) # Web Mercator投影坐标系

# 定义转换方法
transform = osr.CoordinateTransformation(source, target)

# 获取文件路径，此处为示例
filepath = 'example.tif'

# 打开文件并获取数据集
dataset = gdal.Open(filepath, gdal.GA_ReadOnly)

# 获取地理转换
geo_transform = dataset.GetGeoTransform()

# 获取原始数据的坐标系
data_crs = dataset.GetProjection()

# 创建输出数据集，根据输入数据集和转换后的坐标系创建
driver = gdal.GetDriverByName('GTiff')
output_dataset = driver.Create('output.tif', dataset.RasterXSize, dataset.RasterYSize, dataset.RasterCount, gdal.GDT_Float32)
output_dataset.SetProjection(target.ExportToWkt())
output_dataset.SetGeoTransform(geo_transform)

# 循环读取原始数据集，将每个像素转换为新坐标系下的像素值，并写入输出数据集
for i in range(dataset.RasterCount):
    band = dataset.GetRasterBand(i+1)
    band_data = band.ReadAsArray().astype(np.float32)
    new_band_data = np.zeros((dataset.RasterYSize, dataset.RasterXSize), dtype=np.float32)
    for y in range(dataset.RasterYSize):
        for x in range(dataset.RasterXSize):
            px, py, _ = gdal.ApplyGeoTransform(geo_transform, x, y