关于#arcpy#的问题，如何解决？

Question

将1991-1999.tif、2000-2008.tif，这两个栅格文件，需要通过RasterToNumPyArray (arcpy)、numpy将这两个栅格数据在时间维度进行拼接拼接，使其成为1991-2008.tif从而得到土地利用变化情况。
我的尝试：

arr_1991 = arcpy.RasterToNumPyArray("F:\study\luoyang text\CLCD_v01_1991_albert_henan.tif")
arr_1991 = arcpy.RasterToNumPyArray("F:\study\luoyang text\CLCD_v01_1992_albert_henan.tif")
arr_1991_1992 = np.concatenate((arr_1991, arr_1992),axis=0)
new_raster = arcpy.NumPyArrayToRaster(arr_1991_1992, arcpy.Point(x, y),
                                      cell_w, cell_h)
new_raster.save("F:\study\luoyang text\or_1991_1992.tif")  # 保存为一个新的栅格数据

但是结果不是我想要的，我只是在垂直方向按列进行了拼接，但我我需要在在空间维度进行拼接。

我需要的是得到一个图像。我在数组拼接应该出现了问题。

Answer 1

假定灰度像素值0～255表示黑到白共256级灰度
那么
新像素值=(int)((图1像素值+图2像素值)/2)