From 321efec12b98fe1506b13bdc2dce31199b7d5794 Mon Sep 17 00:00:00 2001
From: ralfHielscher <ralf.hielscher@gmail.com>
Date: Sat, 4 Jun 2022 23:13:25 +0200
Subject: [PATCH] remove duplicated vertices at quadruple points

---
 EBSDAnalysis/@EBSD/private/spatialDecomposition.m | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/EBSDAnalysis/@EBSD/private/spatialDecomposition.m b/EBSDAnalysis/@EBSD/private/spatialDecomposition.m
index c5c086de8..1074f6038 100644
--- a/EBSDAnalysis/@EBSD/private/spatialDecomposition.m
+++ b/EBSDAnalysis/@EBSD/private/spatialDecomposition.m
@@ -35,6 +35,18 @@
     [V,D] = voronoin([X;dummyCoordinates],{'Q5','Q6','Qs'}); %,'QbB'
             
   end
+
+  %if check_option(varargin,'quadruplePoints')
+
+  [V,~,ic] = unique(V,'rows');
+  %D = cellfun(@(x) ic(x).',D,'UniformOutput',false);
+
+  for k = 1:length(D)
+    x = ic(D{k}).';              % merge points that coincide
+    D{k} = x(diff([x,x(1)])~=0); % remove dubplicates in D
+  end
+  
+  %end
     
   D = D(1:size(X,1));