Sponsored Content
Top Forums Shell Programming and Scripting Calculate distance and azimuth Post 302395369 by binlib on Monday 15th of February 2010 09:59:43 PM
Old 02-15-2010
I rearranged your formulas and made it awk executable. Not sure if it's right:

Code:
$ cat sphere.awk
function pow(b, p)
{
  return b^p
}

{
  M_PI = 4*atan2(1,1)
  rad = 180 / M_PI
  lat1 = $1 / rad
  lat2 = $3 / rad
  lon1 = $2 / rad
  lon2 = $4 / rad
  latdif = lat1 - lat2
  londif = lon1 - lon2
  meanlat = (lat1 + lat2)/2

  a = 6377276.345
  b = 6356075.4131
  e = sqrt(((a*a)-(b*b))/(a*a))
  mrcurt = (a*(1-(e*e)))/pow((1-((e*e)*pow(sin(meanlat),2))),1.5)
  prcurt =a/sqrt(1-pow((e*sin(meanlat)),2))
  distance=sqrt((pow(londif*prcurt*cos(meanlat),2))+(pow((latdif*mrcurt),2)))

  A=2*atan2(londif*((prcurt/mrcurt)*(cos(meanlat))),latdif)
  B=londif*(sin(meanlat))
  Az=(A-B)/2
  if(londif>0&&latdif<0) Az=Az+M_PI
  if(londif<0&&latdif<0) Az=Az+M_PI
  if(londif<0&&latdif>0) Az=Az+(2*M_PI)

  printf("%8.0f %6.2f\n", distance, Az*rad)
}

$ cat sphere.dat
12.000 25.125 14.235 25.012
14.200 81.000 25.584 25.014
45.023 25.365 25.152 35.222

$ awk -f sphere.awk sphere.dat
  247544 357.15
 5995861 272.61
 2380475 340.65

 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Lat/Long Distance Calculation

I amtrying to write a script that would compute the distance between an "x" number of points. This is what I have come up with so far and it is not working. Can anyone modify it to make it work? A=34.16597 B=-84.33244 C=34.2344 D=-84.29189 test "$A" -eq "$C" -o "$B" -eq "$D" then echo... (3 Replies)
Discussion started by: Ernst
3 Replies

2. UNIX for Dummies Questions & Answers

Long Distance UNIX (Solaris) Cloning ?

Need some advice and guidance for this UNIX beginner. Due to downsizing I have inherited the SysAdmin duties..(sigh). Please excuse and forgive me if I use the wrong terms below.... Situation: We have UNIX ( Solaris 7/8/9( it varies) on Sun Ultra 10's) servers located at several global... (1 Reply)
Discussion started by: HikerLT
1 Replies

3. Shell Programming and Scripting

program to calculate distance between 5 atoms

Hello, I am a beginner with perl. I have a perl program to calculate the distance between 5 atoms or more. i have an array which looks like this: 6.324 32.707 50.379 5.197 32.618 46.826 4.020 36.132 46.259 7.131 38.210 45.919 6.719 38.935 42.270 2.986 39.221 ... (1 Reply)
Discussion started by: annie_singh
1 Replies

4. Programming

Converting distance list to distance matrix in R

Hi power user, I have this type of data (distance list): file1 A B 10 B C 20 C D 50I want output like this # A B C D A 0 10 30 80 B 10 0 20 70 C 30 20 0 50 D 80 70 50 0 Which is a distance matrix I have tried... (0 Replies)
Discussion started by: anjas
0 Replies

5. Shell Programming and Scripting

Calculating distance between two LAT long coordinates

hi, i have a pair of latitude and longitude and i want to calculate the distance between these two points. In vbscript i achieved in the following way...Now i want to implement this in unix shell scripting.... <% Dim lat1, lon1, lat2, lon2 const pi = 3.14159265358979323846 ... (8 Replies)
Discussion started by: aemunathan
8 Replies

6. UNIX for Dummies Questions & Answers

How to make a distance matrix

Hi, I'm trying to generate a distance matrix between sample pairs for use in a tree-drawing program (example below). The example below demonstrates what I'd like to get out of the data - essentially, to calculate the proportion of positions where two samples differ. Any help much appreciated!... (1 Reply)
Discussion started by: auburn
1 Replies

7. Shell Programming and Scripting

Calculate age of a file | calculate time difference

Hello, I'm trying to create a shell script (#!/bin/sh) which should tell me the age of a file in minutes... I have a process, which delivers me all 15 minutes a new file and I want to have a monitoring script, which sends me an email, if the present file is older than 20 minutes. To do... (10 Replies)
Discussion started by: worm
10 Replies

8. Shell Programming and Scripting

finding distance between numbers

Hi, I have a file as ABC 1634230,1634284,1634349,1634468 1634272,1634301,1634356,1634534 What I want is to find distance between the numbers.. column 1 is the gene name and column 2 are starts and column 3 are their respective stops for the starts. So what I want is column 3 which has +1... (2 Replies)
Discussion started by: Diya123
2 Replies

9. Shell Programming and Scripting

Edit distance using perl or awk

Dear all, I am working on a large Sindhi lexicon which I hope to complete by 2017 and place in open source. The database is in Arabic script in two columns delimited by an equal to sign. Column 1 contains a word or words without the short vowel and also some extraneous information which is... (0 Replies)
Discussion started by: gimley
0 Replies

10. Shell Programming and Scripting

Calculate average, azimut and distance

Gents, Please i will to get the distance and azimut from 2 coordinates: Usig excel formula i get the correct values, but i will like to do it using awk. Example A 35089.0 50345.016 9 75 1 2101774 77 70 79 483911.6 2380106.9 137.4 1 1 6 1 A 35089.0 50345.01620 75... (8 Replies)
Discussion started by: jiam912
8 Replies
Geo::GoogleEarth::Pluggable::Plugin::Default(3pm)	User Contributed Perl Documentation	 Geo::GoogleEarth::Pluggable::Plugin::Default(3pm)

NAME
Geo::GoogleEarth::Pluggable::Plugin::Default - Geo::GoogleEarth::Pluggable Default Plugin Methods METHODS
Methods in this package are AUTOLOADed into the Geo::GoogleEarth::Pluggable::Folder namespace at runtime. CONVENTIONS
Plugin Naming Convention: Geo::GoogleEarth::Pluggable::Plugin::CPANID (e.g. "MRDVT") Object Naming Convention: Geo::GoogleEarth::Pluggable::Contrib::"$method" (e.g. Point, CircleByCenterPoint) You only need to have one plugin pointing to all of your contributed objects. The package should be named after the plugin not the objects since there is a many to one relationship. (e.g. Geo-GoogleEarth-Pluggable- Plugin-MRDVT) Point Constructs a new Placemark Point object and appends it to the parent folder object. Returns the object reference if you need to make any setting changes after construction. my $point=$folder->Point(name=>"My Placemark", lat=>38.897607, lon=>-77.036554, alt=>0); LineString $folder->LineString(name=>"My Placemark", coordinates=>[ [lat,lon,alt], {lat=>$lat,lon=>$lon,alt=>$alt}, ]); LinearRing $folder->LinearRing(name=>"My Placemark", coordinates=>[ [lat,lon,alt], {lat=>$lat,lon=>$lon,alt=>$alt}, ]); TODO
Need to determine what methods should be in the Folder package and what should be on the Plugin/Default package and why. BUGS
Please log on RT and send to the geo-perl email list. SUPPORT
DavisNetworks.com supports all Perl applications including this package. AUTHOR
Michael R. Davis (mrdvt92) CPAN ID: MRDVT COPYRIGHT
This program is free software licensed under the... The BSD License The full text of the license can be found in the LICENSE file included with this module. SEE ALSO
Geo::GoogleEarth::Pluggable::Contrib::Point, Geo::GoogleEarth::Pluggable::Contrib::LineString, Geo::GoogleEarth::Pluggable::Contrib::LinearRing perl v5.14.2 2010-12-27 Geo::GoogleEarth::Pluggable::Plugin::Default(3pm)
All times are GMT -4. The time now is 10:56 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy