Skip to content

Production graphique nuit par nuit #3

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
PNPyrenees wants to merge 4 commits into
base: main
Choose a base branch
from
Open

Production graphique nuit par nuit #3

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
a0591f9
Ajout d'un paramètre dans le process_all permetant d'indiquer si on v…
PNPyrenees Jul 25, 2025
59b4334
Update README.md
PNPyrenees Jul 25, 2025
086bcc0
Suppression de print() superflus
PNPyrenees Jul 25, 2025
c8ccfec
Merge remote-tracking branch 'origin/main' into main
PNPyrenees Jul 25, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions .gitignore
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2,6 +2,9 @@
.Rhistory
.Rapp.history

# backup files
*.BAK

# Session Data files
.RData
.RDataTmp
Expand Down
8 changes: 6 additions & 2 deletions README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -111,6 +111,8 @@ sqm_mag_mod <- get_modal_sqm_mag_value(all_data)
# Génération des graphiques
generate_graph(best_night, NOM_SITE, "meilleures nuits", sqm_mag_mod)
generate_graph(all_data, NOM_SITE, "toutes les nuits", sqm_mag_mod)
# Génération d'un graphique de densité par nuit
generate_graph_per_nigth(all_data, NOM_SITE)
```

# Constantes
Expand All @@ -125,12 +127,14 @@ MIN_SQM_MAG_VAL <- 21
MAX_SQM_MAG_VAL <- 22.2
# Différence maximale entre sqm_mag min et max par nuit utilisée pour filtrer les données
DIFF_SQM_MAG <- 1
# Désactivation par défaut du calcul des graphiques de densité par nuit
GET_GRAPH_PER_NIGHT <- FALSE
```

Ces valeurs peuvent être changées lors de l'appel aux fonctions
```R
# Fonction globale
process_all(file_name, nom_site, sun_alt_min = SUN_ALT_MIN, diff_sqm_mag = DIFF_SQM_MAG, min_sqm_mag_val = MIN_SQM_MAG_VAL, max_sqm_mag_val = MAX_SQM_MAG_VAL)
# Fonction globale
process_all(file_name, nom_site, sun_alt_min = SUN_ALT_MIN, diff_sqm_mag = DIFF_SQM_MAG, min_sqm_mag_val = MIN_SQM_MAG_VAL, max_sqm_mag_val = MAX_SQM_MAG_VAL, get_graph_per_night = GET_GRAPH_PER_NIGHT)

# Fonction de calcul des "meilleures nuits"
get_best_night(data_in, sun_alt_min = SUN_ALT_MIN, diff_sqm_mag = DIFF_SQM_MAG, min_sqm_mag_val = MIN_SQM_MAG_VAL, max_sqm_mag_val = MAX_SQM_MAG_VAL)
Expand Down
63 changes: 60 additions & 3 deletions ninox_functions.R
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,8 +12,10 @@ MIN_SQM_MAG_VAL <- 21
MAX_SQM_MAG_VAL <- 22.2
# Différence maximale entre sqm_mag min et max par nuit utilisée pour filtrer les données
DIFF_SQM_MAG <- 1
# Désactivation par défaut du calcul des graphiques de ensité par nuit
GET_GRAPH_PER_NIGHT <- FALSE

process_all <- function(file_name, nom_site, format="csv", year=NULL, sun_alt_min = SUN_ALT_MIN, diff_sqm_mag = DIFF_SQM_MAG, min_sqm_mag_val = MIN_SQM_MAG_VAL, max_sqm_mag_val = MAX_SQM_MAG_VAL) {
process_all <- function(file_name, nom_site, format="csv", year=NULL, sun_alt_min = SUN_ALT_MIN, diff_sqm_mag = DIFF_SQM_MAG, min_sqm_mag_val = MIN_SQM_MAG_VAL, max_sqm_mag_val = MAX_SQM_MAG_VAL, get_graph_per_night = GET_GRAPH_PER_NIGHT) {
# Chargement et prétraitrement du fichier
if (format == "csv") {
all_data <- load_and_process_file_csv(file_name, year)
Expand All @@ -30,19 +32,29 @@ process_all <- function(file_name, nom_site, format="csv", year=NULL, sun_alt_mi
# Selection des nuits dans lune
data_without_moon <- subset(x = all_data, subset = moon_alt < 0)

print(sprintf("nb all_data : %s", nrow(all_data)))
print(sprintf("nb best_night : %s", nrow(best_night)))
print(sprintf("nb flat_night : %s", nrow(flat_night)))
print(sprintf("nb the_best_night : %s", nrow(the_best_night)))

# # #########################################
# # Rapport du lot de données
stats <- calculate_stats_data(all_data)
print_report(stats, file_name, nom_site)

# Génération des graphiques
generate_graph(best_night, nom_site, "", stats$sqm_mag_mod)
generate_graph(best_night, nom_site, "Meilleures nuits", stats$sqm_mag_mod)
generate_graph(all_data, nom_site, "All data", stats$sqm_mag_mod)
generate_graph_density_heatmap(data_without_moon, nom_site, "Nuits sans lune", stats$sqm_mag_mod)
generate_graph_density(the_best_night, nom_site, sprintf("%s", min(as.Date(the_best_night$date))))
generate_graph_density(the_best_night, nom_site, sprintf("La_meilleure_nuit-%s", min(as.Date(the_best_night$date))))

generate_graph_density_all_data(all_data, best_night, flat_night, the_best_night, nom_site )

# LLP : Ajout du traitement générant les graph pour chaque nuit
if (get_graph_per_night) {
generate_graph_per_nigth(all_data, nom_site)
}

}

print_report <- function(stats, file_name, nom_site) {
Expand Down Expand Up @@ -213,13 +225,18 @@ get_best_night <- function(data_in, nb_flat_day = NULL, nb_best_day = NULL, sun_

data_without_sun <- subset(x = data_in, subset = sun_alt < sun_alt_min)


# Ajout de la valeur absolue de la dérivée de sqm_mag (cia fonction diff) en fonction de la nuit
# LLP : Pour chaque ligne il calcul la différence entre le sqm_mag courant et celui de la ligne précédente -> la moyenne de ça par nuit donnant une idée du niveau de variabilité
data_without_sun <- transform(
data_without_sun,
abs_derive = ave(sqm_mag, julian_night, FUN = function(x) c(NA, abs(diff(x))))
)


agg_tbl <- generate_aggrate_per_night(data_without_sun)

write.csv(agg_tbl, "tmp.csv")
# Selection des nuits avec les "meilleurs" valeur :
# diff entre min et max < 1
# et valeurs comprises entre 21 et 22.8
Expand Down Expand Up @@ -247,6 +264,46 @@ get_best_night <- function(data_in, nb_flat_day = NULL, nb_best_day = NULL, sun_
return(best_night)
}


generate_graph_per_nigth <- function(data_in, nom_site) {
# #########################################
# GRAPHIQUE
# Graphique de densité pour chaque nuit

#Création du dossier de destination
dir.create("par_nuits", showWarnings = FALSE)

#Extraction des numéros de nuit unique
unique_julian_night <- unique(data_in$julian_night)

# Loop sur chaque nuit
nigth_idx <- 1
for (d in unique_julian_night) {
# Extraction des données correspondant au numéro de nuit
one_night_data <- data_in[data_in$julian_night == d, ]

sous_titre = sprintf("Nuit-%s", nigth_idx)
title_nb_nights = sprintf("\nDu %s/%s/%s au %s/%s/%s", first(one_night_data$d), first(one_night_data$m), first(one_night_data$y), last(one_night_data$d), last(one_night_data$m), last(one_night_data$y)) # Idélament, changer la valeur en "nuit du jj/MM/yyyy au jj/MM/yyyy"

plot <- ggplot(one_night_data) +
geom_line(aes(x = heure_graph, y = sqm_mag)) +
geom_point(aes(x = heure_graph, y = sqm_mag)) +
scale_y_reverse(breaks = seq(23,16,-1),limits=c(23,16)) +
ylab("NSB (magsqm/arsec²)") +
xlab("Time (UTC)") +
ggtitle(sprintf("NINOX %s %s %s", nom_site, sous_titre, title_nb_nights))

# Sauvegarde le graphique en fichier PNG
ggsave(plot,
filename = gsub(" ", "_", sprintf("par_nuits/%s_%s_densite.jpg", nom_site, sous_titre)),
device = "jpg",
height = 6, width = 5, units = "in"
)

nigth_idx <- nigth_idx + 1
}
}

generate_graph <- function(data_in, nom_site, sous_titre, valeur_modal) {
generate_graph_density(data_in, nom_site, sous_titre)
generate_graph_magnitude(data_in, nom_site, sous_titre, valeur_modal)
Expand Down