texto em 
Inglês (pdf)
Artigo em XML
Referências do artigo
Tradução automática
Enviar este artigo por email
Citado por SciELO 
Similares em
SciELO 
Permalink
Introducción
De acuerdo con CONAGUA (2014) de la superficie de riego en México, 90% tiene problemas de eficiencias de aplicación bajas, que pueden mejorarse con la nivelación de tierras si las pendientes topográficas son adecuadas, lo que implica una mayor uniformidad del riego y una reducción del volumen de agua aplicado (Saucedo et al., 2006; Fuentes y Rendon, 2012).
La nivelación de tierras se justifica en cualquier proyecto de irrigación, ya que generalmente se invierten sumas considerables en obras de captación y distribución, comparativamente, se hacen bajas inversiones en la parcela, que es donde se refleja la bondad de todo un complejo sistema de irrigación (Walker, 1987). Se han desarrollado muchos programas de computadora relativos al tema de la nivelación de terrenos agrícolas, buscando siempre ahorrar tiempo en el diseño y realizar los cálculos óptimos con el menor costo posible (Covarrubias, 1968).
Dentro de los primeros programas en México se tiene a Nivterra v1.0 desarrollado para computadoras personales en Fortran 77, con mejoras posteriores en Pascal para su versión 2.0 (Chávez et al., 1990), dando origen a sistemas que se usan hoy en día en la práctica, como el Sinivet en sus distintas versiones utilizado para levantamientos por cuadrícula y para la nivelación sin divisiones (Hernández y Sánchez, 2008). Los pocos programas existentes en la actualidad para el diseño de proyectos de nivelación de tierras, usa como fuente de datos levantamientos realizados por cuadrícula. Como resultado de varios años de trabajo en el área de desarrollo de software aplicado a la nivelación de tierras se desarrolló el programa Terra 1.0, para procesar datos de campo que se obtienen con levantamientos topográficos con el método tradicional de levantamiento con puntos trazados en una cuadrícula rectangular o con una distribución variable (radiaciones).
Terra, busca adaptarse a las nuevas tecnologías que implican la toma de información de los datos de campo, para el uso y aplicación de la nivelación de tierras. El objetivo fue, el desarrollo de un sistema de cómputo electrónico, que incluye métodos de levantamiento topográfico con equipos modernos, como estaciones totales o GPS de precisión, que auxilie en el cálculo preciso y expedito para la elaboración de proyectos ejecutivos de nivelación de tierras agrícolas. El manejar una masa de puntos tomados por radiaciones en lugar de una cuadrícula, implica generar un modelo digital de elevaciones (MDE) que describa con alta confiabilidad el predio caracterizado. Esto se logra en Terra usando el método de Kriging para el modelado del terreno.
Materiales y métodos
La metodología incluye la aplicación de la interpolación espacial usando el método de Kriging (Villatoro et al., 2008) para la determinación de los pesos o factores de ponderación en el cálculo, estimación y cuantificación de volúmenes de corte y relleno y trazo de curvas de nivel del terreno antes y después de la nivelación, que cubren áreas que el topógrafo no lo hace en una cuadrícula homogéneamente espaciada en el procedimiento tradicional.
La solución del método de Kriging permite definir por el usuario un variograma de elevaciones en función del mejor ajuste realizado a su semivariograma experimental (Clark, 1979). Los modelos utilizados para el ajuste de los variogramas son: lineal, cuadrático, esférico, exponencial y gaussiano, con o sin el efecto pepita (Delhome, 1978). Al igual que Walker (1987); Cano (1997), Terra 1.0 incluye la división del terreno en tablas independientes propuestas por Marr (1957), en caso de que las pendientes del predio sean heterogéneas. En la Figura 1, se describe en forma simplificada el algoritmo usado para el desarrollo del sistema Terra 1.0.
Figura 1 Algoritmo de nivelación, implementado en Terra para calcular el plano proyecto, sin divisiones o en una sola tabla.
Definición de los datos de los levantamientos
Dependiendo del tipo de levantamiento topográfico realizado se obtuvieron coordenadas de campo para cada punto medido, ya sea en cuadrícula rectangular usando la propuesta de Marr (1957), en cuyo caso ya son coordenadas (xi, yj, zij), donde i= 1, 2, …, n; j= 1, 2, …, m; n= número de hileras y m= número de columnas de la cuadrícula o radiaciones, en cuyo caso son coordenadas (Xrk, Yrk, Zrk). Donde: k= 1, 2, …, P y P= número de puntos levantados por radiaciones; las cuales se introducirán el cuadro de diálogo del programa, activando la opción de captura correspondiente según el levantamiento realizado. De la misma forma se introducen las coordenadas de cada uno de los vértices de la poligonal del terreno con sus respectivas coordenadas (Xv, Yv). Donde: v= 1, 2, …, NV y NV= número de vértices de la poligonal del terreno.
Límites de trabajo, área y centroide
Con la información de los puntos de campo capturada, se determinan los límites extremos de una zona rectangular para el trabajo del proyecto de nivelación, definida por (Xmin,Ymin) en la esquina inferior izquierda y (Xmáx,Ymáx) en la esquina superior derecha. Para el caso de un levantamiento por radiaciones se crean los nodos de una cuadricula de dimensiones mxn en la zona rectangular definida, con una distancia entre nodos (Lcuadrícula), así para cada nodo se tiene información (xi, yj), faltando estimar zij. Con las coordenadas de los vértices de la poligonal se calcula el área y el centroide del terreno (xc, yc, zc).
Cálculo del plano proyecto
Éste se realiza por mínimos cuadrados ordinarios mediante una solución matricial como la propuesta por Montgomery y Runger (2008), en la cual es indistinta la forma en que se realizó el levantamiento topográfico, la ecuación del plano proyecto inicial es:
(1).
Donde:
2);
3);
4)
Donde: A= matriz de coordenadas XY; AT= transpuesta de la matriz de coordenadas XY; b= vector de coordenadas Z; β = vector de incógnitas o coeficientes (CA, CB, CC). De la ecuación (2) se despeja el vector de incógnitas o vector de coeficientes ( β ), para obtener los valores de las pendientes del plano (CB y CC) y la ordenada al origen (CA):
5);
6)
Cálculo de cuadrícula homogénea
Si el levantamiento fue mediante cuadrícula, ésta ya fue trazada en campo al hacer el levantamiento topográfico, pero si fue por radiaciones, la cuadricula se crea a partir de los límites calculados previamente y con una equidistancia regular de la cuadrícula asignada para el cálculo de tamaño Lcuadricula (m), que puede ser la misma de la cuadrícula regular u otra según sea el caso. Con esto se definen el número de nodos en X y en Y (Figura 2).
Figura 2 Elementos para determinar la cuadricula homogenea, a partir de un levantamiento topográfico por radiaciones y un MDE.
Los problemas que se detectaron en este punto se dan por tener terrenos muy cortos donde solo pueda definirse una linea, y por tanto, no se alcance a generar ningun cuadro, la solución más práctica es redefinir un valor más pequeño de la Lcuadricula. Con la información anterior, es posible determinar la localización (xi, yj) de cada nodo dentro de la cuadrícula.
Determinación de las coordenadas zij
En cada nodo (i, j) de la cuadrícula se determinan los valores de zij a través un proceso geoestadístico de interpolación espacial, partiendo de la informacion levantada en campo (Xrk, Yrk, Zrk) de los datos de campo definidos como puntos vecinos al nodo (i, j). La definición de puntos vecinos se realiza calculando la distancia de cada nodo (i, j) a todas las coordenadas obtenidas en el levantamiento topográfico (Xrk, Yrk, Zrk), los NC puntos con las distancias más pequeñas conformarán el conjunto de puntos vecinos o cercanos para la estimación de zij, mediante la ecuación (7).
7).
Donde: zij= dato estimado en el nodo (i, j) de la cuadrícula generada de dimesiones mxn; i = 1, 2, 3, n; j = 1, 2, 3, ..., m; n= número de hileras; m= número de columnas. Zq= coordenada Zrk del punto de campo ubicado en la posición q del conjunto de puntos vecinos o cercanos. Wq= peso o factor de ponderación calculado para el sitio q del conjunto de puntos vecinos o cercanos para la estimación de zij en el nodo (i, j). NC= número total de puntos vecinos o cercanos al nodo (i, j), que se utilizarán para estimar el dato zij. La suma de los valores de los pesos o factores de ponderación (Wq) del conjuto de puntos vecinos o cercanos al nodo (i, j) deben cumplir la restricción definida en la ecuación (8):
8).
Cálculo de los pesos o factor de ponderación (Wq)
El cálculo de lo pesos realiza aplicando el método de Kriging. Clayton y André (1998), indican que la escencia de Kriging es la estimación del semivariograma (Figura 3), en éste caso de elevaciones, el cual se construye definiendo 20 clases de distancias con todos los datos de campo (P) y calculando γ h para cada clase ubicada a una distancia (dt+h) el valor del semivariograma γ(h) através de la ecuación 9.
9)
Figura 3 Topografía del terreno natural, aplicando el método de Kriging, en un levantamiento topográfico de distribución variable por radiaciones.
Donde: Zr(dt)= valor de la cota topografía (Z) en el punto t; Zr(dt+h)= valor de la cota topográfica (Z) en el sitio ubicado a la distancia h del punto t; P= número total de puntos del levantamiento por radiaciones.
Calculado el semivariograma para cada una de las distancias definidas en las clases, se ajusta un modelo teorico (Cisneros-lturbe et al., 1999), que pude ser: lineal, cuadrático, esférico, exponencial o gaussiano (Cuadro 1) ecuaciones 10, 11, 12, 13 y 14.
Cuadro 1 Modelos matemáticos (Delhome, 1978; Clark, 1979), disponibles en Terra para el ajuste del variograma teórico, a partir del semivariograma experimental.
| Modelo | Función para (0 < h |
|
| Lineal |
|
(10) |
| Cuadrático |
|
(11) (11a) |
| Esférica |
|
(12) (12a) |
| Exponecial(*) |
|
(13) (13a) (13b) |
| Gaussian (*) |
|
(14) (14a) (14b) |
Cuando h> C2: γ(h)= C1 (para los modelos cuadrático y esférico); (*)= función aproximada a polinomio, mediante series de Taylor y para h= C2.
Michael y Clayton (2014), mencionan que los comportamientos en el origen de los modelos pueden considerar o no el afecto pepita; es decir, la intersección con el eje Y o bien, considerar que salen del origen, como es el caso de la topografía. Inicialmente se probó con modelos de varigrama que incluyen el efecto pepita; sin embargo, la mayoria de soluciones implicaba valores negativos para el coeficiente C0, razon por la cual, se decidió utilizar el modelo que sale del origen (Legrá-Lobaina et al., 2010).
El cálculo de los coeficientes de los modelos (C0, C1 y C2) se hizo con el metodo Newton-Raphson para la solución de sistemas de ecuaciones no lineales (Chapra y Canale, 2015), con resultados poco favorables, por lo que, los modelos se linealizaron (Cuadro 1) ecuaciones 11a, 12a, 13b y 14b, mejorando las soluciones.
Antes de la linealización, los modelos exponencial y Gaussiano se aproximaron a una función polinomial (ecuaciones 13a y 14a), mediante series de Taylor (Purcell y Varberg, 2009). Una vez linealizados los modelos se aplicó la ecuación (5) para determinar los coeficientes β i y para el calculo de los coeficientes C1 y C2 del variograma, las soluciones se listan en el (Cuadro 2) ecuaciones 15 a 20b.
Cuadro 2 Coeficientes C1 y C2, obtenidos para determinar el variograma teórico.
| Model | Coeficientes para variogramas partiendo del origen (C0= 0) | |||
| C2 (rango) | C1 (Sill) | |||
| Cuadrático |
|
(15) |
|
(15a) |
| Esférica |
|
(16) |
|
(16a) |
| Exponecial |
|
(17) (18) |
|
(17a) (18a) |
| Gaussiano |
|
(19) (20) |
|
(19a) (20a) |
Una vez calculados los coeficientes del variograma teórico, la estimación de los pesos o factores de ponderación (Wq) para cada uno de los nodos se hace mediante un proceso matricial (ecuaciones 21 y 22).
21)
22)
Donde: S= es una matriz cuadrada de orden n x n; cada elemento de la matriz se calcula evaluando el variograma teorico con la distancia entre puntos; por ejemplo: el elemento γ 12 es calculado al evaluar la función del variograma con la distancia existente entre el punto 1 y 2 de los NC puntos cercanos o vecinos que se involucraron para el proceso de interpolación. W y b= son los vectores de orden n x 1.
El vector b se obtiene al evaluar la función del variograma con la distancia que hay entre el nodo (i, j) donde se estimará zij y el dato vecino q, w contiene las n icógnitas involucradas en el cálculo; es decir, los pesos de cada uno de los NC datos, más un multiplicador de Lagrange (𝜆) que almacena el error residual calculado en el proceso, así al despejar el vector w se obtendrán los pesos o ponderaciones buscadas (Wq) w= S -1 b (23). Como producto de este paso se obtiene un archivo (*.MLL) con la información zij de cada uno de los nodos de la cuadrícula homogénea.
Ajuste de la ordenada al origen (CA)
El ajuste de la ordenada al origen (CA) del plano proyecto calculado en el segundo paso, se realiza de acuerdo a la relación entre los cortes y rellenos definidos por el coeficiente de abundamiento (R= corte/relleno), la mayoría de autores coiciden que el coeficiente de abundamiento (R) se debe encontrar en el rango: 1.1 ≤ R ≤ 1.5 (Walker 1988).
En la práctica, donde existieron rellenos habrá compactación y quedarán hoyos, estos se presentarán a través del tiempo por causa de la lluvia y compactación del suelo, por lo que es necesario, aumentar los cortes, por lo tanto disminuirán los rellenos esto se hace hasta tener un coeficiente de abundamiento (R) dentro del rango deseado, el ajuste se realiza bajando o subiendo el plano proyecto en ±Δh, h es el valor de incremento o decremento en la elevación del plano en m, dependiendo del suelo y su mayor o menor contenido de materia orgánica.
Es decir, si un suelo es arenoso el coeficiente R será cerca de 1.1, arcilloso el valor del coeficiente será cerca de 1.5. Si se considerara un suelo intermedio (franco arcilloso) los valores del coeficiente variarán entre 1.20 y 1.30 (Walker, 1988). El cálculo de los vólumenes de corte y relleno se hace mediante el método de los 4 puntos; es necesario contar con una cuadrícula con los datos de campo o generada por el método de Kriging.
La incorporación del método de Kriging para la estimación de la cota (zij) en cada nodo (i, j) de la cuadrícula es un procedimiento nuevo en nivelación de tierras y ha mostrado en la práctica resultados adecuados en el cálculo de volúmenes de corte y relleno. Como resultado de paso se obtiene un coeficiente CA ajustado (CAajus), para la ecuación del plano proyecto. Con lo anterior, se calculan valores finales de espesores de cortes o rellenos para cada nodo (i, j) de la cuadrícula y por tanto volúmenes de cortes y rellenos, así como costos de nivelación por hectárea y totales.
Generación y dibujo de isolíneas
Se elaboran planos de curvas de nivel para el terreno natural (configuración antes de la nivelación) y para el proyecto de terreno nivelado (configuración después de la nivelación). Éste proceso se realiza de dos formas: a) se crean curvas de nivel para todo el terreno, sin importar como fue realizado el levantamiento topográfico se crea una nueva cuadrícula más fina (x’i, y’j, z’ij) generada con longitudes entre nodos más pequeña, de aproximadamente un cuarto la longitud de cuadrícula para la nivelación.
Ésta nueva cuadrícula se calcula, a partir de los datos de campo levantados en cuadrícula rectangular (xi, yj, zij) o por radiaciones (Xrk, Yrk, Zrk) aplicando el proceso de interpolación espacial con el método de Kriging explicado en los puntos anteriores o el método del inverso de la distancia. Los límites de la nueva cuadrícula se extienden para cubrir perfectamente la poligonal del terreno y evitar áreas sin curvas. El cálculo y dibujo de isolíneas o curvas de nivel del terreno nivelado, se realiza con los nodos (i, j) de la nueva cuadrícula, cuya cota zij se obtiene con la ecuación final, ecuación (1) definida para el plano proyecto; y b) en caso de que por cambios abruptos de pendientes o tamaño del lote sea necesario dividir el terreno, para hacer su nivelación independiente por tablas, se diseñó un algoritmo que permite dividir el predio en dos partes de forma interactiva con el usuario.
En el ejemplo se muestra la topografía del terreno natural para que sobre ella se trace una línea divisoria que permita separar los puntos levantados en los nuevos polígonos y planos naturales de proyecto, que resultan con la línea divisoria. En caso de necesitar dividir el predio en más de dos partes se tiene una forma de captura para definir los polígonos.
Esto último con el fin de que el proyectista identifique manualmente los linderos de los polígonos y con esto queden dentro todos los puntos levantados en campo. Una vez definidas las tablas se aplicarán los seis puntos descritos anteriores para realizar el cálculo del plano natural de proyecto y final, así como el cálculo de volúmenes de corte y relleno en forma independiente de cada una de las tablas.
Resultados y discusión
Terra 1.0 se desarrolló con el lenguaje de programación Delphi (Pascal), diseñado para funcionar como herramienta en cualquier proyecto de nivelación de tierras donde la topografía sea levantada en cuadricula con trazo rectangular o radiaciones. El sistema permite la nivelación general de todo el predio levantado o la división por tablas.
El cálculo de volúmenes de corte y relleno se realiza con la cuadrícula rectangular aplicando el método de Kriging para el cálculo. La estimación de los valores de la cuadrícula genera valores congruentes con la topografía y por ende el cálculo de volúmenes se obtiene de acuerdo con los procedimientos técnicos preestablecidos.
Captura de la información
Terra, cuenta con tres cuadros de diálogo para la capturar la información de entrada: ‘datos del proyecto y parámetros del diseño’ para los datos generales, ‘datos de la poligonal del terreno y cálculo del centroide’ para capturar las coordenadas (Xv, Yv) de los vértices de la poligonal del terreno y ‘datos topográficos para la nivelación’, para los datos de campo con coordenadas (xi, yj, zij) o (Xrk, Yrk, Zrk) según el levantamiento topográfico realizado y la longitud de equidistancia de la cuadrícula.
Topografía del terreno natural
El programa permite generar y visualizar la topografía del terreno con curvas de nivel (Figura 4). Para poder dibujarlas se ha calculado un modelo digital de elevaciones usando el método de Kriging. Según Delhome (1979) recomienda el uso del método de Kriging en las aplicaciones a la hidrociencia, indicando que describe en forma adecuada el comportamiento de la topografía del terreno natural. La propuesta inicial para el cálculo de los coeficientes C0 y C1 del modelo del variograma teórico, se ajusta automáticamente al modelo lineal; sin embargo, se puede elegir otro en forma personalizada (Figura 3).
Levantamiento topográfico en cuadrícula rectangular
El software desarrollado permite procesar la información obtenida con un levantamiento con una cuadrícula rectangular, por lo que se puede aplicar la metodología tradicional para el cálculo de nivelación de tierras. En este caso ya existe la cuadrícula desde campo por lo que no hay necesidad de crear una nueva. Como es el caso del levantamiento por radiaciones donde se genera la cuadrícula.
División del terreno en tablas
El software desarrollado permite dividir el terreno en dos tablas, trazando una línea divisoria y generando los dos polígonos independientes en forma automática (Figura 5). Cuando se realiza un levantamiento por radiaciones, algunos puntos son tomados cerca de los linderos y con una alta probabilidad pueden caer fuera del terreno, el sistema asigna todos los puntos en alguno de los dos polígonos; sin embargo, cuando la división implica un mayor número de partes, será muy complicado agrupar los puntos, razón por la cual, se cuenta con una hoja de entrada de información donde el usuario debe ingresar las coordenadas de cada polígono y para cada tabla solo se extraerá la información que se encuentre dentro del polígono definido.
Una vez realizada la división del terreno, se analiza por separado cada tabla definida. En el ejemplo se realizó la nivelación por tablas, haciendo una divisoria en la parte media del terreno donde se puede ver claramente que hay un cambio de pendiente brusco, buscando reducir los volúmenes de corte y relleno y por ende el costo total de la nivelación (Figuras 5 y 6).
Resultados del cálculo del proyecto de nivelación
Después de realizar los cálculos respectivos con los datos de entrada, Terra proporciona con características detalladas información de todo el proyecto, agrupado en tres cejas: resultados nivelación, cortes y rellenos y reporte o memoria de cálculo. Las secciones que se describen en los resultados son: datos generales del proyecto, parámetros del diseño, coordenadas de la poligonal y resultados de la nivelación (generales, ajuste al plano, cortes y rellenos, volúmenes y costos).
En la sección de cortes y rellenos, se presenta una tabla que en cada nodo (i, j) contiene, las coordenadas del terreno natural, cota del plano proyecto y el corte o relleno en centímetros de cada punto de la cuadrícula. Finalmente, se presenta la sección del reporte de la memoria con el cálculo desglosado de los resultados: la solución matricial del sistema de ecuaciones para calcular el plano proyecto, el cálculo del área, centroide y cota centroidal, el proyecto inicial con la información topográfica y del plano proyecto, el cálculo de la ecuación del plano y las iteraciones para el ajuste del plano (Figura 6).
Cálculo y generación de isolíneas
Con Terra se pueden calcular y dibujar las curvas de nivel representativas de la cuadrícula analizada. Para ello se crean las isolíneas de dos formas distintas: la primera con los datos de campo obtenidos en el levantamiento topográfico por cuadrícula, la segunda mediante el procedimiento geoestadístico de interpolación espacial detallado en las secciones anteriores, con el cual se crea una cuadrícula con uno de dos métodos: el método de Kriging o el del inverso de la distancia. También se presenta un plano de cortes y rellenos y un diagrama de bloques de cortes y rellenos que sirve como guía al tractorista para el trabajo de la nivelación en campo (Figura 6).
Estimación de costos
El costo total de nivelación para el ejemplo de las figuras 4 y 5, cuyo resumen se muestra en el Cuadro 3, tuvo una reducción de 43%, de $65 114.14 a $ 36 987.43 considerando el costo/m3 de movimiento de tierra en $20.80, al dividir el terreno en dos tablas basado en la topografía del terreno natural, lo cual mejoró el diseño inicial sin divisiones, considerando solamente pendientes naturales promedio.
Cuadro 3 Comparación de volúmenes y costos para el ejemplo, aplicando la nivelación al terreno completo vs el predio dividido en dos tablas según topografía del terreno natural.
| No. Tablas | Área (ha) | Volumen (m3) | Volumen total (m3) | Costo ($) | ||||||
| Corte | Relleno | Corte | Relleno | Por hectárea | Total | |||||
| 1 | 21.371 | 146.48 | 115.987 | 3 130.4874 | 2 478.8105 | 3 046.776 | 65 114.137 | |||
| 2 | 21.371 | 1 778.242 | 1 377.892 | 36 987.43 | ||||||
| 11.215 | 57.144 | 45.057 | 640.859 | 505.309 | 1 188.597 | 13 329.875 | ||||
| 10.157 | 111.984 | 85.912 | 1 137.382 | 872.583 | 2 329.261 | 23 657.555 | ||||
El valor de 1 en núm. tablas, indica que el terreno no fue dividido.
Uno de los aspectos técnicos más importantes del presente trabajo es la generación del modelo digital de elevaciones mediante el método de Kriging, para la determinación de alturas de corte o relleno y cuantificación de los volúmenes de corte y relleno, donde se aplica el método de los cuatro puntos sugerido por el Departamento de Agricultura de Estados Unidos (1970) y que aplica para su estimación una cuadrícula cuya área de influencia sea un cuadrado. Para garantizar que el modelo de Terra, es técnicamente adecuado se comparó el DEM con un programa profesional como Surfer (Golden software, 2011), los modelos generados se muestran en la Figura 7.
Los modelos constan de 39 columnas y 73 hileras, con un total de 2847 datos de elevaciones Z, con medias: Zm- Terra= 20.00739094, Zm-SURFER = 20.00997930 y varianzas: S²- Terra= 0.00373831, S²-SURFER: 0.00388370, al practicarles una prueba estadística para medias diferentes y varianzas distintas se concluyó que los modelos anteriores con 97% de confiabilidad tienen valores de medias iguales y con un valor de 99% de confiabilidad tienen varianzas similares, lo que nos hace concluir que son prácticamente iguales.
Lo anterior, permite aseverar que el modelo generado por Terra, predice el comportamiento en forma adecuada y por tanto al generar la cuadrícula para estimar los valores de cortes y rellenos genera volúmenes congruentes que permiten ajustar el plano proyecto de la nivelación de forma adecuada. Con la aplicación de las tecnologías modernas para la nivelación de tierras, como es la nivelación con rayo láser, pareciera ser innecesario el desarrollo de un proyecto de nivelación.
Sin embargo, no es así; cualquier proyecto de nivelación, independientemente de la forma en que se tomen los datos de campo, necesita un proyecto ejecutivo con datos y resultados que incluyan una memoria de cálculo con la estimación de volúmenes de corte y relleno y por supuesto el costo de la nivelación, el plano con la topografía original y el terreno nivelado. Esto será el sustento del proyecto y el respaldo para que el usuario realice los trámites necesarios en caso de requerir apoyos económicos en instituciones que lo brindan al campo o bien en el pago del movimiento de tierras al tractorista que realice el trabajo, sin duda alguna.
Conclusiones
Terra 1.0, que se presenta, es un sistema de cómputo electrónico desarrollado como herramienta fundamental para auxiliar a los técnicos en la formulación y optimización de proyectos ejecutivos de nivelación de tierras. Incluye algoritmos para el cálculo del proyecto a partir de un levantamiento topográfico realizado en cuadrícula o con distribución variable.
El uso del método de Kriging se integró con éxito para el cálculo de las isolíneas del plano del terreno natural, del terreno nivelado, y la utilización de una cuadrícula levantada o generada que permite realizar el cálculo del volumen de cortes y rellenos. La metodología utilizada, implicó el resolver en forma detallada el ajuste de mínimos cuadrados ordinarios usando un método matricial para el cálculo del plano inicial e implementar un algoritmo eficiente para el ajuste del coeficiente CA de la ecuación del plano.
En el caso de la cuadrícula generada, se incluyen algoritmos para la determinación de semivariogramas y se propusieron ecuaciones algebraicas para el cálculo de los coeficientes C1 y C2 del variograma teórico, para el cálculo de los pesos en la interpolación espacial de los datos de la cuadrícula. En caso de la división en tablas o terrazas, se implementaron algoritmos para la separación de los puntos de campo según los polígonos divididos; se incluyen las rutinas necesarias para que el proyectista pueda generar un diseño técnicamente adecuado de la nivelación, con lo que se logra la disminución de los volúmenes totales de cortes y rellenos y por tanto se logra disminuir los costos de nivelación. El sistema Terra, se podrá utilizar con las nubes de puntos que se generan en levantamientos topográficos que se realizan con VANTs, conocidos como drones.
