Tutorial QGis
Eje mpl o: Conte o de pl antas .
1
[Link]
Índice
1. Cargar archivo
2. Dibujar lotes
3. Crear cuadricula
4. I dentificar
5. Visualización
2
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
1. Cargar ortofoto con valores NDVI
Lo primero que debe hacer es cargar el archivo ortomosaico con los valores de NDVI. Podemos
cargar también el ortomosaico original si lo desea para una comparación visual.
Figura 1. Cargar archivos
Para eso debe hacer clic en la opción administrar fuentes de datos y se desplegara la siguiente
ventana:
Figura 2. Cargar archivos Ráster.
Recuerde que estamos cargando una imagen por lo que debemos seleccionar la opción Raster.
Para seleccionar el archivo deseado debe hacer clic en los en los 3 puntos y seleccionar el
archivo correspondiente a los valores de NDVI. Una vez seleccionado el archivo seleccionamos
añadir y cerramos la venta.
En la figura 3 se observa la vista del archivo. Tenga en cuenta que los valores exportados son
de la banda 1 termal y la banda 2 NDVI. Si no recuerda estos pasos vea el instructivo Agisoft-
Agro.
3
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Figura 3. Ráster.
Para poder visualizar los valores de NDVI, haga clic con el botón derecho sobre la capa de
NDVI, se desplegará un menú, haga clic en propiedades.
Figura 4. Opciones de ráster.
4
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Primero debemos cambiar el tipo de renderizado. Haga clic para desplegar el menú y elija la
opción pseudocolor monobanda.
Figura 5. Opciones de ráster.
Seleccione la banda 2, y elija la rampa de colores que desee desplegando la lista con el
triángulo negro a la derecha de la barra. Haga clic en la opción clasificar y aplicar para que
tome los cambios. Una vez realizado este paso reajuste los valores mínimos al valor de 0.
Aplique y acepte. Debería tener una imagen como se muestra a continuación:
Figura 6. Ráster a valores de NDVI.
5
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
2. Dibujar lotes
Para dibujar los lotes en donde queremos contar la cantidad de plantas debemos ir a la opción:
Capa Crear capa Nueva capa de archivos Shape.
Figura 7. Dibujar lotes.
Para crear los polígonos debe hacer clic en los 3 puntos para desplegar la siguiente ventana en
donde crearemos nuestro archivo vectorial de los lotes. Ponga un nombre y elija en la carpeta
donde guardara los archivos. Se recomienda prolijidad a la hora de almacenar los datos.
En el tipo de geometría, en este caso, elija la opción polígono y corrobore el sistema de
coordenadas en cual esta creando el polígono. Recuerde hacerlo en el mismo sistema que
exporto el archivo ráster.
Figura 8. Nueva capa de poligono.
6
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
En la vista de capas vera una nueva capa llamada como el nombre del archivo. Para comenzar
a dibujar los polígonos de los lotes debe seleccionar la opción de conmutar edición añadir
poligono. Recuerde estar en la capa del lote.
Figura 9. Crear polígono.
Para dibujar el polígono haga clic izquierdo en el primer vértice del polígono que quiera crear,
y un clic por cada lado. Una vez pintado el polígono deseado haga clic derecho para cerrar la
edición del mismo y asignarle un valor. Para dejar de editar deseleccione la opción de
conmutar edición previamente seleccionada. Guarde los cabios.
Una vez realizado el trazado vamos a dejar solo el contorno del polígono. Para ello haga clic
derecho sobre la capa shape que acabamos de crear (lotes) y acceda a las propiedades.
Recuerde estar clic en relleno sencillo
para ver el menú de la imagen
7
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Figura 10. Opciones del polígono.
En la opción de color de relleno, despliegue el menú seleccione el relleno transparente. Haga
clic en aceptar.
Figura 11. Color de relleno
Si saca la capa NDVI vera los lotes dibujados.
Figura 12. Contorno lotes dibujados
8
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
3. Crear cuadricula
Vamos a crear una cuadricula dentro de nuestros lotes para vectorizar las posiciones donde
puede encontrarse lechuga, siendo esta la de mayor índice NDVI en la escala.
Haga clic en la pestaña Vectorial Herramientas de investigación crear cuadricula
Figura 13. Crear cuadricula
En los parámetros para crear la cuadricula vamos a indicar lo siguiente según el orden de las
opciones (recuerde desplegar el menú con el triángulo negro o con los 3 puntos):
• Tipo de cuadricula:
hexagonal.
• Extensión cuadricula: opción
“usar la extensión de la
capa…” y elija la capa
correspondiente al polígono
del lote.
• Espaciado Horizontal y
vertical: elija el espaciado
que considere adecuado. A
menor espaciado más
preciso. Puede elegir 0.5
para probar.
9
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
• Si en las medidas de espaciado tiene la advertencia como en la imagen, es que tiene un
sistema de referencia inadecuado. Elija un sistema de coordenadas proyectado. En este
caso puede usar EPSG:21818.
• Ejecute
• Cerrar
Vera que nuestro ráster origina y nuestros lotes están por debajo de una capa llamada
“cuadricula”, observe la ventana lotes a su izquierda de la pantalla. El paso siguiente es limitar
esta cuadricula al contorno de nuestros lotes, lo que este por fuera de ellos es maleza o suelo
desnudo.
Haga clic en la pestaña Vectorial Herramientas de geoproceso cortar
Figura 15. Opción para unir lote con
cuadricula
Seleccione en capa de entrada la cuadricula y en capa de superposición la capa con los
polígonos del lote. Ejecute y cierre la ventana.
Figura 16. Cortar
10
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Una vez realizado el proceso, destile desde la ventana de las capas, las que no seguiremos
utilizando. En este caso destilde la capa de la cuadricula. Vera que ahora la cuadricula esta
dentro de los limites de mis lotes. Haga clic derecho sobre la capa llamada “cortado” ( puede
asignarle un nombre si lo desea en el menú anterior) y seleccione propiedades, al igual que
con los lotes, cambie el relleno a transparente.
Si hace zoom sobre la imagen vera algo así:
Figura 17. Zoom cuadricula recortada
4. I dentificar
Ahora le diremos al programa que realice una estadística zonal y seleccionaremos a los
polígonos que contengan un promedio de NDVI > 0.7.
En el extremo derecho de la barra superior encontrara el icono , oprímalo y sobre el
menú desplegado busque la herramienta llamada “Estadísticas de zona”.
• Capa ráster: capa con el
ráster NDVI.
• Banda ráster: banda 2,
recuerde banda 1 era
termal para este caso.
• Capa vectorial: la capa
vectorial cortado.
• Ejecutar y cerrar.
Figura 18. Estadísticas zonales
11
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Ahora elegiremos los valores de NDVI
que consideremos que son lechugas,
en este caso vamos a pedir mayores a
0.8. Haga clic con el botón derecho
sobre la capa “cortada” y abra la tabla
de atributos. En el menú superior de la
tabla de atributos abra el
seleccionador de objetos. Lo
encontrara con el siguiente icono:
En este seleccionador puede escribir la expresión para la cual quiere seleccionar datos. En este
caso vamos a escribir lo siguiente:
_mean es el promedio de los valores de NDVI dentro de los hexágonos de la cuadricula. Con esta
expresión solo me seleccionara los hexágonos con promedio de NDVI mayor a 0.8,
asegurándome que son lechugas. Haga clic en la opción inferior derecha “seleccionar objetos
espaciales”. Si Recorre el listado en la tabla de atributos encontrara resaltados en azul los valores
con promedio mayor a 0.7. Cierre la ventana.
Figura 19. Hexágonos promedio
12
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
EL siguiente paso es guardar en una nueva capa estos datos seleccionados. Haga clic derecho
sobre la capa cortada exportar guardar objetos seleccionados como
Figura 20. Capa datos seleccionados
Como en pasos anteriores, deseleccione la capa cortada, y observamos que los puntos están
sobre las lechgas. Necesitamos unificar estos puntos. Para eso haga clic:
Vectorial herramientas
de geoproceso disolver
En la opción de capa de entrada
elija la capa donde guardo los
promedios seleccionados.
Si realizo bien ese paso
encontrara que los hexágonos son
ahora polígonos.
Haga clic en Vectorial
herramientas de geometría
multiparte a
monoparte. Elija la ultima capa procesada en las opciones. (capa disuelta).
13
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Sobre esta capa volvemos a calculas las “estadísticas por zona” calculadas a principio del punto
4 pero esta vez elijiendo la capa de “monoparte”.
Como se necesita contar las lechugas reduciremos estos polígonos a sus centroides, ya que el
dato necesario, en este caso saber que es una lechuga, ya lo sabemos. Para realizar esto debe
Haga clic en Vectorial herramientas de geometría centroides
Elijo nuevamente la ultima capa realizada en este caso “monoparte”, ejecuto y cierro.
Figura 21. Lechugas
El objetivo del trabajo es saber la cantidad de lechugas por lote. Como se vio en el ejercicio se
tomó como parámetro el índice NDVI para discriminar a través de los datos de cámaras
multiespectrales las lechugas frente al suelo y malezas, en este caso identificadas debido a su
índice de NDVI menor al seleccionado para la lechuga. Siempre es recomendable una inspección
visual. En este caso puede cargar el ortomosaico exportado de Metashape con los colores
originales para su comparación.
Para contar las lechugas
automáticamente vamos a hacer clic en:
Vectorial herramientas de
análisis contar puntos en
polígono.
En las opciones el polígono es la capa del
lote y los punto son los centroides.
Ejecute y cierre.
14
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Para concluir destilde todas las capas excepto la que contiene los centroides y la capa llamada
numero.
5. Visualización
Sobre ambas capas realizara una modificación en la visualización. Para ello debe hacer clic
derecho sobre cada capa y elegir la opción propiedades:
Figura 22. Propiedades capas
En la ventana de propiedades debemos cambiar la opción de símbolo único a graduado
15
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
Para la opción de columna vamos a elegir el promedio de los valores de NVDI, puede cambiar el
icono y la graduación de colores. Una ves elegidas estas opciones haga clic en el botón clasificar,
y luego acepte.
Para la capa “numero” debe hacer el mismo procedimiento, pero en la opción de comuna elija
la opción NUMPOINTS, para que la escala represente el conteo de lechugas.
Puede luego variar los nombres de las capas a gusto haciendo siempre clic derecho sobre la capa
a modificar y elegir entre las opciones.
El resultado final sería el siguiente:
16
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
[Link]
17
© 2018 Aurea Imaging SRL Aráoz 2227 - Ciudad de Buenos Aires - Argentina (C1425DGG) +54 (11) 4831-5320
