2nd Degree polynomial fit for the 3D array
33 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
jupiter
am 22 Sep. 2016
Bearbeitet: Jason Stockmann
am 22 Jul. 2020
I am trying to find the 2nd degree polynomial fit for the 3d array which contains the magnetic field distortion information of water in MR imaging. I have two 3D arrays, one having fieldmap values and the other having magnetic field distortion around the MR sample. I am using the expression (Bfieldmap-Xi*Bsample), where Xi is a random value for susceptibility to find the data and try to fit this. But I am not quite sure how to find the fit for the 3D array, since it will have 10 coordinates including all three directions. Please help if someone has already an info on this.
Thanks,
Guru
3 Kommentare
Jason Stockmann
am 22 Jul. 2020
Bearbeitet: Jason Stockmann
am 22 Jul. 2020
Guru, you could try the code below. This worked for me. 'data' is your distortion map. 'polyorder' is the scalar input specifying the order of polynomial you'd like to fit to the data. You need to vectorize both your MRI distortion map (dependent variable) and the coordinate system (independent variable). I am also fitting a 3D polynomial to MRI field map data. I just picked integer indices for the coordinate system (independent variables). I didn't bother scaling them into meaningful values for the image field of view. You could replace them with physically meaningful values if you intend to use them for plotting, etc. later on.
dims=size(data);
[XX,YY,ZZ] = ndgrid(1:dims(1),1:dims(2),1:dims(3));
polymodel = polyfitn([XX(:) YY(:) ZZ(:)],data(:),polyorder);
ypred = polyvaln(polymodel,[XX(:) YY(:) ZZ(:)]);
ypred_array = reshape(ypred,dims);
Akzeptierte Antwort
Steven Lord
am 23 Sep. 2016
If you have Curve Fitting Toolbox it can perform surface fitting. Use the Curve Fitting App to perform the fitting interactively or see "Fit and Plot a Polynomial Surface" on this documentation page for instructions on how to do it programmatically.
If you don't have Curve Fitting Toolbox, see the "Multiple Regression" section on this documentation page.
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Linear and Nonlinear Regression finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!