¿Cálculo de la longitud y anchura de un polígono mediante un SIG de código abierto?

Question

¿Cómo puedo calcular la longitud y con de polígonos utilizando herramientas de código abierto? Por "longitud del polígono" me refiero a la longitud de la línea más larga dentro del polígono (actualización: para ser correcto, sólo necesito la línea más larga dentro del casco convexo del polígono) y la anchura como la medida más larga perpendicular a la medida de longitud.

Intento crear un rectángulo mínimo no alineado con el eje para un polígono. La longitud de este rectángulo mínimo sería el diámetro del polígono.

Answer 1

FWIW, en Postgis hay una función desde 1.5, ST_Longestline que devuelve la longestline entre dos geometrías. Si usted alimenta esa función con la misma geometría dos veces en lugar de dos diferentes, obtendrá la línea más larga en esa geometría. http://postgis.net/docs/ST_LongestLine.html

@Whuber. usted escribe arriba que la línea más larga no tiene que comenzar y terminar en los vértices de los polígonos. Utilizas la letra I como ejemplo. No lo entiendo. Si I está representado por una línea con dos vértices, entonces esa línea tiene que ser la línea más larga en sí, ¿verdad?

Si la letra I está representada por un polígono delgado, la línea más larga será una línea diagonal desde, por ejemplo, el vértice superior izquierdo hasta el vértice inferior derecho.

Pero no puedo entender cómo puede haber una línea más larga que no termine y empiece en vértices del polígono.

¿He entendido algo mal?

/Nicklas

Ok, Whuber, te entendí mal, leí tu respuesta una vez más. Yo estaba hablando de lo que usted llama el diámetro y que estaba hablando de otra cosa. Lo siento

/Nicklas

Answer 2

El STC proporciona un círculo mínimo ( http://tsusiatsoftware.net/jts/javadoc/com/vividsolutions/jts/algorithm/MinimumBoundingCircle.html ) y un Diámetro Mínimo ( http://tsusiatsoftware.net/jts/javadoc/index.html ) - una vez que los tengas, el rectángulo debería ser fácil de calcular.

Answer 3

Utilice Postgis y use ST_envelope(geometry) para obtener las coordenadas mínimas de la caja delimitadora (definida con cuatro puntos). A partir de este resultado, sólo tienes que calcular ahora longitud y anchura

Una solución que utiliza consultas Postgis puras. Hace el trabajo, pero yo no hago la prueba de velocidad, yo uso algunos wkb para la depuración rápida, por lo que este script se puede mejorar.

SELECT line1,
       line2,
       CASE WHEN line1>line2 THEN line1
            WHEN line1<=line2 THEN line2           
       END AS length,
       CASE WHEN line2>line1 THEN line1
            WHEN line2<=line1 THEN line2           
       END AS width FROM (
SELECT ST_length(ST_AsText(ST_MakeLine(ST_AsText(ST_PointN(ST_Boundary(ST_envelope(the_geom)),1)),
ST_AsText(ST_PointN(ST_Boundary(ST_envelope(the_geom)),2))))) AS line1,ST_length(ST_AsText(ST_MakeLine(ST_AsText(ST_PointN(ST_Boundary(ST_envelope(the_geom)),2)),ST_AsText(ST_PointN(ST_Boundary(ST_envelope(the_geom)),3))))) AS line2
from tm_world_borders) AS size;

Answer 4

Si no puedes encontrar una herramienta específica para hacerlo podrías usar pgSQL y PostGIS o conseguir un desarrollador que implemente ese algoritmo por ti.

Hay muchas maneras de hacerlo, pero aquí hay una dividida en 2 etapas:

1) Para cada vértice encontrar la mayor distancia a cada otro vértice dentro del polígono que se puede hacer con un simple pequeño bucle anidado en un SISD como .NET:

for(var i=0; i<count; i++)
 for(var j=i; j<count; j++)
  find the 2 vertices with the greatest distance
construct a segment from those two points
var length = Distance(segment.Begin,segment.End);

2) para la anchura hay algunas maneras, la más fácil de implementar creo que sería esta:

// make the longest segment horizontal
transform the polygon with a rotation of the negative angle of the longest segment
for each vertex in the polygon find the highest and lowest y-coord values
// the "width" would be the max y-coord minus the min y-coord
var width = maxy-miny

Answer 5

Esto resuelve el problema en PostGIS - probado y funciona. Gracias a Mapbox por esta gran y elegante solución. ¡Tenga en cuenta su licencia de uso!

https://raw.githubusercontent.com/mapbox/postgis-vt-util/master/postgis-vt-util.sql