ISO 19157

Aplicación de la norma ISO 19157 en la cartografía de la BTN25 del CNIG

ISO 19157

Aplicación de la norma ISO 19157 en la cartografía BTN25 de CNIG

Harold Mercado Llanos

QGIS plugin


Resumen

La información geográfica es hoy más accesible para los usuarios, contamos con múltiples fuentes y productores, lo que permite que la sociedad tengo muchos datos disponibles, pero no siempre la información es adecuada para una finalidad específica, he aquí la importancia de las herramientas para la evaluación de calidad que deben ser igualmente divulgadas y de fácil acceso, para que la implementación de las normas y estándares, aporten al usuario una valoración sobre la confianza que puede tener sobre el producto cartográfico y que sea una herramienta que le permita estimar la usabilidad (adecuación al uso) y comparar si lo indicado por el productor está conforme a la calidad esperada del conjunto de datos y así pueda tomar mejores decisiones al conocer las fortalezas y debilidades de la información.


En la actualidad el modelo para la evaluación de la calidad de datos geográficos más extendido y que se encuentra vigente es el establecido por la norma ISO 19157, que a pesar del manto de dudas que pueda surgir o las recomendaciones que se quieran aportar, esta norma propone un conjunto de medidas normalizadas de calidad y los procesos básicos de evaluación que son herramientas útiles para obtener medidas cualitativas y cuantitativas sobre el conjunto de datos geográficos.


Como aporte a la discusión, este trabajo pretenden automatizar lo expuesto en la norma ISO 19157, creando una herramienta guiada, estructurada y de carácter gratuito, que recopila los conceptos básicos de la norma, con especial atención a las medidas y métodos definidos en esta, luego se aplica esta herramienta sobre la cartografía oficial del Instituto geográfico Nacional (IGN) de España, para los mapas vectoriales y Bases Cartográficas y Topográficas de la BTN25, donde necesariamente debemos conocer las especificaciones del producto para identificar la estructura y características de los datos geográficos que ofrecen a los usuarios.


La herramienta permite la evaluación de la consistencia lógica, donde se evalúa la consistencia conceptual, la consistencia de formato, la consistencia de dominio y la consistencia topológica, se evalua la compleción (comisión y omisión), la exactitud posicional, la exactitud temporal y la usabilidad, generando para cada una de las evaluaciones un reporte de los resultados asociados a las medidas recomendadas por la norma ISO 19157 en el anexo D.


Esta creada como un complemento (plugin) para Qgis, donde el usuario, mediante la interfase grafica podrá cargar, seleccionar y configurar los datos sobre los que se desea implementar la evaluación y cuenta con ventanas de ayuda para comprender los pasos requeridos en la evaluación de todo el conjunto de datos de una hoja de la BTN25.


Descripción del plugin

El plugin de Qgis sigue el siguiente esquema de desarrollo:


Esquema de desarrollo


Se hace una descripción de las diferentes pestañas con las que cuenta el plugin.

Menu inicial:


menu inicial

Proyecto nuevo:

El plugin requiere una carpeta vacía para crear el proyecto de evaluación de calidad, el sistema verifica esta condición y genera los archivos requeridos en la ubicación local indicada por el usuario. al iniciar un nuevo proyecto el plugin genera los siguientes archivos y directorios:

carpetas

Directorio CD donde se almacenan las CDR que deben ser reproyectados.

Directorio Resultados donde se almacenan los resultados de las evaluaciones de calidad en formato HTML.

Fichero cde.txt Archivo plano donde se relacionan las ubicaciones de los CDE y CDR cargados en el plugin.

Fichero iso19157.qgs Archivo de proyecto QGIS para la interface de la cartografía, los geoprocesos y la visualización de resultados.

Fichero Iso19157.txt Archivo de configuración del proyecto.

Se puede cambiar la ubicación de la carpeta completa del proyecto, pero no se recomienda la manipulación de los archivos o carpetas internas dado que puede llegar a generar perdida de información y errores en el comportamiento del plugin.

Se recomienda no cambiar la ubicación de los conjuntos de datos CDE y CDR, esto genera que el plugin no puede volver a localizar los datos de entrada.

Proyecto existente:

En el caso que ya se tenga iniciado un proyecto de evaluación, el plugin verifica que en la carpeta indicada por el usuario existan los archivos correspondientes a dicho proyecto. Adicionalmente si ya ha realizado alguna evaluación sobre el conjunto de datos se habilita la pestaña Resultados donde encontrará los enlaces correspondientes.

Si se crea un nuevo proyecto o se cargue un proyecto existente de manera satisactoria se genera el mensaje "Puede continuar o revisar el resultado de la evaluación", en este caso se puede pasar a la pestaña "Conjunto de Datos"


Conjunto de datos:


Conjunto de datos

Adicionar objeto:

Permite al usuario localizar en los directorios locales el conjunto de datos a evaluar CDE, solo acepta formato SHP.

Seleccione Conjunto de Datos Evaluación (CDE):

Permite seleccionar mediante un menu desplegable el CDE al que se le desea asociar un conjunto de datos de referencia CDR

Asociar CDR:

Una vez seleccionado el CDE en el menu desplegable, esta opción permite ubicar el archivo CDR en el directorio local, solo acepta formato SHP

Eliminar un CDE cargado:

Dando click sobre la fila del CDE se elimina el conjunto de datos de la evaluación

Validar CD

se valida que el CDE y CDR esten en el mismo sistema de coordenadas, que compartan la misma cobertura espacial y que tengan la misma topología, el resultado de esta validación se presenta en la pestaña "Validar CD"


Validar CD:


Validar CD

Se presenta el resultado de la validación del CDE y CDR se indica el numero de elementos cargados.


Sistema de coordenadas de referencia:

Verifica que el SCR del CDR sea igual que el SCR del CDE, en caso contrario se activa la opcion de "Reproyectar" que genera un nuevo CDR que se guarda en la carpeta del proyecto "CD" con el resultado de la reproyección

Cobertura espacial:

En el caso de la cartografía de la BTN25, se verifica la existencia del Shape "MAR_HOJ" en la carpeta de donde el usuario ha tomado los CDE, si este archivo no es encontrado se genera la cobertura con la extensión máxima de los CDE cargados, igualmente se genera la cobertura con la extensión máxima de los CDR, se registra el area de ambas coberturas y se estima el procentaje de cobertura compartida. Si no tienen cobertura compartida se pide al usuario que modifique los CDR.

Topología:

El plugin solo permite realizar la evaluación de calidad entre conjunto de datos de igual topología (punto-punto, linea-linea o polígono-polígono), en caso que se detecte que no corresponden al mismo tipo, de indica al usuario que debe eliminar el conjunto de datos o modificarlo para que se cumpla la condición.

Continuar:

Esta opción se activa una vez se cumpla el proceso de validación

Elementos de calidad:


Elementos de calidad

Se presentan los CDE cargados y se permite al usuario habilitar o deshabilitar los elementos de calidad a evaluar en cada conjunto de datos, en caso que el CDE no tenga un CDR asociado solo permitirá la evaluacion de la consistencia lógica, la calidad temporal y la usabilidad.


Iniciar:

Una vez el usuario habilite los elementos de calidad a evaluar se puede iniciar con el proceso en la pestaña de Evaluación


Norma ISO 19157

La evaluación de calidad del plugin esta basada en los parametros y medidas establecidos por la norma internacional ISO 19157 que se centra en lo que se denomina elementos cuantitativos de la calidad de los datos espaciales, es decir, aquellos que se pueden expresar de una forma numérica, se organiza en los siguientes elementos de calidad de los datos:



Evaluación de calidad

El plugin implementa la norma ISO 19157 tomando como referencia las medidas y métodos normalizados para cada uno de los elementos y subelementos de calidad, se recomienda seguir el orden de evaluación establecido por la norma ISO 19157 anexo I 2.3.


Para desarrollar el proceso de evaluación de calidad, se guien los pasos indicados por la norma ISO 19157 anexo E.3:


  1. Especificar las unidades de calidad de datos: El plugin en el menu elementos de calidad permite que el usuario especifique las unidades de calidad de datos que evaluará para cada objeto geográfico, esta selección debe estar asociada al tipo de objeto geográfico y a los requerimientos del usuario.

  2. Especificar las medidad de calidad: El plugin desarrolla las medidas de calidad indicadas en la norma ISO 19157 anexo D para cada una de las unidades de calidad seleccionadas por el usuario, no existe filtro previo, dado que se reportan los resultados para todas mas medidas.

  3. Especificar los procedimientos: La norma ISO 19157 determina diferentes métodos de evaluación de calidad de datos:

    El método de evaluación directa se basa en la inspección de los ítems del conjunto de datos, puede clasificarse en interna y externa, el método de evaluación directa interna utiliza solo datos que se encuentran en el propio conjunto de datos que se esta evaluando, la evaluacion directa externa requiere un conjunto de datos externa independientes al conjunto de datos evaluados.
    El método de evaluacion indirecta esta basado en el conocimiento y la experiencia externa sobre el conjunto de datos evaluado, puede ser subjetivo y pueden considerarse el linaje, el productor y el proposito del conjunto de datos.

    El plugin desarrolla diferentes enfoques de inspección para la obtención de la información en el método de evaluación directa:


    El resultado de evaluación de calidad reportado en el plugin incluye el método de evaluación para cada elemento de calidad y el enfoque de inspección.


  4. Determinar la salida de la evaluación de la calidad de datos: EL plugin realiza un informe de calidad independiente, de las medidas de calidad, los procedimientos y los métodos empleados

  5. Consistencia lógica:

    Este elemento de calidad requiere solo los CDE, se verifica la consistencia lógica del 100% de la cobertura de los datos y la evaluación se establece con las características indicadas en las especificaciones de la BTN25. de donde se obtiene el catálogo de objetos

    Para el desarrollo de este modulo se sigue el siguiente esquema:

    Esquema de desarrollo


    (continuación) Para los CDE tipo linea los procesos implementados son:

    Esquema de desarrollo


    (continuación) Para los CDE tipo polígono los procesos implementados son:

    Esquema de desarrollo


    Catálogo de objetos: Fichero en formato plano (txt) donde se relacionan los objetos representados en el CDE, describiendo su nombre, atributos, dominio y tipo de datos, se encuentra cargado en el plugin de manera predeterminada el catálogo de objetos de para la BTN25 con la estructura.

    Catálogo de objetos

    Los usuarios pueden generar el catálogo de objetos personalizado, solo deben cumplir con la siguiente estructura:


    Fichero plano, con extensión "txt", separado por comas ",".

    Encabezado "CODIGO,NOMBRE,TABLA,ATRIBUTO,DOMINIO,TIPO"

    Donde:

    Codigo: Caracteres alfanumericos que representan el CDE.

    Nombre: Nombre del CDE segun las especificaciones.

    Tabla: Nombre de la tabla del CDE

    Atributo: Nombre del campo del CDE

    Dominio: Valores que puede tomar el atributo, se escriben separados por "-", esta caracteristica puede ser nula, en ese caso se deja sin asignar.

    Tipo: Tipo de dato en el que se encuentra almacenado el atributo (campo) segun las especificaciones, debe asignarse uno de estos 3 tipos (String, Integer, Real). Si se define un atributo debe tener un tipo asignado.


    Cuando se tenga el fichero debe ser almacenado en la localización del plugin (por defecto C:\Users\_________\AppData\Roaming\QGIS\QGIS3\profiles\default\python\plugins\iso_19157) en la carpeta "Datos".


    No remueva este fichero de la ubicación dado que afecta el comportamiento del plugin

    Puede modificar el fichero, respetando los parametros antes indicados


    Resultados:

    La evaluación de la consistencia lógica genera para cada CDE, una capa de puntos con la representación de los errores encontrados, se emplea la siguiente simbologia

    CaracterísticaSimbología
    Contornos internos consistentes (poligonos)
    Ordenación de los vertices (capa rios)
    Pt/Vert Duplicados (puntos, lineas y poliginos)◻ - 🛆
    Puntos Superfluos (lineas y poligonos)O
    Bucles (lineas y poligonos)X
    Elementos unificados (lineas)
    Cruces y anclajes en exceso (lineas y poligonos)
    Anclajes por defecto (lineas)>

    Para las características "Entidades atributos duplicados", "Elementos solapados (lineas)" y "Contornos disjuntos (poligonos)", se crean capas adicionales

    Las medidas y métodos de la Norma ISO 19157 (anexo D) implementadas son:

    MedidaAnexo DTipoDescripción
    4Tabla_D.4EnteroNúmero total de duplicaciones exactas de instancias de objetos geográficos en el conjunto de datos
    9Tabla_D.9BooleanoVerdadero indica que un item cumple las reglas del esquema conceptual
    10Tabla_D.10EnteroRecuento de todos los ítem del conjunto de datos que no cumplen las reglas del esquema conceptual
    11Tabla_D.11EnteroNúmero total de superposiciones erróneas en los datos
    12Tabla_D.12RealNúmero de items del conjunto de datos que no cumplen las reglas del esquema conceptual en relación al supuesto número total de estos ítems en el conjunto de datos
    15Tabla_D.15BooleanoVerdadero indica que un ítem no es conforme con su dominio de valores
    16Tabla_D.16EnteroRecuento de todos los ítem del conjunto de datos que no son conformes con su dominio de valores
    18Tabla_D.18RealNúmero de datos que no son conformes con su dominio de valores en relación al número total de ítem
    119Tabla_D.19BooleanoVerdadero indica conflicto en la estructura física
    19Tabla_D.20EnteroRecuento de todos los ítem del conjunto de datos que están almacenados en conflicto con la estructura física del conjunto de datos
    20Tabla_D.21RealNúmero de ítems de datos que están almacenados en conflicto con la estructura física del conjunto de datos, dividido entre el número total de ítems
    23Tabla_D.24EnteroRecuento de ítems del conjunto de datos, dentro del parámetro de la tolerancia*, que no se han combinado debido a subtrazos.
    24Tabla_D.25EnteroRecuento de ítems del conjunto de datos, dentro del parámetro de la tolerancia*, que no se han combinado debido a sobretrazos.
    26Tabla_D.27EnteroRecuento de todos los ítems de los datos que ilegalmente se intersectan con ellos mismos.
    27Tabla_D.28EnteroRecuento de todos los ítems de los datos que ilegalmente se autosuperponen con ellos mismos

    Tomado de la Norma ISO 19157


    Las especificaciones de la BTN25, adicionalmente indica que se deben valorar:


    * Parámetro de tolerancia para la evaluación de la consistencia lógica, se define el parámetro de tolerancia mediante el limite de percepción visual, índice que asume el límite de percepción del ojo humano, para una vista normal entre 0,25 mm, con un error no superior al 0,2 mm, para la observación en condiciones normales, dado que en el proceso de medición de productos cartográficos se requiere dedicación y un gran numero de lecturas, se puede incrementar el valor de percepción normal, en un porcentaje de relajación (Pr), que el usuario define teniendo en cuenta las condiciones externas (ARIZA F), esto permite aumentar el valor de percepción visual debido a la fatiga que se pueda presentar dado que el proceso de medición debe ser continuo.

    Uno de los valores utilizados para Pr es el 25%, considerando un muestreo continuo de 20 a 30 mediciones continuas, por tanto la tolerancia para el límite de percepción visual esta determinado por:

    Tolerancia (metros) = (0,0002 ) . (Pr) . (Me)

    Donde: Me es el modulo de la escala del producto cartográfico.


    En el caso de la BTN25 el limite de percepción visual seria:

    Tolerancia (metros) = 0,0002 x 1,25 x 25.000 = 6,25 m

    Complesión:

    Este elemento de calidad requiere que los CDE tengan asociado con un CDR.

    Para el desarrollo de este modulo se sigue el siguiente esquema:

    Esquema de desarrollo


    Extensión espacial de la evaluación: El usuario debe definir el area de evaluación, el plugin presenta 3 opciones:


    El plugin permite comparar el CDE, con todas las entidades del CDR o realizar un filtro del CDR mediante un valor de atributo específico. Este parámetro debe ser definido por el usuario para cada uno de los CDE a evaluar.


    Resultados:

    Las medidas y métodos de la Norma ISO 19157 (anexo D) implementadas son:

    MedidaAnexo DTipoDescripción
    1Tabla_D.1BooleanoVerdadero indica que el ítem es excedente
    2Tabla_D.2EnteroNúmero de ítems del conjunto o muestra de datos que no deberian haber estado presentes
    3Tabla_D.3RealNúmero de ítems excedentes en el conjunto o muestra de datos en relación al número de ítems que deberían haber estado presentes
    5Tabla_D.5BooleanoVerdadero indica que el ítem es omitido
    6Tabla_D.6EnteroRecuento de todos los ítems que deberían haber estado ene el conjunto o muestra de datos y estan omitidos
    7Tabla_D.7RealNúmero de ítems omitidos en el conjunto o muestra de datos en relación al número de ítems que deberían haber estado presentes

    Tomado de la Norma ISO 19157


    Exactitud posicional:

    Este elemento de calidad requiere que los CDE tengan asociado con un CDR.

    Para el desarrollo de este modulo se sigue el siguiente esquema:

    Esquema de desarrollo


    Extensión espacial de la evaluación: El usuario debe definir el area de evaluación, el plugin presenta 3 opciones:


    El plugin solicita al usuario ingrese el número de puntos a muestrear, como minimo se deben indicar 25 y un máximo de 500 puntos.


    El usuario debe indicar la escala del CD, para calcular las tolerancias, siguiendo los siguientes parámetros:


    Parámetro de tolerancia: (metros)

    limite de percepción visual: índice que asume el límite de percepción del ojo humano, para una vista normal entre 0,25 mm, con un error no superior al 0,2 mm, para la observación en condiciones normales, dado que en el proceso de medición de productos cartográficos se requiere dedicación y un gran numero de lecturas, se puede incrementar el valor de percepción normal, en un porcentaje de relajación (Pr), que el usuario define teniendo en cuenta las condiciones externas (ARIZA F), esto permite aumentar el valor de percepción visual debido a la fatiga que se pueda presentar dado que el proceso de medición debe ser continuo.

    Uno de los valores utilizados para Pr es el 25%, considerando un muestreo continuo de 20 a 30 mediciones continuas, por tanto la tolerancia para el límite de percepción visual esta determinado por:

    Tolerancia (metros) = (0,0002 ) . (Pr) . (Me)


    Error posicionamiento de contornos - EPC: Se basa en la capacidad de posicionar o adquirir información de un producto cartográfico, dadas las características del terreno representado, valorando la incertidumbre en el desplazamiento de un punto debido a la escala de representación. Se utiliza el EPC para determinar la precisión en la medición sobre cartografía, teniendo en cuenta la escala de entrega o para calcular la escala de impresión para “disimular” el desplazamiento del punto representado o las coordenadas dadas. Se utiliza un valor de ponderación de acuerdo al tipo de terreno, denominado Error de contorno – Ec, los valores usuales son: 0.0005 metros para zonas planas y 0.0007 para zonas quebradas o de alta pendiente. La expresión para calcular EPC será:

    EPC = Ec . Me

    Estándar de precisión para datos geoespaciales digitales (ASPRS): El objetivo del estándar de precisión ASPRS para datos geoespaciales digitales reemplazan los estandares existentes para grandes mapas 1990 y las directrices ASPRS, para la precisión vertical Lidar Data 2004, con el nuevo estándar de precisión se aborda mejor la tecnología digital.
    El estandar ASPRS aplicado para la precisión horizontal para mapas planimétricos es aplicable a cualquier escala cartográfica, la implementación se organiza en clases (I, II, III, ... N), donde el RMSE en X o Y (cm) = 0.0125 veces el factor de escala.

    RMSEx o RMSEy = ((Clase) x 0.0125 x (Modulo de la escala))

    (RMSEx / 100) (RMSEy / 100) para indicar la medida en metros




    Precisión Horizontal (radial) al 95% confianza = RMSEr x 1.7308

    Tolerancia del usuario: Se basa en la experiencia del usuario y el conocimiento que tenga de la cartografía evaluada


    Detección de errores groseros (Outliers):

    Se pude detectar la presencia de errores groseros aplicando test estadísticos sobre el residual (R), donde se determina un potencial error grosero se presentan cuando R es mayor que el valor que se calcula de la desvición estandar de las observaciones multiplicado por un factor (M), que depende del tamaño de la muestra. Se puede calcular para 1 dimensión (Test Lineal) o para 2 dimensiones (Test Circular).


    Test Lineal:

    Residual, R =

    R es un potencial cuando:

    y para el calculo de donde v = n - 1 grados de libertad.



    Test Circular:

    Residual, R =

    R es un potencial cuando:

    y para el calculo de donde v = n - 1 grados de libertad.



    Test de significancia del sesgo calculado:

    Para determinar si la presencia de sesgo es significativa, se puede probar contra la distribución t, el sesgo debe considerarse significativo con el nivel de confianza del 90% si el valor cero no se encuentra en el rango definido por:


    donde:


    Se aplica de manera independiente para cada posicion (este, norte, cota):



    Resultados:

    Las medidas y métodos de la Norma ISO 19157 (anexo D) implementadas son:

    MedidaAnexo DTipoDescripción
    28Tabla_D.29MedidaValor medio de las incertidumbres posicionales para un conjunto de posiciones donde la incertidumbre posicional se define como la distancia entre la posición medida y la que se considera como verdadera

    la media de las incertidumbres posicionales de las posiciones horizontales absolutas o externas se calculan como:

    (esta medida de calidad es diferente a la desviación típica).
    128Tabla_D.30MedidaPara un número de puntos (N), se ofrecen las posiciones medidas como coordenadas xmi y ymi dependiendo de las dimensiones en las que se mide la posición del puntos. Se considera que un conjunto correspondiente de coordenadas xti y yti representa las posiciones verdaderas. las desviaciones y sesgos se calculan como:

    sesgo:

    29Tabla_D.31MedidaPara un número de puntos (N), se ofrecen las posiciones medidas como coordenadas xmi y ymi dependiendo de las dimensiones en las que se mide la posición del punto. Se considera que un conjunto correspondiente de coordenadas xti y yti representa las posiciones verdaderas. Todas las incertidumbres posicionales por encima de un umbral predefinido emáx son eliminadas del conjunto. las incertidumbres posicionales se calculan como:

    el cálculo de ei se define en la medida de la calidad "valor medio de las incertidumbres posicionales" en una, dos y tres dimensiones.
    para el número de errores restantes (NR), la la media de las posiciones horizontales absolutas se calcula como:

    30Tabla_D.32Entero Para un número de puntos (N), se ofrecen las posiciones medidas como coordenadas xmi y ymi dependiendo de las dimensiones en las que se mide la posición del punto. Se considera que un conjunto correspondiente de coordenadas xti y yti representa las posiciones verdaderas. El cálculo de ei se define en la medida de calidad "valor medio de las incertidumbres posicionales" en una, dos y tres dimensiones.
    se considera como error todas las incertidumbres posicionales por encima del umbral predefinido.
    31Tabla_D.33RealPara un número de puntos (N), se ofrecen las posiciones medidas como coordenadas xmi y ymi dependiendo de las dimensiones en las que se mide la posición del punto. Se considera que un conjunto correspondiente de coordenadas xti y yti representa las posiciones verdaderas. El cálculo de ei se define en la medida de calidad "valor medio de las incertidumbres posicionales" en una, dos y tres dimensiones.
    se considera como error todas las incertidumbres posicionales por encima del umbral predefinido.

    El número de errores se ofrece en relación al número total de puntos medidos.
    33Tabla_D.35MedidaError lineal probable, LEP (Linear Error probable)
    LE50 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    34Tabla_D.36MedidaError lineal tpico, SLE (Standard Linear Error)
    LE68.3 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    35Tabla_D.37MedidaExactitud lineal al 90% de nivel de significación, LMAS (Linear Map Accuracy Standard)
    LE90 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    36Tabla_D.38MedidaExactitud lineal al 95% de nivel de significación, LMAS (Linear Map Accuracy Standard)
    LE95 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    37Tabla_D.39MedidaExactitud lineal al 99% de nivel de significación, LMAS (Linear Map Accuracy Standard)
    LE99 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    38Tabla_D.40MedidaExactitud lineal casi cierto
    LE99.8 (Ver tabla G.2 de la norma ISO 19157 Anexo G.3.2)
    39Tabla_D.41MedidaError cuadrático medio, ECM (RMSE, Root Mean Square Error)
    El valor verdadero de un observable Z se conocen como xt
    de este modo del estimador

    se obtiene el error cuadrático medio ECM = σz
    41Tabla_D.43MedidaError lineal absoluto al 90% de nivel de significación de datos verticales con sesgo (Alternativa 2)
    Se calcula la comparación de los datos (fuente) y el control (referencia) de la siguiente manera:
    1. 1. Se calcula el error absoluto en la dimensión vertical en cada punto:
    2. 2. Se calcula el error medio vertical:
    3. 3. Se calcula la desviación típica de los errores verticales:
    4. 4. Se calcula el ratio del valor absoluto del error medio entre la desviación típica:
    5. 5. Si ratio > 1,4 entonces se calcula k = 1.2515
    6. 6. si ratio <= 1,4 entonces se calcula k, el ratio del sesgo vertical entre la desviación típica de las altitudes, usando un ajuste polinómico cúbico mediante los valores tabulados definidos en Handbook of Tables for Probability and Statistics
    7. 7. Se calcula LE90 para la fuente:
    8. 8. Se calcula el LE90 absoluto:
    42Tabla_D.44MedidaDesviación típica circular, error puntual de Helmert, CSE (Circular Standard Error)
    CE39.4 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.3)
    43Tabla_D.45MedidaError circular probable, CEP (Circular Error probable)
    CE50 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.3)
    44Tabla_D.46MedidaError circular estandar, CMAS (Circular Map Accuracy Standard)
    CE90 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.3)
    45Tabla_D.47MedidaExactitud de navegación
    CE95 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.3)
    46Tabla_D.48MedidaError circular casi cierto, CNCE (Circular Near Certainty Error)
    CE99.8 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.3)
    47Tabla_D.49MedidaError cuadrático medio planimétrico, ECMP (RMSEP, Root Mean Square Error of Planimetry)
    El valor verdadero de las coordenadas observadas X e Y se conocen como xt y yt
    de este modo del estimador

    se obtiene el error cuadrático medio planimétrico ECMP = σ
    48Tabla_D.50MedidaError circular absoluto al 90% de nivel de significación de datos con sesgo (OTAN)
    Se calcula la comparación de los datos (fuente) y el control (referencia) de la siguiente manera:
    1. 1. Se calcula el error absoluto en la dimensión horizontal en cada punto y en cada coordenada Xt y Yt
    2. 2. Se calcula el error medio horizontal de cada coordenada
    3. 3. Se calcula la desviación típica circular de las diferencias que se han medido entre el producto evaluado y la fuente de referencia:
    4. 4. Se calcula la desviación típica circular de los errores en la fuente de referencia: σR
    5. 5. Se calcula la desviación típica circular de los errores en el producto evaluado:
    6. 6. Se calcula el error circular absoluto al 90% de nivel de confianza de datos con sesgo (CMAS)
    49Tabla_D.51MedidaError circular absoluto al 90% de nivel de significación de datos con sesgo (ACE Absolute circular error)
    Se calcula la comparación de los datos (fuente) y el control (referencia) de la siguiente manera:
    1. 1. Se calcula el error absoluto en la dimensión horizontal en cada punto
    2. 2. Se calcula el error medio horizontal:
    3. 3. Se calcula la desviación típica de los errores horizontales:
    4. 4. Se calcula el ratio del valor absoluto del error medio entre la desviación típica:
    5. 5. Si ratio > 1,4 entonces k = 1,2815
    6. 6. Si ratio <= 1,4 entonces se calcula k, el ratio de la media entre la desviación típica, usando un ajuste polinómico cúbico mediante los valores tabulados definidos en Handbook of Tables for Probability and Statistics
    7. 7. Se calcula CE90 para la fuente:
    8. 8. Se calcula el CE90 absoluto:
    52Tabla_D.54MedidaError vertical relativo Rel LE90
    Se calcula la comparación de los datos (fuente) y el control (referencia) de la siguiente manera:
    1. 1. Se determinan todas las combinaciones de parejas posibles entre los puntos muestreados:
      Las combinaciones de parejas será:
    2. 2. Se calcula el error vertical en cada punto:
    3. 3. Se calcula el error relativo vertical para cada combinación de parejas de puntos:
    4. 4. Se calcula la desviación típica relativa vertical:
    5. 5. se calcula el error vertical relativo pasando el nivel de confianza al 90%
      Rel LE90 = 1,645.μZ rel
    53Tabla_D.55MedidaError horizontal relativo Rel CE90
    Se calcula la comparación de los datos (fuente) y el control (referencia) de la siguiente manera:
    1. 1. Se determinan todas las combinaciones de parejas posibles entre los puntos muestreados:
      Las combinaciones de parejas será:
    2. 2. Se calcula el error en las dimensiones X y Y en cada punto:
    3. 3. Se calcula el error relativo en X y Y para cada combinacion de parejas de puntos:
    4. 4. Se calcula la desviación típica relativa en cada eje:
    5. 5. se calcula la deviación relativa horizontal
    6. 6. se calcula el error horizontal relativo pasando el nivel de significación al 90%
      Rel CE90 = 2,146.μH rel

    Tomado de la Norma ISO 19157


    Tabla G.2 - Relación entre los cuantiles de la distribución normal y el nivel de significación

    Probabilidad PCuantilMedida básicaNombreTipo de valor
    P = 50%μ50% = 0,6745μ50%.σZLE50Medida
    P = 68,3%μ68,3% = 1μ68,3%.σZLE68.3Medida
    P = 90%μ90% = 1,645μ90%.σZLE90Medida
    P = 95%μ95% = 1,960μ95%.σZLE95Medida
    P = 99%μ99% = 2,576μ99%.σZLE99Medida
    P = 99,8%μ99,8% = 3μ99,8%.σZLE99.8Medida

    Tomado de la Norma ISO 19157


    Tabla G.5 - Relación entre la probabilidad P y el radio correspondiente del área circular

    Probabilidad PMedida básica de la calidad de datosNombreTipo de valor
    P = 39,4%CE39.4Medida
    P = 50%CE50Medida
    P = 90%CE90Medida
    P = 95%CE95Medida
    P = 99,8%CE99.8Medida

    Tomado de la Norma ISO 19157


    Se presenta un informe detallado de las medidas de incertidumbre de los puntos muestreados, donde se presentan las coordernadas y las diferencias posicionales para el "error norte" y el "error este"

    Se calculan y visualizan los errores o incertidumbres posicionales, mediante el diferentes estimadores de la estadistica circular


    donde:
    error norte (Yf - Yi)
    error norte (Xf - Xi)

    donde: f = referencia i = evaluado


    Acimut medio:

    donde:

    Modulo medio:

    donde:

    Varianza circular:

    Desviación estandar circular:

    Asimetría (skewness o sesgo):

    donde:

    Curtosis (o elevación):

    donde:

    Dispersión circular:

    donde:

    Desviación estandar angular:

    Desviación angular media (Batschelet, 1981):

    Parámetro Von Mises (aproximaciones):





    Calidad temporal:

    Este elemento de calidad se aplica solo sobre el CDE

    Para el desarrollo de este modulo se sigue el siguiente esquema:

    Esquema de desarrollo


    Resultados:

    Las medidas y métodos de la Norma ISO 19157 (anexo D) implementadas son:

    MedidaAnexo DTipoDescripción
    54Tabla_D.56MedidaExactitud temporal al 68,3% de nivel de significación
    LE68.3 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)
    55Tabla_D.57MedidaExactitud temporal al 50% de nivel de significación
    LE50 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)
    56Tabla_D.58MedidaExactitud temporal al 90% de nivel de significación
    LE90 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)
    57Tabla_D.59MedidaExactitud temporal al 95% de nivel de significación
    LE95 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)
    58Tabla_D.60MedidaExactitud temporal al 99% de nivel de significación
    LE99 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)
    59Tabla_D.61MedidaExactitud temporal al 99.8% de nivel de significación
    LE99.8 (Ver tabla G.5 de la norma ISO 19157 Anexo G.3.2)

    Tomado de la Norma ISO 19157


    Tabla G.2 - Relación entre los cuantiles de la distribución normal y el nivel de significación

    Probabilidad PCuantilMedida básicaNombreTipo de valor
    P = 50%μ50% = 0,6745μ50%.σZLE50Medida
    P = 68,3%μ68,3% = 1μ68,3%.σZLE68.3Medida
    P = 90%μ90% = 1,645μ90%.σZLE90Medida
    P = 95%μ95% = 1,960μ95%.σZLE95Medida
    P = 99%μ99% = 2,576μ99%.σZLE99Medida
    P = 99,8%μ99,8% = 3μ99,8%.σZLE99.8Medida

    Tomado de la Norma ISO 19157


    Exactitud temática:

    Usabilidad: