
function markers = marker_generator(space,tot_markers,error)
%OUTPUT: 
%-  markers = vettore con tutti i marker generati casualmente scegliendo 
%   un numero compreso tra il numero medio di marker visti per telecamera 
%   (=50)e il  n° medio di marker totali(=1500)  (per ora uso solo 150)
%   numero massimo per telecamera

if nargin < 2
    media_n_markers=50;
    max_n_markers=150;
    tot_markers = randint(1,1,[media_n_markers,max_n_markers]);
end

% dovo considereare anche l'incertezza
 if nargin<3
    dx=0.07;
    dy=0.07;
    dz=0.07;
 else
    dx=error(1);
    dy=error(2);
    dz=error(3);
 end

latox=space(1);
latoy=space(2);
latoz=space(3);
    
markers(1,:)= latox*rand(1,tot_markers);
markers(2,:)= latoy*rand(1,tot_markers);
markers(3,:)= latoz*rand(1,tot_markers);

%% plot dei marker con relative incertezze
%% se all'inizio del main mettete dx,dy,dz molto piccoli in relazione
%% ai lati del cubo di visione totale (a=latox,b=latoy etc.) noterete
%% che è inutile calcolare la probab che un voxel contenga parte
%% dell'incertezza....
% figure(11)
% ax = axes('XLim',[0 latox],'YLim',[0 latoy],'ZLim',[0 latoz]);
% view(3); grid on; axis equal
% for i=1:size(markers,2)
%     voxel(markers(:,i),[dx dy dz],'r',1)
%     hold on
% end
% figure(11)
