diff --git a/sort_dist_angle.m b/sort_dist_angle.m
new file mode 100644
index 0000000000000000000000000000000000000000..40d9910752062bed290b7aa0cfb8e75dd42771a0
--- /dev/null
+++ b/sort_dist_angle.m
@@ -0,0 +1,29 @@
+function [sig,order] = sort_dist_angle(dist,angle)
+%SORT_DIST_ANGLE Summary of this function goes here
+%   Detailed explanation goes here
+
+sz = size(dist);
+
+[~,sig] = sort(dist(:)','ascend');
+
+dist  =  dist(sig);
+angle = angle(sig);
+
+for cur_dist = unique(dist)
+    cur_ind = dist==cur_dist;
+    
+      sig_cur_ind =   sig(cur_ind);
+     dist_cur_ind =  dist(cur_ind);
+    angle_cur_ind = angle(cur_ind);
+    
+    [~,cur_sig] = sort(angle_cur_ind,'ascend');
+    
+      sig(cur_ind) =   sig_cur_ind(cur_sig);
+     dist(cur_ind) =  dist_cur_ind(cur_sig);
+    angle(cur_ind) = angle_cur_ind(cur_sig);
+end
+
+order = zeros(sz);
+order(sig) = 1:numel(dist);
+
+end
\ No newline at end of file