python調用GDAL實現幾何校正

引自https://blog.csdn.net/qq_27045589/article/details/81062586

def main():

    infile = "F:\\Temp_Data\\MOD13A2_20180610_H23V04_NDVI.tif"
    outfile = "F:\\Temp_Data\\MOD13A2_20180610_H23V04_NDVI_GEO.tif"

    dataset = gdal.Open(infile, gdal.GA_Update)
    gcps_list = [gdal.GCP(77.760364, 50.00, 0, 0, 0),
                 gdal.GCP(93.330364, 50.00, 0, 1199, 0),
                 gdal.GCP(78.320364, 40.01, 0, 1199, 1199),
                 gdal.GCP(65.270364, 40.01, 0, 0, 1199)]
    sr = osr.SpatialReference()
    sr.SetWellKnownGeogCS('WGS84')
    dataset.SetGCPs(gcps_list, sr.ExportToWkt())

    dst_ds = gdal.Warp(outfile, dataset, format='GTiff', tps=True,
                       xRes=0.01, yRes=0.01, dstNodata=-3000, srcNodata=-3000,
                       resampleAlg=gdal.GRIORA_Bilinear, outputType=gdal.GDT_Float32)
    # print dst_ds
    del dataset, gcps_list, sr


if __name__ == '__main__':
    main()

校正前：

校正后：

 

 

注：envi_register_doit函數校正后圖像范圍發生變化，而GDAL圖像校正不會出現這種問題。