Más

¿Cómo contar polígonos en entidades multiparte?

¿Cómo contar polígonos en entidades multiparte?


Estoy tratando de contar el número de polígonos en una entidad multiparte. Mi objetivo es tener un atributo de 'recuento' para cuantificar la fragmentación de parcelas de tierra multiparte.

He mirado la calculadora de campo y los complementos adecuados, pero no veo nada que pueda ayudar. Hay muchas formas de contar puntos, líneas y vértices, pero no polígonos.


Sé que está buscando una respuesta que use QGIS, pero si está dispuesto a usar las herramientas de análisis geoespacial GIS Whitebox de código abierto en su lugar, aquí hay una respuesta completa que tomará cualquier shapefile y creará un nuevo campo de atributo en su tabla que contiene el número de partes de características. Simplemente necesita abrir la ventana de Scripter en Whitebox, cambiar el lenguaje de script a Groovy y pegar el siguiente código en:

import whitebox.geospatialfiles.ShapeFile import whitebox.geospatialfiles.shapefile. * import whitebox.geospatialfiles.shapefile.attributes. * def inputFileName = "/MY_DIRECTORY/MY_FILE_NAME.shp" ) DBFField field = new DBFField () field.setName ("NUM_PARTS") field.setDataType (DBFField.DBFDataType.NUMERIC) field.setFieldLength (5) field.setDecimalCount (0) table.addField (field) for (int i = 0 ; i 

Deberá actualizar el valor de la variable inputFileName cerca de la parte superior para apuntar a su archivo. También podrías usar Python para lograr lo mismo, pero escribí esto en Groovy. Si lo desea, también puedo hacer que el script excluya partes de los agujeros si está trabajando con polígonos.

A continuación, se muestra un ejemplo de un mapa de un país del mundo que se ha renderizado para mostrar la cantidad de partes que incluye la característica poligonal de cada país:


Cuatro años después, parece que QGIS ha integrado esta función en la calculadora de campo. La sintaxis de la expresión es:

num_geometries ($ geometría)


Si QGIS es la única herramienta disponible El siguiente código de Python en la Consola de Python podría funcionar:

layer = iface.activeLayer () para x en layer.getFeatures (): geometry = x.geometry (). asMultiPolygon () print len ​​(geometría)

Por supuesto, esta no es una solución completa, solo imprimirá un montón de números sin indicar a qué característica pertenecen, pero es un comienzo. También puede hacer que imprima el campo 'NOMBRE' o algún otro atributo de clave preferido. O puede modificar para escribir el entero len (geometría) en algún otro campo de la capa.


Otra solución que no es de QGIS es utilizar OpenJUMP. Tiene una herramienta especial en Herramientas - Estadísticas - Estadísticas de funciones. La herramienta crea una nueva capa con estos atributos:

  • nPts (número de vértices)
  • nHoles (número de agujeros)
  • nComponents (número de componentes)
  • área
  • largo
  • tipo (tipo ot geometría)

Si su objetivo es contar el número de polígonos en multipolígonos, esta herramienta debería funcionar bien. Sin embargo, si desea contar cuántos polígonos hay en GeometryCollections que contienen también otro tipo de geometrías, primero debe descomponer las colecciones.


Puede intentar disociar varios polígonos en un polígono único y luego abrir la tabla de atributos y observar el número de líneas (que representa el número de polígonos).

Para separar multipolígono:Vector-> Herramientas de geometría-> Multiparte a partes individuales


Tenga en cuenta que los datos digitales deben representar con precisión la posición legal del título a registrar.

Cuando se envía una solicitud de registro voluntario que utiliza datos digitales como alternativa a los planes en papel, la solicitud debe incluir un .pdf para ayudar con el mapeo de los sujetos. En estas circunstancias, no proporcionar el .pdf resultará en el rechazo de la solicitud.

Cuando una aplicación utiliza planos / escrituras en papel para definir la extensión, pero también se envían datos digitales, sería beneficioso que también se envíe un .pdf del diseño para usarlo como representación visual. El .pdf puede tener una escala menor que la requerida para los planes en papel, ya que es puramente una representación visual de los sujetos que se registran. Sin embargo, no proporcionar un .pdf en tales circunstancias no resultará en el rechazo de la solicitud.

Cuando se envíen datos digitales como parte de una solicitud, deben enviarse por correo electrónico a d[email protected] Cabe señalar en la sección de información adicional del formulario de solicitud que los datos digitales se envían por correo electrónico y se debe proporcionar la referencia del archivo.


Sintaxis

Las entidades poligonales que se agregarán. Si se trata de una capa que hace referencia a una representación y hay anulaciones de formas presentes en las entidades de entrada, las formas anuladas, no las formas de entidades, se considerarán en el procesamiento de agregación.

La clase de entidad de salida que se creará.

La distancia que se debe satisfacer entre los límites del polígono para que se produzca la agregación. Se debe especificar una distancia y debe ser mayor que cero. Puede elegir una unidad preferida, la predeterminada es la unidad de funciones.

El área mínima para retener un polígono agregado. El valor predeterminado es cero, es decir, para mantener todos los polígonos. Puede especificar una unidad preferida, la predeterminada es la unidad de funciones.

El tamaño mínimo de un agujero poligonal que se conservará. El valor predeterminado es cero, es decir, para mantener todos los huecos del polígono. Puede especificar una unidad preferida, la predeterminada es la unidad de funciones.

Especifica la característica de las entidades de salida al construir los límites agregados.

  • NON_ORTHOGONAL: se crearán entidades de salida con forma orgánica. Esto es adecuado para características naturales, como vegetación o polígonos del suelo. Este es el predeterminado.
  • ORTOGONAL: se crearán entidades de salida de forma ortogonal. Esta opción es adecuada para preservar la característica geométrica de los elementos de entrada antropogénicos, como las huellas de los edificios.

Las capas que contienen las entidades de línea o polígono que son barreras de agregación para las entidades de entrada. Las funciones no se agregarán entre las funciones de barrera. Se ignorarán las entidades de barrera que estén en conflicto geométrico con las entidades de entrada.

Una tabla de relaciones de uno a varios que vincula los polígonos agregados a sus entidades poligonales de origen. Esta tabla contiene dos campos, OUTPUT_FID e INPUT_FID, que almacenan los ID de funciones agregados y sus ID de funciones de origen, respectivamente. Utilice esta tabla para derivar los atributos necesarios para las entidades de salida a partir de sus entidades de origen. El nombre predeterminado para esta tabla es el nombre de la clase de entidad de salida, junto con _tbl. La ruta predeterminada es la misma que la clase de entidad de salida. No se crea ninguna tabla cuando este parámetro se deja en blanco.


Diagrama (izquierda): Polígono hidrográfico antes del proceso Split Narrow Polygons. Diagrama (derecha): Después del proceso Split Narrow Polygons, el polígono hidrográfico se divide en secciones donde el ancho es menor a 15 metros.

Esta herramienta está diseñada para identificar porciones de entidades poligonales que son demasiado estrechas para ser visualmente distinguibles a una escala específica. El ancho mínimo es una indicación de la distancia más corta entre los bordes del polígono que es sensible para mostrar en la escala de salida.

Los campos se agregarán a la clase de entidad de salida que contiene información sobre cada polígono dividido.

  • split_type: identifica qué tipo de área se está identificando.
    • 0: identificado como un área estrecha
    • 1: no identificado como un área estrecha
    • 2: se identifica como un área angosta pero cae por debajo de la longitud mínima
    • 3: no se identifica como un área angosta pero cae por debajo de la longitud mínima
    • 4 — Característica cónica
    • split_length: la longitud del segmento que se está dividiendo.

    El parámetro Longitud mínima es la distancia más corta que debería existir entre polígonos y o la longitud más corta que debería existir para un polígono individual mostrado a escala. Por ejemplo, si decidió reemplazar los polígonos split_type 2 (polígonos estrechos más cortos que la longitud mínima) con una línea central y mostrar todos los polígonos split_type 3 (polígonos anchos más cortos que la longitud mínima), esto podría resultar en una serie de polígonos conectados con segmentos de línea cortos. Dependiendo de los resultados deseados, puede optar por mantener los polígonos split_type 2 como visibles o almacenar los polígonos en búfer para asegurarse de que cumplen con el ancho mínimo.

    El parámetro Taper Length se proporciona para polígonos hidrográficos. Cuando se elimina una sección estrecha de un polígono de río, normalmente se mostrará como la línea central del río a escalas más pequeñas para mantener la conectividad visual de la red hidrográfica. Se puede crear un extremo puntiagudo o ahusado donde el polígono se estrecha para garantizar una transición suave desde el polígono hasta la línea central que se mostrará. Si se proporciona una longitud cónica de 0, no se crearán polígonos con split_type de 4. Todos los demás valores de split_type se crean con extremos rectos.

    La integridad de los resultados de esta herramienta se basa en la integridad topológica de las entradas. Tome nota de los siguientes requisitos y sugerencias de datos de entrada:

    • Geometría vacía o nula: las entidades de entrada deben constar de geometrías válidas. Si es necesario, use la herramienta Reparar geometría defectuosa para reparar estas características.
    • Operaciones de una sola parte: las operaciones de entrada no deben contener operaciones de varias partes. Utilice la herramienta Varias partes a una sola parte o cree una topología con una regla de línea Debe ser una sola parte para convertir entidades en una sola parte.
    • Vértices: si no hay suficientes vértices en un polígono, la calidad de los resultados puede verse comprometida. Se recomienda ejecutar la herramienta Densificar antes de ejecutar esta herramienta si las entidades poligonales tienen una pequeña cantidad de vértices.
    • Segmentos no lineales: los segmentos de curva de Bezier, arco circular y arco elíptico pueden comprometer la calidad de los resultados. Se recomienda ejecutar la herramienta Densificar antes de ejecutar esta herramienta si las entidades poligonales tienen segmentos no lineales, ya que esto agregará vértices a lo largo de las curvas.

    La clase de entidad de salida no contendrá ningún atributo geográfico de las entidades de entrada. Se agregará un campo llamado FID a la clase de entidad de salida. Este campo se completa con el ObjectID de la entidad desde la que se originó el polígono dividido. El campo FID se puede utilizar para vincular los polígonos divididos a sus polígonos de origen.

    Para obtener los mejores resultados, considere ejecutar la herramienta Eliminar parte de polígono para eliminar pequeños agujeros en las entidades de polígono antes de ejecutar Dividir polígonos estrechos. Por ejemplo, eliminar huecos en el medio de polígonos de ríos donde existen polígonos de islas. Si existe un agujero en el polígono, la herramienta considerará la distancia desde el borde del polígono hasta el borde del agujero. Esto podría dar como resultado que el polígono del río en cada lado de la isla se considere estrecho aunque el ancho total del río de borde a borde sea mayor que el ancho mínimo. Rellenar pequeños agujeros antes de ejecutar Dividir polígonos estrechos permitirá que la herramienta ignore la ubicación donde existía un agujero anteriormente y mida el ancho desde los bordes del polígono.


    ¿Cómo contar polígonos en entidades multiparte? - Sistemas de Información Geográfica

    Los shapefiles contienen objetos de datos geográficos en forma de puntos, líneas o polígonos con atributos informativos asociados que describen las características geométricas. Shapefiles es un formato estándar de la industria GIS desarrollado por el Instituto de Investigación de Sistemas Ambientales (ESRI) (consulte http://www.esri.com/). Los shapefiles contienen datos vectoriales geoespaciales y son el tipo principal de datos de entrada para SA Vector Tools.

    Cada shapefile almacena datos para una sola clase de entidad (es decir, un shapefile contiene información para puntos, líneas o polígonos). Las formas de puntos son entidades de una sola coordenada, como chimeneas, escuelas o puertos. Las formas de las líneas pueden ser líneas continuas, como carreteras, y también pueden ser polilíneas, como las ramas o afluentes de un río. Las formas de polígono pueden ser áreas simples o áreas de varias partes como estados y países (por ejemplo, el estado de Michigan es una forma de varias partes). Los shapefiles no almacenan información topológica porque no contienen asociaciones entre puntos, líneas y polígonos.

    Un shapefile consiste en un grupo de archivos que contienen las coordenadas de las formas, atributos asociados y, opcionalmente, información de proyección de mapas. Cada shapefile contiene tres archivos obligatorios:

    .shp & # 8211 coordenadas de formas o características (puntos, líneas o polígonos)

    .shx & # 8211 dato de índice posicional de forma, para respaldar la búsqueda rápida a través de los datos contenidos en el archivo de forma

    .dbf & # 8211 datos de atributos (por ejemplo, recuentos de población, clases de carreteras, capacidades del aeropuerto) para cada forma en formato dBASE III.

    Un shapefile también puede contener muchos archivos opcionales, cuyos detalles no son importantes para esta discusión.

    2.2 Suplentes espaciales

    Los sustitutos espaciales se utilizan con mayor frecuencia para mapear datos de emisiones a nivel de condado en las celdas de la cuadrícula rectangular utilizadas por un modelo de calidad del aire euleriano, como el sistema de modelado de calidad del aire comunitario multiescala (CMAQ) (consulte www.cmascenter.org). Un sustituto espacial es un valor mayor que cero y menor o igual a uno que especifica la fracción de las emisiones en un área (generalmente un condado) que debe asignarse a una celda de la cuadrícula del modelo en particular. Muchos tipos de inventarios de emisiones contienen datos resumidos por condado. Debido a que el área de un condado dado puede caer en varias celdas de la cuadrícula del modelo, se deben usar sustitutos espaciales para indicar la fracción de las emisiones para ese condado que se asignará a cada celda de la cuadrícula que se superponga al condado. Actualmente, la EPA usa alrededor de 65 tipos diferentes de sustitutos que se calculan a partir de alrededor de 24 shapefiles diferentes. Los sustitutos espaciales normalmente se deben calcular para cada cuadrícula que se modelará, a menos que la cuadrícula utilizada sea un subconjunto completamente alineado de una cuadrícula para la que ya se encuentran disponibles sustitutos espaciales.

    Cuando se crean sustitutos, generalmente hay algún tipo de atributo geográfico que se utiliza para "ponderar" las emisiones del condado en celdas de la cuadrícula de una manera que es más específica que una simple distribución uniforme por el condado, que no reflejaría las diferencias en las emisiones entre ciudades. , áreas rurales, puertos, etc. Por ejemplo, los pesos podrían ser información sobre puntos que representan puertos (p. ej., atracaderos) o las millas recorridas por vehículos en líneas que representan la ubicación de carreteras o ferrocarriles o atributos de peso basados ​​en polígono como área geográfica , población, número de hogares y categorías de uso del suelo. Puede crear sustitutos espaciales con el Asignador espacial a partir de shapefiles basados ​​en puntos, líneas o polígonos que contienen información de atributos de peso. Un valor sustituto único (srg) para el condado C y la celda de cuadrícula GC se expresa como:

    donde srg = sustituto, C = condado y GC = celda de la cuadrícula. Tenga en cuenta que un sustituto no tiene que ser para un condado, podría ser para alguna otra región geográfica, como un estado, provincia, nación o tramo censal. Estos polígonos se conocen en el programa como "polígonos de datos".

    Se pueden crear tres tipos de sustitutos con este programa: basados ​​en polígonos, basados ​​en líneas y basados ​​en puntos. Los sustitutos basados ​​en polígono utilizan información de atributos que se basa en el área (por ejemplo, la población en un tramo censal). El valor sustituto se calcula como la relación entre el valor del atributo en la intersección del condado y la celda de la cuadrícula y el valor total del atributo en un "polígono de datos" específico (por ejemplo, condado, estado, tramo censal). Ejemplos de atributos de peso basados ​​en polígono son área, población, número de hogares y uso de la tierra. El numerador (es decir, el valor del atributo de peso en la intersección del condado C y la celda de la cuadrícula GC) y el denominador (es decir, el valor total del atributo en el condado u otro polígono de datos) se calculan de acuerdo con las siguientes ecuaciones:

    donde srg = suplente, C = condado (u otro polígono de datos), GC = celda de cuadrícula, WP = polígono de peso, int. = intersección de, i = polígono de peso, yn = el número de polígonos de peso.


    Registro publicado 2011-03-09
    Última modificación del registro 2021-05-06
    Estado del recurso en marcha

    Descripción del objeto

    Nombre del objeto: WHSE_FOREST_VEGETATION.PEST_INFESTATION_POLY

    Nombre corto: PST_IF_PL
    Comentarios: PEST_INFESTATION_POLY representa una infestación del sitio como un área (polígono). Tipo de entidad: polígono multiparte.


    Registro publicado 2011-03-09
    Última modificación del registro 2021-06-10
    Estado del recurso en marcha

    Descripción del objeto

    Nombre del objeto: WHSE_BASEMAPPING.FWA_WATERSHEDS_POLY

    Nombre corto: FWWTRSHDSP
    Comentarios: Contiene todos los polígonos de cuencas hidrográficas generados a partir de las líneas fronterizas de las cuencas hidrográficas principales, los bordes de los bancos, los bordes delimitadores, los bordes de la línea costera y los bordes de los límites administrativos, siguiendo el modelo de datos 1: 50K. Puede contener geometrías de varias partes. Tipo espacial: geometría multiparte


    Registro publicado 2011-03-09
    Última modificación del registro 2021-06-26
    Estado del recurso en marcha

    Descripción del objeto

    Nombre del objeto: WHSE_BASEMAPPING.FWA_NAMED_WATERSHEDS_POLY

    Nombre corto: FWNMDWTRSH
    Comentarios: Contiene todas las cuencas hidrográficas de los ríos nombrados y los detalles del nombre asociado. Hay aproximadamente 12 000 cuencas hidrográficas nombradas en BC. Contiene geometrías multiparte. Gemetría espacial: polígono multiparte


    GIS 1 (335) Final

    2. Unión de distancia simple: mantiene la cardinalidad de la relación de unión y el criterio espacial es proximidad / distancia
    Resumir la unión espacial AKA a muchos se utiliza en una relación de uno a varios o de varios a varios cuando se viola la cardinalidad de la combinación simple

    3. Unión interna resumida: se usa cuando la característica de destino y la característica de origen exhiben una relación de cardinalidad de uno a muchos y la relación espacial es de contención

    2 Precisión posicional
    -Evaluación de qué tan cerca se corresponden las características de un conjunto de datos con sus posiciones reales en el mundo real

    3 Exactitud de atributo
    -Evalúa qué tan bien los valores de los atributos se corresponden con los valores en el mundo real

    4 Consistencia lógica
    -Describe qué tan bien la función imita las situaciones del mundo real
    -Por lo general, se limita a probar errores de topología.

    5 Lo completo
    -Evalúa si un conjunto de datos capturó cada instancia conocida de un conjunto de objetos del mundo real
    -Evalúa si todos los datos originales se llevaron al producto final


    Registro publicado 2011-03-09
    Última modificación del registro 2021-05-22
    Estado del recurso en marcha

    Descripción del objeto

    Nombre del objeto: WHSE_BASEMAPPING.FWA_ISLANDS_POLY

    Nombre corto: FWSLNDSPL
    Comentarios: Contiene todos los polígonos de islas. Las islas pueden superponerse ya que hay islas dentro de islas. Geometría espacial: polígono multiparte


    Ver el vídeo: Calculating area of polygon shapefile in QGIS