{"version":3,"sources":["jsdelivr-header.js","/npm/nifti-js@1.0.1/nifti.js"],"names":[],"mappings":"AAAA;AACA;AACA;AACA;AACA;AACA;AACA,ACNA,aACA,IAAI,OAAS,QAAQ,UAEjB,iBAAmB,WACnB,IAAI,EAAM,IAAI,YAAY,GACtB,EAAS,IAAI,YAAY,GACzB,EAAU,IAAI,WAAW,GAE7B,OADA,EAAO,GAAK,SACI,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,GACpD,MACY,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,IAAqB,GAAZ,EAAQ,GAC3D,cAEX,QAAQ,KAAK,kCAEhB,CAZsB,GAgBvB,SAAS,iBAAiB,GACxB,IAAI,EAAO,IAAI,WAAW,GACtB,EAAS,EAAK,OACd,EAAO,IAAI,SAAS,GACxB,GAAI,EAAO,WAAW,IACpB,MAAM,IAAI,MAAM,uGAIlB,IAAI,GAAe,EACf,EAAM,IAAI,MAAM,GACpB,EAAI,GAAK,EAAK,SAAS,GAAI,IACvB,EAAE,EAAI,IAAM,EAAI,GAAG,KACrB,GAAgB,EAChB,EAAI,GAAK,EAAK,SAAS,GAAI,KAEzB,EAAE,EAAI,IAAM,EAAI,GAAG,IAErB,QAAQ,KAAK,kHAIf,IAAI,EAAa,EAAK,SAAS,EAAG,GAClC,GAAmB,MAAf,IAAuB,EAAE,EAAI,IAAM,EAAI,GAAG,IAK5C,GAHA,GAAgB,EAChB,EAAI,GAAK,EAAK,SAAS,GAAI,GAER,OADnB,EAAa,EAAK,SAAS,EAAG,IAE5B,MAAM,IAAI,MAAM,2FAEb,IAAI,EAAa,IACtB,MAAM,IAAI,MAAM,qEACQ,MAAf,GACT,QAAQ,KAAK,8GACf,CACA,IAAI,EAAQ,OAAO,aAAa,MAAM,KAAM,EAAK,SAAS,IAAK,MAC/D,GAAc,UAAV,GAA+B,UAAV,EACvB,MAAM,IAAI,MAAM,+FAIlB,IAAI,EAAW,EAAK,QAAQ,IAC5B,EAAI,OAAS,EAAE,KAAK,IAAI,EAAG,EAAI,IAC/B,IAAI,IAAI,EAAE,EAAG,EAAE,EAAI,OAAQ,IACzB,EAAI,GAAK,EAAK,SAAS,GAAG,EAAE,EAAG,GAC3B,EAAI,IAAI,IACV,QAAQ,KAAK,wCACb,EAAI,OAAS,GAGjB,GAAmB,IAAf,EAAI,OAAc,MAAM,IAAI,MAAM,wBAEtC,IAAI,EAAS,IAAI,MAAM,EAAI,QAC3B,IAAQ,EAAE,EAAG,EAAE,EAAO,OAAQ,IAC5B,EAAO,GAAK,EAAK,WAAW,GAAG,EAAE,EAAG,GAGtC,IAAI,EAAO,IAAI,aAAa,IAC5B,IAAQ,EAAE,EAAG,EAAE,GAAI,IACjB,EAAK,GAAK,EAAK,WAAW,IAAI,EAAE,EAAG,GA0DrC,MAtDa,CACX,aAAc,EAEd,WAAY,EACZ,SAAU,EACV,IAAK,EACL,UAAW,EAAK,WAAW,GAAI,GAC/B,UAAW,EAAK,WAAW,GAAI,GAC/B,UAAW,EAAK,WAAW,GAAI,GAC/B,YAAa,EAAK,SAAS,GAAI,GAE/B,SAAU,oBAAoB,EAAK,SAAS,GAAI,IAChD,OAAQ,EAAK,SAAS,GAAI,GAC1B,YAAa,EAAK,SAAS,GAAI,GAC/B,OAAQ,EACR,WAAY,EAAK,WAAW,IAAK,GAEjC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,SAAS,IAAK,GAC9B,WAAY,EAAK,QAAQ,KACzB,WAAY,iBAAiB,EAAK,QAAQ,MAC1C,QAAS,EAAK,WAAW,IAAK,GAC9B,QAAS,EAAK,WAAW,IAAK,GAC9B,eAAgB,EAAK,WAAW,IAAK,GACrC,QAAS,EAAK,WAAW,IAAK,GAE9B,QAAS,OAAO,aAAa,MAAM,KAAM,EAAK,SAAS,IAAK,MAC5D,SAAU,OAAO,aAAa,MAAM,KAAM,EAAK,SAAS,IAAK,MAE7D,WAAY,EAAK,SAAS,IAAK,GAC/B,WAAY,EAAK,SAAS,IAAK,GAE/B,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAChC,UAAW,EAAK,WAAW,IAAK,GAEhC,KAAM,EAEN,YAAa,OAAO,aAAa,MAAM,KAAM,EAAK,SAAS,IAAK,MAEhE,MAAO,EAEP,UAAW,EAAO,WAAa,IAAQ,CAAC,EAAE,EAAE,EAAE,GAAK,CAAC,EAAK,QAAQ,KAAM,EAAK,QAAQ,KAAM,EAAK,QAAQ,KAAM,EAAK,QAAQ,MAS9H,CAGA,SAAS,YAAY,GACnB,IAAI,EAAM,CAAC,EAQX,GAPA,EAAI,UAAY,EAAY,IAAI,GAChC,EAAI,KAAO,EAAY,SACvB,EAAI,SAAW,MACf,EAAI,OAAS,EAAY,aAAe,SAAW,MACnD,EAAI,MAAQ,EAAY,IAAI,MAAM,GAClC,EAAI,YAAY,EAAY,OAAO,MAAM,QAEV,IAA3B,EAAY,WAA0B,CAExC,IADA,EAAI,WAAa,EAAY,WACvB,EAAI,WAAW,OAAS,EAAI,WAChC,EAAI,WAAW,KAAK,IAEtB,EAAI,WAAW,OAAS,EAAI,SAC9B,CAEA,GAA+B,IAA3B,EAAY,WAAkB,CAEhC,IADA,EAAI,SAAW,EAAY,OAAO,MAAM,GAClC,EAAI,SAAS,OAAS,EAAI,WAC9B,EAAI,SAAS,KAAK,KAEpB,EAAI,eAAiB,KAAK,IAAI,EAAI,UAAW,EAC/C,MAAO,GAAI,EAAY,WAAa,EAAG,CAErC,EAAI,MAAQ,0BACZ,IAAI,EAAiC,IAA1B,EAAY,OAAO,GAAa,EAAI,EAAY,OAAO,GAC9D,EAAI,KAAK,KAAK,KAAK,IAAI,EAAI,GAAK,EAAY,UAAU,EAAY,UAAY,EAAY,UAAU,EAAY,UAAY,EAAY,UAAU,EAAY,aAC9J,EAAI,EAAY,UAChB,EAAI,EAAY,UAChB,EAAI,EAAY,UACpB,EAAI,gBAAkB,CACpB,CAAC,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,IAC9G,CAAC,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,IAC9G,CAAC,EAAK,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAK,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAK,EAAY,OAAO,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAC/H,EAAI,YAAc,CAAC,EAAY,UAAU,EAAY,UAAU,EAAY,UAC7E,MACE,QAAQ,KAAK,uBAAyB,EAAY,WAAa,wCAkBjE,OAfI,EAAY,WAAa,GAC3B,QAAQ,KAAK,+CAcR,CACT,CAqCA,SAAS,kBAAkB,EAAQ,EAAM,EAAK,GAC5C,IAAI,EAAG,EAAK,EAAM,EAAW,EAAoC,EAApB,EAAQ,WAAsC,SAAW,MACtG,IAAI,EAAE,EAAG,EAAE,EAAI,OAAQ,IACrB,GAAY,EAAI,GAElB,GAAY,SAAR,EAEF,OAAO,EAAO,MAAM,EAAE,EAAS,EAAQ,WAClC,GAAY,QAAR,GAA0B,SAAR,GAAmB,GAAc,iBAC5D,OAAO,GACP,IAAK,OAEH,OADA,EAAU,GACH,IAAI,UAAU,EAAO,MAAM,EAAE,IACtC,IAAK,QAEH,OADA,EAAU,GACH,IAAI,WAAW,EAAO,MAAM,EAAE,IACvC,IAAK,QAEH,OADA,EAAU,GACH,IAAI,WAAW,EAAO,MAAM,EAAW,EAAT,IACvC,IAAK,SAEH,OADA,EAAU,GACH,IAAI,YAAY,EAAO,MAAM,EAAW,EAAT,IACxC,IAAK,QAEH,OADA,EAAU,GACH,IAAI,WAAW,EAAO,MAAM,EAAW,EAAT,IACvC,IAAK,SAEH,OADA,EAAU,GACH,IAAI,YAAY,EAAO,MAAM,EAAW,EAAT,IAOxC,IAAK,QAEH,OADA,EAAU,GACH,IAAI,aAAa,EAAO,MAAM,EAAW,EAAT,IACzC,IAAK,SAEH,OADA,EAAU,GACH,IAAI,aAAa,EAAO,MAAM,EAAW,EAAT,IACzC,QAEE,YADA,QAAQ,KAAK,2BAA6B,QAK5C,OADA,EAAO,IAAI,SAAS,GACb,GACP,IAAK,OAEH,IADA,EAAM,IAAI,UAAU,GAChB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,QAAQ,GAExB,OAAO,EACT,IAAK,QAEH,IADA,EAAM,IAAI,WAAW,GACjB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,SAAS,GAEzB,OAAO,EACT,IAAK,QAEH,IADA,EAAM,IAAI,WAAW,GACjB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,SAAW,EAAF,GAEzB,OAAO,EACT,IAAK,SAEH,IADA,EAAM,IAAI,YAAY,GAClB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,UAAY,EAAF,GAE1B,OAAO,EACT,IAAK,QAEH,IADA,EAAM,IAAI,WAAW,GACjB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,SAAW,EAAF,GAEzB,OAAO,EACT,IAAK,SAEH,IADA,EAAM,IAAI,YAAY,GAClB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,UAAY,EAAF,GAE1B,OAAO,EAaT,IAAK,QAEH,IADA,EAAM,IAAI,aAAa,GACnB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,WAAa,EAAF,GAE3B,OAAO,EACT,IAAK,SAEH,IADA,EAAM,IAAI,aAAa,GACnB,EAAE,EAAG,EAAE,EAAU,IACnB,EAAI,GAAK,EAAK,WAAa,EAAF,GAE3B,OAAO,EACT,QAEE,YADA,QAAQ,KAAK,0BAA4B,GAI7C,SAAS,EAAU,GACjB,GAAI,EAAO,WAAW,EAAS,EAAY,MAAM,IAAI,MAAM,2CAC7D,CACF,CAEA,SAAS,oBAAoB,GAC3B,OAAO,GACP,KAAK,EACH,MAAO,MACT,KAAK,EACH,MAAO,QACT,KAAK,EACH,MAAO,QACT,KAAK,EACH,MAAO,QACT,KAAK,GACH,MAAO,QACT,KAAK,GACH,MAAO,YACT,KAAK,GACH,MAAO,SACT,KAAK,IACH,MAAO,QACT,KAAK,IACH,MAAO,OACT,KAAK,IACH,MAAO,SACT,KAAK,IACH,MAAO,SACT,KAAK,KACH,MAAO,QACT,KAAK,KACH,MAAO,SACT,KAAK,KACH,MAAO,WACT,KAAK,KACH,MAAO,aACT,KAAK,KACH,MAAO,aACT,KAAK,KACH,MAAO,SACT,QAEE,OADA,QAAQ,KAAK,iCAAmC,GACzC,EAEX,CAEA,SAAS,iBAAiB,GACxB,IAAI,EAAO,EACX,OAAe,EAAR,GACP,KAAK,EACH,EAAQ,GACR,MACF,KAAK,EACH,EAAQ,IACR,MACF,KAAK,EACH,EAAQ,KACR,MACF,KAAK,EACH,EAAQ,KACR,MACF,QACE,QAAQ,KAAK,6BAAqC,EAAN,IAC5C,EAAQ,GAEV,OAAe,GAAR,GACP,KAAK,EACH,EAAO,GACP,MACF,KAAK,EACH,EAAO,IACP,MACF,KAAK,GACH,EAAO,KACP,MACF,KAAK,GACH,EAAO,KACP,MACF,KAAK,GACH,EAAO,KACP,MACF,KAAK,GACH,EAAO,MACP,MACF,KAAK,GACH,EAAO,QACP,MACF,QACE,QAAQ,KAAK,6BAAqC,GAAN,IAC5C,EAAO,GAET,MAAkB,KAAV,GAAyB,KAAT,OAAe,EAAY,CAAC,EAAO,EAAO,EAAO,EAC3E,CAraA,OAAO,QAAQ,iBAAmB,iBAuLlC,OAAO,QAAQ,YAAc,SAAU,GAGrC,OADU,YADQ,iBAAiB,GAGrC,EAIA,OAAO,QAAQ,MAAQ,SAAU,GAC/B,IAAI,EAAc,iBAAiB,GAC/B,EAAM,YAAY,GAOtB,GALiC,IAA7B,EAAY,UAAU,IACxB,QAAQ,KAAK,yFAIW,UAAtB,EAAY,MAAmB,CACjC,IACI,EADO,IAAI,WAAW,GACR,OAClB,GAAI,EAAY,WAAW,KAAO,EAAY,WAAW,EAAO,WAC9D,MAAM,IAAI,MAAM,uBAElB,EAAI,OAAS,EAAO,MAAM,KAAK,MAAM,EAAY,aACpB,IAAzB,EAAY,WAEd,EAAI,KAAO,kBAAkB,EAAI,OAAQ,EAAY,SAAU,EAAY,IAAK,CAAC,WAAY,EAAY,eAE7G,CAEA,OAAO,CACT","file":"/npm/nifti-js@1.0.1/nifti.js","sourceRoot":"","sourcesContent":["/**\n * Minified by jsDelivr using Terser v5.37.0.\n * Original file: /npm/nifti-js@1.0.1/nifti.js\n *\n * Do NOT use SRI with dynamically generated files! More information: https://www.jsdelivr.com/using-sri-with-dynamic-files\n */\n","\"use strict\"\r\nvar assert = require('assert')\r\n\r\nvar systemEndianness = (function() {\r\n    var buf = new ArrayBuffer(4),\r\n        intArr = new Uint32Array(buf),\r\n        byteArr = new Uint8Array(buf)\r\n    intArr[0] = 0x01020304\r\n    if (byteArr[0]==1 && byteArr[1]==2 && byteArr[2]==3 && byteArr[3]==4) {\r\n        return 'big'\r\n    } else if (byteArr[0]==4 && byteArr[1]==3 && byteArr[2]==2 && byteArr[3]==1) {\r\n        return 'little'\r\n    }\r\n    console.warn(\"Unrecognized system endianness!\")\r\n    return undefined\r\n})()\r\n\r\n// Parses a NIfTI header\r\nmodule.exports.parseNIfTIHeader = parseNIfTIHeader\r\nfunction parseNIfTIHeader(buffer_org) {\r\n  var buf8 = new Uint8Array(buffer_org)\r\n  var buffer = buf8.buffer // Make sure we have an ArrayBuffer\r\n  var view = new DataView(buffer)\r\n  if (buffer.byteLength<348) {\r\n    throw new Error(\"The buffer is not even large enough to contain the minimal header I would expect from a NIfTI file!\")\r\n  }\r\n  \r\n  // First read dim[0], to determine byte order\r\n  var littleEndian = true\r\n  var dim = new Array(8)\r\n  dim[0] = view.getInt16(40, littleEndian)\r\n  if (1>dim[0] || dim[0]>7) {\r\n    littleEndian = !littleEndian\r\n    dim[0] = view.getInt16(40, littleEndian)\r\n  }\r\n  if (1>dim[0] || dim[0]>7) {\r\n    // Even if there were other /byte/ orders, we wouldn't be able to detect them using a short (16 bits, so only two bytes).\r\n    console.warn(\"dim[0] is out-of-range, we'll simply try continuing to read the file, but this will most likely fail horribly.\")\r\n  }\r\n  \r\n  // Now check header size and magic\r\n  var sizeof_hdr = view.getInt32(0, littleEndian)\r\n  if (sizeof_hdr !== 348 && (1>dim[0] || dim[0]>7)) {\r\n    // Try to recover from weird dim info\r\n    littleEndian = !littleEndian\r\n    dim[0] = view.getInt16(40, littleEndian)\r\n    sizeof_hdr = view.getInt32(0, littleEndian)\r\n    if (sizeof_hdr !== 348) {\r\n      throw new Error(\"I'm sorry, but I really cannot determine the byte order of the (NIfTI) file at all.\")\r\n    }\r\n  } else if (sizeof_hdr < 348) {\r\n    throw new Error(\"Header of file is smaller than expected, I cannot deal with this.\")\r\n  } else if (sizeof_hdr !== 348) {\r\n    console.warn(\"Size of NIfTI header different from what I expect, but I'll try to do my best anyway (might cause trouble).\")\r\n  }\r\n  var magic = String.fromCharCode.apply(null, buf8.subarray(344, 348))\r\n  if (magic !== \"ni1\\0\" && magic !== \"n+1\\0\") {\r\n    throw new Error(\"Sorry, but this does not appear to be a NIfTI-1 file. Maybe Analyze 7.5 format? or NIfTI-2?\")\r\n  }\r\n  \r\n  // Read some more structured header fields\r\n  var dim_info = view.getInt8(39)\r\n  dim.length = 1+Math.min(7, dim[0])\r\n  for(var i=1; i<dim.length; i++) {\r\n    dim[i] = view.getInt16(40+2*i, littleEndian)\r\n    if (dim[i]<=0) {\r\n      console.warn(\"dim[0] was probably wrong or corrupt\")\r\n      dim.length = i\r\n    }\r\n  }\r\n  if (dim.length === 1) throw new Error(\"No valid dimensions!\")\r\n  \r\n  var pixdim = new Array(dim.length)\r\n  for(var i=0; i<pixdim.length; i++) {\r\n    pixdim[i] = view.getFloat32(76+4*i, littleEndian)\r\n  }\r\n  \r\n  var srow = new Float32Array(12)\r\n  for(var i=0; i<12; i++) {\r\n    srow[i] = view.getFloat32(280+4*i, littleEndian)\r\n  }\r\n  \r\n  // Read simple header fields and build up object representing the header\r\n  var header = {\r\n    littleEndian: littleEndian,\r\n    \r\n    sizeof_hdr: sizeof_hdr,\r\n    dim_info: dim_info,\r\n    dim: dim,\r\n    intent_p1: view.getFloat32(56, littleEndian),\r\n    intent_p2: view.getFloat32(56, littleEndian),\r\n    intent_p3: view.getFloat32(56, littleEndian),\r\n    intent_code: view.getInt16(68, littleEndian),\r\n  \r\n    datatype: decodeNIfTIDataType(view.getInt16(70, littleEndian)),\r\n    bitpix: view.getInt16(72, littleEndian),\r\n    slice_start: view.getInt16(74, littleEndian),\r\n    pixdim: pixdim,\r\n    vox_offset: view.getFloat32(108, littleEndian),\r\n    \r\n    scl_slope: view.getFloat32(112, littleEndian),\r\n    scl_inter: view.getFloat32(116, littleEndian),\r\n    slice_end: view.getInt16(120, littleEndian),\r\n    slice_code: view.getInt8(122),\r\n    xyzt_units: decodeNIfTIUnits(view.getInt8(123)),\r\n    cal_max: view.getFloat32(124, littleEndian),\r\n    cal_min: view.getFloat32(128, littleEndian),\r\n    slice_duration: view.getFloat32(132, littleEndian),\r\n    toffset: view.getFloat32(136, littleEndian),\r\n  \r\n    descrip: String.fromCharCode.apply(null, buf8.subarray(148, 228)),\r\n    aux_file: String.fromCharCode.apply(null, buf8.subarray(228, 252)),\r\n  \r\n    qform_code: view.getInt16(252, littleEndian),\r\n    sform_code: view.getInt16(254, littleEndian),\r\n  \r\n    quatern_b: view.getFloat32(256, littleEndian),\r\n    quatern_c: view.getFloat32(260, littleEndian),\r\n    quatern_d: view.getFloat32(264, littleEndian),\r\n    qoffset_x: view.getFloat32(268, littleEndian),\r\n    qoffset_y: view.getFloat32(272, littleEndian),\r\n    qoffset_z: view.getFloat32(276, littleEndian),\r\n    \r\n    srow: srow,\r\n  \r\n    intent_name: String.fromCharCode.apply(null, buf8.subarray(328, 344)),\r\n    \r\n    magic: magic,\r\n  \r\n    extension: buffer.byteLength < 348+4 ? [0,0,0,0] : [view.getInt8(348), view.getInt8(349), view.getInt8(350), view.getInt8(351)]\r\n  }\r\n  \r\n  // Check bitpix\r\n  \r\n  // \"Normalize\" datatype (so that rgb/complex become several normal floats rather than compound types, possibly also do something about bits)\r\n  // Note that there is actually both an rgb datatype and an rgb intent... (My guess is that the datatype corresponds to sizes = [3,dim[0],...], while the intent might correspond to sizes = [dim[0],...,dim[5]=3].)\r\n  \r\n  return header  \r\n}\r\n\r\n// Converts a NIfTI header to an NRRD-compatible structure\r\nfunction NIfTIToNRRD(niftiHeader) {\r\n  var ret = {}\r\n  ret.dimension = niftiHeader.dim[0]\r\n  ret.type = niftiHeader.datatype // TODO: Check that we do not feed anything incompatible?\r\n  ret.encoding = 'raw'\r\n  ret.endian = niftiHeader.littleEndian ? 'little' : 'big'\r\n  ret.sizes = niftiHeader.dim.slice(1) // Note that both NRRD and NIfTI use the convention that the fastest axis comes first!\r\n  ret.thicknesses=niftiHeader.pixdim.slice(1);\r\n\r\n  if (niftiHeader.xyzt_units !== undefined) {\r\n    ret.spaceUnits = niftiHeader.xyzt_units\r\n    while(ret.spaceUnits.length < ret.dimension) { // Pad if necessary\r\n      ret.spaceUnits.push(\"\")\r\n    }\r\n    ret.spaceUnits.length = ret.dimension // Shrink if necessary\r\n  }\r\n  \r\n  if (niftiHeader.qform_code === 0) { // \"method 1\"\r\n    ret.spacings = niftiHeader.pixdim.slice(1)\r\n    while(ret.spacings.length < ret.dimension) {\r\n      ret.spacings.push(NaN)\r\n    }\r\n    ret.spaceDimension = Math.min(ret.dimension, 3) // There might be non-3D data sets? (Although the NIfTI format does seem /heavily/ reliant on assuming a 3D space.)\r\n  } else if (niftiHeader.qform_code > 0) { // \"method 2\"\r\n    // TODO: Figure out exactly what to do with the different qform codes.\r\n    ret.space = \"right-anterior-superior\" // Any method for orientation (except for \"method 1\") uses this, apparently.\r\n    var qfac = niftiHeader.pixdim[0] === 0.0 ? 1 : niftiHeader.pixdim[0]\r\n    var a = Math.sqrt(Math.max(0.0,1.0-(niftiHeader.quatern_b*niftiHeader.quatern_b + niftiHeader.quatern_c*niftiHeader.quatern_c + niftiHeader.quatern_d*niftiHeader.quatern_d)))\r\n    var b = niftiHeader.quatern_b\r\n    var c = niftiHeader.quatern_c\r\n    var d = niftiHeader.quatern_d\r\n    ret.spaceDirections = [\r\n      [niftiHeader.pixdim[1]*(a*a+b*b-c*c-d*d),niftiHeader.pixdim[1]*(2*b*c+2*a*d),niftiHeader.pixdim[1]*(2*b*d-2*a*c)],\r\n      [niftiHeader.pixdim[2]*(2*b*c-2*a*d),niftiHeader.pixdim[2]*(a*a+c*c-b*b-d*d),niftiHeader.pixdim[2]*(2*c*d+2*a*b)],\r\n      [qfac*niftiHeader.pixdim[3]*(2*b*d+2*a*c),qfac*niftiHeader.pixdim[3]*(2*c*d-2*a*b),qfac*niftiHeader.pixdim[3]*(a*a+d*d-c*c-b*b)]]\r\n    ret.spaceOrigin = [niftiHeader.qoffset_x,niftiHeader.qoffset_y,niftiHeader.qoffset_z]\r\n  } else {\r\n    console.warn(\"Invalid qform_code: \" + niftiHeader.qform_code + \", orientation is probably messed up.\")\r\n  }\r\n  // TODO: Here we run into trouble, because in NRRD we cannot expose two DIFFERENT (not complementary, different!) transformations. Even more frustrating is that sform transformations are actually more compatible with NRRD than the qform methods.\r\n  if (niftiHeader.sform_code > 0) {\r\n    console.warn(\"sform transformation are currently ignored.\")\r\n  }\r\n  /*if (niftiHeader.sform_code > 0) { // \"method 3\"\r\n    ret.space = \"right-anterior-superior\" // Any method for orientation (except for \"method 1\") uses this, apparently.\r\n    ret.spaceDirections = [\r\n      [niftiHeader.srow[0*4 + 0],niftiHeader.srow[1*4 + 0],niftiHeader.srow[2*4 + 0]],\r\n      [niftiHeader.srow[0*4 + 1],niftiHeader.srow[1*4 + 1],niftiHeader.srow[2*4 + 1]],\r\n      [niftiHeader.srow[0*4 + 2],niftiHeader.srow[1*4 + 2],niftiHeader.srow[2*4 + 2]]]\r\n    ret.spaceOrigin = [niftiHeader.srow[0*4 + 3],niftiHeader.srow[1*4 + 3],niftiHeader.srow[2*4 + 3]]\r\n  }*/\r\n  // TODO: Enforce that spaceDirections and so on have the correct size.\r\n  \r\n  // TODO: We're still missing an awful lot of info!\r\n  \r\n  return ret\r\n}\r\n\r\n// Just parses the header\r\n// This expects an ArrayBuffer or (Node.js) Buffer\r\nmodule.exports.parseHeader = function (buffer_org) {\r\n  var niftiHeader = parseNIfTIHeader(buffer_org)\r\n  var ret = NIfTIToNRRD(niftiHeader)\r\n  return ret\r\n}\r\n\r\n// Parses both header and data\r\n// This expects an ArrayBuffer or (Node.js) Buffer\r\nmodule.exports.parse = function (buffer_org) {\r\n  var niftiHeader = parseNIfTIHeader(buffer_org)\r\n  var ret = NIfTIToNRRD(niftiHeader)\r\n\r\n  if (niftiHeader.extension[0] !== 0) {\r\n    console.warn(\"Looks like there are extensions in use in this NIfTI file, which will all be ignored!\")\r\n  }\r\n  \r\n  // Read data if it is here\r\n  if (niftiHeader.magic === \"n+1\\0\") {\r\n    var buf8 = new Uint8Array(buffer_org)\r\n    var buffer = buf8.buffer // Make sure we have an ArrayBuffer\r\n    if (niftiHeader.vox_offset<352 || niftiHeader.vox_offset>buffer.byteLength) {\r\n      throw new Error(\"Illegal vox_offset!\")\r\n    }\r\n    ret.buffer = buffer.slice(Math.floor(niftiHeader.vox_offset))\r\n    if (niftiHeader.datatype !== 0) {\r\n      // TODO: It MIGHT make sense to equate DT_UNKNOWN (0) to 'block', with bitpix giving the block size in bits\r\n      ret.data = parseNIfTIRawData(ret.buffer, niftiHeader.datatype, niftiHeader.dim, {endianFlag: niftiHeader.littleEndian})\r\n    }\r\n  }\r\n  \r\n  return ret\r\n}\r\n\r\nfunction parseNIfTIRawData(buffer, type, dim, options) {\r\n  var i, arr, view, totalLen = 1, endianFlag = options.endianFlag, endianness = endianFlag ? 'little' : 'big'\r\n  for(i=1; i<dim.length; i++) {\r\n    totalLen *= dim[i]\r\n  }\r\n  if (type == 'block') {\r\n    // Don't do anything special, just return the slice containing all blocks.\r\n    return buffer.slice(0,totalLen*options.blockSize)\r\n  } else if (type == 'int8' || type == 'uint8' || endianness == systemEndianness) {\r\n    switch(type) {\r\n    case \"int8\":\r\n      checkSize(1)\r\n      return new Int8Array(buffer.slice(0,totalLen))\r\n    case \"uint8\":\r\n      checkSize(1)\r\n      return new Uint8Array(buffer.slice(0,totalLen))\r\n    case \"int16\":\r\n      checkSize(2)\r\n      return new Int16Array(buffer.slice(0,totalLen*2))\r\n    case \"uint16\":\r\n      checkSize(2)\r\n      return new Uint16Array(buffer.slice(0,totalLen*2))\r\n    case \"int32\":\r\n      checkSize(4)\r\n      return new Int32Array(buffer.slice(0,totalLen*4))\r\n    case \"uint32\":\r\n      checkSize(4)\r\n      return new Uint32Array(buffer.slice(0,totalLen*4))\r\n    //case \"int64\":\r\n    //  checkSize(8)\r\n    //  return new Int64Array(buffer.slice(0,totalLen*8))\r\n    //case \"uint64\":\r\n    //  checkSize(8)\r\n    //  return new Uint64Array(buffer.slice(0,totalLen*8))\r\n    case \"float\":\r\n      checkSize(4)\r\n      return new Float32Array(buffer.slice(0,totalLen*4))\r\n    case \"double\":\r\n      checkSize(8)\r\n      return new Float64Array(buffer.slice(0,totalLen*8))\r\n    default:\r\n      console.warn(\"Unsupported NIfTI type: \" + type)\r\n      return undefined\r\n    }\r\n  } else {\r\n    view = new DataView(buffer)\r\n    switch(type) {\r\n    case \"int8\": // Note that here we do not need to check the size of the buffer, as the DataView.get methods should throw an exception if we read beyond the buffer.\r\n      arr = new Int8Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getInt8(i)\r\n      }\r\n      return arr\r\n    case \"uint8\":\r\n      arr = new Uint8Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getUint8(i)\r\n      }\r\n      return arr\r\n    case \"int16\":\r\n      arr = new Int16Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getInt16(i*2)\r\n      }\r\n      return arr\r\n    case \"uint16\":\r\n      arr = new Uint16Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getUint16(i*2)\r\n      }\r\n      return arr\r\n    case \"int32\":\r\n      arr = new Int32Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getInt32(i*4)\r\n      }\r\n      return arr\r\n    case \"uint32\":\r\n      arr = new Uint32Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getUint32(i*4)\r\n      }\r\n      return arr\r\n    //case \"int64\":\r\n    //  arr = new Int64Array(totalLen)\r\n    //  for(i=0; i<totalLen; i++) {\r\n    //    arr[i] = view.getInt64(i*8)\r\n    //  }\r\n    // return arr\r\n    //case \"uint64\":\r\n    //  arr = new Uint64Array(totalLen)\r\n    //  for(i=0; i<totalLen; i++) {\r\n    //    arr[i] = view.getUint64(i*8)\r\n    //  }\r\n    //  return arr\r\n    case \"float\":\r\n      arr = new Float32Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getFloat32(i*4)\r\n      }\r\n      return arr\r\n    case \"double\":\r\n      arr = new Float64Array(totalLen)\r\n      for(i=0; i<totalLen; i++) {\r\n        arr[i] = view.getFloat64(i*8)\r\n      }\r\n      return arr\r\n    default:\r\n      console.warn(\"Unsupported NRRD type: \" + type)\r\n      return undefined\r\n    }\r\n  }\r\n  function checkSize(sizeOfType) {\r\n    if (buffer.byteLength<totalLen*sizeOfType) throw new Error(\"NIfTI file does not contain enough data!\")\r\n  }\r\n}\r\n\r\nfunction decodeNIfTIDataType(datatype) {\r\n  switch(datatype) {\r\n  case 1:\r\n    return 'bit'\r\n  case 2:\r\n    return 'uint8'\r\n  case 4:\r\n    return 'int16'\r\n  case 8:\r\n    return 'int32'\r\n  case 16:\r\n    return 'float'\r\n  case 32:\r\n    return 'complex64'\r\n  case 64:\r\n    return 'double'\r\n  case 128:\r\n    return 'rgb24'\r\n  case 256:\r\n    return 'int8'\r\n  case 512:\r\n    return 'uint16'\r\n  case 768:\r\n    return 'uint32'\r\n  case 1024:\r\n    return 'int64'\r\n  case 1280:\r\n    return 'uint64'\r\n  case 1536:\r\n    return 'float128'\r\n  case 1792:\r\n    return 'complex128'\r\n  case 2048:\r\n    return 'complex256'\r\n  case 2304:\r\n    return 'rgba32'\r\n  default:\r\n    console.warn(\"Unrecognized NIfTI data type: \" + datatype)\r\n    return datatype\r\n  }\r\n}\r\n\r\nfunction decodeNIfTIUnits(units) {\r\n  var space, time\r\n  switch(units & 7) {\r\n  case 0:\r\n    space = \"\"\r\n    break\r\n  case 1:\r\n    space = \"m\"\r\n    break\r\n  case 2:\r\n    space = \"mm\"\r\n    break\r\n  case 3:\r\n    space = \"um\"\r\n    break\r\n  default:\r\n    console.warn(\"Unrecognized NIfTI unit: \" + (units&7))\r\n    space = \"\"\r\n  }\r\n  switch(units & 56) {\r\n  case 0:\r\n    time = \"\"\r\n    break\r\n  case 8:\r\n    time = \"s\"\r\n    break\r\n  case 16:\r\n    time = \"ms\"\r\n    break\r\n  case 24:\r\n    time = \"us\"\r\n    break\r\n  case 32:\r\n    time = \"Hz\"\r\n    break\r\n  case 40:\r\n    time = \"ppm\"\r\n    break\r\n  case 48:\r\n    time = \"rad/s\"\r\n    break\r\n  default:\r\n    console.warn(\"Unrecognized NIfTI unit: \" + (units&56))\r\n    time = \"\"\r\n  }\r\n  return (space === \"\" && time === \"\") ? undefined : [space, space, space, time]\r\n}\r\n\r\n"]}