El objetivo del presente ejercicio es realizar la extracción por mascara (corte) de un grupo de archivos raster a partir de un shapefile, este usado como máscara, mediante un ciclo en R, es un proceso común cuando se tienen, por ejemplo, la precipitación (o alguna otra variable) para los 12 meses del año de todo Colombia y se quiere obtener solamente la precipitación de una unidad administrativa de menor área, para el presente caso es el departamento del Valle del Cauca; evitando hacer el proceso repetitivo, por ejemplo, en ArcGIS usando la herramienta “Extract by Mask”. Se hace uso de los ciclos, básicamente consiste en decirle al sistema que repita en cada una de las capas de un listado el corte a partir de una capa que representa el límite de nuestra área de estudio.

Luego abrimos el shapefile, que representa la máscara, bajo el comando “shapefile”.

Para poder visualizar ambos archivos hacemos uso del comando “plot”.

Luego hacemos el uso del comando “paste” para indicarle a R que lea todos los archivos que contienen la dirección escrita, esto permite unir la dirección, el 1 a 12 significa cada mes (es como finaliza el archivo) y por último la extensión que es .asc.

Seguido de ello hacemos uso del comando “lapply” para indicar que todos los archivos son raster.

Luego realizamos el ciclo, aspectos a tener en cuenta:

·         El for es el comando que se utiliza

·         La “i” es para indicar que cada vez que esta letra cambie, conforme a la secuencia (1:12) será un archivo distinto al cual se le hace el proceso.

·         Lo que está dentro de las “{}” será lo que repetirá una y otra vez hasta finalizar.

·         El comando “crop” se usa para realizar el corte (equivalente al “extract by mask” de ArcGIS), en la primera posición va el dato de entrada, y en la segunda posición va el dato de máscara crop(rasters, mask),

·         El comando “mask” se usa para indicar que la máscara (extensión del archivo) será equivalente al objeto “mascara”.

·         El comando writeRaster es para guardar los archivos en una dirección de nuestro computador, y el comando “names» se utiliza para que guarde el archivo bajo el mismo nombre que el archivo de entrada.

Luego abrimos el shapefile en R bajo el comando “shapefile”, y el raster base con el comando “raster”, tal como se muestra a continuación.

y = raster(«E:/Blogger/Post_6/_data/ mascara_col.tif»)

Para poder visualizar ambos archivos hacemos uso del comando “plot”.

Luego escogemos un campo de los que tiene el shapefile para hacer sobre este la conversión, para este caso el campo sobre el cual haremos la conversión sera el código del departamento según el código del DANE.

Seguido de ello hacemos la conversión bajo el comando “rasterize”:

Por último, guardamos este archivo raster en nuestro pc bajo el comando “writeRaster”.

plot(conversion)

Luego abrimos la tabla en R bajo el comando “read.table”, como se muestra a continuación.

Seguido de ello convertimos la tabla a formato shapefile, teniendo en cuenta que las columnas 1 y 2 corresponden a la longitud y latitud respectivamente:

Podemos ver el tipo de clase a la que queda convertido el objeto “tabla” bajo el siguiente comando:

Por último, procedemos a realizar el guardado de este archivo.

De este link se puede descargar el dato que se hace uso en este tutorial.

Luego, hacemos un listado de los archivos raster (.tif), seguido de aplicar la función raster para que el programa los reconozca como archivos espaciales, proseguido se realiza un stack que es sobreponer todos los raster, posteriormente se multiplica por 100, y luego se desalinean de la unión que se hizo con el stack.

Por ultimo realizamos un ciclo para escribir los archivos dentro de una carpeta con su respectivo nombre.