SCM Repository
[diderot] / trunk / test / vr-MIP.diderot |
View of /trunk/test/vr-MIP.diderot
Parent Directory
|
Revision Log
Revision 206 -
(download)
(annotate)
Tue Aug 3 16:55:43 2010 UTC (11 years, 10 months ago) by glk
File size: 1269 byte(s)
Tue Aug 3 16:55:43 2010 UTC (11 years, 10 months ago) by glk
File size: 1269 byte(s)
added camera info for vox1 rendering
// vr-MIP.diderot // // maximum intensity projection in Diderot // input string dataFile; // name of dataset input real stepSz = 0.1; // size of steps // e.g. 0.1 input vec3 eye; // location of eye point // e.g. (25,15,10) input vec3 orig; // location of pixel (0,0) // e.g. txs: (8.83877,2.5911,7.65275) // e.g. vox1: (7.45013,2.9202,5.53883) input vec3 cVec; // vector between pixels horizontally // e.g. txs: (-0.0151831,0.0278357,0) // e.g. vox1: (-0.00930055,0.0164548,0) input vec3 rVec; // vector between pixels vertically // e.g. txs: (0.0074887,0.00408474,-0.0305383) // e.g. vox1: (0.00476876,0.00269539,-0.0180902) image(3)[] img = load (dataFile); //field#1(3)[] F = convolve (bspln3, img); field#2(3)[] F = convolve (bspln3, img); actor RayCast (int row, int col) { vec3 pos = orig + real(row)*rVec + real(col)*cVec; vec3 dir = (pos - eye)/|pos - eye|; real t = 0.0; output real maxval = -inf; update { pos = pos + stepSz*dir; if (inside (pos,F)) { real val = F@pos; maxval = max(val, maxval); } if (t > 20.0) stabilize; t = t + stepSz; } /* render: output maxval */ } initially [ RayCast(r, c) | r in 0..199, c in 0..199 ];
root@smlnj-gforge.cs.uchicago.edu | ViewVC Help |
Powered by ViewVC 1.0.0 |