Mise en mémoire tampon en mètres / km à l'aide des couches WGS84 à l'aide de QGIS

12

J'utilise QGIS 1.7.4-Wroclaw et travaille avec une couche SHP dans WGS84. Je veux mettre en mémoire tampon à un certain nombre de kilomètres de cette couche en utilisant l'outil tampon ftools.

Je comprends que l'outil tampon utilise toujours les unités de calque, qui pour WGS84 sont des degrés décimaux. Nous savons tous que les degrés ne se convertissent pas systématiquement en mètres, alors comment dois-je procéder pour créer mes tampons?

Dois-je convertir les fichiers de formes en un autre CRS qui utilise nativement les compteurs?

Si oui, comment en choisir un?

Evan
la source
Quelles sont les étendues de votre fichier de formes? Nous aurions besoin de ces informations pour recommander un CRS.
nhopton
Oui, vous devrez d'abord reprojeter le Shapefile (en utilisant l'outil "Reproject layer" de fTools, par exemple). La meilleure projection dépend de l'endroit où se trouve votre zone (longitude / latitude) et de sa taille.
Jake
Merci, @Jake et @nhopton! Je tamponne depuis plusieurs points dans le nord de Taiwan, mais les tampons seront de l'ordre de 1000-2000 km (couvrant environ la moitié de la Chine ainsi que Taiwan).
Evan
1
Cette carte est définie avec un système de coordonnées géographiques (latitude / longitude), mais l'étendue de la carte suggère qu'il s'agit en fait d'un système de coordonnées projeté (par exemple, Mercator). Si c'est le cas, les résultats des mesures de ligne ou de surface seront incorrects. Pour résoudre ce problème, définissez explicitement un système de coordonnées cartographiques approprié à l'aide du menu Paramètres: Propriétés du projet.
Ger
C'est juste un fichier recadré Natural Earth, donc c'est juste dans WGS84.
Evan

Réponses:

12

Pour cette application, j'utiliserais une projection azimutale équidistante centrée au milieu de vos points sources. Cette projection a la particularité que toutes les distances radiales autour du centre de la projection sont précises.

Cette projection particulière ne fait pas partie des projections standard de QGIS. Vous pouvez définir le vôtre en utilisant Paramètres / CRS personnalisé avec la chaîne de commande +proj=aeqd +lat_0=24.5 +lon_0=121.5, mais malheureusement, les projections personnalisées ne peuvent pas être utilisées par l'outil de reprojection fTools. Au lieu de cela, vous pouvez transformer votre jeu de données sur la ligne de commande à l'aide de la commande

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

Vous pouvez ensuite effectuer la mise en mémoire tampon dans QGIS à l'aide du points_reprojected.shpfichier de formes.

Un tampon de 2000 km autour d'un point dans le nord de Taïwan ressemble à un cercle dans une projection orthographique centrée sur le point:

... et écrasé dans WGS84:

Jake
la source
Merci - reprojeter à Azimuthal Equidistant puis tamponner a fonctionné à merveille! En fait, j'utilisais déjà cette projection pour la projection à la volée, je ne savais tout simplement pas qu'elle aurait les unités dont j'avais besoin. Mais je n'ai pas eu à utiliser la ligne de commande pour reprojeter: Layer> Save as ... a fait l'affaire (après avoir défini le CRS personnalisé).
Evan
@Evan: Heureux de l'entendre! Si vous considérez que votre question a reçu une réponse, vous pouvez la marquer comme telle en cliquant sur la coche à côté de la réponse
Jake
1

Hum, intéressant. Comme un premier essai, je pense que j'aurais un regard sur la génération des (énormes) tampons en KML. Vous pouvez le faire ici:

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

Vous devrez entrer le lat / lon pour un point et une distance tampon. Ce que vous récupérez est un fichier KML contenant le cercle. Jetez un œil à cela dans Google Earth et si cela vous convient, chargez-le dans QGIS.

entrez la description de l'image ici

nhopton
la source
1
Malheureusement, QGIS Lisboa n'aime pas le fichier kml généré.
AndreJ
QGIS n'aime pas le nom de fichier. Essayez de le changer en (par exemple) "cercle1.kml". Mais j'aurais dû vérifier. N.
nhopton
D'accord, chargez le fichier dans Google Earth, puis enregistrez-le au format KML. Modifiez ensuite le nom de fichier à partir du nom par défaut. ogr2ogr convertit également très bien le fichier en SHP. N.
nhopton
La solution de Jake a fini par fonctionner parfaitement pour moi, mais cela pourrait aussi être un bon moyen de le faire. Je m'en souviendrai. Merci!
Evan