# APPORT DE L'IMAGERIE ALSAT2 POUR LA MAINTENANCE DU CADASTRE STEPPIQUE EN ALGERIE

Mohamed Akram SEDDIKI¹,²,³ (Chargé de recherche), Sid Ahmed SOUIAH¹,² (Maitre de Conférence A), Tarek GHODBANI¹,²(Professeur, Directeur du Laboratoire EGEAT), Abd El Halim GUERROUDJ⁴ (Assistant d'études A), Fatiha ZERGUI³ (Ingénieur d'état), Mohamed MILOUA⁵ (Ingénieur d'état).

Akramcnts@yahoo.fr ou aseddiki@cts.asal.dz

N°01, Avenue de la Palestine, BP13, Arzew 31200, Oran, Algérie, Télé : +21341793042, Fax : +21341792176

1. Faculté des Sciences de la terre et de l'Univers, Département de Géographie, Université Oran2, Algérie

2. Laboratoire EGEAT Espace Géographique et Aménagement du Territoire, Université Oran2, Algérie

3. Centre des Techniques Spatiales/Agence Spatiale Algérienne, Algérie

4. Université Ammar THLIDJI, Laghouat, Algérie

5. Institut National de Cartographie et de Télédétection INCT, Algérie

## Résumé

L'établissement du cadastre général en Algérie et la constitution du livre foncier ont été entrepris dès 1977. Cette opération a concerné, en premier lieu, les régions nord du territoire national (les zones rurales et les zones urbaines) et les agglomérations des régions sud.

Ce projet porte sur un vaste territoire couvrant les zones dites 'steppiques' et s'inscrit dans le cadre de recherche de solutions techniques relatives à l'exploitation de l'imagerie du satellite ALSAT2 pour la maintenance des bases de données graphiques du cadastre steppique. Cela par la détermination des capacités informatives que possède cette image en matière d'identification des objets correspondant à une partie ou la totalité des propriétés foncières. Les informations extraites feront partie intégrante d'un Système d'Informations Foncières multi-sources SIF (données du cadastre, d'urbanisme et base de données ortho-image).

**Mots clés:** Imagerie ALSAT2, Algérie, Cadastre Steppique, SIF, classification K-Means, ISODATA, Orientée objet.

## Abstract

The establishment of the general cadastre in Algeria and the constitution of the land register were undertaken, in 1977 and first of all, in the northern regions of the national territory, in rural and urban areas, as well as in the urban areas of the southern regions.

This project covers a vast territory covering the so-called 'steppe' areas and is part of a research framework for technical solutions relating to the exploitation of the imagery of the ALSAT2 satellite for the maintenance of the graphic databases of the steppe land register, this, by determining the informative capacities that this image possesses in terms of identifying objects corresponding to part or all of the land holdings. The information extracted will be an integral part of a multi-source Land Information System SIF (cadastral data, planning data and ortho-image database).

**Keywords:** ALSAT2 Imagery, Algeria, Steppique Cadster, SIF, K-Means classification, ISODATA, Oriented Object.

## 1. Introduction

En Algérie, comme défini dans la loi n°90-25 du 18/11/1990 portant orientation foncière, une terre saharienne est toute terre située en-dessous de l'isohyète de 100mm. Les terres steppiques sont constituées, au sens de cette même loi, par des terres à vocation pastorale ou alfatière (figure 01). Compte tenu de l'étendue des superficies des zones steppiques et de la nature des territoires à traiter, qui constituent dans leur quasi-totalité des zones ayant un statut de propriété publique, les procédés techniques choisis pour la réalisation du cadastre des zones nord du territoire national (délimitation des propriétés basée sur les orthophotoplans et les plans de restitution) sont devenus inadaptés.

Cependant, la délimitation des propriétés en zones steppiques s'effectue à travers des travaux de bureau en se basant sur des documents cartographiques (généralement la couverture cartographique existante est à petite échelle) ce qui donne lieu à la confection d'une matrice cadastrale. Cette dernière comporte la liste des propriétés de l'Etat et un plan cadastral pour constituer le Livre Foncier au profit de l'Etat (des titres de

propriété sous forme de livrets fonciers). Ce genre de délimitation au bureau par des supports cartographiques (sans recourir à des images satellitaires ni à une enquête sur le terrain) présente un défaut d'omission des propriétés privées et des terres 'Arch' [terres dont les tribus jouissaient en permanence, collectivement, suivant les besoins de la communauté. Gérées par un conseil composé de notables de la tribu, elles étaient juridiquement enregistrées au compte de l'Etat depuis la réforme agraire en 1972.] (R. Bensouiah, 2003 et A Ali, 2011).

Ce projet s'inscrit dans un cadre de recherche de solutions techniques relatives à l'exploitation de l'imagerie du satellite ALSAT2 (construit dans le cadre d'un partenariat entre l'ASAL-Algérie et EADS Astrium-France (Herbert J. Kramer, 2002) pour la maintenance du cadastre steppique. Cela par la détermination des capacités informatives que possède cette image en matière d'identification des objets correspondant à une partie ou la totalité des propriétés foncières. Cette question d'identification, à partir des images à haute et très haute résolution spatiale, a été abordée par
plusieurs auteurs et consiste à explorer son contenu et sa précision (M. Ettarid, 2004, S. L-Homme, 2005) ; en précisant en outre la détermination de l'échelle géométrique, la base de données des objets interprétés et sa correspondance ou non aux différentes échelles (D. Holland, 2006).

La comparaison des différentes méthodes d'extraction de l'information par classification à partir des images spatiales et la réduction du taux d'incertitude (A. Masse, 2013), permettent aux services du cadastre (par l'utilisation des capacités du satellite ALSAT2) de mettre à jour les bases de données graphiques géoréférencées après un recensement précis des îlots à l'échelle des sections cadastrales et des communes (enregistrés juridiquement au nom de l'état alors qu'ils appartiennent réellement à des propriétaires privés, des usufruitiers ou des occupants des terres Arch.

Les premières solutions testées dans le cadre de ce travail avaient pour base l'utilisation des algorithmes de classification standards qui sont adaptées aux parcelles de taille grande ou moyennes par rapport à la résolution spatiale de l'image support d'étude. Pour les îlots composés de bâtiments isolés de petites dimensions, une approche d'apprentissage a été développée.

<figure><img src="image_1.png" /><figcaption>**Figure 1** : Localisation du territoire steppique sur la carte de l'Algérie (14% du territoire national).</figcaption></figure>

## 2. Données et zone d'étude

La wilaya d'El Bayadh fait partie intégrante de la région des Hautes Plaines steppiques du Sud-Ouest algérien. Sur le plan physique, elle présente trois grandes zones distinctes :

* Au nord : les Hautes Plaines ;

* Au centre : l'Atlas saharien ;

* Au sud : la frontière présaharienne.

Dans notre étude, pour des raisons de disponibilité des données, on a choisi deux communes : La commune d'El Bayadh d'une superficie de 48783 Ha et la commune de Rogassa qui couvre environ 241600 Ha, délimitées géographiquement par les parallèles 33°,26 et 34°,16 Nord, et les méridiens 0°,73 et 1°,23 Est.

<figure><img src="image_2.png" /><figcaption>**Figure 2** : Ortho-image ALSAT2 drapée sur le MNT de la ville d'El Bayadh (altitude moyenne de 1313 m).</figcaption></figure>

Concernant les données images, sur un nombre de onze (11) scènes, qui couvrent la zone d'étude, sept (07) images ont été utilisées pour produire le support ortho-image dont une scène utilisée pour la détermination de l'échelle par étude du potentiel d'identification des objets. Il s'agit d'une évaluation du potentiel cartographique de l'imagerie satellite qui doit se faire sur la base de critères permettant de conférer au produit les qualités suivantes:
* La possibilité de pouvoir identifier sans ambiguïté, sur l'image, le type et la nature des objets à représenter sur une carte à une échelle donnée.
* La précision de la localisation de ces objets identifiés. (M. Ettarid, 2004).

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 4** : Zoom sur la partie encadrée en bleu sur le segment ALSAT2 (figure 3).</figcaption></figure>

<table><thead><tr><th>Image ALSAT2</th><th>Caractéristiques</th></tr></thead><tbody><tr><td>Nom de l'image</td><td>A1M 2017-06-22 09:58:23.6338960</td></tr><tr><td>Dimension de la scène (lignes et colonnes)</td><td>1750 x 1750</td></tr><tr><td>Format de stockage</td><td>GEOTIFF</td></tr><tr><td>Date de prise de vue</td><td>22-06-2017</td></tr><tr><td>Heure de prise de vue</td><td>09 h 58 '23"</td></tr><tr><td>Système de coordonnées géographiques</td><td>WGS 84</td></tr><tr><td>Nombre de bandes spectrales</td><td>4 Bandes</td></tr><tr><td>Résolution spatiale (mètre)</td><td>2,5 x 2,5</td></tr><tr><td>Localisation géographique</td><td>λmin = 00°51' 40" E<br>φmin= 33°34' 48" N<br>λmax = 01°05' 22"E<br>φmax = 33°46' 07"N</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 5 :** Spatiocarte à partir d'une image ALSAT2 multispectrale de la ville d'El Bayadh.</figcaption></figure>

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 7** : Spatiocarte de la zone d'enquête cadastrale (xmin=307000, ymin= 3731000, xmax=314000, ymax=3735000m). pertinent d'appréciation de l'omission de la présence des propriétés recherchées (propriétés</figcaption></figure>

Pour démontrer l'efficacité des images ALSAT2 en matière d'extraction de l'information (utile pour la mise à jour du cadastre steppique) un jeu de données des bases graphiques Version 1.6 (dernière mise à jour en 2013) a été récupéré auprès de l'Agence Nationale du Cadastre qui couvre la commune d'El Bayadh (chef-lieu de la wilaya) et la commune de Rogassa.

privées, de l'Etat occupées dans le cadre de concessions ou Arch). A cet effet, nous avons opté pour des méthodes qui intègrent, en plus des dimensions spectrales, les informations texturales et géométriques, parmi lesquelles figure l'approche orientée objet (A. Puissant, 2003). La méthodologie suivie repose sur le principe de tester plusieurs classifieurs qui donnent généralement des informations certaines (redondantes) et d'autres en confusion (complémentaires). Le but recherché est la formulation de règles de connaissances permettant l'identification et l'extraction des éléments pertinents à la maintenance du cadastre steppique dans le cadre d'une démarche de classification orientée objet.

Durant tous les traitements, le choix des échantillons était en adéquation avec la nomenclature Spot Théma définie à une résolution spatiale équivalente à l'image support d'étude. (Spot Théma est un produit d'occupation du sol à l'échelle des agglomérations, proposé sur la France métropolitaine et décliné en 8, 29 et 46 classes. Il s'agit d'une base de données vectorielles produite à partir de l'interprétation d'images satellites SPOT 1 à 5 ou équivalent, d'une résolution de 2,5 et 20m, Astrium GEO Information Services, 2011).

Ainsi la notion de zone exclue (omise) du cadastre steppique a été redéfinie dans le cadre de ce travail de manière à prendre en considération toute possibilité d'existence de propriété privée dans le territoire cadastré. Initialement cette notion portait sur les agglomérations, les zones comportant des exploitations agricoles ou tous les établissements (quelle qu'en soit la

<figure><img src="image_8.png" /><figcaption>**Figure 8** : Références cadastrales des sites choisis (Xcentroide=311398, Ycentroide=3732773m).</figcaption></figure>

### 3. Méthodologie et résultats

La disponibilité de données haute résolution à partir d'une plate-forme satellitaire a ouvert de nouvelles possibilités pour la cartographie et la mise à jour du cadastre (P. Jayaprasad et al, 2006). Le paysage steppique est marqué par la présence de détails inférieurs au seuil théorique d'extraction de plusieurs objets, par exemple l'objet bâti isolé qui est un élément
nature) qui sont seulement «délimités» suivant une
vectorisation de leurs contours en raison de la densité
des unités foncières qui la composent. Elles figurent
dans le plan cadastral en tant que périmètres non soumis

aux opérations du cadastre saharien ou steppique. Elles
seront cadastrées ultérieurement selon les modalités
techniques habituelles basées sur les orthophotoplans et
restitutions photogrammétriques.

<figure><img src="image_9.png" /><figcaption>**Figure 9** : Répartition spatiale sur l'image ALSAT2 des échantillons choisis selon la nomenclature d'objets Spot-Théma (Système UTM WGS 84 Zone 31 Nord, xmin=301500, ymin= 3720000, xmax=323000, ymax= 3740000m).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 10** : Exemples de propriétés spectrales de quelques échantillons choisis selon la nomenclature d'objets Spot-Théma (leur identification avec agrandissement à partir de l'image ALSAT2 et sur images THRS type Google-Earth). Les cercles bleus correspondent à des objets bâties qui ne peuvent pas être considérés comme isolés en raison de leur proximité avec des routes ou leurs positions à l'intérieur de parcelles (facilement identifiables). (Système UTM WGS 84 Zone 31 Nord, point en haut x=315650, y= 3724500, point en bas x=316014, y= 3720501m).</figcaption></figure>
La première lecture des images ALSAT2 donne l'impression qu'il s'agit d'un paysage à trois grandes classes thématiques : urbain, sol et végétation. Alors qu'en réalité c'est un milieu (espace) marqué par la dominance du bâti isolé et des cultures qu'on peut faire correspondre à des jardins familiaux de la nomenclature Spot Théma (figure 9).

Notre objectif ici n'est pas de tester les méthodes de classification pour une production cartographique ordinaire, qui doit obéir aux normes de représentation et de définition des objets sur des plans, mais plutôt d'extraire l'information sur ces petits objets 'bâtis isolés (niveau d'identification suffisant) ou parcelles (cartographie à une échelle 1/10000 à 1/25000)' qui fourniront des informations sur la présence des

propriétés privées ou des biens étatiques occupés par des personnes physiques. Par souci d'atteindre une certaine qualité, nous avons veillé à la bonne définition des échantillons d'extraction. Ainsi, nous avons testé les méthodes standards dans le domaine des classifications d'image pour bien cerner la variation spectrale des différents objets. Ensuite, nous les avons intégrées dans un processus de fusion orientée objet pour avoir des résultats avec un taux minimum de confusion thématique.

L'organigramme présenté dans la figure 11 résume toute la chaine de traitement allant de l'image support d'étude en passant par la nomenclature des objets et la définition des zones exclues (omises) jusqu'à l'obtention des résultats et leur validation.

<figure><img src="image_11.png" /><figcaption>**Figure 11** : Approche méthodologique développée pour l'identification des zones d'intérêt (îlots de propriété bâti ou non bâti) dans le cadre de la maintenance du cadastre steppique</figcaption></figure>
## 3.1 Précision géométrique de l'imagerie ALSAT2

Plusieurs normes relatives à la notion de précision de la position existent dans la littérature (ex : P. A. Handbook, 1999) La norme NMAS (United States National Map Accuracy Standards) appliquées aux USA (F.G.D.Committee, 1998), qui est adaptée aux grandes étendues (cas de nos régions steppiques et sahariennes en Algérie). Elle a une relation entre l'indice CE90 et l'échelle de restitution. Elle est donnée pour des échelles plus grandes que 1/20 000 par :

$$E = CE90 \times 100 / (2,54 \times 0,03333)$$

Elle utilise l'indice CE90 dont l'erreur a un intervalle de confiance de 90% pour évaluer la qualité des documents. Cinq points GPS ont servi à la détermination de l'échelle de l'image ALSAT2. Ces points sont exclus du calcul de l'orientation externe.

<table><caption>**Tableau 2 :** Coordonnées des points GPS et Orthoimage ALSAT2 utilisées pour la détermination de l'échelle.</caption><thead><tr><th colspan="2">Coordonnées GPS</th><th colspan="2">Coordonnées Ortho-image ALSAT2</th></tr><tr><th>X</th><th>Y</th><th>X</th><th>Y</th></tr></thead><tbody><tr><td>318892.67</td><td>3733998.59</td><td>318890.00</td><td>3734001.89</td></tr><tr><td>319194.83</td><td>3722897.33</td><td>319197.43</td><td>3722894.64</td></tr><tr><td>309539.50</td><td>3719456.88</td><td>309537.75</td><td>3719456.87</td></tr><tr><td>315693.63</td><td>3716341.13</td><td>315693.79</td><td>3716343.18</td></tr><tr><td>312979.54</td><td>3709923.89</td><td>312980.97</td><td>3709921.24</td></tr></tbody></table>

<table><caption>**Tableau 3 :** Détermination de l'échelle cartographique de l'image ALSAT2 support d'étude.</caption><thead><tr><th>Ecart X</th><th>Ecart Y</th><th>Ecart D</th><th>σX</th><th>σY</th><th>σI</th></tr></thead><tbody><tr><td>2,675</td><td>-3,298</td><td>4,24</td><td>2,18</td><td>2,71</td><td>3,48</td></tr><tr><td>-2,597</td><td>2,697</td><td>3,74</td><td></td><td></td><td></td></tr><tr><td>1,755</td><td>0,014</td><td>1,75</td><td></td><td></td><td></td></tr><tr><td>-0,153</td><td>-2,048</td><td>2,05</td><td></td><td></td><td></td></tr><tr><td>-1,422</td><td>2,651</td><td>3,00</td><td></td><td></td><td></td></tr></tbody><tfoot><tr><td colspan="6">CE90= 7,46 et E=8816</td></tr></tfoot></table>

L'échelle déterminée par application de la norme NMAS est 1/E=1/9000 et qui ne concerne que la précision de position sans relation avec l'identification des objets qui dépend de la résolution spatiale de l'image.

## 3.2 Capacité informative de l'imagerie ALSAT2

Une interprétation spectrale a été faite en examinant plusieurs échantillons présentés :

- **Au niveau des statistiques par pixel** : sont illustrées par le tracé des signatures spectrales dans les quatre bandes (bleu, vert, rouge et proche infrarouge) et ne représentent qu'une variation des pixels isolés (figure10) sans tenir compte des variations du voisinage et de la totalité des pixels qui forme l'objet. Les signatures confirment l'appartenance des pixels à la même classe thématique mais n'indiquent pas la confusion spectrale entre les pixels testés.

- **Au niveau des statistiques par objet** : contrairement au cas précédent les statistiques zonales ou par objet peuvent révéler l'hétérogénéité spectrale de l'objet et son voisinage. Cela dépend d'une part, du nombre de pixels qui le composent, et d'autre part, de la nature de l'objet lui-même.

### 3.2.1 Extraction supervisée Mono bande

Dans le but de traiter les bandes spectrales individuellement en fonction des propriétés intrinsèques de chaque objet, la technique d'extraction est personnalisée. Le travail se fait manuellement et nécessite des connaissances préalables sur le comportement spectral de chaque objet dans chaque bande. Le calcul d'indice tel que le NDVI permet une séparation thématique qui peut servir à une réduction de la probabilité de confusion interclasses.

<figure><img src="image_12.png" /><figcaption>**Figure 12 :** Exemple de confusions thématiques (pixels en jaune) dues à l'extraction mono-bande ($x_{centroide}=312623$, $y_{centroide}=3729383m$).</figcaption></figure>

### 3.2.2 Extraction Multi-bandes

#### Extraction non supervisée

Dans cette phase de traitements deux algorithmes standards sont utilisés (celui d'Isodata et du K-Means) afin de mettre en évidence les variations spectrales des deux objets avec un nombre de classes : 5, sur 10 et 15.
La comparaison des deux algorithmes de classification non supervisée utilisés nous a permis de constater :

- Une extraction moins confuse pour l'objet végétation,

- Une extraction de l'objet bâti avec un taux de confusion plus ou moins important et à des degrés différents,

- Pour les deux algorithmes l'augmentation du nombre de classes améliore légèrement la qualité mais en même temps change les seuils des confusions,

La qualité d'extraction de l'objet bâti par algorithme d'Isodata est relativement meilleure que celle obtenue par l'algorithme K-Means.

<figure><img src="image_13.png" /><figcaption>**Figure 13** : Exemple de confusions thématiques (pixels en jaune) dues à l'extraction par application de l'algorithme K-Means (xcentroide=312819m, ycentroide=3729944m).</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 14** : Exemple de confusions thématiques (pixels en jaune) dues à l'extraction par application de l'algorithme ISODATA (xcentroide=312819m, ycentroide=3729944m).</figcaption></figure>

## Extraction par classification supervisée

Pour choisir les échantillons de classification deux contrôles sont nécessaires :

- Son identité spectrale : par le tracé de la signature qui confirme son appartenance,

- Son histogramme : qui montre s'il est en confusion avec d'autres objets ou non.

L'algorithme choisi pour lancer la classification est celui du maximum de vraisemblance. Contrairement au cas de la végétation, l'objet bâti présente un taux de confusion important. Les opérations de nettoyage se limitent à des simples requêtes sur les superficies.

<figure><img src="image_15.png" /><figcaption>**Figure 15** : Exemple de confusions thématiques (pixels en rouge) dues à l'extraction par application de la classification supervisée par maximum de vraisemblance. (xmin=301600, ymin=3719520, xmax=322926, ymax=3735924m).</figcaption></figure>

## Synthèse

L'analyse des techniques d'extraction testées montre que pour chaque résultat il y a une information certaine sur les deux objets (végétation et bâti). En revanche l'incertitude est fortement liée aux dimensions de l'objet au sol et à sa nature (voir annexe).

L'analyse de l'extraction brute de l'objet bâti (illustrée par le graphe de la figure 16) confirme l'hypothèse que le résultat représente un maximum d'objets avec une forte probabilité de confusion.

Cette hypothèse n'est pas toujours fiable pour le cas de la végétation. Néanmoins, l'hétérogénéité des espèces rend souvent la confusion significative.
<figure><img src="image_16.png" /><figcaption>**Figure 16** : Statistiques des nombres d'objets bâties extraits par différents algorithmes de classification</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 17** : Comparaison des statistiques d'extraction de l'objet végétation</figcaption></figure>

### 3.2.3 Extraction par approche orientée objet

Les objets avec les mêmes revêtements nous fourniront une même signature spectrale sur une image satellitaire mais seront reconnaissables par leurs formes et textures respectives. Ce problème est lié à la variabilité radiométrique, spatiale et spectrale.

Ainsi, la finesse de la résolution spatiale des images haute résolution telle que ALSAT2, exige de travailler

dans un espace à plusieurs dimensions et non pas uniquement limité à la seule dimension spectrale.

Dans cette partie, nous avons opté pour une solution permettant l'intégration d'un certain nombre de règles de connaissances afin de guider l'algorithme à une identification précise de l'objet à extraire. Ces connaissances sont définies pour chaque segment à des échelles différentes.
<figure><img src="image_18.png" /><figcaption>**Figure 18** : Agrandissement de la Spatiocarte ALSAT2 au 1/2500 sur laquelle sont encerclés en bleu des objets bâties et des parcelles omis qui appartiennent à l'îlot 3 et 43 cadastrés au nom de l'état.</figcaption></figure>

Les imagettes illustrées ci-dessous sont extraites à partir de la Spatiocarte pour mettre en exergue les objets extraits (végétaux en couleur verte et bâti en magenta).

<figure><img src="image_19.png" /><figcaption>**Figure 19** : Exemples de parcelles extraites A (parcelles dispersées), B (Exploitations agricoles collectives EAC), C (végétation périurbaine). ($x_{centroide}=311700m, y_{centroide}=3728632m$).</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 20** : Zoom sur la classe A (en couleur verte) niveau d'extraction équivalent à l'échelle 1/10000 ($x_{centroide}=307120m, y_{centroide}=3732058m$).</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Figure 21** : Zoom sur la classe B (en couleur verte) niveau d'extraction équivalent à l'échelle 1/10000 ($x_{centroide}=312164$m, $y_{centroide}=3732177$m).</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 22** : Exemples du bâti identifié (sous forme de groupement de type industriel ou autre, c'est la classe 01, ou isolé classe 02) ($x_{centroide}=315500$m, $y_{centroide}=3724045$m).</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Figure 23** : Exemples de bâtis isolés extraits classe 2 sur la figure 22 (objets encadrés en bleu sur la figure 10). ($x_{centroide}=316200$m, $y_{centroide}=3722751$m).</figcaption></figure>

L'extraction en couleur magenta est qualitative. Le niveau de zoom de l'image est à l'échelle 1/7500 mais les objets extraits ne sont pas représentés dans les vraies dimensions. En agrandissant on perd la qualité visuelle de l'image. Pour les échelles plus grandes que 1/5000 la visualisation devient floue (figures 24). Mais nous pouvons toujours identifier des segments de bâtis ou de routes (figure 25).

Ces résultats corroborent des résultats similaires obtenus dans d'autres études réalisées sur la base de l'imagerie SPOT5. La précision de l'identification dépend directement de la taille, de la forme du bien, de la topographie de la zone, du type de clôtures et de la couverture végétale présente sur la zone d'étude ainsi que de l'échelle de l'ortho-image utilisée dans le processus d'identification.

L'orthoimage SPOT 5 est considérée alors comme une alternative intéressante pour la méthodologie indirecte avec une faible variation topographique. Cependant, il est recommandé d'utiliser des méthodes d'arpentage traditionnelles pour les petites propriétés, le cadastre urbain et l'identification des bâtiments (M. Corlazzoli, 2004).

Pour une échelle de restitution de l'ordre du 10000ème, la résolution spatiale de l'image SPOT5 convient sauf pour la détection de certaines entités de petite taille et de texture particulière. (F. Gourmelon, 2005).

### 3.2.4 Validation et discussion des résultats

En plus des statistiques obtenues par application des différentes méthodes d'extraction (tableaux en annexe) et en dehors des avantages et des limites de chacune des méthodes, la validation des résultats est effectuée par leur confrontation à la réalité lors des sorties sur terrain. Les endroits visités correspondent aux références cadastrales suivantes :

- Commune Rogassa, Section 14 (îlots 05, 06 et 07) : propriétés identifiées, visitées et enquêtées.

- Commune El Bayadh, Section 01 (îlot n°01 et 08) : propriétés identifiées, visitées et enquêtées.

La figure 26 montre des exemples sur la présence des propriétés privées et des propriétés Arch (appartenant à la même tribu) situées dans la commune de Rograssa. De même, la figure 27, montre des exemples sur la présence des propriétés Arch situant dans la commune El Bayadh.
<figure><img src="image_24.png" /><figcaption>**Figure 24** : Dégradation de la qualité visuelle de l'ortho-image à l'échelle 1/2500 (niveau de zoom différent de celui de la figure 18).</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>**Figure 25** : Identification des segments de bâtis et de routes (en couleur bleu) à partir de l'image ALSAT2 affichée à l'échelle 1/2500.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>**Figure 26** : Photos de terrains prises dans la commune de Rogassa (validation des extractions, coordonnées du centre de la parcelle x=307000m, y=3736000, bâti en jaune x=307244m, y=3736145).</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>**Figure 27** : Photos de terrains prises dans la commune d'El Bayadh (validation des extractions, Coordonnées du centre x=312964m, y=3734774m).</figcaption></figure>

**3.3 Conclusion**

La méthodologie suivie dans le cadre de ce travail repose sur le principe que chacun des classifieurs testés livre des informations certaines et d'autres en confusion (A. Masse, 2013), et sur la constitution d'une base de connaissances structurée comme objectif à atteindre avec l'approche orientée objet.

Durant tous les traitements, le choix des échantillons était en adéquation avec la nomenclature Spot Théma définie à une résolution spatiale équivalente à l'image support d'étude. Ainsi la notion de zone à exclure a été redéfinie de manière à intégrer toute possibilité d'existence de propriété privée dans le territoire cadastré.
L'approche suivie pour extraire les informations utiles comprend trois phases :

- Identification des échantillons,

- Extraction et amélioration,

- Une reconnaissance d'objets extraits à partir des images THRS Google-Earth avant validation par les sorties sur terrain. Après enquête et délimitation (croquis de conservation ou document d'arpentage) et ensuite l'application des cas de mise à jour.

L'étude des propriétés spectrales nous a permis d'affiner la précision des choix des échantillons et de réduire la marge d'incertitude.

Contrairement à l'approche géométrique (qui a donné une valeur de 1/9000) l'analyse thématique (en utilisant la nomenclature de 13 objets sélectionnés à partir de Spot Théma) montre qu'une cartographie exhaustive à l'échelle 1/10000 n'est pas possible alors qu'on peut assurer l'obtention du 1/25000 sans aucune difficulté.

<table><thead><tr><th rowspan="2">Nomenclature</th><th colspan="2">Echelles</th></tr><tr><th>1/10000</th><th>1/25000</th></tr></thead><tbody><tr><td>Culture annuelle et prairie</td><td>x</td><td>x</td></tr><tr><td>Jardins familiaux</td><td>x</td><td>x</td></tr><tr><td>Equipement de loisir</td><td>x</td><td>x</td></tr><tr><td>Habitat discontinu collectif</td><td>-</td><td>x</td></tr></tbody></table>

<table><caption>**Tableau 4 :** Exemple d'objets qui peuvent être cartographiés à partir de l'imagerie ALSAT2 aux échelles 1/10000 et 1/20000.</caption><thead><tr><th>Habitat discontinu pavillonnaire</th><th>-</th><th>x</th></tr></thead><tbody><tr><td>Espace libre urbain</td><td>x</td><td>x</td></tr><tr><td>Coupe forestière et jeune plantation</td><td>-</td><td>x</td></tr><tr><td>Culture permanente</td><td>x</td><td>x</td></tr><tr><td>Installation sportive</td><td>x</td><td>x</td></tr><tr><td>Habitat isolé</td><td>-</td><td>x</td></tr><tr><td>Grand équipement</td><td>x</td><td>x</td></tr><tr><td>Tissu urbain continu</td><td>x</td><td>x</td></tr><tr><td>Boisement linéaire</td><td>-</td><td>x</td></tr></tbody></table>

L'imagerie ALSAT2 peut répondre aux attentes du cadastre en zones steppiques en ce qui concerne l'identification des zones à exclure et même la cartographie à l'échelle 1/10000 à 1/25000. A une plus grande échelle il faut vérifier la possibilité d'intégrer des objets extraits directement dans les bases graphiques. Ainsi sur un échantillon de 17 objets de type bâti isolé en calculant les écarts en surface par rapport à la règle du 1/20ème (expliquée ci-dessous) seuls 12 d'entre eux (cf. tableau n°5) ont des détails dont la superficie dépasse 200m² et peuvent être intégrés directement dans les îlots cadastraux.

<figure><img src="image_28.png" /><figcaption>**Figure 28 :** Parfaite superposition des îlots cadastraux (en couleur jaune) sur la Spatiocarte à l'échelle 1/10000, la zone encadrée en bleu correspond à un exemple de la difficulté de l'image ALSAT2 pour la délimitation de l'habitat pavillonnaire discontinu cité en bleu sur le tableau 4.</figcaption></figure>
<table><caption>**Tableau 5** : Contrôle de l'application de la règle du 1/20ème relative à la délivrance des livrets fonciers.</caption><thead><tr><th>Superficie ortho-image (m2)</th><th>Superficie image classifiée (m2)</th><th>Ecart (m2)</th><th>1/20 de la surface</th></tr></thead><tbody><tr><td>444</td><td>418</td><td>26</td><td>22.2</td></tr><tr><td>158</td><td>141</td><td>17</td><td>7.9</td></tr><tr><td>291</td><td>296</td><td>-5</td><td>14.55</td></tr><tr><td>205</td><td>175</td><td>30</td><td>10.25</td></tr><tr><td>278</td><td>319</td><td>-41</td><td>13.9</td></tr><tr><td>332</td><td>339</td><td>-7</td><td>16.6</td></tr><tr><td>99</td><td>82</td><td>17</td><td>4.95</td></tr><tr><td>368</td><td>336</td><td>32</td><td>18.4</td></tr><tr><td>134</td><td>100</td><td>34</td><td>6.7</td></tr><tr><td>147</td><td>100</td><td>47</td><td>7.35</td></tr><tr><td>113</td><td>95</td><td>18</td><td>5.65</td></tr><tr><td>236</td><td>215</td><td>21</td><td>11.8</td></tr><tr><td>278</td><td>221</td><td>57</td><td>13.9</td></tr><tr><td>258</td><td>254</td><td>4</td><td>12.9</td></tr><tr><td>543</td><td>526</td><td>17</td><td>27.15</td></tr><tr><td>548</td><td>540</td><td>8</td><td>27.4</td></tr><tr><td>553</td><td>536</td><td>17</td><td>27.65</td></tr></tbody></table>

En Algérie la règle du 1/20ème doit être respectée pour la délivrance des livrets fonciers aux propriétaires ou les ayants droit. Elle consiste à une vérification de l'écart (qui ne doit pas dépasser 1/20) entre la superficie cadastrée et celle déclarée ou mentionnée sur l'acte.

Les images ALSAT2 permettent non seulement l'identification des parcelles et des objets existants (hangars, serres agricoles, bassins d'irrigation...) mais elles facilitent également l'accès sur le terrain aux propriétés et à leurs délimitations. Elles constituent un support terrain pour effectuer les opérations de maintenance des bases graphiques du cadastre plus facilement, avec une cartographie non exhaustive à l'échelle 1/10000 et complète à l'échelle 1/25000 (grands îlots appartenant à l'état).

Pour le cas des classes « bâti » la capacité de l'image est limitée à l'identification. Elle permet de localiser spatialement des segments indicatifs sur la possibilité de la présence d'une propriété omise (une confirmation est souhaitable avec l'imagerie THRS Google-Earth avant de sortir sur terrain).

La figure 29 donne l'exemple de l'îlot 68 (dont l'îlot primitif porte le numéro 05, qui va faire l'objet d'une mise à jour après enquête foncière. Il faut signaler que l'avantage de l'utilisation de l'imagerie ALSAT2 constitue une adaptation en matière du flux d'information entre le cadastre et la conservation foncière pour mieux gérer l'opération de délivrance des livrets fonciers.

<figure><img src="image_29.png" /><figcaption>**Figure 29** : Endroits de nos sorties sur terrain pour validation des résultats, cinq (05) propriétaires étaient présents sur le terrain et indiqués sur le tableau. Les propriétés enquêtées sont de nature juridique 'privées individuelles'.</figcaption></figure>
Enfin, la méthodologie développée durant ce travail a démontré l'intérêt de la télédétection opérationnelle dans le domaine du foncier à travers l'utilisation des ortho-images ALSAT2. Les résultats obtenus répondent aux exigences du cadastre steppique et saharien en Algérie. Ils permettent la maintenance de ses bases graphiques à des coûts compétitifs et dans les délais requis

(comparativement à la photogrammétrie aérienne). Il ressort de cette expérience que l'intégration de l'imagerie ALSAT2 dans le processus de maintenance engagé par l'Agence Nationale du Cadastre pourra faire l'objet d'un partenariat avec l'Agence Spatiale Algérienne.

<figure><img src="image_30.png" /><figcaption>**Figure 30** : Parcelles restituées automatiquement après extraction par approche orientée objet (à intégrer dans les bases graphiques cadastrales à l'échelle 1/10000, système UTM WGS84 Zone 31 N).</figcaption></figure>
Ahmed Ali A, 2011 : « La législation foncière agricole en Algérie et les formes d'accès à la terre ». In: *Elloumi M. (ed.), Jouve A.-M. (ed.), Napoléone C. (ed.), Paoli J.C. (ed.). Régulation foncière et protection des terres agricoles en Méditerranée. Mon tpellier : CIHEAM, 2011. p. 35-51* (Option s Méditerranéennes: Série B. Etudes et Recherches; n 66)
, http://om.ciheam.org/article.php?IDPDF=801372

Anne Puissant, 2003 : « Information géographique et images à très haute résolution utilité et applications en milieu urbain ». Géographie. Université Louis Pasteur - Strasbourg I, 2003. Français. tel-00467474, https://tel.archives-ouvertes.fr/tel-00467474. P 291

Antoine Masse, 2013: « Développement et automatisation de méthodes de classification à partir de séries temporelles d'images de télédétection - Application aux changements d'occupation des sols et à l'estimation du bilan carbone. Océan, Atmosphère. Université Paul Sabatier - Toulouse III. Français. tel- 00921853v2

Astrium GEO Information Services, 2011 : «SPOT THEMA Base de Données d'Occupation du Sol pour les Acteurs du Territoire ». GEO035-0313, P2 que ce que Spot Théma ?

Caitlin Dempsey, 2000: «GIS Learning National Map Accuracy Standards (NMAS)». https://www.gislounge.com/national-map-accuracy-standards-nmas/.

D.Holland, 2006 : « Updating topographic mapping in Great Britain using imagery from high-resolution satellite sensors », *ISPRS Journal of Photogrammetry and Remote Sensing* 60(3):212 223, DOI: 10.1016/j.isprsjprs.2006.02.002.

F. G. D. Committee, 1998: « Geospatial Positioning Accuracy Standards, Part 3: National Standard for Spatial Data Accuracy». *Subcomm. Base Cartogr. Data*.

Françoise Gourmelon, 2005: « Potentialités de l'imagerie satellitaire spot 5pour la cartographie de la végétation ». *Revue CyberGeo*, https://doi.org/10.4000/cybergeo.3027.P 36.

LEGENDES EN ANGLAIS

Figure 1: Location of the steppe territory on the map of Algeria (14% of the national territory).

Figure 2: AISAT2 ortho-image draped over the DTM of the city of El Bayadh (average altitude of 1313 m).

Figure 3: Example of a segment composed of two ALSAT2 images.

Herbert J. Kramer, 2002: "Observation of the Earth and Its Environment: Survey of Missions and Sensors", (Springer Verlag).
https://directory.eoportal.org/web/eoportal/satellite-missions/a/alsat-2.

Lhomme, Stéphane, 2005 : « Identification du bâti à partir d'images satellitaires à très hautes résolutions spatiales », Thèse en Géomatique appliquée, université de Sherbrooke. http://hdl.handle.net/11143/5824.

Loi n°90-25 du 18/11/1990, modifiée et complétée, portant orientation foncière. http://madrp.gov.dz/wpcontent/uploads/2019/01/L_oi_9025_18_11_1990 trigrid_1990 orientation foncière.pdf

M. Corlazzoli and O. L. Fernandez, 2004: « SPOT 5 Cadastral validation project in Izabal, Guatemala». *Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.*, vol. 35, no. Part 1, pp. 291-296.

M. Ettarid, f. Degaichia, 2004 : « Potentiel cartographique de l'imagerie Ikonos » Commission IV, WGIV/7. https://www.isprs.org/proceedings/xxxv/congress/comm4/papers/519.pdf.

P. A. Handbook, 1999: « Using the National Standards for Spatial Data Accuracy to measure and report geographic data quality». *Minn. Plan. Land Manag. Inf. Cent.*

P. Jayaprasad et al, 2006: « Updation and metrology of cadastre and village boundary using high resolution satellite data». *Int. Arch. Photogramm. REMOTE Sens.*, vol. 36.

Riad BENSOUIAH, 2003: « La lutte contre la désertification dans la steppe algérienne : les raisons de l'échec de la politique environnementale » *Communication présentée aux 15èmes Journées de la Société d'Ecologie Humaine organisée à Marseille (France) sous le thème : « Du Nord au Sud : le recours à l'environnement, le retour des paysans ?*», P 9 et 10.

Figure 4: Zoom in on the blue frame d area on the ALSAT2 segment (Figure 3).

Figure 5: Spatiomap from an ALSAT2 multispectral image of the city of El Bayadh.

Figure 6: Spatiomap of a panchromatic ALSAT2 image of El Bayadh airport (UTM system WGS 84 Zone 31 North, xmin=321000, ymin= 3732000, xmax=325500, ymax=3734000m).
Figure 7: Spatiomap of the cadastral survey area (xmin=307000, ymin= 3731000, xmax=314000, ymax=3735000m).

Figure 8: Cadastral references of the selected sites (UTM system WGS 84 Zone 31 North, x centroid=311398, y centroid=3732773m).

Figure 9: Spatial distribution on the ALSAT2 image of samples selected according to the nomenclature of Spot-Théma objects (UTM system WGS 84 Zone 31 North,xmin=301500,ymin= 3720000,xmax=323000,ymax= 3740000m).

Figure 10: Examples of spectral properties of some samples selected according to the nomenclature of Spot-Théma objects (their identification with enlargement from the ALSAT2 image and on VHSR images type Google-Earth). The blue circles correspond to built objects that cannot be considered as isolated because of their proximity to roads or their positions within parcels (easily identifiable).

Figure 11: Methodological approach developed for the identification of areas of interest (blocks of built or unbuilt property) as part of the maintenance of the steppe cadaster.

Figure 12: Example of thematic confusions (pixels in yellow) due to single-band extraction (x centroid=312623, y centroid=3729383m).

Figure 13: Example of thematic confusions (pixels in yellow) due to extraction by applying the K-Means algorithm (x centroid=312819m, y centroid=3729944m).

Figure 14: Example of thematic confusions (pixels in yellow) due to extraction by applying the ISODATA algorithm (x centroid=312819m, y centroid=3729944m). Figure 15: Example of thematic confusions (pixels in red) due to extraction by applying maximum likelihood supervised classification. (xmin=301600, ymin= 3719520, xmax=322926, ymax= 3735924m).

Figure 16: Statistics of the number of built objects extracted by different classification algorithms.

Figure 17: Comparison of the extraction statistics of the vegetation object.

Figure 18: Enlargement of the ALSAT2 Spatiomap at 1:2500 scale on which are circled in blue omitted built objects and parcels that belong to block 3 and 43 cadastred in the name of the state.

Figure 19: Examples of extracted plots A (scattered plots), B (collective farms EAC), C (peri-urban vegetation). (xcentroide=311700m, ycentroide=3728632m).

Figure 20: Zoom on class A (in green color) equivalent extraction level at 1/10000 scale (xcentroide=307120m, ycentroide=3732058m).

Figure 21: Zoom on class B (in green color) equivalent extraction level at 1/10000 scale (xcentroide=312164m, ycentroide=3732177m).

Figure 22: Examples of the identified building (in the form of grouping of industrial type or other, it is the class 01, or isolated class 02) (xcentroide=315500m, ycentroide=3724045m).

Figure 23: Examples of isolated buildings extracted class 2 on figure 22 (objects framed in blue on figure 10). (xcentroide=316200m, ycentroide=3722751m).

Figure 24: Degradation of the visual quality of the 1:2500 scale orthoimage (different zoom level than in Figure 18).

Figure 25: Identification of building and road segments (in blue) from the ALSAT2 image displayed at 1:2500 scale.

Figure 26: Photos of land taken in the municipality of Rogassa (validation of extractions, coordinates of the center of the parcel x=307000m, y=3736000, frame in yellow x=307244m, y=3736145).

Figure 27: Photos of land taken in the commune of El Bayadh (validation of extractions, center coordinates x=312964m, y=3734774m).

Figure 28: Perfect superposition of the cadastral property (in yellow color) on the Spati at 1/10000 scale, the area framed in blue corresponds to an example of the difficulty of the ALSAT2 image for the delimitation of the discontinuous suburban habitat cited in blue on Table 4.

Figure 29: Locations of our field trips for validation of results, five (05) owners were present in the field and indicated on the table. The properties surveyed are of a 'private individual' legal nature.

Figure 30: Plots of land automatically restored after extraction by object-oriented approach (to be integrated into the cadastral graphic databases at a scale of 1/10000, UTM WGS84 Zone 31 N system).

Table 1: Example of ALSAT2 image metadata

Table 2: GPS and ALSAT2 Orthoimage point coordinates used for scaling.

Table 3: Determination of the map scale of the ALSAT2 image used for the study.

Table 4: Example of objects that can be mapped from ALSAT2 imagery at 1/10000 and 1/20000 scales.

Table 5: Control of the application of the 1/20 meter rule for the issuance of land registers.

Table A1: Statistics of the extraction of the built object by single band analysis

Table A2: Statistics of the extraction of the built object by Isodata algorithm

Table A3: Statistics of the extraction of the built object by K-Means algorithm

Table A4: Statistics of the extraction of the built object by supervised classification (maximum likelihood)

Table A6: Summary of uncertainties
Annexe

<table><caption>**Tableau A1 : Statistiques de l'extraction de l'objet bâti par analyse mono bande**</caption><caption>Extraction supervisée Mono bande : objet bâti isolé</caption><thead><tr><th>Commune</th><th>îlot</th><th>superficie m2</th><th>Extraction brute</th><th>Après nettoyage</th><th>Confusion avec objet route</th><th>Confusion avec d'autres objets</th><th>Objets validés</th></tr></thead><tbody><tr><td rowspan="10">El Bayadh</td><td>2</td><td>353347</td><td>113</td><td>16</td><td>0</td><td>16</td><td>0</td></tr><tr><td>18</td><td>142945</td><td>116</td><td>11</td><td>0</td><td>9</td><td>2</td></tr><tr><td>4</td><td>35398</td><td>23</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td>5</td><td>446761</td><td>52</td><td>10</td><td>0</td><td>6</td><td>4</td></tr><tr><td>6</td><td>446761</td><td>41</td><td>7</td><td>0</td><td>6</td><td>1</td></tr><tr><td>7</td><td>2865137</td><td>353</td><td>21</td><td>0</td><td>13</td><td>8</td></tr><tr><td>12</td><td>148160</td><td>91</td><td>4</td><td>0</td><td>2</td><td>2</td></tr><tr><td>37</td><td>538384</td><td>81</td><td>10</td><td>8</td><td>0</td><td>2</td></tr><tr><td>39</td><td>1523879</td><td>112</td><td>35</td><td>35</td><td>0</td><td>0</td></tr><tr><td>40</td><td>209450</td><td>53</td><td>9</td><td>0</td><td>9</td><td>0</td></tr><tr><td>45</td><td>82432</td><td>1</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td>54</td><td>88663</td><td>23</td><td>2</td><td>0</td><td>1</td><td>1</td></tr><tr><td>241</td><td>198089</td><td>31</td><td>3</td><td>0</td><td>3</td><td>0</td></tr><tr><td rowspan="4">Rogassa</td><td>5</td><td>129395443</td><td>98</td><td>14</td><td>0</td><td>13</td><td>1</td></tr><tr><td>6</td><td>19034257</td><td>614</td><td>62</td><td>42</td><td>8</td><td>12</td></tr><tr><td>65</td><td>109454</td><td>36</td><td>5</td><td>0</td><td>5</td><td>0</td></tr><tr><td>68</td><td>125577517</td><td>2405</td><td>426</td><td>323</td><td>37</td><td>66</td></tr></tbody><tfoot><tr><td colspan="2">Total</td><td></td><td>4243</td><td></td><td></td><td></td><td>101</td></tr></tfoot></table>

<table><caption>**Tableau A2 : Statistiques de l'extraction de l'objet bâti par algorithme IsoData**</caption><caption>Extraction non supervisée IsoData multi bandes: objet bâti isolé</caption><thead><tr><th>Commune</th><th>îlot</th><th>Surface m2</th><th>Extraction brute</th><th>Après nettoyage</th><th>Confusion avec objet route</th><th>Confusion avec d'autres objets</th><th>Objets validés</th></tr></thead><tbody><tr><td rowspan="10">El Bayadh</td><td>1</td><td>12511088</td><td>744</td><td>216</td><td>34</td><td>175</td><td>7</td></tr><tr><td>2</td><td>353347</td><td>88</td><td>21</td><td>0</td><td>18</td><td>3</td></tr><tr><td>3</td><td>445024</td><td>53</td><td>16</td><td>0</td><td>11</td><td>5</td></tr><tr><td>4</td><td>35398</td><td>2</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td>12</td><td>148160</td><td>69</td><td>16</td><td>0</td><td>11</td><td>5</td></tr><tr><td>18</td><td>142945</td><td>80</td><td>23</td><td>0</td><td>23</td><td>0</td></tr><tr><td>37</td><td>538384</td><td>37</td><td>8</td><td>3</td><td>2</td><td>3</td></tr><tr><td>48</td><td>159236</td><td>5</td><td>2</td><td>0</td><td>1</td><td>1</td></tr><tr><td>54</td><td>88663</td><td>10</td><td>5</td><td>0</td><td>1</td><td>4</td></tr><tr><td rowspan="3">Rogassa</td><td>6</td><td>19034257</td><td>569</td><td>129</td><td>40</td><td>71</td><td>18</td></tr><tr><td>5</td><td>129395443</td><td>68</td><td>18</td><td>0</td><td>16</td><td>2</td></tr><tr><td>68</td><td>125577517</td><td>2650</td><td>652</td><td>76</td><td>541</td><td>35</td></tr></tbody><tfoot><tr><td colspan="2">Total</td><td></td><td>4375</td><td></td><td></td><td></td><td>84</td></tr></tfoot></table>
<table><caption>**Tableau A3 :** Statistiques de l'extraction de l'objet bâti par algorithme K-Means</caption><caption>Extraction non supervisée K-Means multi bandes: objet bâti isolé</caption><thead><tr><th>Commune</th><th>Ilot</th><th>Surface m2</th><th>Extraction brute</th><th>Après nettoyage</th><th>Confusion avec objet route</th><th>Confusion avec d'autres objets</th><th>Objets validés</th></tr></thead><tbody><tr><td rowspan="10">El Bayadh</td><td>2</td><td>353347</td><td>37</td><td>8</td><td>0</td><td>8</td><td>0</td></tr><tr><td>18</td><td>142945</td><td>55</td><td>8</td><td>0</td><td>8</td><td>0</td></tr><tr><td>4</td><td>35398</td><td>12</td><td>2</td><td>0</td><td>0</td><td>2</td></tr><tr><td>5</td><td>446761</td><td>5</td><td>1</td><td>0</td><td>1</td><td>0</td></tr><tr><td>7</td><td>2865137</td><td>405</td><td>97</td><td>8</td><td>61</td><td>28</td></tr><tr><td>12</td><td>148160</td><td>78</td><td>18</td><td>11</td><td>1</td><td>6</td></tr><tr><td>37</td><td>538384</td><td>119</td><td>7</td><td>1</td><td>5</td><td>1</td></tr><tr><td>39</td><td>1523879</td><td>44</td><td>8</td><td>0</td><td>8</td><td>0</td></tr><tr><td>40</td><td>209450</td><td>22</td><td>5</td><td>0</td><td>5</td><td>0</td></tr><tr><td>54</td><td>88663</td><td>34</td><td>4</td><td>0</td><td>4</td><td>0</td></tr><tr><td>241</td><td>198089</td><td>59</td><td>7</td><td>0</td><td>7</td><td>0</td></tr><tr><td rowspan="3">Rogassa</td><td>1</td><td>12511088</td><td>1078</td><td>185</td><td>0</td><td>164</td><td>21</td></tr><tr><td>3</td><td>445024</td><td>35</td><td>6</td><td>0</td><td>4</td><td>2</td></tr><tr><td>6</td><td>19034257</td><td>331</td><td>47</td><td>43</td><td>2</td><td>2</td></tr><tr><td rowspan="2">Rogassa</td><td>65</td><td>109454</td><td>23</td><td>5</td><td>0</td><td>4</td><td>1</td></tr><tr><td>68</td><td>125577517</td><td>1487</td><td>74</td><td>12</td><td>53</td><td>9</td></tr></tbody><tfoot><tr><td>Total</td><td></td><td></td><td>3824</td><td></td><td></td><td></td><td>72</td></tr></tfoot></table>

<table><caption>**Tableau A4 :** Statistiques de l'extraction de l'objet bâti par classification supervisée (maximum de vraisemblance)</caption><caption>Extraction supervisée Multi bandes : objet bâti isolé</caption><thead><tr><th>Commune</th><th>Îlot</th><th>Surface m2</th><th>Extraction brute</th><th>Après nettoyage</th><th>Confusion avec objet route</th><th>Confusion avec d'autres objets</th><th>Objets validés</th></tr></thead><tbody><tr><td rowspan="10">El Bayadh</td><td>1</td><td>12511088</td><td>3496</td><td>700</td><td>75</td><td>624</td><td>1</td></tr><tr><td>2</td><td>353347</td><td>219</td><td>42</td><td>0</td><td>38</td><td>4</td></tr><tr><td>18</td><td>142945</td><td>55</td><td>8</td><td>0</td><td>8</td><td>0</td></tr><tr><td>4</td><td>35398</td><td>33</td><td>5</td><td>0</td><td>3</td><td>2</td></tr><tr><td>5</td><td>446761</td><td>155</td><td>29</td><td>0</td><td>21</td><td>8</td></tr><tr><td>6</td><td>446761</td><td>47</td><td>11</td><td>0</td><td>8</td><td>3</td></tr><tr><td>7</td><td>2865137</td><td>588</td><td>120</td><td>0</td><td>106</td><td>14</td></tr><tr><td>12</td><td>142945</td><td>141</td><td>38</td><td>0</td><td>32</td><td>6</td></tr><tr><td>37</td><td>538384</td><td>122</td><td>119</td><td>0</td><td>116</td><td>3</td></tr><tr><td>54</td><td>88663</td><td>58</td><td>20</td><td>0</td><td>15</td><td>5</td></tr><tr><td>241</td><td>198089</td><td>98</td><td>19</td><td>0</td><td>18</td><td>1</td></tr><tr><td rowspan="3">Rogassa</td><td>3</td><td>445024</td><td>180</td><td>33</td><td>0</td><td>30</td><td>3</td></tr><tr><td>48</td><td>159236</td><td>34</td><td>6</td><td>0</td><td>2</td><td>4</td></tr><tr><td>6</td><td>19034257</td><td>7803</td><td>1532</td><td>253</td><td>1251</td><td>28</td></tr><tr><td rowspan="2">Rogassa</td><td>5</td><td>129395443</td><td>358</td><td>100</td><td>6</td><td>89</td><td>5</td></tr><tr><td>68</td><td>125577517</td><td>14138</td><td>2723</td><td>393</td><td>2277</td><td>53</td></tr></tbody><tfoot><tr><td>Total</td><td></td><td></td><td>27525</td><td></td><td></td><td></td><td>140</td></tr></tfoot></table>
<table><caption>Tableau A5 : statistiques de l'extraction de l'objet bâti par approche orientée objet</caption><caption>Extraction du bâti par règles d'apprentissage (sans nettoyage)</caption><thead><tr><th>Commune</th><th>îlot</th><th>Superficie</th><th>Extraction brute</th><th>Avec objet route</th><th>Avec d'autres objets</th><th>Objets validés</th></tr></thead><tbody><tr><td rowspan="6">El Bayadh</td><td>1</td><td>12511087</td><td>42</td><td>10</td><td>24</td><td>8</td></tr><tr><td>2</td><td>353346</td><td>2</td><td>1</td><td>0</td><td>1</td></tr><tr><td>3</td><td>445023</td><td>4</td><td>0</td><td>1</td><td>3</td></tr><tr><td>12</td><td>148159</td><td>4</td><td>0</td><td>1</td><td>3</td></tr><tr><td>18</td><td>142945</td><td>2</td><td>0</td><td>0</td><td>2</td></tr><tr><td>45</td><td>82431</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td>5</td><td>446760</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td rowspan="4">Rogassa</td><td>6</td><td>19034257</td><td>29</td><td>5</td><td>9</td><td>15</td></tr><tr><td>5</td><td>129395443</td><td>3</td><td>0</td><td>2</td><td>1</td></tr><tr><td>68</td><td>125577517</td><td>155</td><td>26</td><td>85</td><td>44</td></tr><tr><td>64</td><td>72575</td><td>2</td><td>0</td><td>1</td><td>1</td></tr><tr><td colspan="2">Total</td><td></td><td>245</td><td></td><td></td><td>80</td></tr></tbody></table>

<table><caption>Tableau A6 : récapitulatif des incertitudes</caption><thead><tr><th>Classification</th><th></th><th>Extraction brute</th><th>Information certaine en%</th><th>Information incertaine en%</th></tr></thead><tbody><tr><td>Segmentation mono bandes</td><td></td><td>637</td><td>15.85</td><td>84.14</td></tr><tr><td rowspan="2">Non supervisée</td><td>ISODATA</td><td>1107</td><td>14.93</td><td>92.41</td></tr><tr><td>K-MEANS</td><td>482</td><td>7.58</td><td>85.06</td></tr><tr><td>Supervisée Maximum de vraisemblance</td><td></td><td>5505</td><td>2.543</td><td>97.45</td></tr><tr><td>Orientée objet</td><td></td><td>245</td><td>33.06</td><td>66.94</td></tr></tbody></table>

NB : Suite au taux d'incertitude de 60% (extraction brute) de l'approche orientée objet pour le cas du bâti isolé, une confirmation de l'identification à partir des images Google-Earth est nécessaire avant d'amorcer les contrôles sur terrain.