Differences

This shows you the differences between two versions of the page.

--- geostats [2014/04/01 17:15]
glaroc [Comment créer un masque dans QGIS]
+++ geostats [2020/02/26 19:26] (current)
qcbs [Numériser un polygone définissant les limites de votre aire d'étude]
@@ Line 1: / Line 1: @@
 ====== Introduction à la géostatistique avec R ======
-Guillaume Larocque - Université Laval, 26 février 2014
+Guillaume Larocque ([email protected]), 27 février 2020. Université de Montréal.
 ===== Présentation =====
@@ Line 8: / Line 8: @@
 [[http://prezi.com/7dqatgpqh4xh/present/?auth_key=70wy1t3&follow=GuillaumeLarocque&kw=present-7dqatgpqh4xh&rc=ref-19631691|Présentation simultanée]]
-===== Ressources =====
+===== Livres =====
   * [[http://link.springer.com/book/10.1007/978-0-387-78171-6/page/1|Applied spatial data analysis with R, Bivand, Pebesma et Gómez-Rubio.]] - Bon survol des outils d'analyse spatiale dans R et un bon chapitre sur la géostatistique.
   * [[http://www.amazon.ca/Geostatistics-Natural-Resources-Evaluation-Goovaerts/dp/0195115384|Geostatistics for Natural Resources Evaluation, Pierre Goovaerts.]] - Très bon livre pour une connaissance approfondie de la géostatistique.
-  * [[http://www.gstat.org|Site web de gstat.]]
-  * [[http://cran.r-project.org/web/packages/gstat/|Package gstat]] - Package de gstat utilisé pour l'atelier.
-  * [[http://cran.r-project.org/web/packages/geoR/|Package geoR]] - Autre package pour la géostatistique en R.
   * [[http://cg.ensmp.fr/bibliotheque/public/MATHERON_Ouvrage_00161.pdf|Georges Matheron - La théorie des variables régionalisées et ses applications (PDF).]] - Text original décrivant les bases de la géostatistique.
   * {{:gisonr.pdf|Analyses spatiales sous R}} - Document préparé par Nicolas Casajus pour un cours offert à l'UQAR à l'hiver 2013.
+===== Articles intéressants =====
+  * [[http://www.esajournals.org/doi/abs/10.2307/2937096|Geostatistical Tools for Modeling and Interpreting Ecological Spatial Dependence]]
+  * [[http://link.springer.com/article/10.1023/A%3A1012734519752|GIS and geostatistics: Essential partners for spatial analysis]]
+  * [[http://www.ecoscience.ulaval.ca/en/paper/spatial-autocorrelation-and-statistical-tests-in-ecology-fait-parti-du-numero-thematique-sur-les-analyses-spatiales-en-ecologie-part-of-the-special-feature-on-spatial-analyses-in-ecology|Spatial autocorrelation and statistical tests in ecology]]
+  * [[http://www.esajournals.org/doi/pdf/10.1890/0012-9658(2003)084%5B1045:SCIPCI%5D2.0.CO%3B2|Spatial covariance in plant communities: intergrating ordination, geostatistics and variance testing]]
+  * [[http://link.springer.com/article/10.1007/s10651-008-0090-z|Coregionalization analysis with a drift for multi-scale assessment of spatial relationships between ecological variables]]
+===== Paquets R =====
+  * [[http://www.gstat.org|Site web de Gstat.]]
+  * [[http://cran.r-project.org/web/packages/gstat/|Gstat package]] - Gstat R package.
+  * [[http://cran.r-project.org/web/packages/geoR/|geoR package geoR]] - Autre package pour la géostatistique dans R.
 ===== Code R =====
-{{::atelier_csbq_geostats_fichiers.zip|Cliquez ici pour télécharger le fichier zip contenant les scripts et autres fichiers qui seront utilisés pour l'atelier.
+{{ ::atelier_geostats.zip | Cliquez ici pour télécharger le fichier zip contenant les scripts et autres fichiers qui seront utilisés pour l'atelier}}. Veuillez extraire ce fichier dans un dossier facile d'accès sur votre ordinateur.
-}}
-Veuillez extraire ce fichier dans un dossier facile d'accès sur votre ordinateur.
 ===== Comment créer un masque dans QGIS =====
 ==== Numériser un polygone définissant les limites de votre aire d'étude ====
-  - Vous pouvez numériser le périmètre de votre aire d'étude à partir d'une couche Google Satellite en suivant les étapes suivantes:
+Vous pouvez numériser le périmètre de votre aire d'étude à partir d'une couche Google Satellite en suivant les étapes suivantes:
-  -   - Activez le 'OpenLayers Plugin' dans les Extensions.
+  - Ouvrez QGIS.
-  -   - Ajoutez une couche 'Google Satellite' en utilisant l'extension OpenLayers (>Extensions>OpenLayers). Notez que le SCR du canevas est maintenant 'Google Mercator' (WGS 84 / Pseudo mercator). Déplacez cette couche en dessous du canevas en cliquant-glissant le nom de la couche vers le bas dans la table des couches à gauche.
+  - Activez l'extension QuickMapServices.
-  -
+  - Ajoutez une couche 'Google Satellite' en utilisant l'extension QuickMapServices (>web>QuickMapServices). Notez que le SCR du canevas est maintenant WGS 84 / Pseudo mercator. Déplacez cette couche en dessous du canevas en cliquant-glissant le nom de la couche vers le bas dans la table des couches à gauche. Si vous avez d'autres couches vecteur ou raster pouvant vous aider à définir les limites de la zone d'étude, ajoutez les au canevas an cliquant sur les icônes appropriés (ajouter couche vecteurs ou raster).
+  - Numérisez la zone définissant votre aire d'étude. Vous devez d'abord ajouter une nouvelle couche shapefile de type polygone sous Couche>Nouveau>Nouvelle couche shapefile. Spécifiez le système de référence approprié selon l'échelle de vos données et donnez à ce fichier shapefile un nom approprié. Il est toujours préférable de travailler avec un système de référence avec des unités en mètres comme UTM, MTM ou Lambert, et non pas en degrés comme latitude/longitude. Activez le mode édition (>Couche>Basculer en mode édition) et numérisez le polygone en cliquant sur l'icône {{::capture_du_2013-09-25_10_58_40.png|}} et spécifiez '1' comme identifiant (ID).
+==== Convertissez le polygone en format raster ====
+  - Trouvez la fonction Rasterize sous Raster>Conversion. Spécifiez le fichier polygone que vous venez de créer à l'étape précédente comme fichier entrant, et spécifiez un nom de fichier avec une extension .tif comme fichier sortant. Choisissez 'Résolution exprimée en unité de la carte par pixel'. Vous devez ensuite spécifier la taille désirée de chaque pixel sur la carte en unités du système de référence du fichier polygone en entrée. Par exemple, si le fichier polygone est en coordonnées UTM, cette résolution est en mètre. Assurez-vous de choisir une résolution raisonnable vous permettant d'obtenir un fichier raster avec un nombre de pixels gérable (par.ex. quelques centaines de lignes et de colonnes, maximum 2000-3000). Par exemple, si votre aire d'étude fait environ 2km x 3km, vous pouvez spécifier 5 (mètres) comme résolution de façon à avoir un fichier raster d'environ 400 pixels x 600 pixels.
+  - Maintenant, ce fichier tif est enregistré sur votre ordinateur. Dans votre aire d'étude, tous les pixels ont une valeur de 1, alors que les pixels à l'extérieur ont une valeur de 0.
+==== Importer le masque dans R ====
+  - Vous pouvez maintenant importer le masque dans R:
+<file>
+mask<-readGDAL(file.choose())
+</file>
+et choisissez le fichier raster en format tif.
+ensuite, vous pouvez faire
+<file>
+mask@data[mask@data==0]=NA
+</file>
+pour changer tous les zéros en valeurs NA (nulles).
-Étape 4: Numérisez la Réserve Molson et l'Arboretum Morgan. Vous devez d'abord ajouter une nouvelle couche shapefile de type polygone (CRS NAD 83/UTM 18N) et ajouter une colonne pour le nom (Donnée Texte, largeur 80). Donnez à ce fichier shapefile un nom approprié (e.g. Reserves_McGill). Activer le mode édition (>Couche>Basculer en mode édition) et numérisez les deux réserves en cliquant sur l'icône {{::capture_du_2013-09-25_10_58_40.png|}}. Pour vous aider à numériser les deux réserves, utilisez les points que vous avez importés à partir du fichier CSV. Pour la réserve Molson, numérisez la forêt située entre le boulevard résidentiel (Boul Perrot) et l'autoroute 20. Pour l'Arboretum, numérisez la région boisée située autour du point central. L'Arboretum fait approximativement 250 hectares. Ajoutez maintenant des polygones adjacents pour indiquer de quelle façon vous prévoyez étendre les réserves. Pour ce faire, il est important de sélectionner l'option 'Éviter intersections' ('Avoid int.) pour cette couche dans >Préférences>Options d'accrochage.
+Votre masque est maintenant prêt à être utilisé!