@d5techs/3dgs-lib 1.4.16 → 1.4.18

package/dist/3dgs-lib.cjs

@@ -4658,7 +4658,7 @@ function createSeededRandom(seed) {
     return ((t ^ t >>> 14) >>> 0) / 4294967296;
   };
 }
-const SH_C0$1 = 0.28209479177387814;
+const SH_C0$2 = 0.28209479177387814;
 function computeImportanceSampling(buffer, dataOffset, stride, totalCount, sampleCount, opacityOffset, opacityType, scale0Offset, scale0Type, scale1Offset, scale1Type, scale2Offset, scale2Type, littleEndian, seed) {
   const dataView = new DataView(buffer, dataOffset);
   const random = createSeededRandom(seed);
@@ -4848,9 +4848,9 @@ async function parsePLYBuffer(buffer, options = {}) {
     const f_dc_0 = offsets.f_dc_0 >= 0 ? readProperty(dataView, base + offsets.f_dc_0, types.f_dc_0, littleEndian) : 0;
     const f_dc_1 = offsets.f_dc_1 >= 0 ? readProperty(dataView, base + offsets.f_dc_1, types.f_dc_1, littleEndian) : 0;
     const f_dc_2 = offsets.f_dc_2 >= 0 ? readProperty(dataView, base + offsets.f_dc_2, types.f_dc_2, littleEndian) : 0;
-    colors[outputIdx * 3 + 0] = 0.5 + SH_C0$1 * f_dc_0;
-    colors[outputIdx * 3 + 1] = 0.5 + SH_C0$1 * f_dc_1;
-    colors[outputIdx * 3 + 2] = 0.5 + SH_C0$1 * f_dc_2;
+    colors[outputIdx * 3 + 0] = 0.5 + SH_C0$2 * f_dc_0;
+    colors[outputIdx * 3 + 1] = 0.5 + SH_C0$2 * f_dc_1;
+    colors[outputIdx * 3 + 2] = 0.5 + SH_C0$2 * f_dc_2;
     const rawOpacity = offsets.opacity >= 0 ? readProperty(dataView, base + offsets.opacity, types.opacity, littleEndian) : 0;
     opacities[outputIdx] = sigmoid(rawOpacity);
     if (shCoeffs && shProps.length > 0) {
@@ -5131,7 +5131,7 @@ var freb = function(eb, start) {
 };
 var _a = freb(fleb, 2), fl = _a.b, revfl = _a.r;
 fl[28] = 258, revfl[258] = 28;
-var _b = freb(fdeb, 0), fd = _b.b;
+var _b = freb(fdeb, 0), fd = _b.b, revfd = _b.r;
 var rev = new u16(32768);
 for (var i = 0; i < 32768; ++i) {
   var x = (i & 43690) >> 1 | (i & 21845) << 1;
@@ -5187,8 +5187,8 @@ for (var i = 280; i < 288; ++i)
 var fdt = new u8(32);
 for (var i = 0; i < 32; ++i)
   fdt[i] = 5;
-var flrm = /* @__PURE__ */ hMap(flt, 9, 1);
-var fdrm = /* @__PURE__ */ hMap(fdt, 5, 1);
+var flm = /* @__PURE__ */ hMap(flt, 9, 0), flrm = /* @__PURE__ */ hMap(flt, 9, 1);
+var fdm = /* @__PURE__ */ hMap(fdt, 5, 0), fdrm = /* @__PURE__ */ hMap(fdt, 5, 1);
 var max = function(a) {
   var m = a[0];
   for (var i = 1; i < a.length; ++i) {
@@ -5382,7 +5382,353 @@ var inflt = function(dat, st, buf, dict) {
   } while (!final);
   return bt != buf.length && noBuf ? slc(buf, 0, bt) : buf.subarray(0, bt);
 };
+var wbits = function(d, p, v) {
+  v <<= p & 7;
+  var o = p / 8 | 0;
+  d[o] |= v;
+  d[o + 1] |= v >> 8;
+};
+var wbits16 = function(d, p, v) {
+  v <<= p & 7;
+  var o = p / 8 | 0;
+  d[o] |= v;
+  d[o + 1] |= v >> 8;
+  d[o + 2] |= v >> 16;
+};
+var hTree = function(d, mb) {
+  var t = [];
+  for (var i = 0; i < d.length; ++i) {
+    if (d[i])
+      t.push({ s: i, f: d[i] });
+  }
+  var s = t.length;
+  var t2 = t.slice();
+  if (!s)
+    return { t: et, l: 0 };
+  if (s == 1) {
+    var v = new u8(t[0].s + 1);
+    v[t[0].s] = 1;
+    return { t: v, l: 1 };
+  }
+  t.sort(function(a, b) {
+    return a.f - b.f;
+  });
+  t.push({ s: -1, f: 25001 });
+  var l = t[0], r = t[1], i0 = 0, i1 = 1, i2 = 2;
+  t[0] = { s: -1, f: l.f + r.f, l, r };
+  while (i1 != s - 1) {
+    l = t[t[i0].f < t[i2].f ? i0++ : i2++];
+    r = t[i0 != i1 && t[i0].f < t[i2].f ? i0++ : i2++];
+    t[i1++] = { s: -1, f: l.f + r.f, l, r };
+  }
+  var maxSym = t2[0].s;
+  for (var i = 1; i < s; ++i) {
+    if (t2[i].s > maxSym)
+      maxSym = t2[i].s;
+  }
+  var tr = new u16(maxSym + 1);
+  var mbt = ln(t[i1 - 1], tr, 0);
+  if (mbt > mb) {
+    var i = 0, dt = 0;
+    var lft = mbt - mb, cst = 1 << lft;
+    t2.sort(function(a, b) {
+      return tr[b.s] - tr[a.s] || a.f - b.f;
+    });
+    for (; i < s; ++i) {
+      var i2_1 = t2[i].s;
+      if (tr[i2_1] > mb) {
+        dt += cst - (1 << mbt - tr[i2_1]);
+        tr[i2_1] = mb;
+      } else
+        break;
+    }
+    dt >>= lft;
+    while (dt > 0) {
+      var i2_2 = t2[i].s;
+      if (tr[i2_2] < mb)
+        dt -= 1 << mb - tr[i2_2]++ - 1;
+      else
+        ++i;
+    }
+    for (; i >= 0 && dt; --i) {
+      var i2_3 = t2[i].s;
+      if (tr[i2_3] == mb) {
+        --tr[i2_3];
+        ++dt;
+      }
+    }
+    mbt = mb;
+  }
+  return { t: new u8(tr), l: mbt };
+};
+var ln = function(n, l, d) {
+  return n.s == -1 ? Math.max(ln(n.l, l, d + 1), ln(n.r, l, d + 1)) : l[n.s] = d;
+};
+var lc = function(c) {
+  var s = c.length;
+  while (s && !c[--s])
+    ;
+  var cl = new u16(++s);
+  var cli = 0, cln = c[0], cls = 1;
+  var w = function(v) {
+    cl[cli++] = v;
+  };
+  for (var i = 1; i <= s; ++i) {
+    if (c[i] == cln && i != s)
+      ++cls;
+    else {
+      if (!cln && cls > 2) {
+        for (; cls > 138; cls -= 138)
+          w(32754);
+        if (cls > 2) {
+          w(cls > 10 ? cls - 11 << 5 | 28690 : cls - 3 << 5 | 12305);
+          cls = 0;
+        }
+      } else if (cls > 3) {
+        w(cln), --cls;
+        for (; cls > 6; cls -= 6)
+          w(8304);
+        if (cls > 2)
+          w(cls - 3 << 5 | 8208), cls = 0;
+      }
+      while (cls--)
+        w(cln);
+      cls = 1;
+      cln = c[i];
+    }
+  }
+  return { c: cl.subarray(0, cli), n: s };
+};
+var clen = function(cf, cl) {
+  var l = 0;
+  for (var i = 0; i < cl.length; ++i)
+    l += cf[i] * cl[i];
+  return l;
+};
+var wfblk = function(out, pos, dat) {
+  var s = dat.length;
+  var o = shft(pos + 2);
+  out[o] = s & 255;
+  out[o + 1] = s >> 8;
+  out[o + 2] = out[o] ^ 255;
+  out[o + 3] = out[o + 1] ^ 255;
+  for (var i = 0; i < s; ++i)
+    out[o + i + 4] = dat[i];
+  return (o + 4 + s) * 8;
+};
+var wblk = function(dat, out, final, syms, lf, df, eb, li, bs, bl, p) {
+  wbits(out, p++, final);
+  ++lf[256];
+  var _a2 = hTree(lf, 15), dlt = _a2.t, mlb = _a2.l;
+  var _b2 = hTree(df, 15), ddt = _b2.t, mdb = _b2.l;
+  var _c = lc(dlt), lclt = _c.c, nlc = _c.n;
+  var _d = lc(ddt), lcdt = _d.c, ndc = _d.n;
+  var lcfreq = new u16(19);
+  for (var i = 0; i < lclt.length; ++i)
+    ++lcfreq[lclt[i] & 31];
+  for (var i = 0; i < lcdt.length; ++i)
+    ++lcfreq[lcdt[i] & 31];
+  var _e = hTree(lcfreq, 7), lct = _e.t, mlcb = _e.l;
+  var nlcc = 19;
+  for (; nlcc > 4 && !lct[clim[nlcc - 1]]; --nlcc)
+    ;
+  var flen = bl + 5 << 3;
+  var ftlen = clen(lf, flt) + clen(df, fdt) + eb;
+  var dtlen = clen(lf, dlt) + clen(df, ddt) + eb + 14 + 3 * nlcc + clen(lcfreq, lct) + 2 * lcfreq[16] + 3 * lcfreq[17] + 7 * lcfreq[18];
+  if (bs >= 0 && flen <= ftlen && flen <= dtlen)
+    return wfblk(out, p, dat.subarray(bs, bs + bl));
+  var lm, ll, dm, dl;
+  wbits(out, p, 1 + (dtlen < ftlen)), p += 2;
+  if (dtlen < ftlen) {
+    lm = hMap(dlt, mlb, 0), ll = dlt, dm = hMap(ddt, mdb, 0), dl = ddt;
+    var llm = hMap(lct, mlcb, 0);
+    wbits(out, p, nlc - 257);
+    wbits(out, p + 5, ndc - 1);
+    wbits(out, p + 10, nlcc - 4);
+    p += 14;
+    for (var i = 0; i < nlcc; ++i)
+      wbits(out, p + 3 * i, lct[clim[i]]);
+    p += 3 * nlcc;
+    var lcts = [lclt, lcdt];
+    for (var it = 0; it < 2; ++it) {
+      var clct = lcts[it];
+      for (var i = 0; i < clct.length; ++i) {
+        var len = clct[i] & 31;
+        wbits(out, p, llm[len]), p += lct[len];
+        if (len > 15)
+          wbits(out, p, clct[i] >> 5 & 127), p += clct[i] >> 12;
+      }
+    }
+  } else {
+    lm = flm, ll = flt, dm = fdm, dl = fdt;
+  }
+  for (var i = 0; i < li; ++i) {
+    var sym = syms[i];
+    if (sym > 255) {
+      var len = sym >> 18 & 31;
+      wbits16(out, p, lm[len + 257]), p += ll[len + 257];
+      if (len > 7)
+        wbits(out, p, sym >> 23 & 31), p += fleb[len];
+      var dst = sym & 31;
+      wbits16(out, p, dm[dst]), p += dl[dst];
+      if (dst > 3)
+        wbits16(out, p, sym >> 5 & 8191), p += fdeb[dst];
+    } else {
+      wbits16(out, p, lm[sym]), p += ll[sym];
+    }
+  }
+  wbits16(out, p, lm[256]);
+  return p + ll[256];
+};
+var deo = /* @__PURE__ */ new i32([65540, 131080, 131088, 131104, 262176, 1048704, 1048832, 2114560, 2117632]);
 var et = /* @__PURE__ */ new u8(0);
+var dflt = function(dat, lvl, plvl, pre, post, st) {
+  var s = st.z || dat.length;
+  var o = new u8(pre + s + 5 * (1 + Math.ceil(s / 7e3)) + post);
+  var w = o.subarray(pre, o.length - post);
+  var lst = st.l;
+  var pos = (st.r || 0) & 7;
+  if (lvl) {
+    if (pos)
+      w[0] = st.r >> 3;
+    var opt = deo[lvl - 1];
+    var n = opt >> 13, c = opt & 8191;
+    var msk_1 = (1 << plvl) - 1;
+    var prev = st.p || new u16(32768), head = st.h || new u16(msk_1 + 1);
+    var bs1_1 = Math.ceil(plvl / 3), bs2_1 = 2 * bs1_1;
+    var hsh = function(i2) {
+      return (dat[i2] ^ dat[i2 + 1] << bs1_1 ^ dat[i2 + 2] << bs2_1) & msk_1;
+    };
+    var syms = new i32(25e3);
+    var lf = new u16(288), df = new u16(32);
+    var lc_1 = 0, eb = 0, i = st.i || 0, li = 0, wi = st.w || 0, bs = 0;
+    for (; i + 2 < s; ++i) {
+      var hv = hsh(i);
+      var imod = i & 32767, pimod = head[hv];
+      prev[imod] = pimod;
+      head[hv] = imod;
+      if (wi <= i) {
+        var rem = s - i;
+        if ((lc_1 > 7e3 || li > 24576) && (rem > 423 || !lst)) {
+          pos = wblk(dat, w, 0, syms, lf, df, eb, li, bs, i - bs, pos);
+          li = lc_1 = eb = 0, bs = i;
+          for (var j = 0; j < 286; ++j)
+            lf[j] = 0;
+          for (var j = 0; j < 30; ++j)
+            df[j] = 0;
+        }
+        var l = 2, d = 0, ch_1 = c, dif = imod - pimod & 32767;
+        if (rem > 2 && hv == hsh(i - dif)) {
+          var maxn = Math.min(n, rem) - 1;
+          var maxd = Math.min(32767, i);
+          var ml = Math.min(258, rem);
+          while (dif <= maxd && --ch_1 && imod != pimod) {
+            if (dat[i + l] == dat[i + l - dif]) {
+              var nl = 0;
+              for (; nl < ml && dat[i + nl] == dat[i + nl - dif]; ++nl)
+                ;
+              if (nl > l) {
+                l = nl, d = dif;
+                if (nl > maxn)
+                  break;
+                var mmd = Math.min(dif, nl - 2);
+                var md = 0;
+                for (var j = 0; j < mmd; ++j) {
+                  var ti = i - dif + j & 32767;
+                  var pti = prev[ti];
+                  var cd = ti - pti & 32767;
+                  if (cd > md)
+                    md = cd, pimod = ti;
+                }
+              }
+            }
+            imod = pimod, pimod = prev[imod];
+            dif += imod - pimod & 32767;
+          }
+        }
+        if (d) {
+          syms[li++] = 268435456 | revfl[l] << 18 | revfd[d];
+          var lin = revfl[l] & 31, din = revfd[d] & 31;
+          eb += fleb[lin] + fdeb[din];
+          ++lf[257 + lin];
+          ++df[din];
+          wi = i + l;
+          ++lc_1;
+        } else {
+          syms[li++] = dat[i];
+          ++lf[dat[i]];
+        }
+      }
+    }
+    for (i = Math.max(i, wi); i < s; ++i) {
+      syms[li++] = dat[i];
+      ++lf[dat[i]];
+    }
+    pos = wblk(dat, w, lst, syms, lf, df, eb, li, bs, i - bs, pos);
+    if (!lst) {
+      st.r = pos & 7 | w[pos / 8 | 0] << 3;
+      pos -= 7;
+      st.h = head, st.p = prev, st.i = i, st.w = wi;
+    }
+  } else {
+    for (var i = st.w || 0; i < s + lst; i += 65535) {
+      var e = i + 65535;
+      if (e >= s) {
+        w[pos / 8 | 0] = lst;
+        e = s;
+      }
+      pos = wfblk(w, pos + 1, dat.subarray(i, e));
+    }
+    st.i = s;
+  }
+  return slc(o, 0, pre + shft(pos) + post);
+};
+var crct = /* @__PURE__ */ function() {
+  var t = new Int32Array(256);
+  for (var i = 0; i < 256; ++i) {
+    var c = i, k = 9;
+    while (--k)
+      c = (c & 1 && -306674912) ^ c >>> 1;
+    t[i] = c;
+  }
+  return t;
+}();
+var crc = function() {
+  var c = -1;
+  return {
+    p: function(d) {
+      var cr = c;
+      for (var i = 0; i < d.length; ++i)
+        cr = crct[cr & 255 ^ d[i]] ^ cr >>> 8;
+      c = cr;
+    },
+    d: function() {
+      return ~c;
+    }
+  };
+};
+var dopt = function(dat, opt, pre, post, st) {
+  if (!st) {
+    st = { l: 1 };
+    if (opt.dictionary) {
+      var dict = opt.dictionary.subarray(-32768);
+      var newDat = new u8(dict.length + dat.length);
+      newDat.set(dict);
+      newDat.set(dat, dict.length);
+      dat = newDat;
+      st.w = dict.length;
+    }
+  }
+  return dflt(dat, opt.level == null ? 6 : opt.level, opt.mem == null ? st.l ? Math.ceil(Math.max(8, Math.min(13, Math.log(dat.length))) * 1.5) : 20 : 12 + opt.mem, pre, post, st);
+};
+var mrg = function(a, b) {
+  var o = {};
+  for (var k in a)
+    o[k] = a[k];
+  for (var k in b)
+    o[k] = b[k];
+  return o;
+};
 var b2 = function(d, b) {
   return d[b] | d[b + 1] << 8;
 };
@@ -5392,9 +5738,30 @@ var b4 = function(d, b) {
 var b8 = function(d, b) {
   return b4(d, b) + b4(d, b + 4) * 4294967296;
 };
+var wbytes = function(d, b, v) {
+  for (; v; ++b)
+    d[b] = v, v >>>= 8;
+};
+function deflateSync(data, opts) {
+  return dopt(data, opts || {}, 0, 0);
+}
 function inflateSync(data, opts) {
   return inflt(data, { i: 2 }, opts && opts.out, opts && opts.dictionary);
 }
+var fltn = function(d, p, t, o) {
+  for (var k in d) {
+    var val = d[k], n = p + k, op = o;
+    if (Array.isArray(val))
+      op = mrg(o, val[1]), val = val[0];
+    if (val instanceof u8)
+      t[n] = [val, op];
+    else {
+      t[n += "/"] = [new u8(0), op];
+      fltn(val, n, t, o);
+    }
+  }
+};
+var te = typeof TextEncoder != "undefined" && /* @__PURE__ */ new TextEncoder();
 var td = typeof TextDecoder != "undefined" && /* @__PURE__ */ new TextDecoder();
 var tds = 0;
 try {
@@ -5418,6 +5785,34 @@ var dutf8 = function(d) {
       r += String.fromCharCode((c & 15) << 12 | (d[i++] & 63) << 6 | d[i++] & 63);
   }
 };
+function strToU8(str, latin1) {
+  var i;
+  if (te)
+    return te.encode(str);
+  var l = str.length;
+  var ar = new u8(str.length + (str.length >> 1));
+  var ai = 0;
+  var w = function(v) {
+    ar[ai++] = v;
+  };
+  for (var i = 0; i < l; ++i) {
+    if (ai + 5 > ar.length) {
+      var n = new u8(ai + 8 + (l - i << 1));
+      n.set(ar);
+      ar = n;
+    }
+    var c = str.charCodeAt(i);
+    if (c < 128 || latin1)
+      w(c);
+    else if (c < 2048)
+      w(192 | c >> 6), w(128 | c & 63);
+    else if (c > 55295 && c < 57344)
+      c = 65536 + (c & 1023 << 10) | str.charCodeAt(++i) & 1023, w(240 | c >> 18), w(128 | c >> 12 & 63), w(128 | c >> 6 & 63), w(128 | c & 63);
+    else
+      w(224 | c >> 12), w(128 | c >> 6 & 63), w(128 | c & 63);
+  }
+  return slc(ar, 0, ai);
+}
 function strFromU8(dat, latin1) {
   if (latin1) {
     var r = "";
@@ -5446,6 +5841,107 @@ var z64e = function(d, b) {
     ;
   return [b8(d, b + 12), b8(d, b + 4), b8(d, b + 20)];
 };
+var exfl = function(ex) {
+  var le = 0;
+  if (ex) {
+    for (var k in ex) {
+      var l = ex[k].length;
+      if (l > 65535)
+        err(9);
+      le += l + 4;
+    }
+  }
+  return le;
+};
+var wzh = function(d, b, f, fn, u, c, ce, co) {
+  var fl2 = fn.length, ex = f.extra, col = co && co.length;
+  var exl = exfl(ex);
+  wbytes(d, b, ce != null ? 33639248 : 67324752), b += 4;
+  if (ce != null)
+    d[b++] = 20, d[b++] = f.os;
+  d[b] = 20, b += 2;
+  d[b++] = f.flag << 1 | (c < 0 && 8), d[b++] = u && 8;
+  d[b++] = f.compression & 255, d[b++] = f.compression >> 8;
+  var dt = new Date(f.mtime == null ? Date.now() : f.mtime), y = dt.getFullYear() - 1980;
+  if (y < 0 || y > 119)
+    err(10);
+  wbytes(d, b, y << 25 | dt.getMonth() + 1 << 21 | dt.getDate() << 16 | dt.getHours() << 11 | dt.getMinutes() << 5 | dt.getSeconds() >> 1), b += 4;
+  if (c != -1) {
+    wbytes(d, b, f.crc);
+    wbytes(d, b + 4, c < 0 ? -c - 2 : c);
+    wbytes(d, b + 8, f.size);
+  }
+  wbytes(d, b + 12, fl2);
+  wbytes(d, b + 14, exl), b += 16;
+  if (ce != null) {
+    wbytes(d, b, col);
+    wbytes(d, b + 6, f.attrs);
+    wbytes(d, b + 10, ce), b += 14;
+  }
+  d.set(fn, b);
+  b += fl2;
+  if (exl) {
+    for (var k in ex) {
+      var exf = ex[k], l = exf.length;
+      wbytes(d, b, +k);
+      wbytes(d, b + 2, l);
+      d.set(exf, b + 4), b += 4 + l;
+    }
+  }
+  if (col)
+    d.set(co, b), b += col;
+  return b;
+};
+var wzf = function(o, b, c, d, e) {
+  wbytes(o, b, 101010256);
+  wbytes(o, b + 8, c);
+  wbytes(o, b + 10, c);
+  wbytes(o, b + 12, d);
+  wbytes(o, b + 16, e);
+};
+function zipSync(data, opts) {
+  if (!opts)
+    opts = {};
+  var r = {};
+  var files = [];
+  fltn(data, "", r, opts);
+  var o = 0;
+  var tot = 0;
+  for (var fn in r) {
+    var _a2 = r[fn], file = _a2[0], p = _a2[1];
+    var compression = p.level == 0 ? 0 : 8;
+    var f = strToU8(fn), s = f.length;
+    var com = p.comment, m = com && strToU8(com), ms = m && m.length;
+    var exl = exfl(p.extra);
+    if (s > 65535)
+      err(11);
+    var d = compression ? deflateSync(file, p) : file, l = d.length;
+    var c = crc();
+    c.p(file);
+    files.push(mrg(p, {
+      size: file.length,
+      crc: c.d(),
+      c: d,
+      f,
+      m,
+      u: s != fn.length || m && com.length != ms,
+      o,
+      compression
+    }));
+    o += 30 + s + exl + l;
+    tot += 76 + 2 * (s + exl) + (ms || 0) + l;
+  }
+  var out = new u8(tot + 22), oe = o, cdl = tot - o;
+  for (var i = 0; i < files.length; ++i) {
+    var f = files[i];
+    wzh(out, f.o, f, f.f, f.u, f.c.length);
+    var badd = 30 + f.f.length + exfl(f.extra);
+    out.set(f.c, f.o + badd);
+    wzh(out, o, f, f.f, f.u, f.c.length, f.o, f.m), o += 16 + badd + (f.m ? f.m.length : 0);
+  }
+  wzf(out, o, files.length, cdl, oe);
+  return out;
+}
 function unzipSync(data, opts) {
   var files = {};
   var e = data.length - 22;
@@ -5480,7 +5976,7 @@ function unzipSync(data, opts) {
   }
   return files;
 }
-const SH_C0 = 0.28209479177387814;
+const SH_C0$1 = 0.28209479177387814;
 const COEFFS_PER_BAND = [3, 8, 15];
 function lerp(a, b, t) {
   return a + (b - a) * t;
@@ -5627,9 +6123,9 @@ async function deserializeSOG(data, onProgress) {
     rotations[i * 4 + 1] = qX;
     rotations[i * 4 + 2] = qY;
     rotations[i * 4 + 3] = qZ;
-    colors[i * 3 + 0] = 0.5 + dcCB[sh0Img.data[off + 0]] * SH_C0;
-    colors[i * 3 + 1] = 0.5 + dcCB[sh0Img.data[off + 1]] * SH_C0;
-    colors[i * 3 + 2] = 0.5 + dcCB[sh0Img.data[off + 2]] * SH_C0;
+    colors[i * 3 + 0] = 0.5 + dcCB[sh0Img.data[off + 0]] * SH_C0$1;
+    colors[i * 3 + 1] = 0.5 + dcCB[sh0Img.data[off + 1]] * SH_C0$1;
+    colors[i * 3 + 2] = 0.5 + dcCB[sh0Img.data[off + 2]] * SH_C0$1;
     opacities[i] = sh0Img.data[off + 3] / 255;
     if (shCoeffs && shPalette && labelsImg) {
       const lOff = (py * labelsImg.width + px) * 4;
@@ -5650,6 +6146,258 @@ async function deserializeSOG(data, onProgress) {
   if (onProgress) onProgress(95, "parse");
   return { count, positions, scales, rotations, colors, opacities, shCoeffs };
 }
+const SH_C0 = 0.28209479177387814;
+function symLog(x) {
+  return Math.sign(x) * Math.log(1 + Math.abs(x));
+}
+function buildCodebook(values, size = 256, iterations = 8) {
+  if (values.length === 0) return new Array(size).fill(0);
+  const sorted = Float32Array.from(values).sort();
+  const codebook = new Float64Array(size);
+  for (let i = 0; i < size; i++) {
+    codebook[i] = sorted[Math.floor(i / (size - 1) * (sorted.length - 1))];
+  }
+  const assignments = new Uint8Array(values.length);
+  for (let iter = 0; iter < iterations; iter++) {
+    for (let i = 0; i < values.length; i++) {
+      let bestIdx = 0;
+      let bestDist = Math.abs(values[i] - codebook[0]);
+      for (let j = 1; j < size; j++) {
+        const d = Math.abs(values[i] - codebook[j]);
+        if (d < bestDist) {
+          bestDist = d;
+          bestIdx = j;
+        }
+      }
+      assignments[i] = bestIdx;
+    }
+    const sums = new Float64Array(size);
+    const counts = new Uint32Array(size);
+    for (let i = 0; i < values.length; i++) {
+      sums[assignments[i]] += values[i];
+      counts[assignments[i]]++;
+    }
+    for (let j = 0; j < size; j++) {
+      if (counts[j] > 0) codebook[j] = sums[j] / counts[j];
+    }
+  }
+  return Array.from(codebook);
+}
+function findNearest(codebook, value) {
+  let bestIdx = 0;
+  let bestDist = Math.abs(value - codebook[0]);
+  for (let j = 1; j < codebook.length; j++) {
+    const d = Math.abs(value - codebook[j]);
+    if (d < bestDist) {
+      bestDist = d;
+      bestIdx = j;
+    }
+  }
+  return bestIdx;
+}
+function calcImageDimensions(count) {
+  const width = Math.ceil(Math.sqrt(count));
+  const height = Math.ceil(count / width);
+  return { width, height };
+}
+async function encodeWebP(width, height, pixelData) {
+  let canvas;
+  let ctx;
+  if (typeof OffscreenCanvas !== "undefined") {
+    canvas = new OffscreenCanvas(width, height);
+    ctx = canvas.getContext("2d");
+  } else {
+    canvas = document.createElement("canvas");
+    canvas.width = width;
+    canvas.height = height;
+    ctx = canvas.getContext("2d");
+  }
+  const imageData = ctx.createImageData(width, height);
+  imageData.data.set(pixelData);
+  ctx.putImageData(imageData, 0, 0);
+  let blob;
+  if (canvas instanceof OffscreenCanvas) {
+    blob = await canvas.convertToBlob({ type: "image/webp", quality: 1 });
+  } else {
+    blob = await new Promise((resolve, reject) => {
+      canvas.toBlob(
+        (b) => b ? resolve(b) : reject(new Error("toBlob failed")),
+        "image/webp",
+        1
+      );
+    });
+  }
+  return new Uint8Array(await blob.arrayBuffer());
+}
+async function serializeSOG(data, coordinateSystem = "blender") {
+  const { count, colors, opacities } = data;
+  let { positions, scales, rotations } = data;
+  if (coordinateSystem === "blender") {
+    positions = new Float32Array(positions);
+    scales = new Float32Array(scales);
+    rotations = new Float32Array(rotations);
+    for (let i = 0; i < count; i++) {
+      const i3 = i * 3, i4 = i * 4;
+      const py = positions[i3 + 1], pz = positions[i3 + 2];
+      positions[i3 + 1] = -pz;
+      positions[i3 + 2] = py;
+      const sy = scales[i3 + 1], sz = scales[i3 + 2];
+      scales[i3 + 1] = sz;
+      scales[i3 + 2] = sy;
+      const ry = rotations[i4 + 2], rz = rotations[i4 + 3];
+      rotations[i4 + 2] = -rz;
+      rotations[i4 + 3] = ry;
+    }
+  }
+  const { width, height } = calcImageDimensions(count);
+  const totalPixels = width * height;
+  const logPositions = new Float32Array(count * 3);
+  for (let i = 0; i < count * 3; i++) {
+    logPositions[i] = symLog(positions[i]);
+  }
+  const mins = [Infinity, Infinity, Infinity];
+  const maxs = [-Infinity, -Infinity, -Infinity];
+  for (let i = 0; i < count; i++) {
+    for (let c = 0; c < 3; c++) {
+      const v = logPositions[i * 3 + c];
+      if (v < mins[c]) mins[c] = v;
+      if (v > maxs[c]) maxs[c] = v;
+    }
+  }
+  const meansLData = new Uint8ClampedArray(totalPixels * 4);
+  const meansUData = new Uint8ClampedArray(totalPixels * 4);
+  meansLData.fill(255);
+  meansUData.fill(255);
+  for (let i = 0; i < count; i++) {
+    const off = i * 4;
+    for (let c = 0; c < 3; c++) {
+      const range = maxs[c] - mins[c];
+      const norm = range < 1e-10 ? 0 : (logPositions[i * 3 + c] - mins[c]) / range;
+      const q16 = Math.round(norm * 65535);
+      meansLData[off + c] = q16 & 255;
+      meansUData[off + c] = q16 >> 8 & 255;
+    }
+    meansLData[off + 3] = 255;
+    meansUData[off + 3] = 255;
+  }
+  const logScales = new Float32Array(count * 3);
+  for (let i = 0; i < count * 3; i++) {
+    logScales[i] = Math.log(Math.max(1e-8, scales[i]));
+  }
+  const scaleCodebook = buildCodebook(logScales);
+  const scalesData = new Uint8ClampedArray(totalPixels * 4);
+  scalesData.fill(255);
+  for (let i = 0; i < count; i++) {
+    const off = i * 4;
+    scalesData[off + 0] = findNearest(scaleCodebook, logScales[i * 3 + 0]);
+    scalesData[off + 1] = findNearest(scaleCodebook, logScales[i * 3 + 1]);
+    scalesData[off + 2] = findNearest(scaleCodebook, logScales[i * 3 + 2]);
+    scalesData[off + 3] = 255;
+  }
+  const quatsData = new Uint8ClampedArray(totalPixels * 4);
+  quatsData.fill(255);
+  const SQRT2 = Math.SQRT2;
+  for (let i = 0; i < count; i++) {
+    const off = i * 4;
+    let qw = rotations[i * 4 + 0];
+    let qx = rotations[i * 4 + 1];
+    let qy = rotations[i * 4 + 2];
+    let qz = rotations[i * 4 + 3];
+    const comps = [qw, qx, qy, qz];
+    let largest = 0;
+    let largestAbs = Math.abs(comps[0]);
+    for (let j = 1; j < 4; j++) {
+      const a = Math.abs(comps[j]);
+      if (a > largestAbs) {
+        largestAbs = a;
+        largest = j;
+      }
+    }
+    if (comps[largest] < 0) {
+      qw = -qw;
+      qx = -qx;
+      qy = -qy;
+      qz = -qz;
+    }
+    let r0, r1, r2;
+    if (largest === 0) {
+      r0 = qx;
+      r1 = qy;
+      r2 = qz;
+    } else if (largest === 1) {
+      r0 = qw;
+      r1 = qy;
+      r2 = qz;
+    } else if (largest === 2) {
+      r0 = qw;
+      r1 = qx;
+      r2 = qz;
+    } else {
+      r0 = qw;
+      r1 = qx;
+      r2 = qy;
+    }
+    quatsData[off + 0] = Math.round(Math.max(0, Math.min(255, (r0 / SQRT2 + 0.5) * 255)));
+    quatsData[off + 1] = Math.round(Math.max(0, Math.min(255, (r1 / SQRT2 + 0.5) * 255)));
+    quatsData[off + 2] = Math.round(Math.max(0, Math.min(255, (r2 / SQRT2 + 0.5) * 255)));
+    quatsData[off + 3] = largest + 252;
+  }
+  const dcValues = new Float32Array(count * 3);
+  for (let i = 0; i < count; i++) {
+    dcValues[i * 3 + 0] = (colors[i * 3 + 0] - 0.5) / SH_C0;
+    dcValues[i * 3 + 1] = (colors[i * 3 + 1] - 0.5) / SH_C0;
+    dcValues[i * 3 + 2] = (colors[i * 3 + 2] - 0.5) / SH_C0;
+  }
+  const dcCodebook = buildCodebook(dcValues);
+  const sh0Data = new Uint8ClampedArray(totalPixels * 4);
+  sh0Data.fill(255);
+  for (let i = 0; i < count; i++) {
+    const off = i * 4;
+    sh0Data[off + 0] = findNearest(dcCodebook, dcValues[i * 3 + 0]);
+    sh0Data[off + 1] = findNearest(dcCodebook, dcValues[i * 3 + 1]);
+    sh0Data[off + 2] = findNearest(dcCodebook, dcValues[i * 3 + 2]);
+    sh0Data[off + 3] = Math.round(Math.max(0, Math.min(255, opacities[i] * 255)));
+  }
+  const [meansLWebP, meansUWebP, scalesWebP, quatsWebP, sh0WebP] = await Promise.all([
+    encodeWebP(width, height, meansLData),
+    encodeWebP(width, height, meansUData),
+    encodeWebP(width, height, scalesData),
+    encodeWebP(width, height, quatsData),
+    encodeWebP(width, height, sh0Data)
+  ]);
+  const meta = {
+    version: 2,
+    asset: { generator: "d5techs/3dgs-lib" },
+    count,
+    antialias: false,
+    means: {
+      mins,
+      maxs,
+      files: ["means_l.webp", "means_u.webp"]
+    },
+    scales: {
+      codebook: scaleCodebook,
+      files: ["scales.webp"]
+    },
+    quats: {
+      files: ["quats.webp"]
+    },
+    sh0: {
+      codebook: dcCodebook,
+      files: ["sh0.webp"]
+    }
+  };
+  const metaBytes = new TextEncoder().encode(JSON.stringify(meta));
+  const zipData = zipSync({
+    "meta.json": metaBytes,
+    "means_l.webp": meansLWebP,
+    "means_u.webp": meansUWebP,
+    "scales.webp": scalesWebP,
+    "quats.webp": quatsWebP,
+    "sh0.webp": sh0WebP
+  }, { level: 0 });
+  return zipData.buffer;
+}
 const WORKGROUP_SIZE$1 = 256;
 const RADIX_BITS = 8;
 const RADIX_SIZE = 256;
@@ -16453,6 +17201,50 @@ class SplatEditor {
     const { positions, scales, rotations, colors, opacities, shCoeffs } = this.compactData;
     return exportEditedPLY(positions, scales, rotations, colors, opacities, shCoeffs ?? null, this.splatState);
   }
+  async exportSOG() {
+    if (!this.compactData) return null;
+    const filtered = this.getFilteredCompactData();
+    if (!filtered) return null;
+    return serializeSOG(filtered);
+  }
+  getFilteredCompactData() {
+    if (!this.compactData) return null;
+    const src = this.compactData;
+    const { count, data } = this.splatState;
+    let keepCount = 0;
+    for (let i = 0; i < count; i++) {
+      if (!(data[i] & State.deleted)) keepCount++;
+    }
+    if (keepCount === 0) return null;
+    if (keepCount === count) return src;
+    const positions = new Float32Array(keepCount * 3);
+    const scales = new Float32Array(keepCount * 3);
+    const rotations = new Float32Array(keepCount * 4);
+    const colors = new Float32Array(keepCount * 3);
+    const opacities = new Float32Array(keepCount);
+    let dst = 0;
+    for (let i = 0; i < count; i++) {
+      if (data[i] & State.deleted) continue;
+      const i3 = i * 3, d3 = dst * 3;
+      positions[d3] = src.positions[i3];
+      positions[d3 + 1] = src.positions[i3 + 1];
+      positions[d3 + 2] = src.positions[i3 + 2];
+      scales[d3] = src.scales[i3];
+      scales[d3 + 1] = src.scales[i3 + 1];
+      scales[d3 + 2] = src.scales[i3 + 2];
+      colors[d3] = src.colors[i3];
+      colors[d3 + 1] = src.colors[i3 + 1];
+      colors[d3 + 2] = src.colors[i3 + 2];
+      const i4 = i * 4, d4 = dst * 4;
+      rotations[d4] = src.rotations[i4];
+      rotations[d4 + 1] = src.rotations[i4 + 1];
+      rotations[d4 + 2] = src.rotations[i4 + 2];
+      rotations[d4 + 3] = src.rotations[i4 + 3];
+      opacities[dst] = src.opacities[i];
+      dst++;
+    }
+    return { count: keepCount, positions, scales, rotations, colors, opacities };
+  }
   downloadPLY(filename = "edited.ply") {
     const buffer = this.exportPLY();
     if (!buffer) return;
@@ -17789,5 +18581,6 @@ exports.loadTextureFromBuffer = loadTextureFromBuffer;
 exports.loadTextureFromURL = loadTextureFromURL;
 exports.mergeBoundingBoxes = mergeBoundingBoxes;
 exports.parsePLYBuffer = parsePLYBuffer;
+exports.serializeSOG = serializeSOG;
 exports.transformBoundingBox = transformBoundingBox;
 //# sourceMappingURL=3dgs-lib.cjs.map