<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 12/03/10 00:38, Donald Tournier wrote:
<blockquote
 cite="mid:325dd9ac1003111638w69a1c285if7561bf40683509d@mail.gmail.com"
 type="cite">
  <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  <style type="text/css">DIV.EVShortcutBanner {PADDING:10px;FONT-WEIGHT:bold;FONT-SIZE:10pt;FONT-FAMILY:verdana;BACKGROUND-COLOR:lightblue}DIV.EVAttachBanner {PADDING:3px;FONT-WEIGHT:bold;FONT-SIZE:10pt;FONT-FAMILY:verdana;BACKGROUND-COLOR:lightblue}DIV.EVShortcutMsgBody {PADDING:10px;FONT-SIZE:10pt;FONT-FAMILY:verdana}TD.EVShortcutAttName {FONT-SIZE:10pt;FONT-FAMILY:verdana}TD.EVShortcutAttSize {FONT-SIZE:7pt;FONT-FAMILY:verdana}</style>
  <div class="EVShortcutBanner">This message has been archived.  <a
 moz-do-not-send="true"
 href="https://kclev.kcl.ac.uk/EnterpriseVault/ViewMessage.asp?VaultId=1CFCE6957FBD5654F8CCD18F833B7A6301110000kcl-evsite01&amp;SavesetId=730000000000000%7E201003120038120000%7E0%7EE7E2A350A97B45AD91183E4EEB3B887">View

the original item </a></div>
  <div class="EVShortcutMsgBody">Hi Philip, Sorry, there is currently
no documentation on that file format... It is quite simple, though. It
consists of a text header, followed by binary data. The first line of
the text header</div>
  <div class="EVAttachBanner">Attachments:</div>
  <table border="0" cellpadding="2" cellspacing="1" width="97%"
 align="center">
    <tbody>
      <tr>
        <td class="EVShortcutAttName" width="80%"><a
 moz-do-not-send="true"
 href="https://kclev.kcl.ac.uk/EnterpriseVault/ViewMessage.asp?VaultId=1CFCE6957FBD5654F8CCD18F833B7A6301110000kcl-evsite01&amp;SavesetId=730000000000000%7E201003120038120000%7E0%7EE7E2A350A97B45AD91183E4EEB3B887&amp;AttachmentId=1">read_mrtrix_tracks.m
        </a></td>
        <td class="EVShortcutAttSize">(2 KB)</td>
      </tr>
    </tbody>
  </table>
</blockquote>
<br>
<br>
Hi<br>
<br>
am trying to process some of the Philips DICOM we acquired, in order to
perform CSD. At the estimate_response stage it crashes, something went
wrong with 'axes specification' but nore sure what. Here is the
sequence of operation I go though:<br>
<br>
    $ mrconvert -info DICOM/ ./mrtrix-data/dwi.mif<br>
    [...]<br>
   <br>
    $ mrconvert dwi.mif -coord 3 0 - | threshold - - | median3D - - |
median3D - mask.mif<br>
     [...]<br>
   <br>
    $ dwi2tensor dwi.mif dti.mif<br>
    [...]<br>
    <br>
    $ tensor2FA dti.mif - | mrmult - mask.mif fa.mif<br>
    $ erode mask.mif - | erode - - | mrmult fa.mif - - | threshold -
-abs 0.6 sf.mif <br>
    $ estimate_response dwi.mif sf.mif  response.txt<br>
<br>
      estimate_response: incorrect number of axes in axes specification
"-0,-1,+2,+3"<br>
      *** glibc detected *** estimate_response: free(): invalid next
size (fast): 0x000000000194c420 ***<br>
      ======= Backtrace: =========<br>
     /lib/libc.so.6(+0x775b6)[0x7f4f4d5065b6]<br>
    /lib/libc.so.6(cfree+0x73)[0x7f4f4d50ce53]<br>
    etc...<br>
<br>
Am trying to diagnose if I'm doing some basic error, or if there's  a
problem with the way the axes are read from the DICOMS?<br>
Any suggestion?<br>
<br>
Ph<br>
</body>
</html>