Grâce à la spectroscopie des vitesses radiales, on connaît en 2003 plus d'une centaines d'exoplanètes géantes autour d'étoiles de type solaire. Les limites de cette technique poussent à utiliser d'autres méthodes pour caractériser complètement les planètes connues, détecter des planètes de taille terrestre, et accéder à des informations statistiques sur une grande population d'objets. Cette thèse porte sur deux de ces techniques complémentaires : d'une part la photométrie des transits et d'autre part l'interférométrie infrarouge à très longue base. La photométrie des transits a été étudiée dans le cadre de la mission spatiale COROT : on a mis en oeuvre un algorithme de détection par filtrage adapté afin de déterminer (i) l'efficacité de détection de l'instrument pour différentes distributions orbitales des planètes, et (ii) la répartition des détections en fonction du type spectral et de la magnitude de l'étoile hôte. En outre, on a étudié l'utilisation des canaux photométriques de COROT pour discriminer les fausses détections et pour réduire la part du bruit de variabilité stellaire. Un traitement global des mesures par interférométrie à très longue base et par spectroscopie permet la caractérisation complète des systèmes doubles. On a appliqué ce traitement à des observations réelles d'étoiles binaires spectroscopiques, puis on a étudié la faisabilité d'observations de couples du type de 51 Pegasi. Par ailleurs, on a contribué à l'effort vers la très haute précision en interférométrie infrarouge par la réalisation d'un catalogue d'étoiles-étalons, et par la caractérisation instrumentale d'un prototype de fibre monomode à 10 µm, un composant essentiel pour les missions spatiales dédiées aux exoplanètes telluriques comme Darwin/TPF.

L'utilisation de collections d'images pour les applications de vision par ordina- teur devient de plus en plus commune des nos jours. L'objectif principal de cette theÌse est d'exploiter et d'extraire des informations importantes d'images de sceÌnes d'exteìrieur aÌ partir de ce type de collections : l'illumination preìsente au moment de la prise, les proprieìteìs de reìflectance des mateìriaux composant les objets dans la sceìne et les proprieìteìs radiomeìtriques des appareils photo utiliseìs. Pour atteindre notre objectif, cette theÌse et composeì de deux parties principales. Dans un pre- mier temps nous allons reìaliser une analyse de diffeìrentes repreìsentations du ciel et une com- paraison des images baseìe sur l'apparence de celui-ci. Une grande partie de l'information visuelle perçue dans les images d'exteìrieures est due aÌ l'illumination en provenance du ciel. Ce facteur est repreìsenteì par les rayons du soleil reìfleìchis et reìfracteìs dans l'atmospheÌre en creìant une illumination globale de l'environnement. En même temps cet environnement deìtermine la façon de percevoir les objets du monde reìel. Eìtant donneì l'importance du ciel comme source d'illumination, nous formulons un processus geìneìrique en trois temps, segmentation, modeìlisation et comparaison des pixels du ciel, pour trouver des images similaires en se basant sur leurs apparences. Diffeìrentes meìthodes sont adopteìes dans les phases de modeìlisation et de comparaison. La performance des algorithmes est valideìe en trouvant des images similaires dans de grandes collections de photos. La deuxieÌme partie de cette theÌse consiste aÌ exploiter l'information geìomeìtrique addition- nelle pour en deìduire les caracteìristiques photomeìtriques de la sceÌne. A partir d'une structure 3D reìcupeìreìe en utilisant des meìthodes disponibles, nous analysons le processus de forma- tion de l'image aÌ partir de modeÌles simples, puis nous estimons les parameÌtres qui les reìgis- sent. Les collections de photos sont geìneìralement captureìes par diffeìrents appareils photos, d'ouÌ l'importance d'insister sur leur calibrage radiomeÌtrique. Notre formulation estime cet eìtalonnage pour tous les appareils photos en même temps, en utilisant une connaissance a priori sur l'espace des fonctions de reìponse des cameìras possibles. Nous proposons ensuite, un cadre d'estimation conjoint pour calculer une repreìsentation de l'illumination globale dans chaque image, l'albedo de la surface qui compose la structure 3D et le calibrage ra- diomeìtrique pour tous les appareils photos.

Dans l'ouvrage [DBB06], Dutré et al. énumèrent dix problèmes d'Illumination Globale non résolus. Parmi ceux-ci, huit sont liés à la mesure et à l'intégration de données réelles dans les simulations. Cette thèse, en quatre parties et trois annexes, s'intéresse à cette problématique. Tout d'abord, les simulations de la propagation de la lumière dans les domaines du rendu physiquement réaliste, de l'éclairage, de la télédétection, de la conception de systèmes optiques... sont envisagées comme résolutions numériques d'un problème d'optique, fédérées par la notion de " photosimulation ". Dans le cadre de l'équation de rendu, les modèles réalistes de sources et le principe de la mesure goniophotométrique en champ lointain sont alors rappelés. La représentation des interactions lumière-matière est ensuite introduite par une exploration de l'apparence macroscopique, qui amène au rappel de la Fonction de Distribution de la Réflectance Bidirectionnelle et Spectrale (SBRDF), de ses principales propriétés et modèles. Le problème de la mesure pratique de la SBRDF d'une surface dans le visible est ensuite exploré. Id est, une taxonomie des méthodes de la littérature est établie ; qui allègue en faveur du développement de nouvelles approches. Un dispositif innovant, multispectral imageant, est alors présenté. Il se fonde sur la capture de la projection de la SBRDF sur un écran lambertien 3/4-sphérique, par une caméra multispectrale grand angle, assemblage d'un objectif fisheye, d'un filtre dynamique LCTF et d'une caméra CCD 12 bits. L'extraction des images capturées de l'information de la SBRDF repose sur un modéle radiométrique, qui explicite la transformation de la lumière en niveaux des pixels, dans le formalisme physique. Ce modèle soulève des problèmes de reconstruction multispectrale et d'interréflexions, pour lesquels de nouveaux algorithmes de résolution sont implantés. Les mesures de SBRDF produites semblent très prometteuses. Dans la troisième partie, le problème de la reconstruction d'une fonction directionnelle, identifié comme fondamental en photosimulation, est traité dans le cadre de la reconstruction de SBRDF discrètes. Pour cela, les propriétés mathématiques souhaitables des fonctions reconstruites sont envisagées. Puis, à l'aune de ce corpus, les approches de la littérature sont discutées ; justifiant la recherche d'algorithmes plus performants. Une nouvelle méthode est alors proposée, fondée sur une double triangulation sphérique des échantillons, et une généralisation à la surface d'un triangle sphérique de l'interpolant spline cubique de Hermite. La fonction reconstruite est interpolante, dérivable, quasi-déterministe, ne présente pas l'artéfact bidirectionnel, intègre la métrique sphérique non euclidienne, et prend en compte le difficile problème du masquage. Qualitativement comme quantitativement, les résultats obtenus semblent plaider en faveur du surcroît de complexité théorique qu'induit cette approche. Dans la dernière partie, cet algorithme de reconstruction est appliqué au problème de la " photométrie en champ proche ", ou, la caractérisation d'une source réelle par un solide photométrique étendu. La supériorité théorique du modèle étendu est d'abord démontrée. Puis, un état de l'art de la photométrie en champ proche est réalisé ; justifiant la mise en œuvre d'une nouvelle approche. Un nouveau dispositif est alors présenté. Il repose sur le déplacement d'un vidéoluminancemètre, couplage original d'une caméra CCD 12 bits avec un luxmètre, le long d'une trajectoire hémisphérique relative à la source. Les procédures de calibrage de l'instrument - géométrique, radiométrique, et photométrique - sont explicitées. Les sources de lumière caractérisées par cette approche ont été jugées suffisamment probantes pour être intégrées aux photosimulations spectrales du CSTB

Les petits corps sont des clés pour comprendre notre système solaire. L'étude de cette population nous donne en effet accès aux informations sur l'état et sur la structure du système solaire primordial et du système solaire actuel, ainsi que sur son évolution et sur les processus de formation des planètes. Connaître la composition de surface des petits corps nous fournit des ingrédients et des proportions pour cette recette cosmique. Cette thèse, qui inclut l'étude des petits corps du système solaire interne et externe, est dédiée à la compréhension de la tendance compositionnelle des corps à travers le système solaire en utilisant des mesures photométriques et spectroscopiques. Je présente une classification (taxonomie) dans les longueurs d'ondes du visible et du proche infrarouge (de 0.4 à 2.4 µm), basée sur les données spectrales de 371 astéroïdes. Cette taxonomie comprend 24 classes qui chacune caractérise au mieux les variations spectrales observées parmi les petits corps du système solaire interne. De part la création de cette taxonomie, nous apprenons qu'en analysant les données dans les longueurs d'ondes du visible uniquement, il reste des incertitudes sur la forme de la bande d'absorption à 1 micron. Bien que la gamme de longueur d'onde du proche infrarouge soit excellente pour interpréter les données incluant les bandes diagnostiques à 1 et 2 microns, les complexes C et X des spectres sans fortes bandes paraissent plutôt dégénérés dans ce régime. J'analyse les couleurs photométriques des 23 objets trans-neptuniens (OTN) et Centaures, parmi lesquels neuf n'avaient jamais été observés précédemment, et je leur assigne une classe taxonomique. Je discute des objets qui ont soit changé de classe depuis les données préalables soit changé considérablement de magnitude absolue. De plus, j'interprète la composition de surfaces de trois petits corps du système solaire externe, l'objet couplé avec Jupiter (52872) Okyrhoe et les OTNs (90482) Orcus et (73480) 2002 PN34 en modélisant des mesures spectroscopiques dans les gammes du visible et du proche infrarouge. Les spectres révèlent des variations de quantité de glace d'eau à la surface de ces corps. Pour Orcus j'apporte des contraintes approximatives sur la présence de matériaux plus volatiles que la glace d'eau. Ensuite, je présente une recherche de l'éthane solide, C2H6, sur les surfaces de Pluton et de Triton. Celle-ci est basée sur les observations spectrales dans les longueurs d'ondes du proche infrarouge. Je modélise chaque surface en utilisant un modèle de transfert radiatif fondé sur la théorie de Hapke (Hapke, 1993) de trois manières : sans éthane, avec de l'éthane pur, et avec de l'éthane dilué dans de l'azote. La présence de moins de quelques pourcents d'éthane sur chaque corps ne permet pas d'exclure ce composant de Triton et Pluton, cependant il n'y a pas non plus de forte détection. Finalement, je reconsidère la connaissance actuelle de la distribution compositionnelle des matériaux de notre système solaire en fournissant une vue globale des petits corps. Je me concentre particulièrement sur la présence de l'eau dans toutes ses phases qui est pertinente surtout pour notre propre planète, la Terre, et la vie. Je compare brièvement la structure générale de notre système solaire aux autres disques d'accrétion, afin de mettre en perspective la vue détaillée mais cependant étroite de notre système solaire avec celle, plus large mais à basse résolution, des autres systèmes planétaires.

Cette thèse, effectuée à l'Institut de Physique Nucléaire de Lyon (IPNL), s'inscrit dans le cadre de l'expérience SNfactory, dédiée à l'observation précise d'un grand nombre de supernovae de type Ia. Ces explosions thermonucléaires de naines blanches, très lumineuses, sont aujourd'hui les meilleures chandelles standards connues et ont de ce fait un grand intérêt pour la cosmologie.<br />Les travaux exposés ici contribuent à l'étalonnage de SNIFS, le spectrographe à champ intégral utilisé par SNfactory pour l'observation des supernovae de type Ia. Une étude des CCD a été effectuée ainsi que la collection de spectres d'étoiles de référence pour l'étalonnage en flux de l'instrument. La contribution principale a consisté en la réalisation d'un code complet pour l'extraction du spectre d'une source ponctuelle à partir du cube pré-étalonné de données. Le procédé choisi est la photométrie de PSF (Point Spread Function), nécessaire ici du fait de la petite taille du champ observé. Le flux perdu peut ainsi être évalué, et l'on atteint une qualité de mesure du flux de l'ordre du pourcent, objectif de l'expérience.

L'utilisation de collections d'images pour les applications de vision par ordinateur devient de plus en plus commune des nos jours. L'objectif principal de cette thèse est d'exploiter et d'extraire des informations importantes d'images de scènes d'extérieur a partir de ce type de collections : l'illumination présente au moment de la prise, les propriétés de reflectance des matériaux composant les objets dans la scène et les propriétés radiométriques des appareils photo utilisés. Pour atteindre notre objectif, cette thèse est composée de deux parties principales. Dans un premier temps nous allons réaliser une analyse de différentes représentations du ciel et une comparaison des images basée sur l'apparence de celui-ci. Une grande partie de l'information visuelle perçue dans les images d'extérieures est due a l'illumination en provenance du ciel. Ce facteur est représenté par les rayons du soleil réfléchis et réfractés dans l'atmosphère en créant une illumination globale de l'environnement. En même temps cet environnement détermine la façon de percevoir les objets du monde réel. Etant donné l'importance du ciel comme source d'illumination, nous formulons un processus générique en trois temps, segmentation, modélisation et comparaison des pixels du ciel, pour trouver des images similaires en se basant sur leurs apparences. Différentes méthodes sont adoptées dans les phases de modélisation et de comparaison. La performance des algorithmes est validée en trouvant des images similaires dans de grandes collections de photos. La deuxième partie de cette thèse consiste a exploiter l'information géométrique additionnelle pour en déduire les caractéristiques photométriques de la scène. A partir d'une structure 3D récupérée en utilisant des méthodes disponibles, nous analysons le processus de formation de l'image a partir de modèles simples, puis nous estimons les paramètres qui les régissent. Les collections de photos sont généralement capturées par différents appareils photos, d'où l'importance d'insister sur leur calibrage radiométrique. Notre formulation estime cet étalonnage pour tous les appareils photos en même temps, en utilisant une connaissance a priori sur l'espace des fonctions de réponse des caméras possibles. Nous proposons ensuite, un cadre d'estimation conjoint pour calculer une représentation de l'illumination globale dans chaque image, l'albedo de la surface qui compose la structure 3D et le calibrage radiométrique pour tous les appareils photos.

No description available.

La structure du milieu interstellaire est très hétérogène dans notre Galaxie. Les observations, ainsi que les lois d'échelle, suggèrent que sa structure pourrait être fractale. Les plus petites structures de gaz ou clumps auraient un rayon de quelques unités astronomiques et une température de quelques kelvins, les rendant très difficilement détectables. Ces nuages de gaz moléculaires froids et denses peuvent donc être des candidats peu lumineux à une matière noire baryonique.L'objectif de cette thèse est de trouver des indices observationnels de l'existence de ces sous-structures gazeuses. La méthodologie consiste à comparer les magnitudes apparentes d'étoiles à différents intervalles de temps et vers plusieurs directions dans la Galaxie; afin de détecter des variations photométriques temporelles compatibles avec un événement d'obscurcissement. En effet, si un nuage moléculaire se trouve en amont d'une étoile, celle-ci sera obscurcie ou occultée par les sous-structures composant le nuage moléculaire, pendant une certaine durée de temps (dépendant de la vitesse de chaque objet et du temps d'observation considéré). Dans un premier temps, nous avons effectué un traitement photométrique complet pour chercher de tels événements dans l'amas globulaire NGC 4833 en comparant des observations photométriques séparées de six mois. Cet amas fut observé par D. Puy et D. Pfenniger en janvier et juillet 2006. Puis, nous avons utilisé des observations de plusieurs amas globulaires observés par le télescope Hubble, dont la photométrie est connue et rendue publique.Finalement, nous avons développé des méthodes d'apprentissage automatique afin de détecter des événements d'obscurcissement dans le Sloan Digital Sky Survey, Stripe 82.Aucun événement d’obscurcissement d'étoiles n'a pu être détecté ni dans la direction des amas globulaires considérés ni dans la région de Stripe 82. Néanmoins, des simulations d'un nuage fractal ont été réalisées pour estimer la probabilité d’occurrence du phénomène d'obscurcissement d'étoiles. Elles montrent que l’événement est extrêmement rare, ce qui permet d'émettre une nouvelle stratégie d'observation pour augmenter les chances de détecter de tels événements. De plus, les méthodes d'apprentissage automatique utilisées pour cette étude sont appliquées à d'autres objets variables tel que les quasars, très nombreux dans Stripe 82, et peuvent être utilisées pour d'autres objets variables tels que les supernovae. / The scaling laws and the very appearance of the interstellar medium suggest that this medium might have a fractal structure. The smallest structure of gas called clumps might have a radius about few astronomical units and a very low temperature of few kelvins, making them very difficult to detect. That is why, cold, dense, molecular clouds could be ideal candidates for baryonic dark matter. Because of their low temperature, they do not radiate, and seem to be dark.The aim of this thesis is to highlight observational clues on the existence of clumps by comparing apparent magnitudes of stars at different times, and towards different directions in the Galaxy. If a molecular fractal cloud lies behind a star, it will be obscured or occulted within a certain time (depending on the velocity of each object and on observational time). Therefore the goal of this thesis is to detect magnitude-time variations of stars compatible with an obscuration event.First, we did a complete photometric data reduction of observations of a globular cluster, NGC 4833, carried out by D. Puy and D. Pfenniger in January and July 2006. We wanted to detect obscuration events of stars in NGC 4833, in a six-month period. Then, we used public data of several globular clusters observed by the Hubble Space telescope.Finally, we developed machine learning methods to detect obscuration events in the Sloan Digital Sky Survey, Stripe 82.No obscuration events were detected either in several globular clusters, or in the Stripe 82 survey. However, we performed simulations of a fractal cloud to estimate probabilities of obscuration events. They show that such events are very rare, allowing us to develop a new observational strategy for increasing the probability of detecting obscuration events. Moreover, machine learning methods used in this study are applied to other variable objects, such as quasars which are numerous inside Stripe 82 survey, and could be used to other variables objects such as supernovae for example.

Plusieurs études ont montré que les étoiles bleues de la branche horizontale (étoiles BHB) chaudes (Teff ≳ 11500 K) des amas globulaires présentent des anomalies d’abondances causées par le processus de la diffusion atomique. Il a été démontré que la diffusion atomique entraîne selon l’étoile, une stratification verticale de certains éléments chimiques. L’objectif de ce projet vise à élargir notre compréhension des étoiles BHB en menant une analyse de leurs propriétés physico-chimiques. L’emphase est mise sur les étoiles BHB du champ qui ont été beaucoup moins étudiées, mais qui, de par leur histoire évolutive et leur environnement différents, devraient permettre d’élargir le domaine des paramètres de base (comme la métallicité). Des modèles d’atmosphère stellaire calculés avec le code atmosphérique PHOENIX sont utilisés pour ce travail. Le spectrographe ESPaDOnS du télescope Canada-France-Hawaï a permis d’obtenir des spectres de très haute résolution (R = 81000) pour quatre étoiles BHB du champ. Le code de transfert radiatif ZEEMAN2 a été utilisé pour une analyse détaillée des raies spectrales. Les résultats de cette analyse indiquent l’existence d’une surabondance de métaux dans les étoiles chaudes HD213781 et Feige 86, marquant ainsi la présence de la diffusion atomique dans ces étoiles. Une stratification verticale du phosphore et une stratification marginale du fer sont détectées dans HD213781 et Feige 86, respectivement. Comme leurs homologues des amas globulaires, les étoiles BHB froides du champ HD128801, HD143459 et HZ27 ne montrent aucune signature de diffusion. Ces variations de l’abondance des éléments soulèvent une question lors de l’évaluation des paramètres atmosphériques (Teff et log g) à partir de modèles de métallicité fixe. De ce fait, l’effet de la métallicité est ici étudié pour 20 étoiles BHB connues dans les amas globulaires M3 et M13. Les paramètres atmosphériques de ces étoiles sont alors extraits en utilisant des modèles d’atmosphère de six métallicités [Fe/H] différentes variant entre -2:0 et +0:5 dex par rapport à la métallicité solaire. Deux méthodes d’analyse, spectroscopique et photométrique, sont utilisées. Seules les valeurs spectroscopiques de Teff et log g obtenues pour les étoiles plus chaudes montrent un changement, bien que léger, avec la métallicité. Par ailleurs, cette analyse souligne qu’un simple ajustement de la métallicité des modèles ne peut régler le problème de la faible gravité spectroscopique observée pour les étoiles BHB chaudes. Enfin, pour augmenter l’échantillon d’étoiles BHB chaudes connues dans le champ, de nouvelles observations de 21 étoiles ont été faites à l’Observatoire du Mont-Mégantic. Entre autres, trois nouvelles étoiles BHB du champ ont été caractérisées pour la première fois. / Several studies have shown that hot (Teff > 11500 K) blue horizontal branch stars (BHB stars) in globular clusters present abundance anomalies caused by the atomic diffusion process. It was demonstrated that atomic diffusion leads in certain stars to vertical stratification of some chemical elements. The aim of this project is to expand our comprehension of BHB stars by analyzing their physico-chemical properties. A special emphasis is put on field BHB stars. These stars are not as well studied, but because of their different evolutive history and environment, they should permit a widening of the domain of the basic parameters (like metallicity). Stellar atmosphere models computed with the atmospheric code PHOENIX were used for this work. The échelle spectrograph ESPaDOnS at the Canada-France-Hawaï telescope allowed to obtain high resolution (R = 81000) spectra of four field BHB stars. The radiative transfer code ZEEMAN2 was used for a detail analysis of the spectral lines. The results of this spectral analysis show the existence of an overabundance of metals in the hot stars HD213781 and Feige 86, indicating the presence of atomic diffusion in these stars. A vertical stratification of phosphorous and a marginal vertical stratification of iron are detected in HD213781 and Feige 86, respectively. Like their counterparts in globular clusters, the cold field BHB stars HD128801, HD143459, and HZ27 present no signature of atomic diffusion. These abundance variations raise a question for the evaluation of the atmospheric parameters from models with a fixed metallicity. Thereby, the metallicity effect is studied here for twenty known BHB stars in the globular clusters M3 and M13. The atmospheric parameters of these stars are extracted using atmosphere models for six different metallicities [Fe/H] in the range between -2:0 and +0:5 dex relative to the solar metallicity. Two analysis methods, based on spectroscopy and photometry, were used. Only the spectroscopic values of Teff and log g obtained for the hottest stars show a variation, although small, with the model metallicity. Otherwise, this analysis demonstrates that a simple adjustment of the model metallicity cannot solve the problem of the low spectroscopic gravity observed for the hot BHB stars. Finally, to increase the sample of known hot field BHB stars, new observations of 21 stars have been made at the Observatoire du Mont-Mégantic. Three new field BHB stars have been characterized for the first time.

Les procédés de combustion sont utilisés dans la plupart des industries chimiques, métallurgiques et manufacturières, pour produire de la vapeur (chaudières), pour sécher des solides ou les transformer dans des fours rotatifs (ou autres). Or, les combustibles fossiles qui les alimentent (ex. : gaz naturel) sont de plus en plus dispendieux, ce qui incite plusieurs compagnies à utiliser d’autres sources de combustibles tels que de la biomasse, des rejets inflammables produits par le procédé lui-même ou des combustibles fossiles de moindre qualité. Ces alternatives sont moins coûteuses, mais de composition, et donc de pouvoir calorifique, plus variable. De telles variations dans la chaleur dégagée par la combustion perturbent l’opération des procédés et la qualité des produits qui dépendent de ces installations. De nouvelles stratégies de contrôle de la combustion doivent donc être élaborées afin de tenir compte de cette nouvelle réalité. Il a été récemment démontré que l’énergie dégagée par la combustion est corrélée à l’aspect visuel de la flamme, principalement sa couleur, ce qui permet d’en quantifier les variations par imagerie numérique. L’objectif de ce projet industriel consiste à faire la démonstration que l’analyse d’images multivariées peut servir à l’identification du comportement d’une chaudière à biomasse. La chaudière à biomasse opérée par Irving Pulp &amp; Paper Ltd (Saint-John, Nouveau-Brunswick) fera office d’exemple. Les résultats montrent qu’un modèle bâtit à partir des informations fournies par les images ainsi que les données de procédé donne de bonnes prédictions de la quantité de vapeur produite (R2modèle=93.6%, R2validation=70.1%) et ce, 2,5 minutes à l’avance. Ce projet est la première étape du développement d’une nouvelle stratégie de contrôle automatique de la combustion de biomasse, capable de stabiliser l’énergie dégagée, malgré les variations imprévisibles dans le pouvoir calorifique et les proportions des combustibles utilisés provenant de différentes sources. / Biomass is increasingly used in the process industry, particularly in utility boilers, as a low cost source of renewable, carbon neutral energy. It is, however, a solid fuel with some degree of moisture which feed rate and heat of combustion is often highly variable and difficult to control. Indeed, the variable bark properties such as its carbon content or its moisture content have an influence on heat released. Moreover, the uncertain and unsteady bark flow rate increases the level of difficulty for predicting heat released. The traditional 3-element boiler control strategy normally used needs to be improved to make sure the resulting heat released remains as steady as possible, thus leading to a more widespread use biomass as a combustible. It has been shown in the past that the flame digital images can be used to estimate the heat released by combustion processes. Therefore, this work investigates the use of Multivariate Image Analysis (MIA) of biomass combustion images for early detection of combustion disturbances. Applied to a bark boiler operated by Irving Pulp &amp; Paper Ltd, it was shown to provide good predictions, 2.5 minutes in advance, of variations in steam flow rate (R2fit=93.6%, R2val=70.1%) when information extracted from images were combined with relevant process data. This project is the first step in the development of a new automatic control scheme for biomass boilers, which would have the ability to take proactive control actions before such disturbances in the manipulated variable (i.e. bark flow and bark properties) could affect steam production and steam header pressure.