用python编写一个将CGCS2000地理坐标系或者投影坐标系的栅格文件转化为WGS84坐标系

在Python中，可以使用 `GDAL`（通过 `osgeo.gdal` 模块）或 `rasterio` 等库来处理栅格数据的坐标系转换。以下是一个使用 `rasterio` 和 `pyproj` 的示例代码，将 CGCS2000（中国2000坐标系）转换为 WGS84（经纬度坐标系）： ```python import rasterio from rasterio.warp import calculate_default_transform, reproject, Resampling # 输入输出文件路径 input_file = "input_cgcs2000.tif" output_file = "output_wgs84.tif" # 定义目标坐标系为 WGS84 dst_crs = 'EPSG:4326' # WGS84 坐标系 with rasterio.open(input_file) as src: # 获取源坐标系和目标坐标系的转换参数 transform, width, height = calculate_default_transform( src.crs, dst_crs, src.width, src.height, *src.bounds) # 创建输出数据集的元数据 kwargs = src.meta.copy() kwargs.update({ 'crs': dst_crs, 'transform': transform, 'width': width, 'height': height }) with rasterio.open(output_file, 'w', **kwargs) as dst: # 对每个波段进行重投影 for i in range(1, src.count + 1): reproject( source=rasterio.band(src, i), destination=rasterio.band(dst, i), src_transform=src.transform, src_crs=src.crs, dst_transform=transform, dst_crs=dst_crs, resampling=Resampling.bilinear) ``` ### 说明： 1. **CGCS2000** 是中国的国家坐标系，通常以 **EPSG:4490** 表示地理坐标系，如果是投影坐标系则可能为 **EPSG:4542（3度带）** 或 **EPSG:4528（6度带）**。 2. **WGS84** 是全球通用的地理坐标系，对应 **EPSG:4326**。 3. 使用 `calculate_default_transform` 自动计算目标栅格的分辨率和变换矩阵。 4. 使用 `reproject` 函数进行重投影。 ### 依赖库安装： ```bash pip install rasterio pyproj ``` ---

创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考