Read in text file, plot, calculate nearest neighbors, and classify based on distance

nndist_layer_code.m

spacinglay=importdata('spacing_layer_grid.csv',',',1);
disp(spacinglay.colheaders)
%% parse the data in matlab
% A tiny piece of code to find which column the given header is in
findColumn=@(idstrut,colName) find(not(cellfun('isempty',strfind(idstrut.colheaders,['"' colName '"']))));
% a function to return all of the data in that column
getColumn=@(idstrut,colName) idstrut.data(:,findColumn(idstrut,colName));

%% calculate the nearest neighbor distance and angle
[nndist,nnangle]=findnn(getColumn(spacinglay,'x'),getColumn(spacinglay,'y'));


%% show the data
figure(1)
subplot(1,2,1);
plot(getColumn(spacinglay,'x'),getColumn(spacinglay,'y'),'r.')
title('Point Location')
subplot(1,2,2);
hist(nndist)
pause(0.5)

%% make two groups based on volume
groupId=kmeans(nndist',2,'EmptyAction','singleton')
% show the groups as different colors
xdata=getColumn(spacinglay,'x');
ydata=getColumn(spacinglay,'y');
figure(2)
plot(xdata(groupId==1),ydata(groupId==1),'r.')
hold on;
plot(xdata(groupId==2),ydata(groupId==2),'bo')
legend({'Group 1','Group 2'});
hold off;