matlab - Calculate great circle distance of all matrix elements to all matrix elements -

-

i have following vector of coordinates

ll = [lat lon]; % 5000 x 2

and want calculate distances of elements elements using:

[dist, ~] = deg2km(distance[lat1, lon1, lat2, lon2]);

so final matrix be, guess, 5000x5000. btw function distance comes in mapping toolbox. tried loop one-by-one, takes centuries run.

update

ll =      -79.49583374               40.029166991     -79.50416707               40.029166991     -79.50416707               40.037500324  d = pdist(ll, @(x2,x1)(deg2km(distance(x1(1),x1(2),x2(1),x2(2))))) d =           2014.58795578131          2014.58795578131         0.168797611011563

while first 2 coordinates:

d = deg2km(distance(ll(1,2),ll(1,1),ll(2,2),ll(2,1))) d =          0.709531880098433

----

any please?

i didn't realise distance matlab function, sorry. try this:

d = pdist(ll, @(x1,x2) (distance(x1(:,1),x1(:,2),x2(:,1),x2(:,2))))

what @(x1,x2) (distance(x1(:,1),x1(:,2),x2(:,1),x2(:,2))) do?

it's anonymous function expects 2 inputs, x1 , x2, , each should n-by-2. let's pretend n equals 1. pdist give our anonymous function every different pair of rows ll inputs x1 , x2. row ll of form [lat, lon] , getting 2 such rows. know distance requires distance(lat1,lon1,lat2,lon2). if x1 row of ll it's [lat1, lon1] x1(1) lat1 , x1(2) lon2. same x2 , [lat2,lon2].

so our anonymous function doing converting distance function taking 4 inputs, more standard formed matlab distance function accepts 2 inputs.

edit:

as swapping latitudes , longitudes, might work you:

d = pdist(ll, @(x1,x2) (distance(x1(:,2),x1(:,1),x2(:,2),x2(:,1))))

Comments

Post a Comment

Popular posts from this blog

Java 3D LWJGL collision -

-
i making 3d java game lwjgl library, , wondering how add collision detection, player not go through models. i using obj models. here objloader class, loads models: package renderengine; import java.io.bufferedreader; import java.io.file; import java.io.filenotfoundexception; import java.io.filereader; import java.util.arraylist; import java.util.list; import models.rawmodel; import org.lwjgl.util.vector.vector2f; import org.lwjgl.util.vector.vector3f; public class objloader { public static rawmodel loadobjmodel(string filename, loader loader){ filereader fr = null; try { fr = new filereader(new file("res/"+filename+".obj")); } catch (filenotfoundexception e) { system.err.println("couldn't load file!"); e.printstacktrace(); } bufferedreader reader = new bufferedreader(fr); string line; list<vector3f> vertices = new arraylist<vector3f...
Read more

methods - python can't use function in submodule -

-
my folder structure this: pythonstuff/ program.py modulea/ __init__.py foo.py bar.py here's code bar.py : def hello(): print("hello, world!") here's code program.py : #!/usr/bin/env python3 modulea import bar bar.hello() i run $ python3 program.py somehow error: file "program.py", line 3, in <module> bar.hello() attributeerror: 'module' object has no attribute 'hello' edit: __init__.py file empty. edit2: after trying realized had bar.py in root directory contained hello() method. bar.py in modulea/ directory empty. add it import os __all__ = [] module in os.listdir(os.path.dirname(__file__)): if module != '__init__.py' , module[-3:] == '.py': __all__.append(module[:-3]) the init .py files required make python treat directories containing packages; done prevent directories common name, such string, unintentionally hiding v...
Read more

c# - ErrorThe type or namespace name 'AxWMPLib' could not be found (are you missing a using directive or an assembly reference?) -

-
what possible reference handler axwmplib._wmpocxevents_playstatechangeeventhandler ? when pass event in method, private void axwindowsmediaplayer1_playstatechange(object sender,axwmplib._wmpocxevents_playstatechangeevent e) compiling throws following error: the type or namespace name 'axwmplib' not found (are missing using directive or assembly reference?) if using visual studio make sure have axwmplib listed under references in solution explorer. if have embedded windows media player control in form, should have been added automatically along wmplib . an easy way add these libraries if missing add windows media player control form (an existing 1 do). right-click on toolbox , pick "choose items..." context menu, tick "windows media player" control on "com components" tab , click ok. should have little windows media player icon attached cursor; click on form add it. add axwmplib , wmplib references solution. can delete w...
Read more