EXRLoader.js 23 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968
  1. /**
  2. * @author Richard M. / https://github.com/richardmonette
  3. *
  4. * OpenEXR loader which, currently, supports reading 16 bit half data, in either
  5. * uncompressed or PIZ wavelet compressed form.
  6. *
  7. * Referred to the original Industrial Light & Magic OpenEXR implementation and the TinyEXR / Syoyo Fujita
  8. * implementation, so I have preserved their copyright notices.
  9. */
  10. // /*
  11. // Copyright (c) 2014 - 2017, Syoyo Fujita
  12. // All rights reserved.
  13. // Redistribution and use in source and binary forms, with or without
  14. // modification, are permitted provided that the following conditions are met:
  15. // * Redistributions of source code must retain the above copyright
  16. // notice, this list of conditions and the following disclaimer.
  17. // * Redistributions in binary form must reproduce the above copyright
  18. // notice, this list of conditions and the following disclaimer in the
  19. // documentation and/or other materials provided with the distribution.
  20. // * Neither the name of the Syoyo Fujita nor the
  21. // names of its contributors may be used to endorse or promote products
  22. // derived from this software without specific prior written permission.
  23. // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
  24. // ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
  25. // WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
  26. // DISCLAIMED. IN NO EVENT SHALL <COPYRIGHT HOLDER> BE LIABLE FOR ANY
  27. // DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
  28. // (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
  29. // LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
  30. // ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  31. // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
  32. // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  33. // */
  34. // // TinyEXR contains some OpenEXR code, which is licensed under ------------
  35. // ///////////////////////////////////////////////////////////////////////////
  36. // //
  37. // // Copyright (c) 2002, Industrial Light & Magic, a division of Lucas
  38. // // Digital Ltd. LLC
  39. // //
  40. // // All rights reserved.
  41. // //
  42. // // Redistribution and use in source and binary forms, with or without
  43. // // modification, are permitted provided that the following conditions are
  44. // // met:
  45. // // * Redistributions of source code must retain the above copyright
  46. // // notice, this list of conditions and the following disclaimer.
  47. // // * Redistributions in binary form must reproduce the above
  48. // // copyright notice, this list of conditions and the following disclaimer
  49. // // in the documentation and/or other materials provided with the
  50. // // distribution.
  51. // // * Neither the name of Industrial Light & Magic nor the names of
  52. // // its contributors may be used to endorse or promote products derived
  53. // // from this software without specific prior written permission.
  54. // //
  55. // // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
  56. // // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
  57. // // LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
  58. // // A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
  59. // // OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
  60. // // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
  61. // // LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
  62. // // DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
  63. // // THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  64. // // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
  65. // // OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  66. // //
  67. // ///////////////////////////////////////////////////////////////////////////
  68. // // End of OpenEXR license -------------------------------------------------
  69. THREE.EXRLoader = function ( manager ) {
  70. this.manager = ( manager !== undefined ) ? manager : THREE.DefaultLoadingManager;
  71. };
  72. THREE.EXRLoader.prototype = Object.create( THREE.DataTextureLoader.prototype );
  73. THREE.EXRLoader.prototype._parser = function ( buffer ) {
  74. const USHORT_RANGE = (1 << 16);
  75. const BITMAP_SIZE = (USHORT_RANGE >> 3);
  76. const HUF_ENCBITS = 16; // literal (value) bit length
  77. const HUF_DECBITS = 14; // decoding bit size (>= 8)
  78. const HUF_ENCSIZE = (1 << HUF_ENCBITS) + 1; // encoding table size
  79. const HUF_DECSIZE = 1 << HUF_DECBITS; // decoding table size
  80. const HUF_DECMASK = HUF_DECSIZE - 1;
  81. const SHORT_ZEROCODE_RUN = 59;
  82. const LONG_ZEROCODE_RUN = 63;
  83. const SHORTEST_LONG_RUN = 2 + LONG_ZEROCODE_RUN - SHORT_ZEROCODE_RUN;
  84. const LONGEST_LONG_RUN = 255 + SHORTEST_LONG_RUN;
  85. const BYTES_PER_HALF = 2;
  86. function reverseLutFromBitmap(bitmap, lut) {
  87. var k = 0;
  88. for (var i = 0; i < USHORT_RANGE; ++i) {
  89. if ((i == 0) || (bitmap[i >> 3] & (1 << (i & 7)))) {
  90. lut[k++] = i;
  91. }
  92. }
  93. var n = k - 1;
  94. while (k < USHORT_RANGE) lut[k++] = 0;
  95. return n;
  96. }
  97. function hufClearDecTable(hdec) {
  98. for (var i = 0; i < HUF_DECSIZE; i++) {
  99. hdec[i] = {}
  100. hdec[i].len = 0;
  101. hdec[i].lit = 0;
  102. hdec[i].p = null;
  103. }
  104. }
  105. function getBits(nBits, c, lc, inBuffer, inOffset) {
  106. while (lc < nBits) {
  107. c = (c << 8) | parseUint8(inBuffer, inOffset);
  108. lc += 8;
  109. }
  110. lc -= nBits;
  111. return { l: (c >> lc) & ((1 << nBits) - 1), c: c, lc: lc };
  112. }
  113. function hufCanonicalCodeTable(hcode) {
  114. var n = new Array(59);
  115. for (var i = 0; i <= 58; ++i) n[i] = 0;
  116. for (var i = 0; i < HUF_ENCSIZE; ++i) n[hcode[i]] += 1;
  117. var c = 0;
  118. for (var i = 58; i > 0; --i) {
  119. var nc = ((c + n[i]) >> 1);
  120. n[i] = c;
  121. c = nc;
  122. }
  123. for (var i = 0; i < HUF_ENCSIZE; ++i) {
  124. var l = hcode[i];
  125. if (l > 0) hcode[i] = l | (n[l]++ << 6);
  126. }
  127. }
  128. function hufUnpackEncTable(inBuffer, inOffset, ni, im, iM, hcode) {
  129. var p = inOffset;
  130. var c = 0;
  131. var lc = 0;
  132. for (; im <= iM; im++) {
  133. if (p.value - inOffset.value > ni) {
  134. return false;
  135. }
  136. var bits = getBits(6, c, lc, inBuffer, p);
  137. var l = bits.l;
  138. c = bits.c;
  139. lc = bits.lc;
  140. hcode[im] = l;
  141. if (l == LONG_ZEROCODE_RUN) {
  142. if (p.value - inOffset.value > ni) {
  143. throw 'Something wrong with hufUnpackEncTable';
  144. }
  145. var bits = getBits(8, c, lc, inBuffer, p);
  146. var zerun = bits.l + SHORTEST_LONG_RUN;
  147. c = bits.c;
  148. lc = bits.lc;
  149. if (im + zerun > iM + 1) {
  150. throw 'Something wrong with hufUnpackEncTable';
  151. }
  152. while (zerun--) hcode[im++] = 0;
  153. im--;
  154. } else if (l >= SHORT_ZEROCODE_RUN) {
  155. var zerun = l - SHORT_ZEROCODE_RUN + 2;
  156. if (im + zerun > iM + 1) {
  157. throw 'Something wrong with hufUnpackEncTable';
  158. }
  159. while (zerun--) hcode[im++] = 0;
  160. im--;
  161. }
  162. }
  163. hufCanonicalCodeTable(hcode);
  164. }
  165. function hufLength(code) { return code & 63; }
  166. function hufCode(code) { return code >> 6; }
  167. function hufBuildDecTable(hcode, im, iM, hdecod) {
  168. for (; im <= iM; im++) {
  169. var c = hufCode(hcode[im]);
  170. var l = hufLength(hcode[im]);
  171. if (c >> l) {
  172. throw 'Invalid table entry';
  173. }
  174. if (l > HUF_DECBITS) {
  175. var pl = hdecod[(c >> (l - HUF_DECBITS))];
  176. if (pl.len) {
  177. throw 'Invalid table entry';
  178. }
  179. pl.lit++;
  180. if (pl.p) {
  181. var p = pl.p;
  182. pl.p = new Array(pl.lit);
  183. for (var i = 0; i < pl.lit - 1; ++i) {
  184. pl.p[i] = p[i];
  185. }
  186. } else {
  187. pl.p = new Array(1);
  188. }
  189. pl.p[pl.lit - 1] = im;
  190. } else if (l) {
  191. var plOffset = 0;
  192. for (var i = 1 << (HUF_DECBITS - l); i > 0; i--) {
  193. var pl = hdecod[(c << (HUF_DECBITS - l)) + plOffset];
  194. if (pl.len || pl.p) {
  195. throw 'Invalid table entry';
  196. }
  197. pl.len = l;
  198. pl.lit = im;
  199. plOffset++;
  200. }
  201. }
  202. }
  203. return true;
  204. }
  205. function getChar(c, lc, inBuffer, inOffset) {
  206. c = (c << 8) | parseUint8(inBuffer, inOffset);
  207. lc += 8;
  208. return { c: c, lc: lc };
  209. }
  210. function getCode(po, rlc, c, lc, inBuffer, inOffset, outBuffer, outBufferOffset, outBufferEndOffset) {
  211. if (po == rlc) {
  212. if (lc < 8) {
  213. var temp = getChar(c, lc, inBuffer, inOffset);
  214. c = temp.c;
  215. lc = temp.lc;
  216. }
  217. lc -= 8;
  218. var cs = (c >> lc);
  219. if (out + cs > oe) {
  220. throw 'Issue with getCode';
  221. }
  222. var s = out[-1];
  223. while (cs-- > 0) {
  224. outBuffer[outBufferOffset.value++] = s;
  225. }
  226. } else if (outBufferOffset.value < outBufferEndOffset) {
  227. outBuffer[outBufferOffset.value++] = po;
  228. } else {
  229. throw 'Issue with getCode';
  230. }
  231. return { c: c, lc: lc };
  232. }
  233. var NBITS = 16;
  234. var A_OFFSET = 1 << (NBITS - 1);
  235. var M_OFFSET = 1 << (NBITS - 1);
  236. var MOD_MASK = (1 << NBITS) - 1;
  237. function wdec14(l, h) {
  238. var ls = (new Int16Array([l]))[0];
  239. var hs = (new Int16Array([h]))[0];
  240. var hi = hs;
  241. var ai = ls + (hi & 1) + (hi >> 1);
  242. var as = ai;
  243. var bs = ai - hi;
  244. return {a: as, b: bs}
  245. }
  246. function wav2Decode(j, buffer, nx, ox, ny, oy, mx) {
  247. var n = (nx > ny) ? ny : nx;
  248. var p = 1;
  249. var p2;
  250. while (p <= n) p <<= 1;
  251. p >>= 1;
  252. p2 = p;
  253. p >>= 1;
  254. while (p >= 1) {
  255. var py = 0;
  256. var ey = py + oy * (ny - p2);
  257. var oy1 = oy * p;
  258. var oy2 = oy * p2;
  259. var ox1 = ox * p;
  260. var ox2 = ox * p2;
  261. var i00, i01, i10, i11;
  262. for (; py <= ey; py += oy2) {
  263. var px = py;
  264. var ex = py + ox * (nx - p2);
  265. for (; px <= ex; px += ox2) {
  266. var p01 = px + ox1;
  267. var p10 = px + oy1;
  268. var p11 = p10 + ox1;
  269. var tmp = wdec14(buffer[px + j], buffer[p10 + j]);
  270. i00 = tmp.a;
  271. i10 = tmp.b;
  272. var tmp = wdec14(buffer[p01 + j], buffer[p11 + j]);
  273. i01 = tmp.a;
  274. i11 = tmp.b;
  275. var tmp = wdec14(i00, i01);
  276. buffer[px + j] = tmp.a;
  277. buffer[p01 + j] = tmp.b;
  278. var tmp = wdec14(i10, i11);
  279. buffer[p10 + j] = tmp.a;
  280. buffer[p11 + j] = tmp.b;
  281. }
  282. if (nx & p) {
  283. var p10 = px + oy1;
  284. var tmp = wdec14(buffer[px + j], buffer[p10 + j]);
  285. i00 = tmp.a;
  286. buffer[p10 + j] = tmp.b;
  287. buffer[px + j] = i00;
  288. }
  289. }
  290. if (ny & p) {
  291. var px = py;
  292. var ex = py + ox * (nx - p2);
  293. for (; px <= ex; px += ox2) {
  294. var p01 = px + ox1;
  295. var tmp = wdec14(buffer[px + j], buffer[p01 + j]);
  296. i00 = tmp.a;
  297. buffer[p01 + j] = tmp.b;
  298. buffer[px + j] = i00;
  299. }
  300. }
  301. p2 = p;
  302. p >>= 1;
  303. }
  304. return py;
  305. }
  306. function hufDecode(encodingTable, decodingTable, inBuffer, inOffset, ni, rlc, no, outBuffer, outOffset) {
  307. var c = 0;
  308. var lc = 0;
  309. var outBufferEndOffset = no;
  310. var inOffsetEnd = parseInt(inOffset.value + (ni + 7) / 8);
  311. while (inOffset.value < inOffsetEnd) {
  312. var temp = getChar(c, lc, inBuffer, inOffset);
  313. c = temp.c;
  314. lc = temp.lc;
  315. while (lc >= HUF_DECBITS) {
  316. var index = (c >> (lc - HUF_DECBITS)) & HUF_DECMASK;
  317. var pl = decodingTable[index];
  318. if (pl.len) {
  319. lc -= pl.len;
  320. var temp = getCode(pl.lit, rlc, c, lc, inBuffer, inOffset, outBuffer, outOffset, outBufferEndOffset);
  321. c = temp.c;
  322. lc = temp.lc;
  323. } else {
  324. if (!pl.p) {
  325. throw 'hufDecode issues';
  326. }
  327. var j;
  328. for (j = 0; j < pl.lit; j++) {
  329. var l = hufLength(encodingTable[pl.p[j]]);
  330. while (lc < l && inOffset.value < inOffsetEnd) {
  331. var temp = getChar(c, lc, inBuffer, inOffset);
  332. c = temp.c;
  333. lc = temp.lc;
  334. }
  335. if (lc >= l) {
  336. if (hufCode(encodingTable[pl.p[j]]) ==
  337. ((c >> (lc - l)) & ((1 << l) - 1))) {
  338. lc -= l;
  339. var temp = getCode(pl.p[j], rlc, c, lc, inBuffer, inOffset, outBuffer, outOffset, outBufferEndOffset);
  340. c = temp.c;
  341. lc = temp.lc;
  342. break;
  343. }
  344. }
  345. }
  346. if (j == pl.lit) {
  347. throw 'hufDecode issues';
  348. }
  349. }
  350. }
  351. }
  352. var i = (8 - ni) & 7;
  353. c >>= i;
  354. lc -= i;
  355. while (lc > 0) {
  356. var pl = decodingTable[(c << (HUF_DECBITS - lc)) & HUF_DECMASK];
  357. if (pl.len) {
  358. lc -= pl.len;
  359. var temp = getCode(pl.lit, rlc, c, lc, inBuffer, inOffset, outBuffer, outOffset, outBufferEndOffset);
  360. c = temp.c;
  361. lc = temp.lc;
  362. } else {
  363. throw 'hufDecode issues';
  364. }
  365. }
  366. return true;
  367. }
  368. function hufUncompress(inBuffer, inOffset, nCompressed, outBuffer, outOffset, nRaw) {
  369. var initialInOffset = inOffset.value;
  370. var im = parseUint32(inBuffer, inOffset);
  371. var iM = parseUint32(inBuffer, inOffset);
  372. inOffset.value += 4;
  373. var nBits = parseUint32(inBuffer, inOffset);
  374. inOffset.value += 4;
  375. if (im < 0 || im >= HUF_ENCSIZE || iM < 0 || iM >= HUF_ENCSIZE) {
  376. throw 'Something wrong with HUF_ENCSIZE';
  377. }
  378. var freq = new Array(HUF_ENCSIZE);
  379. var hdec = new Array(HUF_DECSIZE);
  380. hufClearDecTable(hdec);
  381. var ni = nCompressed - (inOffset.value - initialInOffset);
  382. hufUnpackEncTable(inBuffer, inOffset, ni, im, iM, freq);
  383. if (nBits > 8 * (nCompressed - (inOffset.value - initialInOffset))) {
  384. throw 'Something wrong with hufUncompress';
  385. }
  386. hufBuildDecTable(freq, im, iM, hdec);
  387. hufDecode(freq, hdec, inBuffer, inOffset, nBits, iM, nRaw, outBuffer, outOffset);
  388. }
  389. function applyLut(lut, data, nData) {
  390. for (var i = 0; i < nData; ++i) {
  391. data[i] = lut[data[i]];
  392. }
  393. }
  394. function decompressPIZ(outBuffer, outOffset, inBuffer, inOffset, tmpBufSize, num_channels, exrChannelInfos, dataWidth, num_lines) {
  395. var bitmap = new Uint8Array(BITMAP_SIZE);
  396. var minNonZero = parseUint16(inBuffer, inOffset);
  397. var maxNonZero = parseUint16(inBuffer, inOffset);
  398. if (maxNonZero >= BITMAP_SIZE) {
  399. throw 'Something is wrong with PIZ_COMPRESSION BITMAP_SIZE'
  400. }
  401. if (minNonZero <= maxNonZero) {
  402. for (var i = 0; i < maxNonZero - minNonZero + 1; i++) {
  403. bitmap[i + minNonZero] = parseUint8(inBuffer, inOffset);
  404. }
  405. }
  406. var lut = new Uint16Array(USHORT_RANGE);
  407. var maxValue = reverseLutFromBitmap(bitmap, lut);
  408. var length = parseUint32(inBuffer, inOffset);
  409. hufUncompress(inBuffer, inOffset, length, outBuffer, outOffset, tmpBufSize);
  410. var pizChannelData = new Array(num_channels);
  411. var outBufferEnd = 0
  412. for (var i = 0; i < num_channels; i++) {
  413. var exrChannelInfo = exrChannelInfos[i];
  414. var pixelSize = 2; // assumes HALF_FLOAT
  415. pizChannelData[i] = {};
  416. pizChannelData[i]['start'] = outBufferEnd;
  417. pizChannelData[i]['end'] = pizChannelData[i]['start'];
  418. pizChannelData[i]['nx'] = dataWidth;
  419. pizChannelData[i]['ny'] = num_lines;
  420. pizChannelData[i]['size'] = 1;
  421. outBufferEnd += pizChannelData[i].nx * pizChannelData[i].ny * pizChannelData[i].size;
  422. }
  423. var fooOffset = 0;
  424. for (var i = 0; i < num_channels; i++) {
  425. for (var j = 0; j < pizChannelData[i].size; ++j) {
  426. fooOffset += wav2Decode(
  427. j + fooOffset,
  428. outBuffer,
  429. pizChannelData[i].nx,
  430. pizChannelData[i].size,
  431. pizChannelData[i].ny,
  432. pizChannelData[i].nx * pizChannelData[i].size,
  433. maxValue
  434. );
  435. }
  436. }
  437. applyLut(lut, outBuffer, outBufferEnd);
  438. return true;
  439. }
  440. function parseNullTerminatedString( buffer, offset ) {
  441. var uintBuffer = new Uint8Array( buffer );
  442. var endOffset = 0;
  443. while ( uintBuffer[ offset.value + endOffset ] != 0 ) {
  444. endOffset += 1;
  445. }
  446. var stringValue = new TextDecoder().decode(
  447. new Uint8Array( buffer ).slice( offset.value, offset.value + endOffset )
  448. );
  449. offset.value = offset.value + endOffset + 1;
  450. return stringValue;
  451. }
  452. function parseFixedLengthString( buffer, offset, size ) {
  453. var stringValue = new TextDecoder().decode(
  454. new Uint8Array( buffer ).slice( offset.value, offset.value + size )
  455. );
  456. offset.value = offset.value + size;
  457. return stringValue;
  458. }
  459. function parseUlong( buffer, offset ) {
  460. var uLong = new DataView( buffer.slice( offset.value, offset.value + 4 ) ).getUint32( 0, true );
  461. offset.value = offset.value + 8;
  462. return uLong;
  463. }
  464. function parseUint32( buffer, offset ) {
  465. var Uint32 = new DataView( buffer.slice( offset.value, offset.value + 4 ) ).getUint32( 0, true );
  466. offset.value = offset.value + 4;
  467. return Uint32;
  468. }
  469. function parseUint8( buffer, offset ) {
  470. var Uint8 = new DataView( buffer.slice( offset.value, offset.value + 1 ) ).getUint8( 0, true );
  471. offset.value = offset.value + 1;
  472. return Uint8;
  473. }
  474. function parseFloat32( buffer, offset ) {
  475. var float = new DataView( buffer.slice( offset.value, offset.value + 4 ) ).getFloat32( 0, true );
  476. offset.value += 4;
  477. return float;
  478. }
  479. // https://stackoverflow.com/questions/5678432/decompressing-half-precision-floats-in-javascript
  480. function decodeFloat16( binary ) {
  481. var exponent = ( binary & 0x7C00 ) >> 10,
  482. fraction = binary & 0x03FF;
  483. return ( binary >> 15 ? - 1 : 1 ) * (
  484. exponent ?
  485. (
  486. exponent === 0x1F ?
  487. fraction ? NaN : Infinity :
  488. Math.pow( 2, exponent - 15 ) * ( 1 + fraction / 0x400 )
  489. ) :
  490. 6.103515625e-5 * ( fraction / 0x400 )
  491. );
  492. }
  493. function parseUint16( buffer, offset ) {
  494. var Uint16 = new DataView( buffer.slice( offset.value, offset.value + 2 ) ).getUint16( 0, true );
  495. offset.value += 2;
  496. return Uint16;
  497. }
  498. function parseFloat16( buffer, offset ) {
  499. return decodeFloat16( parseUint16( buffer, offset) );
  500. }
  501. function parseChlist( buffer, offset, size ) {
  502. var startOffset = offset.value;
  503. var channels = [];
  504. while ( offset.value < ( startOffset + size - 1 ) ) {
  505. var name = parseNullTerminatedString( buffer, offset );
  506. var pixelType = parseUint32( buffer, offset ); // TODO: Cast this to UINT, HALF or FLOAT
  507. var pLinear = parseUint8( buffer, offset );
  508. offset.value += 3; // reserved, three chars
  509. var xSampling = parseUint32( buffer, offset );
  510. var ySampling = parseUint32( buffer, offset );
  511. channels.push( {
  512. name: name,
  513. pixelType: pixelType,
  514. pLinear: pLinear,
  515. xSampling: xSampling,
  516. ySampling: ySampling
  517. } );
  518. }
  519. offset.value += 1;
  520. return channels;
  521. }
  522. function parseChromaticities( buffer, offset ) {
  523. var redX = parseFloat32( buffer, offset );
  524. var redY = parseFloat32( buffer, offset );
  525. var greenX = parseFloat32( buffer, offset );
  526. var greenY = parseFloat32( buffer, offset );
  527. var blueX = parseFloat32( buffer, offset );
  528. var blueY = parseFloat32( buffer, offset );
  529. var whiteX = parseFloat32( buffer, offset );
  530. var whiteY = parseFloat32( buffer, offset );
  531. return { redX: redX, redY: redY, greenX, greenY, blueX, blueY, whiteX, whiteY };
  532. }
  533. function parseCompression( buffer, offset ) {
  534. var compressionCodes = [
  535. 'NO_COMPRESSION',
  536. 'RLE_COMPRESSION',
  537. 'ZIPS_COMPRESSION',
  538. 'ZIP_COMPRESSION',
  539. 'PIZ_COMPRESSION'
  540. ];
  541. var compression = parseUint8( buffer, offset );
  542. return compressionCodes[ compression ];
  543. }
  544. function parseBox2i( buffer, offset ) {
  545. var xMin = parseUint32( buffer, offset );
  546. var yMin = parseUint32( buffer, offset );
  547. var xMax = parseUint32( buffer, offset );
  548. var yMax = parseUint32( buffer, offset );
  549. return { xMin: xMin, yMin: yMin, xMax: xMax, yMax: yMax };
  550. }
  551. function parseLineOrder( buffer, offset ) {
  552. var lineOrders = [
  553. 'INCREASING_Y'
  554. ];
  555. var lineOrder = parseUint8( buffer, offset );
  556. return lineOrders[ lineOrder ];
  557. }
  558. function parseV2f( buffer, offset ) {
  559. var x = parseFloat32( buffer, offset );
  560. var y = parseFloat32( buffer, offset );
  561. return [ x, y ];
  562. }
  563. function parseValue( buffer, offset, type, size ) {
  564. if ( type == 'string' || type == 'iccProfile' ) {
  565. return parseFixedLengthString( buffer, offset, size );
  566. } else if ( type == 'chlist' ) {
  567. return parseChlist( buffer, offset, size );
  568. } else if ( type == 'chromaticities' ) {
  569. return parseChromaticities( buffer, offset );
  570. } else if ( type == 'compression' ) {
  571. return parseCompression( buffer, offset );
  572. } else if ( type == 'box2i' ) {
  573. return parseBox2i( buffer, offset );
  574. } else if ( type == 'lineOrder' ) {
  575. return parseLineOrder( buffer, offset );
  576. } else if ( type == 'float' ) {
  577. return parseFloat32( buffer, offset );
  578. } else if ( type == 'v2f' ) {
  579. return parseV2f( buffer, offset );
  580. } else {
  581. throw 'Cannot parse value for unsupported type: ' + type;
  582. }
  583. }
  584. var EXRHeader = {};
  585. var magic = new DataView( buffer ).getUint32( 0, true );
  586. var versionByteZero = new DataView( buffer ).getUint8( 4, true );
  587. var fullMask = new DataView( buffer ).getUint8( 5, true );
  588. // start of header
  589. var offset = { value: 8 }; // start at 8, after magic stuff
  590. var keepReading = true;
  591. while ( keepReading ) {
  592. var attributeName = parseNullTerminatedString( buffer, offset );
  593. if ( attributeName == 0 ) {
  594. keepReading = false;
  595. } else {
  596. var attributeType = parseNullTerminatedString( buffer, offset );
  597. var attributeSize = parseUint32( buffer, offset );
  598. var attributeValue = parseValue( buffer, offset, attributeType, attributeSize );
  599. EXRHeader[ attributeName ] = attributeValue;
  600. }
  601. }
  602. // offsets
  603. var dataWindowHeight = EXRHeader.dataWindow.yMax + 1;
  604. var scanlineBlockSize = 1; // 1 for NO_COMPRESSION
  605. if (EXRHeader.compression == 'PIZ_COMPRESSION') {
  606. scanlineBlockSize = 32;
  607. }
  608. var numBlocks = dataWindowHeight / scanlineBlockSize;
  609. for ( var i = 0; i < numBlocks; i ++ ) {
  610. var scanlineOffset = parseUlong( buffer, offset );
  611. }
  612. // we should be passed the scanline offset table, start reading pixel data
  613. var width = EXRHeader.dataWindow.xMax - EXRHeader.dataWindow.xMin + 1;
  614. var height = EXRHeader.dataWindow.yMax - EXRHeader.dataWindow.yMin + 1;
  615. var numChannels = EXRHeader.channels.length;
  616. var byteArray = new Float32Array( width * height * numChannels );
  617. var channelOffsets = {
  618. R: 0,
  619. G: 1,
  620. B: 2,
  621. A: 3
  622. };
  623. if (EXRHeader.compression == 'NO_COMPRESSION') {
  624. for ( var y = 0; y < height; y ++ ) {
  625. var y_scanline = parseUint32( buffer, offset );
  626. var dataSize = parseUint32( buffer, offset );
  627. for ( var channelID = 0; channelID < EXRHeader.channels.length; channelID ++ ) {
  628. if ( EXRHeader.channels[ channelID ].pixelType == 1 ) {
  629. // HALF
  630. for ( var x = 0; x < width; x ++ ) {
  631. var val = parseFloat16( buffer, offset );
  632. var cOff = channelOffsets[ EXRHeader.channels[ channelID ].name ];
  633. byteArray[ ( ( ( width - y_scanline ) * ( height * numChannels ) ) + ( x * numChannels ) ) + cOff ] = val;
  634. }
  635. } else {
  636. throw 'Only supported pixel format is HALF';
  637. }
  638. }
  639. }
  640. } else if (EXRHeader.compression == 'PIZ_COMPRESSION') {
  641. for ( var scanlineBlockIdx = 0; scanlineBlockIdx < height / scanlineBlockSize; scanlineBlockIdx++ ) {
  642. var line_no = parseUint32( buffer, offset );
  643. var data_len = parseUint32( buffer, offset );
  644. var tmpBufferSize = width * scanlineBlockSize * (EXRHeader.channels.length * BYTES_PER_HALF);
  645. var tmpBuffer = new Uint16Array(tmpBufferSize);
  646. var tmpOffset = { value: 0 };
  647. decompressPIZ(tmpBuffer, tmpOffset, buffer, offset, tmpBufferSize, numChannels, EXRHeader.channels, width, scanlineBlockSize);
  648. for ( var line_y = 0; line_y < scanlineBlockSize; line_y ++ ) {
  649. for ( var channelID = 0; channelID < EXRHeader.channels.length; channelID ++ ) {
  650. if ( EXRHeader.channels[ channelID ].pixelType == 1 ) {
  651. // HALF
  652. for ( var x = 0; x < width; x ++ ) {
  653. var cOff = channelOffsets[ EXRHeader.channels[ channelID ].name ];
  654. var val = decodeFloat16(tmpBuffer[ (channelID * (scanlineBlockSize * width)) + (line_y * width) + x ]);
  655. var true_y = line_y + (scanlineBlockIdx * scanlineBlockSize);
  656. byteArray[ ( ( (height - true_y) * ( width * numChannels ) ) + ( x * numChannels ) ) + cOff ] = val;
  657. }
  658. } else {
  659. throw 'Only supported pixel format is HALF';
  660. }
  661. }
  662. }
  663. }
  664. } else {
  665. throw 'Cannot decompress unsupported compression';
  666. }
  667. return {
  668. header: EXRHeader,
  669. width: width,
  670. height: height,
  671. data: byteArray,
  672. format: THREE.RGBFormat,
  673. type: THREE.FloatType
  674. };
  675. };