¿Obtener metros en lugar de grados de GRASS v.distance?

Question

Estoy tratando de medir la distancia entre puntos a una línea de costa. Tengo un shapefile de la "Alta Mar" que obtuve del sitio web de la NOAA aquí.

http://www.nws.noaa.gov/geodata/catalog/wsom/html/marinezones.htm

También tengo una lista de puntos que representan edificios.

He jugado con diferentes soluciones, pero creo que la v.distance de GRASS es la mejor solución. He podido introducir mis datos en GRASS y puedo trazar el mapa.

El problema es que cuando ejecuto la herramienta v.distance, obtengo grados decimales, pero realmente quiero metros, millas o kilómetros.

Así que la función v.distance parece ser la mejor opción.

Sin embargo, no he podido sacar metros de estos datos. En la guía de ayuda, para la v.distance, encontré esta frase.

In lat-long locations v.distance gives distances (dist and to_along) in meters not in degrees calculated as geodesic distances on a sphere.

No tengo muy claro qué significa esto. Parece que debo configurar la ubicación con un código EPSG que utiliza lat long. He utilizado un puñado de ellos sin éxito. Por ejemplo, he probado con los códigos EPSG 3857, 3786 y 4326, todos sin éxito.

Como nota al margen, he probado la distancia al centro de QGIS que forma parte del plugin MMQGIS. Parece que mapea el centroide más cercano de un polígono. Estoy tratando de obtener la distancia al píxel/borde más cercano de un polígono. O el píxel/borde más cercano de una línea.

Answer 1

A pesar de lo que dice el manual de GRASS para v.distance, las unidades son las mismas que las unidades del sistema de coordenadas, por lo que en un sistema de coordenadas lat/lon obtendrá efectivamente grados decimales. En el manual de GRASS 7 la declaración anterior (incorrecta) no aparece.

Por otro lado, si crea una UBICACIÓN basada en Spherical Mercator, EPSG 3857, como ha mencionado, las unidades son metros, y v.distance da las medidas de distancia en metros. Sin embargo, tenga en cuenta las inexactitudes inherentes a este sistema de coordenadas.

Answer 2

Para obtener las distancias en metros, proyectar los datos. La forma de hacerlo es una cuestión mecánica que se resuelve fácilmente consultando los manuales del software. La parte más difícil de la pregunta se refiere a ¿qué proyección utilizar?

Dado que un comentario a la pregunta indica que América del Norte es la región de interés, la forma de este continente -y la situación de sus costas a ambos lados- aconseja utilizar una proyección cónica o policónica. Una opción razonable es la Cónica Conforme de Lambert. Colocando los paralelos estándar en valores adecuados, esta proyección puede adaptarse para comprometerse decentemente entre la obtención de una buena cobertura de los océanos de las costas orientales y occidentales, minimizando al mismo tiempo el rango de distorsiones.

Dos excelentes herramientas para estudiar las distorsiones en las proyecciones -y por tanto para seleccionar entre las opciones de forma objetiva, inteligente y guiada- son el graticula de latitud-longitud y el Tissot Indicatrix. Esta figura ilustra ambas cosas.

Para lograr una buena precisión de ambas costas de los Estados Unidos, elegí los paralelos estándar de 45 y 25 grados Norte (esencialmente los límites superior e inferior de los estados conterminos). El paralelo central se fija en 35 grados Norte y el meridiano central (que determina la orientación del mapa) en 102 grados Oeste.

• La gratícula es la cuadrícula de líneas grises. Busque alturas iguales de las celdas de la cuadrícula en la dirección norte-sur y anchuras que varíen regularmente (que debe para encogerse hacia los polos). Estos patrones son signos cualitativos de una buena precisión métrica.

• El Indicadores Tissot son las figuras circulares trazadas en puntos estratégicos del océano. Cada una es en realidad dos círculos: uno, en gris, es de tamaño constante en todo el mapa. Otro, en azul, es el actual tamaño del círculo gris. Siempre que los círculos azul y gris coinciden, la precisión métrica es alta. (Se puede ver una parte de una indicatriz situada al oeste de Alaska en la esquina superior izquierda donde el círculo azul es visiblemente mayor que el gris. Esta proyección mostrará una distorsión considerable en las latitudes más allá de los paralelos estándar). Por cierto, los radios rojos siempre apuntan al norte. Su variación a lo largo del mapa muestra una considerable variación en convergencia que son importantes para los análisis que implican el rumbo o la dirección, pero son irrelevantes para las mediciones de distancia.

Las mediciones detalladas de estas indicatrices indican que las mediciones de distancia tendrán una precisión del 1,5% en toda la región que cubren en este mapa.

Como punto de partida para construir esta proyección, considere comenzar con EPSG:102004 o (para una cobertura mucho más amplia de norte a sur) EPSG:102009 .

Answer 3

Usted quiere que sus datos se proyecten en metros, lo que puede proyectar en GRASS o usar QGIS para hacerlo - la proyección EPSG 102008 de North_America_Albers_Equal_Area_Conic funcionaría para América del Norte, por ejemplo.

Entonces no estoy seguro si puedes tener datos en múltiples sistemas de coordenadas en 1 espacio de trabajo, o si debes crear un nuevo espacio de trabajo en EPSG: 102008 para almacenar tus datos proyectados y luego ejecutar tus cálculos de distancia.

Una vez que tenga los datos en metros en su espacio de trabajo de GRASS, sus valores de distancia estarán en metros.