Fixing BBA_empirical coefficient

GEUS-SICE · Oct 18, 2021 · d4ace50 · d4ace50
1 parent 9b50036
commit d4ace50
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/dm.grass.sh b/dm.grass.sh
@@ -176,6 +176,6 @@ parallel "r.out.gdal -m -c input={} output=${outfolder}/${date}/{}.tif ${tifopts
 # Generate some extra rasters
 tifopts='type=Float32 createopt=COMPRESS=DEFLATE,PREDICTOR=2,TILED=YES --q --o'
 r.mapcalc "ndbi = ( r_TOA_01 - r_TOA_21 ) / ( r_TOA_01 + r_TOA_21 )"
-r.mapcalc "bba_emp = (r_TOA_01 + r_TOA_06 + r_TOA_17 + r_TOA_21) / 4.0 * 0.945 + 0.055"
+r.mapcalc "bba_emp = (r_TOA_01 + r_TOA_06 + r_TOA_17 + r_TOA_21) / 4.0 * 1.003 + 0.058" 
 r.out.gdal -f -m -c input=ndbi output=${outfolder}/${date}/NDBI.tif ${tifopts}
 r.out.gdal -f -m -c input=bba_emp output=${outfolder}/${date}/BBA_emp.tif ${tifopts}