[Gvsig_usuarios] Re: Interseccion: Error en calculo de areas?

[Agustin Lobo]

Creo que he encontrado la razon
del problema, aunque parece que
continua habiendo un error. Me explico:
Cuando se hace una interseccion, los campos
AREA y PERIMETRO de la capa original
(en este caso, la "A") no se actualizan, sino
que continuan siendo los valores del poligono
correspondiente en la capa A. Quiza sea solo
cuestion de documentar bien esto de forma
que nadie mas se confunda.

He recalculado las AREAS y PERIMETROS creando
dos campos nuevos en la tabla del vector con los poligonos
de las intersecciones (AREAS2 y PERIMETROS2). Para ello
he utilizado la funcion area() y un procedimiento
equivalente al descrito en Course-gvSIG-1.1-es.pdf p.70 y
sigts. Y aqui continua habiendo un problema:
aunque la mayoria de los poligonos de cada cuadrado
suman los 4000000 m2 de rigor, hay algunos cuya suma
es superior, uno en concreto hasta 5036804 m2. He seleccionado
todos los poligonos correspondientes a ese cuadrado y
no hay ninguno repetido ni ninguno que sobresalga...,
lo que querria decir que hay algun error en el calculo
del area.

Vuelvo a repetir que puedo enviar la capa vectorial
para aclarar mas este asunto.

Saludos,

Agus

Agustin Lobo escribió:
> Tengo la impresion de que hay un error en el calculo
> de las areas resultantes de una interseccion.
> 
> Tengo dos capas de poligonos vectoriales en UTM31N ED50:
> A. Poligonos de un mapa de habitats
> B. Cuadrados de 2km x 2km. Supongamos que tengo solos dos, C1 y C2.
> 
> Interseccion: Resultado, el mapa de habitats extraido para los cuadrados 
> C1 y C2
> 
> Problema: la suma de las areas de cada poligono del mapa dentro de cada
> uno de los cuadrados no es igual y es distinta de 4000000. Sin embargo,
> si selecciono todos los poligonos correspondientes a un cuadrado se
> visualiza la seleccion de todo el cuadrado. Por tanto, el problema no
> es que falte algun poligono, sino que parece que las areas no esten bien 
> calculadas. De hecho, en algunos casos la suma es > 4000000
> 
> Si me indicais un servidor o una direccion, envio el vector resultante
> de la interseccion
> 
> Gracias
> 
> Agus

-- 
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: Agustin.Lobo en ija.csic.es
http://www.ija.csic.es/gt/obster