Introduction à la géostatistique avec R

Guillaume Larocque - Université Laval, 26 février 2014

Présentation sur Prezi

Présentation simultanée

• 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.
• Geostatistics for Natural Resources Evaluation, Pierre Goovaerts. - Très bon livre pour une connaissance approfondie de la géostatistique.
• Site web de gstat.
• Package gstat - Package de gstat utilisé pour l'atelier.
• Package geoR - Autre package pour la géostatistique en R.
• Georges Matheron - La théorie des variables régionalisées et ses applications (PDF). - Text original décrivant les bases de la géostatistique.
• Analyses spatiales sous R - Document préparé par Nicolas Casajus pour un cours offert à l'UQAR à l'hiver 2013.

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.

1. Vous pouvez numériser le périmètre de votre aire d'étude à partir d'une couche Google Satellite en suivant les étapes suivantes:
2. - Activez le 'OpenLayers Plugin' dans les Extensions.
3. - 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.

É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 . 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.