<!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>ecat7r.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">ecat7r.c</div> </div> </div> <div class="contents"> <a href="ecat7r_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) 2003-2010 Turku PET Centre</span> <a name="l00004"></a>00004 <span class="comment"></span> <a name="l00005"></a>00005 <span class="comment"> Library file: ecat7r.c</span> <a name="l00006"></a>00006 <span class="comment"> Description: Functions for reading ECAT 7.x 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"> 2003-07-24 Vesa Oikonen</span> <a name="l00027"></a>00027 <span class="comment"> First created.</span> <a name="l00028"></a>00028 <span class="comment"> 2003-09-08 VO</span> <a name="l00029"></a>00029 <span class="comment"> Added support for 3D sinograms, ecat7ReadScanMatrix().</span> <a name="l00030"></a>00030 <span class="comment"> 2004-05-23 VO</span> <a name="l00031"></a>00031 <span class="comment"> Comments changed into Doxygen format.</span> <a name="l00032"></a>00032 <span class="comment"> 2004-06-21 VO</span> <a name="l00033"></a>00033 <span class="comment"> ecat7ReadScanMatrix():</span> <a name="l00034"></a>00034 <span class="comment"> Before: reads datablocks based on matrix list.</span> <a name="l00035"></a>00035 <span class="comment"> After: if block number based on bin nr is smaller, then read only those</span> <a name="l00036"></a>00036 <span class="comment"> Reason: simulated file with erroneous matrix list.</span> <a name="l00037"></a>00037 <span class="comment"> 2004-09-20 VO</span> <a name="l00038"></a>00038 <span class="comment"> Doxygen style comments are corrected.</span> <a name="l00039"></a>00039 <span class="comment"> 2004-11-10 VO</span> <a name="l00040"></a>00040 <span class="comment"> Calculation of trueblockNr simplified in ecat7ReadScanMatrix().</span> <a name="l00041"></a>00041 <span class="comment"> 2006-02-07 Jarkko Johansson</span> <a name="l00042"></a>00042 <span class="comment"> Comments added in ecat7ReadScanMatrix().</span> <a name="l00043"></a>00043 <span class="comment"> 2007-03-21 VO </span> <a name="l00044"></a>00044 <span class="comment"> ecat7ReadImageheader(): fill_cti[] and fill_user[] are read correctly.</span> <a name="l00045"></a>00045 <span class="comment"> 2007-03-27 VO</span> <a name="l00046"></a>00046 <span class="comment"> Added ecat7ReadPolarmapMatrix().</span> <a name="l00047"></a>00047 <span class="comment"> 2010-08-19 VO</span> <a name="l00048"></a>00048 <span class="comment"> Main header field patient_birth_date can be in two different int formats,</span> <a name="l00049"></a>00049 <span class="comment"> either YYYYMMDD or as seconds from start of year 1970. In latter case</span> <a name="l00050"></a>00050 <span class="comment"> the number can be negative, which is not identified correctly by all C</span> <a name="l00051"></a>00051 <span class="comment"> library versions. Therefore those are converted to YYYYMMDD format.</span> <a name="l00052"></a>00052 <span class="comment"></span> <a name="l00053"></a>00053 <span class="comment"></span> <a name="l00054"></a>00054 <span class="comment">******************************************************************************/</span> <a name="l00055"></a>00055 <span class="preprocessor">#include <locale.h></span> <a name="l00056"></a>00056 <span class="preprocessor">#include <stdio.h></span> <a name="l00057"></a>00057 <span class="preprocessor">#include <stdlib.h></span> <a name="l00058"></a>00058 <span class="preprocessor">#include <math.h></span> <a name="l00059"></a>00059 <span class="preprocessor">#include <ctype.h></span> <a name="l00060"></a>00060 <span class="preprocessor">#include <string.h></span> <a name="l00061"></a>00061 <span class="preprocessor">#include <unistd.h></span> <a name="l00062"></a>00062 <span class="preprocessor">#include <time.h></span> <a name="l00063"></a>00063 <span class="comment">/*****************************************************************************/</span> <a name="l00064"></a>00064 <span class="preprocessor">#include <swap.h></span> <a name="l00065"></a>00065 <span class="preprocessor">#include <datetime.h></span> <a name="l00066"></a>00066 <span class="preprocessor">#include "<a class="code" href="ecat7_8h.html">include/ecat7.h</a>"</span> <a name="l00067"></a>00067 <span class="comment">/*****************************************************************************/</span> <a name="l00068"></a>00068 <a name="l00069"></a>00069 <span class="comment">/*****************************************************************************/</span> <a name="l00078"></a><a class="code" href="ecat7_8h.html#a7d6241dd7a472a4f2aaf3eb16dc03dc1">00078</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a7d6241dd7a472a4f2aaf3eb16dc03dc1">ecat7ReadMainheader</a>(FILE *fp, <a class="code" href="structecat7__mainheader.html">ECAT7_mainheader</a> *h) { <a name="l00079"></a>00079 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00080"></a>00080 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00081"></a>00081 <span class="keyword">struct </span>tm st; <a name="l00082"></a>00082 <a name="l00083"></a>00083 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadMainheader()\n"</span>); <a name="l00084"></a>00084 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00085"></a>00085 little=little_endian(); <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"little=%d\n"</span>, little); <a name="l00086"></a>00086 <a name="l00087"></a>00087 <span class="comment">/* Seek the first block */</span> <a name="l00088"></a>00088 fseek(fp, 0, SEEK_SET); <span class="keywordflow">if</span>(ftell(fp)!=0) <span class="keywordflow">return</span>(2); <a name="l00089"></a>00089 <span class="comment">/* Read the header block */</span> <a name="l00090"></a>00090 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00091"></a>00091 <a name="l00092"></a>00092 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00093"></a>00093 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a71cb64edce706f383f3e68b9548b9ec3">magic_number</a>, buf+0, 14); <a name="l00094"></a>00094 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a229fda43a3abcdb4e055f2b5791823a8">original_file_name</a>, buf+14, 32); <a name="l00095"></a>00095 <span class="keywordflow">if</span>(little) swabip(buf+46, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a413498a0395e39c969bf7bf401cf023c">sw_version</a>, buf+46, 2); <a name="l00096"></a>00096 <span class="keywordflow">if</span>(little) swabip(buf+48, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ae6537726ed92ab81d76c8985fff81244">system_type</a>, buf+48, 2); <a name="l00097"></a>00097 <span class="keywordflow">if</span>(little) swabip(buf+50, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a1914b3a0ebf3077d090752d6ef66c757">file_type</a>, buf+50, 2); <a name="l00098"></a>00098 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a03420554a59f09d747e217011f1859c0">serial_number</a>, buf+52, 10); <a name="l00099"></a>00099 <span class="keywordflow">if</span>(little) swawbip(buf+62, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a1ad956ffc4f0dc8e9bd473bfa20e6b81">scan_start_time</a>, buf+62, 4); <a name="l00100"></a>00100 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a1a6eea4ac705c0f84f755ac6cb12d69d">isotope_name</a>, buf+66, 8); <a name="l00101"></a>00101 <span class="keywordflow">if</span>(little) swawbip(buf+74, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a9389f1b044165caf1a4b79007cb6cdfd">isotope_halflife</a>, buf+74, 4); <a name="l00102"></a>00102 memcpy(&h-><a class="code" href="structecat7__mainheader.html#ab353f8a8b1166323fe7798d81004a225">radiopharmaceutical</a>, buf+78, 32); <a name="l00103"></a>00103 <span class="keywordflow">if</span>(little) swawbip(buf+110, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a5bcbe3f6bdf521d7bc709b79aa8359c0">gantry_tilt</a>, buf+110, 4); <a name="l00104"></a>00104 <span class="keywordflow">if</span>(little) swawbip(buf+114, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a81f5c30a8cfa68da62bd4121b4251e6f">gantry_rotation</a>, buf+114, 4); <a name="l00105"></a>00105 <span class="keywordflow">if</span>(little) swawbip(buf+118, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a8e0caba01dedd5e94dbd36d6db107149">bed_elevation</a>, buf+118, 4); <a name="l00106"></a>00106 <span class="keywordflow">if</span>(little) swawbip(buf+122, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ac39e16ca7d54c7c1ba069f3581700d16">intrinsic_tilt</a>, buf+122, 4); <a name="l00107"></a>00107 <span class="keywordflow">if</span>(little) swabip(buf+126, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a6bb42a878ca6292118b2e591eff6b965">wobble_speed</a>, buf+126, 2); <a name="l00108"></a>00108 <span class="keywordflow">if</span>(little) swabip(buf+128, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a637cb9f894c15a7ff4f147ac1810a80f">transm_source_type</a>, buf+128, 2); <a name="l00109"></a>00109 <span class="keywordflow">if</span>(little) swawbip(buf+130, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a1be6bc903aaef62d6cfac179b7b8a05b">distance_scanned</a>, buf+130, 4); <a name="l00110"></a>00110 <span class="keywordflow">if</span>(little) swawbip(buf+134, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a3dff6c0caf29b90298e542c6e76af2a1">transaxial_fov</a>, buf+134, 4); <a name="l00111"></a>00111 <span class="keywordflow">if</span>(little) swabip(buf+138, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a5682114a0723ab5c3eb65ef673f85779">angular_compression</a>, buf+138, 2); <a name="l00112"></a>00112 <span class="keywordflow">if</span>(little) swabip(buf+140, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a8b370daee9ea0f228870038831684a04">coin_samp_mode</a>, buf+140, 2); <a name="l00113"></a>00113 <span class="keywordflow">if</span>(little) swabip(buf+142, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a97197dc85925d45fc1bf501c2edb94cc">axial_samp_mode</a>, buf+142, 2); <a name="l00114"></a>00114 <span class="keywordflow">if</span>(little) swawbip(buf+144, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ad5be044d585696c1f35bb8b63c26fec9">ecat_calibration_factor</a>, buf+144, 4); <a name="l00115"></a>00115 <span class="keywordflow">if</span>(little) swabip(buf+148, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a56bb68347079aca9062f24fe0fb24a63">calibration_units</a>, buf+148, 2); <a name="l00116"></a>00116 <span class="keywordflow">if</span>(little) swabip(buf+150, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a9ee5c316fc447e305593efc25f51dc67">calibration_units_label</a>, buf+150, 2); <a name="l00117"></a>00117 <span class="keywordflow">if</span>(little) swabip(buf+152, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a11ae37f7e7af8371b4884aff851dec2b">compression_code</a>, buf+152, 2); <a name="l00118"></a>00118 memcpy(&h-><a class="code" href="structecat7__mainheader.html#af24db33af652c26d6b2b8a580f3515c4">study_type</a>, buf+154, 12); <a name="l00119"></a>00119 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a3ba61576ab7c202d8d9ba320ea352fe3">patient_id</a>, buf+166, 16); <a name="l00120"></a>00120 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a077ddc2daa174f716708c2077e9ea838">patient_name</a>, buf+182, 32); <a name="l00121"></a>00121 memcpy(&h-><a class="code" href="structecat7__mainheader.html#ae87bc1bf2265ce23abe05b38a4f40fb8">patient_sex</a>, buf+214, 1); <a name="l00122"></a>00122 memcpy(&h-><a class="code" href="structecat7__mainheader.html#acd68a5818bdac5c729a43959b0cddfb7">patient_dexterity</a>, buf+215, 1); <a name="l00123"></a>00123 <span class="keywordflow">if</span>(little) swawbip(buf+216, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ad308b992472cc20d57261aa792159cdd">patient_age</a>, buf+216, 4); <a name="l00124"></a>00124 <span class="keywordflow">if</span>(little) swawbip(buf+220, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ab0f4df121ec2ed3b4513f25030d7a041">patient_height</a>, buf+220, 4); <a name="l00125"></a>00125 <span class="keywordflow">if</span>(little) swawbip(buf+224, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a5b2180feb84eb042e654666703a068dd">patient_weight</a>, buf+224, 4); <a name="l00126"></a>00126 <span class="keywordflow">if</span>(little) swawbip(buf+228, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ab803645ef6896ad52033966c67ef985c">patient_birth_date</a>, buf+228, 4); <a name="l00127"></a>00127 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a14b5bcdff318fa65d2d95126ebc550fc">physician_name</a>, buf+232, 32); <a name="l00128"></a>00128 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a973da5bf47d9376aff843b93e125156a">operator_name</a>, buf+264, 32); <a name="l00129"></a>00129 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a08d61531e0bfce977dd5542f175925d9">study_description</a>, buf+296, 32); <a name="l00130"></a>00130 <span class="keywordflow">if</span>(little) swabip(buf+328, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a8cb7331fc88e1dbd17e0f403f4d7ae8e">acquisition_type</a>, buf+328, 2); <a name="l00131"></a>00131 <span class="keywordflow">if</span>(little) swabip(buf+330, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a1815f9e0d0083df1ba25af49d033d202">patient_orientation</a>, buf+330, 2); <a name="l00132"></a>00132 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a2db57a2b51c8d87a9ab9ff09a811c556">facility_name</a>, buf+332, 20); <a name="l00133"></a>00133 <span class="keywordflow">if</span>(little) swabip(buf+352, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#aacea3a6d9b88add8710823f1c17acc37">num_planes</a>, buf+352, 2); <a name="l00134"></a>00134 <span class="keywordflow">if</span>(little) swabip(buf+354, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a7df83de0bdabd7ebf190e1a58bfc2e20">num_frames</a>, buf+354, 2); <a name="l00135"></a>00135 <span class="keywordflow">if</span>(little) swabip(buf+356, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a45cd434d5473194ffdbb27c99d9d7e6a">num_gates</a>, buf+356, 2); <a name="l00136"></a>00136 <span class="keywordflow">if</span>(little) swabip(buf+358, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#acd90030d2317a4af5583431115b21a4a">num_bed_pos</a>, buf+358, 2); <a name="l00137"></a>00137 <span class="keywordflow">if</span>(little) swawbip(buf+360, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a8a6bbcaf85225c67f8bfa5e564318e79">init_bed_position</a>, buf+360, 4); <a name="l00138"></a>00138 <span class="keywordflow">if</span>(little) swawbip(buf+364, 15*4); memcpy(h-><a class="code" href="structecat7__mainheader.html#af66557399c7d155de452a077228fea71">bed_position</a>, buf+364, 15*4); <a name="l00139"></a>00139 <span class="keywordflow">if</span>(little) swawbip(buf+424, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a786937fcd3117a9d287da4636e3c0f37">plane_separation</a>, buf+424, 4); <a name="l00140"></a>00140 <span class="keywordflow">if</span>(little) swabip(buf+428, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ae649a9fa654b298d7d06d6cfd1648e9c">lwr_sctr_thres</a>, buf+428, 2); <a name="l00141"></a>00141 <span class="keywordflow">if</span>(little) swabip(buf+430, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ae07a9464fb11913d993a7ce9aaac588f">lwr_true_thres</a>, buf+430, 2); <a name="l00142"></a>00142 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a4f229882ef975488024b05405f56d3ba">upr_true_thres</a>, buf+432, 2); <span class="keywordflow">if</span>(little) swabip(&h-><a class="code" href="structecat7__mainheader.html#a4f229882ef975488024b05405f56d3ba">upr_true_thres</a>, 2); <a name="l00143"></a>00143 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a957502542c1a109b52a59af75ebcb27e">user_process_code</a>, buf+434, 10); <a name="l00144"></a>00144 <span class="keywordflow">if</span>(little) swabip(buf+444, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#ab04f1792d0c7f628dd8e05c5a8f38363">acquisition_mode</a>, buf+444, 2); <a name="l00145"></a>00145 <span class="keywordflow">if</span>(little) swawbip(buf+446, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#abc625a3c1c4052d87e502bfca8599b8f">bin_size</a>, buf+446, 4); <a name="l00146"></a>00146 <span class="keywordflow">if</span>(little) swawbip(buf+450, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#aca100a7dd48d4b56c1ec7028716c6a68">branching_fraction</a>, buf+450, 4); <a name="l00147"></a>00147 <span class="keywordflow">if</span>(little) swawbip(buf+454, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a7d6f47d203ff147bd73cebe41167474c">dose_start_time</a>, buf+454, 4); <a name="l00148"></a>00148 <span class="keywordflow">if</span>(little) swawbip(buf+458, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a002205eea64b6f6be8a8c254c8a6b95f">dosage</a>, buf+458, 4); <a name="l00149"></a>00149 <span class="keywordflow">if</span>(little) swawbip(buf+462, 4); memcpy(&h-><a class="code" href="structecat7__mainheader.html#a0d4893dd1647140f775698b5649223d4">well_counter_corr_factor</a>, buf+462, 4); <a name="l00150"></a>00150 memcpy(&h-><a class="code" href="structecat7__mainheader.html#aa6bbc6d3bc1a60f51dd721db2c3be09d">data_units</a>, buf+466, 32); <a name="l00151"></a>00151 <span class="keywordflow">if</span>(little) swabip(buf+498, 2); memcpy(&h-><a class="code" href="structecat7__mainheader.html#abb92c232a653327ce7a089be23f6df09">septa_state</a>, buf+498, 2); <a name="l00152"></a>00152 memcpy(&h-><a class="code" href="structecat7__mainheader.html#a5766625b9e684946431b919aff7426b9">fill_cti</a>, buf+500, 12); <a name="l00153"></a>00153 <a name="l00154"></a>00154 <span class="comment">/* Patient birth date can have been saved in two different int formats,</span> <a name="l00155"></a>00155 <span class="comment"> either YYYYMMDD or as seconds from start of year 1970. In latter case</span> <a name="l00156"></a>00156 <span class="comment"> the number can be negative. */</span> <a name="l00157"></a>00157 <span class="comment">/* Seconds from start of year 1970 are converted to YYYYMMDD format */</span> <a name="l00158"></a>00158 <span class="keywordflow">if</span>(isdate4(h-><a class="code" href="structecat7__mainheader.html#ab803645ef6896ad52033966c67ef985c">patient_birth_date</a>, NULL, NULL, NULL)!=0) { <a name="l00159"></a>00159 time_to_tm((time_t)h-><a class="code" href="structecat7__mainheader.html#ab803645ef6896ad52033966c67ef985c">patient_birth_date</a>, 12*3600-timezone, &st); <a name="l00160"></a>00160 h-><a class="code" href="structecat7__mainheader.html#ab803645ef6896ad52033966c67ef985c">patient_birth_date</a>=10000*(st.tm_year+1900)+100*(st.tm_mon+1)+st.tm_mday; <a name="l00161"></a>00161 } <a name="l00162"></a>00162 <a name="l00163"></a>00163 <span class="keywordflow">return</span>(0); <a name="l00164"></a>00164 } <a name="l00165"></a>00165 <span class="comment">/*****************************************************************************/</span> <a name="l00166"></a>00166 <a name="l00167"></a>00167 <span class="comment">/*****************************************************************************/</span> <a name="l00177"></a><a class="code" href="ecat7_8h.html#ac96bedc17610a5a05d95ce7742d44b5c">00177</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#ac96bedc17610a5a05d95ce7742d44b5c">ecat7ReadImageheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__imageheader.html">ECAT7_imageheader</a> *h) { <a name="l00178"></a>00178 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00179"></a>00179 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00180"></a>00180 <a name="l00181"></a>00181 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadImageheader()\n"</span>); <a name="l00182"></a>00182 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00183"></a>00183 little=little_endian(); <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"little=%d\n"</span>, little); <a name="l00184"></a>00184 <a name="l00185"></a>00185 <span class="comment">/* Seek the subheader block */</span> <a name="l00186"></a>00186 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00187"></a>00187 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00188"></a>00188 <span class="comment">/* Read the header block */</span> <a name="l00189"></a>00189 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00190"></a>00190 <a name="l00191"></a>00191 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00192"></a>00192 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>, buf+0, 2); <a name="l00193"></a>00193 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#acabca0e942a462b769c59804250b9a1c">num_dimensions</a>, buf+2, 2); <a name="l00194"></a>00194 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a37550f6ce6a114de12aa0d038f573230">x_dimension</a>, buf+4, 2); <a name="l00195"></a>00195 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a217d17f8b1f921fad128ce923fa21c35">y_dimension</a>, buf+6, 2); <a name="l00196"></a>00196 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a7d4bc496c4aa5581d41dff3e5b40a231">z_dimension</a>, buf+8, 2); <a name="l00197"></a>00197 <span class="keywordflow">if</span>(little) swawbip(buf+10, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a340a9f9df62d09e5ea8deefd64fec57d">x_offset</a>, buf+10, 4); <a name="l00198"></a>00198 <span class="keywordflow">if</span>(little) swawbip(buf+14, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ad8fcacd511c334e938839c2a7730bbc3">y_offset</a>, buf+14, 4); <a name="l00199"></a>00199 <span class="keywordflow">if</span>(little) swawbip(buf+18, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ab08e016d3c4fb780aeb1048398308254">z_offset</a>, buf+18, 4); <a name="l00200"></a>00200 <span class="keywordflow">if</span>(little) swawbip(buf+22, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a7290c1acb6f8263f85492e3b5c44cdc9">recon_zoom</a>, buf+22, 4); <a name="l00201"></a>00201 <span class="keywordflow">if</span>(little) swawbip(buf+26, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a2e44d8de4b5809864432376f713968ec">scale_factor</a>, buf+26, 4); <a name="l00202"></a>00202 <span class="keywordflow">if</span>(little) swabip(buf+30, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ad599fd1c46c31b598fb2a6f376d562eb">image_min</a>, buf+30, 2); <a name="l00203"></a>00203 <span class="keywordflow">if</span>(little) swabip(buf+32, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a08e876d5e60f96c5b90fc99eeedde10f">image_max</a>, buf+32, 2); <a name="l00204"></a>00204 <span class="keywordflow">if</span>(little) swawbip(buf+34, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a3f2c13786fd5e912c31f481f1b48cff3">x_pixel_size</a>, buf+34, 4); <a name="l00205"></a>00205 <span class="keywordflow">if</span>(little) swawbip(buf+38, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#aefd3bbd592bedba8d68dbec20d855586">y_pixel_size</a>, buf+38, 4); <a name="l00206"></a>00206 <span class="keywordflow">if</span>(little) swawbip(buf+42, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a404d5327064305a7a9dfe75831b29581">z_pixel_size</a>, buf+42, 4); <a name="l00207"></a>00207 <span class="keywordflow">if</span>(little) swawbip(buf+46, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ae75244f1e5532bb257095e7583ef79a3">frame_duration</a>, buf+46, 4); <a name="l00208"></a>00208 <span class="keywordflow">if</span>(little) swawbip(buf+50, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a5fb0b8e2eec4391ed600f0efa4666631">frame_start_time</a>, buf+50, 4); <a name="l00209"></a>00209 <span class="keywordflow">if</span>(little) swabip(buf+54, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ad5bd2f0c18209cb80d17c9ccabd44753">filter_code</a>, buf+54, 2); <a name="l00210"></a>00210 <span class="keywordflow">if</span>(little) swawbip(buf+56, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a94c84cc00c6100c166167eaf0d4d34c9">x_resolution</a>, buf+56, 4); <a name="l00211"></a>00211 <span class="keywordflow">if</span>(little) swawbip(buf+60, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a88478cdda647f217f81ba4fe753f8602">y_resolution</a>, buf+60, 4); <a name="l00212"></a>00212 <span class="keywordflow">if</span>(little) swawbip(buf+64, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#aa0599cfeedd3708b353ef027f281657d">z_resolution</a>, buf+64, 4); <a name="l00213"></a>00213 <span class="keywordflow">if</span>(little) swawbip(buf+68, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a8649bac62fab217d1b49131a6bbeaffa">num_r_elements</a>, buf+68, 4); <a name="l00214"></a>00214 <span class="keywordflow">if</span>(little) swawbip(buf+72, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a046f5ff250412ec423b6ee43b8cf1496">num_angles</a>, buf+72, 4); <a name="l00215"></a>00215 <span class="keywordflow">if</span>(little) swawbip(buf+76, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ac9bdf6594e830a86e04f92e8c5603225">z_rotation_angle</a>, buf+76, 4); <a name="l00216"></a>00216 <span class="keywordflow">if</span>(little) swawbip(buf+80, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a753cec72803431b0b3ecf130b2fee645">decay_corr_fctr</a>, buf+80, 4); <a name="l00217"></a>00217 <span class="keywordflow">if</span>(little) swawbip(buf+84, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a6e617bc2f7f9ddd225af3946198e2d06">processing_code</a>, buf+84, 4); <a name="l00218"></a>00218 <span class="keywordflow">if</span>(little) swawbip(buf+88, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a5f8654cf5e2a615ee18a0b2226ae2b36">gate_duration</a>, buf+88, 4); <a name="l00219"></a>00219 <span class="keywordflow">if</span>(little) swawbip(buf+92, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#af2652b5447df30fa1198d34a2a6c1a81">r_wave_offset</a>, buf+92, 4); <a name="l00220"></a>00220 <span class="keywordflow">if</span>(little) swawbip(buf+96, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a9562593e2089ae493975e4770c455038">num_accepted_beats</a>, buf+96, 4); <a name="l00221"></a>00221 <span class="keywordflow">if</span>(little) swawbip(buf+100, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#afbe5e1b92387d70242a64e0253c8d196">filter_cutoff_frequency</a>, buf+100, 4); <a name="l00222"></a>00222 <span class="keywordflow">if</span>(little) swawbip(buf+104, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a92cbc1f31172ccde80cfbc74763c6bc7">filter_resolution</a>, buf+104, 4); <a name="l00223"></a>00223 <span class="keywordflow">if</span>(little) swawbip(buf+108, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#aab3e31ea98f55be2dc136b6d38c6084c">filter_ramp_slope</a>, buf+108, 4); <a name="l00224"></a>00224 <span class="keywordflow">if</span>(little) swabip(buf+112, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a09d1048e2e74c28630242ffa98270ed0">filter_order</a>, buf+112, 2); <a name="l00225"></a>00225 <span class="keywordflow">if</span>(little) swawbip(buf+114, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ac8b2021a00f2c1acd43213ae780e9993">filter_scatter_fraction</a>, buf+114, 4); <a name="l00226"></a>00226 <span class="keywordflow">if</span>(little) swawbip(buf+118, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a26a5abf25a440aee0f0adfe0eb5a781b">filter_scatter_slope</a>, buf+118, 4); <a name="l00227"></a>00227 memcpy(&h-><a class="code" href="structecat7__imageheader.html#ad5ed4d561b195452ddb325baab3adbb2">annotation</a>, buf+122, 40); <a name="l00228"></a>00228 <span class="keywordflow">if</span>(little) swawbip(buf+162, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a5ab2348c59531dd764f4b659aaf4f75d">mt_1_1</a>, buf+162, 4); <a name="l00229"></a>00229 <span class="keywordflow">if</span>(little) swawbip(buf+166, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a0b84bc9d478661de5916723009b846c6">mt_1_2</a>, buf+166, 4); <a name="l00230"></a>00230 <span class="keywordflow">if</span>(little) swawbip(buf+170, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ac6112b47c4727d97a00927ec2cd66e08">mt_1_3</a>, buf+170, 4); <a name="l00231"></a>00231 <span class="keywordflow">if</span>(little) swawbip(buf+174, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a9aa22eb676327cd80cc852f1da66b2e9">mt_2_1</a>, buf+174, 4); <a name="l00232"></a>00232 <span class="keywordflow">if</span>(little) swawbip(buf+178, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a0387c09ae6ce5366bd9f557666807db1">mt_2_2</a>, buf+178, 4); <a name="l00233"></a>00233 <span class="keywordflow">if</span>(little) swawbip(buf+182, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a9f12cf3a77cbec5884db289e3f671b88">mt_2_3</a>, buf+182, 4); <a name="l00234"></a>00234 <span class="keywordflow">if</span>(little) swawbip(buf+186, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a93f68b4a932241c32442675541b5eafa">mt_3_1</a>, buf+186, 4); <a name="l00235"></a>00235 <span class="keywordflow">if</span>(little) swawbip(buf+190, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a3ff42a2285fff0e821df591aaf089491">mt_3_2</a>, buf+190, 4); <a name="l00236"></a>00236 <span class="keywordflow">if</span>(little) swawbip(buf+194, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#af571054289d88e95552f0a792667b84f">mt_3_3</a>, buf+194, 4); <a name="l00237"></a>00237 <span class="keywordflow">if</span>(little) swawbip(buf+198, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#af1a33983eb2c35e26f2ca2531b092d64">rfilter_cutoff</a>, buf+198, 4); <a name="l00238"></a>00238 <span class="keywordflow">if</span>(little) swawbip(buf+202, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ab3eaa110df4b571b92924dd94826f3c0">rfilter_resolution</a>, buf+202, 4); <a name="l00239"></a>00239 <span class="keywordflow">if</span>(little) swabip(buf+206, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a2ffc2cb4c75c295a5e58af8592ce9a4c">rfilter_code</a>, buf+206, 2); <a name="l00240"></a>00240 <span class="keywordflow">if</span>(little) swabip(buf+208, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a847e512ab018e42b37430b291f2a6e1b">rfilter_order</a>, buf+208, 2); <a name="l00241"></a>00241 <span class="keywordflow">if</span>(little) swawbip(buf+210, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a574df5d81d5b1ab3b5b0ada42c976dbe">zfilter_cutoff</a>, buf+210, 4); <a name="l00242"></a>00242 <span class="keywordflow">if</span>(little) swawbip(buf+214, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#aacfbdbd324b30a3b319eeca5bbeeee12">zfilter_resolution</a>, buf+214, 4); <a name="l00243"></a>00243 <span class="keywordflow">if</span>(little) swabip(buf+218, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ae90d91c961a0ddaca2e1e3bbd77095c7">zfilter_code</a>, buf+218, 2); <a name="l00244"></a>00244 <span class="keywordflow">if</span>(little) swabip(buf+220, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a8f7008e7dbfe7e4a7121c254375f0f47">zfilter_order</a>, buf+220, 2); <a name="l00245"></a>00245 <span class="keywordflow">if</span>(little) swawbip(buf+222, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ae68e680639ac505c93d21ec447481ebf">mt_1_4</a>, buf+222, 4); <a name="l00246"></a>00246 <span class="keywordflow">if</span>(little) swawbip(buf+226, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#adcd2dc3fbd8ded9a1ed5f8d1e85dbda8">mt_2_4</a>, buf+226, 4); <a name="l00247"></a>00247 <span class="keywordflow">if</span>(little) swawbip(buf+230, 4); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ae3481c71f0ab8a1895d37fc7de71ed80">mt_3_4</a>, buf+230, 4); <a name="l00248"></a>00248 <span class="keywordflow">if</span>(little) swabip(buf+234, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a26e3211da4e7bdee598a912f59bbd80c">scatter_type</a>, buf+234, 2); <a name="l00249"></a>00249 <span class="keywordflow">if</span>(little) swabip(buf+236, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#a060576c7efe995f93a17fb73881b6607">recon_type</a>, buf+236, 2); <a name="l00250"></a>00250 <span class="keywordflow">if</span>(little) swabip(buf+238, 2); memcpy(&h-><a class="code" href="structecat7__imageheader.html#ae67f441f35e6b40499b03fca8a7d3d90">recon_views</a>, buf+238, 2); <a name="l00251"></a>00251 memcpy(&h-><a class="code" href="structecat7__imageheader.html#a70ba01c4d16b5c5e5ad7b6a2143aba4c">fill_cti</a>, buf+240, 174); <a name="l00252"></a>00252 memcpy(&h-><a class="code" href="structecat7__imageheader.html#ac7858af3c6aeab6ee47867cba833e0a9">fill_user</a>, buf+414, 96); <a name="l00253"></a>00253 <a name="l00254"></a>00254 <span class="keywordflow">return</span>(0); <a name="l00255"></a>00255 } <a name="l00256"></a>00256 <span class="comment">/*****************************************************************************/</span> <a name="l00257"></a>00257 <a name="l00258"></a>00258 <span class="comment">/*****************************************************************************/</span> <a name="l00268"></a><a class="code" href="ecat7_8h.html#a610647e39e3d9de0f24b1886a064d6f8">00268</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a610647e39e3d9de0f24b1886a064d6f8">ecat7ReadAttenheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__attenheader.html">ECAT7_attenheader</a> *h) { <a name="l00269"></a>00269 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00270"></a>00270 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00271"></a>00271 <a name="l00272"></a>00272 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadAttenheader()\n"</span>); <a name="l00273"></a>00273 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00274"></a>00274 little=little_endian(); <a name="l00275"></a>00275 <a name="l00276"></a>00276 <span class="comment">/* Seek the subheader block */</span> <a name="l00277"></a>00277 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00278"></a>00278 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00279"></a>00279 <span class="comment">/* Read the header block */</span> <a name="l00280"></a>00280 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00281"></a>00281 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00282"></a>00282 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a1e279bd4817d0a7e940b796ebb901d39">data_type</a>, buf+0, 2); <a name="l00283"></a>00283 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a26c4cf72bd2f93808c53113c45e23552">num_dimensions</a>, buf+2, 2); <a name="l00284"></a>00284 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a559039bce5d7432a2c1b2627400055ed">attenuation_type</a>, buf+4, 2); <a name="l00285"></a>00285 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#ab469110a7579b07996d48b7ec6367cc5">num_r_elements</a>, buf+6, 2); <a name="l00286"></a>00286 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a5b3cb8793edf2f160a1c8893b09ffa4e">num_angles</a>, buf+8, 2); <a name="l00287"></a>00287 <span class="keywordflow">if</span>(little) swabip(buf+10, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a78a581032441626011b5394a43bdf2af">num_z_elements</a>, buf+10, 2); <a name="l00288"></a>00288 <span class="keywordflow">if</span>(little) swabip(buf+12, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#ab27975a4dbecf07f924abec9ca3022d2">ring_difference</a>, buf+12, 2); <a name="l00289"></a>00289 <span class="keywordflow">if</span>(little) swawbip(buf+14, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#ad25d43e88b93dc981f8d3dbc85b246a4">x_resolution</a>, buf+14, 4); <a name="l00290"></a>00290 <span class="keywordflow">if</span>(little) swawbip(buf+18, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a9f559f6c3fdde8d0bf93d110a9057b8b">y_resolution</a>, buf+18, 4); <a name="l00291"></a>00291 <span class="keywordflow">if</span>(little) swawbip(buf+22, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a409675e4d6361e84c616079c2316f061">z_resolution</a>, buf+22, 4); <a name="l00292"></a>00292 <span class="keywordflow">if</span>(little) swawbip(buf+26, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a9d90102a2f5032e260edbb5834dff17f">w_resolution</a>, buf+26, 4); <a name="l00293"></a>00293 <span class="keywordflow">if</span>(little) swawbip(buf+30, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a6c90e753da243d94a24384873b0195b3">scale_factor</a>, buf+30, 4); <a name="l00294"></a>00294 <span class="keywordflow">if</span>(little) swawbip(buf+34, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#adb7c3a7f9ad331e2e2d66ad397d1061e">x_offset</a>, buf+34, 4); <a name="l00295"></a>00295 <span class="keywordflow">if</span>(little) swawbip(buf+38, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#af00ef016e4bec39cccb0566e0d1e1beb">y_offset</a>, buf+38, 4); <a name="l00296"></a>00296 <span class="keywordflow">if</span>(little) swawbip(buf+42, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#aade44803290eb7c042e8db22b647c3f3">x_radius</a>, buf+42, 4); <a name="l00297"></a>00297 <span class="keywordflow">if</span>(little) swawbip(buf+46, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a885d8e8ab2ce770cf193b4273d591de0">y_radius</a>, buf+46, 4); <a name="l00298"></a>00298 <span class="keywordflow">if</span>(little) swawbip(buf+50, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a8ee01a9d388f1116b71395e62618ced1">tilt_angle</a>, buf+50, 4); <a name="l00299"></a>00299 <span class="keywordflow">if</span>(little) swawbip(buf+54, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#acfdb1bdf8255158bf0ce77e6b2d682e3">attenuation_coeff</a>, buf+54, 4); <a name="l00300"></a>00300 <span class="keywordflow">if</span>(little) swawbip(buf+58, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a6734bcf692194b4b2fe3477574549b04">attenuation_min</a>, buf+58, 4); <a name="l00301"></a>00301 <span class="keywordflow">if</span>(little) swawbip(buf+62, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#abc998fff1d231a8c23ab7bd75e980289">attenuation_max</a>, buf+62, 4); <a name="l00302"></a>00302 <span class="keywordflow">if</span>(little) swawbip(buf+66, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#ad5b24b2b096cbddc14a358ba21be77d2">skull_thickness</a>, buf+66, 4); <a name="l00303"></a>00303 <span class="keywordflow">if</span>(little) swabip(buf+70, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#ac4a70cc92f5cd58fcd9c4d10757c04a0">num_additional_atten_coeff</a>, buf+70, 2); <a name="l00304"></a>00304 <span class="keywordflow">if</span>(little) swawbip(buf+72, 8*4); memcpy(h-><a class="code" href="structecat7__attenheader.html#ab71ff5127a01d22c4c376a82a128f926">additional_atten_coeff</a>, buf+72, 8*4); <a name="l00305"></a>00305 <span class="keywordflow">if</span>(little) swawbip(buf+104, 4); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a2ec9a53c0910ba896e41f486b479fd7c">edge_finding_threshold</a>, buf+104, 4); <a name="l00306"></a>00306 <span class="keywordflow">if</span>(little) swabip(buf+108, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#addc2063b495772262146f9385c294323">storage_order</a>, buf+108, 2); <a name="l00307"></a>00307 <span class="keywordflow">if</span>(little) swabip(buf+110, 2); memcpy(&h-><a class="code" href="structecat7__attenheader.html#a45d73fa0fd340f998265aaf15a123ba8">span</a>, buf+110, 2); <a name="l00308"></a>00308 <span class="keywordflow">if</span>(little) swabip(buf+112, 64*2); memcpy(h-><a class="code" href="structecat7__attenheader.html#aeb743233aff66821f8cfbf3ec1f64d3b">z_elements</a>, buf+112, 64*2); <a name="l00309"></a>00309 <span class="keywordflow">if</span>(little) swabip(buf+240, 86*2); memcpy(h-><a class="code" href="structecat7__attenheader.html#a81a71d398c08000665962d87c9c62340">fill_cti</a>, buf+240, 86*2); <a name="l00310"></a>00310 <span class="keywordflow">if</span>(little) swabip(buf+412, 50*2); memcpy(h-><a class="code" href="structecat7__attenheader.html#abb0256f4ba630ed598d257a8a797d449">fill_user</a>, buf+412, 50*2); <a name="l00311"></a>00311 <span class="keywordflow">return</span>(0); <a name="l00312"></a>00312 } <a name="l00313"></a>00313 <span class="comment">/*****************************************************************************/</span> <a name="l00314"></a>00314 <a name="l00315"></a>00315 <span class="comment">/*****************************************************************************/</span> <a name="l00325"></a><a class="code" href="ecat7_8h.html#a498d383ca08a832895af8ead01e8c756">00325</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a498d383ca08a832895af8ead01e8c756">ecat7ReadPolmapheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__polmapheader.html">ECAT7_polmapheader</a> *h) { <a name="l00326"></a>00326 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00327"></a>00327 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00328"></a>00328 <a name="l00329"></a>00329 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadPolarmapheader()\n"</span>); <a name="l00330"></a>00330 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00331"></a>00331 little=little_endian(); <a name="l00332"></a>00332 <a name="l00333"></a>00333 <span class="comment">/* Seek the subheader block */</span> <a name="l00334"></a>00334 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00335"></a>00335 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00336"></a>00336 <span class="comment">/* Read the header block */</span> <a name="l00337"></a>00337 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00338"></a>00338 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00339"></a>00339 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>, buf+0, 2); <a name="l00340"></a>00340 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a1c42cbb82e3d3abdea6eab7893e3151f">polar_map_type</a>, buf+2, 2); <a name="l00341"></a>00341 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a503caad4ac64580ff3508463ae8e2342">num_rings</a>, buf+4, 2); <a name="l00342"></a>00342 <span class="keywordflow">if</span>(little) swabip(buf+6, 32*2); memcpy(h-><a class="code" href="structecat7__polmapheader.html#a7888f1b2bc8adcc5bcbf918bbbb82dbb">sectors_per_ring</a>, buf+6, 32*2); <a name="l00343"></a>00343 <span class="keywordflow">if</span>(little) swawbip(buf+70, 32*4); memcpy(h-><a class="code" href="structecat7__polmapheader.html#a1f8b29deb85540e9d4d30d9a8247c0d6">ring_position</a>, buf+70, 32*4); <a name="l00344"></a>00344 <span class="keywordflow">if</span>(little) swabip(buf+198, 32*2); memcpy(h-><a class="code" href="structecat7__polmapheader.html#af91ac7e401a999aaaaf275c8671ad0ab">ring_angle</a>, buf+198, 32*2); <a name="l00345"></a>00345 <span class="keywordflow">if</span>(little) swabip(buf+262, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#ad59dd3b498ddd78aa65a6cf4425d259a">start_angle</a>, buf+262, 2); <a name="l00346"></a>00346 <span class="keywordflow">if</span>(little) swabip(buf+264, 3*2); memcpy(h-><a class="code" href="structecat7__polmapheader.html#aabc620f253f9b5bc777f1848d8300e5e">long_axis_left</a>, buf+264, 3*2); <a name="l00347"></a>00347 <span class="keywordflow">if</span>(little) swabip(buf+270, 3*2); memcpy(h-><a class="code" href="structecat7__polmapheader.html#af9716ddc7fc52ecbf3980dfd38b4fbd3">long_axis_right</a>, buf+270, 3*2); <a name="l00348"></a>00348 <span class="keywordflow">if</span>(little) swabip(buf+276, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#af4ebf16c597dd432ffb78c267137c40c">position_data</a>, buf+276, 2); <a name="l00349"></a>00349 <span class="keywordflow">if</span>(little) swabip(buf+278, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#ac955d5b0a42d9a35011f0a9687db1ee7">image_min</a>, buf+278, 2); <a name="l00350"></a>00350 <span class="keywordflow">if</span>(little) swabip(buf+280, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a402864bf8af8ab091faf877d831cb7ef">image_max</a>, buf+280, 2); <a name="l00351"></a>00351 <span class="keywordflow">if</span>(little) swawbip(buf+282, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#ad2f57edc57a8a6f328e19ab35502eb20">scale_factor</a>, buf+282, 4); <a name="l00352"></a>00352 <span class="keywordflow">if</span>(little) swawbip(buf+286, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a4a7a70570769b6b97bf67f6358a61f45">pixel_size</a>, buf+286, 4); <a name="l00353"></a>00353 <span class="keywordflow">if</span>(little) swawbip(buf+290, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#af2c0109195bbcc833521aed329d2bd83">frame_duration</a>, buf+290, 4); <a name="l00354"></a>00354 <span class="keywordflow">if</span>(little) swawbip(buf+294, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#ac7757667731f12beb6c6e8a584dd6682">frame_start_time</a>, buf+294, 4); <a name="l00355"></a>00355 <span class="keywordflow">if</span>(little) swabip(buf+298, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a1b6921694ee32fdda5e5c4b1a64f733b">processing_code</a>, buf+298, 2); <a name="l00356"></a>00356 <span class="keywordflow">if</span>(little) swabip(buf+300, 2); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#aa4b2dda11a8f10a58a77307ffbf28805">quant_units</a>, buf+300, 2); <a name="l00357"></a>00357 memcpy(h-><a class="code" href="structecat7__polmapheader.html#a873c092e019a8919b5b82030ed29f6a5">annotation</a>, buf+302, 40); <a name="l00358"></a>00358 <span class="keywordflow">if</span>(little) swawbip(buf+342, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a87bdb3bf6e644c357290c320166e068f">gate_duration</a>, buf+342, 4); <a name="l00359"></a>00359 <span class="keywordflow">if</span>(little) swawbip(buf+346, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#a8a17116a1e9c4a1df60eb9f709b5d03e">r_wave_offset</a>, buf+346, 4); <a name="l00360"></a>00360 <span class="keywordflow">if</span>(little) swawbip(buf+350, 4); memcpy(&h-><a class="code" href="structecat7__polmapheader.html#aae8015e83f96723de1bffab0cf44b2b3">num_accepted_beats</a>, buf+350, 4); <a name="l00361"></a>00361 memcpy(h-><a class="code" href="structecat7__polmapheader.html#ac6a9017cce528db97e827de264aa3142">polar_map_protocol</a>, buf+354, 20); <a name="l00362"></a>00362 memcpy(h-><a class="code" href="structecat7__polmapheader.html#a65ab9e055aedc38c33d62fb6ebf05589">database_name</a>, buf+374, 30); <a name="l00363"></a>00363 <span class="keywordflow">if</span>(little) swabip(buf+404, 27*2); memcpy(h-><a class="code" href="structecat7__polmapheader.html#a164decb72ee26397bf7cba8cd6219b9f">fill_cti</a>, buf+404, 27*2); <a name="l00364"></a>00364 <span class="keywordflow">return</span>(0); <a name="l00365"></a>00365 } <a name="l00366"></a>00366 <span class="comment">/*****************************************************************************/</span> <a name="l00367"></a>00367 <a name="l00368"></a>00368 <span class="comment">/*****************************************************************************/</span> <a name="l00378"></a><a class="code" href="ecat7_8h.html#a53935e3d6fa02ec2414895e205593db3">00378</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a53935e3d6fa02ec2414895e205593db3">ecat7ReadNormheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__normheader.html">ECAT7_normheader</a> *h) { <a name="l00379"></a>00379 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00380"></a>00380 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00381"></a>00381 <a name="l00382"></a>00382 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadNormheader()\n"</span>); <a name="l00383"></a>00383 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00384"></a>00384 little=little_endian(); <a name="l00385"></a>00385 <a name="l00386"></a>00386 <span class="comment">/* Seek the subheader block */</span> <a name="l00387"></a>00387 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00388"></a>00388 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00389"></a>00389 <span class="comment">/* Read the header block */</span> <a name="l00390"></a>00390 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00391"></a>00391 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00392"></a>00392 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#aed9b4cb1c561d080a8336baf7c80fdcc">data_type</a>, buf+0, 2); <a name="l00393"></a>00393 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#ad44f4013947a523f61f406ab26f2e44a">num_r_elements</a>, buf+2, 2); <a name="l00394"></a>00394 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#aa8d1e9e258ad69f1d1af98d30c7a6cf5">num_transaxial_crystals</a>, buf+4, 2); <a name="l00395"></a>00395 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a3cdec5a8adfe16d0d8fcebcfcdbc4d48">num_crystal_rings</a>, buf+6, 2); <a name="l00396"></a>00396 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#aaf63afdf4ccafc47626d7075dbb7727c">crystals_per_ring</a>, buf+8, 2); <a name="l00397"></a>00397 <span class="keywordflow">if</span>(little) swabip(buf+10, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a7ae2012feb0d317925f54c719c6023bb">num_geo_corr_planes</a>, buf+10, 2); <a name="l00398"></a>00398 <span class="keywordflow">if</span>(little) swabip(buf+12, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#aaa65ed681794afc0f8c67b65f05febfb">uld</a>, buf+12, 2); <a name="l00399"></a>00399 <span class="keywordflow">if</span>(little) swabip(buf+14, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a64e2d29a3fd63691aa5fd129a32c32f5">lld</a>, buf+14, 2); <a name="l00400"></a>00400 <span class="keywordflow">if</span>(little) swabip(buf+16, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a7835b4239c66f51f97b0d96c45ce862e">scatter_energy</a>, buf+16, 2); <a name="l00401"></a>00401 <span class="keywordflow">if</span>(little) swawbip(buf+18, 4); memcpy(&h-><a class="code" href="structecat7__normheader.html#a6cb5155206d5f40b4dc6363482364f5a">norm_quality_factor</a>, buf+18, 4); <a name="l00402"></a>00402 <span class="keywordflow">if</span>(little) swabip(buf+22, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a71f98c8e0270a3b738a8228a130e6e4a">norm_quality_factor_code</a>, buf+22, 2); <a name="l00403"></a>00403 <span class="keywordflow">if</span>(little) swawbip(buf+24, 32*4); memcpy(h-><a class="code" href="structecat7__normheader.html#a981b1ecafed3fc8ffb6027ffa4e37d0d">ring_dtcor1</a>, buf+24, 32*4); <a name="l00404"></a>00404 <span class="keywordflow">if</span>(little) swawbip(buf+152, 32*4); memcpy(h-><a class="code" href="structecat7__normheader.html#aca206f1d98356df2cd4ab06bb4b25de0">ring_dtcor2</a>, buf+152, 32*4); <a name="l00405"></a>00405 <span class="keywordflow">if</span>(little) swawbip(buf+280, 8*4); memcpy(h-><a class="code" href="structecat7__normheader.html#acf7f342012459796cfea425e7e9faf7f">crystal_dtcor</a>, buf+280, 8*4); <a name="l00406"></a>00406 <span class="keywordflow">if</span>(little) swabip(buf+312, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a83cc022237d464480c1696fece19729f">span</a>, buf+312, 2); <a name="l00407"></a>00407 <span class="keywordflow">if</span>(little) swabip(buf+314, 2); memcpy(&h-><a class="code" href="structecat7__normheader.html#a7dd4601b4e06fb51df9e6a09d3287154">max_ring_diff</a>, buf+314, 2); <a name="l00408"></a>00408 <span class="keywordflow">if</span>(little) swabip(buf+316, 48*2); memcpy(h-><a class="code" href="structecat7__normheader.html#ad3967ba7c86ba41189c2c6becd34da45">fill_cti</a>, buf+316, 48*2); <a name="l00409"></a>00409 <span class="keywordflow">if</span>(little) swabip(buf+412, 50*2); memcpy(h-><a class="code" href="structecat7__normheader.html#ac9cdc23bf6850c72fca0d38d2e32b640">fill_user</a>, buf+412, 50*2); <a name="l00410"></a>00410 <span class="keywordflow">return</span>(0); <a name="l00411"></a>00411 } <a name="l00412"></a>00412 <span class="comment">/*****************************************************************************/</span> <a name="l00413"></a>00413 <a name="l00414"></a>00414 <span class="comment">/*****************************************************************************/</span> <a name="l00424"></a><a class="code" href="ecat7_8h.html#a37414b52ef395686b6a81d6229d13143">00424</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a37414b52ef395686b6a81d6229d13143">ecat7ReadScanheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__scanheader.html">ECAT7_scanheader</a> *h) { <a name="l00425"></a>00425 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[2*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00426"></a>00426 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00427"></a>00427 <a name="l00428"></a>00428 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadScanheader()\n"</span>); <a name="l00429"></a>00429 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00430"></a>00430 little=little_endian(); <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"little=%d\n"</span>, little); <a name="l00431"></a>00431 <a name="l00432"></a>00432 <span class="comment">/* Seek the subheader block */</span> <a name="l00433"></a>00433 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00434"></a>00434 <span class="comment">/* Read the header block */</span> <a name="l00435"></a>00435 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 2, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00436"></a>00436 <a name="l00437"></a>00437 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00438"></a>00438 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>, buf+0, 2); <a name="l00439"></a>00439 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a4b4361b122c74cfd1c7171357615a2a4">num_dimensions</a>, buf+2, 2); <a name="l00440"></a>00440 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#aee3fd0c2360973b3d53ad1adf6225d7d">num_r_elements</a>, buf+4, 2); <a name="l00441"></a>00441 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a4f58fe54c66427d34390f6530a12f8ed">num_angles</a>, buf+6, 2); <a name="l00442"></a>00442 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a18d65465164971f7dcdfe5b251b6c244">corrections_applied</a>, buf+8, 2); <a name="l00443"></a>00443 <span class="keywordflow">if</span>(little) swabip(buf+10, 64*2); memcpy(h-><a class="code" href="structecat7__scanheader.html#a714bae4e56b451926ceeb0de19222ad7">num_z_elements</a>, buf+10, 64*2); <a name="l00444"></a>00444 <span class="keywordflow">if</span>(little) swabip(buf+138, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ae0d7a41a8f2790fe97d7784367f8daba">ring_difference</a>, buf+138, 2); <a name="l00445"></a>00445 <span class="keywordflow">if</span>(little) swabip(buf+140, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#af0fcf4af875d58acf97d1b412d0299b5">storage_order</a>, buf+140, 2); <a name="l00446"></a>00446 <span class="keywordflow">if</span>(little) swabip(buf+142, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ad0c5b2146bf087faa53cafd30be60aad">axial_compression</a>, buf+142, 2); <a name="l00447"></a>00447 <span class="keywordflow">if</span>(little) swawbip(buf+144, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a3b20e770bbcc8bc038190aab1d777199">x_resolution</a>, buf+144, 4); <a name="l00448"></a>00448 <span class="keywordflow">if</span>(little) swawbip(buf+148, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a377740afa37d1774a1938a37354b5bde">v_resolution</a>, buf+148, 4); <a name="l00449"></a>00449 <span class="keywordflow">if</span>(little) swawbip(buf+152, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a1e5bfef1c77d887f15aa6daada8464d7">z_resolution</a>, buf+152, 4); <a name="l00450"></a>00450 <span class="keywordflow">if</span>(little) swawbip(buf+156, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a3377bb6891117ab290ea776444caecf2">w_resolution</a>, buf+156, 4); <a name="l00451"></a>00451 <span class="keywordflow">if</span>(little) swabip(buf+160, 6*2); memcpy(h-><a class="code" href="structecat7__scanheader.html#add4ebedc9711b45cf37910b50cbfbaf7">fill_gate</a>, buf+160, 6*2); <a name="l00452"></a>00452 <span class="keywordflow">if</span>(little) swawbip(buf+172, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a603f407610a5e2864bfda5c80152f9c9">gate_duration</a>, buf+172, 4); <a name="l00453"></a>00453 <span class="keywordflow">if</span>(little) swawbip(buf+176, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ae2bdd51e3b0c9b838e77793eb71280d7">r_wave_offset</a>, buf+176, 4); <a name="l00454"></a>00454 <span class="keywordflow">if</span>(little) swawbip(buf+180, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ae80300a155f4fcf65cb0851ff8bfa122">num_accepted_beats</a>, buf+180, 4); <a name="l00455"></a>00455 <span class="keywordflow">if</span>(little) swawbip(buf+184, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a77ed7b5edd51bf8feeea367ec4d143bb">scale_factor</a>, buf+184, 4); <a name="l00456"></a>00456 <span class="keywordflow">if</span>(little) swabip(buf+188, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ae8d5fb305463d8a019b452068176c20c">scan_min</a>, buf+188, 2); <a name="l00457"></a>00457 <span class="keywordflow">if</span>(little) swabip(buf+190, 2); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a6cad394ec364530cb33b335f1256fec1">scan_max</a>, buf+190, 2); <a name="l00458"></a>00458 <span class="keywordflow">if</span>(little) swawbip(buf+192, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a421bb8924e8aed2ef38ecc673b71a574">prompts</a>, buf+192, 4); <a name="l00459"></a>00459 <span class="keywordflow">if</span>(little) swawbip(buf+196, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a8f67a369e6ef599551b8271c08a8b5d7">delayed</a>, buf+196, 4); <a name="l00460"></a>00460 <span class="keywordflow">if</span>(little) swawbip(buf+200, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a7d8ff4b790f38d9e8be6bf83c4bd185d">multiples</a>, buf+200, 4); <a name="l00461"></a>00461 <span class="keywordflow">if</span>(little) swawbip(buf+204, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a2e6e0eaff39c88e72fb766c85ebf5c07">net_trues</a>, buf+204, 4); <a name="l00462"></a>00462 <span class="keywordflow">if</span>(little) swawbip(buf+208, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#aaa779c47e6e11a47795aee0c90437adf">tot_avg_cor</a>, buf+208, 4); <a name="l00463"></a>00463 <span class="keywordflow">if</span>(little) swawbip(buf+212, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a80ffd347fbe13b90be3264318d2cb754">tot_avg_uncor</a>, buf+212, 4); <a name="l00464"></a>00464 <span class="keywordflow">if</span>(little) swawbip(buf+216, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#afcd90f2a17b7af7f495e7ffd5a0a52a8">total_coin_rate</a>, buf+216, 4); <a name="l00465"></a>00465 <span class="keywordflow">if</span>(little) swawbip(buf+220, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a73b5c4aae33a122d91703ce6a3330f89">frame_start_time</a>, buf+220, 4); <a name="l00466"></a>00466 <span class="keywordflow">if</span>(little) swawbip(buf+224, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#ae7eafb4bf000a8404407489dbd6b54ac">frame_duration</a>, buf+224, 4); <a name="l00467"></a>00467 <span class="keywordflow">if</span>(little) swawbip(buf+228, 4); memcpy(&h-><a class="code" href="structecat7__scanheader.html#a97b2c1a33ab3c42a21aa9f1d3d9e21e7">deadtime_correction_factor</a>, buf+228, 4); <a name="l00468"></a>00468 <span class="keywordflow">if</span>(little) swabip(buf+232, 90*2); memcpy(h-><a class="code" href="structecat7__scanheader.html#a47f557ddadc32ba65cb7cde63dd959d0">fill_cti</a>, buf+232, 90*2); <a name="l00469"></a>00469 <span class="keywordflow">if</span>(little) swabip(buf+412, 50*2); memcpy(h-><a class="code" href="structecat7__scanheader.html#a87d711f5f18ee88a7ba4f67f65574cee">fill_user</a>, buf+412, 50*2); <a name="l00470"></a>00470 <span class="keywordflow">if</span>(little) swawbip(buf+512, 128*4); memcpy(h-><a class="code" href="structecat7__scanheader.html#a30c149d12f3ee869d8303c8430fca0f0">uncor_singles</a>, buf+512, 128*4); <a name="l00471"></a>00471 <span class="keywordflow">return</span>(0); <a name="l00472"></a>00472 } <a name="l00473"></a>00473 <span class="comment">/*****************************************************************************/</span> <a name="l00474"></a>00474 <a name="l00475"></a>00475 <span class="comment">/*****************************************************************************/</span> <a name="l00485"></a><a class="code" href="ecat7_8h.html#aa37d19fc4f30868ea0b062d39d9d8450">00485</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#aa37d19fc4f30868ea0b062d39d9d8450">ecat7Read2DScanheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__2Dscanheader.html">ECAT7_2Dscanheader</a> *h) { <a name="l00486"></a>00486 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00487"></a>00487 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00488"></a>00488 <a name="l00489"></a>00489 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7Read2DScanheader()\n"</span>); <a name="l00490"></a>00490 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00491"></a>00491 little=little_endian(); <a name="l00492"></a>00492 <a name="l00493"></a>00493 <span class="comment">/* Seek the subheader block */</span> <a name="l00494"></a>00494 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00495"></a>00495 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00496"></a>00496 <span class="comment">/* Read the header block */</span> <a name="l00497"></a>00497 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00498"></a>00498 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00499"></a>00499 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>, buf+0, 2); <a name="l00500"></a>00500 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a592447ad964ba95e2a29a7d3d2284576">num_dimensions</a>, buf+2, 2); <a name="l00501"></a>00501 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#adfe07cff44e3145ca65c234079ff255a">num_r_elements</a>, buf+4, 2); <a name="l00502"></a>00502 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ad0d1958b231f10aa68d60c653ec30627">num_angles</a>, buf+6, 2); <a name="l00503"></a>00503 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a8c126cf20e0b59765720d52fbfbf92d9">corrections_applied</a>, buf+8, 2); <a name="l00504"></a>00504 <span class="keywordflow">if</span>(little) swabip(buf+10, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a697bb7a95d88bce42bf80e2d6bd69037">num_z_elements</a>, buf+10, 2); <a name="l00505"></a>00505 <span class="keywordflow">if</span>(little) swabip(buf+12, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a6a37718cdfb000b2724062232ecadcf6">ring_difference</a>, buf+12, 2); <a name="l00506"></a>00506 <span class="keywordflow">if</span>(little) swawbip(buf+14, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a77e48781a4d1eda3cd828fdc2c3b1a62">x_resolution</a>, buf+14, 4); <a name="l00507"></a>00507 <span class="keywordflow">if</span>(little) swawbip(buf+18, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ab4d6c7e3cdd71cf5edb222956600083b">y_resolution</a>, buf+18, 4); <a name="l00508"></a>00508 <span class="keywordflow">if</span>(little) swawbip(buf+22, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#aed6ef090a59741c23fc97dec8eb6e17a">z_resolution</a>, buf+22, 4); <a name="l00509"></a>00509 <span class="keywordflow">if</span>(little) swawbip(buf+26, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a946ef333c9627293f30db28447a09f55">w_resolution</a>, buf+26, 4); <a name="l00510"></a>00510 <span class="keywordflow">if</span>(little) swabip(buf+30, 6*2); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#a4a79a87fe7768f023c873b9d38bb75d5">fill_gate</a>, buf+30, 6*2); <a name="l00511"></a>00511 <span class="keywordflow">if</span>(little) swawbip(buf+42, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a790afb788d414578a5c3ab75ee33e6ef">gate_duration</a>, buf+42, 4); <a name="l00512"></a>00512 <span class="keywordflow">if</span>(little) swawbip(buf+46, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ac6159922811bbc4faa5227977aa1aecc">r_wave_offset</a>, buf+46, 4); <a name="l00513"></a>00513 <span class="keywordflow">if</span>(little) swawbip(buf+50, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a268204e52f918f4738cbd4b6b41c0010">num_accepted_beats</a>, buf+50, 4); <a name="l00514"></a>00514 <span class="keywordflow">if</span>(little) swawbip(buf+54, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ae27dc6ec3582a4e2a0e69ea9377993a3">scale_factor</a>, buf+54, 4); <a name="l00515"></a>00515 <span class="keywordflow">if</span>(little) swabip(buf+58, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a6a644bf6cbb45523522f20ebf730af2c">scan_min</a>, buf+58, 2); <a name="l00516"></a>00516 <span class="keywordflow">if</span>(little) swabip(buf+60, 2); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#abaf0931401015351496475447bdf2e4a">scan_max</a>, buf+60, 2); <a name="l00517"></a>00517 <span class="keywordflow">if</span>(little) swawbip(buf+62, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ac4abf37f9a786b4133fe6526cf29463c">prompts</a>, buf+62, 4); <a name="l00518"></a>00518 <span class="keywordflow">if</span>(little) swawbip(buf+66, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a01d95cfc32cb38a15067f6f2b6835176">delayed</a>, buf+66, 4); <a name="l00519"></a>00519 <span class="keywordflow">if</span>(little) swawbip(buf+70, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a875117198383dd1ede79e8e0c7e9e98d">multiples</a>, buf+70, 4); <a name="l00520"></a>00520 <span class="keywordflow">if</span>(little) swawbip(buf+74, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#af370e602f761d19cc1663c84f8fb0fae">net_trues</a>, buf+74, 4); <a name="l00521"></a>00521 <span class="keywordflow">if</span>(little) swawbip(buf+78, 16*4); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#af2e871980537eb5bf534ff9ed4d7346f">cor_singles</a>, buf+78, 16*4); <a name="l00522"></a>00522 <span class="keywordflow">if</span>(little) swawbip(buf+142, 16*4); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#ae09638e910e990877f98f8f5bf0450fe">uncor_singles</a>, buf+142, 16*4); <a name="l00523"></a>00523 <span class="keywordflow">if</span>(little) swawbip(buf+206, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#ab18f6c6f465371a84b37b909e3294dda">tot_avg_cor</a>, buf+206, 4); <a name="l00524"></a>00524 <span class="keywordflow">if</span>(little) swawbip(buf+210, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a3e4245e6136c30d208c36e2e60432ea1">tot_avg_uncor</a>, buf+210, 4); <a name="l00525"></a>00525 <span class="keywordflow">if</span>(little) swawbip(buf+214, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a2f154fd294d938039e15dbc452aae7a0">total_coin_rate</a>, buf+214, 4); <a name="l00526"></a>00526 <span class="keywordflow">if</span>(little) swawbip(buf+218, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a8930d980328907f75c2518b36e67b36f">frame_start_time</a>, buf+218, 4); <a name="l00527"></a>00527 <span class="keywordflow">if</span>(little) swawbip(buf+222, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#afe892b507d85054ad72ce5ac289b3b05">frame_duration</a>, buf+222, 4); <a name="l00528"></a>00528 <span class="keywordflow">if</span>(little) swawbip(buf+226, 4); memcpy(&h-><a class="code" href="structecat7__2Dscanheader.html#a8788db06847593186cbeaef45e03b838">deadtime_correction_factor</a>, buf+226, 4); <a name="l00529"></a>00529 <span class="keywordflow">if</span>(little) swabip(buf+230, 8*2); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#a32d5dfdc807697d412704c5ac53adced">physical_planes</a>, buf+230, 8*2); <a name="l00530"></a>00530 <span class="keywordflow">if</span>(little) swabip(buf+246, 83*2); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#ad3b13c4044cd9de84db192d700ebdca4">fill_cti</a>, buf+246, 83*2); <a name="l00531"></a>00531 <span class="keywordflow">if</span>(little) swabip(buf+412, 50*2); memcpy(h-><a class="code" href="structecat7__2Dscanheader.html#a07da82cdd5f206c0c52eeb20576f2a77">fill_user</a>, buf+412, 50*2); <a name="l00532"></a>00532 <span class="keywordflow">return</span>(0); <a name="l00533"></a>00533 } <a name="l00534"></a>00534 <span class="comment">/*****************************************************************************/</span> <a name="l00535"></a>00535 <a name="l00536"></a>00536 <span class="comment">/*****************************************************************************/</span> <a name="l00546"></a><a class="code" href="ecat7_8h.html#ad9fd9a7a7dc9cabfddf398c32b01bdce">00546</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#ad9fd9a7a7dc9cabfddf398c32b01bdce">ecat7Read2DNormheader</a>(FILE *fp, <span class="keywordtype">int</span> blk, <a class="code" href="structecat7__2Dnormheader.html">ECAT7_2Dnormheader</a> *h) { <a name="l00547"></a>00547 <span class="keywordtype">unsigned</span> <span class="keywordtype">char</span> buf[<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>]; <a name="l00548"></a>00548 <span class="keywordtype">int</span> little; <span class="comment">/* 1 if current platform is little endian (i386), else 0 */</span> <a name="l00549"></a>00549 <a name="l00550"></a>00550 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7Read2Dnormheader()\n"</span>); <a name="l00551"></a>00551 <span class="keywordflow">if</span>(fp==NULL || h==NULL) <span class="keywordflow">return</span>(1); <a name="l00552"></a>00552 little=little_endian(); <a name="l00553"></a>00553 <a name="l00554"></a>00554 <span class="comment">/* Seek the subheader block */</span> <a name="l00555"></a>00555 fseek(fp, (blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00556"></a>00556 <span class="keywordflow">if</span>(ftell(fp)!=(blk-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(2); <a name="l00557"></a>00557 <span class="comment">/* Read the header block */</span> <a name="l00558"></a>00558 <span class="keywordflow">if</span>(fread(buf, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, 1, fp)<1) <span class="keywordflow">return</span>(3); <a name="l00559"></a>00559 <span class="comment">/* Copy the header fields and swap if necessary */</span> <a name="l00560"></a>00560 <span class="keywordflow">if</span>(little) swabip(buf+0, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#aed37b336ee5cf6446aa8f14a8a8589e9">data_type</a>, buf+0, 2); <a name="l00561"></a>00561 <span class="keywordflow">if</span>(little) swabip(buf+2, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a2513b09534273845652deb6b9ab94e25">num_dimensions</a>, buf+2, 2); <a name="l00562"></a>00562 <span class="keywordflow">if</span>(little) swabip(buf+4, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a9445d293a44e6b33f1114ed77a8aae28">num_r_elements</a>, buf+4, 2); <a name="l00563"></a>00563 <span class="keywordflow">if</span>(little) swabip(buf+6, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a56ea7bd534112992fbd8380d9358bb67">num_angles</a>, buf+6, 2); <a name="l00564"></a>00564 <span class="keywordflow">if</span>(little) swabip(buf+8, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a83e70d1ac4fa6dd8421c1302370529ed">num_z_elements</a>, buf+8, 2); <a name="l00565"></a>00565 <span class="keywordflow">if</span>(little) swabip(buf+10, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#ac43014ca13110db8929ecfa9ae0141b0">ring_difference</a>, buf+10, 2); <a name="l00566"></a>00566 <span class="keywordflow">if</span>(little) swawbip(buf+12, 4); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a41e0cea1c8ef72712a6f44c1634fa5fc">scale_factor</a>, buf+12, 4); <a name="l00567"></a>00567 <span class="keywordflow">if</span>(little) swawbip(buf+16, 4); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a27bc628ff44c68378fa60074c442acf1">norm_min</a>, buf+16, 4); <a name="l00568"></a>00568 <span class="keywordflow">if</span>(little) swawbip(buf+20, 4); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#ad8025a3f59f73ba95c2412bd0a72ff9c">norm_max</a>, buf+20, 4); <a name="l00569"></a>00569 <span class="keywordflow">if</span>(little) swawbip(buf+24, 4); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a7bbbf33025bb7abbe325825cb7ede630">fov_source_width</a>, buf+24, 4); <a name="l00570"></a>00570 <span class="keywordflow">if</span>(little) swawbip(buf+28, 4); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a52b683dd331c8c3d3f9f12d354e83522">norm_quality_factor</a>, buf+28, 4); <a name="l00571"></a>00571 <span class="keywordflow">if</span>(little) swabip(buf+32, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a8642a53226b65eaf30bd7205432871a0">norm_quality_factor_code</a>, buf+32, 2); <a name="l00572"></a>00572 <span class="keywordflow">if</span>(little) swabip(buf+34, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a0a45487c98abca585bb1865c00d8bbfb">storage_order</a>, buf+34, 2); <a name="l00573"></a>00573 <span class="keywordflow">if</span>(little) swabip(buf+36, 2); memcpy(&h-><a class="code" href="structecat7__2Dnormheader.html#a863a208786a3728549c8f8e3f1aa7082">span</a>, buf+36, 2); <a name="l00574"></a>00574 <span class="keywordflow">if</span>(little) swabip(buf+38, 64*2); memcpy(h-><a class="code" href="structecat7__2Dnormheader.html#a020286ea0efc6ef7946ca7120f193a57">fill_cti</a>, buf+38, 64*2); <a name="l00575"></a>00575 <span class="keywordflow">if</span>(little) swabip(buf+166, 123*2); memcpy(h-><a class="code" href="structecat7__2Dnormheader.html#a020286ea0efc6ef7946ca7120f193a57">fill_cti</a>, buf+166, 123*2); <a name="l00576"></a>00576 <span class="keywordflow">if</span>(little) swabip(buf+412, 50*2); memcpy(h-><a class="code" href="structecat7__2Dnormheader.html#a10847d8706567154e91057fc98012254">fill_user</a>, buf+412, 50*2); <a name="l00577"></a>00577 <span class="keywordflow">return</span>(0); <a name="l00578"></a>00578 } <a name="l00579"></a>00579 <span class="comment">/*****************************************************************************/</span> <a name="l00580"></a>00580 <a name="l00581"></a>00581 <span class="comment">/*****************************************************************************/</span> <a name="l00595"></a><a class="code" href="ecat7_8h.html#a4bf3a2271415eb36df5e2cee74170414">00595</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a4bf3a2271415eb36df5e2cee74170414">ecat7ReadMatrixdata</a>(FILE *fp, <span class="keywordtype">int</span> start_block, <span class="keywordtype">int</span> block_nr, <span class="keywordtype">char</span> *data, <span class="keywordtype">int</span> dtype) { <a name="l00596"></a>00596 <span class="keywordtype">int</span> i, n, little, err=0; <a name="l00597"></a>00597 <span class="keywordtype">char</span> *cptr; <a name="l00598"></a>00598 <span class="keywordtype">float</span> f; <a name="l00599"></a>00599 <a name="l00600"></a>00600 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadMatrixdata(fp, %d, %d, data, %d)\n"</span>, <a name="l00601"></a>00601 start_block, block_nr, dtype); <a name="l00602"></a>00602 <span class="comment">/* Check the arguments */</span> <a name="l00603"></a>00603 <span class="keywordflow">if</span>(block_nr<=0 || start_block<1 || data==NULL) <span class="keywordflow">return</span>(1); <a name="l00604"></a>00604 <span class="comment">/* Seek the first data block */</span> <a name="l00605"></a>00605 fseek(fp, (start_block-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, SEEK_SET); <a name="l00606"></a>00606 <span class="keywordflow">if</span>(ftell(fp)!=(start_block-1)*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>) <span class="keywordflow">return</span>(9); <a name="l00607"></a>00607 <span class="comment">/* Read the data blocks */</span> <a name="l00608"></a>00608 <span class="keywordflow">if</span>(fread(data, <a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>, block_nr, fp) < block_nr) <span class="keywordflow">return</span>(2); <a name="l00609"></a>00609 <span class="comment">/* Translate data if necessary */</span> <a name="l00610"></a>00610 little=little_endian(); <a name="l00611"></a>00611 <span class="keywordflow">switch</span>(dtype) { <a name="l00612"></a>00612 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>: <span class="comment">/* byte format...no translation necessary */</span> <a name="l00613"></a>00613 <span class="keywordflow">break</span>; <a name="l00614"></a>00614 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a>: <span class="comment">/* byte conversion necessary on big endian platform */</span> <a name="l00615"></a>00615 <span class="keywordflow">if</span>(!little) {cptr=data; swabip(cptr, block_nr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>);} <a name="l00616"></a>00616 <span class="keywordflow">break</span>; <a name="l00617"></a>00617 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a>: <a name="l00618"></a>00618 <span class="keywordflow">for</span>(i=0, cptr=data; i<block_nr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>; i+=4, cptr+=4) { <a name="l00619"></a>00619 n=<a class="code" href="ecat7r_8c.html#a61ea7d463c2bdbdae837a54d6fe44f4c">ecat7rInt</a>(cptr, 1, little); memcpy(cptr, &n, 4); <a name="l00620"></a>00620 } <a name="l00621"></a>00621 <span class="keywordflow">break</span>; <a name="l00622"></a>00622 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a>: <a name="l00623"></a>00623 <span class="keywordflow">for</span>(i=0, cptr=data; i<block_nr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>; i+=4, cptr+=4) { <a name="l00624"></a>00624 f=<a class="code" href="ecat7r_8c.html#a460559c9d89abe660cd05157ce8f1fa2">ecat7rFloat</a>(cptr, 1, little); memcpy(cptr, &f, 4); <a name="l00625"></a>00625 } <a name="l00626"></a>00626 <span class="keywordflow">break</span>; <a name="l00627"></a>00627 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>: <span class="comment">/* IEEE float ; byte conversion necessary on little endian platforms */</span> <a name="l00628"></a>00628 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>: <span class="comment">/* SUN int ; byte conversion necessary on little endian platforms */</span> <a name="l00629"></a>00629 <span class="keywordflow">if</span>(little) swawbip(data, block_nr*MatBLKSIZE); <a name="l00630"></a>00630 <span class="keywordflow">break</span>; <a name="l00631"></a>00631 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>: <span class="comment">/* SUN short ; byte conversion necessary on little endian platforms */</span> <a name="l00632"></a>00632 <span class="keywordflow">if</span>(little) swabip(data, block_nr*MatBLKSIZE); <a name="l00633"></a>00633 <span class="keywordflow">break</span>; <a name="l00634"></a>00634 <span class="keywordflow">default</span>: <span class="comment">/* if something else, for now think it as an error */</span> <a name="l00635"></a>00635 err=2; <a name="l00636"></a>00636 <span class="keywordflow">break</span>; <a name="l00637"></a>00637 } <a name="l00638"></a>00638 <span class="keywordflow">return</span>(err); <a name="l00639"></a>00639 } <a name="l00640"></a>00640 <span class="comment">/*****************************************************************************/</span> <a name="l00641"></a>00641 <a name="l00642"></a>00642 <span class="comment">/*****************************************************************************/</span> <a name="l00656"></a><a class="code" href="ecat7_8h.html#ac62693d7c451818000deb09e2f59cc68">00656</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#ac62693d7c451818000deb09e2f59cc68">ecat7ReadImageMatrix</a>(FILE *fp, <span class="keywordtype">int</span> first_block, <span class="keywordtype">int</span> last_block, <a class="code" href="structecat7__imageheader.html">ECAT7_imageheader</a> *h, <span class="keywordtype">float</span> **fdata) { <a name="l00657"></a>00657 <span class="keywordtype">int</span> i, ret, blockNr, pxlNr; <a name="l00658"></a>00658 <span class="keywordtype">char</span> *mdata, *mptr; <a name="l00659"></a>00659 <span class="keywordtype">float</span> *_fdata, *fptr; <a name="l00660"></a>00660 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sptr; <a name="l00661"></a>00661 <span class="keywordtype">int</span> *iptr; <a name="l00662"></a>00662 <a name="l00663"></a>00663 <a name="l00664"></a>00664 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadImageMatrix(fp, %d, %d, hdr, fdata)\n"</span>, <a name="l00665"></a>00665 first_block, last_block); <a name="l00666"></a>00666 <span class="keywordflow">if</span>(fp==NULL || first_block<=<a class="code" href="ecat63_8h.html#a08ed2cc1d5f72492fb6b4eceb66665ad">MatFirstDirBlk</a> || h==NULL) { <a name="l00667"></a>00667 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid function parameter.\n"</span>); <a name="l00668"></a>00668 <span class="keywordflow">return</span>(1); <a name="l00669"></a>00669 } <a name="l00670"></a>00670 *fdata=(<span class="keywordtype">float</span>*)NULL; <a name="l00671"></a>00671 <a name="l00672"></a>00672 <span class="comment">/* Read subheader */</span> <a name="l00673"></a>00673 ret=<a class="code" href="ecat7r_8c.html#ac96bedc17610a5a05d95ce7742d44b5c">ecat7ReadImageheader</a>(fp, first_block, h); <a name="l00674"></a>00674 <span class="keywordflow">if</span>(ret) { <a name="l00675"></a>00675 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read subheader (%d).\n"</span>, ret); <a name="l00676"></a>00676 <span class="keywordflow">return</span>(5); <a name="l00677"></a>00677 } <a name="l00678"></a>00678 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>>4) <a class="code" href="ecat7p_8c.html#a0854c4b26b52aba47971deb210c1a5a7">ecat7PrintImageheader</a>(h, stdout); <a name="l00679"></a>00679 pxlNr=h-><a class="code" href="structecat7__imageheader.html#a37550f6ce6a114de12aa0d038f573230">x_dimension</a>*h-><a class="code" href="structecat7__imageheader.html#a217d17f8b1f921fad128ce923fa21c35">y_dimension</a>; <a name="l00680"></a>00680 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__imageheader.html#acabca0e942a462b769c59804250b9a1c">num_dimensions</a>>2) pxlNr*=h-><a class="code" href="structecat7__imageheader.html#a7d4bc496c4aa5581d41dff3e5b40a231">z_dimension</a>; <a name="l00681"></a>00681 <span class="keywordflow">if</span>(pxlNr<=0) { <a name="l00682"></a>00682 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid matrix dimension.\n"</span>); <a name="l00683"></a>00683 <span class="keywordflow">return</span>(6); <a name="l00684"></a>00684 } <a name="l00685"></a>00685 <a name="l00686"></a>00686 <span class="comment">/* Read matrix data */</span> <a name="l00687"></a>00687 blockNr=last_block-first_block; <span class="keywordflow">if</span>(blockNr<1) <span class="keywordflow">return</span>(0); <a name="l00688"></a>00688 mdata=(<span class="keywordtype">char</span>*)malloc(blockNr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>); <a name="l00689"></a>00689 <span class="keywordflow">if</span>(mdata==NULL) { <a name="l00690"></a>00690 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00691"></a>00691 <span class="keywordflow">return</span>(8); <a name="l00692"></a>00692 } <a name="l00693"></a>00693 mptr=mdata; <a name="l00694"></a>00694 ret=<a class="code" href="ecat7r_8c.html#a4bf3a2271415eb36df5e2cee74170414">ecat7ReadMatrixdata</a>(fp, first_block+1, blockNr, mptr, h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>); <a name="l00695"></a>00695 <span class="keywordflow">if</span>(ret || mdata==NULL) { <a name="l00696"></a>00696 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read matrix data (%d).\n"</span>, ret); <a name="l00697"></a>00697 free(mdata); <span class="keywordflow">return</span>(9); <a name="l00698"></a>00698 } <a name="l00699"></a>00699 <a name="l00700"></a>00700 <span class="comment">/* Allocate memory for float data */</span> <a name="l00701"></a>00701 _fdata=(<span class="keywordtype">float</span>*)malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)); <a name="l00702"></a>00702 <span class="keywordflow">if</span>(_fdata==NULL) { <a name="l00703"></a>00703 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00704"></a>00704 free(mdata); <span class="keywordflow">return</span>(11); <a name="l00705"></a>00705 } <a name="l00706"></a>00706 <a name="l00707"></a>00707 <span class="comment">/* Convert matrix data to floats */</span> <a name="l00708"></a>00708 fptr=_fdata; mptr=mdata; <a name="l00709"></a>00709 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>) { <a name="l00710"></a>00710 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr++, fptr++) <a name="l00711"></a>00711 *fptr=h-><a class="code" href="structecat7__imageheader.html#a2e44d8de4b5809864432376f713968ec">scale_factor</a>*(<span class="keywordtype">float</span>)(*mptr); <a name="l00712"></a>00712 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a> || h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>) { <a name="l00713"></a>00713 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=2, fptr++) { <a name="l00714"></a>00714 sptr=(<span class="keywordtype">short</span> <span class="keywordtype">int</span>*)mptr; <a name="l00715"></a>00715 *fptr=h-><a class="code" href="structecat7__imageheader.html#a2e44d8de4b5809864432376f713968ec">scale_factor</a>*(float)(*sptr); <a name="l00716"></a>00716 } <a name="l00717"></a>00717 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a> || h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>) { <a name="l00718"></a>00718 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=4, fptr++) { <a name="l00719"></a>00719 iptr=(<span class="keywordtype">int</span>*)mptr; <a name="l00720"></a>00720 *fptr=h-><a class="code" href="structecat7__imageheader.html#a2e44d8de4b5809864432376f713968ec">scale_factor</a>*(float)(*iptr); <a name="l00721"></a>00721 } <a name="l00722"></a>00722 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a> || h-><a class="code" href="structecat7__imageheader.html#a9b90dc8511beaf487a57afc61044d09a">data_type</a>==<a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>) { <a name="l00723"></a>00723 memcpy(fptr, mptr, pxlNr*4); <a name="l00724"></a>00724 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, fptr++) *fptr *= h-><a class="code" href="structecat7__imageheader.html#a2e44d8de4b5809864432376f713968ec">scale_factor</a>; <a name="l00725"></a>00725 } <a name="l00726"></a>00726 free(mdata); <a name="l00727"></a>00727 *fdata=_fdata; <a name="l00728"></a>00728 <a name="l00729"></a>00729 <span class="keywordflow">return</span>(0); <a name="l00730"></a>00730 } <a name="l00731"></a>00731 <span class="comment">/*****************************************************************************/</span> <a name="l00732"></a>00732 <a name="l00733"></a>00733 <span class="comment">/*****************************************************************************/</span> <a name="l00749"></a><a class="code" href="ecat7_8h.html#a2cdc756d0ec7d7b4d22f7d6c8b73fd7b">00749</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a2cdc756d0ec7d7b4d22f7d6c8b73fd7b">ecat7Read2DScanMatrix</a>(FILE *fp, <span class="keywordtype">int</span> first_block, <span class="keywordtype">int</span> last_block, <a name="l00750"></a>00750 <a class="code" href="structecat7__2Dscanheader.html">ECAT7_2Dscanheader</a> *h, <span class="keywordtype">float</span> **fdata) { <a name="l00751"></a>00751 <span class="keywordtype">int</span> i, ret, blockNr, pxlNr; <a name="l00752"></a>00752 <span class="keywordtype">char</span> *mdata, *mptr; <a name="l00753"></a>00753 <span class="keywordtype">float</span> *_fdata, *fptr; <a name="l00754"></a>00754 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sptr; <a name="l00755"></a>00755 <span class="keywordtype">int</span> *iptr; <a name="l00756"></a>00756 <a name="l00757"></a>00757 <a name="l00758"></a>00758 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7Read2DScanMatrix(fp, %d, %d, hdr, fdata)\n"</span>, <a name="l00759"></a>00759 first_block, last_block); <a name="l00760"></a>00760 <span class="keywordflow">if</span>(fp==NULL || first_block<=<a class="code" href="ecat63_8h.html#a08ed2cc1d5f72492fb6b4eceb66665ad">MatFirstDirBlk</a> || h==NULL) { <a name="l00761"></a>00761 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid function parameter.\n"</span>); <a name="l00762"></a>00762 <span class="keywordflow">return</span>(1); <a name="l00763"></a>00763 } <a name="l00764"></a>00764 *fdata=(<span class="keywordtype">float</span>*)NULL; <a name="l00765"></a>00765 <a name="l00766"></a>00766 <span class="comment">/* Read subheader */</span> <a name="l00767"></a>00767 ret=<a class="code" href="ecat7r_8c.html#aa37d19fc4f30868ea0b062d39d9d8450">ecat7Read2DScanheader</a>(fp, first_block, h); <a name="l00768"></a>00768 <span class="keywordflow">if</span>(ret) { <a name="l00769"></a>00769 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read subheader (%d).\n"</span>, ret); <a name="l00770"></a>00770 <span class="keywordflow">return</span>(5); <a name="l00771"></a>00771 } <a name="l00772"></a>00772 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>>4) <a class="code" href="ecat7p_8c.html#aaa6831e8494e81833098b42b37e3f96a">ecat7Print2DScanheader</a>(h, stdout); <a name="l00773"></a>00773 pxlNr=h-><a class="code" href="structecat7__2Dscanheader.html#adfe07cff44e3145ca65c234079ff255a">num_r_elements</a>*h-><a class="code" href="structecat7__2Dscanheader.html#ad0d1958b231f10aa68d60c653ec30627">num_angles</a>; <a name="l00774"></a>00774 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__2Dscanheader.html#a592447ad964ba95e2a29a7d3d2284576">num_dimensions</a>>2) pxlNr*=h-><a class="code" href="structecat7__2Dscanheader.html#a697bb7a95d88bce42bf80e2d6bd69037">num_z_elements</a>; <a name="l00775"></a>00775 <span class="keywordflow">if</span>(pxlNr<=0) { <a name="l00776"></a>00776 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid matrix dimension.\n"</span>); <a name="l00777"></a>00777 <span class="keywordflow">return</span>(6); <a name="l00778"></a>00778 } <a name="l00779"></a>00779 <a name="l00780"></a>00780 <span class="comment">/* Read matrix data */</span> <a name="l00781"></a>00781 blockNr=last_block-first_block; <span class="keywordflow">if</span>(blockNr<1) <span class="keywordflow">return</span>(0); <a name="l00782"></a>00782 mdata=(<span class="keywordtype">char</span>*)malloc(blockNr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>); <a name="l00783"></a>00783 <span class="keywordflow">if</span>(mdata==NULL) { <a name="l00784"></a>00784 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00785"></a>00785 <span class="keywordflow">return</span>(8); <a name="l00786"></a>00786 } <a name="l00787"></a>00787 mptr=mdata; <a name="l00788"></a>00788 ret=<a class="code" href="ecat7r_8c.html#a4bf3a2271415eb36df5e2cee74170414">ecat7ReadMatrixdata</a>(fp, first_block+1, blockNr, mptr, h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>); <a name="l00789"></a>00789 <span class="keywordflow">if</span>(ret || mdata==NULL) { <a name="l00790"></a>00790 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read matrix data (%d).\n"</span>, ret); <a name="l00791"></a>00791 free(mdata); <span class="keywordflow">return</span>(9); <a name="l00792"></a>00792 } <a name="l00793"></a>00793 <a name="l00794"></a>00794 <span class="comment">/* Allocate memory for float data */</span> <a name="l00795"></a>00795 _fdata=(<span class="keywordtype">float</span>*)malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)); <a name="l00796"></a>00796 <span class="keywordflow">if</span>(_fdata==NULL) { <a name="l00797"></a>00797 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00798"></a>00798 free(mdata); <span class="keywordflow">return</span>(11); <a name="l00799"></a>00799 } <a name="l00800"></a>00800 <a name="l00801"></a>00801 <span class="comment">/* Convert matrix data to floats */</span> <a name="l00802"></a>00802 fptr=_fdata; mptr=mdata; <a name="l00803"></a>00803 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>) { <a name="l00804"></a>00804 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr++, fptr++) <a name="l00805"></a>00805 *fptr=h-><a class="code" href="structecat7__2Dscanheader.html#ae27dc6ec3582a4e2a0e69ea9377993a3">scale_factor</a>*(<span class="keywordtype">float</span>)(*mptr); <a name="l00806"></a>00806 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a> || h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>) { <a name="l00807"></a>00807 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=2, fptr++) { <a name="l00808"></a>00808 sptr=(<span class="keywordtype">short</span> <span class="keywordtype">int</span>*)mptr; <a name="l00809"></a>00809 *fptr=h-><a class="code" href="structecat7__2Dscanheader.html#ae27dc6ec3582a4e2a0e69ea9377993a3">scale_factor</a>*(float)(*sptr); <a name="l00810"></a>00810 } <a name="l00811"></a>00811 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a> || h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>) { <a name="l00812"></a>00812 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=4, fptr++) { <a name="l00813"></a>00813 iptr=(<span class="keywordtype">int</span>*)mptr; <a name="l00814"></a>00814 *fptr=h-><a class="code" href="structecat7__2Dscanheader.html#ae27dc6ec3582a4e2a0e69ea9377993a3">scale_factor</a>*(float)(*iptr); <a name="l00815"></a>00815 } <a name="l00816"></a>00816 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a> || h-><a class="code" href="structecat7__2Dscanheader.html#aa13b8697f5c56f7e31e942cafb91bddc">data_type</a>==<a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>) { <a name="l00817"></a>00817 memcpy(fptr, mptr, pxlNr*4); <a name="l00818"></a>00818 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, fptr++) *fptr *= h-><a class="code" href="structecat7__2Dscanheader.html#ae27dc6ec3582a4e2a0e69ea9377993a3">scale_factor</a>; <a name="l00819"></a>00819 } <a name="l00820"></a>00820 free(mdata); <a name="l00821"></a>00821 *fdata=_fdata; <a name="l00822"></a>00822 <a name="l00823"></a>00823 <span class="keywordflow">return</span>(0); <a name="l00824"></a>00824 } <a name="l00825"></a>00825 <span class="comment">/*****************************************************************************/</span> <a name="l00826"></a>00826 <a name="l00827"></a>00827 <span class="comment">/*****************************************************************************/</span> <a name="l00844"></a><a class="code" href="ecat7_8h.html#aeed5aa16cd7ac3d07994e039ede3142d">00844</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#aeed5aa16cd7ac3d07994e039ede3142d">ecat7ReadScanMatrix</a>(FILE *fp, <span class="keywordtype">int</span> first_block, <span class="keywordtype">int</span> last_block, <a class="code" href="structecat7__scanheader.html">ECAT7_scanheader</a> *h, <span class="keywordtype">float</span> **fdata) { <a name="l00845"></a>00845 <span class="keywordtype">int</span> i, ret, blockNr, trueblockNr, pxlNr, dimz; <a name="l00846"></a>00846 <span class="keywordtype">char</span> *mdata, *mptr; <a name="l00847"></a>00847 <span class="keywordtype">float</span> *_fdata, *fptr; <a name="l00848"></a>00848 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sptr; <a name="l00849"></a>00849 <span class="keywordtype">int</span> *iptr; <a name="l00850"></a>00850 <a name="l00851"></a>00851 <a name="l00852"></a>00852 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadScanMatrix(fp, %d, %d, hdr, fdata)\n"</span>, <a name="l00853"></a>00853 first_block, last_block); <a name="l00854"></a>00854 <span class="keywordflow">if</span>(fp==NULL || first_block<=<a class="code" href="ecat63_8h.html#a08ed2cc1d5f72492fb6b4eceb66665ad">MatFirstDirBlk</a> || h==NULL) { <a name="l00855"></a>00855 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid function parameter.\n"</span>); <a name="l00856"></a>00856 <span class="keywordflow">return</span>(1); <a name="l00857"></a>00857 } <a name="l00858"></a>00858 *fdata=(<span class="keywordtype">float</span>*)NULL; <a name="l00859"></a>00859 <a name="l00860"></a>00860 <span class="comment">/* Read subheader */</span> <a name="l00861"></a>00861 ret=<a class="code" href="ecat7r_8c.html#a37414b52ef395686b6a81d6229d13143">ecat7ReadScanheader</a>(fp, first_block, h); <a name="l00862"></a>00862 <span class="keywordflow">if</span>(ret) { <a name="l00863"></a>00863 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read subheader (%d).\n"</span>, ret); <a name="l00864"></a>00864 <span class="keywordflow">return</span>(5); <a name="l00865"></a>00865 } <a name="l00866"></a>00866 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>>4) <a class="code" href="ecat7p_8c.html#ade61cec0f7e15550290193ef0d6e26cc">ecat7PrintScanheader</a>(h, stdout); <a name="l00867"></a>00867 pxlNr=h-><a class="code" href="structecat7__scanheader.html#aee3fd0c2360973b3d53ad1adf6225d7d">num_r_elements</a>*h-><a class="code" href="structecat7__scanheader.html#a4f58fe54c66427d34390f6530a12f8ed">num_angles</a>; <a name="l00868"></a>00868 <span class="keywordflow">for</span>(i=dimz=0; i<64; i++) dimz+=h-><a class="code" href="structecat7__scanheader.html#a714bae4e56b451926ceeb0de19222ad7">num_z_elements</a>[i]; pxlNr*=dimz; <a name="l00869"></a>00869 <span class="keywordflow">if</span>(pxlNr<=0) { <a name="l00870"></a>00870 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"invalid matrix dimension.\n"</span>); <a name="l00871"></a>00871 <span class="keywordflow">return</span>(6); <a name="l00872"></a>00872 } <a name="l00873"></a>00873 trueblockNr=pxlNr*<a class="code" href="ecat7r_8c.html#a17795a6c58bfa4c1159e0fe642098678">ecat7pxlbytes</a>(h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>); <a name="l00874"></a>00874 trueblockNr=(trueblockNr+<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>-1)/<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>; <a name="l00875"></a>00875 <a name="l00876"></a>00876 <span class="comment">/* Read matrix data; note that header takes 2 blocks */</span> <a name="l00877"></a>00877 blockNr=last_block-first_block-1; <span class="keywordflow">if</span>(blockNr<1) <span class="keywordflow">return</span>(0); <a name="l00878"></a>00878 <span class="keywordflow">if</span>(blockNr<trueblockNr) trueblockNr=blockNr; <a name="l00879"></a>00879 mdata=(<span class="keywordtype">char</span>*)malloc(blockNr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>); <a name="l00880"></a>00880 <span class="keywordflow">if</span>(mdata==NULL) { <a name="l00881"></a>00881 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00882"></a>00882 <span class="keywordflow">return</span>(8); <a name="l00883"></a>00883 } <a name="l00884"></a>00884 mptr=mdata; <span class="comment">/* note that only true block nr is read! */</span> <a name="l00885"></a>00885 ret=<a class="code" href="ecat7r_8c.html#a4bf3a2271415eb36df5e2cee74170414">ecat7ReadMatrixdata</a>(fp, first_block+2, trueblockNr, mptr, h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>); <a name="l00886"></a>00886 <span class="keywordflow">if</span>(ret || mdata==NULL) { <a name="l00887"></a>00887 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read matrix data (%d).\n"</span>, ret); <a name="l00888"></a>00888 free(mdata); <span class="keywordflow">return</span>(9); <a name="l00889"></a>00889 } <a name="l00890"></a>00890 <a name="l00891"></a>00891 <span class="comment">/* Allocate memory for float data */</span> <a name="l00892"></a>00892 _fdata=(<span class="keywordtype">float</span>*)malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)); <a name="l00893"></a>00893 <span class="keywordflow">if</span>(_fdata==NULL) { <a name="l00894"></a>00894 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00895"></a>00895 free(mdata); <span class="keywordflow">return</span>(11); <a name="l00896"></a>00896 } <a name="l00897"></a>00897 <a name="l00898"></a>00898 <span class="comment">/* Convert matrix data to floats */</span> <a name="l00899"></a>00899 fptr=_fdata; mptr=mdata; <a name="l00900"></a>00900 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>) { <a name="l00901"></a>00901 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr++, fptr++) <a name="l00902"></a>00902 *fptr=h-><a class="code" href="structecat7__scanheader.html#a77ed7b5edd51bf8feeea367ec4d143bb">scale_factor</a>*(<span class="keywordtype">float</span>)(*mptr); <a name="l00903"></a>00903 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a> || h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>) { <a name="l00904"></a>00904 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=2, fptr++) { <a name="l00905"></a>00905 sptr=(<span class="keywordtype">short</span> <span class="keywordtype">int</span>*)mptr; <a name="l00906"></a>00906 *fptr=h-><a class="code" href="structecat7__scanheader.html#a77ed7b5edd51bf8feeea367ec4d143bb">scale_factor</a>*(float)(*sptr); <a name="l00907"></a>00907 } <a name="l00908"></a>00908 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a> || h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>) { <a name="l00909"></a>00909 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=4, fptr++) { <a name="l00910"></a>00910 iptr=(<span class="keywordtype">int</span>*)mptr; <a name="l00911"></a>00911 *fptr=h-><a class="code" href="structecat7__scanheader.html#a77ed7b5edd51bf8feeea367ec4d143bb">scale_factor</a>*(float)(*iptr); <a name="l00912"></a>00912 } <a name="l00913"></a>00913 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a> || h-><a class="code" href="structecat7__scanheader.html#ac1dfbebf9227256db4dd3ec96a95277f">data_type</a>==<a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>) { <a name="l00914"></a>00914 memcpy(fptr, mptr, pxlNr*4); <a name="l00915"></a>00915 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, fptr++) *fptr *= h-><a class="code" href="structecat7__scanheader.html#a77ed7b5edd51bf8feeea367ec4d143bb">scale_factor</a>; <a name="l00916"></a>00916 } <a name="l00917"></a>00917 free(mdata); <a name="l00918"></a>00918 *fdata=_fdata; <a name="l00919"></a>00919 <a name="l00920"></a>00920 <span class="keywordflow">return</span>(0); <a name="l00921"></a>00921 } <a name="l00922"></a>00922 <span class="comment">/*****************************************************************************/</span> <a name="l00923"></a>00923 <a name="l00924"></a>00924 <span class="comment">/*****************************************************************************/</span> <a name="l00939"></a><a class="code" href="ecat7_8h.html#aa07f0b9089fda6bc1f5e8ce179f62686">00939</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#aa07f0b9089fda6bc1f5e8ce179f62686">ecat7ReadPolarmapMatrix</a>(FILE *fp, <span class="keywordtype">int</span> first_block, <span class="keywordtype">int</span> last_block, <a class="code" href="structecat7__polmapheader.html">ECAT7_polmapheader</a> *h, <span class="keywordtype">float</span> **fdata) { <a name="l00940"></a>00940 <span class="keywordtype">int</span> i, ret, blockNr, pxlNr; <a name="l00941"></a>00941 <span class="keywordtype">char</span> *mdata, *mptr; <a name="l00942"></a>00942 <span class="keywordtype">float</span> *_fdata, *fptr; <a name="l00943"></a>00943 <span class="keywordtype">short</span> <span class="keywordtype">int</span> *sptr; <a name="l00944"></a>00944 <span class="keywordtype">int</span> *iptr; <a name="l00945"></a>00945 <a name="l00946"></a>00946 <a name="l00947"></a>00947 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>) printf(<span class="stringliteral">"ecat7ReadPolarmapMatrix(fp, %d, %d, hdr, fdata)\n"</span>, <a name="l00948"></a>00948 first_block, last_block); <a name="l00949"></a>00949 <span class="keywordflow">if</span>(fp==NULL || first_block<=<a class="code" href="ecat63_8h.html#a08ed2cc1d5f72492fb6b4eceb66665ad">MatFirstDirBlk</a> || h==NULL) <span class="keywordflow">return</span> 1; <a name="l00950"></a>00950 *fdata=(<span class="keywordtype">float</span>*)NULL; <a name="l00951"></a>00951 <a name="l00952"></a>00952 <span class="comment">/* Read subheader */</span> <a name="l00953"></a>00953 ret=<a class="code" href="ecat7r_8c.html#a498d383ca08a832895af8ead01e8c756">ecat7ReadPolmapheader</a>(fp, first_block, h); <a name="l00954"></a>00954 <span class="keywordflow">if</span>(ret) { <a name="l00955"></a>00955 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot read subheader (%d).\n"</span>, ret); <a name="l00956"></a>00956 <span class="keywordflow">return</span> 2; <a name="l00957"></a>00957 } <a name="l00958"></a>00958 <span class="keywordflow">if</span>(<a class="code" href="ecat7_8h.html#aeba6b1646ad1dd69809d9ba3a8732939">ECAT7_TEST</a>>4) <a class="code" href="ecat7p_8c.html#a7b20d17081345d4f4ede5d4791e18e31">ecat7PrintPolmapheader</a>(h, stdout); <a name="l00959"></a>00959 <span class="keywordflow">for</span>(i=pxlNr=0; i<h-><a class="code" href="structecat7__polmapheader.html#a503caad4ac64580ff3508463ae8e2342">num_rings</a>; i++) pxlNr+=h-><a class="code" href="structecat7__polmapheader.html#a7888f1b2bc8adcc5bcbf918bbbb82dbb">sectors_per_ring</a>[i]; <a name="l00960"></a>00960 <span class="keywordflow">if</span>(pxlNr<=0) <span class="keywordflow">return</span> 3; <a name="l00961"></a>00961 <a name="l00962"></a>00962 <span class="comment">/* Read matrix data */</span> <a name="l00963"></a>00963 blockNr=last_block-first_block; <span class="keywordflow">if</span>(blockNr<1) <span class="keywordflow">return</span> 0; <a name="l00964"></a>00964 mdata=(<span class="keywordtype">char</span>*)malloc(blockNr*<a class="code" href="ecat63_8h.html#abdf75c7aade8a596594e8e7ce8b16872">MatBLKSIZE</a>); <a name="l00965"></a>00965 <span class="keywordflow">if</span>(mdata==NULL) <span class="keywordflow">return</span> 4; <a name="l00966"></a>00966 mptr=mdata; <a name="l00967"></a>00967 ret=<a class="code" href="ecat7r_8c.html#a4bf3a2271415eb36df5e2cee74170414">ecat7ReadMatrixdata</a>(fp, first_block+1, blockNr, mptr, h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>); <a name="l00968"></a>00968 <span class="keywordflow">if</span>(ret || mdata==NULL) { <a name="l00969"></a>00969 <span class="keywordflow">if</span>(mdata!=NULL) free(mdata); <a name="l00970"></a>00970 <span class="keywordflow">return</span> 5; <a name="l00971"></a>00971 } <a name="l00972"></a>00972 <a name="l00973"></a>00973 <span class="comment">/* Allocate memory for float data */</span> <a name="l00974"></a>00974 _fdata=(<span class="keywordtype">float</span>*)malloc(pxlNr*<span class="keyword">sizeof</span>(<span class="keywordtype">float</span>)); <a name="l00975"></a>00975 <span class="keywordflow">if</span>(_fdata==NULL) { <a name="l00976"></a>00976 sprintf(<a class="code" href="ecat7_8h.html#a66c8b21d11aea27d874d191a12312207">ecat7errmsg</a>, <span class="stringliteral">"cannot allocate memory.\n"</span>); <a name="l00977"></a>00977 free(mdata); <span class="keywordflow">return</span> 4; <a name="l00978"></a>00978 } <a name="l00979"></a>00979 <a name="l00980"></a>00980 <span class="comment">/* Convert matrix data to floats */</span> <a name="l00981"></a>00981 fptr=_fdata; mptr=mdata; <a name="l00982"></a>00982 <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>) { <a name="l00983"></a>00983 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr++, fptr++) <a name="l00984"></a>00984 *fptr=h-><a class="code" href="structecat7__polmapheader.html#ad2f57edc57a8a6f328e19ab35502eb20">scale_factor</a>*(<span class="keywordtype">float</span>)(*mptr); <a name="l00985"></a>00985 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a> || h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>) { <a name="l00986"></a>00986 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=2, fptr++) { <a name="l00987"></a>00987 sptr=(<span class="keywordtype">short</span> <span class="keywordtype">int</span>*)mptr; <a name="l00988"></a>00988 *fptr=h-><a class="code" href="structecat7__polmapheader.html#ad2f57edc57a8a6f328e19ab35502eb20">scale_factor</a>*(float)(*sptr); <a name="l00989"></a>00989 } <a name="l00990"></a>00990 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a> || h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>) { <a name="l00991"></a>00991 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, mptr+=4, fptr++) { <a name="l00992"></a>00992 iptr=(<span class="keywordtype">int</span>*)mptr; <a name="l00993"></a>00993 *fptr=h-><a class="code" href="structecat7__polmapheader.html#ad2f57edc57a8a6f328e19ab35502eb20">scale_factor</a>*(float)(*iptr); <a name="l00994"></a>00994 } <a name="l00995"></a>00995 } <span class="keywordflow">else</span> <span class="keywordflow">if</span>(h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a> || h-><a class="code" href="structecat7__polmapheader.html#a85514ffcdb6dc106fe8d742bcb248947">data_type</a>==<a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>) { <a name="l00996"></a>00996 memcpy(fptr, mptr, pxlNr*4); <a name="l00997"></a>00997 <span class="keywordflow">for</span>(i=0; i<pxlNr; i++, fptr++) *fptr *= h-><a class="code" href="structecat7__polmapheader.html#ad2f57edc57a8a6f328e19ab35502eb20">scale_factor</a>; <a name="l00998"></a>00998 } <a name="l00999"></a>00999 free(mdata); <a name="l01000"></a>01000 *fdata=_fdata; <a name="l01001"></a>01001 <a name="l01002"></a>01002 <span class="keywordflow">return</span> 0; <a name="l01003"></a>01003 } <a name="l01004"></a>01004 <span class="comment">/*****************************************************************************/</span> <a name="l01005"></a>01005 <a name="l01006"></a>01006 <span class="comment">/*****************************************************************************/</span> <a name="l01015"></a><a class="code" href="ecat7_8h.html#a460559c9d89abe660cd05157ce8f1fa2">01015</a> <span class="keywordtype">float</span> <a class="code" href="ecat7r_8c.html#a460559c9d89abe660cd05157ce8f1fa2">ecat7rFloat</a>(<span class="keywordtype">void</span> *bufi, <span class="keywordtype">int</span> isvax, <span class="keywordtype">int</span> islittle) { <a name="l01016"></a>01016 <span class="keyword">union </span>{<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> ul; <span class="keywordtype">float</span> f;} t; <a name="l01017"></a>01017 <a name="l01018"></a>01018 memcpy(&t.ul, bufi, 4); <span class="keywordflow">if</span>(t.ul==0) {<span class="keywordflow">return</span>(0.0);} <a name="l01019"></a>01019 <span class="keywordflow">if</span>(isvax) { <span class="comment">/* if input is in VAX format */</span> <a name="l01020"></a>01020 <span class="comment">/* Swap words on i386 and bytes on SUN */</span> <a name="l01021"></a>01021 <span class="keywordflow">if</span>(islittle) swawip(&t.ul, 4); <span class="keywordflow">else</span> swabip(&t.ul, 4); <a name="l01022"></a>01022 t.ul-=(2L<<23); <span class="comment">/* subtract 2 from exp */</span> <a name="l01023"></a>01023 } <span class="keywordflow">else</span> { <span class="comment">/* input is in i386 format */</span> <a name="l01024"></a>01024 <span class="keywordflow">if</span>(!islittle) swawbip(&t.ul, 4); <span class="comment">/* Switch words and bytes on SUN */</span> <a name="l01025"></a>01025 } <a name="l01026"></a>01026 <span class="keywordflow">return</span>(t.f); <a name="l01027"></a>01027 } <a name="l01028"></a>01028 <a name="l01038"></a><a class="code" href="ecat7_8h.html#a61ea7d463c2bdbdae837a54d6fe44f4c">01038</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a61ea7d463c2bdbdae837a54d6fe44f4c">ecat7rInt</a>(<span class="keywordtype">void</span> *bufi, <span class="keywordtype">int</span> isvax, <span class="keywordtype">int</span> islittle) { <a name="l01039"></a>01039 <span class="keywordtype">int</span> i; <a name="l01040"></a>01040 <a name="l01041"></a>01041 <span class="comment">/* Swap both words and bytes on SUN */</span> <a name="l01042"></a>01042 memcpy(&i, bufi, 4); <span class="keywordflow">if</span>(!islittle) swawbip(&i, 4); <a name="l01043"></a>01043 <span class="keywordflow">return</span>(i); <a name="l01044"></a>01044 } <a name="l01045"></a>01045 <span class="comment">/*****************************************************************************/</span> <a name="l01046"></a>01046 <a name="l01047"></a>01047 <span class="comment">/*****************************************************************************/</span> <a name="l01055"></a><a class="code" href="ecat7_8h.html#a17795a6c58bfa4c1159e0fe642098678">01055</a> <span class="keywordtype">int</span> <a class="code" href="ecat7r_8c.html#a17795a6c58bfa4c1159e0fe642098678">ecat7pxlbytes</a>(<span class="keywordtype">short</span> <span class="keywordtype">int</span> data_type) { <a name="l01056"></a>01056 <span class="keywordtype">int</span> byteNr=0; <a name="l01057"></a>01057 <span class="keywordflow">switch</span>(data_type) { <a name="l01058"></a>01058 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#aa620685ae5b24a3e4410c109eb6a7b2a">ECAT7_BYTE</a>: byteNr=1; <span class="keywordflow">break</span>; <a name="l01059"></a>01059 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#ab6ce9728e8574761efa8303ff6a34fe9">ECAT7_VAXI2</a>: <a name="l01060"></a>01060 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#aee0f2ad8e2d9d538fda30aec053a642a">ECAT7_SUNI2</a>: byteNr=2; <span class="keywordflow">break</span>; <a name="l01061"></a>01061 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a3256675b42164c36df6e3810ac8c82fa">ECAT7_VAXI4</a>: <a name="l01062"></a>01062 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#acff598c85518f15df8aec8175730cbe4">ECAT7_VAXR4</a>: <a name="l01063"></a>01063 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a2b1d3c2ff1a9c241f41fe68189b95570">ECAT7_IEEER4</a>: <a name="l01064"></a>01064 <span class="keywordflow">case</span> <a class="code" href="ecat7_8h.html#a81ca91463b7e64939e4970bb48b4b26f">ECAT7_SUNI4</a>: byteNr=4; <span class="keywordflow">break</span>; <a name="l01065"></a>01065 } <a name="l01066"></a>01066 <span class="keywordflow">return</span>(byteNr); <a name="l01067"></a>01067 } <a name="l01068"></a>01068 <span class="comment">/*****************************************************************************/</span> <a name="l01069"></a>01069 <a name="l01070"></a>01070 <span class="comment">/*****************************************************************************/</span> <a name="l01071"></a>01071 </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>