clear
close all

load GLpoints.dat;
GLP = GLpoints;

load MFcoordinates.dat;
MFP=MFcoordinates;
MFN=length(MFP(:,1));


load GoCcoordinates.sorted.dat
GoCP=GoCcoordinates_sorted;
GoCN=length(GoCP(:,1));


load GCcoordinates.sorted.dat
GCP=GCcoordinates_sorted;
GCN=length(GCP(:,1));

%%% Plot the positions of MFs, GCs and GoCs.

figure() 

plot3(MFP(:,1),MFP(:,2),ones(MFN,1),'.y','MarkerSize',12),hold on;
plot3(GLP(:,2),GLP(:,3),GLP(:,4),'.m','MarkerSize',4),hold on;
plot3(GoCP(:,1),GoCP(:,2),GoCP(:,3),'.r','MarkerSize',12),hold on;
plot3(GCP(:,1),GCP(:,2),GCP(:,3),'.b','MarkerSize',4);

axis([0, 1500, 0, 1000, 0, 200]);