Extraction Des Linéaments Géologique À Partir Des Images Landsat: Cas de La Plaine de Tamlilt (Haute Atlas Oriental)
Extraction Des Linéaments Géologique À Partir Des Images Landsat: Cas de La Plaine de Tamlilt (Haute Atlas Oriental)
Extraction Des Linéaments Géologique À Partir Des Images Landsat: Cas de La Plaine de Tamlilt (Haute Atlas Oriental)
Présenté par:
LAMZOURI Hala
Encadrée par:
1
Université Sidi Mohammed Ben Abdellah
Faculté des Sciences et Techniques
www.fst-usmba.ac.ma
-----------------------------------------------------------------------------------------------------------
LAMZOURI Hala
Titre : Extraction des linéaments géologique à partir des images Landsat : Cas de la
région de Tamlilt (Haute Atlas Oriental)
Résumé
Les éléments de base de cette étude sont l’analyse et l’interprétation d’une image satellitaire
Landsat en vue de la cartographie géologique.
L’apport essentiel de la télédétection est la vue synoptique. Dans cette analyse on a fait la
comparaison entre les différents documents (Carte et image traitée).
Lors des missions du terrain, l’objectif est de mieux cibler et repérer géographiquement les
structures identifiées. Le report des informations dans une base de données géographique,
permettent d’établir un schéma compréhensif dans lequel les éléments prennent leur
signification géologique.
Il existe de nos jours une panoplie de satellites qui permettent de distinguer plusieurs types
de minéraux pour une cartographie géologique très exacte et permettent de gagner plus de
temps tout en ciblant les campagnes de vérité terrain.
2
Liste des figures
Traitement utilisés
Figure 9 : 5x5……………………………………………………………………..………..…23
Figure 14 : carte des nouvelles fractures obtenues à partir des images satellite …….….……26
3
Dédicace
Ma chère mère et mon cher père
Dieu seul sait combien je vous aime et combien je vous chérie.
L’éducation que vous m’avez donnée est un trésor inépuisable.
Le chemin que je suis en train de tracer aujourd’hui est le couronnement de vos sacrifices, de votre
soutien incomparable, de votre patience et de vos efforts.
Grâce à vos conseils sages, vos encouragements pour moi, votre compréhension, vous m’avez
épargné beaucoup de souffrances et vous m’avez aidé à surpasser tous les obstacles qui se sont
dressés sur ma route.
Je vous aime, et j’espère que vous trouveriez dans ce travail l’expression de ma gratitude et le fruit
de vos sacrifices.
Puisse le tout puissant vous procurer santé et bonheur, puisse-t-il m’aider pour jamais vous
décevoir.
Mes chères sœurs & frères
Najlae, Hajar, Khaoula, , Ayman ,Abdelwahed , Otmane, Hiba et Sara
En témoignage de mon amour fraternel, je vous remercie du fond du cœur pour votre soutien moral
et vous encouragements permanents. Je vous souhaite une vie pleine de bonheur, réussite et joie.
Que Dieu vous protège et vous procure longue et heureuse vie.
A toute ma famille
Votre amour pour moi m’a toujours poussé à faire de mon mieux pour que vous soyez fiers de moi.
En gage de mon profond attachement, à l’aide précieuse que vous m’avez toujours portée.
Veuillez trouver ici l’expression de mon dévouement, ma reconnaissance et mon grand respect.
A tous mes très chers amis :
Pour l’amitié sincère et l’affection profonde que nous partageons, pour les moments de joies que
nous avons passées ensemble. Je vous dédie ce travail en vous souhaitant une vie pleine de santé et
de réussite.
Je vous remercie sincèrement pour tout ce que vous m’avez apporté et appris.
Que Dieu vous garde et vous comble de joie.
A tous mes collègues :
Que ce projet soit le témoignage de notre respect.
A toute la famille de Géo-ressource et environnement
Veuillez trouver ici l’expression de mon dévouement, ma reconnaissance et mon vif attachement.
4
Remerciement
Nos remerciements vont particulièrement :
Aux membres de jury (Mme Rais et Mr. Jabrane) qui ont bien voulu faire l’honneur
d’évaluer et d’apprécier ce travail.
A tous les enseignants du département environnement pour les efforts énormes qu’ils
ont fournit tout au long de notre formation.
Mme S. MNISSAR à tous ses conseils, ses efforts pour effectuer ce projet.
Enfin, il est pratiquement impossible d’énumérer tous ceux qui ont apporté leurs aides
et à qui nous exprimons également toute notre gratitude.
5
Plan :
Introduction
Chapitre I : présentation de l’entreprise
1) Présentation de l’ONHYM……………………. ………………………..…..…6
1-1) Création et organisation de l’ONHYM ……………………. ……………6
1.1.1) Historique ……………………. ………………..…..……….…7
1.1.2) Mission et domaine d’activité de l’ONYM…………………….7
Chapitre II : la géologie
1) Géologie de Maroc …………………………………………………………....8
1.1) domaine saharien et anti atlassique ……………………………………..8
1.2) domaine atlassique mésetien ……………………………………………8
1.3) domaine rifain …………………………………………………………..8
2) Géologie du Maroc Oriental ………………………… ……………………......9
3) Géologie local la bouteniére de Tamlelt …………… ………………………..12
3.1) cadre géologique ………………………………………………………..13
3.2) Géologie de la plaine de Tamelelt …………………………..…………14
3.3) Mines ……………………………………………………………………14
Chapitre III : Généralité sur la télédétection
1) Définition ……………………. ………………..…..…………………15
2) Principe ……………………. ………………..…..…………...………15
3) Intérêt ……………………. ………………..…..…………………..…15
4) résolution spatiale ……………………. ………………..…..………...16
5) résolution spectrale ……………………. ………………..…..………16
6) images satellite ……………………. ………………..…..…………...16
7) images Landsat ……………………. ………………..…..………...…17
Chapitre VI : Extraction des structures par l’utilisation d’ERDAS et ArcGis
1) Présentation d’ERDAS……………………. ………………..………..18
2) étapes de cette étude d’ERDAS ……………………. ……………….18
3) traitements utilisés……………………. ………………..…..……...…22
4) conclusion……………………. ………………..…..…………………24
Chapitre V : Analyse et des résultats de l’interprétation par Imagerie Landsat
1) Interprétation des Linéaments……………………. ……………….....26
2) Confrontation avec les cartes géologiques existantes………………...27
3) cartographie des altérations……………………. ………………..…...27
4) conclusion …………………. ………………..…..…….……………..28
6
Introduction Générale
Le stage dans une entreprise, représente une étape importante, puisqu’il permet de découvrir le
milieu professionnel et de compléter la formation théorique par une expérience pratique et de vivre
de proche le climat de la responsabilité.
Ce rapport présente l’application que j’ai effectuée au cours de mon stage au siège de l’ONHYM
(office national des hydrocarbures et des mines) à rabat .
Les objectifs de ce travail sont multiples : premièrement comprendre le travail de base et les
mécanismes de fonctionnement de l’entreprise, utiliser les capacités d’observation, d’analyse et
d’intégration dans une organisation, et intégrer les connaissances acquises tout au long de la
formation universitaire dans le cadre professionnel.
7
Chapitre I :
Présentation Générale de l’ONHYM :
1) Présentation d’ONHYM
L’ONHYM- L’Office National des Hydrocarbures et des Mines- est un office au sein du
Ministère de l’énergie et des mines à caractère semi-public et s’occupe de toutes études et
recherches visant à l’exploitation et à la valorisation des substances à l’exclusion des
hydrocarbures (pétrole).Il est soumis à la tutelle de l'Etat, laquelle a pour objet d'assurer le
respect, par les organes compétents de l'office, des dispositions de la loi 33-01.
8
Les Laboratoires Mines réalisent des prestations aussi bien pour les géologues de l’ONHYM
que pour les clients externes. Ces prestations concernent les trois domaines suivants :
Minéralurgie, Chimie Minérale et Minéralogie et pétrographie
1.1.1) Historique :
L’office Nationale des Hydrocarbure et des Mines (ONHYM) créé le 17 aout 2005,
par la fusion du bureau de Recherche et de Participation Minière (BRPM) et de l’Office
National de Recherche et d’Exploitations pétrolières (ONAREP).
De puis leur création, en 1928 pour le BRPM et en 1891 pour l’ONAREP, ces deux
organismes ont été les leaders au Maroc dans leurs domaines d’activité respectifs, et ont
bénéficié d’une grande renommée auprès des opérateurs internationaux .
C’est donc un nouvel élan, porté par une stratégie réactualisée dans le droit fil des
orientations gouvernementales, que l’ONHYM repend et développe les acquis des organismes
fusionnés.
9
Chapitre II :
La géologie
1) Géologie du Maroc
Situé à l’extrémité occidentale de l’Afrique du nord, le Maroc est bordé par la méditerranée au
nord et l’Atlantique à l’ouest. Sa configuration actuelle résulte des effets des différents cycles
orogéniques qui s’y sont succédés au cours des temps géologiques à savoir les cycles précambrien,
hercynien et alpin. Trois grands domaines structuraux y sont reconnus ce sont du Sud au Nord :
10
PRINCIPAUX DOMAINES
STRUCTURAUX
DU MAROC Pays des
Horsts
Hauts
0 200Km Massif Hercynien Plateaux
Central
Haut Atlas
Haut Atlas oriental
Jebilet central
Haut Atlas
occidental Antis Atlas
occidental
Antis Atlas
central
Bassin de Tindouf
Quaternaire et
Cénozoïque
Mésozoïque
Paléozoïque
Précambrien
Le Haut Atlas Oriental s’étend depuis l’est de la ville de Bou Arfa jusqu’à la frontière maroco-
algérienne (figure 1). Les formations géologiques correspondent, à l’ouest à la plaine de Tamlelt et
12
Figure 2: Localisation de la zone d’étude sur la carte du Maroc
La morphologie en dépression de cette zone est causée par une subsidence active et prolongée qui
poursuit actuellement à la faveur de failles quaternaires longitudinales au N et au S de cette plaine.
Au cours des années (1997-1998), des travaux de recherches géologiques à des échelles détaillées
ont concerné surtout les minéralisations aurifères. Celles-ci sont encaissées dans des terrains d’âge
indéterminé et situées dans la bordure Nord de la boutonnière
13
3.2) Géologie de la plaine de Tamlelt
La plaine de Tamlelt est formée essentiellement par une surface quaternaire. Les
affleurements y sont rares et apparaissent surtout prés des bordures où on note les formations
suivantes :
Un petit affleurement de schiste satiné métamorphique du PII juste au Sud du Douar de Ain
Chair.
Des formations volcaniques du PIII près de Menhoulou et de Guelb bou salem
Des affleurements de schiste et quartzite du cambrien qui apparaissent surtout dans la bordure
Ouest de la plaine.
Les schistes épimétaporphiques attribués au Primaire indéterminé qui apparaissent dans la
bordure Nord
3.3) Mine
L’activité minière occupe la seconde place dans l’économie de la province, et favorisée par
l’importance des substances minières du plomb, de Zinc et de la Barytine. C’est une activité
artisanale, lancée en sens du Dahir N° 1-60-019 du 1er Décembre 1960 par la Centrale D’achat et de
Développement de la région Minière de Tafilat et de Figuig (CADETAF) qui encadrent les artisans
sur les plans administratif et technique et assure la commercialisation de la production.
L’exploitation minière se fait à titre individuel ou dans le cadre de coopératives qui sont aux
nombres de 15. Egalement, il y a à signaler l’exploitation d’un des plus grands gisements de
Barytine au Maroc (Mine de Zelmou à la commune rurale de Bouanane) par la Compagnie
Marocaine de Barytine (COMABAR). Il Existe d’autres gisements renfermant les minerais de
Manganèse et de cuivre.(Mnissar, 2001)
14
Chapitre III :
Généralité sur la télédétection
1) Définition :
La télédétection est l’ensemble des techniques qui permettent, par l’acquisition d’images, d’obtenir
de l’information sur la surface de la Terre, sans contact direct avec celle-ci. La télédétection englobe
tout le processus qui consiste à capter et enregistrer l’énergie d’un rayonnement électromagnétique
émis ou réfléchi, à traiter et analyser l’information qu’il représente, pour ensuite mettre en
application cette information.
2) Principe :
3) Intérêt :
La télédétection géologique concerne de nombreux champs d'application parmi lesquels vient en
premier ordre la cartographie géologique, suivie de la recherche de ressources du sous-sol et le suivi
de certains phénomènes naturels. On utilise les propriétés spectrales des objets pour les reconnaître:
c'est ce qui a donné lieu à une nouvelle discipline de la télédétection : la géologie spectrale.
Dans cette section on a fait l'analyse de l'imagerie satellite est réalisée afin de localiser les structures
géologiques qui peuvent contenir des hydrocarbures. Durant cette analyse de la surface, le photo-
interprète peut estimer l'étendue et la géométrie ainsi que l'origine des structures, ce qui donne une
esquisse primaire de ce qui se passe en profondeur. Le choix de la résolution est important si on
compte avoir plus de détails sur les structures géologiques de surface. (Mnissar, 2001).
4) résolution spatiale :
La plupart des images de télédétection sont composées d’une matrice d’éléments d’image appelés
« pixels », soit les plus petits éléments d’une image. Le pixel est habituellement carré et représente
une partie de l’image. Il importe, cependant, de faire la distinction entre la taille des pixels et la
15
résolution spatiale, les deux notions n’étant pas interchangeables. Lorsque l’image provenant d’un
capteur ayant une résolution spatiale de 30 m est affichée à l’écran à pleine résolution, chaque pixel
correspond à une superficie au sol de 30 m x 30 m. Dans notre image, la taille des pixels et la
résolution sont identiques. Il est toutefois possible d’afficher une image sur laquelle la taille des
pixels sera différente de la résolution spatiale. Sur de nombreuses affiches montrant des images de la
Terre provenant d’un satellite, les pixels ont été réduits proportionnellement pour montrer de plus
vastes superficies, mais cela ne change en rien la résolution spatiale du capteur qui a pris l’image.(
nilerak.hatfieldgroup.com)
5) Résolution spectrale :
On peut distinguer différente classes de caractéristiques et de détails dans une image en comparant
leurs réponses sur un ensemble précis de longueurs d’onde. On distingue habituellement les grandes
classes, comme l’eau et la végétation, en utilisant un grand intervalle de longueurs d’onde, par
exemple le visible ou le proche infrarouge .d’autres classes plus spécifiques, notamment les divers
types de roche, ne sont pas toujours aussi faciles à différencier par l’un ou l’autre de ces intervalles
de longueurs d’onde. Elles exigent l’utilisation d’un intervalle de longueurs d’onde beaucoup plus fin
et, par conséquent, un capteur ayant une résolution spectrale beaucoup plus élevée. La résolution
spectrale renvoie à la capacité d’un instrument à capter de petits intervalles de longueurs d’onde. Plus
la résolution spectrale est élevée, plus l’intervalle de longueurs d’onde est fin pour une bande ou un
canal donné. (nilerak.hatfieldgroup.com)
6) images satellite :
Une image satellite est comme une photo d'une partie de la Terre, prise de l'espace. Cette image est
prise par un satellite placé en orbite autour de la planète. A la différence des images obtenues d'un
appareil photo, ou dessinées sur du papier, une image satellite est une image numérique, traité
informatiquement, et celle-ci est mise en place à partir d'ondes transmises par le satellite à
l'équipement au sol. L'image numérique est fabriquée par un assemble de pixels, représentés par une
série de bits, ayant chacun comme caractéristique un niveau de couleur ou de gris prélevé à
l'emplacement correspondant dans l'image réelle, cet emplacement étant défini par ses coordonnées
cartésiennes dans un plan d'ordonné y, d’abscisse x et de centre O.
Grâce à leur altitude, les satellites fournissent des images de grande taille. Dans le jargon, l'image
prise est appelée « scène ». Effectivement, une image satellite représente en moyenne une « scène »
16
de 60km de côté, mais tout dépend du satellite, cela peut varier. La résolution des images joue un rôle
important : la précision d'un pixel peut aller de 20 à 30 mètres en couleur et jusqu'à 1 ou 2 mètres en
noir et blanc. Ceci dépend des capacités et de l'altitude du satellite. Les résolutions les plus hautes
sont malheureusement encore réservées aux services militaires. Grâce à la technologie qui s'améliore
de jour en jour, les satellites les plus récents fournissent des images comparables à des photos
aériennes.
7) Image Landsat
Le programme d’observation de la surface terrestre Landsat est le plus ancien des Etats Unis. Sept
satellites Landsat ont été lancés de puis juillet 1972. Le dernier d’entre eux, Landsat 7.
Les satellites Landsat 4,5 et 7 évoluent à une altitude moyenne de 705 km, sur des orbites circulaires
quasi polaire caractérisées par une inclinaison de 98,2° (ce qui permet l'héliosynchronisme). Un tour
de la Terre leur prend 98,9 min, si bien qu'ils décrivent 14,5 révolutions par jour. Un cycle orbital
complet dure 16 jours.
Ces caractéristiques orbitales ont des conséquences sur l'acquisition des images :
- la couverture est complète entre les parallèles 81° nord et sud.
- le recouvrement des images s'accroît avec la latitude (7% à l’équateur, 54% par 60° de latitude)
La résolution spatiale est passée de 80 m à 30 m (ou 15 m en mode panchromatique) et les
domaines spectraux explorés concernent le visible, l'infrarouge proche et moyen ainsi que
l'infrarouge thermique. (www.futura-sciences.com)
17
Chapitre VI :
Extraction des structures par l’utilisation des logiciels
ERDAS et ArcGIS
1) Présentation d’ERDAS
ERDAS IMAGINE vous pouvez traiter des images et en extraire des informations tel un
professionnel chevronné, indépendamment de votre expérience. Vous pouvez analyser les
données de n'importe quelle source. ERDAS IMAGINE simplifie la classification,
l'orthorectification, le mosaïquage, la projection et l'interprétation des images, en maintenant
l'intégrité des données géospatiales pour la mise à jour de votre SIG.
ERDAS Imagine combine la télédétection et les possibilités SIG, permettant ainsi de créer
des données géospatiales, d'extraire le maximum d'informations et de mettre à jour vos données
SIG. (www.quelsoft.com)
2) étapes de cette étude d’ERDAS
18
Figure 4 : sélection de couches téléchargées
b) la mosaïque d’images
Le collage des images de chaque (row-path)
dataprep image mosaic mosaic tool
19
D) la superposition de la zone d’étude sur les images Landsat
20
E) le coupage de la zone de travail par “ AOI“
Les traitements spécifiques sont utilisés pour obtenir une différenciation des fractures. Ces
traitements sont rendus possibles grâce au format numérique des images satellitaires qui offre de
multiples possibilités d’analyse.
Les images numériques sont riches en informations de diverses natures sur des objets au sol et ne
peuvent de ce fait être directement exploitées. Pour faciliter les opérations d’interprétation et mieux
discerner les éléments d’analyse structurale, des traitements spécifiques (filtres directionnels) à la
cartographie structurale ont été utilisés
21
Les traitements utilisés :
1x5 5x5
22
Figure 12: Texture Figure13: shaded
On a utilisé quelques traitements pour mettre en évidence les structures qui existent sur
l’image satellitaire.
Pour l’application d’un filtre 11x11 on’ a pas la possibilité de montrer les structures
présentent sur l’image, (un peu flow).
Le filtre 1x5 est un bon traitement de voir même à une petite échelle sur l’image.
Les images de Texture et Shaded n’ont pas donné de bons résultats.
Le traitement qui montre bien les structures sur l’image est l’application du filtre 7x7
Conclusion
Les traitements d’images satellisables font ressortir des systèmes structuraux et des systèmes
hydrologiques qui ne pouvaient pas être observés aussi complètement par les anciennes techniques.
Les applications économiques de ce genre d’études sont importantes particulièrement pour les
recherches minérales. Cette approche permet de renforcer les connaissances théoriques, en
particulier dans le domaine structural grâce à l’application de filtres directionnels.
Ces données peuvent être complétées par intégration des données des campagnes géophysiques et
des relevés géochimiques dans les systèmes d’information géographique.
Les capteurs utilisés pour cette étude le TM (Thematic Mapper) de Landsat 5 de 30m de résolution.
23
Chapitre V :
Analyse des résultats de l’interprétation d’image
Landsat TM
Les linéaments sont déduits suivant des critères bien précis : il s’agit de discontinuités
correspondant soit à des fractures montrant un décalage entre les formations géologiques ou des
prolongements des failles représentées sur la carte existante, soit à un contact anormal entre
deux formations géologiques différentes sans décalage. Les linéaments qui ne présentent pas de
décalage entre les structures ont fait l’objet d’une vérification sur le terrain. La figure 14
représente la carte structurale regroupant tout le réseau de fractures nouvelles, cartographiées de
la région.
Figure 14: Carte des nouvelles fractures obtenues à partir des images satellite (MNISSAR.S
2002)
La carte structurale élaborée à partir de l’image satellitaire montre que la fracturation s’organise selon
quatre directions principales : E-O à ENE-OSO, NNE-SSO, NNO-SSE et N-S.
24
- La direction E-O à ENE-OSO : elle correspond à des failles plurikilométriques (figure 15) qui
affectent aussi bien les terrains triasiques jurassiques que crétacés. Parmi les plus importantes,
citons celles qui affleurent parallèlement à l’Oued Tisserfine (F1, F2 et F3 de la figure 15).
Sur le terrain, cette direction E-O correspond soit à des failles à jeu normal au cours du Trias
et du Lias inférieur (Piqué et al., 1998) soit à des failles réactivées avec un jeu inverse.
- La direction NNE-SSO correspond à des failles métriques à kilométriques qui affectent les
terrains du Jurassique supérieur et Crétacé. Parmi les plus importantes citons celles de Jbel
Maïz (failles F4 et F5 de la figure 14). Sur le terrain, cette direction correspond à un premier
jeu normal, auquel succède un second jeu en décrochement senestre attesté par le décalage des
failles normales de direction E-O (figure 15).
- La direction NNO-SSE correspond à des failles décamétriques qui affectent les terrains du
Trias et du Jurassique. Sur la figure 15, la faille F6 partant de Jbel Maïz décale en
décrochement dextre les failles de direction E-O.
Une vision nouvelle des réseaux de failles qui ont affecté le Haut Atlas au cours du Mésozoïque et
du Cénozoïque. Les principales directions de failles mises en évidence E-O à ENE-OSO, NNE-
SSO, NNO-SSE et N-S ont guidé la structuration de cette région dont l’évolution géodynamique est
caractérisée par une phase extensive au Mésozoïque et une phase compressive Cénozoïque. Ces
25
résultats sont conformes à celles obtenues par plusieurs auteurs (Du Dresnay, 1963 ; Aït Brahim,
1991 ; EL Kochri et Chorowicz, 1996).
5. conclusion
Ces résultats sont conformes aux orientations cartographiées par du Dresnay (1963) mais
avec un nombre nettement plus important que la carte existante. Une analyse approfondie de ces
résultats donne une nouvelle vision de l’interprétation de la fracturation dans le Haut Atlas Oriental.
Une analyse approfondie de ces résultats donne une nouvelle vision de l’interprétation de la
fracturation dans le Haut Atlas Oriental. En effet, la formation de cette partie de la chaîne serait le
résultat d’un jeu de failles normales au Mésozoïque et de leur rejeu au Cénozoïque en failles
inverses couplées à des jeux en décrochement dextre.
Conclusion générale :
Cette approche a permis de réaliser une carte de la géométrie des structures qui
complète les travaux anciens réalisés dans la zone.
Les différents traitements réalisés lors de cette étude ont permis d’extraire de
manière directe et automatique les différents linéaments présents. En effet, nous
avons obtenu une cartographie numérique précise plus détaillée et qui peut être à
26
chaque fois corrigée de manière plus souple et facile à manipuler. Ceci grâce au
logiciel SIG utilisé dans ce travail.
27
Bibliographie :
Aït Brahim L. et Chotin P. (2001) - Paleostress evolution of the moroccan African margin from
triassic to present. Per-Tethys, Memoir 6: Mém. Mus. nat. Hist. Nat, Paris (à paraître).
Dresnay (1966) -
Houari (2003) -
http://www.quelsoft.com/fiche/erdas-imagine-c32-612-2022.html
http://www.futura-sciences.com/magazines/espace/infos/dico/d/univers-landsat-3738/
http://nilerak.hatfieldgroup.com/French/NRAK/EO/html/rsbch7.html
http://nilerak.hatfieldgroup.com/French/NRAK/EO/html/rsbch6.html
http://www.onhym.com/emploi-et-carriere.html
28