Il s'agit de la commande Planimeter 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
Planimètre -- calcule l'aire des polygones géodésiques
SYNOPSIS
Planimètre [ -r ] [ -s ] [ -l ] [ -e a f ] [ -w ] [ -p prec ] [ -G | -E | -Q | -R ] [
--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
Mesurer l'aire d'un polygone géodésique. Lit les sommets des polygones à partir de l'entrée standard, un
par ligne. Les sommets peuvent être donnés sous forme de latitude et de longitude, de coordonnées UTM/UPS ou MGRS,
interprété de la même manière que GéoConvertir(1). (Les coordonnées MGRS signifient le centre du
carré MGRS correspondant.) La fin de l'entrée, une ligne vide ou une ligne qui ne peut pas être
interprété comme un sommet signale la fin d'un polygone et le début du suivant. Pour
chaque polygone imprime une ligne récapitulative avec le nombre de points, le périmètre (en mètres),
et la superficie (en mètres^2).
Les arêtes du polygone sont données par le le plus court géodésique entre des sommets consécutifs.
Dans certains cas, il peut y avoir deux ou plusieurs géodésiques les plus courtes, et dans ce cas, la
polygone n'est pas spécifié de manière unique par ses sommets. Cela ne se produit qu'avec des bords très longs
(pour l'ellipsoïde WGS84, toute arête plus courte que 19970 km est uniquement spécifiée par son extrémité
points). Dans de tels cas, insérez un sommet supplémentaire près du milieu du bord long pour
définir la limite du polygone.
Par défaut, les polygones traversés dans le sens inverse des aiguilles d'une montre renvoient une zone positive et
ceux parcourus dans le sens des aiguilles d'une montre renvoient une zone négative. Cette convention de signe est
inversé si le -r option est donnée.
Bien sûr, encercler une zone dans le sens des aiguilles d'une montre équivaut à encercler la
reste de l'ellipsoïde dans le sens inverse des aiguilles d'une montre. L'interprétation par défaut utilisée
by Planimètre est celui qui se traduit par une plus petite ampleur de la zone ; c'est-à-dire l'ampleur
de l'aire est inférieure ou égale à la moitié de l'aire totale de l'ellipsoïde. Si la -s
option est donnée, alors l'interprétation retenue est celle qui aboutit à une zone positive ;
c'est-à-dire que l'aire est positive et inférieure à l'aire totale de l'ellipsoïde.
Seuls les polygones simples (c'est-à-dire non auto-sécants) sont pris en charge pour le calcul de l'aire.
Les polygones peuvent inclure un ou les deux pôles. Il n'est pas nécessaire de fermer le polygone.
OPTIONS
-r basculer si la traversée dans le sens inverse des aiguilles d'une montre du polygone renvoie un résultat positif (le
par défaut) ou résultat négatif.
-s basculer s'il faut retourner un résultat signé (par défaut) ou non.
-l basculer si les sommets représentent un polygone (par défaut) ou une polyligne. Pour un
polyligne, le nombre de points et la longueur du chemin les joignant sont renvoyés ;
le chemin n'est pas fermé et la zone n'est pas signalée.
-e spécifier l'ellipsoïde via a f; le rayon équatorial est a et l'aplatissement est f.
Paramètres f = 0 donne une sphère. Spécifier f < 0 pour un ellipsoïde allongé. Un simple
fraction, par exemple, 1/297, est autorisé pour f. Par défaut, l'ellipsoïde WGS84 est utilisé, a =
6378137 m, f = 1/298.257223563. Si vous saisissez des sommets sous forme de coordonnées UTM/UPS ou MGRS,
utiliser l'ellipsoïde par défaut, puisque la conversion de ces coordonnées en latitude et
longitude utilise toujours les paramètres WGS84.
-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).
-p régler la précision de sortie sur prec (par défaut 6); le périmètre est donné (en mètres) avec
prec chiffres après la virgule ; la surface est donnée (en mètres^2) avec (prec - 5)
chiffres après la virgule.
-G utiliser la formulation en série pour les géodésiques. Il s'agit de l'option par défaut et
recommandé pour les applications terrestres. Cette option, -G, et les trois suivants
options, -E, -Q et -R, s'excluent mutuellement.
-E utiliser des algorithmes "exacts" (basés sur des intégrales elliptiques) pour les calculs géodésiques.
Celles-ci sont plus précises que les extensions de série (par défaut) pour |f| > 0.02. (Mais
notez que l'implémentation des zones dans GeodesicExact utilise une série d'ordre élevé et
ceci n'est exact que pour les aplatissements modestes.)
-Q effectuer le calcul sur la sphère authalique. Le calcul de l'aire est précis même
si l'aplatissement est important, à condition de les bords sont suffisamment courts. Le périmètre
le calcul n'est pas précis.
-R Les lignes joignant les sommets sont des lignes de rhumb au lieu de géodésiques.
--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. Pour un polygone donné, le dernier de ces
La chaîne trouvée sera ajoutée à la ligne de sortie (séparée par un espace).
--version
version imprimée et quitter.
-h utilisation de l'impression et sortie.
--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.
EXEMPLES
Exemple (la superficie du carré MGRS de 100 km 18SWK)
Planimètre <
18n500000
18n600000
18n600000
18n500000
EOF
=> 4 400139.53295860 10007388597.1913
Le code suivant prend la sortie de gdalinfo et rapporte la zone couverte par les données
(en supposant que les bords de l'image sont des géodésiques).
#! / Bin / sh
egrep '^((Supérieur|Inférieur) (Gauche|Droite)|Centre) ' |
sed -e 's/d /d/g' -e "s/' /'/g" | tr -s '(),\r\t' ' ' | awk '{
si ($1 $2 == "UpperLeft")
ul = 6 $ " " 5 $;
else if ($1 $2 == "LowerLeft")
ll = 6 $ " " 5 $ ;
else if ($1 $2 == "UpperRight")
ur = 6 $ " " 5 $;
else if ($1 $2 == "LowerRight")
lr = 6 $ " " 5 $;
else if ($1 == "Centre") {
printf "%s\n%s\n%s\n%s\n\n", ul, ll, lr, ur;
ul = ll = ur = lr = "" ;
}
}
' | Planimètre | couper -f3 -d' '
Utilisez Planimeter en ligne en utilisant les services onworks.net