Lines Matching refs:eigen
120 void insert_fiber(const CImg<T>& fiber, const CImg<te>& eigen, const CImg<tc>& palette, in insert_fiber() argument
125 float x0 = fiber(0,0), y0 = fiber(0,1), z0 = fiber(0,2), fa0 = eigen.linear_atXYZ(x0,y0,z0,12); in insert_fiber()
128 float x1 = fiber(l,0), y1 = fiber(l,1), z1 = fiber(l,2), fa1 = eigen.linear_atXYZ(x1,y1,z1,12); in insert_fiber()
144 CImg<> get_fibertrack(CImg<T>& eigen, in get_fibertrack() argument
148 { T &u = eigen(i,j,k,3), &v = eigen(i,j,k,4), &w = eigen(i,j,k,5); \ in get_fibertrack()
156 pu = eigen(X0,Y0,Z0,3), in get_fibertrack()
157 pv = eigen(X0,Y0,Z0,4), in get_fibertrack()
158 pw = eigen(X0,Y0,Z0,5); in get_fibertrack()
163 if (X<0 || X>eigen.width() - 1 || Y<0 || Y>eigen.height() - 1 || Z<0 || Z>eigen.depth() - 1 || in get_fibertrack()
164 eigen((int)X,(int)Y,(int)Z,12)<FAmin || l>lmax) stopflag = true; in get_fibertrack()
169 … cx = (int)X, px = (cx - 1<0)?0:cx - 1, nx = (cx + 1>=eigen.width())?eigen.width() - 1:cx + 1, in get_fibertrack()
170 … cy = (int)Y, py = (cy - 1<0)?0:cy - 1, ny = (cy + 1>=eigen.height())?eigen.height() - 1:cy + 1, in get_fibertrack()
171 … cz = (int)Z, pz = (cz - 1<0)?0:cz - 1, nz = (cz + 1>=eigen.depth())?eigen.depth() - 1:cz + 1; in get_fibertrack()
172 const T cu = eigen(cx,cy,cz,3), cv = eigen(cx,cy,cz,4), cw = eigen(cx,cy,cz,5); in get_fibertrack()
185 u0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,3), in get_fibertrack()
186 v0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,4), in get_fibertrack()
187 w0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,5), in get_fibertrack()
188 u1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,3), in get_fibertrack()
189 v1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,4), in get_fibertrack()
190 w1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,5), in get_fibertrack()
191 u2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,3), in get_fibertrack()
192 v2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,4), in get_fibertrack()
193 w2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,5), in get_fibertrack()
194 u3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,3), in get_fibertrack()
195 v3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,4), in get_fibertrack()
196 w3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,5); in get_fibertrack()
214 pu = eigen(X0,Y0,Z0,3); in get_fibertrack()
215 pv = eigen(X0,Y0,Z0,4); in get_fibertrack()
216 pw = eigen(X0,Y0,Z0,5); in get_fibertrack()
221 if (X<0 || X>eigen.width() - 1 || Y<0 || Y>eigen.height() - 1 || Z<0 || Z>eigen.depth() - 1 || in get_fibertrack()
222 eigen((int)X,(int)Y,(int)Z,12)<FAmin || l>lmax) stopflag = true; in get_fibertrack()
226 … cx = (int)X, px = (cx - 1<0)?0:cx - 1, nx = (cx + 1>=eigen.width())?eigen.width() - 1:cx + 1, in get_fibertrack()
227 … cy = (int)Y, py = (cy - 1<0)?0:cy - 1, ny = (cy + 1>=eigen.height())?eigen.height() - 1:cy + 1, in get_fibertrack()
228 … cz = (int)Z, pz = (cz - 1<0)?0:cz - 1, nz = (cz + 1>=eigen.depth())?eigen.depth() - 1:cz + 1; in get_fibertrack()
229 const T cu = eigen(cx,cy,cz,3), cv = eigen(cx,cy,cz,4), cw = eigen(cx,cy,cz,5); in get_fibertrack()
242 u0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,3), in get_fibertrack()
243 v0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,4), in get_fibertrack()
244 w0 = 0.5f*dl*eigen.linear_atXYZ(X,Y,Z,5), in get_fibertrack()
245 u1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,3), in get_fibertrack()
246 v1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,4), in get_fibertrack()
247 w1 = 0.5f*dl*eigen.linear_atXYZ(X + u0,Y + v0,Z + w0,5), in get_fibertrack()
248 u2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,3), in get_fibertrack()
249 v2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,4), in get_fibertrack()
250 w2 = 0.5f*dl*eigen.linear_atXYZ(X + u1,Y + v1,Z + w1,5), in get_fibertrack()
251 u3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,3), in get_fibertrack()
252 v3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,4), in get_fibertrack()
253 w3 = 0.5f*dl*eigen.linear_atXYZ(X + u2,Y + v2,Z + w2,5); in get_fibertrack()
340 CImg<> eigen(tensors.width(),tensors.height(),tensors.depth(),13); in main() local
345 eigen(x,y,z,0) = val[0]; eigen(x,y,z,1) = val[1]; eigen(x,y,z,2) = val[2]; in main()
350 eigen(x,y,z,3) = vec(0,0); eigen(x,y,z,4) = vec(0,1); eigen(x,y,z,5) = vec(0,2); in main()
351 eigen(x,y,z,6) = vec(1,0); eigen(x,y,z,7) = vec(1,1); eigen(x,y,z,8) = vec(1,2); in main()
352 eigen(x,y,z,9) = vec(2,0); eigen(x,y,z,10) = vec(2,1); eigen(x,y,z,11) = vec(2,2); in main()
354 eigen(x,y,z,12) = fa; in main()
364 …std::fprintf(stderr,"\n- Maximum diffusivity = %g, Maximum FA = %g",eigmax,eigen.get_shared_channe… in main()
367 eigen.get_shared_channels(0,2)/=eigmax; in main()
376 XYZ[0] = eigen.width()/2; XYZ[1] = eigen.height()/2; XYZ[2] = eigen.depth()/2; in main()
392 …else { xm = ym = zm = 0; xM = eigen.width() - 1; yM = eigen.height() - 1; zM = eigen.height() - 1;… in main()
393 const CImg<> img = eigen.get_crop(xm,ym,zm,xM,yM,zM); in main()
411 const CImg<> veigen = eigen.get_crop(xm,ym,zm,xM,yM,zM); in main()
414 if (fiber.width()) insert_fiber(fiber,eigen,fiber_palette, in main()
470 …else { xm = ym = zm = 0; xM = eigen.width() - 1; yM = eigen.height() - 1; zM = eigen.height() - 1;… in main()
471 const CImg<> img = eigen.get_crop(xm,ym,zm,xM,yM,zM); in main()
473 eigen.get_shared_channel(0).mean(), in main()
474 eigen.get_shared_channel(12).mean()); in main()
486 const CImg<> fiber = get_fibertrack(eigen,x,y,z,lmax,dl,famin,cmin); in main()
508 const CImg<> fiber = get_fibertrack(eigen,x,y,z,lmax,dl,famin,cmin); in main()
525 const CImg<> fiber = get_fibertrack(eigen,x,y,z,lmax,dl,famin,cmin); in main()