Example of use of dcoordinates
Revision as of 18:09, 12 February 2019 by Daniel Castaño (talk | contribs) (Created page with "<pre>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],...")
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)));