Introducción

El agave azul (Agave tequilana Weber, var. azul) se utiliza, principalmente, para la producción del Tequila, una bebida espirituosa (^{Bautista-Justo et al., 2001}); y recientemente para la producción de miel de agave e inulina. En México, el agave azul se cultiva en la región conocida como Denominación de Origen del Tequila (DOT) reconocimiento que otorgó la Unión Europea a México en 1997 y que comprende a los estados de Jalisco (126 municipios), Michoacán (30 municipios), Tamaulipas (11 municipios), Nayarit (8 municipios) y Guanajuato (7 municipios), contribuyendo Jalisco con el 90% de la producción (^{González et al., 2007}).

Desde finales de la década de los 90s, las plantaciones de agave han presentado problemas debido a factores que restringen el desarrollo normal de la planta. En los primeros años de la década pasada, a través del Consejo Regulador del Tequila, instituciones de investigación, enseñanza y productores, desarrollaron trabajos para determinar las causas del problema, enfocándose a enfermedades y plagas; entre los más estudiados están el daño por el picudo del agave (Scyphophorus acupunctatus) y la marchitez del agave atribuida a la bacteria Erwinia caratovora (^{González et al., 2007}).

Una restricción es un factor o grupo de factores que afectan negativamente el desarrollo de un organismo; por ejemplo, la luz solar como forma de radiación electromagnética influye en procesos fisiológicos de plantas y animales (^{Valladares, 2004}). Al considerar una restricción como factor negativo en el desarrollo de la planta que impide su manifestación como planta vigorosa, y al estar afectada la plantación por algún tipo de restricción, ésta se manifiesta físicamente con características que el observador puede distinguir a simple vista.

En la presente investigación, se usó el concepto restricción, según ^{Bozinovic y Canals (2009)} en el sentido de limitación. Es decir, según estos autores, “ningún organismo puede hacer de todo pues está limitado por sus genes, su historia, sus características biofísicas y mecánicas, su fisiología y conducta, estilo de vida y por el ambiente. En este sentido existirán restricciones internas (históricas) y externas (ecológicas). Dentro de estas restricciones existirán compromisos y plasticidad que dictarán el conjunto posible de características”.

Las tecnologías de información geográfica (TIG), incluyen a la percepción remota, que actualmente utiliza imágenes de satélite multiespectrales que se aplican para estimar tipos de cubiertas sobre la superficie terrestre y la condición de las mismas (^{Lencinas y Mohr-Bell, 2007}). De esta manera, se puede aplicar esta tecnología para identificar las superficies con plantaciones de agave, diferenciar sus condiciones de desarrollo y discriminar aquellas superficies cuyas características en la imagen sugieran alguna restricción.

^{Flores et al. (2008)}, mencionan que para el agave existen estudios concernientes a su fisiología, enfermedades y condiciones ambientales, más no hay estudios documentados en el uso de tecnología de satélite para su localización. El Consejo Regulador del Tequila, presenta un artículo titulado “El Agave Azul a través de la Imagen de Satélite” (“Blue Agave Through a Satellite Image”) donde se puede visualizar la metodología para cuantificar predios con plantaciones de agave (^{CRT, 2004}). Los trabajos publicados con respecto a este cultivo están enfocados a producción de tequila, aspectos ambientales, de manejo, de la composición propia del vegetal o de la organización para la producción, y en el caso de la aplicación de tecnología satelital se refieren a cultivos “comunes” (^{Manjarrez-Domínguez et al., 2007}; ^{Lencinas y Mohr-Bell, 2007}; ^{Vega et al., 2008}).

Dado que cualquier fenómeno que ocurra sobre la superficie de la tierra: vegetación, cultivos, cuerpos de agua, suelos etc. se puede detectar y analizar mediante la tecnología de la percepción remota (^{González, 2004}), es factible inferir la manifestación de la “salud” del cultivo o de su condición de desarrollo a través de la interpretación de imágenes de satélite.

Una imagen de satélite es un arreglo de columnas y renglones que conforman una matriz de datos numéricos que representan la intensidad de la energía electromagnética reflejada o emitida por los objetos en la superficie de la tierra (^{Jensen, 2005}; ^{Smith, 2012}). Las imágenes se pueden registrar en bandas individuales del espectro electromagnético, como ocurre en los satélites del programa LANDSAT; es decir, la misma escena capturada en diferentes bandas; Banda 1, Banda 2, Banda 3, etc., que en conjunto conforman una imagen multiespectral. Estas imágenes son adquiridas transformando la energía recibida por el sensor a valores de 8 bits; el pixel toma un valor de 0 a 255, dependiendo de la energía reflejada por los objetos. Estos valores, que se conocen como números digitales (ND), son los datos originales (^{Chuvieco, 1995}). Las imágenes del programa de satélites LANDSAT son consideradas como de resolución media con un tamaño de pixel de 30 m en multiespectral y 120 m en infrarrojo térmico (^{Smith, 2012}).

Se utilizó la imagen del satélite LANDSAT 5 adquirida el 2 de abril del año 2011; a la cual se le realizaron los procesos de calibración radiométrica, corrección atmosférica y georreferenciación.

El análisis con fines de obtener información se puede hacer mediante procesos digitales aplicando las estadísticas de la imagen (clasificación supervisada, clasificación no supervisada) o mediante el análisis visual a través de un compuesto a color (^{Lillesand y Kieffer, 1987}; ^{López y Caselles, 1991}; ^{López et al., 2008}).

El objetivo de la presente investigación fue identificar a través de la imagen LANDSAT 5, la superficie con plantaciones de agave, la presencia de agave con restricciones en su desarrollo, su cuantificación y ubicación.

Materiales y Métodos

El área de influencia de la Denominación de Origen del Tequila (DOT), es de carácter interestatal y las superficies ocupadas por agave se distribuyen de manera dispersa aunque en ciertas localidades tiene una mayor concentración (^{Valenzuela, 1994}; ^{González et al., 2007}).

Esta investigación se llevó a cabo en los municipios de El Arenal, Amatitlán y Tequila, zona tradicional del cultivo de agave y producción de tequila. El área de trabajo se distribuye dentro de las coordenada extremas Noroeste: 21° 00’ 10.7’’ N, 103° 54’ 10.3’’ O y Sureste: 20° 40’ 17.1’’N, 103° 33’ 40.6’’ O, para el sistema de referencia WGS84. En este polígono es donde se localiza la mayor concentración de predios ocupados con agave, para estos municipios. Presenta una fuerte influencia volcánica (Volcán de Tequila), se incluye en la provincia fisiográfica conocida como Eje Transversal Neovolcánico (^{SPP, 1981}), que determina en gran medida las características de superficie y las condiciones de suelo, principalmente pedregosidad y relieve ondulado que orientan el uso hacia el cultivo del agave (^{Demant, 1978}).

Para el reconocimiento por percepción remota de factores que restringen el desarrollo del cultivo del agave, se utilizó la imagen LANDSAT 5 del 2 de abril de 2011, adquirida durante la estación seca enero-junio, que permite un mejor contraste entre los predios con agave y otras cubiertas.

El satélite LANDSAT 5 está equipado con el sensor denominado Mapeador Temático (Thematic Mapper TM) que captura la radiación proveniente de la superficie terrestre en siete bandas espectrales, cada una para registrar características de objetos en la superficie: suelos, vegetación, agua, etc. (^{López y Caselles, 1991})

En el aspecto cartográfico, se utilizaron archivos imágenes de cartas topográficas a escala 1:50000, editadas por el Instituto Nacional de Estadística Geografía e Informática (INEGI) que abarcan el área de proyecto; estas cartas fueron F13D54, Tequila; F13D55, Tesistán; F13D64, Tala y F13D65, Guadalajara oeste, así como, archivos vectoriales para estas mismas cartas. Los archivos de imagen topográfica se georreferenciaron y unieron en mosaico mediante la herramienta Geocoding Wizard y Mosaiking (^{ER Mapper, 1999}).

Como un recurso abierto de las tecnologías de la información geográfica, se ha hecho uso del sitio Google Earth, tanto para obtener escenas de interés como para ubicar predios con coordenadas geográficas en datum WGS84, obtenidas en campo. Las imágenes que maneja Google Earth, para el área de estudio son imágenes de alta resolución que permiten hacer acercamientos a nivel de líneas de agave.

Para la ubicación de los predios en campo se usó un navegador satelital marca Garmin Etrex de 12 canales. Este navegador ofrece posiciones con un rango de precisión de hasta 4 m, que para fines de ubicar y registrar el predio que se visita para la obtención de datos, es suficiente. Para el procesamiento de datos de la imagen, se usó una computadora y el software ER MAPPER 6 (^{ER Mapper, 1999}).

Actualmente, la localización geográfica de cualquier punto sobre la superficie de la tierra se hace con base en sistemas de referencia. Estos son modelos que establecen la relación entre la forma de la Tierra y sus dimensiones. El primero de estos modelos utilizado para México por la Comisión de Estudios del Territorio Nacional (CETENAL) creada en 1968 y actualmente Instituto Nacional de Estadística, Geografía e Informática, fue el sistema de referencia conocido como North American Datum, 1927 (NAD27) (^{Hansen, 2005}).

Se generó el mapa base para localizar cualquier punto o rasgo sobre la superficie del área de proyecto y a la vez referir a este mapa base los materiales de imágenes que se fueran obteniendo. Esto se hizo georreferenciando las cartas topográficas de INEGI a escala 1:50000 del área de estudio, a su sistema de referencia inicial NAD27. Posteriormente, se procesaron estos materiales para unirlos en mosaico y por último hacer la reproyección del mosaico obtenido al sistema de referencia WGS84. Estos procesos se hicieron según el Manual de Trabajo de ER MAPPER 6 (^{ER Mapper, 1999}).

Las bandas de la imagen se procesaron para su calibración radiométrica, corrección atmosférica, georreferenciación. Esto se realizó para toda la escena 29/46 (path/route) en su conjunto, banda por banda y ya corregida por estos conceptos se le realizó un recorte para dejar solo el área de estudio.

Se generó el compuesto a color usando la Banda 4 (infrarrojo cercano), Banda 5 (infrarrojo medio) y la Banda 6 (infrarrojo térmico) asignando canal rojo, verde y azul, respectivamente. Este manejo de banda dio como salida una imagen en compuesto a color (falso color), que fue la que mejor separó las cubiertas o usos del suelo y que a su vez ayudó a diferenciar el cultivo del agave azul de otras cubiertas y de sí mismo, incluyendo edades y estados de desarrollo (^{Flores et al., 2008}).

Los sitios en los que los patrones de color sugirieron cambios en la cubierta se marcaron con puntos para reconocimiento en campo, los que fueron examinados para la captura de información. La información obtenida hace referencia al estado del cultivo o cubierta y condiciones de sitio, registrándose en fotografías, audio y video. El trabajo de campo se realizó con apoyo del Consejo Regulador del Tequila (CRT) y su supervisor de campo para la zona. Se trabajó durante los años 2011, 2012 y 2013.

La imagen del LANDSAT 5, está definida por su resolución espacial (30 m el tamaño del pixel para las bandas del visible e infrarrojo y 120 m para la banda térmica), resolución espectral (siete bandas del espectro electromagnético), resolución radiométrica (imagen de 8 bits, el valor del pixel tiene un rango de 0 a 255) y su resolución temporal (lapso de revisita de 16 días).

La calibración radiométrica (^{Chander et al., 2009}) se realizó banda por banda utilizando el software ENVI 4.7. El número digital o nivel digital (ND) a través del proceso de calibración radiométrica es transformado a cantidad de energía (radianza) y ésta a valor de reflectancia (^{Chander et al., 2009}).

El siguiente proceso en el tratamiento de la imagen fue el de corrección atmosférica. La energía recibida por el sensor es la suma de la energía reflejada por los objetos en la superficie y la energía dispersada y reflejada por la atmósfera. Fue necesario corregir estos efectos de la atmósfera para dejar sólo los correspondientes a los objetos en la superficie.

El método empleado fue el de pixel oscuro; mínimo del histograma. Este se basa en que la reflectividad en el océano (cuerpos de agua) es prácticamente nula, toda la energía recibida desde los pixeles del océano proviene de la atmósfera; estos pixeles del océano coincidirían con los valores mínimos del histograma de las imágenes (^{Chávez, 1988}).

De las estadísticas de cada banda de la imagen se obtuvo el valor mínimo del pixel, diferente de cero. Ese valor mínimo corresponde al cuerpo de agua que absorbe la energía incidente y la que llega al sensor corresponde a la reflejada por la atmósfera.

Este valor mínimo se restó a los valores de los pixeles de toda la banda y por cada banda, obteniendo una nueva imagen con las siete bandas corregidas por el efecto de la atmósfera. Respecto a la georreferenciación, la imagen ya se obtiene corregida geométricamente (nivel 1). La corrección geométrica consistió en la verificación de este proceso. Para este caso se comparó la imagen de trabajo, con imágenes vectoriales editadas por INEGI; cartas F13D54 Tequila y F14D55 Guadalajara oeste a escala 1:50 000 coincidiendo rasgos vectoriales con rasgos de la imagen de manera fiel.

Con estos tres procesos asegurados: calibración radiométrica, corrección atmosférica y rectificación geométrica, el siguiente proceso fue el recorte del área de estudio. Este consistió en tres etapas: 1) Hacer compilación de las 7 bandas de LANDSAT agrupadas en un sólo archivo (Layer Stacking), para tener la escena completa de la imagen en un solo archivo. 2) Obtener una subescena que incluye el área de interés. Esta subescena se distribuye en cuadrante noroeste entre las coordenadas extremas noroeste 21° 00’ 10.7’’ N, 103° 54’ 10.3’’ O y sureste: 20° 40’ 17.1’’ N, 103° 33’ 40.6’’ O. 3) Trazar el polígono que delimita el área de estudio, y su conversión a REGIÓN para segregarlo del resto y dejarlo como área de estudio. Estos procesos se realizaron con las herramientas del programa ER MAPPER 6 (^{ER Mapper, 1999}).

La interpretación digital de la imagen se hizo sobre el área de estudio segregada y considerando las 7 bandas de la imagen. El método utilizado fue el de Clasificación Supervisada (^{ER Mapper 6, 1999}) que opera con campos de entrenamiento (muestras) de datos reales para aplicarlos a los datos de toda la imagen a clasificar. En la formación de grupos o clases, éstos deberán corresponder a los diferentes tipos de usos o cubiertas asignados en los campos de entrenamiento.

Se utilizaron cuatro campos de entrenamiento para cada una de las cubiertas identificadas: bosque de pino, bosque de encino, matorral, etc. y, por supuesto, agave en diferentes edades y estados de producción.

En este proceso se utilizó la opción de máxima probabilidad. Se obtuvo un mapa de clases en base a los tipos de uso o cubiertas indicadas en los campos de entrenamiento, así como, el listado de las estadísticas para cada clase de cubierta de donde se obtiene la superficie total del área de proyecto y la superficie de los diferentes usos de suelo o cubierta. Para la evaluación de la precisión de esta información se trabajó una matriz de error (^{Conglaton, 1991}; ^{Jensen, 2005}), utilizando puntos de distribución aleatoria y comparándolos con la verdad de tierra.

Adicionalmente se consideró conveniente mostrar los valores de reflectancia para cada una de las cubiertas o uso del suelo reportados por el proceso de clasificación supervisada para la banda 4 (infrarrojo cercano); así mismo, se muestran los valores del índice de vegetación (NDVI), para las mimas cubiertas.

Resultados y Discusión

A las restricciones identificadas por las observaciones de campo, los patrones de color reflejados en la imagen, y el resultado de la clasificación supervisada, se les agrupó en tres clases:

En el compuesto a color de la imagen de trabajo, y con un contraste de 0.98, los predios con Restricciones Tipo 1 que son los desatendidos, los asociados con maíz o los invadidos por maleza, presentan tonos claros con variantes tenues a verde, azul, o gris rosáceo que en el sistema RGB corresponden a las mezclas 181, 244, 253 para el azul tenue y 190,143, 179 para tonos gris rosáceo. No debe ignorarse que estas condiciones del cultivo pueden confundirse con otras cubiertas como pastizales o vegetación arbustiva debido a que el agave desatendido se comporta como una especie arbustiva más. Por otro lado, la interpretación se hace sobre una imagen adquirida en periodo seco en el que la actividad fotosintética solo se presenta en especies forestales, cultivos de riego, agave y vegetación acuática. (^{Pimienta et al., 2006}).

En el caso de predios de agave afectados por alguna enfermedad, las hojas medias y cogollo tienden a secarse, lo que se observa principalmente en predios mayores a 4 años y en la imagen se registran con tonalidades tendientes al naranja y al azul púrpura, que en el sistema RGB corresponden a las mezclas 250, 196, 84 para el tono naranja y 117, 50, 190 para el tono azul púrpura. Igualmente, estos tonos se pueden confundir con pastos, céspedes o caña de azúcar y con agave de 2 ó 3 años, respectivamente; la discriminación corresponderá al responsable técnico del área.

Para el grupo 3, o restricciones tipo 3, en el compuesto a color se manifiesta en un azul-púrpura; advirtiéndose que para este caso pudiera confundirse con agaves de 2 a 3 años de edad.

Este análisis visual, en primera instancia, permite focalizar los predios con problemas, así como los predios con buen desarrollo y estimar cualitativamente una expectativa de rendimiento o bien dirigir acciones a predios con restricciones susceptibles de ser beneficiados.

Del proceso de clasificación supervisada, se obtuvo un mapa temático de distribución de las cubiertas clasificadas como agave y no agave (Figura 1). Así mismo, se generó otro mapa que muestra las restricciones en clases, donde el color azul representa a las restricciones tipo 1 asociadas con manejo; el naranja las restricciones tipo 2 asociadas a la sanidad de la planta, y en rojo, restricciones tipo 3 en donde el deterioro por desecación de la planta está asociado a exceso de humedad en terrenos planos (Figura 2).

Figura 1: Resultado de la clasificación supervisada que muestra la distribución de superficies ocupadas por agave contra otros tipos de usos del suelo (no agave). 

Figura 2: Resultado de la clasificación supervisada que muestra la distribución de las superficies con restricciones en el cultivo. Restricción tipo 1 = aspectos económicos; Restricción tipo 2 = plagas y enfermedades; Restricción tipo 3 = efectos por exceso de humedad del suelo. 

La Figura 3 representa el mapa temático de restricciones sobre el mosaico de carta topográfica a escala 1:50000; este mapa es de especial interés ya que los predios con restricciones son ubicables en campo.

Figura 3: Mapa de clases con restricciones y su base cartográfica. Restricción tipo 1 = aspectos económicos; Restricción tipo 2 = plagas y enfermedades; Restricción tipo 3 = efectos por exceso de humedad del suelo. 

Se obtuvo la superficie en hectáreas por uso del suelo o cubierta terrestre, por medio de la clasificación supervisada (Cuadro 1). El proceso reportó para las Restricciones Tipo 1 asociadas al manejo del cultivo, 4601 ha; para las Restricciones Tipo 2 asociadas a la sanidad del cultivo, 729 ha y para Restricciones Tipo 3 asociadas a excesos de humedad en terrenos planos, 630 ha que en total suman 5962 ha que representan el 45% del total ocupado por agave en la zona de estudio. ^{Uvalle-Bueno et al. (2006)} mencionan que en estudio fitosanitario del CRT de 1997, al menos el 23% del inventario del agave tequilero tenía algún grado de enfermedad, valor que prácticamente se duplicó en 2005. Este dato comparable con lo obtenido en la presente investigación (año 2011).

Restricción tipo 1 = aspectos económicos; Restricción tipo 2 = plagas y enfermedades; Restricción tipo 3 = efectos por exceso de humedad del suelo.

Cuadro 1: Superficies obtenidas de la clasificación supervisada para diferentes usos del suelo. 

Los datos de superficie acorde con el uso de suelo (Cuadro 1), junto con el mapa de clases donde se muestra la distribución del agave con restricciones en su desarrollo, permitirán hacer evaluaciones sobre lo que se espera en rendimiento del cultivo y tomar las acciones que esto sugiera.

La información publicada sobre el cultivo del agave azul y su relación con imágenes de satélite, ya sea para fines de identificación, ubicación, monitoreo o estado de desarrollo, es escasa o prácticamente nula. Las experiencias son de la segunda mitad de los años 90’s y principios de los años del S. XXI, no documentadas. No obstante el uso de imágenes satelitales para el estudio de recursos naturales como vegetación, suelos, agua y minerales si está documentada (^{Valladares, 2004}; ^{Lencinas y Mohr-Bell, 2007}; ^{Manjarrez-Domínguez et al., 2007}), utilizando imágenes en compuesto a color, procesos de clasificación no supervisada y supervisada e índices de desarrollo vegetal, así como procesos estadísticos. Estos trabajos han mostrado su utilidad para separar cubiertas o usos del suelo que son diferentes entre sí en la realidad.

Esta diferenciación de cubiertas fue congruente con el mapa de clases obtenido por la clasificación supervisada (Figura 2), logrando separar superficies ocupadas con agave con problemas de desarrollo, al contrastar el mapa de clase con el compuesto a color. La evaluación de la precisión, a través del análisis de la matriz de error (^{Conglaton, 1991}; ^{Jensen, 2005}) se muestra en el Cuadro 2, cuya precisión general resultó del 73%.

Restricción tipo 1 = aspectos económicos; Restricción tipo 2 = plagas y enfermedades; Restricción tipo 3 = efectos por exceso de humedad del suelo.

Cuadro 2: Resumen de la matriz de error: clasificación supervisada. 

De este análisis, para las clases de agave con restricciones se concluye que el resultado de la clasificación supervisada es de utilidad, considerando que los errores o confusiones se dan entre clases cuyo aspecto en campo son similares. Cabe señalar que prácticamente no se presenta confusión con ninguna de las otras clases de agave en diferentes edades de desarrollo, salvo la clase Agave Restricciones 3 que tiende a reflejar la misma matiz de color que la clase Agave 2008 debido a la influencia del suelo.

Dadas las características del cultivo, objeto de esta investigación, su manejo y desarrollo, es de esperarse que en el procesamiento de la imagen siempre haya un riesgo de confusión; así se apreció en la presente investigación, pero esta confusión puede ser aminorada con la interpretación del personal a cargo del manejo de esta tecnología y con conocimiento del hábito del cultivo. No obstante, el resultado es una herramienta básica para verificar su corroboración en campo y establecer criterios de decisión sobre el manejo del cultivo según las condiciones y restricciones encontradas.

El Cuadro 3 muestra los valores de reflectancia para cada una de las cubiertas reportadas. Se presentan los valores de la Banda 4 (infrarrojo cercano) por ser la que mejor responde a la vegetación (^{Jackson y Huete, 1991}), asimismo, se incluyen los valores de NDVI; también adimensionales, obtenidos en el proceso de datos de la imagen de trabajo con el propósito de mostrar la respuesta, especialmente de las plantas de agave al estímulo de la energía incidente.

Restricción tipo 1 = aspectos económicos; Restricción tipo 2 = plagas y enfermedades; Restricción tipo 3 = efectos por exceso de humedad del suelo.

Cuadro 3: Valores de reflectancia para la Banda 4 y valores de NDVI por uso del suelo o cubiertas registradas. 

La vegetación de agave de 2004 a 2007 varió en el rango de 0.241 a 0.294. Para esta banda, B4, en este rango ocurren el agave con restricciones 1, el agave con restricciones 2, el bosque de pino y la vegetación secundaria. Las cubiertas de agave 2008, 2009 y 2010 presentaron un rango de 0.142 a 0.195, menor a las cubiertas de agave estudiadas, lo cual se explica por ser plantas menores a 3 años. El valor más alto lo reporta la vegetación bajo riego con 0.434. Con respecto al índice de vegetación (NDVI), los valores que comparativamente presentan un valor más alto de 0.301 y 0.357 correspondieron al agave 2008 y 2007 respectivamente, esto para las edades de 3 a 4 años. En el resto de cubiertas, la vegetación de riego con 0.769 y el bosque de pino con 0.609 son las que reportaron los valores más altos de NDVI, lo cual se debe a que la vegetación de tipo forestal cubre totalmente el suelo y la vegetación bajo riego se mantiene activa en la época de sequía (noviembre - junio).

Conclusiones

El uso de imágenes LANDSAT por su resolución espectral y radiométrica es aceptable para el estudio del cultivo del agave azul según lo observado en campo y lo expresado por el procesamiento de la imagen. El proceso permite la identificación puntual del cultivo y su condición de desarrollo, lo que ayuda a focalizar la toma de decisiones.

En el presente estudio el 45% de la superficie ocupada con agave reflejó problemas en su desarrollo; esto representa, aproximadamente 5962 ha, con una confiablidad obtenida de la matriz de error del 73%. Este hallazgo justifica la aplicación de la metodología usada en acciones de seguimiento del cultivo en su desarrollo.

En la presente investigación la Banda 4 y el Índice de Vegetación del agave con restricciones reporta valores de 0.192 a 0.282 y de 0.235 a 0.261, respectivamente. Por la escasez de información sobre el tema y los resultados reportados, este trabajo representa un esfuerzo que retoma acción realizadas en el pasado reciente sobre esta temática y permite abrir una etapa de trabajos que amplíen el estudio del agave con esta herramienta metodológica.