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()