Skip to content
Snippets Groups Projects
Commit a669b7ea authored by Mark Jenkinson's avatar Mark Jenkinson
Browse files

Replaced all origin and signed pixdims with sforms

parent cd6e986e
No related branches found
No related tags found
No related merge requests found
...@@ -99,7 +99,6 @@ namespace MISCMATHS { ...@@ -99,7 +99,6 @@ namespace MISCMATHS {
{ {
Time_Tracer ts(string("Volume::writeAsFloat" + fname).c_str()); Time_Tracer ts(string("Volume::writeAsFloat" + fname).c_str());
float fmin, fmax;
const ColumnVector& outputvol = *this; const ColumnVector& outputvol = *this;
FSLIO* OP = FslOpen(fname.c_str(), "wb"); FSLIO* OP = FslOpen(fname.c_str(), "wb");
...@@ -115,20 +114,11 @@ namespace MISCMATHS { ...@@ -115,20 +114,11 @@ namespace MISCMATHS {
float *qv = new float[sizeVol]; float *qv = new float[sizeVol];
fmin = outputvol(1);
fmax = outputvol(1);
for (int i=1; i<=sizeVol; i++) 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); qv[i-1] = outputvol(i);
} }
FslSetMinMax(OP, (short)fmin, (short)fmax);
FslWriteHeader(OP); FslWriteHeader(OP);
FslWriteVolumes(OP, qv, 1); FslWriteVolumes(OP, qv, 1);
...@@ -140,7 +130,6 @@ namespace MISCMATHS { ...@@ -140,7 +130,6 @@ namespace MISCMATHS {
void Volume::writeAsInt(const string& fname) void Volume::writeAsInt(const string& fname)
{ {
Time_Tracer ts("Volume::writeAsInt"); Time_Tracer ts("Volume::writeAsInt");
int fmin, fmax;
const ColumnVector& outputvol = *this; const ColumnVector& outputvol = *this;
FSLIO* OP = FslOpen(fname.c_str(), "wb"); FSLIO* OP = FslOpen(fname.c_str(), "wb");
...@@ -157,20 +146,12 @@ namespace MISCMATHS { ...@@ -157,20 +146,12 @@ namespace MISCMATHS {
short *qv = new short[sizeVol]; short *qv = new short[sizeVol];
fmin = (int)outputvol(1);
fmax = (int)outputvol(1);
for (int i=1; i<=sizeVol; i++) 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); qv[i-1] = (short)outputvol(i);
} }
FslSetMinMax(OP, (short)fmin, (short)(fmax+0.9999));
FslWriteHeader(OP); FslWriteHeader(OP);
FslWriteVolumes(OP, qv, 1); FslWriteVolumes(OP, qv, 1);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment