diff --git a/volume.cc b/volume.cc index 97d8917bc507b2a2b5ce19a4df22689e9dddd8a6..11a2f8838d55b9a1460b44895bf47f27e0bf20bd 100644 --- a/volume.cc +++ b/volume.cc @@ -99,7 +99,6 @@ namespace MISCMATHS { { Time_Tracer ts(string("Volume::writeAsFloat" + fname).c_str()); - float fmin, fmax; const ColumnVector& outputvol = *this; FSLIO* OP = FslOpen(fname.c_str(), "wb"); @@ -115,20 +114,11 @@ namespace MISCMATHS { float *qv = new float[sizeVol]; - fmin = outputvol(1); - fmax = outputvol(1); for (int i=1; i<=sizeVol; i++) { - if(outputvol(i) < fmin) - fmin = outputvol(i); - else if(outputvol(i) > fmax) - fmax = outputvol(i); - qv[i-1] = outputvol(i); } - FslSetMinMax(OP, (short)fmin, (short)fmax); - FslWriteHeader(OP); FslWriteVolumes(OP, qv, 1); @@ -140,7 +130,6 @@ namespace MISCMATHS { void Volume::writeAsInt(const string& fname) { Time_Tracer ts("Volume::writeAsInt"); - int fmin, fmax; const ColumnVector& outputvol = *this; FSLIO* OP = FslOpen(fname.c_str(), "wb"); @@ -157,20 +146,12 @@ namespace MISCMATHS { short *qv = new short[sizeVol]; - fmin = (int)outputvol(1); - fmax = (int)outputvol(1); + for (int i=1; i<=sizeVol; i++) { - if(outputvol(i) < fmin) - fmin = (int)outputvol(i); - else if(outputvol(i) > fmax) - fmax = (int)outputvol(i); - qv[i-1] = (short)outputvol(i); } - FslSetMinMax(OP, (short)fmin, (short)(fmax+0.9999)); - FslWriteHeader(OP); FslWriteVolumes(OP, qv, 1);