Threshold or somehow segment your images to get a mask of the "black" parts, then use
rgbImage = imoverlay(grayImage, mask, yourColor);
where yourColor is a 3 element array where each number of the fraction from black to bright for red, green, and blue, like [0.3, 0.6, 0.9] or whatever.
