Mini Club Astronomie-Python

[sebb916]

tiens je vais peut être faire quelque chose de ça :

[pejive]

Oui mais c'est moins performant et il ne faut pas avoir de gros doigts
viewtopic.php?f=34&t=5204
J'essaie de ne pas l'abandonner dans un fond de tiroir mais je m' éclate davantage sur un pc avec Pyzo

[sebb916]

je comprends, moi j'ai surtout utilisé le tableur sur le CX cas, en pleine manip avec les élèves, c'est plus facile à déplacer qu'un ordi ou même une tablette, mais les fonctionnalités restent limités et il faut de petits doigts, je suis d'accord

[pejive]

conversion d'une image RVB en niveaux de gris

image-py06.png (258.71 Kio) Consulté 3756 fois

[pejive]

Remplacement du bleu par un bleu synthétique réalisé avec 40% de rouge et 60% de vert
J'ai adapté un script proposé par ChatGPT
Il n'écrit pas les canaux tout à fait comme je l'ai fait au dessus; pour subplots également son écriture est plus compacte; va falloir creuser cette syntaxe

image-py07.png (83.98 Kio) Consulté 3748 fois

Figure_1BS.jpeg (27.58 Kio) Consulté 3748 fois

test sur une de mes anciennes images de M13 très chromatisée

fini les vilaines auréoles bleues !

M13sb.png (221.12 Kio) Consulté 3737 fois

Autre exemple

ngc1977sb.png (527.14 Kio) Consulté 3712 fois

[kaelig]

@pejive tu es un grand malade

[pejive]

Sur CloudyNight ils donnent plutôt B= 2*G-R pour le bleu synthétique (mais pour la Lune et les planètes)
Essai sur NGC1977

SB-2G-R.png (584.61 Kio) Consulté 3700 fois

pas idéal pour le fond du ciel
Essai avec M13

M13sb2gr.jpeg (55.15 Kio) Consulté 3695 fois

Par contre c'est effectivement meilleur sur Jupiter
(image initiale sans filtre au mak 127)

jupiter-sb.png (277.41 Kio) Consulté 3695 fois

Ma formule donne sensiblement le même résultat

jupiter-sb1.jpg (15.34 Kio) Consulté 3687 fois

[clouzot]

Sait-on jamais, s'il y a des utilisateurs de SkySafari dans la salle, voici un petit script Python qui permet de générer une image d'horizon à partir de mesures des parties du ciel visibles ou invisibles. Dans mon cas, coincé entre deux immeubles et sous le balcon du voisin, c'était indispensable...

Pour faire le relevé de points, une petite app en ligne utilisable sur n'importe quel smartphone sans installation : https://rkinnett.github.io/gyrocam/
Ca sort un fichier texte avec des couples azimuth - élévation qu'il suffit ensuite de recopier dans le script ci-dessous

Code : Tout sélectionner

import numpy as np
from PIL import Image, ImageDraw

# Image dimensions
WIDTH = 2048
HEIGHT = 1024

# Horizon measurement points (azimuth, elevation)
horizon_points = [
    (0, 89), (37, 89), (79, 89), (80, 54), (85, 50),
    (101, 40), (112, 25), (113, 17), (114, 0), (118, 0),
    (120, 0), (116, 0), (131, 0), (160, 0), (167, 0),
    (170, 0), (175, 0), (176, 0), (178, 15), (181, 25),
    (181, 30), (188, 35), (290, 89), (355, 89), (359, 89)
]


def az_el_to_pixel(az, el):
    """Convert azimuth/elevation to pixel coordinates"""
    # Azimuth: 0° at left edge, 360° at right edge
    x = (az / 360.0) * WIDTH

    # Elevation: +90° at top (y=0), 0° at middle (y=HEIGHT/2), -90° at bottom (y=HEIGHT)
    # y = HEIGHT/2 - (el / 90.0) * (HEIGHT / 2)
    y = HEIGHT / 2 * (1 - el / 90.0)

    return (x, y)


# Create RGBA image (transparent background)
img = Image.new('RGBA', (WIDTH, HEIGHT), (0, 0, 0, 0))
draw = ImageDraw.Draw(img)

# Convert horizon points to pixel coordinates
pixel_points = [az_el_to_pixel(az, el) for az, el in horizon_points]

# Add wrap-around point at the end (360° = 0°)
if horizon_points[0][0] == 0:
    wrap_point = (WIDTH, pixel_points[0][1])
    pixel_points.append(wrap_point)

# Create polygon points for the terrain (from horizon down to bottom of image)
polygon_points = []

# Add all horizon points
polygon_points.extend(pixel_points)

# Close the polygon by going along the bottom edge
polygon_points.append((WIDTH, HEIGHT))  # Bottom-right
polygon_points.append((0, HEIGHT))  # Bottom-left

# Draw the terrain as a filled polygon (opaque black)
draw.polygon(polygon_points, fill=(0, 0, 0, 255))

# Draw the horizon line for visibility (optional, can be removed)
draw.line(pixel_points, fill=(64, 64, 64, 255), width=2)

# Save the image
output_filename = 'skysafari_horizon.png'
img.save(output_filename, 'PNG')
print(f"Horizon image saved as '{output_filename}'")
print(f"Image size: {WIDTH}x{HEIGHT} pixels")
print(f"Format: PNG with alpha channel")
print(f"Number of horizon points: {len(horizon_points)}")

Résultat (vous pourrez vérifier que je n'y vois pas grand chose depuis mon balcon...)

[darki]

En fait @clouzot , tu n’as pas besoin de goto, tu pointes au milieu et tu attend que les objets passent, c’est plus simple

[clouzot]

En fait @clouzot , tu n’as pas besoin de goto, tu pointes au milieu et tu attend que les objets passent, c’est plus simple

C'est à peu près l'ambiance ici, oui. Heureusement qu'on a les lampadaires, ça rend les choses un peu plus pimentées

[pejive]

[pejive]

2e tentative (faut déjà comprendre comment et quand cliquer )
Az,El
99,7
104,34
79,42
95,19
109,14
128,8
150,3
je ne sais pas si Python en sortira quelque chose; avec un papier et un crayon çà doit être plus simple
A première vue c'est n'importe quoi

[pejive]

Premiers pas dans l'interface graphique tkinter
La présentation et la programmation sont à améliorer, mais çà semble marcher

my1stTkcalc.jpg (28.17 Kio) Consulté 1571 fois

my1stTkcalc1.jpg (24.98 Kio) Consulté 1568 fois

work in progress...

my1stTkcalc2.jpg (20.6 Kio) Consulté 1558 fois

[FredWeber]

Pour faire le relevé de points, une petite app en ligne utilisable sur n'importe quel smartphone sans installation

Super je vais tester et voir pour integrer ça dans kstars...

En tout cas c'est déja intégré dans le paramétrage, l'ajout d'un ficheir "horizon" sous forme de script texte avec les couples Az/El .... Il n'y a plus qu'à tester l'appli sur le smartphone.

Script Testé .... ça génère un fichier png au format equirectangulaire

Mais bizarrement, le Nord de Kstars n'est pas le même que celui de SkySafary