I'm trying to replicate the stretched color scheme employed by ArcGIS for displaying rasters in R. I want to use the method that defines the min/max of the color scale as mean(raster_values) +/- 2*sd(raster_values), and then anything outside this range will get set to the colors of the min or max depending upon the inequality.I'm trying to replicate the stretched color sch