Les images de James Webb

Créer sa propre image

Un exemple pour mettre en couleurs les images du télescope James Webb

Utiliser les fonctions filtres

import matplotlib.pyplot as plt
from PIL import Image
import numpy as np

#Import de l'image:
img = Image.open("cat_nvg.png")

#Création d'une nouvelle image vide à la bonne taille:
new_img = Image.new("I", img.size)

#Numpy 2D array :
data = np.array(img.getdata())

#Formatage 0-1
data = data/max(data)

#Fonctions filtres:
def resserrer(x):
    return 4*(x-0.5)**3+0.5

def racine(x):
    return np.sqrt(x)

def carre(x):
    return x*x

def contraste(x):
    return x*(-2*x**2+3*x)

def renverser(x):
    return 1-x

#Vos modifications ici:
new_data = contraste(data)

#Remplissage de l'ancienne image et de la nouvelle en 0-255:
img.putdata(255*data)
new_img.putdata(255*new_data)

#Affichage
fig = plt.figure(figsize=(7, 5))
fig.add_subplot(2, 2, 1)
plt.imshow(img)
plt.title("Originale")

fig.add_subplot(2, 2, 3)
plt.imshow(new_img)
plt.title("Contraste")

#Pour le prof
fig.add_subplot(2, 2, 2)
t = np.arange(0., 1., 0.1)
plt.plot(t, t, 'r--')
plt.axis('square')

fig.add_subplot(2, 2, 4)
t = np.arange(0., 1., 0.1)
plt.plot(t, contraste(t), 'r--', t, t, 'b--')
plt.axis('square')

plt.show()

Modifier votre Nébuleuse de la Carène ou vos Piliers de la Création

Modifier la ligne 65 ou la ligne 47 par exemple...

import matplotlib.pyplot as plt
from PIL import Image
import numpy as np

#Import des images filtres:
R_img = Image.open("Carene_RED_small.png")
G_img = Image.open("Carene_GREEN_small.png")
B_img = Image.open("Carene_BLUE_small.png")

#R_img = Image.open("Piliers_RED_small_i.png")
#G_img = Image.open("Piliers_GREEN_small_i.png")
#B_img = Image.open("Piliers_BLUE_small_i.png")

#Création de la nouvelle image couleur vide:
RGB_img = Image.new("RGB", R_img.size)

#Numpy 2D array :
R_data = np.array(R_img.getdata())
G_data = np.array(G_img.getdata())
B_data = np.array(B_img.getdata())

#Formatage 0-1
R_data = R_data/max(R_data)
G_data = G_data/max(G_data)
B_data = B_data/max(B_data)

#Fonctions filtres:
def resserrer(x):
    return 4*(x-0.5)**3+0.5

def racine(x):
    return np.sqrt(x)

def carre(x):
    return x*x

def contraste(x):
    return x*(-2*x**2+3*x)

def renverser(x):
    return 1-x

#Modifications:
#R_data = carre(R_data)
#R_data = racine(R_data)
#G_data = resserrer(G_data)
B_data = contraste(racine(B_data))

#Remplissage des images filtres modifiées sur 0-255:
R_img.putdata(255*R_data)
G_img.putdata(255*G_data)
B_img.putdata(255*B_data)

#Remplissage de l'image couleur:
#On parcourt tous les pixels de l'image RGB:
for x in range(RGB_img.size[0]):
    for y in range(RGB_img.size[1]):
        #On recupère les valeurs de pixels de chaque image filtre:
        R = R_img.getpixel((x, y))
        G = G_img.getpixel((x, y))
        B = B_img.getpixel((x, y))
        if R > 4*B:
            B = R
        #pix = (R, G, B)
        pix = ((6*R+2*G+B)//9, (R+3*G+B)//5, (G+2*B+G)//4)
        #On ajoute le pixel dans la nouvelle image
        RGB_img.putpixel((x, y), pix)
plt.imshow(RGB_img)
plt.show()

Cours1 : Déplacer la tortue Cours2 : Répéter Cours3 : Procédures Cours4 : On mélange tout