r.projgrass - En ligne dans le Cloud

PROGRAMME:

Nom

r.proj - Re-projette une carte raster de l'emplacement donné à l'emplacement actuel.

MOTS-CLÉS

raster, projection, transformation

SYNOPSIS

r.proj
r.proj --Aidez-moi
r.proj [-GNL] Localisation=prénom [jeu de cartes=prénom] [contribution=prénom] [dbase=chemin] [sortie=prénom]
[méthode=un magnifique] [Mémoire=entier] [RAPIDE=flotter] [--écraser] [--vous aider]
[--verbeux] [--calme] [--ui]

Drapeaux:
-l
Répertorier les cartes raster dans le jeu de cartes en entrée et quitter

-n
Ne pas effectuer d'optimisation de recadrage de région

-p
Imprimer les limites de la carte d'entrée dans la projection actuelle et quitter

-g
Imprimer les limites de la carte d'entrée dans la projection actuelle et la sortie (style shell)

--écraser
Autoriser les fichiers de sortie à écraser les fichiers existants

--Aidez-moi
Imprimer le récapitulatif d'utilisation

--verbeux
Sortie du module verbeux

--silencieux
Sortie module silencieuse

--interface utilisateur
Forcer le lancement de la boîte de dialogue GUI

Paramètres:
Localisation=prénom [obligatoire]
Emplacement contenant la carte raster en entrée
Nom de l'emplacement (pas le chemin de l'emplacement)

jeu de cartes=prénom
Jeu de cartes contenant une carte raster en entrée
Par défaut : nom du jeu de cartes actuel

contribution=prénom
Nom de la carte raster en entrée à reprojeter

dbase=chemin
Chemin d'accès à la base de données GRASS de l'emplacement d'entrée
Par défaut : chemin vers la base de données SIG GRASS actuelle

sortie=prénom
Nom de la carte raster en sortie (par défaut : identique à « entrée »)

méthode=un magnifique
Méthode d'interpolation à utiliser
Options: la plus proche, bilinéaire, bicubique, les lanczos, bilinéaire_f, bicubique_f, lanczos_f
Valeur par défaut: le plus proche
le plus proche: voisin le plus proche
bilinéaire: interpolation bilinéaire
bicubique: interpolation bicubique
lanczos: filtre lanczos
bilinéaire_f: interpolation bilinéaire avec repli
bicubique_f: interpolation bicubique avec repli
lanczos_f: filtre lanczos avec repli

Mémoire=entier
Mémoire maximale à utiliser (en Mo)
Taille du cache pour les lignes raster
Valeur par défaut: 300

RAPIDE=flotter
Résolution de la carte raster en sortie

DESCRIPTION

r.proj projette une carte raster dans un jeu de cartes spécifié d'un emplacement spécifié à partir du
projection de l'emplacement d'entrée sur une carte raster à l'emplacement actuel. La projection
les informations sont extraites des fichiers PROJ_INFO actuels, tels qu'ils sont définis et visualisés avec g.proj.

Introduction
Carte projections
Les projections cartographiques sont une méthode de représentation d'informations à partir d'une surface courbe (généralement une
sphéroïde) en deux dimensions, généralement pour permettre l'indexation par coordonnées cartésiennes.
Il existe une grande variété de projections, les communes étant divisées en un certain nombre de
classes, y compris cylindrique et pseudo-cylindrique, conique et pseudo-conique, et
méthodes azimutales, dont chacune peut être conforme, à aire égale ou aucune.

La projection particulière choisie dépend de l'objectif du projet et de la taille,
la forme et l'emplacement de la zone d'intérêt. Par exemple, des projections cylindriques normales
sont bonnes pour les cartes qui sont plus étendues est-ouest que nord-sud et en équatorial
régions, tandis que les projections coniques sont meilleures aux latitudes moyennes ; cylindrique transversale
les projections sont utilisées pour les cartes qui sont plus étendues nord-sud qu'est-ouest;
les projections azimutales sont utilisées pour les régions polaires. Les versions obliques de l'un d'entre eux peuvent
également être utilisé. Les projections conformes préservent les relations angulaires, et mieux préservent
longueur d'arc, tandis que les projections à aire égale sont plus appropriées pour les études statistiques et
travail dans lequel la quantité de matière est importante.

Les projections sont définies par des relations mathématiques précises, de sorte que la méthode de projection
coordonnées d'un référentiel géographique (latitude-longitude) dans une projection
Le référentiel cartésien (par exemple les mètres) est régi par ces équations. Projections inverses
peut également être atteint. Le progiciel Unix du domaine public PROJ.4 [1] a été
conçu pour effectuer ces transformations, et le manuel de l'utilisateur contient un
description de plus de 100 projections utiles. Cela comprend également une bibliothèque de programmeurs de
les méthodes de projection pour soutenir le développement d'autres logiciels.

Ainsi, la conversion d'une carte vectorielle - dans laquelle les objets sont localisés avec des paramètres spatiaux arbitraires
précision - d'une projection dans une autre est généralement accomplie par une simple
processus : d'abord, l'emplacement de tous les points de la carte est converti à partir de la source
par une projection inverse en latitude-longitude, puis par une projection vers l'avant
projection dans la cible. (Bien sûr, la procédure sera en une étape si soit le
la source ou la cible est en coordonnées géographiques.)

Cependant, la conversion d'une carte raster ou d'une image entre différentes projections implique
Considérations supplémentaires. Un raster peut être considéré comme représentant un échantillonnage d'un
processus à un ensemble régulier et ordonné d'emplacements. L'ensemble des emplacements qui se trouvent au
les intersections d'une grille cartésienne dans une projection ne coïncideront généralement pas avec
les points d'échantillonnage dans une autre projection. Ainsi, la conversion de cartes raster implique une
étape d'interpolation dans laquelle les valeurs des points à des emplacements intermédiaires par rapport à la
grille source sont estimées.

En saillie vecteur cartes dans les le GRASS SIG
La capture, l'importation et le transfert de données SIG nécessitent souvent une étape de projection, car la source
ou le client sera fréquemment dans une projection différente de la projection de travail.

Dans certains cas, il est pratique d'effectuer la conversion en dehors du package, avant l'importation
ou après exportation, à l'aide de logiciels tels que PROJ.4's cs2cs [1]. C'est une méthode facile pour
convertir un fichier ASCII contenant une liste de points de coordonnées, car il n'y a pas
topologie à conserver et cs2cs peut être utilisé pour traiter des listes simples en utilisant une ligne
commander. le m.proj module fournit une interface pratique à cs2cs.

Les cartes vectorielles sont généralement plus complexes, car des parties des données stockées dans les fichiers seront
décrire la topologie, et pas seulement les coordonnées. Dans GRASS GIS, le v.proj module est fourni à
reprojetez les cartes vectorielles, transférez la topologie et les attributs ainsi que les coordonnées des nœuds.
Ce programme utilise la définition de projection et les paramètres qui sont stockés dans le
Fichiers PROJ_INFO et PROJ_UNITS dans le répertoire de mapset PERMANENT pour chaque emplacement GRASS.

Design of r.proj
Comme indiqué brièvement ci-dessus, l'étape fondamentale de la re-projection d'un raster est le rééchantillonnage
la grille source aux emplacements correspondant aux intersections d'une grille dans la cible
projection. La procédure de base pour y parvenir est donc la suivante :

r.proj convertit une carte en une nouvelle projection géographique. Il lit une carte d'un autre
emplacement, le projette et l'écrit à l'emplacement actuel. Les données projetées sont
rééchantillonné avec l'une des quatre méthodes différentes : plus proche voisin, bilinéaire, bicubique
iterpolation ou lanczos.

Votre méthode=la plus proche méthode, qui effectue une affectation de voisin le plus proche, est la plus rapide des
les trois méthodes de rééchantillonnage. Il est principalement utilisé pour les données catégorielles telles que l'utilisation des terres
classification, car cela ne changera pas les valeurs des cellules de données. Les méthode=bilinéaire
méthode détermine la nouvelle valeur de la cellule sur la base d'une moyenne de distance pondérée des 4
cellules environnantes dans la carte d'entrée. Les méthode=bicubique méthode détermine la nouvelle valeur de
la cellule sur la base d'une moyenne de distance pondérée des 16 cellules environnantes dans l'entrée
carte. le méthode=lanzcos La méthode détermine la nouvelle valeur de la cellule en fonction d'une valeur pondérée
distance moyenne des 25 cellules environnantes dans la carte d'entrée. Par rapport au bicubique,
lanczos met un poids plus élevé sur les cellules proches du centre et un poids plus faible sur les cellules éloignées
du centre, résultant en un contraste légèrement meilleur.

Les méthodes d'interpolation bilinéaire, bicubique et lanczos sont les plus appropriées pour
données continues et provoquer un certain lissage. La quantité de lissage diminue de bilinéaire
à bicubique à lanczos. Ces options ne doivent pas être utilisées avec des données catégorielles, car le
les valeurs des cellules seront modifiées.

Dans les méthodes bilinéaire, bicubique et lanczos, si l'une des cellules environnantes est utilisée pour
interpoler la nouvelle valeur de cellule sont nulles, la cellule résultante sera nulle, même si le
la cellule la plus proche n'est pas nulle. Cela entraînera un certain amincissement le long des frontières nulles, telles que le
côtes des terres émergées d'un DEM. L'interpolation bilinear_f, bicubic_f et lanczos_f
Les méthodes peuvent être utilisées si l'amincissement le long des bords nuls n'est pas souhaité. Ces méthodes « tombent
back" à des méthodes d'interpolation plus simples le long de frontières nulles. C'est-à-dire de lanczos à
bicubique à bilinéaire au plus proche.

Si l'affectation du voisin le plus proche est utilisée, la carte en sortie a le même format raster que le
carte d'entrée. Si l'une des interpolations est utilisée, la carte de sortie est écrite comme flottante
point.

Notez que, suivant les conventions normales de GRASS, la couverture et la résolution du
la grille résultante est définie par les paramètres de région actuels, qui peuvent être ajustés à l'aide
g.région. Le raster cible sera relativement non biaisé dans tous les cas si sa grille a un
résolution similaire à la source, de sorte que l'étape de rééchantillonnage/interpolation n'est qu'une
fonctionnement local. Si la résolution est modifiée de manière significative, le comportement du
la généralisation ou le raffinement dépendra du modèle du processus représenté.
Ce sera très différent pour les données catégoriques par rapport aux données numériques. Notez que trois
des méthodes pour l'étape d'interpolation locale sont fournies.

r.proj prend en charge les transformations de données générales, en utilisant le PROJ.4 système de coordonnées
bibliothèque de traduction.

NOTES

If sortie n'est pas spécifié, il est défini pour être le même que le nom de la carte d'entrée.
If jeu de cartes n'est pas spécifié, son nom est supposé être le même que celui de l'ensemble de cartes actuel
nom.
If dbase n'est pas spécifié, il est supposé être la base de données actuelle. L'utilisateur n'a qu'à
spécifier dbase si l'emplacement source est stocké dans une autre base de données GRASS distincte.

Pour éviter une consommation de temps excessive lors de la reprojection d'une carte, la région et la résolution de
l'emplacement cible doit être défini de manière appropriée au préalable.

Un moyen simple de le faire est de vérifier les limites projetées de la carte d'entrée dans le courant
projection de l'emplacement à l'aide du -p drapeau. Les -g flag signale la même chose, mais sous une forme
qui peut être directement coupé et collé dans un g.région commander. Après avoir défini la région dans
de cette façon, vous pouvez vérifier la résolution de la cellule avec "g.région -p" puis enclenchez-le sur un
grille avec g.région's -a drapeau. Par exemple g.région -a res=5 -p. Notez qu'il ne s'agit que d'une approximation
guider.

Une manière plus complexe, mais plus précise, de procéder consiste à générer une carte vectorielle "boîte" de
la région dans l'emplacement source en utilisant v.dans.région -d. Cette carte "boîte" est alors
reprojeté dans l'emplacement cible avec v.proj. Suivant la région dans l'emplacement cible
est fixé à l'étendue de la nouvelle carte vectorielle avec g.région avec le raster souhaité
résolution (g.région -m peut être utilisé dans des emplacements de latitude/longitude pour mesurer la
longueur géodésique d'un pixel). r.proj est ensuite exécuté pour la carte raster que l'utilisateur souhaite
reprojeter. Dans ce cas, un peu de préparation va un long chemin.

Lors de la reprojection de cartes du monde entier, l'utilisateur doit désactiver le découpage de carte avec le -n drapeau.
Le rognage n'est pas utile ici car le module a de toute façon toute la carte en mémoire. outre
que les "bords" du monde sont difficiles (ou impossibles) à trouver dans des projections autres que
latitude-longitude de sorte que les résultats peuvent être étranges avec le rognage.

EXEMPLES

En ligne méthode
Avec GRASS en cours d'exécution dans l'emplacement de destination, utilisez le -g indicateur pour afficher la carte d'entrée
limites une fois projetées dans la projection de travail actuelle, puis utilisez-les pour définir la région
limites avant d'effectuer la reprojection :
# calcule où sera la carte de sortie
r.proj input=emplacement d'élévation=ll_wgs84 mapset=user1 -p
Cols sources : 8162
Lignes sources : 12277
Local nord : -4265502.30382993
Local sud : -4473453.15255565
Local ouest : 14271663.19157564
Local est : 14409956.2693866
# même calcul, mais sous une forme qui peut être coupée et collée dans un appel g.region
r.proj input=emplacement d'élévation=ll_wgs84 mapset=user1 -g
n=-4265502.30382993 s=-4473453.15255565 w=14271663.19157564 e=14409956.2693866 lignes=12277 cols=8162
g.région n=-4265502.30382993 s=-4473453.15255565 \
w=14271663.19157564 e=14409956.2693866 lignes=12277 cols=8162 -p
projection : 99 (Mercator)
zone : 0
donnée : wgs84
ellipsoïde : wgs84
nord : -4265502.30382993
sud : -4473453.15255565
ouest : 14271663.19157564
est : 14409956.2693866
nbre : 16.93824621
aigrettes : 16.94352828
lignes : 12277
cols : 8162
cellules : 100204874
# résolution arrondie à quelque chose de plus propre
g.région res=17 -a -p
projection : 99 (Mercator)
zone : 0
donnée : wgs84
ellipsoïde : wgs84
nord : -4265487
sud : -4473465
ouest : 14271653
est : 14409965
nbre : 17
aigrettes : 17
lignes : 12234
cols : 8136
cellules : 99535824
# enfin, effectuer la reprojection
r.proj input=emplacement d'élévation=ll_wgs84 mapset=user1 memory=800

v.dans.région méthode
# Dans l'emplacement source, utilisez v.in.region pour générer un cadre de délimitation autour du
# région d'intérêt:
v.in.region -d sortie = type de limites = zone
# Maintenant, passez à l'emplacement cible et importez la boîte englobante vectorielle
# (vous pouvez exécuter v.proj -l pour obtenir une liste de cartes vectorielles à l'emplacement source) :
v.proj input=limites emplacement=source_location_name sortie=bounds_reprojected
# Définir la région dans l'emplacement cible avec celle du vecteur nouvellement importé
# carte des limites et alignez la résolution sur la résolution de cellule souhaitée du
# carte raster finale reprojetée :
g.region vector=bounds_reprojected res=5 -a
# Maintenant, reprojetez le raster dans l'emplacement cible
r.proj entrée=élévation.dem sortie=élévation.dem.reproj \
location=source_location_name mapset=PERMANENT res=5 method=bicubique

Références

1 Evenden, GI (1990) Procédures de projection cartographique pour l'environnement UNIX -
un manuel d'utilisation. USGS Open-File Report 90-284 (OF90-284.pdf) Voir aussi ici :
Rapport intérimaire et 2e rapport intérimaire sur la version 4, Evenden 1994).

2 Richards, John A. (1993), Analyse d'images numériques par télédétection, Springer-Verlag,
Berlin, 2e édition.

PROJ 4: Bibliothèque de support de projection/datum.

Plus loin lire

· Grilles et données ASPRS

· Liste de transformation des projections (PROJ.4)

· MapRef - La collection de projections cartographiques et de systèmes de référence pour l'Europe

· Système d'information et de service pour les systèmes européens de coordonnées de référence - CRS

· Projections cartographiques de Carlos A. Furuti

Utilisez r.projgrass en ligne en utilisant les services onworks.net