Almacenamiento intermedio en metros / km usando capas WGS84 usando QGIS

12

Estoy usando QGIS 1.7.4-Wroclaw y estoy trabajando con una capa SHP en WGS84. Quiero amortiguar a un cierto número de kilómetros desde esta capa usando la herramienta de protección ftools.

Entiendo que la herramienta de búfer siempre usa las unidades de capa, que para WGS84 son grados decimales. Todos sabemos que los grados no se convierten constantemente en metros, entonces, ¿cómo debo hacer para hacer mis amortiguadores?

¿Es necesario para mí convertir los archivos de forma a un CRS diferente que usa medidores de forma nativa?

Si es así, ¿cómo elijo uno?

qgis buffer units Evan
fuente
¿Cuáles son las extensiones de tu shapefile? Necesitaríamos esta información para recomendar un CRS.
nhopton
Sí, primero deberá volver a proyectar el Shapefile (por ejemplo, utilizando la herramienta "Reproyectar capa" de fTools). La mejor proyección depende de dónde se encuentre su área (longitud / latitud) y qué tan grande sea.
Jake
Gracias, @Jake y @nhopton! Estoy almacenando en búfer desde varios puntos en el norte de Taiwán, pero los búferes estarán en el rango de 1000-2000 km (cubriendo aproximadamente la mitad de China y Taiwán).
Evan
1
Este mapa se define con un sistema de coordenadas geográficas (latitud / longitud) pero la extensión del mapa sugiere que en realidad es un sistema de coordenadas proyectadas (por ejemplo, Mercator). Si es así, los resultados de las mediciones de línea o área serán incorrectos. Para solucionar esto, establezca explícitamente un sistema de coordenadas de mapa apropiado utilizando el menú Configuración: Propiedades del proyecto
Ger
Es solo un archivo recortado de Natural Earth, en realidad, así que está solo en WGS84.
Evan

Respuestas:

12

Para esta aplicación, usaría una proyección equidistante azimutal centrada en el centro de los puntos de origen. Esta proyección tiene la característica agradable de que todas las distancias radiales alrededor del centro de la proyección son precisas.

Esa proyección particular no es parte de las proyecciones estándar de QGIS. Puede definir su propio uso de Configuración / CRS personalizado con la cadena de comando +proj=aeqd +lat_0=24.5 +lon_0=121.5, pero desafortunadamente, la herramienta de Reproyección de fTools no puede usar proyecciones personalizadas. En cambio, puede transformar su conjunto de datos en la línea de comando usando el comando

ogr2ogr points_reprojected.shp points.shp -t_srs "+proj=aeqd +lat_0=24.5 +lon_0=121.5"

Luego puede hacer el almacenamiento en búfer en QGIS usando el points_reprojected.shparchivo de forma.

Un búfer de 2000 km alrededor de un punto en el norte de Taiwán parece un círculo en una proyección ortográfica centrada en el punto:

... y aplastado en WGS84:

Jake
fuente
Gracias, ¡la reproyección al Equidistante azimutal y luego el almacenamiento en búfer funcionaron maravillosamente! De hecho, ya estaba usando esa proyección para la proyección sobre la marcha, simplemente no me di cuenta de que tendría las unidades que necesitaba. Pero no tuve que usar la línea de comando para reproyectar: ​​Capa> Guardar como ... hizo el truco (después de definir el CRS personalizado).
Evan
@Evan: ¡Me alegra escucharlo! Si considera respondida su pregunta, puede marcarla como tal haciendo clic en la marca de verificación junto a la respuesta
Jake
1

Hum, interesante. Como primer intento, creo que podría echar un vistazo a generar los buffers (enormes) como KML. Puedes hacer esto aquí:

http://www.thesamestory.com/kmlcircle/

Deberá ingresar el lat / lon para un punto y una distancia de búfer. Lo que obtienes es un archivo KML que contiene el círculo. Eche un vistazo a esto en Google Earth y, si está bien, cárguelo en QGIS.

ingrese la descripción de la imagen aquí

nhopton
fuente
1
Desafortunadamente, a QGIS Lisboa no le gusta el archivo kml generado.
AndreJ
QGIS no le gusta el nombre del archivo. Intente cambiarlo a (por ejemplo) "circle1.kml". Pero debería haberlo comprobado. N.
nhopton
Bien, cargue el archivo en Google Earth y luego guárdelo como KML. Luego cambie el nombre del archivo por el nombre predeterminado. ogr2ogr convierte el archivo a SHP bastante bien también. N.
nhopton
La solución de Jake terminó funcionando perfectamente para mí, pero esta podría ser una buena manera de hacerlo también. Lo tendré en cuenta. ¡Gracias!
Evan