Sorry your browser is not supported!

Give PLEASE LEARN HOW distance from point to point on the Earth\\\\\\\'s surface
If the point given in the FORMAT

Latitude:
Denomination - the degree, minutes & seconds (hDDD, MM, SS.sssss) OR
Packed - small second-degrees (hDDDMMSS.sssss) OR
Decimal number - degree (hDDD.ddddddddd) OR

Longitude:
Denomination - the degree, minutes & seconds (hDDD, MM, SS.sssss) OR
Packed - small second-degrees (hDDDMMSS.sssss) OR
Decimal number - degree (hDDD.ddddddddd) OR

A RESULT OF NEED IN METRE

forall

Back to top

Profile PM Website

TheComet

17

Years of Service

User Offline

Joined: 18th Oct 2007

Location: I`m under ur bridge eating ur goatz.

Posted: 4th Aug 2009 17:27

Sorry, I don't understand...

Please elaborate.

TheComet

Make the path of your enemies easier with Waypoint Pro!

Back to top

Profile PM Email Website

NeX the Fairly Fast Ferret

20

Years of Service

User Offline

Joined: 10th Apr 2005

Location: The Fifth Plane of Oblivion

Posted: 4th Aug 2009 18:00

He wants to know how to figure out the distance in metres between two points given in some packed format resembling longitude and latitude.

Sorry, can't help you there.

Back to top

Profile PM Website

Phaelax

DBPro Master

22

Years of Service

User Offline

Joined: 16th Apr 2003

Location: Metropia

Posted: 4th Aug 2009 20:17

The "great-circle" distance. Some use the haversine formula.
Here's what I used recently for a website I'm doing which requires me to find the distances between zip codes(given I have the zip code's long and lat.

R = 6371 `km radius of earth
dLat = toRad(lat2-lat1)
dLon = toRad(lon2-lon1) 
a = sin(dLat/2) * sin(dLat/2) +
        cos(toRad(lat1)) * cos(toRad(lat2)) * 
        sin(dLon/2) * sin(dLon/2); 
c = 2 * atanfull(sqrt(a), sqrt(1-a)); 
d = R * c;

[url="http://dbcc.zimnox.com"]

[/url]

Back to top

Profile PM Email Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 4th Aug 2009 21:03

THANKS Phaelax
toRad() - WHAT IS FUNCTION ?

forall

Back to top

Profile PM Website

Phaelax

DBPro Master

22

Years of Service

User Offline

Joined: 16th Apr 2003

Location: Metropia

Posted: 4th Aug 2009 21:12

DB doesn't have one, it just converts from degrees into radians. A quick google search will show you the formula.

[url="http://dbcc.zimnox.com"]

[/url]

Back to top

Profile PM Email Website

kaedroho

17

Years of Service

User Offline

Joined: 21st Aug 2007

Location: Oxford,UK

Posted: 4th Aug 2009 21:16

You dont need to convert to radians, DBs trig commands use degrees.

Back to top

Profile PM Email Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 4th Aug 2009 21:39 Edited at: 4th Aug 2009 23:25

Thank you very much once again - very helpful


R# = 6371.0 `km radius of earth

lat1# = 47.93049125
lat2# = 47.91645709
lon2# = 39.81548309
lon1# = 39.81548309
dLat# =TORAD(lat2#-lat1#)
dLon# =TORAD(lon2#-lon1#)

a# = sin(dLat#/2.0) * sin(dLat#/2.0) + cos(TORAD(lat1#)) * cos(TORAD(lat2#)) * sin(dLon#/2) * sin(dLon#/2.0)
c# = 2.0 * atanfull(sqrt(a#), sqrt(1.0-a#))
d# = R# * c#


DO

TEXT 100,100,STR$(D#)

SYNC
LOOP

FUNCTION TORAD(SS#)
radians# = ( pi() / 180.0 ) * SS#
ENDFUNCTION   radians#

forall

Back to top

Profile PM Website

kaedroho

17

Years of Service

User Offline

Joined: 21st Aug 2007

Location: Oxford,UK

Posted: 6th Aug 2009 18:18

You dont need to convert to radians! DBPro uses degrees!

Back to top

Profile PM Email Website

prasoc

16

Years of Service

User Offline

Joined: 8th Oct 2008

Location:

Posted: 6th Aug 2009 21:24

He wants one: radians = (degree) * (pi/180.0f)

Your signature has been erased by a mod

Back to top

Profile PM Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 7th Aug 2009 00:24

I solved this problem.
required plug STYX for DBpro as atan (), tan (), atanfull () not working correctly
Here's the code:

REM Проект: dist
REM Создан: 06.08.2009 20:56:08
REM
REM ***** Главный исходный файл *****
REM
l1#= 40
lo1#= 0
l2#= -40
lo2#= 0

pia# = 3.14159265358979
rad# = 6378
`'получение координат точек в радианах
lat1# = l1#
lat2# = l2#
long1# = lo1#
long2# = lo2#
`'косинусы и синусы широт и разниц долгот
cl1# =cos( lat1# )
cl2# =cos( lat2# )
sl1# =sin( lat1# )
sl2# =sin( lat2# )
delta# = long2# - long1#
cdelta# =cos( delta# )
sdelta# =sin( delta# )
`'вычисления длины большого круга
p1# = (cl2#*sdelta#)^2.0
p2# = ((cl1#*sl2#) - (sl1#*cl2#*cdelta#))^2.0
p3# = (p1# + p2#)^0.5
p4# = sl1#*sl2#
p5# = cl1#*cl2#*cdelta#
p6# = p4# + p5#
p7# = p3#/p6#
ADD EXPRESSION VAR "aa",p7#
expr$ ="aTAN( aa )"
anglerad#  = EVALUATE EXPRESSION( expr$ )
dist# = anglerad#*rad#

do
         text 0,0,"distance = " + str$(dist#)
         text 20, 50 ,"lat1# = "+ str$(lat1#  )
         text 20, 70 ,"lat2# = "+ str$(lat2#  )
         text 20, 90 ,"long1# = "+ str$(long1#  )
         text 20, 110 ,"long2# = "+ str$(long2#  )
         text 20, 130 ,"cos( lat1# ) = "+ str$( cl1# )
         text 20, 150 ,"cos( lat2# ) = "+ str$( cl2# )
         text 20, 170 ,"sin( lat1# ) = "+ str$( sl1# )
         text 20, 190 ,"sin( lat2# ) = "+ str$( sl2# )

sync
loop

REM Проект: dist
REM Создан: 06.08.2009 20:56:08
REM
REM ***** Главный исходный файл *****
REM
l1#= 40
lo1#= 0
l2#= -40
lo2#= 0

pia# = 3.14159265358979
rad# = 6378
`'получение координат точек в радианах
lat1# = l1#
lat2# = l2#
long1# = lo1#
long2# = lo2#
`'косинусы и синусы широт и разниц долгот
cl1# =cos( lat1# )
cl2# =cos( lat2# )
sl1# =sin( lat1# )
sl2# =sin( lat2# )
delta# = long2# - long1#
cdelta# =cos( delta# )
sdelta# =sin( delta# )
`'вычисления длины большого круга
p1# = (cl2#*sdelta#)^2.0
p2# = ((cl1#*sl2#) - (sl1#*cl2#*cdelta#))^2.0
p3# = (p1# + p2#)^0.5
p4# = sl1#*sl2#
p5# = cl1#*cl2#*cdelta#
p6# = p4# + p5#
p7# = p3#/p6#
ADD EXPRESSION VAR "aa",p7#
expr$ ="aTAN( aa )"
anglerad#  = EVALUATE EXPRESSION( expr$ )
dist# = anglerad#*rad#

do
         text 0,0,"distance = " + str$(dist#)
         text 20, 50 ,"lat1# = "+ str$(lat1#  )
         text 20, 70 ,"lat2# = "+ str$(lat2#  )
         text 20, 90 ,"long1# = "+ str$(long1#  )
         text 20, 110 ,"long2# = "+ str$(long2#  )
         text 20, 130 ,"cos( lat1# ) = "+ str$( cl1# )
         text 20, 150 ,"cos( lat2# ) = "+ str$( cl2# )
         text 20, 170 ,"sin( lat1# ) = "+ str$( sl1# )
         text 20, 190 ,"sin( lat2# ) = "+ str$( sl2# )

sync
loop

forall

Back to top

Profile PM Website

Mahoney

17

Years of Service

User Offline

Joined: 14th Apr 2008

Location: The Interwebs

Posted: 7th Aug 2009 04:30

Apollo: Perfect English NOT requiring.

Windows Vista Home Premium Intel Pentium Dual-Core 1.6 Ghz 1GB DDR2 RAM GeForce 8600GT Twin Turbo

Back to top

Profile PM

kaedroho

17

Years of Service

User Offline

Joined: 21st Aug 2007

Location: Oxford,UK

Posted: 7th Aug 2009 18:19 Edited at: 7th Aug 2009 18:20

Quote: "atan (), tan (), atanfull () not working correctly "

For the third time. THEY USE DEGREES!!

Styx is probably working because it uses radians.

DBPro trig commands are 100s of times faster than styxs/

this code should work

R# = 6371.0 `km radius of earth

lat1# = 47.93049125
lat2# = 47.91645709
lon2# = 39.81548309
lon1# = 39.81548309
dLat# =lat2#-lat1#
dLon# =lon2#-lon1#

a# = sin(dLat#/2.0) * sin(dLat#/2.0) + cos(lat1#) * cos(lat2#) * sin(dLon#/2) * sin(dLon#/2.0)
c# = 2.0 * atanfull(sqrt(a#), sqrt(1.0-a#))
d# = R# * c#


DO

TEXT 100,100,STR$(D#)

SYNC
LOOP

Back to top

Profile PM Email Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 7th Aug 2009 21:22

R# = 6371.0 `km radius of earth

lat1# = 47.93049125
lat2# = 47.91645709
lon2# = 39.81548309
lon1# = 39.81548309
dLat# =lat2#-lat1#
dLon# =lon2#-lon1#

a# = sin(dLat#/2.0) * sin(dLat#/2.0) + cos(lat1#) * cos(lat2#) * sin(dLon#/2) * sin(dLon#/2.0)
c# = 2.0 * atanfull(sqrt(a#), sqrt(1.0-a#))
d# = R# * c#


DO

TEXT 100,100,STR$(D#)

SYNC
LOOP

this is not working
atanfull() is a problem

different results

forall

Back to top

Profile PM Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 7th Aug 2009 21:45

lat1# = 0
lat2# = 0
lon2# = 45
lon1# = -45

should be 10018 km
check out? if it means I get something from DBPro

forall

Back to top

Profile PM Website

kaedroho

17

Years of Service

User Offline

Joined: 21st Aug 2007

Location: Oxford,UK

Posted: 7th Aug 2009 21:51

Alright.

Try your code with DBPros atan command.

DBPros atan command is alot faster as styxes gets sent to styx in a string, and needs to be ran through loads of code just to find out what you want it to do, then finds the variables then sends the atan command to the CPU.

DBPro compiles the atan command into machine code to be sent directly to the CPU, without having to process any strings.

Back to top

Profile PM Email Website

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 7th Aug 2009 21:57

not work.
with DBPros atan command result is -574019....

forall

Back to top

Profile PM Website

Neuro Fuzzy

17

Years of Service

User Offline

Joined: 11th Jun 2007

Location:

Posted: 7th Aug 2009 22:08 Edited at: 7th Aug 2009 22:18

are you absolutely, 100% sure that your using radians? I've used ALL of dbpro's trig commands without any problem whatsoever.

change phalaex's example to this:

R = 6371 `km radius of earth
dLat = lat2-lat1
dLon = lon2-lon1)
a = sin(dLat/2) * sin(dLat/2) +
        cos(lat1) * cos(lat2) * 
        sin(dLon/2) * sin(dLon/2); 
c = 2 * atanfull(sqrt(a), sqrt(1-a)); 
d = R * c;

[edit]

uhh, also, the example you posted does work... but the answer is in kilometers.

Back to top

Profile PM Email

OldPMan

TGC Store Seller

16

Years of Service

User Offline

Joined: 10th Aug 2008

Location:

Posted: 7th Aug 2009 22:16

sorry
this is work