siemens_ima_read.html

mri_toolbox是一个工具用来MRI. 来自于SourceForge, 我上传这个软件,希望能结识对医疗软件感兴趣的兄弟.

0118             <span class="keyword">for</span> y = 1:RowDim,
0119                 x = 1:PixelDim;
0120                 ima.img(x,y,z) = tmp(n:n+(PixelDim-1));
0121                 n = n + PixelDim;
0122             <span class="keyword">end</span>
0123         <span class="keyword">end</span>
0124         
0125     <span class="keyword">end</span>
0126 <span class="keyword">return</span>
0127 
0128 
0129 <a name="_sub2" href="#_subfunctions" class="code">function [ ima, machine ] = ima_hdr_read(fileprefix,machine)</a>
0130 
0131 
0132 
0133 <span class="comment">% Siemens Magnetom Vision Native Format</span>
0134 <span class="comment">%</span>
0135 <span class="comment">% The exact details of the format are not known, but</span>
0136 <span class="comment">% a little guess work has determined what follows. The</span>
0137 <span class="comment">% data types are Sun, hence the byte order is big-endian</span>
0138 <span class="comment">% and the all the floats I have found are doubles. Offsets</span>
0139 <span class="comment">% here are in bytes from the start of the header. The</span>
0140 <span class="comment">% uncompressed image data starts at offset 6144.</span>
0141 <span class="comment">%</span>
0142 <span class="comment">%         0         u_int      SiemensStudyDateYYYY</span>
0143 <span class="comment">%         4         u_int      SiemensStudyDateMM</span>
0144 <span class="comment">%         8         u_int      SiemensStudyDateDD</span>
0145 <span class="comment">%         12        u_int      AcquisitionDateYYYY</span>
0146 <span class="comment">%         16        u_int      AcquisitionDateMM</span>
0147 <span class="comment">%         20        u_int      AcquisitionDateDD</span>
0148 <span class="comment">%         24        u_int      ImageDateYYYY</span>
0149 <span class="comment">%         28        u_int      ImageDateMM</span>
0150 <span class="comment">%         32        u_int      ImageDateDD</span>
0151 <span class="comment">%         36        u_int      SiemensStudyTimeHH</span>
0152 <span class="comment">%         40        u_int      SiemensStudyTimeMM</span>
0153 <span class="comment">%         44        u_int      SiemensStudyTimeSS</span>
0154 <span class="comment">%         52        u_int      AcquisitionTimeHH</span>
0155 <span class="comment">%         56        u_int      AcquisitionTimeMM</span>
0156 <span class="comment">%         60        u_int      AcquisitionTimeSS</span>
0157 <span class="comment">%         68        u_int      ImageTimeHH</span>
0158 <span class="comment">%         72        u_int      ImageTimeMM</span>
0159 <span class="comment">%         76        u_int      ImageTimeSS</span>
0160 <span class="comment">%         96        char[7]    Manufacturer</span>
0161 <span class="comment">%         105       char[25]   InstitutionName</span>
0162 <span class="comment">%         186       char[4]    Annotation</span>
0163 <span class="comment">%         281       char[15]   ModelName</span>
0164 <span class="comment">%         412       u_int      LastMoveDateYYYY</span>
0165 <span class="comment">%         416       u_int      LastMoveDateMM</span>
0166 <span class="comment">%         420       u_int      LastMoveDateDD</span>
0167 <span class="comment">%         424       u_int      LastMoveTimeHH</span>
0168 <span class="comment">%         428       u_int      LastMoveTimeMM</span>
0169 <span class="comment">%         432       u_int      LastMoveTimeSS</span>
0170 <span class="comment">%         768       char[25]   PatientName</span>
0171 <span class="comment">%         795       char[12]   PatientID</span>
0172 <span class="comment">%         808       u_int      DOBYYYY</span>
0173 <span class="comment">%         812       u_int      DOBMM</span>
0174 <span class="comment">%         816       u_int      DOBDD</span>
0175 <span class="comment">%         851       char[3]    PatientAge</span>
0176 <span class="comment">%         854       char       PatientAgeUnits      ('Y'=years)</span>
0177 <span class="comment">%         1052      u_int      RegistrationDateYYYY</span>
0178 <span class="comment">%         1056      u_int      RegistrationDateMM</span>
0179 <span class="comment">%         1060      u_int      RegistrationDateDD</span>
0180 <span class="comment">%         1064      u_int      RegistrationTimeHH</span>
0181 <span class="comment">%         1068      u_int      RegistrationTimeMM</span>
0182 <span class="comment">%         1072      u_int      RegistrationTimeSS</span>
0183 <span class="comment">%         1544      double     SliceThickness</span>
0184 <span class="comment">%         1560      double     RepetitionTime</span>
0185 <span class="comment">%         1568      double     EchoTime</span>
0186 <span class="comment">%         1592      double     FrequencyMHz</span>
0187 <span class="comment">%         1639      char[5]    Station</span>
0188 <span class="comment">%         1712      u_int      CalibrationDateYYYY</span>
0189 <span class="comment">%         1716      u_int      CalibrationDateMM</span>
0190 <span class="comment">%         1720      u_int      CalibrationDateDD</span>
0191 <span class="comment">%         1724      u_int      CalibrationTimeHH</span>
0192 <span class="comment">%         1728      u_int      CalibrationTimeMM</span>
0193 <span class="comment">%         1732      u_int      CalibrationTimeSS</span>
0194 <span class="comment">%         1767      char[16]   ReceivingCoil</span>
0195 <span class="comment">%         1828      char[4]    ImagedNucleus</span>
0196 <span class="comment">%         2112      double     FlipAngle</span>
0197 <span class="comment">%         2560      double     MagneticFieldStrength</span>
0198 <span class="comment">%         2864      u_int      DisplayMatrixSize</span>
0199 <span class="comment">%         2944      char[65]   SequencePrgName</span>
0200 <span class="comment">%         3009      char[65]   SequenceWkcName</span>
0201 <span class="comment">%         3074      char[9]    SequenceAuthor</span>
0202 <span class="comment">%         3083      char[8]    SequenceType</span>
0203 <span class="comment">%         3744      double     FOVRow</span>
0204 <span class="comment">%         3752      double     FOVColumn</span>
0205 <span class="comment">%         3768      double     CenterPointX</span>
0206 <span class="comment">%         3776      double     CenterPointY</span>
0207 <span class="comment">%         3784      double     CenterPointZ</span>
0208 <span class="comment">%         3792      double     NormalVectorX</span>
0209 <span class="comment">%         3800      double     NormalVectorY</span>
0210 <span class="comment">%         3808      double     NormalVectorZ</span>
0211 <span class="comment">%         3816      double     DistanceFromIsocenter</span>
0212 <span class="comment">%         3832      double     RowVectorX</span>
0213 <span class="comment">%         3840      double     RowVectorY</span>
0214 <span class="comment">%         3848      double     RowVectorZ</span>
0215 <span class="comment">%         3856      double     ColumnVectorX</span>
0216 <span class="comment">%         3864      double     ColumnVectorY</span>
0217 <span class="comment">%         3872      double     ColumnVectorZ</span>
0218 <span class="comment">%         3880      char[3]    OrientationSet1Top</span>
0219 <span class="comment">%         3884      char[3]    OrientationSet1Left</span>
0220 <span class="comment">%         3888      char[3]    OrientationSet1Back</span>
0221 <span class="comment">%         3892      char[3]    OrientationSet2Down</span>
0222 <span class="comment">%         3896      char[3]    OrientationSet2Right</span>
0223 <span class="comment">%         3900      char[3]    OrientationSet2Front</span>
0224 <span class="comment">%         3904      char[32]   SequenceName</span>
0225 <span class="comment">%</span>
0226 <span class="comment">%         4136      double     GapRatio (see below)</span>
0227 <span class="comment">%</span>
0228 <span class="comment">%         5000      double     PixelSizeRow</span>
0229 <span class="comment">%         5008      double     PixelSizeColumn</span>
0230 <span class="comment">%</span>
0231 <span class="comment">%         5504      char[12]   TextPatientID</span>
0232 <span class="comment">%         5517      char       TextPatientSex</span>
0233 <span class="comment">%         5518      char[3]    TextPatientAge</span>
0234 <span class="comment">%         5521      char       TextPatientAgeUnits       ('Y'=years)</span>
0235 <span class="comment">%         5529      char[7]    TextPatientPosition</span>
0236 <span class="comment">%         5541      char[5]    TextImageNumberFlag       ('IMAGE'=image)</span>
0237 <span class="comment">%         5546      char[3]    TextImageNumber</span>
0238 <span class="comment">%         5559      char[2]    TextDateDD</span>
0239 <span class="comment">%         5562      char[3]    TextDateMM</span>
0240 <span class="comment">%         5566      char[4]    TextDateYYYY</span>
0241 <span class="comment">%         5571      char[2]    TextTimeHH</span>
0242 <span class="comment">%         5574      char[2]    TextTimeMM</span>
0243 <span class="comment">%         5577      char[2]    TextAcquisitionTimeFlag   ('TA'=acquisition time)</span>
0244 <span class="comment">%         5583      char[2]    TextAcquisitionTimeMM</span>
0245 <span class="comment">%         5586      char[2]    TextAcquisitionTimeSS</span>
0246 <span class="comment">%         5601      char[4]    TextAnnotation</span>
0247 <span class="comment">%         5655      char[25]   TextOrganization</span>
0248 <span class="comment">%         5682      char[5]    TextStation</span>
0249 <span class="comment">%         5695      char[3]    TextAcquisitionMatrixPhase</span>
0250 <span class="comment">%         5698      char       TextAcquisitionMatrixPhaseAxis  ('h'=horizontal,' '=vertical)</span>
0251 <span class="comment">%         5700      char[3]    TextAcquisitionMatrixFreq</span>
0252 <span class="comment">%         5703      char       TextAcquisitionMatrixFreqO      ('o'=o,' '=blank)</span>
0253 <span class="comment">%         5704      char       TextAcquisitionMatrixFreqS      ('s'=s,' '=blank)</span>
0254 <span class="comment">%         5706      char[8]    TextSequence</span>
0255 <span class="comment">%         5714      char[3]    TextFlipAngle</span>
0256 <span class="comment">%         5718      char[4]    TextScanNumberFlag        ('SCAN'=scan)</span>
0257 <span class="comment">%         5723      char[3]    TextScanNumberA</span>
0258 <span class="comment">%         5726      char[3]    TextScanNumberB</span>
0259 <span class="comment">%         5730      char[2]    TextRepetitionTimeFlag    ('TR'=tr)</span>
0260 <span class="comment">%         5734      char[7]    TextRepetitionTime</span>
0261 <span class="comment">%         5742      char[2]    TextEchoTimeFlag          ('TE'=te)</span>
0262 <span class="comment">%         5746      char[5]    TextEchoTime</span>
0263 <span class="comment">%         5752      char       TextEchoNumber</span>
0264 <span class="comment">%         5790      char[2]    TextSliceThicknessFlag    ('SL'=slice thickness)</span>
0265 <span class="comment">%         5794      char[7]    TextSliceThickness</span>
0266 <span class="comment">%         5802      char[2]    TextSlicePositionFlag     ('SP'=slice position)</span>
0267 <span class="comment">%         5806      char[7]    TextSlicePosition</span>
0268 <span class="comment">%         5814      char[3]    TextAngleFlag1            ('Sag'=sagittal,'Cor'=coronal,'Tra'=transverse)</span>
0269 <span class="comment">%         5817      char       TextAngleFlag2            ('&gt;'=gt,'&lt;'=lt)</span>
0270 <span class="comment">%         5818      char[3]    TextAngleFlag3            ('Sag'=sagittal,'Cor'=coronal,'Tra'=transverse)</span>
0271 <span class="comment">%         5821      char[4]    TextAngle</span>
0272 <span class="comment">%         5838      char[3]    TextFOVFlag               ('FoV'=field of view)</span>
0273 <span class="comment">%         5842      char[3]    TextFOVH</span>
0274 <span class="comment">%         5846      char[3]    TextFOVV</span>
0275 <span class="comment">%         5874      char[2]    TextTablePositionFlag     ('TP'=table position)</span>
0276 <span class="comment">%         5878      char[7]    TextTablePosition</span>
0277 <span class="comment">%         5938      char[5]    TextStudyNumberFlag       ('STUDY'=study)</span>
0278 <span class="comment">%         5943      char[2]    TextStudyNumber</span>
0279 <span class="comment">%         5956      char[2]    TextDOBDD</span>
0280 <span class="comment">%         5959      char[3]    TextDOBMM</span>
0281 <span class="comment">%         5963      char[4]    TextDOBYYYY</span>
0282 <span class="comment">%         5992      char[3]    TextStudyNumberFlag2      ('STU'=study)</span>
0283 <span class="comment">%         5996      char[3]    TextImageNumberFlag2      ('IMA'=study)</span>
0284 <span class="comment">%         5999      char[2]    TextStudyNumber2</span>
0285 <span class="comment">%         6002      char[2]    TextImageNumber2</span>
0286 <span class="comment">%         6013      char[5]    TextStudyImageNumber3</span>
0287 <span class="comment">%         6031      char[15]   TextModelName</span>
0288 <span class="comment">%         6058      char[25]   TextPatientName</span>
0289 <span class="comment">%         6085      char[2]    TextScanStartTimeHH</span>
0290 <span class="comment">%         6088      char[2]    TextScanStartTimeMM</span>
0291 <span class="comment">%         6091      char[2]    TextScanStartTimeSS</span>
0292 
0293 
0294 <span class="comment">% I have asked David Clunie to update his FAQ to include information for</span>
0295 <span class="comment">% finding size of the gap given a Siemens Vision format image. For the record</span>
0296 <span class="comment">% the answer is:  GapRatio is at byte 4136 (measured from the beginning of</span>
0297 <span class="comment">% the file) and it is of type double (8-byte float, Big Endian). For example,</span>
0298 <span class="comment">% if SliceThickness is 8mm and the GapRatio is 0.25, then the Gap is</span>
0299 <span class="comment">% (8x0.25=) 2mm, and the interslice distance is (8+2=) 10mm.</span>
0300 
0301 
0302 <span class="keyword">return</span></pre></div>
<hr><address>Generated on Fri 21-May-2004 12:38:21 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> &copy; 2003</address>
</body>
</html>