Suite

Mosaic to new raster : même nombre de bandes émises

Mosaic to new raster : même nombre de bandes émises


J'essaie de mosaïquer 32 fichiers raster (valeur minimale 0, valeur maximale 13) dans un fichier raster afin que je puisse utiliser ultérieurement ce nouveau raster pour extraire plusieurs valeurs en points. Cependant, lorsque j'essaie de mosaïquer les rasters, cette erreur apparaît :

<>.

Je suis nouveau sur ArcGIS, donc je ne sais pas ce que cela signifie, d'autant plus que je ne vois pas de jeu de données cible ? De plus, je ne sais pas quelle devrait être la valeur de la bande? Lorsque j'ai enquêté sur la source de données, il a dit:

Les images composites de référence sont des observations médianes d'un ensemble d'observations de saison de croissance évaluées par la qualité dans quatre bandes spectrales, en particulier les bandes Landsat 3, 4, 5 et 7

Par conséquent, logiquement, je supposerais que la valeur de la bande est de quatre ? Je ne comprends tout simplement pas pourquoi cela ne fonctionne pas. J'utilise ArcMap 10.2.


Voici quelques informations générales sur les bandes raster. Pensez-y comme ceci : une image RVB aurait 3 bandes, une pour le rouge, une pour le vert et une pour le bleu. Chaque cellule raster aurait 3 valeurs, une pour chaque bande. L'imagerie Landsat en a 4, car elles incluent également une bande supplémentaire pour les valeurs infrarouges.

Dans votre cas, il n'y a qu'un seul canal car chaque cellule n'a qu'une valeur, une valeur de 0 à 13. Pour le confirmer, ouvrez les propriétés du raster et regardez dans l'onglet source. Consultez également l'onglet symbologie.

Lorsque le message d'erreur indique "raster cible", il fait simplement référence à la sortie que vous essayez de créer. Vous mettez un tas de rasters à bande unique et vous vous attendez à obtenir un raster à 4 bandes.


Il s'est avéré que je ne pouvais pas reprojeter à la volée car je mosaïque tous les fichiers raster. Je devais d'abord les mosaïquer et une fois cela fait, alors seulement je pouvais le re-projeter. Intéressant.


Mosaïque à nouveau raster : problème du même nombre de bandes - Systèmes d'Information Géographique

Une grande variété d'images satellitaires et de photographies aériennes est disponible pour une utilisation dans les systèmes d'information géographique (SIG). Bien que ces produits soient essentiellement des graphiques raster, ils sont fondamentalement différents dans leur utilisation au sein d'un SIG. L'imagerie satellitaire et la photographie aérienne fournissent des informations contextuelles importantes pour un SIG et sont souvent utilisées pour effectuer une numérisation tête haute (chapitre 5 « Geospatial Data Management », Section 5.1.4 « Secondary Data Capture ») où les caractéristiques de l'image sont converties en vecteurs ensembles de données.


Ajout de données raster à une mosaïque

Lorsque vous ajoutez des données raster à une mosaïque, un pointeur vers l'emplacement des données est stocké dans la mosaïque. Les données raster réelles ne sont pas transférées ou stockées dans la mosaïque.

Les données raster sont ajoutées à une mosaïque en spécifiant un type de raster. Le type raster identifie les métadonnées, telles que le géoréférencement, la date d'acquisition, le type de capteur et les longueurs d'onde de bande, ainsi qu'un format raster. Un format raster définit la manière dont les pixels sont stockés, comme le nombre de lignes et de colonnes, le nombre de canaux, les valeurs de pixels réelles et d'autres paramètres spécifiques au format raster. Cependant, en ajoutant des données raster selon un type de raster, les métadonnées appropriées sont lues et utilisées pour définir tout traitement devant être appliqué. Par exemple, lors de l'ajout d'une scène QuickBird Standard, une scène peut être définie par un fichier .imd. Le fichier .imd contient des informations de métadonnées sur le jeu de données raster et peut pointer vers un ou plusieurs fichiers .tif. Pour ajouter correctement ces données, utilisez le type raster QuickBird car il recherche cette combinaison de types de fichiers. Si vous ajoutez les données raster à l'aide du type raster Raster Dataset, seuls les fichiers .tif seront reconnus et ajoutés, ils seront ajoutés en tant que fichiers TIFF, et toutes les informations de métadonnées pouvant affecter les fonctions nécessaires ou le géoréférencement seront manquantes.

Le type raster Jeu de données raster fait référence à tout format raster pris en charge par ArcGIS. Les autres types de raster sont spécifiques à un produit (capteur, fournisseur de données ou fournisseur) ou aux informations de métadonnées associées. Vous avez la possibilité de modifier n'importe quel type de raster lors de l'ajout de données. Par exemple, vous souhaiterez peut-être définir une combinaison de bandes spécifique, un algorithme de netteté panoramique ou ajouter un filtre lorsque les données sont ajoutées, plutôt que plus tard. Vous pouvez modifier n'importe quel type de raster en changeant ses propriétés par défaut, telles que la combinaison de canaux, ou en définissant le modèle d'altitude à utiliser lors de l'orthorectification des données. Vous pouvez même ajouter ou modifier les chaînes de fonctions. Chaque fois que vous modifiez le type de raster lors de l'ajout de rasters à une mosaïque, vous pouvez l'enregistrer dans un nouveau fichier .art afin de pouvoir charger des données supplémentaires à un autre moment en utilisant les mêmes paramètres modifiés. En option, vous pouvez l'enregistrer en tant que .art.xml si vous souhaitez le lire ou le modifier directement.

Les données raster peuvent être ajoutées à partir de nombreuses sources différentes, y compris un fichier raster unique, un catalogue d'images, une table, une mosaïque, un ou plusieurs dossiers, ArcSDE, une couverture WCS, un service d'imagerie, un service de carte mis en cache ou un service WMS .

Si vous comptez utiliser l'outil Ajouter des rasters à la mosaïque plusieurs fois pour ajouter des données à votre mosaïque, vous pouvez décocher Mettre à jour les plages de taille de cellule et Mettre à jour la limite jusqu'à ce que vous ayez ajouté toutes les données raster que vous souhaitez inclure. dans la mosaïque.

Si vous ajoutez des données raster à une mosaïque dont les aperçus sont générés, vous devrez mettre à jour les aperçus ou en ajouter de nouveaux. Cela peut être fait en même temps que les données raster sont ajoutées en cochant l'option Mettre à jour les aperçus dans la boîte de dialogue de l'outil Ajouter des rasters à la mosaïque. Ou vous pouvez le faire plus tard à l'aide de l'outil Build Overviews.

  1. Ouvrez l'outil Ajouter des rasters à la mosaïque.
    • Pour ce faire, dans ArcCatalog ou dans la fenêtre Catalogue, cliquez avec le bouton droit sur une mosaïque et cliquez sur Ajouter des rasters .
    • Faites glisser et déposez un dossier ou un jeu de données sur une mosaïque.
    • Naviguez dans les boîtes à outils du système jusqu'à l'outil Ajouter des rasters à un jeu de données en mosaïque.
    • Utilisez la fenêtre de recherche pour localiser l'outil Ajouter des rasters à un jeu de données en mosaïque.
  2. Vérifiez ou saisissez le chemin et le nom de la mosaïque à laquelle les données seront ajoutées.
  3. Cliquez sur la flèche déroulante Type de raster et sélectionnez le type de raster ou cliquez sur le bouton Parcourir et accédez à un fichier de type raster (.art ou .art.xml).
  4. Cliquez éventuellement sur le bouton Propriétés du type raster .

Dans cette boîte de dialogue, vous pouvez modifier divers paramètres, tels que le modèle de traitement, la chaîne de fonctions, la combinaison de bandes et le type d'étirement.

Il existe de nombreuses options sur cet outil. Pour plus d'informations, consultez Outil Ajouter des rasters à la mosaïque et Création d'une mosaïque.


IllegalArgumentException : les nombres de canaux raster source et les composants de l'espace colorimétrique source ne correspondent pas. Pour une exception d'image couleur

La réponse ci-dessus que quelqu'un a suggérée convertit mon image colorée en une image en noir et blanc. Donc ce n'est pas approprié pour ma question.

voici le code et ci-dessous est l'image. Téléchargez l'image et enregistrez-la sur votre PC. Et essayez d'exécuter le code ci-dessus avec une valeur correcte de chemin, cela lèvera une exception dans le sujet

Simplement si quelqu'un peut obtenir un objet java.awt.image.BufferedImage à partir de l'image étant donné que cela suffit (ne devrait pas convertir l'image en échelle de gris).

Tu es un génie si tu peux répondre à ça :D. Aide SVP.


Mosaïque à nouveau raster : problème du même nombre de bandes - Systèmes d'Information Géographique

J'ai besoin de convertir un DEM en contours.

Mon problème :
Le DEM fait plus de 30 Go. J'ai essayé de découper le DEM de deux manières différentes en fonction d'une couche de polygones qui supprimerait 90% du DEM. Cependant, les deux fois, le DEM écrêté résultant est toujours proche ou supérieur à 30 Go et le traitement prend des jours.

Ce que j'ai essayé :
1. Découpage du DEM via l'analyse d'images. J'ai sélectionné tous les polygones via la table attributaire de la couche de polygones. Sélectionnez ensuite le DEM répertorié dans l'analyse d'image et utilisez l'outil Clip dans la section de traitement. La zone découpée apparaît rapidement. Lorsque je vais enregistrer/exporter le DEM tronqué, le processus est lent et le résultat est que le DEM tronqué est toujours (ou proche) de la taille du fichier d'origine, même si le DEM est tronqué. Lorsque j'essaie de créer les contours de ce DEM écrêté, j'obtiens le message par défaut que le fichier dépasse la taille limite de 2 Go.

2. L'utilisation de l'outil ArcToolbox Clip Raster donne le même résultat.

3. J'ai également utilisé ETSurface et j'ai obtenu le même résultat (mais un temps de traitement beaucoup plus rapide).

Mes questions:
1. Pourquoi le tampon écrêté a-t-il presque la même taille que l'original ?
2. Une durée de fonctionnement de 2-3 jours est-elle normale ? 3. Je dois me retrouver avec des contours basés sur la zone du polygone.
4. Y a-t-il une autre façon que je peux essayer de faire cela?


Format de données raster pour cet atelier

Les données raster peuvent se présenter sous de nombreux formats différents. Pour cet atelier, nous utiliserons le format GeoTIFF qui porte l'extension .tif . Un fichier .tif stocke les métadonnées ou les attributs du fichier sous forme de balises tif intégrées. Par exemple, votre appareil photo peut stocker une étiquette décrivant la marque et le modèle de l'appareil photo ou la date à laquelle la photo a été prise lorsqu'il enregistre un .tif . Un GeoTIFF est un format d'image .tif standard avec des informations spatiales supplémentaires (géoréférencement) intégrées dans le fichier sous forme de balises. Ces balises doivent inclure les métadonnées raster suivantes :

  1. Le degré
  2. Résolution
  3. Système de référence de coordonnées (CRS) - nous présenterons ce concept dans un épisode ultérieur
  4. Valeurs qui représentent des données manquantes ( NoDataValue ) - nous présenterons ce concept dans une leçon ultérieure.

Nous discuterons de ces attributs plus en détail dans une leçon ultérieure. Dans cette leçon, nous apprendrons également à utiliser R pour extraire les attributs raster d'un fichier GeoTIFF.


Définition Source de la définition

Aperçu des entités et des attributs MELCD est un raster dont les valeurs de pixel représentent la classification des zones en 23 classifications : route/piste aménagée à haute intensité aménagée à moyenne intensité aménagée à faible intensité aménagée espace ouvert terre cultivée champ de bleuets pâturage/foin prairie/forêt herbacée à feuilles caduques forêt à feuilles persistantes mixte forêt broussailles arbustives forêt humide zone humide rivage non consolidé terre nue eau libre coupe à blanc coupe partielle légère (après 1995) coupe partielle lourde (après 1995) régénération de la forêt (après 1995) alpine/toundra.


Gdal_calc.py¶

Calculatrice raster de ligne de commande avec syntaxe numpy. Utilisez n'importe quelle arithmétique de base prise en charge par les tableaux numpy tels que + , - , * et avec des opérateurs logiques tels que > . Notez que tous les fichiers doivent avoir les mêmes dimensions (sauf si l'option d'étendue est utilisée), mais aucune vérification de projection n'est effectuée (sauf si l'option projectionCheck est utilisée).

Afficher ce message d'aide et quitter

Calcul en syntaxe numpy en utilisant + , - , / , * , ou toute fonction de tableau numpy (c'est-à-dire log10() ). Plusieurs options --calc peuvent être répertoriées pour produire un fichier multibande (GDAL >= 3.2).

Fichier raster gdal d'entrée, vous pouvez utiliser n'importe quelle lettre (a-z, A-Z). (minuscule supporté depuis GDAL 3.3)

Une lettre peut être répétée (GDAL >= 3.3). L'effet sera de créer un tableau numpy 3-dim. Dans un tel cas, la formule de calcul doit utiliser cette entrée comme un tableau à 3 dimensions et doit renvoyer un tableau 2D (voir exemples ci-dessous). Si le calcul ne renvoie pas un tableau 2D, une erreur serait générée.

Numéro de bande raster pour le fichier A (par défaut 1).

Fichier de sortie à générer ou à remplir.

Sortie NoDataValue (valeur spécifique au type de données par défaut). Pour indiquer de ne pas définir de NoDataValue, utilisez –NoDataValue=none (GDAL >= 3.3)

Utilisation de l'API Python : aucune valeur indiquera la valeur par défaut spécifique au type de données. La valeur 'aucun' indiquera que vous n'avez pas défini de NoDataValue.

Ignore les bandes d'entrée NoDataValue. Par défaut, les bandes d'entrée NoDataValue ne participent pas au calcul. En définissant ce paramètre, aucun traitement spécial ne sera effectué sur l'entrée NoDataValue. et ils participeront au calcul comme toute autre valeur. La sortie n'aura pas de valeur NoDataValue définie, sauf si vous avez explicitement spécifié une valeur spécifique en définissant –NoDataValue=<value>.

Le type de données de sortie doit être l'un des [ Int32 , Int16 , Float64 , UInt16 , Byte , UInt32 , Float32 ].

Malgré le type de données défini à l'aide de --type , lors d'opérations arithmétiques intermédiaires utilisant des opérandes du même type, le résultat de l'opération respectera le type de données d'origine. Cela peut conduire à des résultats inattendus dans le résultat final.

Format GDAL pour le fichier de sortie.

Permet de spécifier un nom de fichier d'une table de couleurs (ou d'un objet ColorTable) (avec interprétation de l'index de palette) à utiliser pour le raster en sortie. Formats pris en charge : txt (c'est-à-dire comme gdaldem, mais les noms de couleurs ne sont pas pris en charge), qlr, qml (c'est-à-dire exportés depuis QGIS)

cette option détermine comment gérer les rasters avec différentes étendues. cette option est mutuellement exclusive avec l'option projwin, qui est utilisée pour fournir une extension personnalisée. pour toutes les options ci-dessous, la taille en pixels (résolution) et le SRS (Spatial Reference System) de tous les rasters en entrée doivent être les mêmes. ignorer (par défaut) - seules les dimensions des rasters sont comparées. si les dimensions ne concordent pas, l'opération échouera. fail - les dimensions et l'étendue (limites) des rasters doivent concorder, sinon l'opération échouera. union - l'étendue (limites) de la sortie sera le rectangle minimal qui contient toutes les étendues d'entrée. intersection - l'étendue (limites) de la sortie sera le rectangle maximal contenu dans toutes les étendues d'entrée.

cette option fournit une étendue personnalisée pour la sortie, elle est mutuellement exclusive avec l'option d'étendue.

Par défaut, aucune vérification de projection ne sera effectuée. En définissant cette option, si la projection n'est pas la même pour toutes les bandes, l'opération échouera.

Passe une option de création au pilote de format de sortie. Plusieurs options peuvent être répertoriées. Voir la documentation spécifique au format pour les options de création légale pour chaque format.

Traiter toutes les bandes d'un raster donné (a-z, A-Z). Nécessite un seul calcul pour toutes les bandes.

Écraser le fichier de sortie s'il existe déjà.

Imprimer les informations de débogage.

Supprimer les messages de progression.


Mosaïque à nouveau raster : problème du même nombre de bandes - Systèmes d'Information Géographique

L'API raster ArcSDE fournit un mécanisme permettant de charger des données raster à partir d'une application sur le serveur ArcSDE et de les récupérer. L'API est une interface de bas niveau pour l'accès de base aux données raster en lecture-écriture au serveur ArcSDE. Ainsi, il permet aux organisations de stocker leurs données raster dans l'un des SGBD relationnels commerciaux pris en charge (IBM DB2, Informix, Microsoft SQL Server et Oracle) ou open source (PostgreSQL) pour un accès rapide, en ligne et multi-utilisateurs aux données raster continues.

Les valeurs stockées dans un raster représentent des données discrètes ou continues. Les données raster discrètes représentent des entités géographiques qui ont des limites définissables. Les entités géographiques discrètes comprennent les bâtiments, les puits et les parcelles de terrain. Les données raster continues représentent des informations géographiques qui peuvent varier à chaque emplacement, comme l'altitude, la température ou la force du signal radio. Les phénomènes géographiques continus n'ont pas de frontières distinctes.

Avec les données raster, chaque pixel est indépendant des autres pixels, enregistrant des informations uniquement sur la zone qu'il couvre. Par exemple, un pixel ne sait pas que le pixel à côté peut faire partie du même arbre ou bâtiment.

Un raster est parfois aussi appelé ligne de balayage. Une image matricielle est donc une image constituée de telles lignes de balayage.


Préparation des données spatiales

En général, quatre étapes doivent être mises en œuvre pour délimiter les zones de disponibilité des aquifères en utilisant les approches FR et SE. Ces étapes sont les suivantes : (a) Collecter des données sur les emplacements des forages et des données arbitraires divisées en deux ensembles : formation et test selon des critères spécifiques tels que (70/30) ou (80/20). L'ensemble de données de formation est uniquement utilisé pour étudier la relation statistique entre les emplacements des forages et les eaux souterraines influençant le facteur d'occurrence. L'ensemble de test est utilisé pour valider les résultats et montrer la capacité du modèle à prédire les zones de disponibilité (test à l'aveugle). (b) Construire la base de données géospatiale. Dans cette étape de l'analyse, les couches thématiques raster probantes des facteurs d'occurrence des eaux souterraines sont préparées à l'aide de différentes ressources telles que l'enquête conventionnelle, l'enquête sur le terrain et la RS. Toutes les couches thématiques doivent être converties au format raster pour être utilisées dans une analyse ultérieure. (c) La relation entre les emplacements des forages (pour l'ensemble de données de formation) et les facteurs d'occurrence des eaux souterraines est étudiée à l'aide de FR et de l'indice d'entropie dans la troisième étape de l'analyse. Le calcul du rapport de vraisemblance pour chaque classe de chaque facteur est effectué et le poids approprié pour chaque facteur utilisé est également calculé. Le GAI est ensuite calculé et classé en différentes classes en utilisant un schéma de classification approprié tel que Natural Break, Geometric, etc. en fonction des conditions de la zone d'étude et de l'expérience personnelle. (d) La quatrième étape implique la validation des résultats et la comparaison de l'efficacité du modèle dans la prédiction des zones de disponibilité des eaux souterraines. En général, la validation des résultats est effectuée pour la formation (appelée taux de réussite) et les tests (taux de prédiction) en utilisant la technique bien connue des caractéristiques d'exploitation relatives (ROC), ou directement en comparant les emplacements des forages avec les zones prospectives d'eau souterraine prévues. Parfois, si le GAI est estimé à l'aide de différentes techniques, les capacités de ces techniques sont également comparées et le modèle sélectionné est le meilleur avec une erreur de prédiction minimale. Un organigramme est présenté sur la figure 4 pour clarifier la procédure mentionnée précédemment.

Organigramme de la délimitation des zones de disponibilité des eaux souterraines

Inventaire des emplacements de forage

Les données de forage ont été obtenues auprès de la Commission générale des eaux souterraines/Ministère des ressources en eau, Irak. L'enregistrement des données impliquait de nombreuses données pertinentes telles que l'emplacement géographique (UTM), le débit du forage (l/s), la profondeur du forage (m), le type d'aquifère (confiné, non confiné) et l'analyse physico-chimique des constituants des eaux souterraines. Au total, il existe 211 forages productifs dans la zone d'étude. À partir de ces forages, seuls les forages à débit relativement élevé (8 l/s) (137) ont été cartographiés et utilisés pour construire les modèles de potentiel des eaux souterraines. Le débit a été choisi après examen de la littérature et obtention d'avis d'experts sur les eaux souterraines (Jabar Al-Syadi, Expert, Commission générale des eaux souterraines/Ministère des ressources en eau, Irak, communication personnelle). Sur les 211 forages, 95 (70 %) ont été sélectionnés au hasard comme données de formation et les 43 autres (30 %) ont été conservés à des fins de validation. Le logiciel commercial statistique MINITAB v.16 a été utilisé pour le fractionnement des données.

Sélection des facteurs d'occurrence des eaux souterraines

Dans cette étude, dix facteurs d'occurrence des eaux souterraines ont été sélectionnés sur la base de l'opinion d'experts et d'une revue de la littérature. Ces facteurs étaient l'altitude (m), la pente (°), la courbure, l'aspect, le TWI, le SPI, la géologie, le sol, l'utilisation des terres/couverture des terres (LULC) et la distance aux failles. L'altitude est un facteur important pour l'occurrence des eaux souterraines parce que les conditions météorologiques et climatiques varient considérablement à différentes altitudes, ce qui a entraîné des différences dans le sol et la végétation (Aniya 1985). La pente est une montée ou une descente de la surface du sol. C'est un facteur important pour les études de cartographie de la disponibilité des eaux souterraines car il contrôle l'accumulation d'eau dans une zone et améliore ainsi la recharge des eaux souterraines (Ozdemir 2011a). La courbure est la dérivée seconde d'une surface, ou la pente de la pente (Kimerling et al. 2011). Il représente la morphologie de la topographie. Il existe trois types de courbure différents : total, profil et plan. La courbure du profil est parallèle à la direction de la pente maximale et affecte principalement l'accélération de la décélération de l'écoulement à travers la surface. Une valeur négative de la courbure du profil indique que la surface est convexe vers le haut, et une courbure négative implique que la surface est concave vers le haut, tandis que la valeur zéro indique que la surface est linéaire (Oh et Lee 2010). La courbure du plan est perpendiculaire à la direction de la pente maximale et affecte principalement la convergence et la divergence de l'écoulement à travers la surface. Une valeur positive indique que la surface est convexe latéralement, une valeur négative indique que la surface est concave latéralement, une valeur de zéro implique que la surface est linéaire (Kimerling et al. 2011). La combinaison du profil et de la courbure en plan est appelée courbure totale. La prise en compte des courbures en plan et en profil permet de comprendre plus précisément l'écoulement sur une surface. La direction de l'aspect ou de la pente identifie la direction de la pente descendante du taux maximum de changement de valeur de chaque cellule d'un raster à ses voisines (Burrough et McDonnell 1998). Les valeurs d'aspect indiquent la direction de la boussole à laquelle la surface fait face à cet endroit. Il est mesuré dans le sens des aiguilles d'une montre en degrés de 0 (nord) à 360 (à nouveau plein nord), en faisant un cercle complet. En raison du fait que les zones plates n'ont pas de direction de pente descendante, une valeur de -1 est donnée. L'aspect affecte fortement les processus hydrologiques via l'évapotranspiration, la direction des précipitations frontales, et affecte donc le processus d'altération et le développement de la végétation et des racines, en particulier dans les environnements plus secs, il est donc pris en compte dans cette étude. Les indices topographiques tels que TWI et SPI ont un rôle vital dans la variation spatiale des conditions hydrologiques telles que l'humidité du sol, l'écoulement des eaux souterraines et la stabilité des pentes (Lee et Kim 2011). Ces indices topographiques ont été utilisés pour décrire les modèles spatiaux d'humidité du sol (Moore et al. 1991). La géologie est un facteur important dans la cartographie du potentiel des eaux souterraines (Ozdemir 2011a, b Lee et al. 2012 Manap et al. 2011 Nampak et al. 2014). La géologie influence l'occurrence des eaux souterraines car la variation lithologique conduit souvent à une différence de porosité et de conductivité hydraulique des roches et des sols. Le sol fait référence à la partie supérieure de la zone non saturée caractérisée par des activités biologiques importantes. Le sol a un impact sur la quantité de recharge, qui peut s'infiltrer dans les eaux souterraines, et donc augmenter le stockage des eaux souterraines d'un aquifère. L'autre facteur considéré, le LULC, définit l'état biologique de la surface de la terre et comment les gens utilisent la terre et l'activité socio-économique. Elle est considérée comme un facteur de délimitation de la disponibilité des eaux souterraines par de nombreux chercheurs (Manap et al. 2011 Nampak et al. 2014 Gumma et Pavelic 2013 Moghaddam et al. 2013). La variation des catégories LULC contribue à la variation des conditions du sol et, par conséquent, à l'occurrence des eaux souterraines. Le réglage de la structure contrôle l'apparition et le mouvement des eaux souterraines. La plupart des roches possèdent des fractures et d'autres discontinuités qui facilitent le stockage et le mouvement des fluides à travers elles. Certaines discontinuités, par ex. failles et digues, peuvent également agir comme barrières à l'eau (Singhal et Gupta 1999). Les deux sont pris en considération dans cette étude en tant que principaux facteurs contribuant à la disponibilité des eaux souterraines dans la zone d'étude.

Préparation des couches raster thématiques

Tous les facteurs d'occurrence des eaux souterraines mentionnés ci-dessus ont été préparés au format raster avec une taille de cellule de 30 × 30 m à l'aide du logiciel commercial ArcGIS 10.2 et de ses extensions Spatial Analyst, Geostastical Analyst, Image Analyst, Arc Hydro et ArcTool box. Pour la classification des valeurs continues des couches raster influentes, différents schémas de classification tels que Jenks, equal et manual ont été utilisés. Le processus de classification des données combine les données brutes dans des classes prédéfinies, ou des bacs. Ces classes peuvent être représentées sur une carte par des symboles uniques ou, dans le cas de cartes thématiques, par une couleur ou une teinte unique. Le schéma de classification de Jenks (également appelé méthode de classification des ruptures naturelles) est une méthode de regroupement de données conçue pour déterminer la meilleure disposition des valeurs en différentes classes. La méthode cherche à réduire la variance au sein des classes et à maximiser la variance entre les classes (Jenks et al. 1967). Essentiellement, la méthode de Jenks minimise les variances au sein de la classe (les rend aussi similaires que possible) et maximise la variance entre les groupes (rend les classes de données aussi différentes que possible). L'avantage de la classification Natural Breaks (Jenks) est qu'elle identifie des classes réelles au sein des données. Ceci est utile car il crée des cartes de couches thématiques qui ont des représentations précises des tendances dans les données. La sélection de ce schéma de classification est basée sur des revues de la littérature et l'expérience de l'auteur de la zone d'étude et de ses conditions (Al-Abadi 2015b).

Pour créer des facteurs topographiques, c'est-à-dire l'élévation, l'angle de pente, la courbure, l'aspect, le TWI et le SPI, le modèle numérique mondial d'élévation du radiomètre spatial avancé d'émission et de réflexion (ASTER-GDEM) (http://gdem.ersdac.jspacesystems.or .jp/search.jsp) a été utilisé. L'ASTER-GDEM a été développé par le ministère de l'Économie du Japon et la National Aeronautics and Space Administration (NASA) des États-Unis. La résolution spatiale de la tuile ASTER-GDEM est d'environ 30 m. Six tuiles DEM brutes ont été téléchargées à partir de l'emplacement Web précédent et fusionnées pour créer un nouveau raster. Le nouveau raster est ensuite reprojeté dans le système de coordonnées projetées UTM (référence WGS, 38 N), et les cuvettes sont remplies pour créer une couche raster d'altitude remplie de la zone d'étude (Fig. 2a). Une zone sans puits, également appelée MNT sans dépression, est l'entrée préférable pour calculer la direction d'écoulement dans le bassin. Le processus cherche à remplir les puits d'un réseau DEM. Par conséquent, si des cellules avec une altitude plus élevée entourent une cellule du réseau DEM, l'eau est piégée et ne peut pas s'écouler. Les puits de remplissage modifient les valeurs d'altitude autour de la cellule pour éliminer ces problèmes. La couche raster d'altitude a été dérivée directement du DEM rempli et classée à l'aide d'un schéma de classification des ruptures en cinq classes (Fig. 2) à utiliser dans l'analyse ultérieure. L'angle de pente (°) de la zone d'étude a été dérivé du MNT rempli à l'aide de l'extension Spatial Analyst d'ArcMap et présenté sur la figure 5 après avoir été classé manuellement en cinq catégories : pente douce et plate <5°, pente moyenne (5°–15° ), pente modérée (15°–30°), pente raide (30°–50°) et pente très raide >50°, (Pourghasemi et al. 2013). La carte raster de courbure a également été dérivée directement du MNT rempli et le raster résultant a été classé manuellement en trois classes : convexe <0, plat 0 et concave >0 (Fig. 6). La carte d'aspect a également été préparée à partir du DEM de remplissage et classée en dix classes (Fig. 7) : Plat (−1), Nord (0–22,5) (337,5–360), Nord-est (22,5–67,5), Est (67,5–112,5 ), sud-est (112,5–157,5), sud (157,5–202,5), sud-ouest (202,5–247,5), ouest (247,5–292,5) et nord-ouest (292,5–337,5). Le TWI et le SPI sont définis mathématiquement comme (Moore et al. 1991) :

où, une est la zone locale non inclinée s'écoulant à travers un certain point par unité de longueur de contour et de tan β est la pente locale en degrés, et UNE s est la zone de chalandise spécifique. Pour calculer les facteurs TWI et SPI, la direction du flux et les couches d'accumulation du flux doivent d'abord être calculées, qui sont considérées comme les principales étapes de l'analyse du terrain et de la délimitation des bassins versants. La fonction de direction d'écoulement est calculée pour chaque cellule du DEM rempli dans la direction dans laquelle l'eau s'écoulerait à travers elle. La valeur de la cellule dans un raster de direction de flux est un nombre compris entre 1 et 128 qui représente une direction cardinale. La grille de direction du flux est utilisée comme entrée pour créer la couche d'accumulation de flux. La fonction d'accumulation de flux calcule pour chaque cellule du tableau DEM rempli, le nombre de cellules qui y coulent. L'accumulation de débit permet de déterminer la zone de drainage vers n'importe quel point spécifié dans un MNT. L'extension Arc Hydro d'ArcGIS a été utilisée pour dériver ces couches. La couche d'accumulation de flux a été utilisée pour dériver les couches raster de TWI et SPI en utilisant les équations. 10 et 11, puis classés en quatre classes pour les deux couches en utilisant le système de classification de Jenks (Figs. 8, 9).