Mesh creation in nirfast

I am doing 3d reconstruction in which I have used cylindrical mesh with 12695 node points are there. Input mua=0.02, but the reconstructed mua=0.0159. I want to create a more dense cylinder mesh. How do I use these files?

function mesh = create_mesh(outputfn,shape,sizevar,type)

% create_mesh(outputfn,shape,sizevar,type,link,source,meas)
%
% Creates a mesh with the desired properties
%
% outputfn is the filename to save the mesh to
% shape is the shape of the mesh
% sizevar is a structured variable containing the shape sizes
% type is the mesh type


%% error checking
if (strcmp(shape,'circle') || strcmp(shape,'rectangle')) && ...
        (strcmp(type,'stnd_bem') || strcmp(type,'fluor_bem') || strcmp(type,'spec_bem'))
    errordlg('BEM works for 3D shapes only','NIRFAST Error');
    error('BEM works for 3D shapes only');
end

sizevar.outputfn = outputfn;
%% node, elements, bndvtx
if strcmp(type,'stnd_bem') || strcmp(type,'fluor_bem') || strcmp(type,'spec_bem')
    eval(['mesh = make_' lower(shape) '(sizevar,1);']);
else
    eval(['mesh = make_' lower(shape) '(sizevar);']);
end

%% type, name
mesh.type = type;
mesh.name = outputfn;

%% dimension
if strcmp(shape,'Circle') || strcmp(shape,'Rectangle')
    mesh.dimension = 2;
else
    mesh.dimension = 3;
end

%% region
if strcmp(type,'stnd_bem') || strcmp(type,'fluor_bem') || strcmp(type,'spec_bem')
    mesh.region = [ones(size(mesh.elements,1),1) zeros(size(mesh.elements,1),1)];
else
    mesh.region = zeros(size(mesh.nodes,1),1);
end

%% param
mesh = set_mesh_type(mesh,type);

%% save the mesh
save_mesh(mesh,outputfn);

____________________________________________________________________
function mesh = make_cylinder(sizevar,surf)

% mesh = make_cylinder(sizevar)
%
% Creates a 3D cylindrical mesh
%
% sizevar.xc - x coordinate of the center
% sizevar.yc - y coordinate of the center
% sizevar.zc - z coordinate of the center
% sizevar.r - radius
% sizevar.height - height
% sizevar.dist - distance between nodes
% surf - optional, if 1 surface mesh is created


% error checking
if ~isfield(sizevar,'xc') || ~isfield(sizevar,'yc') || ~isfield(sizevar,'height') ...
       || ~isfield(sizevar,'zc') || ~isfield(sizevar,'r') || ~isfield(sizevar,'dist')
    errordlg('Missing input(s) for cylinder creation','NIRFAST Error');
    error('Missing input(s) for cylinder creation');
end
if sizevar.r <= 0
    errordlg('Radius must be positive','NIRFAST Error');
    error('Radius must be positive');
end
if sizevar.dist <= 0
    errordlg('Distance between nodes must be positive','NIRFAST Error');
    error('Distance between nodes must be positive');
end
if sizevar.height <= 0
    errordlg('Height must be positive','NIRFAST Error');
    error('Height must be positive');
end

hf = waitbar(0,'Creating boundary nodes');

x=[];
y=[];
z=[];
x(1) = 0;
y(1) = 0;
z(1) = 0;
i=2;

delta = sizevar.dist;
h = sizevar.height;
ri = sizevar.r;

hi = 0-h/2:delta:h/2;
hmax = hi(end); hmin = -h/2;
hi = hi(2:end-1);
thi=0:2*3.1415/floor(2*3.1415*ri/sizevar.dist):2*3.141;
[theta, R, Z] = ndgrid(thi, ri, hi);
[x y z] = pol2cart(theta, R, Z);
x = x(:); y = y(:); z = z(:);

clear theta;
i = 1;
for ri=sizevar.dist:sizevar.dist:sizevar.r
    for thi=0:2*3.1415/floor(2*3.1415*ri/sizevar.dist):2*3.141
        r(i) = ri; theta(i) = thi;
        i = i + 1;
    end
end
[t1 t2 t3] = pol2cart(theta, r, ones(1,length(r))*hmax);
x = [x; t1']; y = [y; t2']; z = [z; t3'];
[t1 t2 t3] = pol2cart(theta, r, ones(1,length(r))*hmin);
x = [x; t1']; y = [y; t2']; z = [z; t3'];

mesh.bndvtx = ones(length(x),1);
x = x + sizevar.xc;
y = y + sizevar.yc;
z = z + sizevar.zc;
mesh.nodes(:,1) = x(:);
mesh.nodes(:,2) = y(:);
mesh.nodes(:,3) = z(:);
    
if isfield(sizevar,'outputfn') && ~isempty(sizevar.outputfn)
    outputdir = fileparts(sizevar.outputfn);
else
    outputdir = pwd;
end

waitbar(0.1,hf,'Creating surface');
    
if exist('surf','var') && surf == 1
    mesh.elements = MyRobustCrust(mesh.nodes);
    mesh.elements = FixPatchOrientation(mesh.nodes,mesh.elements,[],1);
else
    mesh.elements = MyRobustCrust(mesh.nodes);
    writenodelm_nod_elm(fullfile(outputdir,'test_node_ele'),...
        mesh.elements,mesh.nodes);
    waitbar(0.6,hf,'Creating volume');
    mesh = checkerboard3d_mm(fullfile(outputdir,...
        'test_node_ele.ele'),'stnd');
    delete(fullfile(outputdir,'test_node_ele.node'),...
        fullfile(outputdir,'test_node_ele.ele'));
end

waitbar(1.0,hf,'Done');
close(hf);

Answer 1

[darova]

Looks like element size is coded here