From 98f6a46533b5961de6fdc99327975bf790d5b82c Mon Sep 17 00:00:00 2001
From: Tim Behrens <behrens@fmrib.ox.ac.uk>
Date: Wed, 7 Jul 2004 10:37:07 +0000
Subject: [PATCH] *** empty log message ***
---
particle.h | 7 ++++++-
probtrack.cc | 5 +++--
2 files changed, 9 insertions(+), 3 deletions(-)
diff --git a/particle.h b/particle.h
index cf1a7c8..b5220fd 100644
--- a/particle.h
+++ b/particle.h
@@ -178,7 +178,12 @@ namespace PARTICLE{
}
}
-
+
+ void set_dir(const float& rx,const float& ry,const float& rz){
+ m_rx=rx;m_ry=ry;m_rz=rz;m_has_jumped=true;
+ }
+
+
bool check_dir(const float& theta,const float& phi, const float& thr){
if(m_has_jumped){
float rx_new=cos(phi)*sin(theta);
diff --git a/probtrack.cc b/probtrack.cc
index e69de89..20301ae 100644
--- a/probtrack.cc
+++ b/probtrack.cc
@@ -1923,7 +1923,7 @@ void meshtrack(){
x=xst;y=yst;z=zst;
part.change_xyz(x,y,z);
-
+ part.set_dir((*i)->local_normal().X,(*i)->local_normal().Y,(*i)->local_normal().Z);//Set the start dir so that we track inwards from cortex
for( int it = 1 ; it <= nsteps/2; it++){
if( (mask( round(part.x()), round(part.y()), round(part.z())) == 1) ){
if(opts.loopcheck.value()){
@@ -1989,7 +1989,8 @@ void meshtrack(){
tmp2=rand(); tmp2/=RAND_MAX;
if(th_ph_f(3)>tmp2){
- if(!part.check_dir(th_ph_f(1),th_ph_f(2),opts.c_thr.value())){
+ if(!part.check_dir(th_ph_f(1),th_ph_f(2),opts.c_thr.value()) && it!=1){
+ //Don't do curvature checking on the first step as we have set the old direction to the surface normal
break;
}
--
GitLab