function vData = GenerateSampleCoronaData( cNumberOfPoints, cInnerRadius, cOuterRadius, cNoiseSize )

%
% function vData = GenerateSampleCoronaData( cNumberOfPoints, cInnerRadius, cOuterRadius, cNoiseSize )
%
 
vData = zeros(cNumberOfPoints,2);

% Create the blob in the middle
lHalfPoints = floor(cNumberOfPoints/2);
vData(1:lHalfPoints,:) = cInnerRadius*randn(lHalfPoints,2);

% Create the circle around
lRandRadia = randn(1,cNumberOfPoints-lHalfPoints);
for lk = 1:cNumberOfPoints-lHalfPoints
    lNewPointRadius = cOuterRadius+ ( cNoiseSize*lRandRadia(lk) );
    vData(lHalfPoints+lk,1) = lNewPointRadius * cos(2*pi*lk/(cNumberOfPoints-lHalfPoints));
    vData(lHalfPoints+lk,2) = lNewPointRadius * sin(2*pi*lk/(cNumberOfPoints-lHalfPoints));    
end;

return;
