matlab - Creating 3D model from cross-sectional images and normal vectors -


i'm trying map series of binary images 3d space based on pixel position q,p , image's location , normal vector. each image (829x829), have x,y,z position of center , dx,dy,dz across transducer (the normal vector out of image, values in mm). know physical diameter of images (30 mm), i'm able interpolate each pixel's position in reference known centerpoint. important note cannot measure rolling of image (rotation center), assume image oriented height axis in xy plane.

my current method create rotation matrix , local x,y,z each pixel of each frame, transform them multiplying matrix this:

for frame=1:numberofframes      %create unit normal vector      vz = [dx,dy,dz] / norm([dx,dy,dz]);      %define x unit vector upwards, y unit vector orthogonal both     vx = cross([0,1,0],vz);     vy = cross(vz,vx);      %create rotation matrix     r = [vx(1) vy(1) vz(1); vx(2) vy(2) vz(2); vx(3) vy(3) vz(3)];        z = zposition(frame);      q=1:height          x = xposition(frame) + (millimeters/pixel*(height/2 - q));          p=1:width              %counter n makes easier multiply r , re-assign values              (rather cycling through q , p again)             n = n + 1;              y = yposition(frame) + (millimeters/pixel*(width/2 - p));              %create position array correlates positions each pixel             positions(n,:) = [x,y,z];         end     end transformedpositions(:,1:3) = round(positions(:,1:3)*r); end

i create 4d volume matrix indices transformedpositions , proper binary value. create model patch(isosurface(1:maxvaluey,1:maxvaluex,1:maxvaluez,volume));

my resulting models rotated seem flat, every image plane parallel , close together, when in reality images taken several cm apart. think rotation not orienting them differently, , pixel interpolation heavily dependent on image orientation(it should independent).

questions:

  • is matrix sensible? using euler's rotation matrix (with roll=0) gives me looking models different matrix.
  • is order of operations correct here, or should translate pixels , image last (i.e. add x,y,z of centerpoint or add in pixel location after transforming - have tried former , doesn't seem alter image)
  • is setting z values same each image proper method?

any advice appreciated! here's sample image, cross-section looking 360* around transducer in center:sample


Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more