Wiki de Martin Aubé | A11Cours02GroupeAEquipe02 / Methof1

Comme expliqué précédemment, nous utiliserons deux techniques distinctes et novatrices afin de mesurer la pollution lumineuse. La première méthode repose sur l'usage du spectromètre, tandis que la seconde nécessite le radiomètre.

Mesure de la pollution lumineuse avec le spectromètre

D'abord, le spectromètre récolte des données durant un temps de pose de quatre heures. On génère ainsi, un spectre illustrant la luminance énergétique par nanomètre (W×m-2×sr-1×nm-1) en fonction de la longueur d'onde (nm). L'ensemble de l'environnement du spectromètre ainsi que ses composantes électroniques dégagent de la chaleur. Cette énergie thermique influence les valeurs enregistrées par les capteurs CCD, puisqu'une température chaude va fausser la luminance à la hausse. Afin de supprimer ce bruit thermique, on enregistre un spectre durant quatre heures avec l'obturateur fermé. Puisqu'aucune lumière ne parvient au capteur CCD, les données enregistrées par le capteur correspondent uniquement à l'agitation thermique. À partir de là, on soustrait ce spectre de celui obtenu avec le ciel, de sorte qu'on élimine en grande partie le bruit thermique. Aussi, il est important de s'assurer que le ciel était sans nuages et sans Lune durant l'enregistrement des données. Afin de s'en assurer, on regarde l'aspect général du spectre. Si on remarque que le spectre épouse une courbe en cloche telle qu'à la figure 1a, on sait que la Lune était présente durant le captage des données. De plus, si le spectre a le même aspect que celui de la figure 1boù est cette figure?, on sait que le ciel était trop nuageux. Dans ces deux situations, on doit éliminer les spectres, puisque les nuages et la Lune empêchent l'analyse de nos résultats. Préalablement à nos prises de mesure du ciel, on a produit le spectre des lampes les plus utilisées la nuit dans la ville de Sherbrooke. On a donc le spectre d'émission d'une lampe à halogénure métallique (ML) et d'une lampe au sodium à haute pression (HPS). On peut trouver les pics d'émission les plus importants de ces lampes et noter leur longueur d'onde caractéristique. À partir de là, on peut retrouver ces pics sur les spectres du ciel. La figure 2 présente les spectres des lampes ML et HPS tandis que le tableau 1 décrit certaines raies facilement visibles sur les spectres du ciel (indiquer la longueur d'onde et la LE c.-à-d. l'aire sous la courbe du pic). Avec le programme Fityk, on fait une analyse rapide et efficace de nos spectres. D'abord, on supprime les pics très pointus causés par les rayons cosmiques. Aussi, on ajuste de larges courbes gaussiennes à notre nuage de points, afin de modéliser le rayonnement de corps noir des étoiles. Enfin, on modélise chaque raie spectrale à l'aide d'une courbe gaussienne. L'aire sous la courbe de ces raies indique la luminance énergétique (W×m-2×sr-1). Sur le spectre du ciel, on mesure la luminance des pics se trouvant aux longueurs d'onde du tableau 1, ces données sont notées LE'. Comme énoncé dans le cadre théorique, un rayon de luminance énergétique LE et de longueur d'onde λ pointé vers le ciel devrait en partie réfléchir vers le sol selon une loi exponentielle de la forme LE' = LE × cλα. À partir des données du tableau 1 (LE) et des spectres du ciel (LE'), on peut déterminer les meilleurs coefficients c et α. En multipliant chaque spectre des lampes par les facteurs de correction cλα, on obtient la luminance énergétique réfléchie par le ciel. En intégrant ces spectres transformés, on obtient une bonne mesure de la pollution lumineuse notée PLSAND.

Mesure de la pollution lumineuse avec le radiomètre

Le radiomètre nous permet d'obtenir des images du ciel nocturne. Celles-ci sont obtenues avec un temps de pose de 100 secondes. Comme dans le cas du spectromètre, on doit d'abord s'assurer que le ciel était sans nuages et sans Lune au moment de la pose. Cette vérification est faite de façon visuelle avec les images du ciel. La figure 3a présente un ciel nuageux, tandis que la figure 3b présente une image avec la Lune. Aussi, comme avec le spectromètre, il faut supprimer le bruit thermique de l'image. Pour ce faire, on prend une image avec l'obturateur fermé (un Dark), puis on supprime cette image de celles du ciel. On élimine ainsi le bruit thermique. Dans l'analyse, nous sélectionnons les endroits du ciel les plus sombres, puisque c'est dans ces zones que l'importance relative de la pollution lumineuse est la plus grande. Sur l'image, on ne retient donc que les 50 pixels qui ont la plus faible intensité lumineuse. Cependant, certains pixels ont une intensité lumineuse anormalement basse. Ces aberrations sont causées par des rayons cosmiques ayant percuté le capteur CCD durant la prise du Dark. Lorsqu'un de ces rayons de grande intensité traverse l'obturateur et frappe le CCD, des hautes intensités seront mesurées sur les pixels atteints. Lorsque le Dark est soustrait de l'image, on obtient ainsi des intensités anormalement petites sur certains pixels. Afin de ne pas tenir compte de ces points erronés, ils sont systématiquement exclus de notre échantillon des 50 pixels les plus sombres. Pour obtenir cet échantillon, nous appliquons un algorithme programmé en Fortran. Les étapes du processus sont définies ci-dessous :

On fait d'abord l'histogramme de l'intensité lumineuse des pixels (on exprime la fréquence f en fonction de l'intensité i pour obtenir une fonction notée f(i)), comme à la figure 4.
Depuis l'intensité i = 1, on parcourt l'histogramme jusqu'à une intensité ayant une fréquence non nulle.
À partir de cette intensité, que l'on note k, on fait la somme de sa fréquence et avec celles des 20 intensités suivantes (f(k) + f(k+1) + ... + f(k+20) = S).
Si cette somme est égale à la fréquence de l'intensité k (si S = f(k)), cela signifie que la fréquence des 20 intensités suivantes est nulle. Dans ce cas, l'intensité k ainsi que les 100 intensités précédentes sont exclues (on pose {#f(k-100 \gt i \le k) = 0)#}, car elles sont marginales. Au contraire, si cette somme est supérieure à la fréquence de l'intensité k, cela signifie que l'intensité k est proche d'au moins une autre valeur. L'intensité k est alors conservée.
Ensuite, on parcourt l'histogramme jusqu'à une autre intensité ayant une fréquence non nulle.
On répète les étapes 3, 4 et 5 jusqu'à la fin de l'histogramme.
Pour finir, on fait la moyenne des 50 intensités les plus basses de l'histogramme corrigé. Cette moyenne est notée Ix (où x est le nom du filtre).

L'intensité moyenne de ces pixels est la luminance énergétique des zones du ciel sans étoiles (W×m-2×sr-1). Puisque les différents filtres ne laissent pas passer la lumière aux mêmes longueurs d'onde, les valeurs de Ix varient. Une fois qu'on a obtenu Ix pour chaque filtre, on utilise un programme permettant de trouver une combinaison linéaire entre tous les Ix et chaque valeur de référence PLSAND. On peut écrire cette combinaison linéaire comme suit :

PLSAND ≈ a IClear + b IIR + c IHα + d IDS + e IComet

La combinaison linéaire est donc une fonction à plusieurs variables ayant pour image les valeurs PLSAND enregistrées dans différentes zones du ciel. Notre programme permettra de trouver les meilleurs coefficients a, b, c, d et e permettant d'optimiser le coefficient de corrélation linéaire de cette relation. Une fois que les meilleurs coefficients ont été trouvés pour notre combinaison, on fait une analyse statistique pour déterminer la précision de notre approximation en évaluant la variabilité des résultats sur plusieurs nuits. On évalue aussi son exactitude en comparant notre approximation aux données du spectromètre (PLSAND).