Fundamentos de análisis geográfico con SEXTANTE - La Salle

332 CAPÍTULO 20. PROGRAMACIÓN DE MÓDULOS DE SEXTANTE
Igualmente, para comprobar si una celda contiene valor de sin datos, no compruebes directamente
el valor, sino utiliza el siguiente método
virtual bool is_NoData (int x, int y)
En ocasiones, puedes necesitar algún grid accesorio además de los que introduce el usuario
(por ejemplo, para usarlo como grid intermedio donde almacenar información temporal
de una operación. Para ello, puedes usar algunas de las diferentes formas de la función
API Create Grid() (no es un método de CGrid. Puesto que lo más habitual es que necesites un
grid de las mismas dimensiones y características que uno que ya tienes (uno de los introducidos
por el usuario), la forma más frecuente es la definida en el siguiente prototipo
CGrid* API_Create_Grid (const CGrid &Grid);
20.4. Manejando capas vectoriales. La clase CShapes
Si se conoce bien la jerarquía de las capas vectoriales tal y como en su momento fue
explicada, manejar internamente capas vectoriales es sencillo. Para acceder a cada uno de los
elementos de una capa vectorial existen métodos apropiados de uso intuitivo, y el trabajo
tanto con las propias capas como con sus tablas de atributos no es en absoluto complejo.
Empezando por el nivel más alto de la jerarquía encontramos la clase CShapes, que almacena
una serie de objetos de la clase CShape, siguiendo esa jerarquía que ya conocemos. Existen
4 clases que heredan de CShape, y que sirven para los cuatro tipos de entidades vectoriales
que ya hemos visto: CShape Polygon, CShape Line y CShape Points. Todas ellas no dejan de
ser sino conjuntos de puntos, y para el manejo de éstos tienen una serie de funciones comunes.
Para conocer las coordenadas de el punto i de la parte j de una entidad (es decir, un objeto
de clase CShape), usa el método Get Point() en una llamada como la siguiente.
TGEO_Point pt = Get_Point(i,j);
La estructura TGEO Point simplemente tiene dos valores, x e y, que puedes consultar para
saber las coordenadas del punto en ambos ejes.
Para ver un ejemplo práctico de cómo crear entidades y cómo leer la información de otra
ya existente, el siguiente bloque de código toma una capa de polígonos (a partir de un puntero
pInput) y crea otra (en pOutput) igual. Creo que esta es una buena manera de ir aprendiendo
los métodos más útiles de estas clases. Por el momento, sólo prestamos atención a la geometría
y no a las tablas de atributos.
int i,j,k;
TGEO_Point Point;
CShape *pShape,*pShape2;
for(i=0; iGet_Count(); i++){
pShape = pInput->Get_Shape(i);
pShape2 = pOutput->Add_Shape();
for(int j=0; jGet_Part_Count(); j++){
for(int k=0; kGet_Point_Count(j); k++){
Point = pShape->Get_Point(k,j);
pShape2->Add_Point(Point.x,Point.y,j);
}
}
}