¿Cómo puedo calcular la orientación de cada lado de un polígono en ArcGIS 10?

Question

Quiero examinar la orientación de cada línea en un polígono para que yo pueda calcular su exposición a la radiación solar. Cada polígono representa un edificio, y ha sido asociada a una altura. Por el momento, sólo quiero considerar la orientación, y después de considerar el sombreado de los problemas.

Un enfoque que yo pensaba que era para dividir el polígono en las líneas, y para calcular la orientación de cada línea, pero la dificultad es que tengo que identificar la cara exterior de la línea. Aunque la mayoría de los polígonos son simples cuatro del lado de figuras con líneas rectas, hay un pequeño número donde este no es el caso (un problema que sólo quiero a considerar, pero no tiene que resolver todavía).

Estoy familiarizado con python y planeado hacer todo esto desde una secuencia de comandos.

Answer 1

Si sólo desea la mayoría de la orientación, echa un vistazo @Mapperz 's respuesta anterior.

De lo contrario, como usted dice que usted podría dividir los polígonos en líneas usando el Polígono A la Línea de la herramienta. Esto añade una a la izquierda y a la derecha campo FID, donde el campo es -1 si no hay ningún polígono exterior - esto puede causar un poco de retirada de estiércol sobre si sus edificios son adyacentes o se superponen.

Desde allí se podía dividir sus líneas en cada vértice (tal vez de usar Dividido en COGO líneas) y, a continuación, calcular los ángulos en cada una de las líneas (posiblemente, por la Actualización de atributos COGO).

Asumiendo que usted tiene su ángulo de campo calculado a partir del Norte el aspecto será correcto cuando el left_FID es -1, y para obtener el aspecto cuando el right_FID es -1 acaba de añadir 180°. A continuación, basado en el original de la FID, usted puede agregar, conseguir que la mayoría de aspecto en función de la duración, etc.

El Polígono de la línea de herramienta de secuencias de comandos, (que yo sepa) COGO herramientas no son, por lo que tendría que venir con algo allí mismo.

Espero que esto ayude!

Answer 2

Calcular El Polígono Principal Ángulo (Cartografía)

Calcula el dominante de los ángulos de un polígono de entrada características y asigna los valores de un campo especificado en la clase de entidad

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//007000000028000000.htm

Answer 3

Búsqueda de orientación

Dentro de una secuencia de comandos, el polígono estará disponible como un conjunto de anillos-un anillo exterior y cero o más anillos interiores-con cada anillo representa por un cíclicamente ordenó tupla de vectores (v[0], v1, ..., v[m-1], v[m] = v[0]). Cada vector da las coordenadas de un vértice (no hay dos sucesivos vértices coincidentes). A partir de esto es sencillo, como otros han señalado, para obtener los vectores normales (es decir, los vectores perpendiculares al borde direcciones):

n[i] = t(v[i+1] - v[i]).

La "t" de la operación a la que gira un vector de 90 grados hacia la izquierda:

t( (x, y) ) = (y, -x).

Sólo las direcciones de estos vectores normales de la materia, por lo reescala a tener unidad de longitud: un vector (x, y) escala a (x/s, s/s), donde s = Sqrt(x^2 + y^2) (que es la longitud de su borde correspondiente). A partir de ahora, vamos a suponer que esto se ha hecho. Escribir las componentes de la resultante de la unidad de vectores normales como

n[i] = (u[i], v[i]), i = 0, 1, ..., m-1.

Discriminar fuera desde el interior

Como comentario, esto deja una direccional de la ambigüedad: debemos ser el uso de n[i] o -n[i]? Que uno de los puntos hacia el exterior? Esta pregunta es equivalente a encontrar el grado de Gauss mapa. Para calcular esto, usted necesita a la suma de los ángulos por los que la dirección de las normales de cambiar a medida que la marcha alrededor de un anillo. Debido a que los vectores normales tienen unidad de longitud, el coseno del ángulo entre dos sucesivas, de bordes

Cos(q_i) = n[i] . n[i+1] = u[i]*u[i+1] + v[i]*v[i+1], i = 0, 1, ..., m-1.

(Definir n[m] = n[0].)

El seno del ángulo entre dos sucesivas, de bordes

El pecado(q_i) = n[i] . t(n[i+1]) = u[i]*v[i+1] - v[i]*u[i+1].

(Tenga en cuenta que estos cálculos requieren sólo sumas, diferencias y productos hasta el momento.) La aplicación de los principales inversa de la función tangente (ATan2) para cualquier (coseno, seno) par da el ángulo q_i entre -180 y 180 grados. Sumando los ángulos para i = 0, 1, ..., n-1 produce (hasta el punto flotante de error) el total de la curvatura del anillo, que debe ser un múltiplo de 360 grados; para un cierre no-auto-intersección del anillo es wil +360 o -360. En el primer caso, el grado es 1 y en el segundo caso, el grado es -1. Las normales son todo exterior orientado a cuando el grado de que el anillo exterior es de +1 y los grados de los anillos interiores son de -1. Re-orientar el anillo por anillo, según sea necesario, de acuerdo a esta regla. Es decir, si el grado de cualquier anillo es el opuesto de lo que es necesario, anular todas las normales para que el anillo. Ahora puedes continuar con su insolación cálculos.

Answer 4

Puede esto ayudar?