def return_johnson_uvk_img(filename,
camera=0,
scale_min=0.1,
scale_max=50,
size_scale=1.0):
if (not os.path.exists(filename)):
print(("file not found:", filename))
sys.exit()
b_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(
filename, band="U_Johnson.res")
g_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(
filename, band="V_Johnson.res")
r_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(
filename, band="K_Johnson.res")
b_image = sunpy__load.load_broadband_image(
filename, band='U_Johnson.res',
camera=camera) * b_effective_wavelength / g_effective_wavelength * 2.5
g_image = sunpy__load.load_broadband_image(
filename, band='V_Johnson.res',
camera=camera) * g_effective_wavelength / g_effective_wavelength
r_image = sunpy__load.load_broadband_image(
filename, band='K_Johnson.res',
camera=camera) * r_effective_wavelength / g_effective_wavelength * 1.5
n_pixels
img = np.zeros((n_pixels, n_pixels, 3), dtype=float)
img[:, :, 0] = asinh(r_image,
scale_min=scale_min,
scale_max=scale_max,
non_linear=0.5)
img[:, :, 1] = asinh(g_image,
scale_min=scale_min,
scale_max=scale_max,
non_linear=0.5)
img[:, :, 2] = asinh(b_image,
scale_min=scale_min,
scale_max=scale_max,
non_linear=0.5)
img[img < 0] = 0
return img
def return_h_band_img(filename,camera=0,scale_min=0.1,scale_max=50,size_scale=1.0):
if (not os.path.exists(filename)):
print "file not found:", filename
sys.exit()
image = sunpy__load.load_broadband_image(filename,band='H_Johnson.res', camera=camera)
n_pixels = image.shape[0]
img = np.zeros((n_pixels, n_pixels), dtype=float)
img[:,:] = asinh(image, scale_min=scale_min, scale_max=scale_max,non_linear=0.5)
img[img<0] = 0
return img
def return_sdss_gri_img(filename,camera=0,scale_min=0.1,scale_max=50,size_scale=1.0, non_linear=0.5):
if (not os.path.exists(filename)):
print "file not found:", filename
sys.exit()
b_image = sunpy__load.load_broadband_image(filename,band='g_SDSS.res',camera=camera) * 0.7
g_image = sunpy__load.load_broadband_image(filename,band='r_SDSS.res',camera=camera) * 1.0
r_image = sunpy__load.load_broadband_image(filename,band='i_SDSS.res',camera=camera) * 1.4
n_pixels = r_image.shape[0]
img = np.zeros((n_pixels, n_pixels, 3), dtype=float)
img[:,:,0] = asinh(r_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[:,:,1] = asinh(g_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[:,:,2] = asinh(b_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[img<0] = 0
del b_image, g_image, r_image
gc.collect()
return img
def return_johnson_uvk_img(filename,camera=0,scale_min=0.1,scale_max=50,size_scale=1.0):
if (not os.path.exists(filename)):
print "file not found:", filename
sys.exit()
b_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(filename,band="U_Johnson.res")
g_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(filename,band="V_Johnson.res")
r_effective_wavelength = sunpy__load.load_broadband_effective_wavelengths(filename,band="K_Johnson.res")
b_image = sunpy__load.load_broadband_image(filename,band='U_Johnson.res',camera=camera) * b_effective_wavelength / g_effective_wavelength * 2.5
g_image = sunpy__load.load_broadband_image(filename,band='V_Johnson.res',camera=camera) * g_effective_wavelength / g_effective_wavelength
r_image = sunpy__load.load_broadband_image(filename,band='K_Johnson.res',camera=camera) * r_effective_wavelength / g_effective_wavelength * 1.5
n_pixels
img = np.zeros((n_pixels, n_pixels, 3), dtype=float)
img[:,:,0] = img_scale.asinh(r_image, scale_min=scale_min, scale_max=scale_max,non_linear=0.5)
img[:,:,1] = img_scale.asinh(g_image, scale_min=scale_min, scale_max=scale_max,non_linear=0.5)
img[:,:,2] = img_scale.asinh(b_image, scale_min=scale_min, scale_max=scale_max,non_linear=0.5)
img[img<0] = 0
return img
def return_h_band_img(filename,camera=0,scale_min=0.1,scale_max=50,size_scale=1.0):
if (not os.path.exists(filename)):
print "file not found:", filename
sys.exit()
image = sunpy__load.load_broadband_image(filename,band='H_Johnson.res', camera=camera)
n_pixels = image.shape[0]
img = np.zeros((n_pixels, n_pixels), dtype=float)
img[:,:] = img_scale.asinh(image, scale_min=scale_min, scale_max=scale_max,non_linear=0.5)
img[img<0] = 0
return img
def return_sdss_gri_img(filename,camera=0,scale_min=0.1,scale_max=50,size_scale=1.0, non_linear=0.5):
if (not os.path.exists(filename)):
print "file not found:", filename
sys.exit()
b_image = sunpy__load.load_broadband_image(filename,band='g_SDSS.res',camera=camera) * 0.7
g_image = sunpy__load.load_broadband_image(filename,band='r_SDSS.res',camera=camera) * 1.0
r_image = sunpy__load.load_broadband_image(filename,band='i_SDSS.res',camera=camera) * 1.4
n_pixels = r_image.shape[0]
img = np.zeros((n_pixels, n_pixels, 3), dtype=float)
img[:,:,0] = img_scale.asinh(r_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[:,:,1] = img_scale.asinh(g_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[:,:,2] = img_scale.asinh(b_image, scale_min=scale_min, scale_max=scale_max,non_linear=non_linear)
img[img<0] = 0
del b_image, g_image, r_image
gc.collect()
return img