Error running dwi2response dhollander

Hello, when I use the HCP dataset, the following error occurred when generating the response function through dhollander. Can you help me look at it.

root@zzc-virtual-machine://home/zzc/SC_Martix# dwi2response dhollander DWI.mif wm.txt gm.txt csf.txt
dwi2response: 
dwi2response: Note that this script makes use of commands / algorithms that have relevant articles for citation. Please consult the help page (-help option) for more information.
dwi2response: 
dwi2response: Generated scratch directory: /home/zzc/SC_Martix/dwi2response-tmp-16NMRP/
dwi2response: Importing DWI data (/home/zzc/SC_Martix/DWI.mif)...
dwi2response: Changing to scratch directory (/home/zzc/SC_Martix/dwi2response-tmp-16NMRP/)
dwi2response: Computing brain mask (dwi2mask)...
dwi2response: -------
dwi2response: 4 unique b-value(s) detected: 5,998,1998,2996 with 18,90,90,90 volumes
dwi2response: -------
dwi2response: Preparation:
dwi2response: * Eroding brain mask by 3 pass(es)...
dwi2response:   [ mask: 880236 -> 728255 ]
dwi2response: * Computing signal decay metric (SDM):
dwi2response:  * b=5...
dwi2response:  * b=998...
dwi2response:  * b=1998...
dwi2response:  * b=2996...
dwi2response: * Removing erroneous voxels from mask and correcting SDM...

dwi2response: [ERROR] Error trying to calculate statistics from image 'safe_mask.mif'

This is bval and bvec of DWI file. It seems that there are two problems. One is that there is no file with b=0, and the other is that there is no bvec vector of 0, 0, 0.

 0.54218611650394    0.672049144404883   -0.504365108403665                    5
  -0.918106027041497   -0.306174009013839   -0.251672007411375                 1000
   0.382921084996627   -0.843687187192569   -0.376249083496686                 1995
 -0.0979770807989744   -0.442573364995367    0.891363735090669                 3005
  -0.204410162696602    0.907238722184919   -0.367606292593889                  995
   0.420482043989621    0.790796082680479    0.444788046489021                 2995
   0.713198097813056    0.520708071409533   -0.469267064308591                 2005
  -0.337916186702255    0.287233158701917    0.896275495205981                  990
  -0.595529494166168   0.0984770817044055   -0.797274661554707                 1990
   -0.90288223833564    0.417035110116462    0.104333027504118                 3000
  -0.172680031702285   -0.874402160711568   -0.453434083405999                 1000
  -0.258874105604443   0.0773700315713278    0.962807392816524                 1985
  -0.845580276880818   -0.460928150889544     0.26933108819389                 2995
  -0.409899258109698     0.27768217490657   -0.868835547120556                 1005
  -0.976788698148888    0.205441936510282  -0.0606419812630351                 1995
  -0.293249274591162    0.947821887671434   -0.125054117096231                 2995
   0.542752291688461     0.67024836018575   -0.506149271989239                    5
  -0.902899827878606    0.288457944993165    0.318690939292449                  995
  -0.424887074412828   -0.882829154626655   -0.200209035106045                 2000
  -0.104117993199621    0.382692975098606    0.917989940196657                 3010
  -0.754053365600972   0.0812180393801047    -0.65177231600084                 3005
  -0.542995147214054   -0.609858165315784    0.577260156514941                  995
  -0.188334072100899      0.7961743050038   -0.575010220202745                 2005
  -0.791947431774366    0.542407295682443   -0.280381152890925                  995
  -0.304833998998809   -0.544402998197872   -0.781473997396946                 1990
   -0.69799928137964     0.10319904159699     0.70862328567933                 2985
  -0.438156189814333   -0.391789169812817   -0.809024350526466                 1005
  -0.793025835629173    0.251160947909239    0.555002884920417                 2000
  -0.307747173705584    -0.89438350481623     0.32460718320589                 3005
  -0.486607052789893    0.807463087683228    0.333492036193073                  995
  -0.867470080106658   -0.429637039703297   -0.250814023201925                 1995
  -0.516189404504821    0.601352471305617   -0.609855477905696                 2990
  -0.543239049701324   -0.668560061101629    0.507857046401237                    5
  -0.299073868499974    0.793319651199931    0.530281766799954                 3005
   -0.59307624999097   -0.803392338687767   0.0531160223891912                  995
  -0.619700396386662    0.753802482083775   -0.218525139795297                 2000
  -0.220479205006467   -0.289159268808481    0.931544866027322                  990
  -0.235762270090666   -0.520906596779376    0.820409939867518                 1990
  -0.844673202015499   -0.377372090206924   -0.379628090806966                 2990
  -0.900420863430002   -0.214813967407158    0.378281942612604                  995
  -0.284420181992201  -0.0382440244789514   -0.957936613073734                 2985
 -0.0317900183400865   -0.995056574002709   0.0940840542702561                 1995
  -0.571379341114204   -0.422917252410513    0.703325419817484                 3005
  -0.833420123387469   0.0974280144285351   -0.543984080591821                  995
  -0.244201914612837    0.562210803529554    0.790116723841534                 2005
 0.00376600041994212   -0.999348111484641  -0.0359050040094482                 1000
  -0.795861221597321    0.347206096698831    -0.49603713809833                 1990
  -0.999893301458346 -0.00990900298658721  -0.0107330032395529                 2995
  -0.543458310277987   -0.668341381572928    0.507910289979427                    5
 -0.0404320094397021    0.596294139195606    0.801747187194092                 1005
  -0.948777265248752   -0.197271055189344     0.24679106898667                 2005
 -0.0746470431175321   -0.975369563467753   -0.207562119893138                 2995
  -0.643577914809253    0.749577900810777    0.154726979502225                 3000
 -0.0580860279098485    0.661451317798274   -0.747735359298049                 1005
   -0.16639709079375    0.183086099893123   -0.968912528663606                 1985
  -0.665435329233247 -0.00505800250325271    0.746438369337294                 1005
  -0.448908958790446   -0.477353956189841   -0.755389930583924                 2985
  -0.442191314495712   -0.834387593491909    0.329035233996809                 2005
  -0.649559034414025   -0.635016033613711   -0.418124022109028                  995
    -0.6784017246993   -0.366664851199622   -0.636653741699343                 1990
 -0.0511820312104192    0.783971478006421   -0.618683377305067                 3005
  -0.582544002218093    0.810136003025162  -0.0657430002520419                 1000
  -0.585145879388032    0.476836901690247    0.655919864786585                 2990
  -0.750801237987625     0.60789319268998    0.258386081895741                 2000
  -0.543328464487833   -0.667284570385058    0.509436435488592                    5
  -0.131552009999449  -0.0654380049697259   -0.989147075095856                  990
  -0.601377068628482   -0.157044017907438    0.783379089437102                 2005
   -0.67644637810541   -0.736237411605888  -0.0193590108201548                 2995
  -0.992240602906261    0.104311958300658  -0.0676579729204269                 1000
  -0.069078023558802   -0.870552296884902    -0.48720316609155                 2005
  -0.835667252699957    0.422089127699978   -0.351427106299982                 3000
  -0.258221160010953   -0.876172543037165    0.406993252217264                  995
  -0.412554306188299     0.52357638858515   -0.745430553278857                 2005
  -0.913091205751194  -0.0270510060985541    0.406857091678253                 3005
  -0.614132944112913    0.523251952411003    0.590802946212423                 1005
   0.757527339081152    0.651709291683785    -0.03778001690906                 2000
  -0.267939044791979  -0.0508470084984779      0.9620931607712                 3010
  -0.572262235014347    0.542237222613594   -0.615219252615424                  990
  -0.128238094597249   -0.252472186294584   -0.959068707779427                 2005
  -0.347430140499695    0.937583379199178  -0.0151560061299867                 1995
  -0.541318275673902   -0.668358340367777    0.510168259775404                    5
  -0.847558373491777   -0.515673227194997    0.125403055298783                 1000
  -0.337305154801889    0.357488164002002   -0.870877399604877                 3010
  -0.542391447494371   -0.700061577492734     0.46446238319518                 2985
  -0.663109871071164   -0.149311970993507   -0.733478857368104                 1005
  -0.930698900385316  -0.0446329952192958   -0.363052961094272                 1995
 -0.0782209762476072   -0.720001781377975   -0.689549790578907                 3000
  -0.172307041197301    0.839073200686858    0.516010123391918                 1000
  -0.536620165900573    0.343118106100367    0.770914238300824                 2005
  -0.560124983480788      0.7784279770733    -0.28338999169028                 3000
  -0.526839094326192   -0.325222058216169    0.785284140639041                  990
  -0.177254095408779   -0.799087430239579    0.574491309328454                 2005
  -0.700595846988757   -0.245110946496067   -0.670138853589246                 2990
  -0.348816983401123   -0.936938955403017    -0.02172799897007                 3000
    -0.8283810510807    0.547281033787249    0.119450007397217                 1000
  -0.532177321318868   -0.677373408924016    0.507890306618007                    5
  -0.814988073270315    0.389781035085802    0.428795038584381                 2990
  -0.346746254408215   -0.927152680221967   -0.141968104203364                  995
   -0.94617988778803    0.173319979397807    0.273319967596542                 1995
  -0.517321784114036   -0.701339707319029   -0.490408795413306                 1995
  -0.112762077895147     0.35153024268487   -0.929360641659999                 1005
  -0.840123462042673    0.525170288826675   -0.135605074606888                 2000
  -0.298144139687845    0.926328434062235    0.230273107890612                 2995
  -0.337641993492075   -0.655881987484606   -0.675141987084154                 1005
   -0.14881402119752   -0.255703036395738    0.955233136084079                 2010
  -0.930077920522362   -0.361758969108698  -0.0639179945415368                 2995
  -0.472352013212713    0.776553021720901   -0.416952011711222                 1005
  -0.049598019790546    0.973122388310713    0.224884089702476                 2000
 -0.0273890081586775      0.1138940339945   -0.993115296052046                 3010
  -0.979063729935171   -0.178594950706416   0.0976629730635084                 1000
  -0.459848259306567   -0.167778094602396   -0.872003491612453                 2005
  -0.940229811180774   0.0824369834483143   -0.330411933693244                 2990
  -0.533880156404884   -0.676014198106185    0.507914148804647                    5
  -0.180477034392347   0.0496910094578929    0.982323186958345                  990
  -0.246591051406178   -0.674366140516895    0.696005145017437                 2985
  -0.806299190989417   -0.196798046597417    0.557810132092679                 2005
  -0.300615225590592    0.950508713170253    0.078509058907543                 1000
  -0.514723201596613    0.703482275595371   -0.490074191996775                 2005
  -0.409931147805028    0.258913093403176    0.874597315410727                 2985
  -0.742665352090999      0.2309981094972    0.628560297992382                  995
  -0.529633874894963    0.625927852194047    0.572452864794556                 2005
 -0.0422350188685365    0.952386425566999   -0.301954134889537                 2990
  -0.790286309879031   -0.341599133990936   -0.508682199486503                  995
  -0.468778186319099   -0.611190242924901    0.637725253525982                 2005
  -0.767804993463307   -0.246474997888221    0.591375994971738                 2990
  -0.227243106201486   -0.713119333304664    0.663190309904338                 1005
  -0.975122361018939   -0.211698078404112  -0.0657290243412766                 2000
  -0.208906105908866   -0.380788193016161   -0.900754456538229                 3005
  -0.535606357005268   -0.673278448706622    0.509727339705013                    5
  -0.678534153754295   -0.658244149155661    0.326046073878038                  995
  -0.682611038463893    -0.64164203616606   -0.349768019681499                 2995
  -0.036472990351478   -0.668433823127086   -0.742876803430103                 1990
  -0.336122002811862    0.588045004920752    0.735680006125962                 1005
  -0.912497235948654    0.268853069514335    -0.30832907971644                  995
  -0.125563079102896    0.540663340412472   -0.831815523719188                 2005
  -0.625957354490272    0.337120190894761   -0.703226398289071                 2990
  -0.342901997307079   -0.936119992719326   0.0780869993916121                 2000
  -0.774936555102518    0.625351448002032  -0.0917000656902979                 3000
  -0.202806045991185   -0.967807219457933    0.149060033793521                 1000
   -0.78848901270744  -0.0430550006904062   -0.613540009905789                 2005
  -0.221658080301966    0.616475223305467     -0.7555302737067                 2985
  -0.437429225392914  -0.0251720129695922   -0.898900463085438                  990
  -0.217475114998019    0.933837493891495   -0.283993150197413                 1995
  -0.243177062492935    0.615362158182122    0.749796192678216                 2985
  -0.535443201890281   -0.675215254487744    0.507331191290792                    5
  -0.109931042594911   -0.373968144882689   -0.920903356757372                  990
  -0.690661075325484   -0.678676074025042   -0.249772027209216                 2000
  -0.491969275293566   -0.181523101597626    0.851478476388864                 3010
  -0.740172893886437    -0.33999195129377    0.580128916889369                  995
  -0.517600996421408   0.0770259994731858    0.852147994135245                 1990
   -0.96894423924152    0.156768038706718    0.191235047208195                 2995
  -0.332028217996791    0.638833419493826   -0.694009455693293                 1005
  -0.590350316311876    0.795578426216004    0.136167073002739                 2000
  -0.543119233898825  -0.0384150165499169   -0.838776361298185                 3005
  -0.762914031144165   -0.625884025536232   -0.161961006609376                 1000
  -0.132806009797253    0.323149023893316    0.936983069280619                 1990
   -0.11533107769638   -0.881882594172318    -0.45714530798565                 2995
  -0.856434023133257   -0.470951012718288    0.211485005708212                 2000
  -0.704297283912938    0.624523251811473    0.337544136106201                  995
  -0.659114090112972   -0.717668098114124    0.224769030704424                 2990
   -0.53641853907571   -0.674530677869456    0.507211509677032                    5
   -0.68715000950305    0.271184003701204   -0.674006009302991                  995
   -0.58483185249482    0.685442827093929    0.433750890596158                 3005
  -0.630792266421301     0.39177816551323   -0.669784282922618                 1990
  -0.326996209098638    0.838801536396506   -0.435299278398187                 3005
  -0.953233772730456   0.0450279892614386    0.298860928709549                 1000
  -0.891991358413171    0.414531166506121    0.180320072402663                 2000
 -0.0565350136680607    0.944877228467588    0.322507077988937                 1000
  -0.295192027093887     -0.7822520718838   -0.548583050388639                 1990
  -0.616057135817788   -0.538810118815558   -0.574593126716591                 2990
 -0.0804870060289366   -0.995467074586848   0.0506670037993306                 2990
  -0.430746140191944  -0.0423420137792081    0.901479293483141                  990
  -0.424362267313138   -0.361697227911198      0.8301155229257                 2010
  -0.461756075099316   -0.807025131298804   -0.368092059899454                  995
  -0.103807047698811   -0.931956428289323     0.34739215969602                 2005
  -0.951770287611756   -0.238695072102948    0.192764058302381                 3000
  -0.535279233697402   -0.676437295396717    0.505874220897545                    5
  0.0163910142891597   -0.857938747756016    0.513490447573674                 1005
 -0.0187230040108281    0.144270030906381    0.989361211743761                 2980
  -0.218271992095099    0.772560972082653    0.596243978486612                 2005
  -0.954694974392554  -0.0502959986496077   -0.293304992097712                 1000
  -0.373026127804849   0.0943190323212262   -0.923014316311999                 1985
  -0.928780813679188    0.212599957295236   -0.303590939093197                 2005
  -0.921936325324869    0.386457136410425  -0.0261590092307056                 1000
  -0.329746262802453   -0.433583345503226    0.838613668306239                 2985
  -0.948229407608583     0.28506212250258   -0.140002060201267                 2995
  -0.165879077103063   -0.523905243409675    0.835468388215428                  990
  -0.634706192581732   -0.722912219379194    0.273031082792142                 1995
   -0.88491126655139   -0.141315042592237   -0.443804133675621                 3005
   -0.49145009328647    0.870418165276036   0.0291380055291978                 2995
  -0.372169159704515    0.915764392811109   -0.151214064901834                  995
  -0.532328190304126   -0.676463241805243    0.508944181903945                    5
  -0.865439164196562   0.0275480052298906    0.500256094898013                 1995
 -0.0644610183226942   -0.750342213231362   -0.657899186927498                 1005
  -0.499489089819598   -0.844067151833118   -0.195093035107655                 3000
  -0.743362112221737   -0.486362073414222   -0.459200069313428                 1990
  -0.545650923666897    0.279127960983066    0.790159889452064                  995
     0.0115049999003  0.00190999998304981    0.999931991126076                 2010
  -0.821346013594986    0.173852002898939    0.543292008996684                 2995
  -0.485191095172684   -0.823235161453652    0.294743057783406                  995
  -0.721736620102218    0.580299498501783   -0.377291324101159                 2005
  -0.183108042598608   -0.840140195393614    0.510525118696119                 2990
  -0.601801256689062   -0.475599202891356   -0.641592273688339                  995
  -0.186655049106081   -0.960631252931295   -0.205785054206704                 2000
  -0.275143915694631   -0.634136805687626     -0.7226107785859                 2985
  -0.229626216004976    0.147136138403188   -0.962092904920847                 1005
  -0.696849482085681    0.546618378188768   -0.464337321290459                 2990
  -0.712007199596499    0.482381135297628    0.510249143097491                 1995
  -0.531759196215453   -0.678061250219705    0.507410187314746                    5
  -0.149391109694924    0.311409228789419   -0.938459689368115                 3010
   -0.82765412805363  -0.0437320067675499    0.559532086568651                 1005
  0.0134330080997712    0.699289421688089   -0.714712430987827                 1990
  -0.658153651790673     0.69582068919014   -0.287519284795926                 1000
  -0.496585912084385   -0.391075930687702   -0.774894862775633                 2005
   -0.72430238549253   -0.484491257895003    0.490565261094941                 2990
  -0.416511076797789     0.72162113299617    0.552975101897065                  995
  -0.996057111937569  -0.0189580021307151   0.0866650097432688                 2000
  -0.134866982303014    0.895557882120017    0.424012944209477                 2995
  -0.907887416481526   -0.416074732591534   -0.051207967088958                 1000
  -0.428624269620105    0.866071544740623   -0.257296161812069                 2005
  -0.395566931393427    0.453313921392468    0.798769861386727                 2990
  -0.824912699369194   -0.557734796779171  -0.0919289664965669                 3000
  -0.309442087816951   -0.242589068813289   -0.919454260950368                 1005
  -0.640239235201521   -0.394739145000938    0.658995242101565                 1990
  -0.531366393394744   -0.679470503093279    0.505934374594995                    5
 -0.0476310407491453    0.980997839382397   -0.188081160896625                  995
  -0.504097224500187    0.175789078300065   -0.845567376500313                 3010
 -0.0706150246480543   -0.465479162487174   -0.882237307975691                 2005
  -0.765970952073046    0.416704973885337   -0.489535969382774                  995
  -0.270784138890875     0.94015348216832     0.20685110609303                 2000
  -0.783812994715697    0.554823996211111    0.278939998105586                 3000
  -0.489252879996745   0.0341839916097726    0.871471786194202                 2985
  -0.798514326483142   -0.466324190690155    0.380679155691963                  995
  -0.303403046181566    0.387049058876484   -0.870712132447099                 1990
  -0.472938162776739   -0.865574297857427    0.164653056691902                 2995
  -0.144390048307502     0.39322413142043    0.908034303547177                  995
  -0.567747010883857   -0.823179015776595 -0.00629600012082099                 2000
  -0.393167971718266   -0.535973961424901     0.74709494623471                  995
  -0.864794048289507   -0.242045013497063   -0.439938024594662                 2005
  -0.431671175400669   -0.254370103400394    -0.86542235170134                 3010
  -0.530896988106243   -0.679789984807994     0.50599798870595                    5
  -0.968537360976696   -0.145106054096509   -0.202187075395135                 2995
  -0.677158990476306    0.724826989774638    0.126812998195563                 1000
  -0.703782324882484   0.0470290217088295    0.708857327182358                 1990
   -0.83440842666783   -0.121585062195312   -0.537568274879274                  995
  -0.349045167510771    0.351375168610843    0.868736416826807                 1990
  -0.149397987801021    0.987395919406747    0.052244995740357                 3000
   -0.57718788318007   -0.797103838672477   -0.177424964093874                 1000
  -0.758778410426177    0.650348351822436    0.036088019521245                 2000
-0.00389600317994624   -0.629139513491318    0.777282634389274                 2985
  -0.285735097386808    0.486046165777561    -0.82590228166187                  990
  -0.374027132399387   -0.771789273298735    0.514242182099157                 1990
  -0.825886989778906    0.263202996793277   -0.498632993887264                 2990
  -0.779066824202405    0.406557908201255    0.477247892301473                 1000
-0.00841799881333982   -0.542350923521894   -0.840109881533914                 3005
  -0.770621020959692    0.172471004690979    -0.61351201666791                 1990
  -0.530226987199898   -0.681430983499869    0.504490987799903                    5
 -0.0559090488981408   -0.148206129595072    0.987374863667166                 1005
  -0.896094410971699   -0.443799708285984 -0.00752499505376234                 2000
  -0.347689025510691    0.727060053422357   -0.592027043518204                 2990
  -0.145146009294737   -0.960087061665187   -0.239093015391331                 1000
  -0.339880015484212   -0.156524007092729    0.927352042256923                 1990
  -0.289948956002616   -0.904983862508164   -0.311341952702809                 3000
  -0.970956938928738    0.202494987305993    0.127429992003772                 1000
  -0.570855840300793    0.726667796801009    0.382200893100531                 2000
  -0.696158872405727   -0.125273977001031    0.706872870505816                 2985
  -0.639837016808101   0.0750630019709503   -0.764836020109683                  995
  -0.333102138101926    0.689642286003987   -0.642990266603718                 1990
  -0.465677954707704    0.821652920113593    0.328679968005437                 3000
  -0.214555938201927   -0.544082843204888   -0.811134766307287                 1005
  -0.428309165987752   -0.648099251081467    0.629697243981993                 3005
  -0.469145190688719   -0.588413239185851   -0.658538267684165                 2005
  -0.529874112786429   -0.682837145382512    0.502958107087119                    5
  -0.269314248403491    0.794952733110306   -0.543617501307047                  995
  -0.156433122592337    0.986892773351658  -0.0396400310580583                 1995
  -0.876319662271815   -0.276220893491116    0.394671847887306                 3005
  -0.302263207200713    0.894873613402112    0.328387225100775                 1000
  -0.412825441777745   -0.771916826058387    0.483445517373938                  995
  -0.730055236098461    0.633929205098663   -0.255251082599462                 2990
  -0.221878992401386    0.163462994401021    0.961274967006006                 2985
  -0.848303010359843    -0.34876400428349    0.398429004881139                 2005
  -0.274947076710469  -0.0932330260035501   -0.956928266836438                 1990
  -0.803743729025936    -0.48623583611569   -0.342885884411065                 1000
  -0.742515057023319   -0.388279029812194   -0.545812041917141                 2990
   -0.29900595731035   -0.877100874830361   -0.375884946313012                 2000
  -0.651566104005555   -0.211416033801802    0.728536116306211                 1005
  -0.987863008700787  0.00887300007800707   -0.155074001400123                 1995
  -0.462363280023938    0.626741379532448    0.627228379832473                 3000

This is the DWI information displayed by the dirstat command. It can be seen that the b value information displayed by dirstat is inconsistent with the b value information provided by bval.

root@zzc-virtual-machine://home/zzc/SC_Martix# dirstat DWI.mif 
DWI.mif (b=0) [ 18 volumes ]

DWI.mif (b=998.11111111111109) [ 90 directions ]

  Bipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 13.1877, range [ 10.2194 - 16.6271 ]
    energy: total = 7430.09, mean = 165.113, range [ 162.787 - 167.897 ]

  Unipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 13.3211, range [ 10.2194 - 16.9701 ]
    energy: total = 4514.97, mean = 100.333, range [ 79.8884 - 113.111 ]

  Spherical Harmonic fit:
    condition numbers for lmax = 2 -> 10: [ 1.02808 1.11594 1.25888 1.41302 2.36592 ]

  Asymmetry of sampling:
    norm of mean direction vector = 0.498569
    WARNING: sampling is strongly asymmetric - this may affect resiliance to eddy-current distortions

DWI.mif (b=1998.2777777777778) [ 90 directions ]

  Bipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 13.133, range [ 11.0274 - 16.6268 ]
    energy: total = 7433.84, mean = 165.196, range [ 162.217 - 167.931 ]

  Unipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 13.6647, range [ 11.0274 - 26.9589 ]
    energy: total = 4413.53, mean = 98.0784, range [ 57.76 - 113.806 ]

  Spherical Harmonic fit:
    condition numbers for lmax = 2 -> 10: [ 1.04592 1.12854 1.27451 1.49217 2.19035 ]

  Asymmetry of sampling:
    norm of mean direction vector = 0.464999
    WARNING: sampling is strongly asymmetric - this may affect resiliance to eddy-current distortions

DWI.mif (b=2996) [ 90 directions ]

  Bipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 12.8808, range [ 9.74136 - 16.2665 ]
    energy: total = 7435.1, mean = 165.224, range [ 162.071 - 168.576 ]

  Unipolar electrostatic repulsion model:
    nearest-neighbour angles: mean = 13.2226, range [ 9.74136 - 32.8503 ]
    energy: total = 4489.95, mean = 99.7766, range [ 64.4733 - 112.631 ]

  Spherical Harmonic fit:
    condition numbers for lmax = 2 -> 10: [ 1.04974 1.14677 1.26432 1.48579 2.11973 ]

  Asymmetry of sampling:
    norm of mean direction vector = 0.489402
    WARNING: sampling is strongly asymmetric - this may affect resiliance to eddy-current distortions

Hi @zzcpage,

Thanks for all the debugging information – always useful. I’ve moved your posts to a fresh topic since the original one is 3 years old… The problem and causes are likely to be very different.

The gradient table is unlikely to be a problem here, nothing seems out of the ordinary. I’ll clear up the couple of issues you’ve mentioned below, but I don’t see that as being problematic.

To figure this out, can you try again, running with the -debug and the -nocleanup options, so we can get a clearer picture of which step is failing?


Regarding the gradient table:

It’s not uncommon, since the imaging gradients will themselves contribute to the diffusion weighting, and some manufacturers will take this into account when computing the actual DW directions and b-values. You will therefore generally end up with non-zero b=0 volumes, and slightly different gradients & b-values than you may have specified – you’ll note for example that your b-values are not exactly 1000, 2000, or 3000. These are typically minor variations though, and shouldn’t have any noticeable impact on the end result.

See above, and refer to our documentation page on how we handle the DW gradient table. This will hopefully clear up any confusion. In particular, MRtrix will (by default) consider that any volume with b<10 as a b=0 volume, so that data such as yours can be processed as expected.


Also: I note you’re running as root, which I wouldn’t recommend as a general rule… You’re in a VM, but I still get nervous seeing user commands running as super-user! Up to you though, it’s doubt it’s going anything to do with the problems you’re seeing here.

All the best,
Donald.

Thank you for your help. When I run the

dwi2Response dhollander DWI.mif wm.txt gm.txt csf.txt - noclean - debug

command, the following error prompt appears:

dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:145): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'safe_mask.mif', 'full_sdm.mif', '0', '-if', '10', '-min', 'safe_sdm.mif']]
Command:  mrcalc safe_mask.mif full_sdm.mif 0 -if 10 -min safe_sdm.mif
          mrcalc: [INFO] opening image "safe_mask.mif"...
          mrcalc: [INFO] image "safe_mask.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "full_sdm.mif"...
          mrcalc: [INFO] image "full_sdm.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "10"...
          mrcalc: [INFO] creating image "safe_sdm.mif"...
          mrcalc: [INFO] image "safe_sdm.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: min ((safe_mask.mif ? full_sdm.mif : 0), 10)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrstats -> mrstats
dwi2response: [DEBUG] run.version_match() (from image.py:233): Version-matched executable for mrstats: /home/zzc/mrtrix3/bin/mrstats
dwi2response: [DEBUG] run.exe_name() (from image.py:233): /home/zzc/mrtrix3/bin/mrstats -> /home/zzc/mrtrix3/bin/mrstats
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrstats safe_mask.mif -output mean -output median -output std -output std_rv -output min -output max -output count -mask safe_mask.mif' (piping data to local storage)

dwi2response: [ERROR] Error trying to calculate statistics from image 'safe_mask.mif'

All input logs of the command are as follows:

dwi2response: 
dwi2response: Note that this script makes use of commands / algorithms that have relevant articles for citation. Please consult the help page (-help option) for more information.
dwi2response: 
dwi2response: [DEBUG] path.script_subdir_name() (from algorithm.py:69): dwi2response
dwi2response: [DEBUG] path.from_user() (from dwi2response:83): DWI.mif -> /home/zzc/SC_Martix/DWI.mif
dwi2response: [DEBUG] path.name_temporary() (from image.py:36): /home/zzc/SC_Martix/mrtrix-tmp-ZmYOTG.json
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrinfo -> mrinfo
dwi2response: [DEBUG] run.version_match() (from image.py:37): Version-matched executable for mrinfo: /home/zzc/mrtrix3/bin/mrinfo
dwi2response: [DEBUG] run.exe_name() (from image.py:37): /home/zzc/mrtrix3/bin/mrinfo -> /home/zzc/mrtrix3/bin/mrinfo
dwi2response: Loading header for image file '/home/zzc/SC_Martix/DWI.mif'
dwi2response: [DEBUG] image.__init__() (from dwi2response:83): ['/home/zzc/mrtrix3/bin/mrinfo', '/home/zzc/SC_Martix/DWI.mif', '-json_all', '/home/zzc/SC_Martix/mrtrix-tmp-ZmYOTG.json']
dwi2response: [DEBUG] image.__init__() (from dwi2response:83): {'_name': '/home/zzc/SC_Martix/DWI.mif', '_size': [145, 174, 145, 288], '_spacing': [1.25, 1.25, 1.25, 1.0], '_strides': [-2, 3, 4, 1], '_format': 'MRtrix', '_datatype': 'Float32LE', '_intensity_offset': 0.0, '_intensity_scale': 1.0, '_transform': [[1.0, 0.0, 0.0, -90.0], [-0.0, 1.0, 0.0, -126.0], [-0.0, 0.0, 1.0, -72.0], [0.0, 0.0, 0.0, 1.0]], '_keyval': {'command_history': "mrconvert data.nii.gz DWI.mif -fslgrad bvecs bvals -datatype float32 -strides '0,0,0,1'  (version=3.0.3-120-gbb1c2fd4)", 'comments': 'FSL5.0', 'dw_scheme': [[0.5421861165, 0.6720491444, -0.5043651084, 5.0], [-0.918106027, -0.306174009, -0.2516720074, 1000.0], [0.382921085, -0.8436871872, -0.3762490835, 1995.0], [-0.0979770808, -0.442573365, 0.8913637351, 3005.0], [-0.2044101627, 0.9072387222, -0.3676062926, 995.0], [0.420482044, 0.7907960827, 0.4447880465, 2995.0], [0.7131980978, 0.5207080714, -0.4692670643, 2005.0], [-0.3379161867, 0.2872331587, 0.8962754952, 990.0], [-0.5955294942, 0.09847708171, -0.7972746616, 1990.0], [-0.9028822383, 0.4170351101, 0.1043330275, 3000.0], [-0.1726800317, -0.8744021607, -0.4534340834, 1000.0], [-0.2588741056, 0.07737003157, 0.9628073928, 1985.0], [-0.8455802769, -0.4609281509, 0.2693310882, 2995.0], [-0.4098992581, 0.2776821749, -0.8688355471, 1005.0], [-0.9767886981, 0.2054419365, -0.06064198126, 1995.0], [-0.2932492746, 0.9478218877, -0.1250541171, 2995.0], [0.5427522917, 0.6702483602, -0.506149272, 5.0], [-0.9028998279, 0.288457945, 0.3186909393, 995.0], [-0.4248870744, -0.8828291546, -0.2002090351, 2000.0], [-0.1041179932, 0.3826929751, 0.9179899402, 3010.0], [-0.7540533656, 0.08121803938, -0.651772316, 3005.0], [-0.5429951472, -0.6098581653, 0.5772601565, 995.0], [-0.1883340721, 0.796174305, -0.5750102202, 2005.0], [-0.7919474318, 0.5424072957, -0.2803811529, 995.0], [-0.304833999, -0.5444029982, -0.7814739974, 1990.0], [-0.6979992814, 0.1031990416, 0.7086232857, 2985.0], [-0.4381561898, -0.3917891698, -0.8090243505, 1005.0], [-0.7930258356, 0.2511609479, 0.5550028849, 2000.0], [-0.3077471737, -0.8943835048, 0.3246071832, 3005.0], [-0.4866070528, 0.8074630877, 0.3334920362, 995.0], [-0.8674700801, -0.4296370397, -0.2508140232, 1995.0], [-0.5161894045, 0.6013524713, -0.6098554779, 2990.0], [-0.5432390497, -0.6685600611, 0.5078570464, 5.0], [-0.2990738685, 0.7933196512, 0.5302817668, 3005.0], [-0.59307625, -0.8033923387, 0.05311602239, 995.0], [-0.6197003964, 0.7538024821, -0.2185251398, 2000.0], [-0.220479205, -0.2891592688, 0.931544866, 990.0], [-0.2357622701, -0.5209065968, 0.8204099399, 1990.0], [-0.844673202, -0.3773720902, -0.3796280908, 2990.0], [-0.9004208634, -0.2148139674, 0.3782819426, 995.0], [-0.284420182, -0.03824402448, -0.9579366131, 2985.0], [-0.03179001834, -0.995056574, 0.09408405427, 1995.0], [-0.5713793411, -0.4229172524, 0.7033254198, 3005.0], [-0.8334201234, 0.09742801443, -0.5439840806, 995.0], [-0.2442019146, 0.5622108035, 0.7901167238, 2005.0], [0.00376600042, -0.9993481115, -0.03590500401, 1000.0], [-0.7958612216, 0.3472060967, -0.4960371381, 1990.0], [-0.9998933015, -0.009909002987, -0.01073300324, 2995.0], [-0.5434583103, -0.6683413816, 0.50791029, 5.0], [-0.04043200944, 0.5962941392, 0.8017471872, 1005.0], [-0.9487772653, -0.1972710552, 0.246791069, 2005.0], [-0.07464704312, -0.9753695635, -0.2075621199, 2995.0], [-0.6435779148, 0.7495779008, 0.1547269795, 3000.0], [-0.05808602791, 0.6614513178, -0.7477353593, 1005.0], [-0.1663970908, 0.1830860999, -0.9689125287, 1985.0], [-0.6654353292, -0.005058002503, 0.7464383693, 1005.0], [-0.4489089588, -0.4773539562, -0.7553899306, 2985.0], [-0.4421913145, -0.8343875935, 0.329035234, 2005.0], [-0.6495590344, -0.6350160336, -0.4181240221, 995.0], [-0.6784017247, -0.3666648512, -0.6366537417, 1990.0], [-0.05118203121, 0.783971478, -0.6186833773, 3005.0], [-0.5825440022, 0.810136003, -0.06574300025, 1000.0], [-0.5851458794, 0.4768369017, 0.6559198648, 2990.0], [-0.750801238, 0.6078931927, 0.2583860819, 2000.0], [-0.5433284645, -0.6672845704, 0.5094364355, 5.0], [-0.13155201, -0.06543800497, -0.9891470751, 990.0], [-0.6013770686, -0.1570440179, 0.7833790894, 2005.0], [-0.6764463781, -0.7362374116, -0.01935901082, 2995.0], [-0.9922406029, 0.1043119583, -0.06765797292, 1000.0], [-0.06907802356, -0.8705522969, -0.4872031661, 2005.0], [-0.8356672527, 0.4220891277, -0.3514271063, 3000.0], [-0.25822116, -0.876172543, 0.4069932522, 995.0], [-0.4125543062, 0.5235763886, -0.7454305533, 2005.0], [-0.9130912058, -0.0270510061, 0.4068570917, 3005.0], [-0.6141329441, 0.5232519524, 0.5908029462, 1005.0], [0.7575273391, 0.6517092917, -0.03778001691, 2000.0], [-0.2679390448, -0.0508470085, 0.9620931608, 3010.0], [-0.572262235, 0.5422372226, -0.6152192526, 990.0], [-0.1282380946, -0.2524721863, -0.9590687078, 2005.0], [-0.3474301405, 0.9375833792, -0.01515600613, 1995.0], [-0.5413182757, -0.6683583404, 0.5101682598, 5.0], [-0.8475583735, -0.5156732272, 0.1254030553, 1000.0], [-0.3373051548, 0.357488164, -0.8708773996, 3010.0], [-0.5423914475, -0.7000615775, 0.4644623832, 2985.0], [-0.6631098711, -0.149311971, -0.7334788574, 1005.0], [-0.9306989004, -0.04463299522, -0.3630529611, 1995.0], [-0.07822097625, -0.7200017814, -0.6895497906, 3000.0], [-0.1723070412, 0.8390732007, 0.5160101234, 1000.0], [-0.5366201659, 0.3431181061, 0.7709142383, 2005.0], [-0.5601249835, 0.7784279771, -0.2833899917, 3000.0], [-0.5268390943, -0.3252220582, 0.7852841406, 990.0], [-0.1772540954, -0.7990874302, 0.5744913093, 2005.0], [-0.700595847, -0.2451109465, -0.6701388536, 2990.0], [-0.3488169834, -0.9369389554, -0.02172799897, 3000.0], [-0.8283810511, 0.5472810338, 0.1194500074, 1000.0], [-0.5321773213, -0.6773734089, 0.5078903066, 5.0], [-0.8149880733, 0.3897810351, 0.4287950386, 2990.0], [-0.3467462544, -0.9271526802, -0.1419681042, 995.0], [-0.9461798878, 0.1733199794, 0.2733199676, 1995.0], [-0.5173217841, -0.7013397073, -0.4904087954, 1995.0], [-0.1127620779, 0.3515302427, -0.9293606417, 1005.0], [-0.840123462, 0.5251702888, -0.1356050746, 2000.0], [-0.2981441397, 0.9263284341, 0.2302731079, 2995.0], [-0.3376419935, -0.6558819875, -0.6751419871, 1005.0], [-0.1488140212, -0.2557030364, 0.9552331361, 2010.0], [-0.9300779205, -0.3617589691, -0.06391799454, 2995.0], [-0.4723520132, 0.7765530217, -0.4169520117, 1005.0], [-0.04959801979, 0.9731223883, 0.2248840897, 2000.0], [-0.02738900816, 0.113894034, -0.9931152961, 3010.0], [-0.9790637299, -0.1785949507, 0.09766297306, 1000.0], [-0.4598482593, -0.1677780946, -0.8720034916, 2005.0], [-0.9402298112, 0.08243698345, -0.3304119337, 2990.0], [-0.5338801564, -0.6760141981, 0.5079141488, 5.0], [-0.1804770344, 0.04969100946, 0.982323187, 990.0], [-0.2465910514, -0.6743661405, 0.696005145, 2985.0], [-0.806299191, -0.1967980466, 0.5578101321, 2005.0], [-0.3006152256, 0.9505087132, 0.07850905891, 1000.0], [-0.5147232016, 0.7034822756, -0.490074192, 2005.0], [-0.4099311478, 0.2589130934, 0.8745973154, 2985.0], [-0.7426653521, 0.2309981095, 0.628560298, 995.0], [-0.5296338749, 0.6259278522, 0.5724528648, 2005.0], [-0.04223501887, 0.9523864256, -0.3019541349, 2990.0], [-0.7902863099, -0.341599134, -0.5086821995, 995.0], [-0.4687781863, -0.6111902429, 0.6377252535, 2005.0], [-0.7678049935, -0.2464749979, 0.591375995, 2990.0], [-0.2272431062, -0.7131193333, 0.6631903099, 1005.0], [-0.975122361, -0.2116980784, -0.06572902434, 2000.0], [-0.2089061059, -0.380788193, -0.9007544565, 3005.0], [-0.535606357, -0.6732784487, 0.5097273397, 5.0], [-0.6785341538, -0.6582441492, 0.3260460739, 995.0], [-0.6826110385, -0.6416420362, -0.3497680197, 2995.0], [-0.03647299035, -0.6684338231, -0.7428768034, 1990.0], [-0.3361220028, 0.5880450049, 0.7356800061, 1005.0], [-0.9124972359, 0.2688530695, -0.3083290797, 995.0], [-0.1255630791, 0.5406633404, -0.8318155237, 2005.0], [-0.6259573545, 0.3371201909, -0.7032263983, 2990.0], [-0.3429019973, -0.9361199927, 0.07808699939, 2000.0], [-0.7749365551, 0.625351448, -0.09170006569, 3000.0], [-0.202806046, -0.9678072195, 0.1490600338, 1000.0], [-0.7884890127, -0.04305500069, -0.6135400099, 2005.0], [-0.2216580803, 0.6164752233, -0.7555302737, 2985.0], [-0.4374292254, -0.02517201297, -0.8989004631, 990.0], [-0.217475115, 0.9338374939, -0.2839931502, 1995.0], [-0.2431770625, 0.6153621582, 0.7497961927, 2985.0], [-0.5354432019, -0.6752152545, 0.5073311913, 5.0], [-0.1099310426, -0.3739681449, -0.9209033568, 990.0], [-0.6906610753, -0.678676074, -0.2497720272, 2000.0], [-0.4919692753, -0.1815231016, 0.8514784764, 3010.0], [-0.7401728939, -0.3399919513, 0.5801289169, 995.0], [-0.5176009964, 0.07702599947, 0.8521479941, 1990.0], [-0.9689442392, 0.1567680387, 0.1912350472, 2995.0], [-0.332028218, 0.6388334195, -0.6940094557, 1005.0], [-0.5903503163, 0.7955784262, 0.136167073, 2000.0], [-0.5431192339, -0.03841501655, -0.8387763613, 3005.0], [-0.7629140311, -0.6258840255, -0.1619610066, 1000.0], [-0.1328060098, 0.3231490239, 0.9369830693, 1990.0], [-0.1153310777, -0.8818825942, -0.457145308, 2995.0], [-0.8564340231, -0.4709510127, 0.2114850057, 2000.0], [-0.7042972839, 0.6245232518, 0.3375441361, 995.0], [-0.6591140901, -0.7176680981, 0.2247690307, 2990.0], [-0.5364185391, -0.6745306779, 0.5072115097, 5.0], [-0.6871500095, 0.2711840037, -0.6740060093, 995.0], [-0.5848318525, 0.6854428271, 0.4337508906, 3005.0], [-0.6307922664, 0.3917781655, -0.6697842829, 1990.0], [-0.3269962091, 0.8388015364, -0.4352992784, 3005.0], [-0.9532337727, 0.04502798926, 0.2988609287, 1000.0], [-0.8919913584, 0.4145311665, 0.1803200724, 2000.0], [-0.05653501367, 0.9448772285, 0.322507078, 1000.0], [-0.2951920271, -0.7822520719, -0.5485830504, 1990.0], [-0.6160571358, -0.5388101188, -0.5745931267, 2990.0], [-0.08048700603, -0.9954670746, 0.0506670038, 2990.0], [-0.4307461402, -0.04234201378, 0.9014792935, 990.0], [-0.4243622673, -0.3616972279, 0.8301155229, 2010.0], [-0.4617560751, -0.8070251313, -0.3680920599, 995.0], [-0.1038070477, -0.9319564283, 0.3473921597, 2005.0], [-0.9517702876, -0.2386950721, 0.1927640583, 3000.0], [-0.5352792337, -0.6764372954, 0.5058742209, 5.0], [0.01639101429, -0.8579387478, 0.5134904476, 1005.0], [-0.01872300401, 0.1442700309, 0.9893612117, 2980.0], [-0.2182719921, 0.7725609721, 0.5962439785, 2005.0], [-0.9546949744, -0.05029599865, -0.2933049921, 1000.0], [-0.3730261278, 0.09431903232, -0.9230143163, 1985.0], [-0.9287808137, 0.2125999573, -0.3035909391, 2005.0], [-0.9219363253, 0.3864571364, -0.02615900923, 1000.0], [-0.3297462628, -0.4335833455, 0.8386136683, 2985.0], [-0.9482294076, 0.2850621225, -0.1400020602, 2995.0], [-0.1658790771, -0.5239052434, 0.8354683882, 990.0], [-0.6347061926, -0.7229122194, 0.2730310828, 1995.0], [-0.8849112666, -0.1413150426, -0.4438041337, 3005.0], [-0.4914500933, 0.8704181653, 0.02913800553, 2995.0], [-0.3721691597, 0.9157643928, -0.1512140649, 995.0], [-0.5323281903, -0.6764632418, 0.5089441819, 5.0], [-0.8654391642, 0.02754800523, 0.5002560949, 1995.0], [-0.06446101832, -0.7503422132, -0.6578991869, 1005.0], [-0.4994890898, -0.8440671518, -0.1950930351, 3000.0], [-0.7433621122, -0.4863620734, -0.4592000693, 1990.0], [-0.5456509237, 0.279127961, 0.7901598895, 995.0], [0.0115049999, 0.001909999983, 0.9999319911, 2010.0], [-0.8213460136, 0.1738520029, 0.543292009, 2995.0], [-0.4851910952, -0.8232351615, 0.2947430578, 995.0], [-0.7217366201, 0.5802994985, -0.3772913241, 2005.0], [-0.1831080426, -0.8401401954, 0.5105251187, 2990.0], [-0.6018012567, -0.4755992029, -0.6415922737, 995.0], [-0.1866550491, -0.9606312529, -0.2057850542, 2000.0], [-0.2751439157, -0.6341368057, -0.7226107786, 2985.0], [-0.229626216, 0.1471361384, -0.9620929049, 1005.0], [-0.6968494821, 0.5466183782, -0.4643373213, 2990.0], [-0.7120071996, 0.4823811353, 0.5102491431, 1995.0], [-0.5317591962, -0.6780612502, 0.5074101873, 5.0], [-0.1493911097, 0.3114092288, -0.9384596894, 3010.0], [-0.8276541281, -0.04373200677, 0.5595320866, 1005.0], [0.0134330081, 0.6992894217, -0.714712431, 1990.0], [-0.6581536518, 0.6958206892, -0.2875192848, 1000.0], [-0.4965859121, -0.3910759307, -0.7748948628, 2005.0], [-0.7243023855, -0.4844912579, 0.4905652611, 2990.0], [-0.4165110768, 0.721621133, 0.5529751019, 995.0], [-0.9960571119, -0.01895800213, 0.08666500974, 2000.0], [-0.1348669823, 0.8955578821, 0.4240129442, 2995.0], [-0.9078874165, -0.4160747326, -0.05120796709, 1000.0], [-0.4286242696, 0.8660715447, -0.2572961618, 2005.0], [-0.3955669314, 0.4533139214, 0.7987698614, 2990.0], [-0.8249126994, -0.5577347968, -0.0919289665, 3000.0], [-0.3094420878, -0.2425890688, -0.9194542609, 1005.0], [-0.6402392352, -0.394739145, 0.6589952421, 1990.0], [-0.5313663934, -0.6794705031, 0.5059343746, 5.0], [-0.04763104075, 0.9809978394, -0.1880811609, 995.0], [-0.5040972245, 0.1757890783, -0.8455673765, 3010.0], [-0.07061502465, -0.4654791625, -0.882237308, 2005.0], [-0.7659709521, 0.4167049739, -0.4895359694, 995.0], [-0.2707841389, 0.9401534822, 0.2068511061, 2000.0], [-0.7838129947, 0.5548239962, 0.2789399981, 3000.0], [-0.48925288, 0.03418399161, 0.8714717862, 2985.0], [-0.7985143265, -0.4663241907, 0.3806791557, 995.0], [-0.3034030462, 0.3870490589, -0.8707121325, 1990.0], [-0.4729381628, -0.8655742979, 0.1646530567, 2995.0], [-0.1443900483, 0.3932241314, 0.9080343035, 995.0], [-0.5677470109, -0.8231790158, -0.006296000121, 2000.0], [-0.3931679717, -0.5359739614, 0.7470949462, 995.0], [-0.8647940483, -0.2420450135, -0.4399380246, 2005.0], [-0.4316711754, -0.2543701034, -0.8654223517, 3010.0], [-0.5308969881, -0.6797899848, 0.5059979887, 5.0], [-0.968537361, -0.1451060541, -0.2021870754, 2995.0], [-0.6771589905, 0.7248269898, 0.1268129982, 1000.0], [-0.7037823249, 0.04702902171, 0.7088573272, 1990.0], [-0.8344084267, -0.1215850622, -0.5375682749, 995.0], [-0.3490451675, 0.3513751686, 0.8687364168, 1990.0], [-0.1493979878, 0.9873959194, 0.05224499574, 3000.0], [-0.5771878832, -0.7971038387, -0.1774249641, 1000.0], [-0.7587784104, 0.6503483518, 0.03608801952, 2000.0], [-0.00389600318, -0.6291395135, 0.7772826344, 2985.0], [-0.2857350974, 0.4860461658, -0.8259022817, 990.0], [-0.3740271324, -0.7717892733, 0.5142421821, 1990.0], [-0.8258869898, 0.2632029968, -0.4986329939, 2990.0], [-0.7790668242, 0.4065579082, 0.4772478923, 1000.0], [-0.008417998813, -0.5423509235, -0.8401098815, 3005.0], [-0.770621021, 0.1724710047, -0.6135120167, 1990.0], [-0.5302269872, -0.6814309835, 0.5044909878, 5.0], [-0.0559090489, -0.1482061296, 0.9873748637, 1005.0], [-0.896094411, -0.4437997083, -0.007524995054, 2000.0], [-0.3476890255, 0.7270600534, -0.5920270435, 2990.0], [-0.1451460093, -0.9600870617, -0.2390930154, 1000.0], [-0.3398800155, -0.1565240071, 0.9273520423, 1990.0], [-0.289948956, -0.9049838625, -0.3113419527, 3000.0], [-0.9709569389, 0.2024949873, 0.127429992, 1000.0], [-0.5708558403, 0.7266677968, 0.3822008931, 2000.0], [-0.6961588724, -0.125273977, 0.7068728705, 2985.0], [-0.6398370168, 0.07506300197, -0.7648360201, 995.0], [-0.3331021381, 0.689642286, -0.6429902666, 1990.0], [-0.4656779547, 0.8216529201, 0.328679968, 3000.0], [-0.2145559382, -0.5440828432, -0.8111347663, 1005.0], [-0.428309166, -0.6480992511, 0.629697244, 3005.0], [-0.4691451907, -0.5884132392, -0.6585382677, 2005.0], [-0.5298741128, -0.6828371454, 0.5029581071, 5.0], [-0.2693142484, 0.7949527331, -0.5436175013, 995.0], [-0.1564331226, 0.9868927734, -0.03964003106, 1995.0], [-0.8763196623, -0.2762208935, 0.3946718479, 3005.0], [-0.3022632072, 0.8948736134, 0.3283872251, 1000.0], [-0.4128254418, -0.7719168261, 0.4834455174, 995.0], [-0.7300552361, 0.6339292051, -0.2552510826, 2990.0], [-0.2218789924, 0.1634629944, 0.961274967, 2985.0], [-0.8483030104, -0.3487640043, 0.3984290049, 2005.0], [-0.2749470767, -0.093233026, -0.9569282668, 1990.0], [-0.803743729, -0.4862358361, -0.3428858844, 1000.0], [-0.742515057, -0.3882790298, -0.5458120419, 2990.0], [-0.2990059573, -0.8771008748, -0.3758849463, 2000.0], [-0.651566104, -0.2114160338, 0.7285361163, 1005.0], [-0.9878630087, 0.008873000078, -0.1550740014, 1995.0], [-0.46236328, 0.6267413795, 0.6272283798, 3000.0]], 'mrtrix_version': '3.0.3-120-gbb1c2fd4'}}
dwi2response: Generated scratch directory: /home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/
dwi2response: [DEBUG] path.from_user() (from dwi2response:93): DWI.mif -> /home/zzc/SC_Martix/DWI.mif
dwi2response: Importing DWI data (/home/zzc/SC_Martix/DWI.mif)...
dwi2response: [DEBUG] path.from_user() (from dwi2response:94): DWI.mif -> /home/zzc/SC_Martix/DWI.mif
dwi2response: [DEBUG] path.to_scratch() (from dwi2response:94): dwi.mif -> /home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/dwi.mif
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrconvert -> mrconvert
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrconvert: /home/zzc/mrtrix3/bin/mrconvert
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrconvert": Not a text file
dwi2response: [DEBUG] run.command() (from dwi2response:94): To execute: [['/home/zzc/mrtrix3/bin/mrconvert', '/home/zzc/SC_Martix/DWI.mif', '/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/dwi.mif', '-strides', '0,0,0,1']]
Command:  mrconvert /home/zzc/SC_Martix/DWI.mif /home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/dwi.mif -strides 0,0,0,1
          mrconvert: [INFO] opening image "/home/zzc/SC_Martix/DWI.mif"...
          mrconvert: [INFO] image "/home/zzc/SC_Martix/DWI.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrconvert: [INFO] found 288x4 diffusion gradient table
          mrconvert: [INFO] opening image "0,0,0,1"...
          mrconvert: [INFO] creating image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/dwi.mif"...
          mrconvert: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/dwi.mif" created with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrconvert: [100%] copying from "/home/zzc/SC_Martix/DWI.mif" to "/home/zzc/...response-tmp-NX2EQ2/dwi.mif"
dwi2response: Changing to scratch directory (/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/)
dwi2response: Computing brain mask (dwi2mask)...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): dwi2mask -> dwi2mask
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for dwi2mask: /home/zzc/mrtrix3/bin/dwi2mask
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/dwi2mask": Not a text file
dwi2response: [DEBUG] run.command() (from dwi2response:113): To execute: [['/home/zzc/mrtrix3/bin/dwi2mask', 'dwi.mif', 'mask.mif']]
Command:  dwi2mask dwi.mif mask.mif
          dwi2mask: [INFO] opening image "dwi.mif"...
          dwi2mask: [INFO] image "dwi.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwi2mask: [INFO] found 288x4 diffusion gradient table
          dwi2mask: [INFO] Diffusion gradient encoding data clustered into 3 non-zero shells and 18 b=0 volumes
          dwi2mask: [done] computing dwi brain mask
          dwi2mask: [INFO] creating image "mask.mif"...
          dwi2mask: [INFO] image "mask.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          dwi2mask: [done] applying mask cleaning filter
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrstats -> mrstats
dwi2response: [DEBUG] run.version_match() (from image.py:233): Version-matched executable for mrstats: /home/zzc/mrtrix3/bin/mrstats
dwi2response: [DEBUG] run.exe_name() (from image.py:233): /home/zzc/mrtrix3/bin/mrstats -> /home/zzc/mrtrix3/bin/mrstats
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrstats mask.mif -output mean -output median -output std -output std_rv -output min -output max -output count -mask mask.mif' (piping data to local storage)
dwi2response: Result: ImageStatistics(mean=1.0, median=1.0, std=0.0, std_rv=0.0, min=1.0, max=1.0, count=880236)
dwi2response: -------
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrinfo -> mrinfo
dwi2response: [DEBUG] run.version_match() (from image.py:153): Version-matched executable for mrinfo: /home/zzc/mrtrix3/bin/mrinfo
dwi2response: [DEBUG] run.exe_name() (from image.py:153): /home/zzc/mrtrix3/bin/mrinfo -> /home/zzc/mrtrix3/bin/mrinfo
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrinfo dwi.mif -shell_bvalues' (piping data to local storage)
dwi2response: Result: 5 998.111 1998.28 2996
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrinfo -> mrinfo
dwi2response: [DEBUG] run.version_match() (from image.py:153): Version-matched executable for mrinfo: /home/zzc/mrtrix3/bin/mrinfo
dwi2response: [DEBUG] run.exe_name() (from image.py:153): /home/zzc/mrtrix3/bin/mrinfo -> /home/zzc/mrtrix3/bin/mrinfo
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrinfo dwi.mif -shell_sizes' (piping data to local storage)
dwi2response: Result: 18 90 90 90
dwi2response: 4 unique b-value(s) detected: 5,998,1998,2996 with 18,90,90,90 volumes
dwi2response: -------
dwi2response: Preparation:
dwi2response: * Eroding brain mask by 3 pass(es)...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): maskfilter -> maskfilter
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for maskfilter: /home/zzc/mrtrix3/bin/maskfilter
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/maskfilter": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:110): To execute: [['/home/zzc/mrtrix3/bin/maskfilter', 'mask.mif', 'erode', 'eroded_mask.mif', '-npass', '3']]
Command:  maskfilter mask.mif erode eroded_mask.mif -npass 3
          maskfilter: [INFO] opening image "mask.mif"...
          maskfilter: [INFO] image "mask.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          maskfilter: [INFO] creating image "eroded_mask.mif"...
          maskfilter: [INFO] image "eroded_mask.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          maskfilter: [100%] applying erode filter to image mask.mif
          dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrstats -> mrstats
dwi2response: [DEBUG] run.version_match() (from image.py:233): Version-matched executable for mrstats: /home/zzc/mrtrix3/bin/mrstats
dwi2response: [DEBUG] run.exe_name() (from image.py:233): /home/zzc/mrtrix3/bin/mrstats -> /home/zzc/mrtrix3/bin/mrstats
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrstats mask.mif -output mean -output median -output std -output std_rv -output min -output max -output count -mask mask.mif' (piping data to local storage)
dwi2response: Result: ImageStatistics(mean=1.0, median=1.0, std=0.0, std_rv=0.0, min=1.0, max=1.0, count=880236)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrstats -> mrstats
dwi2response: [DEBUG] run.version_match() (from image.py:233): Version-matched executable for mrstats: /home/zzc/mrtrix3/bin/mrstats
dwi2response: [DEBUG] run.exe_name() (from image.py:233): /home/zzc/mrtrix3/bin/mrstats -> /home/zzc/mrtrix3/bin/mrstats
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrstats eroded_mask.mif -output mean -output median -output std -output std_rv -output min -output max -output count -mask eroded_mask.mif' (piping data to local storage)
dwi2response: Result: ImageStatistics(mean=1.0, median=1.0, std=0.0, std_rv=0.0, min=1.0, max=1.0, count=728255)
dwi2response:   [ mask: 880236 -> 728255 ]
dwi2response: * Computing signal decay metric (SDM):
dwi2response:  * b=5...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): dwiextract -> dwiextract
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for dwiextract: /home/zzc/mrtrix3/bin/dwiextract
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/dwiextract": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrmath -> mrmath
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrmath: /home/zzc/mrtrix3/bin/mrmath
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrmath": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:127): To execute: [['/home/zzc/mrtrix3/bin/dwiextract', 'dwi.mif', '-shells', '5', '-'], ['/home/zzc/mrtrix3/bin/mrcalc', '-', '0', '-max', '-'], ['/home/zzc/mrtrix3/bin/mrmath', '-', 'mean', 'mean_b5.mif', '-axis', '3']]
Command:  dwiextract dwi.mif -shells 5 - | mrcalc - 0 -max - | mrmath - mean mean_b5.mif -axis 3
          dwiextract: [INFO] opening image "dwi.mif"...
          dwiextract: [INFO] image "dwi.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [INFO] found 288x4 diffusion gradient table
          dwiextract: [INFO] Diffusion gradient encoding data clustered into 3 non-zero shells and 18 b=0 volumes
          dwiextract: [INFO] creating image "-"...
          dwiextract: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-w8YaC9.mif" created with dimensions 145x174x145x18, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [100%] extracting volumes
                    mrcalc: [INFO] opening image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-w8YaC9.mif" opened with dimensions 145x174x145x18, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-HtrLdn.mif" created with dimensions 145x174x145x18, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [100%] computing: max (/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-w8YaC9.mif, 0)
                    mrmath: [INFO] opening image "-"...
          mrmath: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-HtrLdn.mif" opened with dimensions 145x174x145x18, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrmath: [INFO] creating image "mean_b5.mif"...
          mrmath: [INFO] image "mean_b5.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrmath: [100%] computing mean along axis 3...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:129): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b5.mif', '-finite', 'mean_b5.mif', '0', '-if', '0', '-le', 'err_b5.mif', '-datatype', 'bit']]
Command:  mrcalc mean_b5.mif -finite mean_b5.mif 0 -if 0 -le err_b5.mif -datatype bit
          mrcalc: [INFO] opening image "mean_b5.mif"...
          mrcalc: [INFO] image "mean_b5.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "err_b5.mif"...
          mrcalc: [INFO] image "err_b5.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: ((finite (mean_b5.mif) ? mean_b5.mif : 0) <= 0)
dwi2response:  * b=998...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): dwiextract -> dwiextract
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for dwiextract: /home/zzc/mrtrix3/bin/dwiextract
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/dwiextract": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrmath -> mrmath
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrmath: /home/zzc/mrtrix3/bin/mrmath
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrmath": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:127): To execute: [['/home/zzc/mrtrix3/bin/dwiextract', 'dwi.mif', '-shells', '998', '-'], ['/home/zzc/mrtrix3/bin/mrcalc', '-', '0', '-max', '-'], ['/home/zzc/mrtrix3/bin/mrmath', '-', 'mean', 'mean_b998.mif', '-axis', '3']]
Command:  dwiextract dwi.mif -shells 998 - | mrcalc - 0 -max - | mrmath - mean mean_b998.mif -axis 3
          dwiextract: [INFO] opening image "dwi.mif"...
          dwiextract: [INFO] image "dwi.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [INFO] found 288x4 diffusion gradient table
          dwiextract: [INFO] Diffusion gradient encoding data clustered into 3 non-zero shells and 18 b=0 volumes
          dwiextract: [INFO] creating image "-"...
          dwiextract: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-vbeZH8.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [100%] extracting volumes
          mrcalc: [INFO] opening image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-vbeZH8.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-dV4aEj.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [100%] computing: max (/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-vbeZH8.mif, 0)
                    mrmath: [INFO] opening image "-"...
          mrmath: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-dV4aEj.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrmath: [INFO] creating image "mean_b998.mif"...
          mrmath: [INFO] image "mean_b998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrmath: [100%] computing mean along axis 3...
          dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:129): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b998.mif', '-finite', 'mean_b998.mif', '0', '-if', '0', '-le', 'err_b998.mif', '-datatype', 'bit']]
Command:  mrcalc mean_b998.mif -finite mean_b998.mif 0 -if 0 -le err_b998.mif -datatype bit
          mrcalc: [INFO] opening image "mean_b998.mif"...
          mrcalc: [INFO] image "mean_b998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "err_b998.mif"...
          mrcalc: [INFO] image "err_b998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: ((finite (mean_b998.mif) ? mean_b998.mif : 0) <= 0)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:134): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b5.mif', 'mean_b998.mif', '-divide', '-log', 'sdm_b998.mif']]
Command:  mrcalc mean_b5.mif mean_b998.mif -divide -log sdm_b998.mif
          mrcalc: [INFO] opening image "mean_b5.mif"...
          mrcalc: [INFO] image "mean_b5.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "mean_b998.mif"...
          mrcalc: [INFO] image "mean_b998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] creating image "sdm_b998.mif"...
          mrcalc: [INFO] image "sdm_b998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: log ((mean_b5.mif / mean_b998.mif))
dwi2response:  * b=1998...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): dwiextract -> dwiextract
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for dwiextract: /home/zzc/mrtrix3/bin/dwiextract
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/dwiextract": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrmath -> mrmath
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrmath: /home/zzc/mrtrix3/bin/mrmath
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrmath": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:127): To execute: [['/home/zzc/mrtrix3/bin/dwiextract', 'dwi.mif', '-shells', '1998', '-'], ['/home/zzc/mrtrix3/bin/mrcalc', '-', '0', '-max', '-'], ['/home/zzc/mrtrix3/bin/mrmath', '-', 'mean', 'mean_b1998.mif', '-axis', '3']]
Command:  dwiextract dwi.mif -shells 1998 - | mrcalc - 0 -max - | mrmath - mean mean_b1998.mif -axis 3
          dwiextract: [INFO] opening image "dwi.mif"...
          dwiextract: [INFO] image "dwi.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [INFO] found 288x4 diffusion gradient table
          dwiextract: [INFO] Diffusion gradient encoding data clustered into 3 non-zero shells and 18 b=0 volumes
          dwiextract: [INFO] creating image "-"...
          dwiextract: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-76GK89.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [100%] extracting volumes
                    mrcalc: [INFO] opening image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-76GK89.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-vWh1ol.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [100%] computing: max (/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-76GK89.mif, 0)
                    mrmath: [INFO] opening image "-"...
          mrmath: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-vWh1ol.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrmath: [INFO] creating image "mean_b1998.mif"...
          mrmath: [INFO] image "mean_b1998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrmath: [100%] computing mean along axis 3...
          dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:129): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b1998.mif', '-finite', 'mean_b1998.mif', '0', '-if', '0', '-le', 'err_b1998.mif', '-datatype', 'bit']]
Command:  mrcalc mean_b1998.mif -finite mean_b1998.mif 0 -if 0 -le err_b1998.mif -datatype bit
          mrcalc: [INFO] opening image "mean_b1998.mif"...
          mrcalc: [INFO] image "mean_b1998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "err_b1998.mif"...
          mrcalc: [INFO] image "err_b1998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: ((finite (mean_b1998.mif) ? mean_b1998.mif : 0) <= 0)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:134): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b5.mif', 'mean_b1998.mif', '-divide', '-log', 'sdm_b1998.mif']]
Command:  mrcalc mean_b5.mif mean_b1998.mif -divide -log sdm_b1998.mif
          mrcalc: [INFO] opening image "mean_b5.mif"...
          mrcalc: [INFO] image "mean_b5.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "mean_b1998.mif"...
          mrcalc: [INFO] image "mean_b1998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] creating image "sdm_b1998.mif"...
          mrcalc: [INFO] image "sdm_b1998.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: log ((mean_b5.mif / mean_b1998.mif))
dwi2response:  * b=2996...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): dwiextract -> dwiextract
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for dwiextract: /home/zzc/mrtrix3/bin/dwiextract
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/dwiextract": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrmath -> mrmath
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrmath: /home/zzc/mrtrix3/bin/mrmath
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrmath": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:127): To execute: [['/home/zzc/mrtrix3/bin/dwiextract', 'dwi.mif', '-shells', '2996', '-'], ['/home/zzc/mrtrix3/bin/mrcalc', '-', '0', '-max', '-'], ['/home/zzc/mrtrix3/bin/mrmath', '-', 'mean', 'mean_b2996.mif', '-axis', '3']]
Command:  dwiextract dwi.mif -shells 2996 - | mrcalc - 0 -max - | mrmath - mean mean_b2996.mif -axis 3
          dwiextract: [INFO] opening image "dwi.mif"...
          dwiextract: [INFO] image "dwi.mif" opened with dimensions 145x174x145x288, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [INFO] found 288x4 diffusion gradient table
          dwiextract: [INFO] Diffusion gradient encoding data clustered into 3 non-zero shells and 18 b=0 volumes
          dwiextract: [INFO] creating image "-"...
          dwiextract: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-SXZbBS.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          dwiextract: [100%] extracting volumes
          mrcalc: [INFO] opening image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-SXZbBS.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "-"...
          mrcalc: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-H2Gv9m.mif" created with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrcalc: [100%] computing: max (/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-SXZbBS.mif, 0)
                    mrmath: [INFO] opening image "-"...
          mrmath: [INFO] image "/home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/mrtrix-tmp-H2Gv9m.mif" opened with dimensions 145x174x145x90, voxel spacing 1.25x1.25x1.25x1, datatype Float32LE
          mrmath: [INFO] creating image "mean_b2996.mif"...
          mrmath: [INFO] image "mean_b2996.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrmath: [100%] computing mean along axis 3...
          dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:129): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b2996.mif', '-finite', 'mean_b2996.mif', '0', '-if', '0', '-le', 'err_b2996.mif', '-datatype', 'bit']]
Command:  mrcalc mean_b2996.mif -finite mean_b2996.mif 0 -if 0 -le err_b2996.mif -datatype bit
          mrcalc: [INFO] opening image "mean_b2996.mif"...
          mrcalc: [INFO] image "mean_b2996.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "err_b2996.mif"...
          mrcalc: [INFO] image "err_b2996.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: ((finite (mean_b2996.mif) ? mean_b2996.mif : 0) <= 0)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:134): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'mean_b5.mif', 'mean_b2996.mif', '-divide', '-log', 'sdm_b2996.mif']]
Command:  mrcalc mean_b5.mif mean_b2996.mif -divide -log sdm_b2996.mif
          mrcalc: [INFO] opening image "mean_b5.mif"...
          mrcalc: [INFO] image "mean_b5.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "mean_b2996.mif"...
          mrcalc: [INFO] image "mean_b2996.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] creating image "sdm_b2996.mif"...
          mrcalc: [INFO] image "sdm_b2996.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: log ((mean_b5.mif / mean_b2996.mif))
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:140): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'sdm_b998.mif', '90', '-mult', 'sdm_b1998.mif', '90', '-mult', '-add', 'sdm_b2996.mif', '90', '-mult', '-add', '270', '-divide', 'full_sdm.mif']]
Command:  mrcalc sdm_b998.mif 90 -mult sdm_b1998.mif 90 -mult -add sdm_b2996.mif 90 -mult -add 270 -divide full_sdm.mif
          mrcalc: [INFO] opening image "sdm_b998.mif"...
          mrcalc: [INFO] image "sdm_b998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "90"...
          mrcalc: [INFO] opening image "sdm_b1998.mif"...
          mrcalc: [INFO] image "sdm_b1998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "90"...
          mrcalc: [INFO] opening image "sdm_b2996.mif"...
          mrcalc: [INFO] image "sdm_b2996.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "90"...
          mrcalc: [INFO] opening image "270"...
          mrcalc: [INFO] creating image "full_sdm.mif"...
          mrcalc: [INFO] image "full_sdm.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: ((((sdm_b998.mif * 90) + (sdm_b1998.mif * 90)) + (sdm_b2996.mif * 90)) / 270)
dwi2response: * Removing erroneous voxels from mask and correcting SDM...
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:142): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'full_sdm.mif', '-finite', 'full_sdm.mif', '0', '-if', '0', '-le', 'err_sdm.mif', '-datatype', 'bit']]
Command:  mrcalc full_sdm.mif -finite full_sdm.mif 0 -if 0 -le err_sdm.mif -datatype bit
          mrcalc: [INFO] opening image "full_sdm.mif"...
          mrcalc: [INFO] image "full_sdm.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] creating image "err_sdm.mif"...
          mrcalc: [INFO] image "err_sdm.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: ((finite (full_sdm.mif) ? full_sdm.mif : 0) <= 0)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:144): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'err_b5.mif', 'err_b998.mif', '-add', 'err_b1998.mif', '-add', 'err_b2996.mif', '-add', 'err_sdm.mif', '-add', '0', 'eroded_mask.mif', '-if', 'safe_mask.mif', '-datatype', 'bit']]
Command:  mrcalc err_b5.mif err_b998.mif -add err_b1998.mif -add err_b2996.mif -add err_sdm.mif -add 0 eroded_mask.mif -if safe_mask.mif -datatype bit
          mrcalc: [INFO] opening image "err_b5.mif"...
          mrcalc: [INFO] image "err_b5.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "err_b998.mif"...
          mrcalc: [INFO] image "err_b998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "err_b1998.mif"...
          mrcalc: [INFO] image "err_b1998.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "err_b2996.mif"...
          mrcalc: [INFO] image "err_b2996.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "err_sdm.mif"...
          mrcalc: [INFO] image "err_sdm.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "eroded_mask.mif"...
          mrcalc: [INFO] image "eroded_mask.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] creating image "safe_mask.mif"...
          mrcalc: [INFO] image "safe_mask.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [100%] computing: (((((err_b5.mif + err_b998.mif) + err_b1998.mif) + err_b2996.mif) + err_sdm.mif) ? 0 : eroded_mask.mif)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrcalc -> mrcalc
dwi2response: [DEBUG] run.version_match() (from run.py:348): Version-matched executable for mrcalc: /home/zzc/mrtrix3/bin/mrcalc
dwi2response: [DEBUG] run._shebang() (from run.py:355): File "/home/zzc/mrtrix3/bin/mrcalc": Not a text file
dwi2response: [DEBUG] run.command() (from dhollander.py:145): To execute: [['/home/zzc/mrtrix3/bin/mrcalc', 'safe_mask.mif', 'full_sdm.mif', '0', '-if', '10', '-min', 'safe_sdm.mif']]
Command:  mrcalc safe_mask.mif full_sdm.mif 0 -if 10 -min safe_sdm.mif
          mrcalc: [INFO] opening image "safe_mask.mif"...
          mrcalc: [INFO] image "safe_mask.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Bit
          mrcalc: [INFO] opening image "full_sdm.mif"...
          mrcalc: [INFO] image "full_sdm.mif" opened with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [INFO] opening image "0"...
          mrcalc: [INFO] opening image "10"...
          mrcalc: [INFO] creating image "safe_sdm.mif"...
          mrcalc: [INFO] image "safe_sdm.mif" created with dimensions 145x174x145, voxel spacing 1.25x1.25x1.25, datatype Float32LE
          mrcalc: [100%] computing: min ((safe_mask.mif ? full_sdm.mif : 0), 10)
dwi2response: [DEBUG] run.exe_name() (from run.py:561): mrstats -> mrstats
dwi2response: [DEBUG] run.version_match() (from image.py:233): Version-matched executable for mrstats: /home/zzc/mrtrix3/bin/mrstats
dwi2response: [DEBUG] run.exe_name() (from image.py:233): /home/zzc/mrtrix3/bin/mrstats -> /home/zzc/mrtrix3/bin/mrstats
dwi2response: Command: '/home/zzc/mrtrix3/bin/mrstats safe_mask.mif -output mean -output median -output std -output std_rv -output min -output max -output count -mask safe_mask.mif' (piping data to local storage)

dwi2response: [ERROR] Error trying to calculate statistics from image 'safe_mask.mif'
dwi2response: Scratch directory retained; location: /home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/

OK, I’m not sure what’s going on, but my best guess is that the ‘safe’ mask it’s trying to generate ends up empty - no voxels at all. You can check with:

/home/zzc/mrtrix3/bin/mrstats /home/zzc/SC_Martix/dwi2response-tmp-NX2EQ2/safe_mask.mif

If that shows both min & max of zero, then that’s the problem.

As to what might be causing that, it’s hard to tell without more information about your data. Are these adult human in vivo data, acquired with a product sequence?

Thank you for your reply. My data is from dMRI data pre-processed by HCP in HCP 1200

As you said, safe_ Mask. mif voxel is empty

Thank you for your help. I think I may have found the answer to the question. I guess this may be caused by insufficient memory of my server. I extracted the first 32 volumes from the DWI.mif file to estimate the response function, and successfully ran the command dwi2response dhollander DWI.mif wm.txt gm.txt csf.txt. My server memory is 8GB, the size of DWI.mif before extraction is 4G, and the size of DWI.mif after extraction of the first 32 volumes is 447M, so I think the previous error may be caused by insufficient server memory. I want to know, if I run 4G DWI.mif, how much memory does the server need?

Ok, that makes sense. These are indeed huge datasets. I’m not sure how much would be needed exactly, but I would expect something like 4× the input (uncompressed) data should be sufficient – though you may find some of the subsequent steps will probably require substantially more than that…