3.1. Présentation de la spécificité de MicMac pour le traitement des élévations

1MicMac est un logiciel très performant, qui nous permet de proposer ici une alternative à Metashape, qui est le logiciel actuellement le plus utilisé pour la production d’orthophotographies en archéologie. Seulement, ce dernier logiciel présente un inconvénient majeur à son utilisation dans des contexte scientifique, celui de produire des résultats difficilement voire même nullement reproductibles (particulièrement en cas d’utilisation du GPU) (Pamart et al. 2019). De même, pour produire des orthophotographies verticales, l’utilisateur doit faire lui-même un grand nombre de manipulations directement sur le nuage de points ou en post-traitement pour obtenir un résultat (masques, orientation manuelle du modèle, mise à l’échelle manuelle du résultat, etc.). La production d’orthophotographies passe nécessairement par des étapes de maillage et de texturage du modèle 3D, qui impliquent des interpolations induisant de potentielles altérations et inexactitudes dans les données produites. Ces aspects sont des limites à l’évaluation de la fiabilité et la reproductibilité de ses résultats, qui dépendent autant de la qualité du set de données que de la rigueur des utilisateurs. MicMac, malgré une interface en lignes de commande peu engageante, permet d’obtenir des relevés fiables et reproductibles : les lacunes de l’acquisition ne sont pas compensées de façon opaque par le logiciel, et il n’a besoin d’aucune intervention de l’utilisateur sur le modèle.

2Les orthophotographies verticales sont produites selon un calcul mathématique de géométrie dans l’espace, que nous verrons plus en détail dans la suite du chapitre. Dans les grandes lignes, cela consiste en le calcul d’un nouveau repère spatial local et cartésien temporaire, sur la base du géoréférencement : l’utilisateur renseigne le plan (c’est-à-dire le mur, la coupe, la face, etc.) dont il souhaite produire l’orthophotographie, qui permettra de basculer du géoréférencement, avec une orthophotographie calculée à partir du plan horizontal du modèle, à une orthophotographie verticale, non géoréférencée.

3.2. Cadre d’expérimentation

3Le protocole de photogrammétrie présenté ici est issu de nos travaux de recherche sur la restitution architecturale de l’église paléochrétienne de Mirine en Croatie (fig. 29). Elle fait partie du site archéologie de Mirine-Fulfinum, situé près du village d’Omišalj sur l’île de Krk. L’église de Mirine (fig. 30) a été construite en périphérie de la cité antique de Fulfinum, fondée ex nihilo à partir du i^er siècle de n. è. Le plan de l’église est en croix latine, à nef unique et chevet plat, flanqué, au nord et au sud, de deux modestes annexes liturgiques. La nef est précédée d’un vestibule flanqué d’une tour au sud. De proportion carrée, les bras de transept sont peu saillants. Construite en opus incertum de moellons de calcaire local, elle mesure environ 30 m de long par 12 m de large – hors vestibule et bras de transept – et est préservée en élévation sur 5 à 9 m de haut. L’atrium, ajouté dans une seconde phase contre la tour et l’église, moins bien conservé que cette dernière, mesure environ 25 m de côté.

Fig. 29 : photographie aérienne du complexe paléochrétien de Mirine (cl. M. Thivet, 2016).

Fig. 30 : situation du site de Mirine-Fulfinum (extraits de Čaušević-Bully et al. 2017 et Čaušević-Bully, Bully 2019).

4Partiellement sondée après la Seconde Guerre mondiale, la première étude architecturale de l’église est entreprise dans les années 1970 par Ante Šonje (Šonje 1978). Il réalise une étude synthétique formelle et stylistique, permettant de situer la construction de l’église à la moitié du v^e siècle. L’église n’est fouillée qu’à partir de 1993, jusqu’en 2005, sous la direction de Nino Novak (Bully, Čaušević-Bully 2011). Malheureusement, ces fouilles n’ont pas fait l’objet de publications régulières, et il n’existe pas de rapport de fouilles. Seules deux publications nous renseignent sur les problématiques de recherche de Novak (Novak, Bložić 1995 ; Novak 1999). Il s’est d’abord concentré sur les questions liées à la fondation de la cité de Fulfinum et ses relations avec Curicum (Krk) avant son abandon et la construction du complexe de Mirine. La première publication mentionne 6 sondages, ouverts dans la nef et dans la tour. La seconde se concentre uniquement sur la restitution du chœur de l’église, et ne s’attarde pas sur les fouilles qu’il dirige pourtant encore à ce moment-là.

5De nos jours, l’église ne fait plus l’objet de fouilles systématiques. Néanmoins, les secteurs à proximité, et particulièrement l’église aux trois absides, sont fouillés par Morana Čaušević-Bully dans le cadre du programme KVARNER de l’EFR depuis 2016. Les problématiques soulevées par les études des constructions autour de l’église ont conduit l’équipe à se réintéresser à Mirine : outre les prospections géophysiques qui ont été menées dans la nef et l’atrium en 2012, une acquisition de données aéroportées en vue d’un relevé photogrammétrique complet a été réalisée en 2016 par Matthieu Thivet. Cette campagne de relevés en photogrammétrie avait pour but de dresser un plan des vestiges et d’obtenir un modèle 3D. Étant donné les objectifs de ces premiers relevés, ils étaient inexploitables pour notre étude puisque les photos ont été prises en plongée et trop loin des vestiges. Cette même année, un dernier sondage a été ouvert dans l’atrium, sur la base des prospections géophysiques, mais il s’est avéré être négatif.

6Le travail de master dont est issu ce protocole se proposait de réactualiser les connaissances acquises sur l’église, malgré une bibliographie abondante la concernant (Gnirs 1914 ; Szabo 1930 ; Šonje 1978 ; Novak, Bložić 1995 ; Chevalier 1995a ; Chevalier 1995b ; Novak 1999 ; Bully, Čaušević-Bully 2011 ; Čaušević-Bully 2013), et notamment des travaux proposant des restitutions. Certains points de l’architecture n’ont jamais été véritablement abordés : les restitutions en volume existantes ne se sont concentrées que sur l’aspect des façades extérieures, ou sur le chœur. La question de l’évolution du bâtiment dans son architecture ou ses aménagements ne s’est jamais véritablement posée. Toute la difficulté de ce travail réside dans le fait que l’église a été restaurée et partiellement reconstruite peu de temps après les fouilles, d’autant plus que celles-ci n’ont jamais été publiées et qu’il n’existe aucun véritable relevé de l’état du bâtiment avant restauration. L’une des premières étapes de ce travail consistait donc à réaliser une critique d’authenticité des maçonneries, qui n’aurait pas été possible sans relevés des élévations.

7Compte tenu de ses dimensions conséquentes, nous nous sommes naturellement tournés vers des méthodes de relevés numériques, dont la photogrammétrie.

3.3. Spécificité de l’acquisition pour les élévations

8Les précautions à prendre concernant le choix du matériel, les grandes lignes du protocole de préparation et d’acquisition des photos et de la topographie varient peu de ce qui est déjà présenté dans le chapitre précédent. Nous présenterons ici nos propres conditions d’acquisition.

9Étant donné nos problématiques de recherches, nous avions prévu une résolution d’environ 1 cm/pixel. Sachant que les murs sont plutôt homogènes et exempts de décors sophistiqués, il n’était pas nécessaire d’être plus précis, même si dans les faits nous avons atteint une résolution de 1 mm/pixel. L’acquisition des clichés s’est faite, pour la plupart des élévations, avec un drone Dji Phantom 4 Pro. D’autres endroits, comme la façade occidentale, ont été acquis avec un appareil photo numérique reflex Canon EOS 6D équipé d’un objectif 24 mm à focale fixe et une perche. Ce choix s’explique par l’accessibilité tant au sol qu’aérienne, notamment en raison de la végétation abondante qui entoure l’église, du faible recul par rapport aux parois par endroits, et à cause de la grande hauteur des élévations.

10Cependant, tout l’édifice n’a pas été relevé en photogrammétrie, pour trois raisons : quelques élévations étaient très difficilement accessibles, quel que soit le matériel prévu ; celles-ci n’étaient pas indispensables à notre étude ; et l’ensemble du bâtiment a aussi été relevé au scanner laser terrestre. L’utilisation des deux méthodes était particulièrement pertinente dans notre cas, car nous nous intéressions aux maçonneries pour en évaluer l’authenticité, et il nous fallait disposer d’un relevé entier qui nous permette d’appréhender les dimensions et le volume de l’édifice, comme support pour la restitution.

Conditions d’acquisition

11Nous avons mené l’acquisition de chaque élévation de manière à les traiter indépendamment les unes des autres. Nous le verrons plus loin, il est possible d’extraire plusieurs orthophotomosaïques d’un même nuage de points, mais cela demande une grande puissance de calcul et capacité de stockage, et une préparation plus fastidieuse des données. Cela présente également l’avantage de pouvoir tester le set fraîchement acquis en cas de doute, les ordinateurs portables ayant une capacité de calcul et de stockage moindre qu’une machine fixe. Dans notre cas, nos sets ne dépassaient pas la cinquantaine de photos, sauf exception. Nous avons ainsi pu vérifier si nos traitements étaient satisfaisants, tant que nous étions sur le site, quitte à les compléter.

12Si les dimensions de l’élévation à relever vous le permettent, l’installation des GCP ne diffère pas vraiment du protocole d’orthophotographies en plan. À Mirine, nous n’avions pas connaissance des recommandations pour le placement des GCP. Nous en avons utilisé 4, que nous avons essayé de répartir sur toute la surface qui était accessible à pied et avec une échelle. Les auteurs spécialisés recommandent au moins 5 GCP, et jusqu’à 10 maximum s’il y a des paliers topographiques à couvrir (Girod 2018, p. 31-32). Il est suggéré d’en placer 4 aux angles et 1 au centre de l’objet. Si les élévations sont trop hautes, placez les 2 du dessus le plus haut possible.

13Comme pour une orthophotographie plane, privilégiez un cheminement en S, en déplaçant le capteur latéralement (droite, gauche, haut ou bas), tout en veillant à prendre les clichés à un axe le plus perpendiculaire possible par rapport à la paroi et ainsi obtenir un angle d’incidence faible (fig. 31). Il convient alors d’avoir recours à l’utilisation d’un drone ou d’une perche pour les parties hautes inaccessibles. Lorsque l’angle d’incidence de la caméra contre l’objet relevé est trop important, cela crée des anomalies sur le traitement qui gênent la lecture des résultats.

Fig. 31 : schéma en plan d’une acquisition d'élévation. Haut : déplacement du capteur par rapport à un point de rotation, à éviter. Bas : cheminement latéral entre les clichés, à privilégier.

Contraintes environnementales

14Il est essentiel d’éviter la présence de végétation sur les clichés : les mouvements dus au vent empêchent la corrélation entre les images de se faire correctement et les feuilles entraînent des anomalies lors du traitement. Si les zones végétalisées sont inévitables, contournez-les, quitte à prendre plus de photos. Par exemple, à Mirine, la paroi nord du mur gouttereau nord a demandé la prise de 150 clichés à cause de la végétation puisque nous n’avions pas suffisamment de recul pour n’en prendre qu’une cinquantaine, comme pour les autres élévations. Aussi, nous n’avons pas pu traiter la paroi extérieure de la tour, car les branches d’un arbre juste devant gênaient la prise de photos : soit parce que des feuilles étaient dans le champ, soit parce que nous ne pouvions pas manipuler l’appareil correctement. D’une manière générale, tout élément sensible au vent ou susceptible de bouger doit être soigneusement évité, fixé ou extrait de la scène, de même que les personnes, les véhicules en mouvement, les animaux, etc. (Van Dongen 2011)

15La paroi nord du mur gouttereau nord a été particulièrement retorse. Outre la végétation, le ciel était dégagé et l’église complètement ensoleillée le jour de l’acquisition, donc nous avons attendu la fin de journée pour éviter les ombres portées. Mais la paroi étant orientée nord-nord-est, certaines photos présentent un contre-jour inévitable. Le contre-jour est un éclairage plutôt défavorable pour le traitement et doit être proscrit si les clichés présentent des phénomènes de reflet et de dispersion de la lumière dans les lentilles de l’objectif, qui prennent la forme d’artefacts lumineux sur l’image (tels que des halos). Nous n’avons donc pas utilisé les photos trop touchées, quitte à avoir des lacunes dans le traitement à cause d’un moins bon recouvrement.

3.4. Traitement

16Le protocole présenté ici reprend en grande partie celui présenté plus haut, par conséquent, nous éviterons les redites pour la majeure partie des commandes et y renvoyons le lecteur. Nous présenterons le protocole à partir du premier changement dans la chaîne de commandes. Les commandes sont communes jusqu’à l’étape «  2.4.12. Visualisation du nuage de points homologues géoréférencés  » avec l’outil AperiCloud. Le reste du processus de traitement est directement issu d’un tutoriel du wiki de MicMac, dédié aux élévations, consultable à l’adresse suivante : https://micmac.ensg.eu/index.php/Pierrerue_tutorial. S’il est illustré avec le calcul d’orthoimages d’élévations d’un édifice, ce protocole est aussi adapté aux coupes stratigraphiques, avec une particularité que nous expliciterons dans le sous-chapitre suivant (cf. 3.5).

3.4.1. Troncage du nuage de points avant densification du nuage de points orienté

Commande

mm3d SaisieMasqQT sparse_cloud_geo.ply

Description

17SaisieMasqQT – Cet outil est une interface graphique permettant de créer un masque. Ce masque peut s’appliquer à un nuage de points ou à une photo.

18sparse_cloud_geo.ply – C’est le nom du fichier sur lequel on veut appliquer le masque. Il n’est pas obligatoire d’écrire le nom du fichier dans la commande, puisqu’un fichier peut être directement glissé-déposé dans la fenêtre de l’outil. Cependant, par souci de reproductibilité des résultats, nous recommandons qu’il soit inscrit dans le process.

À prendre en compte

19À cette étape, un masque est créé pour alléger le calcul du nuage dense. Ce n’est pas nécessaire pour le calcul d’une orthophotomosaïque en plan puisque les débords de l’objet relevé sont peu gênants. En revanche, lorsque l’on relève des élévations ou une coupe, il faut concentrer le calcul exclusivement sur la surface d’intérêt et éviter de calculer les éléments des autres plans. Le traitement est plus léger et le résultat plus lisible.

20Pour créer un masque, activez le mode édition avec F9 et dessinez un polygone, puis appliquez une action. Tous les raccourcis et actions sont consultables en tapant sur F1. Attention, le masque correspond à la surface à traiter, et non la surface à ignorer. Si l’on n’est pas attentif à cela, cela peut donner des résultats totalement aberrants ensuite. Le polygone apparaît en vert sur les images. Sur un nuage de points, les points ne faisant pas partie du masque deviennent gris/violets après validation (fig. 32). Quelques paramètres, dont la langue (Anglais, Français ou Espagnol), sont modifiables et sont conservés pour les prochaines utilisations.

Fig. 32 : interface graphique de saisie de masque 3D. Haut : sélection des points. Bas : validation des points à ne pas prendre en compte (points et caméras grisés). Mur gouttereau Nord de l’église paléochrétienne de Mirine, Omišalj, Croatie, 2021.

Production

21Cette commande produit deux fichiers XML : sparse_cloud_geo_polyg3d.xml et sparse_cloud_geo_selectionInfo.xml les informations 3D et 2D du polygone de sélection.

3.4.2. Densification du nuage de points

Commande

mm3d C3DC MicMac ".*JPG" CAL-Terrain Masq3D=sparse_cloud_geo.ply Out=dense_cloud.ply

Description

22C3DC – Outil de densification automatique du nuage de points à partir des images orientées.

23MicMac – Cet argument définit le mode de reconstruction du nuage de points. Celui-ci est un mode multistéréoscopique à moyenne densité. Pour chaque image, les calculs de corrélation sont effectués sur tous les couples possibles avec une définition divisée par 4 (ZoomF=4). Pour une densité élevée, il faut utiliser l’argument BigMac (ZoomF=2) ; pour une densité faible, le mode QuickMac (ZoomF=8).

24".*JPG" – Toutes les images sont prises en compte dans le calcul.

25CAL-Terrain – La commande se base sur le répertoire de calibration et d’orientations définitif.

Production

26La commande produit le fichier dense_cloud.ply, il s’agit du nuage de points dense au format PLY, visualisable avec CloudCompare. MMByPairFile.xml et MMByPairCAWSI.xml apparaissent également dans le répertoire de travail, ils enregistrent les couples stéréoscopiques utilisés pour la reconstruction. Deux nouveaux répertoires sont créés : PIMs-MicMac et Pyram. En effet, C3DC calcule les cartes de profondeur pour chaque image grâce à l’outil PIMs, qu’il stocke dans le premier. Le second conserve des copies de chaque photo sous-échantillonnée, quel que soit le mode de C3DC choisi.

3.4.3. Saisie du masque définissant le plan de l’orthophotographie

Commande

mm3d SaisieMasqQT 1.JPG Attr=plan

Description

27SaisieMasqQT – Même interface graphique qu’appelée plus haut, utilisée ici pour réaliser un masque 2D (fig. 33). N’importe quelle photo du set peut faire l’affaire, tant que la surface qui servira à calculer le nouveau plan y figure. Par analogie avec des projections géométrales d’architecture, ce plan est celui en face duquel la « vue » finale doit être. Il n’est pas nécessaire de dessiner une surface recouvrant tout le plan, un polygone quelconque suffit pour que MicMac en détermine le plan moyen.

Fig. 33 : interface graphique de saisie de masque 2D. Haut : sélection du polygone. Bas : validation du masque à appliquer (surface verte). Mur gouttereau Nord de l’église paléochrétienne de Mirine, Mirine-Fulfinum, Croatie, 2021.

281.JPG – Une seule photo est à appeler, la plus pertinente pour renseigner le plan de référence.

29Attr=plan – Cet argument définit le suffixe à appliquer à l’image correspondant au masque. Par défaut, la commande ajoute simplement _Masq à tous les masques, le suffixe permet alors de les différencier si plusieurs masques sont produits.

À prendre en compte

30Pour calculer des orthophotographies d’élévations avec MicMac, l’axe vertical du repère doit être soit face à la surface, tel que nous le serions nous-mêmes en se tenant face à un mur (même s’il présente un fruit ou un contre-fruit), soit perpendiculaire au plan dont on souhaite obtenir l’orthophotographie. Cette étape permet donc de définir le plan par rapport auquel la « vue » va se placer, à l’aide d’une image. Lors d’un processus d’orthophotographie en plan, l’axe vertical est déjà défini par le géoréférencement et n’a pas besoin d’être modifié, donc aucune intervention n’est nécessaire.

3.4.4. Calcul du repère temporaire

31C’est cette étape, combinée à la précédente, qui fait la spécificité du traitement des orthophotographies verticales sur MicMac, puisqu’elle permet d’obtenir un résultat mathématiquement fiable et reproductible sans aucune interaction manuelle sur le modèle. Elle consiste à recalculer un repère cartésien temporaire à partir du géoréférencement, sans modifier l’orientation du modèle. Vous avez le choix entre deux solutions, selon le résultat que vous souhaitez obtenir et/ou de l’objet dont vous faites le relevé. Le résultat dépendra de l’argument définissant l’orientation : HOR ou HORVy.

Commande

mm3d RepLocBascule ".*JPG" Ori-CAL-Terrain HOR Repere-elevation.xml PostPlan=_Masqplan

ou

mm3d RepLocBascule ".*JPG" Ori-CAL-Terrain HORVy Repere-elevation.xml PostPlan=_Masqplan

Description

32RepLocBascule – Cet outil permet de créer un repère temporaire sans changer l’orientation du modèle.

33Ori-CAL-Terrain – La commande se base sur le répertoire de calibration et d’orientations définitif.

34HOR – Cet argument rabat l’axe Ox à l’intersection entre le plan horizontal défini par le géoréférencement et le plan moyen défini par le masque à l’étape précédente. L’axe Oy devient perpendiculaire au plan moyen défini par le masque. Cet argument est pertinent si vous avez besoin que le « point de vue » se positionne face à la surface définie dans le masque. Cet argument est pertinent pour le relevé d’une coupe stratigraphique, ou toute autre surface qui ne serait pas parfaitement verticale (fig. 34).

Fig. 34 : schéma du fonctionnement théorique de l’argument HOR de l’outil RepLocBascule et exemple de résultat obtenu avec cette option.

35HORVy – Cet argument bascule l’axe Oy, qui devient l’axe vertical, tout en préservant l’horizontalité. Ox est aussi rabattu à l’intersection entre le plan horizontal défini par le géoréférencement et le plan moyen défini par le masque comme avec l’argument HOR, mais HORVy modifie les coordonnées de Oy qui a alors les caractéristiques d’un axe strictement vertical. On peut en effet voir dans le fichier Repere-elevation.xml que ses coordonnées X Y et Z sont 0 0 et 1. Préférez cet argument pour dessiner un pierre-à-pierre ou pour toute documentation d’architecture technique normalisée (fig. 35).

36Repere-elevation.xml – Nom du fichier XML qui stockera les informations relatives au repère temporaire.

37PostPlan=_Masqplan – Cet argument permet d’appeler le masque produit à l’étape précédente en ne renseignant que le suffixe appliqué au nom de l’image.

Fig. 35 : schéma du fonctionnement théorique de l’argument HORVy de l’outil RepLocBascule et exemple de résultat obtenu avec cette option.

Production

38Cette commande produit un fichier XML : Repere-elevation.xml, qui recueille les informations du repère temporaire cartésien qui servira ensuite à produire les mosaïques de l’orthophotographie et dont la vue sera orientée selon vos besoin.

3.4.5. Production des orthophotographies

Commande

mm3d PIMs2MNT MicMac DoOrtho=1 Repere=Repere-elevation.xml

Description

39PIMs2MNT – Cet outil calcule dans un premier temps des cartes de profondeur pour chaque image (à partir de l’outil PIMs, pour Per Image Matching), permettant ainsi de produire un MNT.

40MicMac – Cet argument doit absolument correspondre à la définition choisie à l’étape 3 .4.2. Les différentes options possibles sont les mêmes que proposées pour l’outil C3DC.

41DoOrtho=1 – Cet argument demande à l’outil de recomposer une orthophotographie pour chaque photo, car il ne le fait pas par défaut.

42Repere=Repere-elevation.xml – Sans cet argument, le nouveau repère local précédemment calculé n’est pas pris en compte.

Production

43Cette commande produit quatre nouveaux répertoires. Les fichiers utiles au calcul suivant sont directement stockés dans PIMs-ORTHO. Il réunit toutes les mosaïques de l’orthophotographie et la carte de profondeur de chaque image ainsi que le fichier de géoréférencement World (.tfw) et XML correspondant à chaque mosaïque. PIMs-TmpBasc, PIMs-TmpMnt et PIMs-TmpMntOrtho sont quant à eux des répertoires temporaires.

3.4.6. Si plusieurs nuages de points sont modélisés en même temps (optionnel)

Description

44Pat="elevation1.*JPG" – Cet argument, ajouté à la commande précédente, permet de calculer l’orthophotomosaïque pour une partie seulement du nuage de points, en ne prenant en compte que les images concernées.

À prendre en compte

45Cette option nécessite une préparation du set de données un peu plus laborieuse en amont : les photos de chaque partie doivent être nommées avec un préfixe. Il faut là encore être attentif à la casse, éviter les accents, espaces et autres caractères spéciaux. Nous avons mis ici le préfixe elevation1 à titre d’exemple, mais le nom choisi pour chaque mur lors de l’enregistrement est sûrement plus pertinent pour ne pas se perdre. La commande Tapas doit être lancée pour chaque partie si les réglages de l’appareil photo ont été changés, pour une meilleure profondeur de champ lors de l’acquisition d’angles sortants, par exemple. Il est aussi conseillé, pour de meilleurs résultats, de ne pas prendre les photos des angles en considération dans le calcul. En effet, l’orthophotomosaïque produite risque alors d’être étirée à cause de l’angle d’incidence du point de vue contre la surface de l’objet.

46Nous déconseillons, dans le cadre d’un grand nombre d’orthophotomosaïques d’élévations à calculer, de modéliser l’ensemble de l’objet d’étude. Cela demandera une grande puissance de calcul, pour un gain de temps incertain : traiter chaque élévation individuellement reste la solution la plus optimale si l’on ne dispose pas d’une station de travail (très) performante et d’une bonne capacité de stockage. Cela peut néanmoins être envisageable si le bâtiment est petit ou s’il faut relever un petit groupe d’élévations contiguës et que ces élévations présentent une surface réduite acquise avec peu de clichés (moins de 50 photos). Si le site nécessite impérativement une modélisation 3D complète, il est alors sûrement plus pertinent de se tourner vers une acquisition au scanner laser terrestre plutôt que photogrammétrique.

47Cet argument est indispensable si vous modélisez plusieurs élévations en même temps, sinon, PIMs2MNT ne fonctionnera pas. En effet, le calcul de l’orthophotographie n’est possible qu’avec des photos présentant un angle d’incidence maximum de la caméra par rapport au plan moyen défini à l’étape 3.4.3. Il est prédéfini dans PIMs2MNT à 1 radian (environ 57°) et est non modifiable dans la commande. Par conséquent, toute photo du set dépassant cet angle d’incidence maximum empêchera PIMs2MNT de calculer l’orthophotographie.

3.4.7. Production de l’orthophotomosaïque

Commande

mm3d Tawny PIMs-ORTHO/

Description

48Tawny – Cet outil assemble automatiquement les orthophotographies pour produire une orthophotomosaïque.

49PIMs-ORTHO/ – Indique dans quel répertoire se trouvent les orthophotographies à assembler.

Production

50Cette commande produit deux fichiers dans le répertoire PIMs-ORTHO : Orthophotomosaic.tif, le résultat final du traitement (fig. 36), et Orthophotomosaic.tfw. Ce dernier est le fichier de coordonnées World de l’orthophotomosaïque. La première ligne indique la résolution au sol exprimée en mètres. Les deux dernières lignes correspondent aux coordonnées tronquées X et Y du pixel dans l’angle bas gauche de l’image.

Fig. 36 : résultat brut du traitement. Mur gouttereau Nord de l’église paléochrétienne de Mirine, Omišalj, Croatie, 2021.

À prendre en compte

51Ici, avec la bascule du repère, ces coordonnées ne revêtent plus aucune réalité géographique, mais attention tout de même : le fichier TIFF doit toujours être accompagné, dans le même répertoire, du fichier World correspondant et posséder exactement le même nom (sauf l’extension). En effet, afin de pouvoir post-traiter l’orthophotomosaïque dans un SIG, le fichier World permet d’obtenir une image directement à taille réelle sans aucune manipulation de la part de l’utilisateur.

3.4.8. Assemblage des tuiles de l’orthophotomosaïque (optionnel)

Commande

mm3d ConvertIm PIMs-ORTHO/Orthophotomosaic.tif

Description

52ConvertIm – Lorsque l’orthophotomosaïque fait plus de 1 Go, MicMac découpe l’image en plusieurs tuiles. Cette commande permet de les réassembler automatiquement.

53PIMs-ORTHO/Orthophotomosaic.tif – Cet argument permet la recherche automatique des tuiles du fichier Orthophotomosaic.tif dans le répertoire PIMs-ORTHO.

À prendre en compte

54Cette commande produit le fichier Orthophotomosaic_Out.tif dans le répertoire PIMs-ORTHO. Il ne faut pas oublier de renommer le fichier World Orthophotomosaic.tfw en Orthophotomosaic_Out.tfw pour pouvoir exploiter les résultats du traitement par la suite.

3.5. Le cas particulier des coupes archéologiques

55Concernant la production des orthophotographies de coupes archéologiques, toutes les étapes précédentes sont communes au traitement pour la production d’orthophotographies en élévation et nous renvoyons au chapitre précédent. Nous reprenons ici seulement l’étape « 3.4.4. Calcul du repère temporaire », pour le cas des coupes archéologiques. C’est cette étape, combinée à la précédente (cf. 3.4.3), qui fait la spécificité du traitement des orthophotographies verticales sur MicMac, puisqu’elle permet d’obtenir un résultat mathématiquement fiable et reproductible sans aucune interaction manuelle sur le modèle. Elle consiste à recalculer un repère cartésien temporaire à partir du géoréférencement, sans modifier l’orientation du modèle.

Commande

mm3d RepLocBascule ".*JPG" Ori-CAL-Terrain HOR Repere-elevation.xml PostPlan=_Masqplan

Description

56RepLocBascule – Cet outil permet de créer un repère temporaire sans changer l’orientation du modèle.

57Ori-CAL-Terrain – La commande se base sur le répertoire de calibration et d’orientations définitif.

58HOR – Cet argument rabat l’axe Ox à l’intersection entre le plan horizontal défini par le géoréférencement et le plan moyen défini par le masque à l’étape précédente. L’axe Oy devient perpendiculaire au plan moyen défini par le masque. Cet argument est pertinent si vous avez besoin que le « point de vue » se positionne face à la surface définie dans le masque. Cet argument est pertinent pour le relevé d’une coupe stratigraphique, ou toute autre surface qui ne serait pas parfaitement verticale (fig. 30).

59Repere-elevation.xml – Nom du fichier XML qui stockera les informations relatives au repère temporaire.

60PostPlan=_Masqplan – Cet argument permet d’appeler le masque produit à l’étape précédente en ne renseignant que le suffixe appliqué au nom de l’image.

Production

61Cette commande produit un fichier XML : Repere-elevation.xml, qui recueille les informations du repère temporaire cartésien qui servira ensuite à produire les mosaïques de l’orthophotographie et dont la vue sera orientée selon vos besoin.

À prendre en compte

62La suite de la chaîne de commande se poursuit normalement à partir de l’étape 3.4.5. Pour plus de renseignements sur leur fonctionnement, les commandes et les arguments optionnels possibles, nous vous renvoyons au chapitre précédent.

mm3d PIMs2MNT MicMac DoOrtho=1 Repere=Repere-elevation.xml

mm3d Tawny PIMs-ORTHO/