<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<title>img_ana.c Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="doxygen.css" rel="stylesheet" type="text/css"/>
</head>
<body>
<!-- Generated by Doxygen 1.7.4 -->
<div id="top">
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li><a href="dirs.html"><span>Directories</span></a></li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>Globals</span></a></li>
</ul>
</div>
<div class="header">
<div class="headertitle">
<div class="title">img_ana.c</div> </div>
</div>
<div class="contents">
<a href="img__ana_8c.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/******************************************************************************</span>
<a name="l00002"></a>00002 <span class="comment"></span>
<a name="l00003"></a>00003 <span class="comment"> Copyright (c) 2007,2008,2009 Turku PET Centre</span>
<a name="l00004"></a>00004 <span class="comment"></span>
<a name="l00005"></a>00005 <span class="comment"> Library: img_ana.c</span>
<a name="l00006"></a>00006 <span class="comment"> Description: I/O routines for IMG data from/to Analyze 7.5 format.</span>
<a name="l00007"></a>00007 <span class="comment"></span>
<a name="l00008"></a>00008 <span class="comment"> This library is free software; you can redistribute it and/or</span>
<a name="l00009"></a>00009 <span class="comment"> modify it under the terms of the GNU Lesser General Public</span>
<a name="l00010"></a>00010 <span class="comment"> License as published by the Free Software Foundation; either</span>
<a name="l00011"></a>00011 <span class="comment"> version 2.1 of the License, or (at your option) any later version.</span>
<a name="l00012"></a>00012 <span class="comment"></span>
<a name="l00013"></a>00013 <span class="comment"> This library is distributed in the hope that it will be useful,</span>
<a name="l00014"></a>00014 <span class="comment"> but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
<a name="l00015"></a>00015 <span class="comment"> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.</span>
<a name="l00016"></a>00016 <span class="comment"> See the GNU Lesser General Public License for more details:</span>
<a name="l00017"></a>00017 <span class="comment"> http://www.gnu.org/copyleft/lesser.html</span>
<a name="l00018"></a>00018 <span class="comment"></span>
<a name="l00019"></a>00019 <span class="comment"> You should have received a copy of the GNU Lesser General Public License</span>
<a name="l00020"></a>00020 <span class="comment"> along with this library/program; if not, write to the Free Software</span>
<a name="l00021"></a>00021 <span class="comment"> Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA </span>
<a name="l00022"></a>00022 <span class="comment"></span>
<a name="l00023"></a>00023 <span class="comment"> Turku PET Centre, Turku, Finland, http://www.turkupetcentre.fi/</span>
<a name="l00024"></a>00024 <span class="comment"></span>
<a name="l00025"></a>00025 <span class="comment"> Modification history:</span>
<a name="l00026"></a>00026 <span class="comment"> 2007-01-31 Vesa Oikonen</span>
<a name="l00027"></a>00027 <span class="comment"> Functions moved from imgfile.c.</span>
<a name="l00028"></a>00028 <span class="comment"> Prompts and randoms are read from SIF data to IMG when reading Analyze.</span>
<a name="l00029"></a>00029 <span class="comment"> 2007-02-27 VO</span>
<a name="l00030"></a>00030 <span class="comment"> Bug correction in imgWriteAnalyze(): fp was not closed in all errors.</span>
<a name="l00031"></a>00031 <span class="comment"> 2007-03-25 VO</span>
<a name="l00032"></a>00032 <span class="comment"> Added functions imgReadAnalyzeHeader(), imgGetAnalyzeHeader(),</span>
<a name="l00033"></a>00033 <span class="comment"> imgSetAnalyzeHeader(), imgReadAnalyzeFrame(), imgReadAnalyzeFirstFrame(),</span>
<a name="l00034"></a>00034 <span class="comment"> and imgWriteAnalyzeFrame().</span>
<a name="l00035"></a>00035 <span class="comment"> 2007-17-07 Harri Merisaari</span>
<a name="l00036"></a>00036 <span class="comment"> Modified for optional ANSi compatibility </span>
<a name="l00037"></a>00037 <span class="comment"> 2007-09-10 VO</span>
<a name="l00038"></a>00038 <span class="comment"> Return value of localtime() is checked.</span>
<a name="l00039"></a>00039 <span class="comment"> 2008-07-07 VO</span>
<a name="l00040"></a>00040 <span class="comment"> If information on decay correction is not found in Analyze header, then</span>
<a name="l00041"></a>00041 <span class="comment"> it is assumed that image data is corrected for decay; previously assumed</span>
<a name="l00042"></a>00042 <span class="comment"> that image was NOT corrected for decay.</span>
<a name="l00043"></a>00043 <span class="comment"> 2009-12-10 VO</span>
<a name="l00044"></a>00044 <span class="comment"> strcpy() replaced with strncpy() in filling header fields.</span>
<a name="l00045"></a>00045 <span class="comment"></span>
<a name="l00046"></a>00046 <span class="comment"></span>
<a name="l00047"></a>00047 <span class="comment">******************************************************************************/</span>
<a name="l00048"></a>00048 <span class="preprocessor">#include <stdio.h></span>
<a name="l00049"></a>00049 <span class="preprocessor">#include <stdlib.h></span>
<a name="l00050"></a>00050 <span class="preprocessor">#include <unistd.h></span>
<a name="l00051"></a>00051 <span class="preprocessor">#include <math.h></span>
<a name="l00052"></a>00052 <span class="preprocessor">#include <string.h></span>
<a name="l00053"></a>00053 <span class="preprocessor">#include <time.h></span>
<a name="l00054"></a>00054 <span class="comment">/*****************************************************************************/</span>
<a name="l00055"></a>00055 <span class="preprocessor">#include "petc99.h"</span>
<a name="l00056"></a>00056 <span class="preprocessor">#include "swap.h"</span>
<a name="l00057"></a>00057 <span class="preprocessor">#include "halflife.h"</span>
<a name="l00058"></a>00058 <span class="comment">/*****************************************************************************/</span>
<a name="l00059"></a>00059 <span class="preprocessor">#include "<a class="code" href="img_8h.html">include/img.h</a>"</span>
<a name="l00060"></a>00060 <span class="preprocessor">#include "<a class="code" href="analyze_8h.html">include/analyze.h</a>"</span>
<a name="l00061"></a>00061 <span class="preprocessor">#include "<a class="code" href="imgmax_8h.html">include/imgmax.h</a>"</span>
<a name="l00062"></a>00062 <span class="preprocessor">#include "<a class="code" href="imgdecay_8h.html">include/imgdecay.h</a>"</span>
<a name="l00063"></a>00063 <span class="preprocessor">#include "<a class="code" href="sif_8h.html">include/sif.h</a>"</span>
<a name="l00064"></a>00064 <span class="preprocessor">#include "<a class="code" href="imgfile_8h.html">include/imgfile.h</a>"</span>
<a name="l00065"></a>00065 <span class="comment">/*****************************************************************************/</span>
<a name="l00066"></a>00066
<a name="l00067"></a>00067 <span class="comment">/*****************************************************************************/</span>
<a name="l00083"></a><a class="code" href="imgfile_8h.html#aefc54f0812e16d43a17edcdc342be35b">00083</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#aefc54f0812e16d43a17edcdc342be35b">imgReadAnalyze</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *dbname, <a class="code" href="structIMG.html">IMG</a> *img) {
<a name="l00084"></a>00084 FILE *fp;
<a name="l00085"></a>00085 <span class="keywordtype">int</span> ret, fi, pi, xi, yi;
<a name="l00086"></a>00086 <span class="keywordtype">float</span> *fdata=NULL, *fptr;
<a name="l00087"></a>00087 <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> dsr;
<a name="l00088"></a>00088 <span class="keywordtype">char</span> datfile[FILENAME_MAX], hdrfile[FILENAME_MAX], siffile[FILENAME_MAX];
<a name="l00089"></a>00089 <span class="keywordtype">char</span> buf[128], *cptr;
<a name="l00090"></a>00090 <span class="keywordtype">int</span> dimNr, dimx, dimy, dimz=1, dimt=1, pxlNr=0;
<a name="l00091"></a>00091 <a class="code" href="structSIF.html">SIF</a> sif;
<a name="l00092"></a>00092 <span class="keyword">struct </span>tm *st;
<a name="l00093"></a>00093
<a name="l00094"></a>00094
<a name="l00095"></a>00095 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"imgReadAnalyze(%s, *img)\n"</span>, dbname);
<a name="l00096"></a>00096
<a name="l00097"></a>00097 <span class="comment">/* Check the arguments */</span>
<a name="l00098"></a>00098 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00099"></a>00099 <span class="keywordflow">if</span>(img==NULL || img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a7aeab2fb426b858c78ec515aaa8f0a67">IMG_STATUS_INITIALIZED</a>) {
<a name="l00100"></a>00100 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>); <span class="keywordflow">return</span>(2);}
<a name="l00101"></a>00101 <span class="keywordflow">if</span>(dbname==NULL || !dbname[0]) {<a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>); <span class="keywordflow">return</span>(1);}
<a name="l00102"></a>00102
<a name="l00103"></a>00103 <span class="comment">/* Make the image and header filenames */</span>
<a name="l00104"></a>00104 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a962221ab7adb0bed2050284912d473ca">anaExists</a>(dbname)==0) {
<a name="l00105"></a>00105 <span class="comment">/* Check if filename was given accidentally with extension */</span>
<a name="l00106"></a>00106 strcpy(datfile, dbname); cptr=strrchr(datfile, <span class="charliteral">'.'</span>);
<a name="l00107"></a>00107 <span class="keywordflow">if</span>(cptr!=NULL && (strcmp(cptr, <span class="stringliteral">".img"</span>)==0 || strcmp(cptr, <span class="stringliteral">".hdr"</span>)==0)) {
<a name="l00108"></a>00108 *cptr=(char)0; strcpy(hdrfile, datfile);
<a name="l00109"></a>00109 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a962221ab7adb0bed2050284912d473ca">anaExists</a>(datfile)==0) { <span class="comment">/* still not found */</span>
<a name="l00110"></a>00110 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba201ea0defe535c9b0449047fd6337370">STATUS_NOHEADERFILE</a>); <span class="keywordflow">return</span>(3);}
<a name="l00111"></a>00111 strcat(datfile, <span class="stringliteral">".img"</span>); strcat(hdrfile, <span class="stringliteral">".hdr"</span>);
<a name="l00112"></a>00112 } <span class="keywordflow">else</span> {
<a name="l00113"></a>00113 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba201ea0defe535c9b0449047fd6337370">STATUS_NOHEADERFILE</a>); <span class="keywordflow">return</span>(3);
<a name="l00114"></a>00114 }
<a name="l00115"></a>00115 } <span class="keywordflow">else</span> {
<a name="l00116"></a>00116 <span class="comment">/* Database name was given and img and hdr files were found */</span>
<a name="l00117"></a>00117 strcpy(datfile, dbname); strcat(datfile, <span class="stringliteral">".img"</span>);
<a name="l00118"></a>00118 strcpy(hdrfile, dbname); strcat(hdrfile, <span class="stringliteral">".hdr"</span>);
<a name="l00119"></a>00119 }
<a name="l00120"></a>00120
<a name="l00121"></a>00121 <span class="comment">/* Read Analyze header file */</span>
<a name="l00122"></a>00122 ret=<a class="code" href="analyze_8c.html#aa7bf1965dcdbff0686193675a3175a1c">anaReadHeader</a>(hdrfile, &dsr);
<a name="l00123"></a>00123 <span class="keywordflow">if</span>(ret) {
<a name="l00124"></a>00124 <span class="keywordflow">if</span>(ret==1) <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00125"></a>00125 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(ret==2) <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba201ea0defe535c9b0449047fd6337370">STATUS_NOHEADERFILE</a>);
<a name="l00126"></a>00126 <span class="keywordflow">else</span> <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baef2b214bb9b6ccc009f27c76e81b3fef">STATUS_UNSUPPORTED</a>);
<a name="l00127"></a>00127 <span class="keywordflow">return</span>(3);
<a name="l00128"></a>00128 }
<a name="l00129"></a>00129 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) <a class="code" href="analyze_8c.html#a5c764c78d790dabed25b1f7971717a05">anaPrintHeader</a>(&dsr, stdout);
<a name="l00130"></a>00130
<a name="l00131"></a>00131 <span class="comment">/* Open image datafile */</span>
<a name="l00132"></a>00132 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) fprintf(stdout, <span class="stringliteral">"reading image data %s\n"</span>, datfile);
<a name="l00133"></a>00133 <span class="keywordflow">if</span>((fp=fopen(datfile, <span class="stringliteral">"rb"</span>)) == NULL) {<a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba524fb2fd0efecb30fe42b109b3d28fe8">STATUS_NOIMGDATA</a>); <span class="keywordflow">return</span>(5);}
<a name="l00134"></a>00134
<a name="l00135"></a>00135 <span class="comment">/* Prepare IMG for Analyze image */</span>
<a name="l00136"></a>00136 <span class="comment">/* Get the image dimensions from header */</span>
<a name="l00137"></a>00137 dimNr=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[0];
<a name="l00138"></a>00138 <span class="keywordflow">if</span>(dimNr<2) {fclose(fp); <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7d0a6ddf0a60e475f14f6b0e9b6487a1">STATUS_INVALIDHEADER</a>); <span class="keywordflow">return</span>(4);}
<a name="l00139"></a>00139 dimx=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[1]; dimy=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[2];
<a name="l00140"></a>00140 <span class="keywordflow">if</span>(dimNr>2) {dimz=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[3]; <span class="keywordflow">if</span>(dimNr>3) dimt=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4];}
<a name="l00141"></a>00141 pxlNr=dimx*dimy*dimz;
<a name="l00142"></a>00142 <span class="keywordflow">if</span>(pxlNr<1) {fclose(fp); <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7d0a6ddf0a60e475f14f6b0e9b6487a1">STATUS_INVALIDHEADER</a>); <span class="keywordflow">return</span>(4);}
<a name="l00143"></a>00143 <span class="comment">/* Allocate memory for IMG */</span>
<a name="l00144"></a>00144 ret=<a class="code" href="img_8c.html#a9b2670da9c7685e6d53058cdf1db1199">imgAllocate</a>(img, dimz, dimy, dimx, dimt);
<a name="l00145"></a>00145 <span class="keywordflow">if</span>(ret) {fclose(fp); <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>); <span class="keywordflow">return</span>(11);}
<a name="l00146"></a>00146 <span class="comment">/* Copy information from Analyze header */</span>
<a name="l00147"></a>00147 img-><a class="code" href="structIMG.html#a0d9199e1c3c4bd94d56eb30e421b309c">type</a>=<a class="code" href="img_8h.html#ae8b5bb16b84adcbcd2bba3e10a36d62d">IMG_TYPE_IMAGE</a>;
<a name="l00148"></a>00148 strncpy(img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>, dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ab9d251da82c987e8fe2eb23c6cf83f57">patient_id</a>, MAX_STUDYNR_LEN);
<a name="l00149"></a>00149 strcpy(img-><a class="code" href="structIMG.html#a29013e391052918bd5aa6974eeec811d">patientName</a>, dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ab9d251da82c987e8fe2eb23c6cf83f57">patient_id</a>);
<a name="l00150"></a>00150 img-><a class="code" href="structIMG.html#a02cc239e917e7ebf85e34f21448ff080">sizex</a>=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[1];
<a name="l00151"></a>00151 img-><a class="code" href="structIMG.html#a80b0609c63f33934bf6e6678c6121cb3">sizey</a>=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[2];
<a name="l00152"></a>00152 img-><a class="code" href="structIMG.html#afc2a1d2417a8a3afe2c500173a05c472">sizez</a>=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[3];
<a name="l00153"></a>00153 <span class="comment">/*if(dsr.dime.funused2>1.E-5) img->zoom=dsr.dime.funused2;*/</span>
<a name="l00154"></a>00154 <span class="keywordflow">if</span>(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>>1.E-5) img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>;
<a name="l00155"></a>00155 <span class="keywordflow">for</span>(pi=0; pi<dimz; pi++) img-><a class="code" href="structIMG.html#a78e3164d578cf4a570e930dceb10fb51">planeNumber</a>[pi]=pi+1;
<a name="l00156"></a>00156 <span class="keywordflow">if</span>(dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>) img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>=<a class="code" href="img_8h.html#a38e113a0101e5e683653263ddabe0fc8">IMG_ANA_L</a>; <span class="keywordflow">else</span> img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>=<a class="code" href="img_8h.html#a9145a82aa093b7e52a693ec5c38094ef">IMG_ANA</a>;
<a name="l00157"></a>00157 <span class="comment">/* Decay correction */</span>
<a name="l00158"></a>00158 <span class="keywordflow">if</span>(strstr(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"Decay corrected."</span>)!=NULL)
<a name="l00159"></a>00159 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=1;
<a name="l00160"></a>00160 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(strstr(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"No decay correction."</span>)!=NULL)
<a name="l00161"></a>00161 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=0;
<a name="l00162"></a>00162 <span class="keywordflow">else</span>
<a name="l00163"></a>00163 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=1;
<a name="l00164"></a>00164
<a name="l00165"></a>00165 <span class="comment">/* Allocate memory for one image frame */</span>
<a name="l00166"></a>00166 fdata=malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>));
<a name="l00167"></a>00167 <span class="keywordflow">if</span>(fdata==NULL) {fclose(fp); <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>); <span class="keywordflow">return</span>(12);}
<a name="l00168"></a>00168
<a name="l00169"></a>00169 <span class="comment">/* Read one image frame at a time */</span>
<a name="l00170"></a>00170 <span class="keywordflow">for</span>(fi=0; fi<dimt; fi++) {
<a name="l00171"></a>00171 fptr=fdata;
<a name="l00172"></a>00172 ret=<a class="code" href="analyze_8c.html#abfb69549a8e8c6a4eab3c52045262d22">anaReadImagedata</a>(fp, &dsr, fi+1, fptr);
<a name="l00173"></a>00173 <span class="keywordflow">if</span>(ret) {free(fdata); fclose(fp); <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba524fb2fd0efecb30fe42b109b3d28fe8">STATUS_NOIMGDATA</a>); <span class="keywordflow">return</span>(7);}
<a name="l00174"></a>00174 <span class="comment">/* Copy pixel values to IMG */</span>
<a name="l00175"></a>00175 fptr=fdata;
<a name="l00176"></a>00176 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a55febab802306e82d162793558f5fbce">anaFlipping</a>()==0) { <span class="comment">/* no flipping in z-direction */</span>
<a name="l00177"></a>00177 <span class="keywordflow">for</span>(pi=0; pi<img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>; pi++)
<a name="l00178"></a>00178 <span class="keywordflow">for</span>(yi=dimy-1; yi>=0; yi--)
<a name="l00179"></a>00179 <span class="keywordflow">for</span>(xi=dimx-1; xi>=0; xi--)
<a name="l00180"></a>00180 img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][fi]=*fptr++;
<a name="l00181"></a>00181 } <span class="keywordflow">else</span> {
<a name="l00182"></a>00182 <span class="keywordflow">for</span>(pi=dimz-1; pi>=0; pi--)
<a name="l00183"></a>00183 <span class="keywordflow">for</span>(yi=dimy-1; yi>=0; yi--)
<a name="l00184"></a>00184 <span class="keywordflow">for</span>(xi=dimx-1; xi>=0; xi--)
<a name="l00185"></a>00185 img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][fi]=*fptr++;
<a name="l00186"></a>00186 }
<a name="l00187"></a>00187 } <span class="comment">/* next frame */</span>
<a name="l00188"></a>00188 free(fdata);
<a name="l00189"></a>00189 fclose(fp);
<a name="l00190"></a>00190
<a name="l00191"></a>00191 <span class="comment">/* Try to read frame time information from SIF file */</span>
<a name="l00192"></a>00192 <span class="comment">/* Make filename from database or image data file */</span>
<a name="l00193"></a>00193 strcpy(siffile, dbname); strcat(siffile, <span class="stringliteral">".sif"</span>);
<a name="l00194"></a>00194 <span class="keywordflow">if</span>(access(siffile, 0) == -1) {
<a name="l00195"></a>00195 strcpy(siffile, datfile); strcat(siffile, <span class="stringliteral">".sif"</span>);
<a name="l00196"></a>00196 }
<a name="l00197"></a>00197 <span class="comment">/* Check if SIF file is found */</span>
<a name="l00198"></a>00198 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"reading SIF file %s\n"</span>, siffile);
<a name="l00199"></a>00199 <span class="keywordflow">if</span>(access(siffile, 0) == -1) {
<a name="l00200"></a>00200 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">" No SIF file; therefore unknown frame times.\n"</span>);
<a name="l00201"></a>00201 <span class="keywordflow">return</span>(0);
<a name="l00202"></a>00202 }
<a name="l00203"></a>00203 <span class="comment">/* If found, then read it */</span>
<a name="l00204"></a>00204 <a class="code" href="sif_8c.html#a61c0218e380da8ace65d78f82a8f9f69">sifInit</a>(&sif); ret=<a class="code" href="sifio_8c.html#ae6c3f31d7cbe31ca6b74760978fd39c7">sifRead</a>(siffile, &sif);
<a name="l00205"></a>00205 <span class="keywordflow">if</span>(ret) {<a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba289e951f1be6b96d46b55e834db00f78">STATUS_NOSIFDATA</a>); <span class="keywordflow">return</span>(21);}
<a name="l00206"></a>00206 <span class="comment">/* Copy scan start time */</span>
<a name="l00207"></a>00207 <span class="keywordflow">if</span>(sif.<a class="code" href="structSIF.html#a98b5ff616033e728f0edc9bdc7eddce2">scantime</a>>0) {
<a name="l00208"></a>00208 st=localtime(&sif.<a class="code" href="structSIF.html#a98b5ff616033e728f0edc9bdc7eddce2">scantime</a>);
<a name="l00209"></a>00209 <span class="keywordflow">if</span>(st!=NULL) strftime(buf, 128, <span class="stringliteral">"%Y-%m-%d %H:%M:%S"</span>, st);
<a name="l00210"></a>00210 <span class="keywordflow">else</span> strcpy(buf, <span class="stringliteral">"1900-01-01 00:00:00"</span>);
<a name="l00211"></a>00211 img-><a class="code" href="structIMG.html#a5db78f8b423fc42b80bba42e27c7b2c3">scanStart</a>=sif.<a class="code" href="structSIF.html#a98b5ff616033e728f0edc9bdc7eddce2">scantime</a>;
<a name="l00212"></a>00212 }
<a name="l00213"></a>00213 <span class="comment">/* Copy frame times */</span>
<a name="l00214"></a>00214 <span class="keywordflow">if</span>(sif.<a class="code" href="structSIF.html#a441bd4d1d3825ace8f21d02e9f5d650d">frameNr</a>!=img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>) {
<a name="l00215"></a>00215 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba5c5cd324e637bab1f54039528e2790a7">STATUS_WRONGSIFDATA</a>); <a class="code" href="sif_8c.html#ad2378bec7fec69900af6f52737c89784">sifEmpty</a>(&sif); <span class="keywordflow">return</span>(22);}
<a name="l00216"></a>00216 <span class="keywordflow">for</span>(fi=0; fi<sif.<a class="code" href="structSIF.html#a441bd4d1d3825ace8f21d02e9f5d650d">frameNr</a>; fi++) {
<a name="l00217"></a>00217 img-><a class="code" href="structIMG.html#a8ff22fc209b6f44ad14341c8644788d2">start</a>[fi]=sif.<a class="code" href="structSIF.html#a085f3f1a4b2a77eabe13f52c2010a3c0">x1</a>[fi]; img-><a class="code" href="structIMG.html#aa2b629a51ca62bcafad0a66d32e1ca8b">end</a>[fi]=sif.<a class="code" href="structSIF.html#a45001ad49c925e6ae266357b3ab7d8a4">x2</a>[fi];
<a name="l00218"></a>00218 img-><a class="code" href="structIMG.html#af7c69142be89fac081e7f8bc354b14c8">mid</a>[fi]=0.5*(img-><a class="code" href="structIMG.html#a8ff22fc209b6f44ad14341c8644788d2">start</a>[fi]+img-><a class="code" href="structIMG.html#aa2b629a51ca62bcafad0a66d32e1ca8b">end</a>[fi]);
<a name="l00219"></a>00219 }
<a name="l00220"></a>00220 <span class="comment">/* Copy prompts and randoms */</span>
<a name="l00221"></a>00221 <span class="keywordflow">for</span>(fi=0; fi<sif.<a class="code" href="structSIF.html#a441bd4d1d3825ace8f21d02e9f5d650d">frameNr</a>; fi++) {
<a name="l00222"></a>00222 img-><a class="code" href="structIMG.html#a76d63a001b8569d9f80ccb325e064be1">prompts</a>[fi]=sif.<a class="code" href="structSIF.html#a20d6fecbd751d872300c661e19903eed">prompts</a>[fi];
<a name="l00223"></a>00223 img-><a class="code" href="structIMG.html#ae47ebbbdc0c308f473fec801a2df513d">randoms</a>[fi]=sif.<a class="code" href="structSIF.html#af2fa6ee9e91796baa37261bd866b641f">randoms</a>[fi];
<a name="l00224"></a>00224 }
<a name="l00225"></a>00225
<a name="l00226"></a>00226 <span class="comment">/* Set isotopeHalflife, if isotope is found */</span>
<a name="l00227"></a>00227 <span class="keywordflow">if</span>(strlen(sif.<a class="code" href="structSIF.html#aaf18449475434d1cf5ab11b2c36780b1">isotope_name</a>)>1) {
<a name="l00228"></a>00228 img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>=60.0*hlFromIsotope(sif.<a class="code" href="structSIF.html#aaf18449475434d1cf5ab11b2c36780b1">isotope_name</a>);
<a name="l00229"></a>00229 }
<a name="l00230"></a>00230 <a class="code" href="sif_8c.html#ad2378bec7fec69900af6f52737c89784">sifEmpty</a>(&sif);
<a name="l00231"></a>00231
<a name="l00232"></a>00232 <span class="keywordflow">return</span>(0);
<a name="l00233"></a>00233 }
<a name="l00234"></a>00234 <span class="comment">/*****************************************************************************/</span>
<a name="l00235"></a>00235
<a name="l00236"></a>00236 <span class="comment">/*****************************************************************************/</span>
<a name="l00253"></a><a class="code" href="imgfile_8h.html#a380113542ddd5040b383221f1377d86f">00253</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#a380113542ddd5040b383221f1377d86f">imgWriteAnalyze</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *dbname, <a class="code" href="structIMG.html">IMG</a> *img) {
<a name="l00254"></a>00254 FILE *fp;
<a name="l00255"></a>00255 <span class="keywordtype">int</span> ret, fi, pi, xi, yi, little;
<a name="l00256"></a>00256 <span class="keywordtype">float</span> g;
<a name="l00257"></a>00257 <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> dsr;
<a name="l00258"></a>00258 <span class="keywordtype">char</span> datfile[FILENAME_MAX], hdrfile[FILENAME_MAX];
<a name="l00259"></a>00259 <span class="keyword">const</span> <span class="keywordtype">char</span> *cptr;
<a name="l00260"></a>00260 <span class="keywordtype">int</span> pxlNr=0;
<a name="l00261"></a>00261 <span class="keyword">struct </span>tm *st;
<a name="l00262"></a>00262 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sdata, *sptr, smin, smax;
<a name="l00263"></a>00263
<a name="l00264"></a>00264
<a name="l00265"></a>00265 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"imgWriteAnalyze(%s, *img)\n"</span>, dbname);
<a name="l00266"></a>00266
<a name="l00267"></a>00267 <span class="comment">/* Check the arguments */</span>
<a name="l00268"></a>00268 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00269"></a>00269 <span class="keywordflow">if</span>(img==NULL || img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a>) {
<a name="l00270"></a>00270 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>); <span class="keywordflow">return</span>(2);}
<a name="l00271"></a>00271 <span class="keywordflow">if</span>(dbname==NULL || !dbname[0]) {<a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>); <span class="keywordflow">return</span>(1);}
<a name="l00272"></a>00272
<a name="l00273"></a>00273 <span class="comment">/* Make the image and header filenames */</span>
<a name="l00274"></a>00274 strcpy(datfile, dbname); strcat(datfile, <span class="stringliteral">".img"</span>);
<a name="l00275"></a>00275 strcpy(hdrfile, dbname); strcat(hdrfile, <span class="stringliteral">".hdr"</span>);
<a name="l00276"></a>00276
<a name="l00277"></a>00277
<a name="l00278"></a>00278 <span class="comment">/*</span>
<a name="l00279"></a>00279 <span class="comment"> * Fill Analyze header</span>
<a name="l00280"></a>00280 <span class="comment"> */</span>
<a name="l00281"></a>00281 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>==<a class="code" href="img_8h.html#a38e113a0101e5e683653263ddabe0fc8">IMG_ANA_L</a>) dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>=1; <span class="keywordflow">else</span> dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>=0;
<a name="l00282"></a>00282 <span class="comment">/* Header key */</span>
<a name="l00283"></a>00283 memset(&dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__KEY.html">ANALYZE_HEADER_KEY</a>));
<a name="l00284"></a>00284 memset(&dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__IMGDIM.html">ANALYZE_HEADER_IMGDIM</a>));
<a name="l00285"></a>00285 memset(&dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__HISTORY.html">ANALYZE_HEADER_HISTORY</a>));
<a name="l00286"></a>00286 dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a6d7253a777a1f2c86fa14e77a5c47515">sizeof_hdr</a>=348;
<a name="l00287"></a>00287 strcpy(dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a96f4c147959660ddbd2a579e19e24c3b">data_type</a>, <span class="stringliteral">""</span>);
<a name="l00288"></a>00288 cptr=strrchr(dbname, <span class="charliteral">'/'</span>); <span class="keywordflow">if</span>(cptr==NULL) cptr=strrchr(dbname, <span class="charliteral">'\\'</span>);
<a name="l00289"></a>00289 <span class="keywordflow">if</span>(cptr!=NULL) cptr++; <span class="keywordflow">if</span>(cptr==NULL) cptr=dbname;
<a name="l00290"></a>00290 strncpy(dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a86cb24f71952964e9ebed40133da6725">db_name</a>, cptr, 17);
<a name="l00291"></a>00291 dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a09532be9a54a2e27170c8776c0ae1b3a">extents</a>=16384;
<a name="l00292"></a>00292 dsr.<a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a6be3dcc756c59c1ddec21ca6311f2c5a">regular</a>=<span class="charliteral">'r'</span>;
<a name="l00293"></a>00293 <span class="comment">/* Image dimension */</span>
<a name="l00294"></a>00294 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[0]=4;
<a name="l00295"></a>00295 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[1]=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>;
<a name="l00296"></a>00296 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[2]=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>;
<a name="l00297"></a>00297 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[3]=img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>;
<a name="l00298"></a>00298 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]=img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>;
<a name="l00299"></a>00299 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a2a05d2cb24eeac5ac09098974f6b0fc9">datatype</a>=<a class="code" href="analyze_8h.html#a57a617cb3056b972eabcb94cdc289298">ANALYZE_DT_SIGNED_SHORT</a>;
<a name="l00300"></a>00300 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#ad0ca8c2b6ac4f211ff570c186979bbf9">bitpix</a>=16;
<a name="l00301"></a>00301 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[0]=0.0;
<a name="l00302"></a>00302 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[1]=img-><a class="code" href="structIMG.html#a02cc239e917e7ebf85e34f21448ff080">sizex</a>;
<a name="l00303"></a>00303 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[2]=img-><a class="code" href="structIMG.html#a80b0609c63f33934bf6e6678c6121cb3">sizey</a>;
<a name="l00304"></a>00304 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[3]=img-><a class="code" href="structIMG.html#afc2a1d2417a8a3afe2c500173a05c472">sizez</a>;
<a name="l00305"></a>00305 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[4]=0.0;
<a name="l00306"></a>00306 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>=0.0; <span class="comment">/* Scale factor is set later */</span>
<a name="l00307"></a>00307 <span class="comment">/* dsr.dime.funused2=img->zoom; */</span> <span class="comment">/* Reconstruction zoom */</span>
<a name="l00308"></a>00308 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>=img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>;
<a name="l00309"></a>00309 <span class="comment">/* Data history */</span>
<a name="l00310"></a>00310 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>==1) strcpy(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"Decay corrected."</span>);
<a name="l00311"></a>00311 <span class="keywordflow">else</span> strcpy(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"No decay correction."</span>);
<a name="l00312"></a>00312 strncpy(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a0a750bf7ef6504c76d5c933debe56959">scannum</a>, img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>, 10);
<a name="l00313"></a>00313 st=localtime(&img-><a class="code" href="structIMG.html#a5db78f8b423fc42b80bba42e27c7b2c3">scanStart</a>);
<a name="l00314"></a>00314 <span class="keywordflow">if</span>(st!=NULL) {
<a name="l00315"></a>00315 strftime(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a682bbce3bd0941ee16b62069e4cd1f15">exp_date</a>, 10, <span class="stringliteral">"%Y-%m-%d"</span>, st);
<a name="l00316"></a>00316 strftime(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a3938bc6e0b7a1770366b838d1d96c37e">exp_time</a>, 10, <span class="stringliteral">"%H:%M:%S"</span>, st);
<a name="l00317"></a>00317 } <span class="keywordflow">else</span> {
<a name="l00318"></a>00318 strncpy(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a682bbce3bd0941ee16b62069e4cd1f15">exp_date</a>, <span class="stringliteral">"1900-01-01"</span>, 10);
<a name="l00319"></a>00319 strncpy(dsr.<a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a3938bc6e0b7a1770366b838d1d96c37e">exp_time</a>, <span class="stringliteral">"00:00:00"</span>, 10);
<a name="l00320"></a>00320 }
<a name="l00321"></a>00321
<a name="l00322"></a>00322 <span class="comment">/*</span>
<a name="l00323"></a>00323 <span class="comment"> * Scale data to short int range</span>
<a name="l00324"></a>00324 <span class="comment"> * Determine and set scale factor and cal_min & cal_max</span>
<a name="l00325"></a>00325 <span class="comment"> */</span>
<a name="l00326"></a>00326 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"scaling data to short ints\n"</span>);
<a name="l00327"></a>00327 ret=<a class="code" href="imgmax_8c.html#a97c48d1ed44d81b2d96ce1a139b99009">imgMinMax</a>(img, &dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>, &dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>);
<a name="l00328"></a>00328 <span class="keywordflow">if</span>(ret) {<a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>); <span class="keywordflow">return</span>(3);}
<a name="l00329"></a>00329 <span class="keywordflow">if</span>(fabs(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>)>fabs(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>)) g=fabs(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>);
<a name="l00330"></a>00330 <span class="keywordflow">else</span> g=fabs(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>);
<a name="l00331"></a>00331 <span class="keywordflow">if</span>(g<1E-20) g=1.0; <span class="keywordflow">else</span> g=32767./g; dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>=1.0/g;
<a name="l00332"></a>00332 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"min=%g max=%g scale_factor=%g\n"</span>,
<a name="l00333"></a>00333 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>, dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>, dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>);
<a name="l00334"></a>00334
<a name="l00335"></a>00335 <span class="comment">/* Allocate memory for short int array */</span>
<a name="l00336"></a>00336 pxlNr=(img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>)*(img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>)*(img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>);
<a name="l00337"></a>00337 sdata=malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>));
<a name="l00338"></a>00338 <span class="keywordflow">if</span>(sdata==NULL) {
<a name="l00339"></a>00339 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>);
<a name="l00340"></a>00340 <span class="keywordflow">return</span> 12;
<a name="l00341"></a>00341 }
<a name="l00342"></a>00342
<a name="l00343"></a>00343 <span class="comment">/* Open image data file for write */</span>
<a name="l00344"></a>00344 <span class="keywordflow">if</span>((fp=fopen(datfile, <span class="stringliteral">"wb"</span>)) == NULL) {
<a name="l00345"></a>00345 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baa2334ba6708d894ba27671de890f0e0a">STATUS_CANTWRITEIMGFILE</a>);
<a name="l00346"></a>00346 free(sdata);
<a name="l00347"></a>00347 <span class="keywordflow">return</span> 14;
<a name="l00348"></a>00348 }
<a name="l00349"></a>00349
<a name="l00350"></a>00350 <span class="comment">/* Copy and write image matrix data to short int array */</span>
<a name="l00351"></a>00351 <span class="comment">/* Data is written one frame at a time */</span>
<a name="l00352"></a>00352 smin=smax=temp_roundf(g*img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[0][0][0][0]);
<a name="l00353"></a>00353 <span class="keywordflow">for</span>(fi=0; fi<img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>; fi++) {
<a name="l00354"></a>00354 sptr=sdata;
<a name="l00355"></a>00355 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a55febab802306e82d162793558f5fbce">anaFlipping</a>()==0) {
<a name="l00356"></a>00356 <span class="keywordflow">for</span>(pi=0; pi<img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>; pi++)
<a name="l00357"></a>00357 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--)
<a name="l00358"></a>00358 <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--) {
<a name="l00359"></a>00359 *sptr=temp_roundf(g*img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][fi]);
<a name="l00360"></a>00360 <span class="keywordflow">if</span>(*sptr>smax) smax=*sptr; <span class="keywordflow">else</span> <span class="keywordflow">if</span>(*sptr<smin) smin=*sptr;
<a name="l00361"></a>00361 sptr++;
<a name="l00362"></a>00362 }
<a name="l00363"></a>00363 } <span class="keywordflow">else</span> {
<a name="l00364"></a>00364 <span class="keywordflow">for</span>(pi=img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>-1; pi>=0; pi--)
<a name="l00365"></a>00365 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--)
<a name="l00366"></a>00366 <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--) {
<a name="l00367"></a>00367 *sptr=temp_roundf(g*img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][fi]);
<a name="l00368"></a>00368 <span class="keywordflow">if</span>(*sptr>smax) smax=*sptr; <span class="keywordflow">else</span> <span class="keywordflow">if</span>(*sptr<smin) smin=*sptr;
<a name="l00369"></a>00369 sptr++;
<a name="l00370"></a>00370 }
<a name="l00371"></a>00371 }
<a name="l00372"></a>00372 <span class="comment">/* Change byte order if necessary */</span>
<a name="l00373"></a>00373 little=little_endian();
<a name="l00374"></a>00374 <span class="keywordflow">if</span>(little!=dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>)
<a name="l00375"></a>00375 swabip(sdata, pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>));
<a name="l00376"></a>00376 <span class="comment">/* Write image data */</span>
<a name="l00377"></a>00377 <span class="keywordflow">if</span>(fwrite(sdata, 2, pxlNr, fp) != pxlNr) {
<a name="l00378"></a>00378 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baa2334ba6708d894ba27671de890f0e0a">STATUS_CANTWRITEIMGFILE</a>);
<a name="l00379"></a>00379 free(sdata); fclose(fp);
<a name="l00380"></a>00380 <span class="keywordflow">return</span> 15;
<a name="l00381"></a>00381 }
<a name="l00382"></a>00382 }
<a name="l00383"></a>00383 <span class="comment">/* Done writing */</span>
<a name="l00384"></a>00384 fclose(fp);
<a name="l00385"></a>00385 free(sdata);
<a name="l00386"></a>00386
<a name="l00387"></a>00387 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"smin=%d smax=%d\n"</span>, smin, smax);
<a name="l00388"></a>00388
<a name="l00389"></a>00389 <span class="comment">/* Set header glmin & glmax */</span>
<a name="l00390"></a>00390 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aa92cd82aae233799feeea86dfdd0beb4">glmin</a>=smin; dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a57093cdaa14482a67fcc4a29ae7c0cb9">glmax</a>=smax;
<a name="l00391"></a>00391
<a name="l00392"></a>00392 <span class="comment">/* Write Analyze header */</span>
<a name="l00393"></a>00393 ret=<a class="code" href="analyze_8c.html#a3ec6f97c394c1bafffc57bd45d81c078">anaWriteHeader</a>(hdrfile, &dsr);
<a name="l00394"></a>00394 <span class="keywordflow">if</span>(ret) {
<a name="l00395"></a>00395 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bac4469979d9ee5c8d7ed3f63b82fcdf41">STATUS_CANTWRITEHEADERFILE</a>);
<a name="l00396"></a>00396 <span class="keywordflow">return</span> 21;
<a name="l00397"></a>00397 }
<a name="l00398"></a>00398 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00399"></a>00399 <span class="keywordflow">return</span> 0;
<a name="l00400"></a>00400 }
<a name="l00401"></a>00401 <span class="comment">/*****************************************************************************/</span>
<a name="l00402"></a>00402
<a name="l00403"></a>00403 <span class="comment">/*****************************************************************************/</span>
<a name="l00414"></a><a class="code" href="imgfile_8h.html#aca1bc7e49593e3137706287dfb9707c7">00414</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#aca1bc7e49593e3137706287dfb9707c7">imgReadAnalyzeHeader</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *dbname, <a class="code" href="structIMG.html">IMG</a> *img) {
<a name="l00415"></a>00415 <span class="keywordtype">char</span> hdrfile[FILENAME_MAX], siffile[FILENAME_MAX];
<a name="l00416"></a>00416 <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> ana_header;
<a name="l00417"></a>00417 <a class="code" href="structSIF.html">SIF</a> sif;
<a name="l00418"></a>00418 <span class="keywordtype">double</span> f;
<a name="l00419"></a>00419 <span class="keywordtype">int</span> ret;
<a name="l00420"></a>00420
<a name="l00421"></a>00421 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgReadAnalyzeHeader(%s, *img)\n"</span>, dbname);
<a name="l00422"></a>00422
<a name="l00423"></a>00423 <span class="comment">/* Check the input */</span>
<a name="l00424"></a>00424 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00425"></a>00425 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a7aeab2fb426b858c78ec515aaa8f0a67">IMG_STATUS_INITIALIZED</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00426"></a>00426 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00427"></a>00427 <span class="keywordflow">if</span>(dbname==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00428"></a>00428
<a name="l00429"></a>00429 <span class="comment">/* Determine the names of hdr and sif files */</span>
<a name="l00430"></a>00430 ret=<a class="code" href="analyze_8c.html#a6ecb9b0f8fb5f1fa4920a378d2ce0534">anaDatabaseExists</a>(dbname, hdrfile, NULL, siffile);
<a name="l00431"></a>00431 <span class="keywordflow">if</span>(ret==0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba9f298b5339db644dbc2b9dc01ef12006">STATUS_NOFILE</a>;
<a name="l00432"></a>00432
<a name="l00433"></a>00433 <span class="comment">/* Read Analyze header file */</span>
<a name="l00434"></a>00434 ret=<a class="code" href="analyze_8c.html#aa7bf1965dcdbff0686193675a3175a1c">anaReadHeader</a>(hdrfile, &ana_header);
<a name="l00435"></a>00435 <span class="keywordflow">if</span>(ret!=0) {
<a name="l00436"></a>00436 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>>1) printf(<span class="stringliteral">"anaReadHeader() return value := %d\n"</span>, ret);
<a name="l00437"></a>00437 <span class="keywordflow">if</span>(ret==1) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00438"></a>00438 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(ret==2) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba201ea0defe535c9b0449047fd6337370">STATUS_NOHEADERFILE</a>;
<a name="l00439"></a>00439 <span class="keywordflow">else</span> <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baef2b214bb9b6ccc009f27c76e81b3fef">STATUS_UNSUPPORTED</a>;
<a name="l00440"></a>00440 <span class="keywordflow">return</span>(<a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00441"></a>00441 }
<a name="l00442"></a>00442 <span class="comment">/* and set IMG contents */</span>
<a name="l00443"></a>00443 ret=<a class="code" href="img__ana_8c.html#a825dd53b92d0a2e0ec0e38308132be9f">imgGetAnalyzeHeader</a>(img, &ana_header);
<a name="l00444"></a>00444 <span class="keywordflow">if</span>(ret!=0) {
<a name="l00445"></a>00445 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, ret);
<a name="l00446"></a>00446 <span class="keywordflow">return</span>(ret);
<a name="l00447"></a>00447 }
<a name="l00448"></a>00448
<a name="l00449"></a>00449 <span class="comment">/* If SIF does not exist, then that's it */</span>
<a name="l00450"></a>00450 <span class="keywordflow">if</span>(!siffile[0]) {
<a name="l00451"></a>00451 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00452"></a>00452 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00453"></a>00453 }
<a name="l00454"></a>00454
<a name="l00455"></a>00455 <span class="comment">/* SIF is available, so read that too */</span>
<a name="l00456"></a>00456 <a class="code" href="sif_8c.html#a61c0218e380da8ace65d78f82a8f9f69">sifInit</a>(&sif); ret=0;
<a name="l00457"></a>00457 <span class="keywordflow">if</span>(<a class="code" href="sifio_8c.html#ae6c3f31d7cbe31ca6b74760978fd39c7">sifRead</a>(siffile, &sif)!=0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00458"></a>00458 <span class="comment">/* Copy scan time */</span>
<a name="l00459"></a>00459 img-><a class="code" href="structIMG.html#a5db78f8b423fc42b80bba42e27c7b2c3">scanStart</a>=sif.<a class="code" href="structSIF.html#a98b5ff616033e728f0edc9bdc7eddce2">scantime</a>;
<a name="l00460"></a>00460 <span class="comment">/* Study number, if not yet defined */</span>
<a name="l00461"></a>00461 <span class="keywordflow">if</span>(!img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>[0] && strlen(sif.<a class="code" href="structSIF.html#a2d3e6df8c97d55d0c28d83501d12f523">studynr</a>)>1 )
<a name="l00462"></a>00462 strncpy(img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>, sif.<a class="code" href="structSIF.html#a2d3e6df8c97d55d0c28d83501d12f523">studynr</a>, MAX_STUDYNR_LEN);
<a name="l00463"></a>00463 <span class="comment">/* Isotope half-life, if not yet defined */</span>
<a name="l00464"></a>00464 f=hlFromIsotope(sif.<a class="code" href="structSIF.html#aaf18449475434d1cf5ab11b2c36780b1">isotope_name</a>);
<a name="l00465"></a>00465 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a><=0.0 && f>0.0) img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>=60.0*f;
<a name="l00466"></a>00466 <a class="code" href="sif_8c.html#ad2378bec7fec69900af6f52737c89784">sifEmpty</a>(&sif);
<a name="l00467"></a>00467
<a name="l00468"></a>00468 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00469"></a>00469 }
<a name="l00470"></a>00470 <span class="comment">/*****************************************************************************/</span>
<a name="l00471"></a>00471
<a name="l00472"></a>00472 <span class="comment">/*****************************************************************************/</span>
<a name="l00481"></a><a class="code" href="imgfile_8h.html#a825dd53b92d0a2e0ec0e38308132be9f">00481</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#a825dd53b92d0a2e0ec0e38308132be9f">imgGetAnalyzeHeader</a>(<a class="code" href="structIMG.html">IMG</a> *img, <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> *h) {
<a name="l00482"></a>00482 <span class="keywordtype">int</span> dimNr, dimx, dimy, dimz=1, dimt=1, pxlNr=0;
<a name="l00483"></a>00483
<a name="l00484"></a>00484 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgGetAnalyzeHeader(*img, *dsr)\n"</span>);
<a name="l00485"></a>00485
<a name="l00486"></a>00486 <span class="comment">/* Check the input */</span>
<a name="l00487"></a>00487
<a name="l00488"></a>00488 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00489"></a>00489 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a7aeab2fb426b858c78ec515aaa8f0a67">IMG_STATUS_INITIALIZED</a> && img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a>)
<a name="l00490"></a>00490 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00491"></a>00491 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00492"></a>00492 <span class="keywordflow">if</span>(h==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00493"></a>00493
<a name="l00494"></a>00494 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7d0a6ddf0a60e475f14f6b0e9b6487a1">STATUS_INVALIDHEADER</a>);
<a name="l00495"></a>00495
<a name="l00496"></a>00496 <span class="comment">/* Get the image dimensions from header */</span>
<a name="l00497"></a>00497 dimNr=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[0];
<a name="l00498"></a>00498 <span class="keywordflow">if</span>(dimNr<2) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7d0a6ddf0a60e475f14f6b0e9b6487a1">STATUS_INVALIDHEADER</a>;
<a name="l00499"></a>00499 dimx=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[1]; dimy=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[2];
<a name="l00500"></a>00500 <span class="keywordflow">if</span>(dimNr>2) {dimz=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[3]; <span class="keywordflow">if</span>(dimNr>3) dimt=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4];}
<a name="l00501"></a>00501 pxlNr=dimx*dimy*dimz;
<a name="l00502"></a>00502 <span class="keywordflow">if</span>(pxlNr<1) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7d0a6ddf0a60e475f14f6b0e9b6487a1">STATUS_INVALIDHEADER</a>;
<a name="l00503"></a>00503 img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>=dimx; img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>=dimy; img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>=dimz; img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>=dimt;
<a name="l00504"></a>00504
<a name="l00505"></a>00505 <span class="comment">/* Copy information from Analyze header */</span>
<a name="l00506"></a>00506 img-><a class="code" href="structIMG.html#a0d9199e1c3c4bd94d56eb30e421b309c">type</a>=<a class="code" href="img_8h.html#ae8b5bb16b84adcbcd2bba3e10a36d62d">IMG_TYPE_IMAGE</a>;
<a name="l00507"></a>00507 strncpy(img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>, h-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ab9d251da82c987e8fe2eb23c6cf83f57">patient_id</a>, MAX_STUDYNR_LEN);
<a name="l00508"></a>00508 strcpy(img-><a class="code" href="structIMG.html#a29013e391052918bd5aa6974eeec811d">patientName</a>, h-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ab9d251da82c987e8fe2eb23c6cf83f57">patient_id</a>);
<a name="l00509"></a>00509 img-><a class="code" href="structIMG.html#a02cc239e917e7ebf85e34f21448ff080">sizex</a>=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[1];
<a name="l00510"></a>00510 img-><a class="code" href="structIMG.html#a80b0609c63f33934bf6e6678c6121cb3">sizey</a>=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[2];
<a name="l00511"></a>00511 img-><a class="code" href="structIMG.html#afc2a1d2417a8a3afe2c500173a05c472">sizez</a>=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[3];
<a name="l00512"></a>00512 <span class="comment">/*if(h->dime.funused2>1.E-5) img->zoom=h->dime.funused2;*/</span>
<a name="l00513"></a>00513 <span class="keywordflow">if</span>(h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>>1.E-5) img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>=h-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>;
<a name="l00514"></a>00514 <span class="keywordflow">if</span>(h-><a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>) img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>=<a class="code" href="img_8h.html#a38e113a0101e5e683653263ddabe0fc8">IMG_ANA_L</a>; <span class="keywordflow">else</span> img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>=<a class="code" href="img_8h.html#a9145a82aa093b7e52a693ec5c38094ef">IMG_ANA</a>;
<a name="l00515"></a>00515
<a name="l00516"></a>00516 <span class="comment">/* Decay correction */</span>
<a name="l00517"></a>00517 <span class="keywordflow">if</span>(strstr(h-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"Decay corrected."</span>)!=NULL)
<a name="l00518"></a>00518 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=1;
<a name="l00519"></a>00519 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(strstr(h-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"No decay correction."</span>)!=NULL)
<a name="l00520"></a>00520 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=0;
<a name="l00521"></a>00521 <span class="keywordflow">else</span>
<a name="l00522"></a>00522 img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>=1;
<a name="l00523"></a>00523
<a name="l00524"></a>00524 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00525"></a>00525 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00526"></a>00526 }
<a name="l00527"></a>00527 <span class="comment">/*****************************************************************************/</span>
<a name="l00528"></a>00528
<a name="l00529"></a>00529 <span class="comment">/*****************************************************************************/</span>
<a name="l00542"></a><a class="code" href="imgfile_8h.html#a234e8f67805da1e1c6f9f5361387b603">00542</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#adac51c408702b1d888a7056d2e764f9e">imgSetAnalyzeHeader</a>(<a class="code" href="structIMG.html">IMG</a> *img, <span class="keyword">const</span> <span class="keywordtype">char</span> *dbname, <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> *dsr, <span class="keywordtype">float</span> fmin, <span class="keywordtype">float</span> fmax) {
<a name="l00543"></a>00543 <span class="keyword">struct </span>tm *st;
<a name="l00544"></a>00544 <span class="keywordtype">char</span> *cptr;
<a name="l00545"></a>00545 <span class="keywordtype">float</span> g;
<a name="l00546"></a>00546
<a name="l00547"></a>00547 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgSetAnalyzeHeader(*img, *dsr)\n"</span>);
<a name="l00548"></a>00548
<a name="l00549"></a>00549 <span class="comment">/* Check the input */</span>
<a name="l00550"></a>00550 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00551"></a>00551 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a7aeab2fb426b858c78ec515aaa8f0a67">IMG_STATUS_INITIALIZED</a> && img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a>)
<a name="l00552"></a>00552 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00553"></a>00553 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00554"></a>00554 <span class="keywordflow">if</span>(dsr==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00555"></a>00555
<a name="l00556"></a>00556 <span class="comment">/* Byte order */</span>
<a name="l00557"></a>00557 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>==<a class="code" href="img_8h.html#a38e113a0101e5e683653263ddabe0fc8">IMG_ANA_L</a>) dsr-><a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>=1; <span class="keywordflow">else</span> dsr-><a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>=0;
<a name="l00558"></a>00558 <span class="comment">/* Header key */</span>
<a name="l00559"></a>00559 memset(&dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__KEY.html">ANALYZE_HEADER_KEY</a>));
<a name="l00560"></a>00560 memset(&dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__IMGDIM.html">ANALYZE_HEADER_IMGDIM</a>));
<a name="l00561"></a>00561 memset(&dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>, 0, <span class="keyword">sizeof</span>(<a class="code" href="structANALYZE__HEADER__HISTORY.html">ANALYZE_HEADER_HISTORY</a>));
<a name="l00562"></a>00562 dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a6d7253a777a1f2c86fa14e77a5c47515">sizeof_hdr</a>=348;
<a name="l00563"></a>00563 strcpy(dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a96f4c147959660ddbd2a579e19e24c3b">data_type</a>, <span class="stringliteral">""</span>);
<a name="l00564"></a>00564 cptr=strrchr(dbname, <span class="charliteral">'/'</span>); <span class="keywordflow">if</span>(cptr==NULL) cptr=strrchr(dbname, <span class="charliteral">'\\'</span>);
<a name="l00565"></a>00565 <span class="keywordflow">if</span>(cptr!=NULL) cptr++; <span class="keywordflow">if</span>(cptr==NULL) cptr=(<span class="keywordtype">char</span>*)dbname;
<a name="l00566"></a>00566 strncpy(dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a86cb24f71952964e9ebed40133da6725">db_name</a>, cptr, 17);
<a name="l00567"></a>00567 dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a09532be9a54a2e27170c8776c0ae1b3a">extents</a>=16384;
<a name="l00568"></a>00568 dsr-><a class="code" href="structANALYZE__DSR.html#a5477228c28d719c72bc01fc847cb4796">hk</a>.<a class="code" href="structANALYZE__HEADER__KEY.html#a6be3dcc756c59c1ddec21ca6311f2c5a">regular</a>=<span class="charliteral">'r'</span>;
<a name="l00569"></a>00569 <span class="comment">/* Image dimension */</span>
<a name="l00570"></a>00570 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[0]=4;
<a name="l00571"></a>00571 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[1]=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>;
<a name="l00572"></a>00572 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[2]=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>;
<a name="l00573"></a>00573 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[3]=img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>;
<a name="l00574"></a>00574 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]=img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>;
<a name="l00575"></a>00575 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a2a05d2cb24eeac5ac09098974f6b0fc9">datatype</a>=<a class="code" href="analyze_8h.html#a57a617cb3056b972eabcb94cdc289298">ANALYZE_DT_SIGNED_SHORT</a>;
<a name="l00576"></a>00576 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#ad0ca8c2b6ac4f211ff570c186979bbf9">bitpix</a>=16;
<a name="l00577"></a>00577 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[0]=0.0;
<a name="l00578"></a>00578 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[1]=img-><a class="code" href="structIMG.html#a02cc239e917e7ebf85e34f21448ff080">sizex</a>;
<a name="l00579"></a>00579 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[2]=img-><a class="code" href="structIMG.html#a80b0609c63f33934bf6e6678c6121cb3">sizey</a>;
<a name="l00580"></a>00580 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[3]=img-><a class="code" href="structIMG.html#afc2a1d2417a8a3afe2c500173a05c472">sizez</a>;
<a name="l00581"></a>00581 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a02e4cc5653e6ca7c740fa2f5cf0156d5">pixdim</a>[4]=0.0;
<a name="l00582"></a>00582 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>=0.0; <span class="comment">/* Scale factor is set later */</span>
<a name="l00583"></a>00583 <span class="comment">/* dsr.dime.funused2=img->zoom; */</span> <span class="comment">/* Reconstruction zoom */</span>
<a name="l00584"></a>00584 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a34b35e530052db74f16fb51001d8a27b">funused3</a>=img-><a class="code" href="structIMG.html#ac9c048a9dc8205b22518b05486d8030b">isotopeHalflife</a>;
<a name="l00585"></a>00585 <span class="comment">/* Data history */</span>
<a name="l00586"></a>00586 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a6ce801dea67ec8224e8cfc9aebb8a909">decayCorrected</a>==1) strcpy(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"Decay corrected."</span>);
<a name="l00587"></a>00587 <span class="keywordflow">else</span> strcpy(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#ad52cfe59c4eb1beb8b90bfc73d2fc571">descrip</a>, <span class="stringliteral">"No decay correction."</span>);
<a name="l00588"></a>00588 strncpy(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a0a750bf7ef6504c76d5c933debe56959">scannum</a>, img-><a class="code" href="structIMG.html#ad76f79021c93175ec98839613ba4dcd2">studyNr</a>, 10);
<a name="l00589"></a>00589 st=localtime(&img-><a class="code" href="structIMG.html#a5db78f8b423fc42b80bba42e27c7b2c3">scanStart</a>);
<a name="l00590"></a>00590 <span class="keywordflow">if</span>(st!=NULL) {
<a name="l00591"></a>00591 strftime(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a682bbce3bd0941ee16b62069e4cd1f15">exp_date</a>, 10, <span class="stringliteral">"%Y-%m-%d"</span>, st);
<a name="l00592"></a>00592 strftime(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a3938bc6e0b7a1770366b838d1d96c37e">exp_time</a>, 10, <span class="stringliteral">"%H:%M:%S"</span>, st);
<a name="l00593"></a>00593 } <span class="keywordflow">else</span> {
<a name="l00594"></a>00594 strcpy(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a682bbce3bd0941ee16b62069e4cd1f15">exp_date</a>, <span class="stringliteral">"1900-01-01"</span>);
<a name="l00595"></a>00595 strcpy(dsr-><a class="code" href="structANALYZE__DSR.html#a0cfb3442d9752e0bae02673441aadf60">hist</a>.<a class="code" href="structANALYZE__HEADER__HISTORY.html#a3938bc6e0b7a1770366b838d1d96c37e">exp_time</a>, <span class="stringliteral">"00:00:00"</span>);
<a name="l00596"></a>00596 }
<a name="l00597"></a>00597 <span class="comment">/* Determine and set scale factor and cal_min & cal_max */</span>
<a name="l00598"></a>00598 <span class="keywordflow">if</span>(fmin<fmax) {
<a name="l00599"></a>00599 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>=fmin; dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>=fmax;
<a name="l00600"></a>00600 } <span class="keywordflow">else</span> { <span class="comment">/* not given in function call, try to find those here */</span>
<a name="l00601"></a>00601 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>==<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a> &&
<a name="l00602"></a>00602 <a class="code" href="imgmax_8c.html#a97c48d1ed44d81b2d96ce1a139b99009">imgMinMax</a>(img, &dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>, &dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>)==0) {}
<a name="l00603"></a>00603 <span class="keywordflow">else</span> <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00604"></a>00604 }
<a name="l00605"></a>00605 <span class="keywordflow">if</span>(fabs(dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>) > fabs(dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>)) g = fabs(dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a75debe576328454d2b266a4f881ccccd">cal_min</a>);
<a name="l00606"></a>00606 <span class="keywordflow">else</span> g = fabs(dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a4b6e880b44f61ec68395f8686ab81f1e">cal_max</a>);
<a name="l00607"></a>00607 <span class="comment">/* if(fabs(dsr->dime.cal_min)>fabs(dsr->dime.cal_max)) g=fabs(dsr->dime.cal_min); */</span>
<a name="l00608"></a>00608 <span class="comment">/* else g=fabs(dsr->dime.cal_max); */</span>
<a name="l00609"></a>00609 <span class="keywordflow">if</span>(g<1E-20) g=1.0; <span class="keywordflow">else</span> g=32767./g; dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>=1.0/g;
<a name="l00610"></a>00610 <span class="comment">/* Set header glmin & glmax */</span>
<a name="l00611"></a>00611 dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aa92cd82aae233799feeea86dfdd0beb4">glmin</a>=temp_roundf(fmin*g); dsr-><a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a57093cdaa14482a67fcc4a29ae7c0cb9">glmax</a>=temp_roundf(fmax*g);
<a name="l00612"></a>00612 <span class="comment">/* printf("glmin=%d\n", dsr->dime.glmin); */</span>
<a name="l00613"></a>00613 <span class="comment">/* printf("glmax=%d\n", dsr->dime.glmax); */</span>
<a name="l00614"></a>00614
<a name="l00615"></a>00615 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00616"></a>00616 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00617"></a>00617 }
<a name="l00618"></a>00618 <span class="comment">/*****************************************************************************/</span>
<a name="l00619"></a>00619
<a name="l00620"></a>00620 <span class="comment">/*****************************************************************************/</span>
<a name="l00629"></a><a class="code" href="imgfile_8h.html#a8b134310e111af2b0577f5e95dc22e21">00629</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#a8b134310e111af2b0577f5e95dc22e21">imgReadAnalyzeFirstFrame</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *fname, <a class="code" href="structIMG.html">IMG</a> *img) {
<a name="l00630"></a>00630 <span class="keywordtype">int</span> ret=0;
<a name="l00631"></a>00631
<a name="l00632"></a>00632 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgReadAnalyzeFirstFrame(%s, *img)\n"</span>, fname);
<a name="l00633"></a>00633 <span class="comment">/* Check the input */</span>
<a name="l00634"></a>00634 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00635"></a>00635 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a7aeab2fb426b858c78ec515aaa8f0a67">IMG_STATUS_INITIALIZED</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00636"></a>00636 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00637"></a>00637 <span class="keywordflow">if</span>(fname==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00638"></a>00638
<a name="l00639"></a>00639 <span class="comment">/* Read header information from file */</span>
<a name="l00640"></a>00640 ret=<a class="code" href="img__ana_8c.html#aca1bc7e49593e3137706287dfb9707c7">imgReadAnalyzeHeader</a>(fname, img); <span class="keywordflow">if</span>(ret) <span class="keywordflow">return</span>(ret);
<a name="l00641"></a>00641 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>>3) <a class="code" href="img_8c.html#a32edd1b78039f6cc22b7bbe413aa3a56">imgInfo</a>(img);
<a name="l00642"></a>00642
<a name="l00643"></a>00643 <span class="comment">/* Allocate memory for one frame */</span>
<a name="l00644"></a>00644 img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>=1;
<a name="l00645"></a>00645 ret=<a class="code" href="img_8c.html#a9b2670da9c7685e6d53058cdf1db1199">imgAllocate</a>(img, img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>, img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>, img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>, img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>);
<a name="l00646"></a>00646 <span class="keywordflow">if</span>(ret) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>;
<a name="l00647"></a>00647
<a name="l00648"></a>00648 <span class="comment">/* Read the first frame */</span>
<a name="l00649"></a>00649 ret=<a class="code" href="img__ana_8c.html#a9094bdc56d1564b7ad412f0dfc114d9b">imgReadAnalyzeFrame</a>(fname, 1, img, 0); <span class="keywordflow">if</span>(ret) <span class="keywordflow">return</span>(ret);
<a name="l00650"></a>00650
<a name="l00651"></a>00651 <span class="comment">/* All went well */</span>
<a name="l00652"></a>00652 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>);
<a name="l00653"></a>00653 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00654"></a>00654 }
<a name="l00655"></a>00655 <span class="comment">/*****************************************************************************/</span>
<a name="l00656"></a>00656
<a name="l00657"></a>00657 <span class="comment">/*****************************************************************************/</span>
<a name="l00674"></a><a class="code" href="imgfile_8h.html#ac9d774f38a77bf88db92bb943c3ab324">00674</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#a9094bdc56d1564b7ad412f0dfc114d9b">imgReadAnalyzeFrame</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *fname, <span class="keywordtype">int</span> frame_to_read, <a class="code" href="structIMG.html">IMG</a> *img, <span class="keywordtype">int</span> frame_index) {
<a name="l00675"></a>00675 FILE *fp;
<a name="l00676"></a>00676 <span class="keywordtype">int</span> ret, pi, xi, yi;
<a name="l00677"></a>00677 <span class="keywordtype">float</span> *fdata=NULL, *fptr;
<a name="l00678"></a>00678 <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> dsr;
<a name="l00679"></a>00679 <span class="keywordtype">char</span> datfile[FILENAME_MAX], hdrfile[FILENAME_MAX], siffile[FILENAME_MAX];
<a name="l00680"></a>00680 <a class="code" href="structSIF.html">SIF</a> sif;
<a name="l00681"></a>00681
<a name="l00682"></a>00682
<a name="l00683"></a>00683 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgReadAnalyzeFrame(%s, %d, *img, %d)\n"</span>,
<a name="l00684"></a>00684 fname, frame_to_read, frame_index);
<a name="l00685"></a>00685
<a name="l00686"></a>00686 <span class="comment">/* Check the input */</span>
<a name="l00687"></a>00687 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00688"></a>00688 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00689"></a>00689 <span class="keywordflow">if</span>(fname==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00690"></a>00690 <span class="keywordflow">if</span>(frame_index<0 || frame_index>img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>-1) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00691"></a>00691 <span class="keywordflow">if</span>(frame_to_read<1) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00692"></a>00692 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00693"></a>00693
<a name="l00694"></a>00694 <span class="comment">/* Determine the names of hdr, data and sif files */</span>
<a name="l00695"></a>00695 ret=<a class="code" href="analyze_8c.html#a6ecb9b0f8fb5f1fa4920a378d2ce0534">anaDatabaseExists</a>(fname, hdrfile, datfile, siffile);
<a name="l00696"></a>00696 <span class="keywordflow">if</span>(ret==0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba9f298b5339db644dbc2b9dc01ef12006">STATUS_NOFILE</a>;
<a name="l00697"></a>00697
<a name="l00698"></a>00698 <span class="comment">/* Read Analyze header file */</span>
<a name="l00699"></a>00699 ret=<a class="code" href="analyze_8c.html#aa7bf1965dcdbff0686193675a3175a1c">anaReadHeader</a>(hdrfile, &dsr);
<a name="l00700"></a>00700 <span class="keywordflow">if</span>(ret!=0) {
<a name="l00701"></a>00701 <span class="keywordflow">if</span>(ret==1) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00702"></a>00702 <span class="keywordflow">else</span> <span class="keywordflow">if</span>(ret==2) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba201ea0defe535c9b0449047fd6337370">STATUS_NOHEADERFILE</a>;
<a name="l00703"></a>00703 <span class="keywordflow">else</span> <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baef2b214bb9b6ccc009f27c76e81b3fef">STATUS_UNSUPPORTED</a>;
<a name="l00704"></a>00704 <span class="keywordflow">return</span>(<a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>);
<a name="l00705"></a>00705 }
<a name="l00706"></a>00706
<a name="l00707"></a>00707 <span class="comment">/* Open image datafile */</span>
<a name="l00708"></a>00708 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>>2) fprintf(stdout, <span class="stringliteral">"reading image data %s\n"</span>, datfile);
<a name="l00709"></a>00709 <span class="keywordflow">if</span>((fp=fopen(datfile, <span class="stringliteral">"rb"</span>)) == NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba524fb2fd0efecb30fe42b109b3d28fe8">STATUS_NOIMGDATA</a>;
<a name="l00710"></a>00710
<a name="l00711"></a>00711 <span class="comment">/* Allocate memory for one image frame */</span>
<a name="l00712"></a>00712 fdata=malloc(img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>*img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>*img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>));
<a name="l00713"></a>00713 <span class="keywordflow">if</span>(fdata==NULL) {fclose(fp); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>;}
<a name="l00714"></a>00714
<a name="l00715"></a>00715 <span class="comment">/* Read the required image frame */</span>
<a name="l00716"></a>00716 fptr=fdata;
<a name="l00717"></a>00717 ret=<a class="code" href="analyze_8c.html#abfb69549a8e8c6a4eab3c52045262d22">anaReadImagedata</a>(fp, &dsr, frame_to_read, fptr);
<a name="l00718"></a>00718 fclose(fp);
<a name="l00719"></a>00719 <span class="keywordflow">if</span>(ret==3) {free(fdata); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55babe9c1bcf366cd905b79148ba91b34756">STATUS_NOMATRIX</a>;} <span class="comment">/* no more frames */</span>
<a name="l00720"></a>00720 <span class="keywordflow">if</span>(ret!=0) {free(fdata); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baef2b214bb9b6ccc009f27c76e81b3fef">STATUS_UNSUPPORTED</a>;}
<a name="l00721"></a>00721
<a name="l00722"></a>00722 <span class="comment">/* Copy pixel values to IMG */</span>
<a name="l00723"></a>00723 fptr=fdata;
<a name="l00724"></a>00724 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a55febab802306e82d162793558f5fbce">anaFlipping</a>()==0) { <span class="comment">/* no flipping in z-direction */</span>
<a name="l00725"></a>00725 <span class="keywordflow">for</span>(pi=0; pi<img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>; pi++)
<a name="l00726"></a>00726 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--)
<a name="l00727"></a>00727 <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--)
<a name="l00728"></a>00728 img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][frame_index]=*fptr++;
<a name="l00729"></a>00729 } <span class="keywordflow">else</span> {
<a name="l00730"></a>00730 <span class="keywordflow">for</span>(pi=img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>-1; pi>=0; pi--)
<a name="l00731"></a>00731 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--)
<a name="l00732"></a>00732 <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--)
<a name="l00733"></a>00733 img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[pi][yi][xi][frame_index]=*fptr++;
<a name="l00734"></a>00734 }
<a name="l00735"></a>00735 free(fdata);
<a name="l00736"></a>00736
<a name="l00737"></a>00737 <a class="code" href="img_8c.html#ad124d2d496934d677156cc8779a1e567">imgSetStatus</a>(img, <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>); <span class="comment">/* If the rest is failed, no problem */</span>
<a name="l00738"></a>00738
<a name="l00739"></a>00739 <span class="comment">/* </span>
<a name="l00740"></a>00740 <span class="comment"> * Try to read frame time information from SIF file</span>
<a name="l00741"></a>00741 <span class="comment"> */</span>
<a name="l00742"></a>00742 <a class="code" href="sif_8c.html#a61c0218e380da8ace65d78f82a8f9f69">sifInit</a>(&sif);
<a name="l00743"></a>00743 <span class="keywordflow">if</span>(<a class="code" href="sifio_8c.html#ae6c3f31d7cbe31ca6b74760978fd39c7">sifRead</a>(siffile, &sif)!=0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00744"></a>00744 <span class="comment">/* Frame information */</span>
<a name="l00745"></a>00745 <span class="keywordflow">if</span>(sif.<a class="code" href="structSIF.html#a441bd4d1d3825ace8f21d02e9f5d650d">frameNr</a>>=frame_to_read) {
<a name="l00746"></a>00746 img-><a class="code" href="structIMG.html#a8ff22fc209b6f44ad14341c8644788d2">start</a>[frame_index]=sif.<a class="code" href="structSIF.html#a085f3f1a4b2a77eabe13f52c2010a3c0">x1</a>[frame_to_read-1];
<a name="l00747"></a>00747 img-><a class="code" href="structIMG.html#aa2b629a51ca62bcafad0a66d32e1ca8b">end</a>[frame_index]=sif.<a class="code" href="structSIF.html#a45001ad49c925e6ae266357b3ab7d8a4">x2</a>[frame_to_read-1];
<a name="l00748"></a>00748 img-><a class="code" href="structIMG.html#af7c69142be89fac081e7f8bc354b14c8">mid</a>[frame_index]=0.5*(img-><a class="code" href="structIMG.html#a8ff22fc209b6f44ad14341c8644788d2">start</a>[frame_index]+img-><a class="code" href="structIMG.html#aa2b629a51ca62bcafad0a66d32e1ca8b">end</a>[frame_index]);
<a name="l00749"></a>00749 img-><a class="code" href="structIMG.html#a76d63a001b8569d9f80ccb325e064be1">prompts</a>[frame_index]=sif.<a class="code" href="structSIF.html#a20d6fecbd751d872300c661e19903eed">prompts</a>[frame_to_read-1];
<a name="l00750"></a>00750 img-><a class="code" href="structIMG.html#ae47ebbbdc0c308f473fec801a2df513d">randoms</a>[frame_index]=sif.<a class="code" href="structSIF.html#af2fa6ee9e91796baa37261bd866b641f">randoms</a>[frame_to_read-1];
<a name="l00751"></a>00751 }
<a name="l00752"></a>00752 <a class="code" href="sif_8c.html#ad2378bec7fec69900af6f52737c89784">sifEmpty</a>(&sif);
<a name="l00753"></a>00753
<a name="l00754"></a>00754 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00755"></a>00755 }
<a name="l00756"></a>00756 <span class="comment">/*****************************************************************************/</span>
<a name="l00757"></a>00757
<a name="l00758"></a>00758 <span class="comment">/*****************************************************************************/</span>
<a name="l00781"></a><a class="code" href="imgfile_8h.html#aa69ada7b39ad67366fb5b358198feec0">00781</a> <span class="keywordtype">int</span> <a class="code" href="img__ana_8c.html#add3b38cb5c7cb556ee55dee3b7d66130">imgWriteAnalyzeFrame</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *dbname, <span class="keywordtype">int</span> frame_to_write, <a class="code" href="structIMG.html">IMG</a> *img, <span class="keywordtype">int</span> frame_index, <span class="keywordtype">float</span> fmin, <span class="keywordtype">float</span> fmax) {
<a name="l00782"></a>00782 <a class="code" href="structIMG.html">IMG</a> test_img;
<a name="l00783"></a>00783 <span class="keywordtype">int</span> ret=0, pxlNr, zi, xi, yi, little;
<a name="l00784"></a>00784 FILE *fp;
<a name="l00785"></a>00785 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sdata=NULL, *sptr;
<a name="l00786"></a>00786 <span class="keywordtype">char</span> datfile[FILENAME_MAX], hdrfile[FILENAME_MAX], siffile[FILENAME_MAX];
<a name="l00787"></a>00787 <a class="code" href="structANALYZE__DSR.html">ANALYZE_DSR</a> dsr;
<a name="l00788"></a>00788 <span class="keywordtype">float</span> scale_factor=1.0;
<a name="l00789"></a>00789
<a name="l00790"></a>00790
<a name="l00791"></a>00791 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"\nimgWriteAnalyzeFrame(%s, %d, *img, %d, %g, %g)\n"</span>,
<a name="l00792"></a>00792 dbname, frame_to_write, frame_index, fmin, fmax);
<a name="l00793"></a>00793
<a name="l00794"></a>00794 <span class="comment">/*</span>
<a name="l00795"></a>00795 <span class="comment"> * Check the input </span>
<a name="l00796"></a>00796 <span class="comment"> */</span>
<a name="l00797"></a>00797 <span class="keywordflow">if</span>(dbname==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00798"></a>00798 <span class="keywordflow">if</span>(img==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00799"></a>00799 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5b8bdc61088df4a3c3f89c31895b14a4">status</a>!=<a class="code" href="img_8h.html#a3117b1b0649d3761dcda1395202129e5">IMG_STATUS_OCCUPIED</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00800"></a>00800 <span class="keywordflow">if</span>(frame_to_write<0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00801"></a>00801 <span class="keywordflow">if</span>(frame_index<0 || frame_index>=img-><a class="code" href="structIMG.html#a77605d4aa64b86294590c79e9a761b15">dimt</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00802"></a>00802 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>!=<a class="code" href="img_8h.html#a38e113a0101e5e683653263ddabe0fc8">IMG_ANA_L</a> && img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>!=<a class="code" href="img_8h.html#a9145a82aa093b7e52a693ec5c38094ef">IMG_ANA</a>) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baf8f3d89c25483639f6cbce4f6b08c1f8">STATUS_FAULT</a>;
<a name="l00803"></a>00803
<a name="l00804"></a>00804 <span class="comment">/*</span>
<a name="l00805"></a>00805 <span class="comment"> * If database does not exist, then create it with new header,</span>
<a name="l00806"></a>00806 <span class="comment"> * and if it does exist, then read and check header information.</span>
<a name="l00807"></a>00807 <span class="comment"> * Create or edit header to contain correct frame nr.</span>
<a name="l00808"></a>00808 <span class="comment"> * Determine the global scaling factor. </span>
<a name="l00809"></a>00809 <span class="comment"> */</span>
<a name="l00810"></a>00810 <a class="code" href="img_8c.html#af95041292a0f9e4a82a10b49a074a95c">imgInit</a>(&test_img);
<a name="l00811"></a>00811 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a6ecb9b0f8fb5f1fa4920a378d2ce0534">anaDatabaseExists</a>(dbname, hdrfile, datfile, siffile)==0) { <span class="comment">/* not existing */</span>
<a name="l00812"></a>00812
<a name="l00813"></a>00813 <span class="comment">/* Create database filenames */</span>
<a name="l00814"></a>00814 sprintf(hdrfile, <span class="stringliteral">"%s.hdr"</span>, dbname);
<a name="l00815"></a>00815 sprintf(datfile, <span class="stringliteral">"%s.img"</span>, dbname);
<a name="l00816"></a>00816 sprintf(siffile, <span class="stringliteral">"%s.sif"</span>, dbname);
<a name="l00817"></a>00817
<a name="l00818"></a>00818 <span class="comment">/* Set main header */</span>
<a name="l00819"></a>00819 <a class="code" href="img__ana_8c.html#adac51c408702b1d888a7056d2e764f9e">imgSetAnalyzeHeader</a>(img, dbname, &dsr, fmin, fmax);
<a name="l00820"></a>00820 <span class="keywordflow">if</span>(frame_to_write==0) frame_to_write=1;
<a name="l00821"></a>00821 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]=frame_to_write;
<a name="l00822"></a>00822 scale_factor=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>;
<a name="l00823"></a>00823 <span class="keywordflow">if</span>(fabs(scale_factor)>1.0E-20) scale_factor=1.0/scale_factor;
<a name="l00824"></a>00824
<a name="l00825"></a>00825 <span class="comment">/* Write Analyze header */</span>
<a name="l00826"></a>00826 ret=<a class="code" href="analyze_8c.html#a3ec6f97c394c1bafffc57bd45d81c078">anaWriteHeader</a>(hdrfile, &dsr); <span class="keywordflow">if</span>(ret && <a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>) printf(<span class="stringliteral">"anaWriteHeader() := %d\n"</span>, ret);
<a name="l00827"></a>00827 <span class="keywordflow">if</span>(ret) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bac4469979d9ee5c8d7ed3f63b82fcdf41">STATUS_CANTWRITEHEADERFILE</a>;
<a name="l00828"></a>00828
<a name="l00829"></a>00829 <span class="comment">/* Remove datafile if necessary */</span>
<a name="l00830"></a>00830 <span class="keywordflow">if</span>(access(datfile, 0) != -1) <span class="keyword">remove</span>(datfile);
<a name="l00831"></a>00831
<a name="l00832"></a>00832 } <span class="keywordflow">else</span> { <span class="comment">/* database does exist */</span>
<a name="l00833"></a>00833
<a name="l00834"></a>00834 <span class="comment">/* Read header information for checking */</span>
<a name="l00835"></a>00835 ret=<a class="code" href="img__ana_8c.html#aca1bc7e49593e3137706287dfb9707c7">imgReadAnalyzeHeader</a>(dbname, &test_img); <span class="keywordflow">if</span>(ret!=0) <span class="keywordflow">return</span> ret;
<a name="l00836"></a>00836 <span class="comment">/* Check that file format is the same */</span>
<a name="l00837"></a>00837 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a>!=test_img.<a class="code" href="structIMG.html#a5e25c84c23270f974f1eb7c32a01bac7">_fileFormat</a> || img-><a class="code" href="structIMG.html#a0d9199e1c3c4bd94d56eb30e421b309c">type</a>!=test_img.<a class="code" href="structIMG.html#a0d9199e1c3c4bd94d56eb30e421b309c">type</a>)
<a name="l00838"></a>00838 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba8860e85936fb99f4d0f248a9ca7fd4f7">STATUS_WRONGFILETYPE</a>;
<a name="l00839"></a>00839 <span class="comment">/* Check that matrix sizes are the same */</span>
<a name="l00840"></a>00840 <span class="keywordflow">if</span>(img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>!=test_img.<a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a> || img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>!=test_img.<a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a> ||
<a name="l00841"></a>00841 img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>!=test_img.<a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>)
<a name="l00842"></a>00842 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bae8ce9dbf165d2c5fdb9c4609c01d94e5">STATUS_VARMATSIZE</a>;
<a name="l00843"></a>00843 <a class="code" href="img_8c.html#ada573d4b40da77be4455bb9fab9a3033">imgEmpty</a>(&test_img);
<a name="l00844"></a>00844
<a name="l00845"></a>00845 <span class="comment">/* Read the header, set new frame number, and write it back */</span>
<a name="l00846"></a>00846 <span class="comment">/* Get also the scale factor */</span>
<a name="l00847"></a>00847 <span class="keywordflow">if</span>((ret=<a class="code" href="analyze_8c.html#aa7bf1965dcdbff0686193675a3175a1c">anaReadHeader</a>(hdrfile, &dsr))!=0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafa1051360f88591f14aa0383d8087132">STATUS_NOMAINHEADER</a>;
<a name="l00848"></a>00848 scale_factor=1.0/dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#a162fbe1e4494d6fa56091b447afaf4a0">funused1</a>;
<a name="l00849"></a>00849 <span class="keywordflow">if</span>(frame_to_write==0) frame_to_write=dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]+1;
<a name="l00850"></a>00850 <span class="keywordflow">if</span>(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]<frame_to_write) {
<a name="l00851"></a>00851 <span class="keywordflow">if</span>(dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]+1<frame_to_write) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bad6498890650ae8d01e69f43514eb4fed">STATUS_MISSINGMATRIX</a>;
<a name="l00852"></a>00852 dsr.<a class="code" href="structANALYZE__DSR.html#a6599d263ca760b52abba8f87657a88ba">dime</a>.<a class="code" href="structANALYZE__HEADER__IMGDIM.html#aaffbef7b5025d4ccb95161217424f84f">dim</a>[4]=frame_to_write;
<a name="l00853"></a>00853 }
<a name="l00854"></a>00854 <span class="keywordflow">if</span>((ret=<a class="code" href="analyze_8c.html#a3ec6f97c394c1bafffc57bd45d81c078">anaWriteHeader</a>(hdrfile, &dsr))!=0) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba89d8cfdc7190aa78887b412e66e292e4">STATUS_NOWRITEPERM</a>;
<a name="l00855"></a>00855 }
<a name="l00856"></a>00856 <span class="keywordflow">if</span>(<a class="code" href="img_8h.html#a0d66ed136e5992b7b00ce6f43bad1439">IMG_TEST</a>>2) {
<a name="l00857"></a>00857 printf(<span class="stringliteral">"frame_to_write := %d\n"</span>, frame_to_write);
<a name="l00858"></a>00858 printf(<span class="stringliteral">"hdrfile := %s\n"</span>, hdrfile);
<a name="l00859"></a>00859 printf(<span class="stringliteral">"datfile := %s\n"</span>, datfile);
<a name="l00860"></a>00860 printf(<span class="stringliteral">"siffile := %s\n"</span>, siffile);
<a name="l00861"></a>00861 }
<a name="l00862"></a>00862
<a name="l00863"></a>00863 <span class="comment">/* Allocate memory for matrix short int data (one plane) */</span>
<a name="l00864"></a>00864 pxlNr=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>*img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>;
<a name="l00865"></a>00865 sdata=(<span class="keywordtype">short</span> <span class="keywordtype">int</span>*)malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>));
<a name="l00866"></a>00866 <span class="keywordflow">if</span>(sdata==NULL) <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bafbe470f6663c748e405bf509d2bb57bb">STATUS_NOMEMORY</a>;
<a name="l00867"></a>00867
<a name="l00868"></a>00868 <span class="comment">/* Open datafile, not removing possible old contents */</span>
<a name="l00869"></a>00869 <span class="keywordflow">if</span>(frame_to_write==1) fp=fopen(datfile, <span class="stringliteral">"wb"</span>); <span class="keywordflow">else</span> fp=fopen(datfile, <span class="stringliteral">"r+b"</span>);
<a name="l00870"></a>00870 <span class="keywordflow">if</span>(fp==NULL) {free(sdata); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baa2334ba6708d894ba27671de890f0e0a">STATUS_CANTWRITEIMGFILE</a>;}
<a name="l00871"></a>00871 <span class="comment">/* Move file pointer to the place of current frame */</span>
<a name="l00872"></a>00872 <span class="keywordflow">if</span>(fseek(fp, (frame_to_write-1)*pxlNr*img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>), SEEK_SET)!=0) {
<a name="l00873"></a>00873 free(sdata); fclose(fp); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55bad6498890650ae8d01e69f43514eb4fed">STATUS_MISSINGMATRIX</a>;}
<a name="l00874"></a>00874 little=little_endian();
<a name="l00875"></a>00875 <span class="comment">/* Copy, scale and write data plane-by-plane */</span>
<a name="l00876"></a>00876 <span class="keywordflow">if</span>(<a class="code" href="analyze_8c.html#a55febab802306e82d162793558f5fbce">anaFlipping</a>()==0) {
<a name="l00877"></a>00877 <span class="keywordflow">for</span>(zi=0; zi<img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>; zi++) {
<a name="l00878"></a>00878 sptr=sdata; <span class="comment">/*printf("plane := %d\n scale_factor := %g\n", zi+1, scale_factor);*/</span>
<a name="l00879"></a>00879 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--) <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--) {
<a name="l00880"></a>00880 *sptr=temp_roundf(scale_factor*img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[zi][yi][xi][frame_index]); sptr++;
<a name="l00881"></a>00881 }
<a name="l00882"></a>00882 <span class="comment">/* Change byte order if necessary */</span>
<a name="l00883"></a>00883 sptr=sdata; <span class="keywordflow">if</span>(little!=dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>) swabip(sptr, pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>));
<a name="l00884"></a>00884 <span class="comment">/* Write image data */</span>
<a name="l00885"></a>00885 sptr=sdata;
<a name="l00886"></a>00886 <span class="keywordflow">if</span>(fwrite(sptr, <span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>), pxlNr, fp) != pxlNr) {
<a name="l00887"></a>00887 free(sdata); fclose(fp); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baa2334ba6708d894ba27671de890f0e0a">STATUS_CANTWRITEIMGFILE</a>;
<a name="l00888"></a>00888 }
<a name="l00889"></a>00889 }
<a name="l00890"></a>00890 } <span class="keywordflow">else</span> {
<a name="l00891"></a>00891 <span class="keywordflow">for</span>(zi=img-><a class="code" href="structIMG.html#a9c04780b47e6a627a2b700b76352c3fc">dimz</a>-1; zi>=0; zi--) {
<a name="l00892"></a>00892 sptr=sdata;
<a name="l00893"></a>00893 <span class="keywordflow">for</span>(yi=img-><a class="code" href="structIMG.html#aa17e87812503d3c142f31c08abcbf79d">dimy</a>-1; yi>=0; yi--) <span class="keywordflow">for</span>(xi=img-><a class="code" href="structIMG.html#a0818b82059710627a94282a793ff7030">dimx</a>-1; xi>=0; xi--) {
<a name="l00894"></a>00894 *sptr=temp_roundf(scale_factor*img-><a class="code" href="structIMG.html#a43e2c98dc0728731601a76abd51bb935">m</a>[zi][yi][xi][frame_index]); sptr++;
<a name="l00895"></a>00895 }
<a name="l00896"></a>00896 <span class="comment">/* Change byte order if necessary */</span>
<a name="l00897"></a>00897 sptr=sdata; <span class="keywordflow">if</span>(little!=dsr.<a class="code" href="structANALYZE__DSR.html#a4bc9f7d8293b4ac8ff45331cf58ca5df">little</a>) swabip(sptr, pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>));
<a name="l00898"></a>00898 <span class="comment">/* Write image data */</span>
<a name="l00899"></a>00899 sptr=sdata;
<a name="l00900"></a>00900 <span class="keywordflow">if</span>(fwrite(sptr, <span class="keyword">sizeof</span>(<span class="keywordtype">short</span> <span class="keywordtype">int</span>), pxlNr, fp) != pxlNr) {
<a name="l00901"></a>00901 free(sdata); fclose(fp); <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55baa2334ba6708d894ba27671de890f0e0a">STATUS_CANTWRITEIMGFILE</a>;
<a name="l00902"></a>00902 }
<a name="l00903"></a>00903 }
<a name="l00904"></a>00904 }
<a name="l00905"></a>00905 free(sdata);
<a name="l00906"></a>00906 fclose(fp);
<a name="l00907"></a>00907
<a name="l00908"></a>00908 <span class="keywordflow">return</span> <a class="code" href="img_8h.html#a06fc87d81c62e9abb8790b6e5713c55ba7e4a42e3b6dd63708c64cf3db6f69566">STATUS_OK</a>;
<a name="l00909"></a>00909 }
<a name="l00910"></a>00910 <span class="comment">/*****************************************************************************/</span>
<a name="l00911"></a>00911
<a name="l00912"></a>00912 <span class="comment">/*****************************************************************************/</span>
</pre></div></div>
</div>
<hr class="footer"/><address class="footer"><small>Generated on Mon Aug 8 2011 by 
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.4 </small></address>
</body>
</html>
