Optic_flow_RF_estimation_Zhang_et_al_2022/code/General/sphDensity.m

function [sample_pnt,sampleDen] = sphDensity(sphpoint,sample_pntNum,kernel_sig)
if ~exist('kernel_sig','var'), kernel_sig = [];end
if isempty(kernel_sig),kernel_sig = .01;end

gaussFunc = @(x,sig) 1./(sig*sqrt(2*pi)).*exp(-x.^2./(2*sig^2));
if ~isa(sample_pntNum,'floQ')
    sample_pnt = CPS_sampling('pnt_num',sample_pntNum);
    sample_pnt = floQ(sample_pnt);
else
    sample_pnt = sample_pntNum;
end
pnt_dist   = real(asin(norm(sample_pnt+conj(sphpoint))./2)).*2;  
weightMat  = gaussFunc(pnt_dist,kernel_sig);
sampleDen  = mean(weightMat,2);
end