Il s'agit de la commande GeoidEval qui peut être exécutée dans le fournisseur d'hébergement gratuit OnWorks en utilisant l'un de nos multiples postes de travail en ligne gratuits tels que Ubuntu Online, Fedora Online, l'émulateur en ligne Windows ou l'émulateur en ligne MAC OS
PROGRAMME:
Nom
GeoidEval -- recherche les hauteurs du géoïde
SYNOPSIS
ÉvalGéoïde [ -n prénom ] [ -d dir ] [ -l ] [ -a | -c sud ouest nord est ] [ -w ] [ -z Zone
] [ --msltohae ] [ --haetomsl ] [ -v ] [ --comment-delimiter commenterdelim ] [ --version |
-h | --Aidez-moi ] [ --fichier-d'entrée dans le fichier | --chaîne-d'entrée chaîne ] [ --séparateur de ligne saut de ligne
] [ --fichier de sortie fichier de sortie ]
DESCRIPTION
ÉvalGéoïde lit les positions sur l'entrée standard et imprime les hauteurs correspondantes de
le géoïde au-dessus de l'ellipsoïde WGS84 en sortie standard.
Les positions sont données en latitude et longitude, UTM/UPS ou MGRS, dans l'un des formats
accepté par GéoConvertir(1). (Les coordonnées MGRS signifient centre du MGRS correspondant
carré.) Si le -z option est spécifiée, la zone spécifiée est ajoutée au début de chaque ligne
d'entrée (qui doit être en coordonnées UTM/UPS). Cela permet un fichier avec des abscisses UTM et
les ordonnées dans une seule zone à utiliser comme entrée standard.
Des résultats plus précis pour la hauteur du géoïde sont fournis par Gravité(1). Cet utilitaire peut
calculer également la direction de la gravité avec précision.
La hauteur du géoïde au-dessus de l'ellipsoïde, N, est parfois appelée l'ondulation du géoïde.
Il peut être utilisé pour convertir une hauteur au-dessus de l'ellipsoïde, h, à la hauteur correspondante
au-dessus du géoïde (la hauteur orthométrique, à peu près la hauteur au-dessus du niveau moyen de la mer), H,
en utilisant les relations
h = N + H, H =-N + h.
OPTIONS
-n utiliser le géoïde prénom au lieu du "egm96-5" par défaut. Voir "GEOIDES".
-d lire les données du géoïde à partir de dir au lieu de la valeur par défaut. Voir "GEOIDES".
-l utiliser l'interpolation bilinéaire au lieu de cubique. Voir "INTERPOLATION".
-a mettre en cache l'ensemble des données en mémoire. Voir "CACHE".
-c cache les données délimitées par sud ouest nord est en mémoire. Les deux premiers arguments
spécifiez le coin SW du cache et les deux derniers arguments spécifient le coin NE.
Notre -w flag spécifie que la longitude précède la latitude pour ces coins, à condition
qu'il apparaisse avant -c. Voir "CACHE".
-w lors de la lecture des coordonnées géographiques, la longitude précède la latitude (cela peut être
remplacé par un indicateur d'hémisphère, N, S, E, W).
-z préfixer chaque ligne d'entrée par Zone, par exemple, "38n". Cela doit être utilisé lorsque l'entrée
se compose d'abscisse et d'ordonnée UTM/UPS.
--msltohae
l'entrée standard doit inclure un jeton final sur chaque ligne qui est traité comme une hauteur
(en mètres) au-dessus du géoïde et la sortie fait écho à la ligne d'entrée avec la hauteur
converti en hauteur au-dessus de l'ellipsoïde (HAE). Si -z Zone est spécifié alors le troisième
le jeton est traité comme la hauteur ; cela permet de convertir les données LIDAR où
chaque ligne se compose de : est-nord hauteur intensité.
--haetomsl
c'est similaire à --msltohae sauf que le jeton de hauteur est traité comme une hauteur (dans
mètres) au-dessus de l'ellipsoïde et la sortie fait écho à la ligne d'entrée avec la hauteur
convertie en hauteur au-dessus du géoïde (MSL).
-v imprimer des informations sur le géoïde sur l'erreur standard avant de traiter l'entrée.
--comment-delimiter
définir le délimiteur de commentaire sur commenterdelim (par exemple, "#" ou "//"). S'il est défini, l'entrée
les lignes seront scannées pour ce délimiteur et, s'il est trouvé, le délimiteur et le reste de
la ligne sera supprimée avant le traitement et ensuite ajoutée à la sortie
ligne (séparée par un espace).
--version
version imprimée et quitter.
-h imprimer l'utilisation, le chemin et le nom du géoïde par défaut, et quitter.
--Aidez-moi
imprimez la documentation complète et quittez.
--fichier-d'entrée
lire l'entrée du fichier dans le fichier au lieu de l'entrée standard ; un nom de fichier de "-"
signifie entrée standard.
--chaîne-d'entrée
lire l'entrée de la chaîne chaîne au lieu de l'entrée standard. Toutes les occurrences
du caractère de séparation de ligne (la valeur par défaut est un point-virgule) dans chaîne sont convertis en
sauts de ligne avant le début de la lecture.
--séparateur de ligne
définir le caractère de séparation de ligne sur saut de ligne. Par défaut, il s'agit d'un point-virgule.
--fichier de sortie
écrire la sortie dans le fichier fichier de sortie au lieu de la sortie standard ; un nom de fichier de "-"
signifie sortie standard.
GÉOIDES
ÉvalGéoïde calcule les hauteurs du géoïde en interpolant sur les données dans une table régulièrement espacée
(voir "INTERPOLATION"). Les tables du géoïde suivantes sont disponibles (cependant, certaines peuvent ne pas être
installée):
erreur bilinéaire erreur cubique
nom grille du géoïde max rms max rms
egm84-30 EGM84 30' 1.546 m 70 mm 0.274 m 14 mm
egm84-15 EGM84 15' 0.413 m 18 mm 0.021 m 1.2 mm
egm96-15 EGM96 15' 1.152 m 40 mm 0.169 m 7.0 mm
egm96-5 EGM96 5' 0.140 m 4.6 mm 0032 m 0.7 mm
egm2008-5 EGM2008 5' 0.478 m 12 mm 0.294 m 4.5 mm
egm2008-2_5 EGM2008 2.5' 0.135 m 3.2 mm 0.031 m 0.8 mm
egm2008-1 EGM2008 1' 0.025 m 0.8 mm 0022 m 0.7 mm
Par défaut, le géoïde "egm96-5" est utilisé. Cela peut changer en définissant l'environnement
variable "GEOGRAPHICLIB_GEOID_NAME" ou avec le -n option. Les erreurs répertoriées ici sont
estimations des erreurs de quantification et d'interpolation dans les hauteurs signalées par rapport à
le géoïde spécifié.
Les données du géoïde seront chargées à partir d'un répertoire spécifié au moment de la compilation. Ceci peut
modifié en définissant les variables d'environnement "GEOGRAPHICLIB_GEOID_PATH" ou
"GEOGRAPHICLIB_DATA", ou avec le -d option. La -h l'option imprime le chemin du géoïde par défaut
et nom. Utilisez le -v option pour connaître le nom de chemin complet du fichier de données.
Les instructions pour télécharger et installer les données du géoïde sont disponibles sur
<http://geographiclib.sf.net/html/geoid.html#geoidinst>.
REMARQUE: tous les géoïdes ci-dessus s'appliquent à l'ellipsoïde WGS84 (a = 6378137 m, f =
1/298.257223563) uniquement.
INTERPOLATION
L'interpolation cubique est utilisée pour calculer la hauteur du géoïde à moins que -l est spécifié dans lequel
cas, une interpolation bilinéaire est utilisée. L'interpolation cubique est basée sur les moindres carrés
ajustement d'un polynôme cubique à un pochoir à 12 points
. 1 1 .
1 2 2 1 XNUMX
1 2 2 1 XNUMX
. 1 1 .
Le cube est contraint d'être indépendant de la longitude lors de l'évaluation de la hauteur à un
des pôles. L'interpolation cubique est considérablement plus précise que bilinéaire ; Cependant, il
entraîne de petites discontinuités dans la hauteur renvoyée sur les limites des cellules.
CACHE
Par défaut, le fichier de données est lu aléatoirement pour calculer les hauteurs du géoïde en entrée
postes. Cela suffit généralement pour une utilisation interactive. Si de nombreuses hauteurs doivent être
calculé, utiliser -c sud ouest nord est prévenir ÉvalGéoïde lire un rectangle de données
dans la mémoire ; les hauteurs à l'intérieur de ce rectangle peuvent alors être calculées sans aucun disque
accès. Si -a est spécifié que toutes les données du géoïde sont lues ; dans le cas de "egm2008-1", ce
nécessite environ 0.5 Go de RAM. L'évaluation des hauteurs en dehors de la zone cachée provoque
les données nécessaires à lire à partir du disque. Utilisez le -v possibilité de vérifier la taille du
cache
Indépendamment du fait qu'un cache soit demandé (avec le -a or -c options), les données de
dernière cellule de la grille en cache. Cela permet à la hauteur du géoïde le long d'un chemin continu d'être
retourné avec peu de surcharge de disque.
ENVIRONNEMENT
GEOGRAPHICLIB_GEOID_NAME
Remplacez le nom de géoïde par défaut au moment de la compilation de "egm96-5". Les -h l'option signale le
valeur de GEOGRAPHICLIB_GEOID_NAME, s'il est défini, sinon il rapporte le temps de compilation
valeur. Si la -n prénom l'option est utilisée, alors prénom a la priorité.
GEOGRAPHICLIB_GEOID_PATH
Remplacez le chemin du géoïde par défaut au moment de la compilation. C'est typiquement
"/usr/local/share/GeographicLib/geoids" sur les systèmes de type Unix et
"C:/ProgramData/GeographicLib/geoids" sur les systèmes Windows. Les -h l'option signale le
valeur de GEOGRAPHICLIB_GEOID_PATH, s'il est défini, sinon il rapporte le temps de compilation
valeur. Si la -d dir l'option est utilisée, alors dir a la priorité.
GEOGRAPHICLIB_DATA
Une autre façon de remplacer le chemin du géoïde par défaut au moment de la compilation. S'il est défini (et si
GEOGRAPHICLIB_GEOID_PATH n'est pas défini), alors $GEOGRAPHICLIB_DATA/geoids est utilisé.
LES ERREURS
Une ligne d'entrée illégale imprimera un message d'erreur sur la sortie standard commençant par
« ERREUR : » et les causes ÉvalGéoïde pour renvoyer un code de sortie de 1. Cependant, une erreur ne
causer ÉvalGéoïde Terminer; les lignes suivantes seront converties.
ABRÉVIATIONS
Le géoïde est généralement approximé par un « modèle de gravité terrestre ». Les modèles publiés par le
Les NGA sont :
EGM84
Un modèle de gravité terrestre publié par la NGA en 1984,
<http://earth-info.nga.mil/GandG/wgs84/gravitymod/wgs84_180/wgs84_180.html>.
EGM96
Un modèle de gravité terrestre publié par la NGA en 1996,
<http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm96/egm96.html>.
EGM2008
Un modèle de gravité terrestre publié par la NGA en 2008,
<http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm2008>.
WGS84
Système géodésique mondial 1984, .
AOH Hauteur au-dessus de l'ellipsoïde WGS84.
MSL Niveau moyen de la mer, utilisé comme raccourci pratique pour le géoïde. (Cependant, généralement,
le géoïde diffère de quelques mètres du niveau moyen de la mer.)
EXEMPLES
La hauteur du géoïde EGM96 à Tombouctou
écho 16:46:33N 3:00:34W | ÉvalGéoïde
=> 28.7068 -0.02e-6 -1.73e-6
Le premier nombre renvoyé est la hauteur du géoïde et les 2e et 3e sont ses pentes en
les directions nord et est.
Convertir un point dans la zone UTM 18n de MSL en HAE
écho 531595 4468135 23 | GeoidEval --msltohae -z 18n
=> 531595 4468135 -10.842
Utilisez GeoidEval en ligne en utilisant les services onworks.net