Find Nearest Neighbors on Sphere
- 6 files
- 25 downloads
A geodesic grid is an even grid over the surface of a sphere. The algorithm is optimized for a grid generated by GridSphere and won't work on an arbitrary geodesic grid. All code is compatible with GNU Octave. The GridSphere function is available on the MATLAB file exchange at http://www.mathworks.com/matlabcentral/fileexchange/28842-gridsphere.
FindNearestNeighbors and GridSphere share some functions in common. Each package contains a copy of these functions so that both can stand alone. To eliminate duplicates, simply move all the files into a single folder and replace the shared files when prompted.
Cite As
Kurt von Laven (2020). Find Nearest Neighbors on Sphere (https://www.mathworks.com/matlabcentral/fileexchange/28844-find-nearest-neighbors-on-sphere), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (5)
Updates
| 1.7 | Added a link to another FX submission that this submission uses code from. |
|
| 1.6.0.0 | Made minor corrections to comments. Removed write and execute permissions from group/others on all *.m files. Removed trailing whitespace and switched from 2 spaces after each period to 1 in the interest of readability. |
|
| 1.5.0.0 | Made code available as a toolbox and removed some unnecessary hidden files. |
|
| 1.4.0.0 | Replaced tabs with spaces so that the source code displays consistently in all text editors. |
|
| 1.3.0.0 | Replaced the ElementWiseMax function with an equivalent built-in function. |
|
| 1.2.0.0 | Added missing Logarithm function file. |
|
| 1.1.0.0 | Added missing function file. |
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
FindNearestNeighbors/
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Jonathan Wolf (view profile)
Smart code, well documented and straightforward to understand.
Only one thing that could be improved: The code does not work reliably for a high density of points (e.g. for a very dense mesh): The reason is that NUM_LINES_OF_LAT_TO_SEARCH in FindNearestNeighborsWithinNearbyLatitudes is too small (by default 3). It is also not possible to just increase it until the neighbor search works reliably because one then easily gets indexes that are out of bounds. Had to add some if statements to make it work for any point density and automatically determine a useful NUM_LINES_OF_LAT_TO_SEARCH.
Brendan Carter (view profile)
_
Brendan Carter (view profile)
Fast and found it easy to use. Thanks!
Paul (view profile)
-- Wrong page. Comment intended for GridSphere code.
Paul (view profile)
Blazingly fast and well documented. Only an example would not have been unbecoming, e.g. to illustrate the fact that returned values are in degrees (whereas Matlab's sph2cart requires radians).
[lat,long] = GridSphere(1000);
latrad = deg2rad(lat);
longrad = deg2rad(long);
[x,y,z] = sph2cart(longrad, latrad, 1);
scatter3(x, y, z, 22);
axis equal vis3d;