Suite

Découpez raster par raster avec extraction de données et changement de résolution


J'ai besoin d'empiler plusieurs cartes (tif) ensemble et coupez-les tous en fonction de l'étendue d'une petite carte raster (deuxième raster en haut de l'image ci-dessous). Cependant, le la taille des pixels dans chaque carte est différente. Existe-t-il un moyen de visualiser le processus ci-dessous ?

S'il vous plaît, répondez uniquement à la solution R ou à la solution QGIS.

bibliothèque(raster) e1<-extent(c(0,6,0,6)) r1<-raster(nrows=3,ncols=3,ext=e1) valeurs(r1)<-c(0,4,1 ,0,1,1,1,2,3) plot(r1) e2<-extent(c(2,5,2,5)) r2<-raster(nrows=4,ncols=3,ext=e2) valeurs(r2)<-c(0,8,2,1,4,9,9,4,0,8,0,0) tracé(r2)

Vous pouvez:

Chargez les bibliothèques requises :

bibliothèque(raster) bibliothèque(rgdal)

Lire les rasters :

r1 = raster("./dir/r1.tif") r2 = raster("./dir/r2.tif")

Rééchantillonner à la grille plus fine

r.nouveau = resample(r1, r2, "bilinéaire")

Si nécessaire (pour le masquage), définissez les étendues pour qu'elles correspondent

ex = étendue(r1) r2 = recadrage(r2, ex)

Suppression des données inutiles

r.nouveau = masque(r.nouveau, r2)


Voir la vidéo: Arcgis: Merge the raster datasets fusionner raster arcgis Mosaic raster dataset ArcGIS (Octobre 2021).