Arcgis之栅格数据转投影转换（84转2000）

Question

如题所述

Answer 1

案例：将栅格（.tif）格式数据的投影信息由84坐标系转换为2000坐标系

场景： 想要将data.tif格式数据的投影信息转换为2000坐标系并输出为data_srs.tif

脚本：Python脚本示范

def ProjectionTransformation_WGS84To2000(src_data, dst_data):

    arcpy.env.workspace = r'C:/Users/Administrator/Documents/ArcGIS/Default.gdb'

    # src_sr = arcpy.SpatialReference("WGS 1984")

    src_sr = arcpy.SpatialReference("WGS 1984 UTM Zone 48N")

    # dst_sr = r"prj/China Geodetic Coordinate System 2000.prj"

    dst_sr = arcpy.SpatialReference("CGCS2000 3 Degree GK CM 114E")

    print(src_sr.name, dst_sr.name)

    # 创建投影转换

    CustomGeoTransformation = "GEOGTRAN[METHOD['Geocentric_Translation'],PARAMETER['X_Axis_Translation',''],PARAMETER['Y_Axis_Translation',''],PARAMETER['Z_Axis_Translation','']]"

    arcpy.CreateCustomGeoTransformation_management("WGS84_TO_CGCS2000",

                                                   src_sr,

                                                   dst_sr,

                                                   CustomGeoTransformation)

    arcpy.ProjectRaster_management(src_data,

                                   dst_data,

                                   dst_sr,

                                   "BILINEAR",

                                   "#",

                                   "WGS84_TO_CGCS2000",

                                   "#",

                                   src_sr)

ProjectionTransformation_WGS84To2000(r'data.tif', r'data_srs.tif')