grdimage assigns the colors to each cell but the radiances/shadings are computed in grdgradient. The RGB color is converted to HSV where shadings are applied and the result is converted back to RGB. This is a different procedure from gdal_calc.py that blends images … derived from gradients.