# Visualize an RGB image in all potential colorspaces as one glance

import os
import cv2
import numpy as np
from plantcv.plantcv import params
from plantcv.plantcv import resize
from plantcv.plantcv import plot_image
from plantcv.plantcv import fatal_error
from plantcv.plantcv import print_image
from plantcv.plantcv import rgb2gray_hsv
from plantcv.plantcv import rgb2gray_lab


def colorspaces(rgb_img):
    """ Visualize an RGB image in all potential colorspaces

    Inputs:
    rgb_img      = RGB image data

    Returns:
    plotting_img = Plotting image containing the original image and L,A,B,H,S, and V colorspaces

    :param segmented_img: numpy.ndarray
    :return labeled_img: numpy.ndarray

    """

    if not len(np.shape(rgb_img)) == 3:
        fatal_error("Input image is not RGB!")

    # Store and disable debug mode
    debug = params.debug
    params.debug = None

    # Initialize grayscale images list, rgb images list, plotting coordinates
    colorspace_names = ["H", "S", "V", "L", "A", "B"]
    all_colorspaces = []
    labeled_imgs = []
    y = int(np.shape(rgb_img)[0] / 2)
    x = int(np.shape(rgb_img)[1] / 2)

    # Loop through and create grayscale imgs from each colorspace
    for i in range(0, 3):
        channel = colorspace_names[i]
        all_colorspaces.append(rgb2gray_hsv(rgb_img=rgb_img, channel=channel))
    for i in range(3, 6):
        channel = colorspace_names[i]
        all_colorspaces.append(rgb2gray_lab(rgb_img=rgb_img, channel=channel))

    # Plot labels of each colorspace on the corresponding img
    for i, colorspace in enumerate(all_colorspaces):
        converted_img = cv2.cvtColor(colorspace, cv2.COLOR_GRAY2RGB)
        labeled = cv2.putText(img=converted_img, text=colorspace_names[i], org=(x, y),
                              fontFace=cv2.FONT_HERSHEY_SIMPLEX,
                              fontScale=params.text_size, color=(255, 0, 255), thickness=params.text_thickness)
        labeled_imgs.append(labeled)

    # Compile images together, including a larger version of the original image
    all_labeled = np.vstack([np.hstack([labeled_imgs[0], labeled_imgs[1], labeled_imgs[2]]),
                          np.hstack([labeled_imgs[3], labeled_imgs[4], labeled_imgs[5]])])
    plotting_img = resize(np.hstack([resize(img=rgb_img, resize_x=2, resize_y=2), all_labeled]), x=.5, y=.5)

    # Reset debug mode
    params.debug = debug

    if params.debug == "print":
        # If debug is print, save the image to a file
        print_image(plotting_img, os.path.join(params.debug_outdir, str(params.device) + "_vis_colorspaces.png"))
    elif params.debug == "plot":
        # If debug is plot, print to the plotting device
        plot_image(plotting_img)

    return plotting_img

@@ -398,19 +398,38 @@
 (rgb_img
+, original_img=True
 ):%0A
-
     %22%22%22 
@@ -526,16 +526,103 @@
 age data
+%0A    original_img = Whether or not to include the original image the the debugging plot
 %0A%0A    Re
@@ -722,24 +722,24 @@
 olorspaces%0A%0A
-
     :param s
@@ -766,16 +766,46 @@
 ndarray%0A
+    :param original_img: bool%0A
     :ret
@@ -2377,20 +2377,132 @@
 %5B5%5D%5D)%5D)%0A
+%0A
     
+# If original_img is True then also plot the original image with the rest of them%0A    if original_img:%0A        
 plotting
@@ -2508,23 +2508,16 @@
 g_img = 
-resize(
 np.hstac
@@ -2577,17 +2577,57 @@
 abeled%5D)
-,
+%0A    plotting_img = resize(plotting_img, 
  x=.5, y

add option to exclude ori image in colorspaces function