原始碼如下，其中pInput為輸入影像路徑，pOut為輸出影像路徑，gcs為輸出影像的地理座標系

int DoReproject(string pInput, string pOut,string gcs)
{
    GDALAllRegister();
    CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
    GDALDataset* pDTIn = (GDALDataset*)GDALOpen(pInput.c_str(), GA_ReadOnly);
    if (pDTIn == NULL)
    {
        printf("檔案開啟失敗: %s\n"
 
, pInput.c_str());
        return -1;
    }
    ///判斷是否與目標座標系統一直
    OGRSpatialReference obj_OSRS;
    obj_OSRS.SetWellKnownGeogCS(gcs.c_str());
    char *pszWKT = NULL;
    obj_OSRS.exportToWkt(&pszWKT);

    ///建立從UTM到84座標系的座標轉換關係
    void *  hTransformArg = GDALCreateGenImgProjTransformer((GDALDatasetH)pDTIn, pDTIn->GetProjectionRef(),
        NULL, pszWKT, false
 
, 0, 1);
    if (hTransformArg == NULL)
    {
        printf("無法建立轉換方程！\n");
        GDALClose((GDALDatasetH)pDTIn);
        return -1;
    }

    //計算輸出影像範圍
    int nWidth = 0;
    int nHeight = 0;
    double dOutTransform[6];
    ///hTransformArg：座標轉換關係
    CPLErr eErr = GDALSuggestedWarpOutput((GDALDatasetH)pDTIn,
        GDALGenImgProjTransform, hTransformArg,
        dOutTransform, &nWidth, &nHeight);
    /*if (eErr == CE_None)
    {
        GDALDestroyGenImgProjTransformer(hTransformArg);
        return -1;
    }*/
 


    //建立輸出影像
    GDALDataType eType = pDTIn->GetRasterBand(1)->GetRasterDataType();
    int nBandCnt = pDTIn->GetRasterCount();
    const char* format = "GTiff";
    GDALDriver* pDriver = GetGDALDriverManager()->GetDriverByName(format);
    GDALDataset* pDTReproj = (GDALDataset*)GDALCreate((GDALDriverH)pDriver,
        pOut.c_str(), nWidth, nHeight, nBandCnt, eType, NULL);
    if (pDTReproj == NULL)
    {
        printf("建立失敗: %s\n", pOut.c_str());
        return -1;
    }
    pDTReproj->SetProjection(pszWKT);
    pDTReproj->SetGeoTransform(dOutTransform);

    //psWarpOptions set
    GDALWarpOptions *psWarpOptions = GDALCreateWarpOptions();
    psWarpOptions->hDstDS = pDTReproj;///輸出影象指標
    psWarpOptions->hSrcDS = pDTIn;///原始影象指標
    psWarpOptions->nBandCount = 0;///處理的波段資料，0為處理所有波段
    psWarpOptions->eResampleAlg = GRA_Bilinear;///重取樣方式
    psWarpOptions->panSrcBands = NULL;
    psWarpOptions->panDstBands = NULL;
    psWarpOptions->pfnProgress = GDALTermProgress;
    psWarpOptions->pTransformerArg =
        GDALCreateGenImgProjTransformer((GDALDatasetH)pDTIn,
        GDALGetProjectionRef(pDTIn),
        (GDALDatasetH)pDTReproj,
        GDALGetProjectionRef(pDTReproj),
        FALSE, 0.0, 1);   //建立從原始到輸出結果的投影變換引數
    psWarpOptions->pfnTransformer = GDALGenImgProjTransform;  //座標變換函式指標


    //真正的重投影
    GDALWarpOperation oOperation;
    oOperation.Initialize(psWarpOptions);
    oOperation.ChunkAndWarpImage(0, 0, nWidth, nHeight);

    //銷燬轉換
    GDALDestroyGenImgProjTransformer(psWarpOptions->pTransformerArg);
    GDALDestroyWarpOptions(psWarpOptions);

    GDALClose(pDTIn);
    GDALClose(pDTReproj);

    return true;
}