Hi, thank you very much Gabri,<br><br>the thing is that I am not able to find this class in the code. I have cheked several times both the code and the description given in <a href="http://www.itk.org/Doxygen/html/classitk_1_1MRIBiasFieldCorrectionFilter.html#itk_1_1MRIBiasFieldCorrectionFiltera60">
http://www.itk.org/Doxygen/html/classitk_1_1MRIBiasFieldCorrectionFilter.html#itk_1_1MRIBiasFieldCorrectionFiltera60</a> and don't see where the p coefficients are given after each iteration.<br><br>Furthermore, I tried to simulate it by stopping the program and starting it again with the p_coeffs given by the last execution (in steps of 100 or 1000 iterations) and I must say that I'm pretty sure there is a bug in the program in this aspect: given coefficients are not taken into account and they're set at zero.
<br><br>The command line used is (for example):<br><br>BiasCorrector.exe --input Image.mhd --output Im_Corr_01.mhd --input-mask Mask_V2.mhd --class-mean 1000 --class-sigma 100 --degree 3 --volume-max-iteration 50 --coefficients -
0.055574 -0.050695 -0.061534 -0.077325 0.016309 -0.031734 -0.032184 -0.039354 -0.012286 -0.018609 -0.021895 -0.007437 -0.016401 0.014835 0.005662 -0.028097 -0.008361 0.025609 -0.020652 -0.043081<br><br>After hundreds of iterations, the result still got stacked at the same correction (as if they were independent corrections not using a previous vector).
<br><br>Thanks,<br><br> Xavier<br><br><br><br><div><span class="gmail_quote">2008/1/14, gabri <<a href="mailto:tartuz@gmail.com">tartuz@gmail.com</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>try to connect an observer to the class you are using. Look some<br>examples in the "Registration" folder , or read the "Registration" topic<br>in the ITK software on the manual.<br><br>best ,<br>
<br>Gabriele<br><br>Xavier Penya ha scritto:<br>> Hi all,<br>><br>> I am working with this ITK feature, and I would like to take out the<br>> so called "p vectors" for each one of the iterations (and not just the
<br>> final one, as usual) in order to control the correction process. I<br>> would like to find out where I have to change the code.<br>><br>> Thank you very much,<br>><br>> Xavier<br>> ------------------------------------------------------------------------
<br>><br>> _______________________________________________<br>> Insight-users mailing list<br>> <a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><br>> <a href="http://www.itk.org/mailman/listinfo/insight-users">
http://www.itk.org/mailman/listinfo/insight-users</a><br>><br><br><br><br></blockquote></div><br>