Example of use of dcoordinates

vol  = dsphere(30,256);
v = dmodels.vesicle();
v.center = 128.5*[1,1,1];
v.separation = 20;
v.radius = 40;
v.updateCrop;
t = v.grepTable();
myTemplate = dcylinder([4,12],24);
mask      = myTemplate;
for i=1:size(t,1);
   rotTemplate = dynamo_rot(myTemplate,t(i,7:9));
   rotMask     = dynamo_rot(mask,t(i,7:9));  
   vol = dpkvol.embed(rotTemplate,vol,t(i,24:26),rotMask);     
end 

We can quickly check how our " virus" looks like:

disoview(vol); 

Remember that disoview is a very basic renderer of isosurfaces.

vols = dcoordinates(vol,'cart2sph');

figure;
dshow(squeeze(sum(vols(30:32,:,:),1))); 

figure;
dshow(squeeze(sum(vols(35:40,:,:),1)));