gdal_calc

python gdal_calc.py -A input.tif -B input2.tif --outfile=result.tif --calc="(A+B)/2"

python gdal_merge.py -o ouutfile input1.tif input2.tif

gdal_translate result.tif result_cut.tif -projwin 136 37 139 34