vmin=np.nanpercentile(hcop, 10),
#vmax=np.nanpercentile(hcop,99.0),
vmax=0.06,
clip=True)(hcop)
rgb_hcop = np.zeros(hc3n.shape + (3, ))
rgb_hcop[:, :, blue] = np.nan_to_num(monochrome_hcop)
hsv_hcop = rgb_to_hsv(rgb_hcop)
hsv_hcop[:, :, hue] = 0 / 360.
rgb_hcop = hsv_to_rgb(hsv_hcop)
rgb_im[:, :, :alpha] += rgb_hcop
vmin = 0.0002
monochrome_cont = LogNorm(vmin=vmin,
vmax=np.nanpercentile(cont, 99.9995),
clip=True)(cont * (cont > vmin))
monochrome_cont = monochrome_cont.filled(np.nan)
# vmin = 0.0001
# monochrome_cont = Normalize(vmin=vmin, vmax=np.nanpercentile(cont,99.9),
# clip=True)(cont*(cont > vmin))
#monochrome_cont[monochrome_cont.mask] = 0.0
#monochrome_cont.mask[:] = False
# rgb_cont = monochrome_cont[:,:,None]
# rgb_cont = np.zeros(hc3n.shape + (3,))
# rgb_cont[:,:,blue] = np.nan_to_num(monochrome_cont)
# hsv_cont = rgb_to_hsv(rgb_cont)
# hue, saturation, value = 0,1,2
# hsv_cont[:,:,hue] = 150/360.
# rgb_cont = hsv_to_rgb(hsv_cont)
# rgb_im[:,:,:alpha] += rgb_cont
#rgb_im[:,:,:alpha] += np.nan_to_num(monochrome_cont[:,:,None])
nir_im = nir_hdu[0].data
# BEGIN IMAGE MAKING HERE
red, green, blue, alpha = 0, 1, 2, 3
hue, saturation, value = 0, 1, 2
# FIRST LAYER: NIR = white
rgb_im = np.zeros(nir_im.shape + (4, ))
rgb_im[:, :, alpha] = 1.0
nir_im_rezero = nir_im - np.nanmin(nir_im)
colorized_nir = LogNorm(vmin=np.nanpercentile(nir_im_rezero, 30),
vmax=np.nanpercentile(nir_im_rezero, 99.9995),
clip=True)(nir_im_rezero)
colorized_nir.fill_value = 0.0
rgb_nir_im = np.nan_to_num(colorized_nir.filled())[:, :, None] * np.ones(3)
rgb_im[:, :, :alpha] += rgb_nir_im
cm_im_rezero = cm_im - np.nanmin(cm_im)
monochrome_cm_im = LogNorm(
vmin=0.0003, #np.nanpercentile(cm_im-np.nanmin(cm_im),30),
vmax=np.nanpercentile(cm_im_rezero, 99.9995),
clip=True)(cm_im_rezero)
monochrome_cm_im.fill_value = 0.0
monochrome_cm_im = monochrome_cm_im.filled()
#rgb_cm_im = np.nan_to_num(monochrome_cm_im)
#rgb_im[:,:,red] += rgb_cm_im
rgb_cm_im = np.zeros(nir_im.shape + (3, ))
rgb_cm_im[:, :, blue] = np.nan_to_num(monochrome_cm_im)
#rgb_cm_im[:,:,blue] = monochrome_cm_im.filled()
