6 votos

¿Zonal de Estadísticas de cambio de valor cuando se aplica a un subconjunto de datos?

Me han generado zonal estadísticas de uso de la característica de los datos de la zona en un valor de trama. Este shapefile contiene alrededor de 20000 polígonos y el análisis fue sin un tirón.

Mientras que las pruebas de calidad, me corrió la misma operación en 10 polígonos tomadas a partir de los mismos datos. Los valores en la zonal max cambiado esta vez. He comparado los valores mediante la superposición de los polígonos en el valor de trama y se encontró que los valores de los píxeles que pertenecen en el polígono, pero no sólo la máxima de píxeles.

Soy consciente de que un interno de la trama lleva a cabo la conversión para calcular las estadísticas. Quiero saber exactamente lo que sucede durante el proceso de la zonal de estadísticas y si el tamaño de la muestra afecta a la salida.

Tengo ArcGIS 10 con una licencia de Arcview y arcgis spatial analyst. Cualquier orientación es muy apreciada.

17voto

cjstehno Puntos 131

Zonal estadísticas en última instancia, es una trama en trama de la operación. Cuando los polígonos se proporcionan en formato vectorial, que se convertirá automáticamente en formato raster. Espero que este sea el mismo cellsize y registrada con el valor de trama, pero no estoy seguro de eso. (Llevaría algo difícil de ingeniería inversa para comprobar. Posiblemente los polígonos se convierten utilizando el análisis espacial de las propiedades son, en efecto, y luego se vuelven a muestrear para el zonal estadísticas de operación. Eso sería una tontería, pero no es difícil ver cómo el software puede ser diseñado de esa manera.)

(La conversión a formato raster explica por qué, por el camino, no se puede realizar directamente zonal estadísticas con la superposición de polígonos. Cada celda ráster puede representar a más de un polígono en un momento.)

Históricamente, Spatial Analyst ha tenido un tiempo duro con el zonal de estadísticas. Algunos relativamente recientes versiones (como recientemente como 9.3, recuerdo) plantea un montón de preguntas sobre el ESRI foros relativos a las respuestas incorrectas, falta de valores, y francamente errores cuando más de 2000 polígonos estaban involucrados. Algunos de estos (si no todos) sospecho que puede haber sido debido a los malentendidos acerca de cómo utilizar zonal estadísticas y lo que hace. En particular, la conversión de vector polígono de trama asigna un celular a un polígono si y sólo si la celda del centro está contenida en el polígono. (¿Cómo se toma la decisión, cuando el centro se encuentra en el polígono de la frontera es otro de los indocumentados misterio que requieren algún tedioso de la ingeniería inversa.) Cuando un tortuoso polígono logra hilo alrededor de toda la célula cercana de los centros, que se polígono simplemente desaparece de los resultados en conjunto!

Incluso hay más cosas bajo el capó y que, también, ha sido el cambio de versión del software. Cuando sobre la marcha de proyección está involucrado, los polígonos pueden o no ser que se vuelve a proyectar antes de la conversión de raster. Por lo tanto, puede que la trama en sí, para que la materia. Precisamente lo que sucede puede depender del entorno: si esto es en la Trama de la Calculadora, la línea de comandos, una secuencia de comandos de Python, una herramienta, o lo que sea.

Siendo conscientes de todo esto es que más de la mitad de la batalla, porque podemos protegernos a nosotros mismos desde el software. La mejor manera es crear un ráster de los polígonos que tiene el mismo cellsize, con domicilio a las células, como el valor de trama, y tiene la misma proyección. Deshabilitar cualquier remuestreo automático o en el de la mosca de la proyección. De esta manera minimizar el invisible "ayuda" del software es dar a usted y tener una posibilidad razonable de saber exactamente lo que las entradas para el cálculo. Si haces eso y todavía encontrar discrepancias, a continuación, (en mi humilde opinión) tiene amplia base formal de informe de error.

(Mediante los siguientes procedimientos religiosamente, no he tenido ningún problema con la zonal estadísticas en cualquier versión de Spatial Analyst (1.0 a través de 9.3), incluso con redes muy grandes y un gran número de polígonos. Por ejemplo, he obtenido los valores correctos para un conjunto de bloques Censales a lo largo de un gran estado de los EE.UU. en una resolución de 10 metros. (Que habría sido de más de 10^5 polígonos, recuerdo.) Sin embargo, yo a veces tomar medidas extraordinarias para impedir problemas: he logrado que el Censo del bloque de cálculo de forma automática la ruptura del estado en sus 50 o así condados, realizar el cálculo para cada condado, y, a continuación, volver a montar los resultados. Entre otras cosas, esto permite la paralelización de la computación. También me permitió manejar los pequeños polígonos, que no fueron capturados durante el proceso de rasterización, en un paso independiente mediante el cual el valor de cuadrícula fue consultado en el polígono de los centroides. Estos valores sirven como sustitutos de la zonal significa ser calculada.)

Por cierto, el tamaño de celda ("tamaño de muestra"?) afecta a algunas de salida en formas importantes. Si aumenta el valor de trama,

  1. El número de celdas de valor dentro de un polígono va a cambiar.

  2. Valores individuales son interpolados, causando un ligero cambio en su distribución estadística. Esto puede incluir cambios en los máximos y mínimos, que tienden a desplazarse hacia la media.

  3. Las estadísticas relacionadas con la uniformidad y la frecuencia (tal como un modo) puede cambiar por completo.

  4. Algunos geométricas resúmenes, especialmente perímetro, puede cambiar arbitrariamente.

  5. Muchos de los más llamativos son los cálculos estadísticos, tales como los errores estándar de las estimaciones o regresiones, puede ser completamente engañosa debido a la arbitrariedad de los cambios en el número de células (ver #1).

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X