Python中光栅图像的比较

Python中光栅图像的比较,python,gis,raster,Python,Gis,Raster,我有两张两年来拍摄的一个国家的光栅.tif图像()。我还有三个有行政边界的文件。我想创建一个代码,当从()获得管理边界文件时,它将创建一个CSV文件,其中包含比较两个图像的信息。这就是我到目前为止所做的: def night_lights(x,y): table( admin_level:"Admin Level" admin_id:"Admin ID" admin_name:"Admin Name"

我有两张两年来拍摄的一个国家的光栅.tif图像()。我还有三个有行政边界的文件。我想创建一个代码,当从()获得管理边界文件时,它将创建一个CSV文件,其中包含比较两个图像的信息。这就是我到目前为止所做的:

def night_lights(x,y):
    table(
    admin_level:"Admin Level"
    admin_id:"Admin ID"
    admin_name:"Admin Name"
    lum_1993:"Luminosity in 1993"
    lum_2013:"Luminosity in 2013"
    mean_change:"Mean Change between Years")

while ID_0 =! NULL
    admin_level.append(x.ID_)
    admin_id.append()
    admin_name.append()
    lum_1993.append(.mean)
    lum_2013.append(.mean)
    mean_change.append((lum_2013+lum_1993)/2)
    return table.to_csv('china_'admin_level)
我需要从文件中提取光栅信息和变量的帮助。我如何将所有这些放在一个表中并导出一个csv,而无需编写三次代码