@ortho-earth/voxel 0.4.0 → 0.6.0

package/dist/voxel-gl.js

@@ -142,7 +142,7 @@ void main() {
 	}
 	${pick ? `discard;` : `if (u_mode != 1 || acc.a < 1.0 / 255.0) discard; o = acc;`}
 }`;
-let glc = null, host = null, prog = null, pickProg = null, boxVao = null, wantDepth = false, dummyTex = null, max3d = 256;
+let budget = Infinity, glc = null, host = null, prog = null, pickProg = null, boxVao = null, wantDepth = false, dummyTex = null, max3d = 256;
 let pickFbo = null, pickTex = null, pickW = 0, pickH = 0;
 const layers = new Map(), probes = [];
 const R_SPHERE = 6371000, R_ELL = 6378137, D2R = Math.PI / 180;
@@ -235,33 +235,15 @@ export function normalMats(A, Ainv) {
     return { nrm, gradM };
 }
 const upCmap = (t, cmap) => { const gl = glc; gl.bindTexture(gl.TEXTURE_2D, t); gl.texImage2D(gl.TEXTURE_2D, 0, gl.RGBA8, 256, 1, 0, gl.RGBA, gl.UNSIGNED_BYTE, cmap); gl.bindTexture(gl.TEXTURE_2D, null); };
-function build(d) {
-    const gl = glc;
-    const [nx, ny, nz] = d.dims;
+/** 元のセルの大きさ（m）＝box は A の列の短い方・region は cellM の短い方 */
+export const cellOf = (A, geo) => geo ? Math.min(...geo.cellM) : Math.min(...[0, 1, 2].map(c => Math.hypot(A[c], A[3 + c], A[6 + c])));
+function mkLvl(b, comps, region, single) {
+    const gl = glc, [nx, ny, nz] = b.dims;
     if (nx > max3d || ny > max3d || nz > max3d)
-        throw new Error(`voxel: dims ${d.dims.join("×")} exceed MAX_3D_TEXTURE_SIZE ${max3d}`);
-    const mk = (data) => {
-        const tex = gl.createTexture();
-        gl.bindTexture(gl.TEXTURE_3D, tex);
-        gl.pixelStorei(gl.UNPACK_ALIGNMENT, 1);
-        if (d.comps === 4)
-            gl.texImage3D(gl.TEXTURE_3D, 0, gl.RGBA8, nx, ny, nz, 0, gl.RGBA, gl.UNSIGNED_BYTE, data);
-        else
-            gl.texImage3D(gl.TEXTURE_3D, 0, gl.R8, nx, ny, nz, 0, gl.RED, gl.UNSIGNED_BYTE, data);
-        for (const [p, v] of [[gl.TEXTURE_MIN_FILTER, gl.LINEAR], [gl.TEXTURE_MAG_FILTER, gl.LINEAR], [gl.TEXTURE_WRAP_S, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_T, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_R, gl.CLAMP_TO_EDGE]])
-            gl.texParameteri(gl.TEXTURE_3D, p, v);
-        gl.bindTexture(gl.TEXTURE_3D, null);
-        return tex;
-    };
-    const texs = (d.frames?.length ? d.frames : [d.tex]).map(mk), frame = Math.min(texs.length - 1, Math.max(0, d.frame ?? 0)), tex = texs[frame];
-    const cmap = gl.createTexture();
-    gl.bindTexture(gl.TEXTURE_2D, cmap);
-    for (const [p, v] of [[gl.TEXTURE_MIN_FILTER, gl.LINEAR], [gl.TEXTURE_MAG_FILTER, gl.LINEAR], [gl.TEXTURE_WRAP_S, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_T, gl.CLAMP_TO_EDGE]])
-        gl.texParameteri(gl.TEXTURE_2D, p, v);
-    upCmap(cmap, d.cmap);
-    const A = Float64Array.from(d.A), Ainv = inv3(A), { nrm, gradM } = normalMats(A, Ainv);
+        throw new Error(`voxel: dims ${b.dims.join("×")} exceed MAX_3D_TEXTURE_SIZE ${max3d}`);
+    const A = Float64Array.from(b.A), Ainv = inv3(A), { nrm, gradM } = normalMats(A, Ainv);
     let vao = null, vb = null, nVerts = 36;
-    const geo = d.shape === "region" ? d.geo ?? null : null;
+    const geo = region ? b.geo ?? null : null;
     if (geo) { // region＝殻（東北上 m）
         vao = gl.createVertexArray();
         gl.bindVertexArray(vao);
@@ -273,13 +255,131 @@ function build(d) {
         gl.bindVertexArray(null);
         nVerts = geo.shell.length / 3;
     }
-    return { shape: geo ? "region" : "box", dims: [nx, ny, nz], comps: d.comps, n: nx * ny * nz, anchor: d.anchor, alt: d.alt || "relative", A, Ainv, radius: d.radius || 0, mode: d.mode, thr: d.thr || 0, step: d.step || 0.5, opacity: d.opacity ?? 1, occl: d.depthTest !== false, clip: d.clip ?? null, properties: d.properties ?? null, range: d.range, tex, cmap, on: d.on !== false, base: 0, dist: 0, nrm, gradM, geo, vao, vb, nVerts, layer: d.layer ?? null, off: d.off ?? null, own: d.own ?? null, full: d.full ?? null, vmax: d.vmax ?? 255, texs, frame };
+    const frames = b.frames?.length ? b.frames : [b.tex];
+    return { s: b.s, dims: [nx, ny, nz], off: single ? null : b.off, own: single ? null : b.own, anchor: b.anchor, A, Ainv, radius: b.radius || 0, nrm, gradM, geo, vao, vb, nVerts, frames, tex: null, up: -1, bytes: nx * ny * nz * comps };
+}
+/** 段 i を層の「いま」にする（置き方・範囲・殻） */
+function useLevel(L, i) {
+    const V = L.levels[i];
+    L.lv = i;
+    L.s = V.s;
+    L.dims = V.dims;
+    L.off = V.off;
+    L.own = V.own;
+    L.anchor = V.anchor;
+    L.A = V.A;
+    L.Ainv = V.Ainv;
+    L.radius = V.radius;
+    L.nrm = V.nrm;
+    L.gradM = V.gradM;
+    L.geo = V.geo;
+    L.vao = V.vao;
+    L.vb = V.vb;
+    L.nVerts = V.nVerts;
+    L.full = V.off ? L.fullFine : null;
+}
+/** いまの段の、見せるこまを GPU に（無ければ作る・こまが違えば上げ直す） */
+function ensure(L) {
+    const gl = glc, V = L.levels[L.lv], f = Math.min(V.frames.length - 1, Math.max(0, L.frame)), [nx, ny, nz] = V.dims, data = V.frames[f];
+    gl.pixelStorei(gl.UNPACK_ALIGNMENT, 1);
+    if (!V.tex) {
+        V.tex = gl.createTexture();
+        gl.bindTexture(gl.TEXTURE_3D, V.tex);
+        if (L.comps === 4)
+            gl.texImage3D(gl.TEXTURE_3D, 0, gl.RGBA8, nx, ny, nz, 0, gl.RGBA, gl.UNSIGNED_BYTE, data);
+        else
+            gl.texImage3D(gl.TEXTURE_3D, 0, gl.R8, nx, ny, nz, 0, gl.RED, gl.UNSIGNED_BYTE, data);
+        for (const [p, v] of [[gl.TEXTURE_MIN_FILTER, gl.LINEAR], [gl.TEXTURE_MAG_FILTER, gl.LINEAR], [gl.TEXTURE_WRAP_S, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_T, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_R, gl.CLAMP_TO_EDGE]])
+            gl.texParameteri(gl.TEXTURE_3D, p, v);
+        gl.bindTexture(gl.TEXTURE_3D, null);
+        V.up = f;
+    }
+    else if (V.up !== f) {
+        gl.bindTexture(gl.TEXTURE_3D, V.tex);
+        gl.texSubImage3D(gl.TEXTURE_3D, 0, 0, 0, 0, nx, ny, nz, L.comps === 4 ? gl.RGBA : gl.RED, gl.UNSIGNED_BYTE, data);
+        gl.bindTexture(gl.TEXTURE_3D, null);
+        V.up = f;
+    }
+    L.tex = V.tex;
+}
+const freeTex = (V) => { if (V.tex) {
+    glc.deleteTexture(V.tex);
+    V.tex = null;
+    V.up = -1;
+} };
+function build(d) {
+    const gl = glc, region = d.shape === "region", single = !d.off;
+    const levels = [mkLvl({ s: 1, off: d.off ?? [0, 0, 0], own: d.own ?? [0, 0, 0, ...d.dims], dims: d.dims, tex: d.tex, frames: d.frames ?? null, vmax: d.vmax ?? 255, anchor: d.anchor, A: d.A, radius: d.radius, geo: d.geo ?? null }, d.comps, region, single), ...(d.lods ?? []).map(b => mkLvl(b, d.comps, region, false))];
+    const cmap = gl.createTexture();
+    gl.bindTexture(gl.TEXTURE_2D, cmap);
+    for (const [p, v] of [[gl.TEXTURE_MIN_FILTER, gl.LINEAR], [gl.TEXTURE_MAG_FILTER, gl.LINEAR], [gl.TEXTURE_WRAP_S, gl.CLAMP_TO_EDGE], [gl.TEXTURE_WRAP_T, gl.CLAMP_TO_EDGE]])
+        gl.texParameteri(gl.TEXTURE_2D, p, v);
+    upCmap(cmap, d.cmap);
+    const V = levels[0], fullFine = d.full ?? d.dims;
+    const L = { shape: region ? "region" : "box", dims: V.dims, comps: d.comps, n: V.dims[0] * V.dims[1] * V.dims[2], anchor: V.anchor, alt: d.alt || "relative", A: V.A, Ainv: V.Ainv, radius: V.radius, mode: d.mode, thr: d.thr || 0, step: d.step || 0.5, opacity: d.opacity ?? 1, occl: d.depthTest !== false, clip: d.clip ?? null, properties: d.properties ?? null, range: d.range,
+        tex: null, cmap, on: d.on !== false, base: 0, dist: 0, nrm: V.nrm, gradM: V.gradM, geo: V.geo, vao: V.vao, vb: V.vb, nVerts: V.nVerts, layer: d.layer ?? null, off: V.off, own: V.own, full: d.full ?? null, vmax: d.vmax ?? 255,
+        frame: Math.max(0, d.frame ?? 0), lv: 0, s: 1, levels, cell: cellOf(V.A, V.geo), fullFine: [fullFine[0], fullFine[1], fullFine[2]] };
+    useLevel(L, 0);
+    return L;
 }
 const drop = (L) => { if (!L)
-    return; const gl = glc; for (const t of L.texs)
-    gl.deleteTexture(t); gl.deleteTexture(L.cmap); if (L.vb)
-    gl.deleteBuffer(L.vb); if (L.vao)
-    gl.deleteVertexArray(L.vao); };
+    return; const gl = glc; for (const V of L.levels) {
+    freeTex(V);
+    if (V.vb)
+        gl.deleteBuffer(V.vb);
+    if (V.vao)
+        gl.deleteVertexArray(V.vao);
+} gl.deleteTexture(L.cmap); };
+/** LoD の割り振り（第 5 段 ③）＝予算（bytes）に収まる中で一番細かく。遠い（元のセルが小さく見える）層から粗く：
+ *  lv_i = clamp(floor(x + log2(dist_i / cell_i)), 0, 段の数 − 1) の x を二分で（x が大きいほど粗い・合計は x に単調）。一番粗くしても越えるなら全部一番粗く */
+export function planLevels(items, budget) {
+    if (!items.length)
+        return [];
+    const g = items.map(it => Math.log2(Math.max(1e-9, it.dist) / Math.max(1e-9, it.cell)));
+    const at = (x) => items.map((it, i) => Math.min(it.bytes.length - 1, Math.max(0, Math.floor(x + g[i]))));
+    const total = (lv) => lv.reduce((a, l, i) => a + items[i].bytes[l], 0);
+    let lo = -Math.max(...g) - 1, hi = -Math.min(...g) + 9; // lo＝全部 0・hi＝全部一番粗い
+    const fine = at(lo);
+    if (total(fine) <= budget)
+        return fine;
+    const coarse = at(hi);
+    if (total(coarse) > budget)
+        return coarse;
+    for (let k = 0; k < 40; k++) {
+        const m = (lo + hi) / 2;
+        if (total(at(m)) <= budget)
+            hi = m;
+        else
+            lo = m;
+    }
+    return at(hi);
+}
+/** 描く層の段を決めてテクスチャを揃え、予算を越えれば描かない層のテクスチャを捨てる */
+function settle(vis) {
+    const lv = planLevels(vis.map(L => ({ dist: L.dist, cell: L.cell, bytes: L.levels.map(V => V.bytes) })), budget);
+    vis.forEach((L, i) => { if (L.lv !== lv[i])
+        useLevel(L, lv[i]); L.levels.forEach((V, k) => { if (k !== L.lv)
+        freeTex(V); }); ensure(L); });
+    let used = 0;
+    for (const L of layers.values())
+        for (const V of L.levels)
+            if (V.tex)
+                used += V.bytes;
+    if (used <= budget)
+        return;
+    const on = new Set(vis);
+    for (const L of layers.values()) {
+        if (on.has(L))
+            continue;
+        for (const V of L.levels)
+            if (V.tex) {
+                used -= V.bytes;
+                freeTex(V);
+            }
+        if (used <= budget)
+            return;
+    }
+}
 function rebase() { let b = 1; for (const L of layers.values()) {
     L.base = b;
     b += L.n;
@@ -290,6 +390,10 @@ export function message(d) {
         host?.requestDraw?.();
         return;
     }
+    if (d.type === "budget" && !glc) {
+        budget = d.bytes > 0 ? d.bytes : Infinity;
+        return;
+    }
     if (!glc)
         return;
     if (d.type === "layer") {
@@ -308,10 +412,8 @@ export function message(d) {
         const L = layers.get(d.q);
         if (!L)
             return;
-        if (d.frame != null && L.texs.length > 1) {
-            L.frame = Math.min(L.texs.length - 1, Math.max(0, d.frame));
-            L.tex = L.texs[L.frame];
-        }
+        if (d.frame != null)
+            L.frame = Math.max(0, d.frame); // 次に描く時に上げ直す（ensure）
         if (d.mode != null)
             L.mode = d.mode;
         if (d.thr != null)
@@ -348,6 +450,10 @@ export function message(d) {
         rebase();
         host?.requestDraw?.();
     }
+    else if (d.type === "budget") {
+        budget = d.bytes > 0 ? d.bytes : Infinity;
+        host?.requestDraw?.();
+    }
 }
 // 層の mvp（box＝セル空間 → clip＝mvp × [E N U o] × A・region＝東北上 m → clip）と、目の位置（同じ空間）。倍精度で掛けてから f32
 const mvpL = new Float64Array(16), mvpC = new Float32Array(16), eyeC = new Float32Array(3), geoU = { o: new Float32Array(3), E: new Float32Array(3), N: new Float32Array(3), U: new Float32Array(3), R: 0, rAx: 1 };
@@ -436,11 +542,11 @@ export function stepsOf(L) {
 }
 /** 断面 → セル空間の範囲 */
 export const clipBounds = (dims, clip) => ({ lo: [0, 1, 2].map(i => (clip ? clip[i * 2] : 0) * dims[i]), hi: [0, 1, 2].map(i => (clip ? clip[i * 2 + 1] : 1) * dims[i]) });
-/** 描く範囲（テクスチャのセル座標）＝断面（親の割合）と持ち分の交わり。空なら null（そのブリックは描かない） */
+/** 描く範囲（テクスチャのセル座標）＝断面（親の割合）と持ち分の交わり。空なら null（そのブリックは描かない）。粗い段（s）は親もその段のセルで測る */
 export function drawBounds(L) {
     if (!L.full || !L.off || !L.own)
         return clipBounds(L.dims, L.clip);
-    const c = clipBounds(L.full, L.clip), lo = [0, 1, 2].map(i => Math.max(L.own[i], c.lo[i] - L.off[i])), hi = [0, 1, 2].map(i => Math.min(L.own[3 + i], c.hi[i] - L.off[i]));
+    const s = L.s ?? 1, c = clipBounds([L.full[0] / s, L.full[1] / s, L.full[2] / s], L.clip), lo = [0, 1, 2].map(i => Math.max(L.own[i], c.lo[i] - L.off[i])), hi = [0, 1, 2].map(i => Math.min(L.own[3 + i], c.hi[i] - L.off[i]));
     return lo.every((v, i) => v < hi[i]) ? { lo, hi } : null;
 }
 /** 閾より下しか無い層（ブリック）＝描かない（vmax＝R8 の 0〜255。RGBA は α 0 の時だけ） */
@@ -456,12 +562,13 @@ export function centerDist(L, enu) {
     }
     return Math.hypot(enu[0] - c[0], enu[1] - c[1], enu[2] - c[2]);
 }
-/** 当たり（層の中の番号）→ 親の層の id・番号・セル（ブリックは off と full で戻す） */
+/** 当たり（層の中の番号）→ 親の層の id・番号・セル（ブリックは off と full で戻す・粗い段は元のセルの真ん中へ） */
 export function pickedCell(q, L, index) {
     const [nx, ny] = L.dims, i = index % nx, j = Math.floor(index / nx) % ny, k = Math.floor(index / (nx * ny));
     if (!L.off || !L.full)
         return { layer: L.layer ?? q, index, cell: [i, j, k] };
-    const cell = [i + L.off[0], j + L.off[1], k + L.off[2]];
+    const s = L.s ?? 1, f = (v, a) => s === 1 ? v : Math.min(L.full[a] - 1, Math.floor((v + 0.5) * s));
+    const cell = [f(i + L.off[0], 0), f(j + L.off[1], 1), f(k + L.off[2], 2)];
     return { layer: L.layer ?? q, index: cell[0] + L.full[0] * (cell[1] + L.full[1] * cell[2]), cell };
 }
 function drawAll(pick, s, api, w, h) {
@@ -486,6 +593,7 @@ function drawAll(pick, s, api, w, h) {
     for (const L of layers.values())
         if (L.on && L.n && !isEmpty(L) && drawBounds(L) && place(L, s, api))
             vis.push(L);
+    settle(vis);
     vis.sort((a, b) => b.dist - a.dist); // 層は奥から
     for (const L of vis) {
         place(L, s, api);
@@ -533,7 +641,7 @@ function answerProbes(w, h) {
         return;
     const gl = glc;
     const ids = probes.splice(0);
-    const ls = [...layers].map(([q, L]) => ({ q, on: L.on, n: L.n, dims: L.dims, mode: L.mode, shape: L.shape, dist: +L.dist.toFixed(4), anchor: L.anchor, alt: L.alt }));
+    const ls = [...layers].map(([q, L]) => ({ q, on: L.on, n: L.n, dims: L.dims, mode: L.mode, shape: L.shape, dist: +L.dist.toFixed(4), anchor: L.anchor, alt: L.alt, lv: L.lv, s: L.s, gpu: L.levels.reduce((a, V) => a + (V.tex ? V.bytes : 0), 0) }));
     try {
         const data = new Uint8Array(w * h * 4);
         gl.readPixels(0, 0, w, h, gl.RGBA, gl.UNSIGNED_BYTE, data);

package/dist/voxel-gpu.js

@@ -209,7 +209,7 @@ const stepsOf = (L) => {
 };
 const align256 = (v) => Math.ceil(v / 256) * 256;
 function drawBounds(L) {
-    const dims = L.full && L.off && L.own ? L.full : L.dims, c = { lo: [0, 1, 2].map(i => (L.clip ? L.clip[i * 2] : 0) * dims[i]), hi: [0, 1, 2].map(i => (L.clip ? L.clip[i * 2 + 1] : 1) * dims[i]) };
+    const s = L.s ?? 1, dims = L.full && L.off && L.own ? [L.full[0] / s, L.full[1] / s, L.full[2] / s] : L.dims, c = { lo: [0, 1, 2].map(i => (L.clip ? L.clip[i * 2] : 0) * dims[i]), hi: [0, 1, 2].map(i => (L.clip ? L.clip[i * 2 + 1] : 1) * dims[i]) };
     if (!L.full || !L.off || !L.own)
         return c;
     const lo = [0, 1, 2].map(i => Math.max(L.own[i], c.lo[i] - L.off[i])), hi = [0, 1, 2].map(i => Math.min(L.own[3 + i], c.hi[i] - L.off[i]));
@@ -230,12 +230,36 @@ function pickedCell(q, L, index) {
     const [nx, ny] = L.dims, i = index % nx, j = Math.floor(index / nx) % ny, k = Math.floor(index / (nx * ny));
     if (!L.off || !L.full)
         return { layer: L.layer ?? q, index, cell: [i, j, k] };
-    const cell = [i + L.off[0], j + L.off[1], k + L.off[2]];
+    const s = L.s ?? 1, f = (v, a) => s === 1 ? v : Math.min(L.full[a] - 1, Math.floor((v + 0.5) * s));
+    const cell = [f(i + L.off[0], 0), f(j + L.off[1], 1), f(k + L.off[2], 2)];
     return { layer: L.layer ?? q, index: cell[0] + L.full[0] * (cell[1] + L.full[1] * cell[2]), cell };
 }
+const cellOf = (A, geo) => geo ? Math.min(...geo.cellM) : Math.min(...[0, 1, 2].map(c => Math.hypot(A[c], A[3 + c], A[6 + c])));
+function planLevels(items, budget) {
+    if (!items.length)
+        return [];
+    const g = items.map(it => Math.log2(Math.max(1e-9, it.dist) / Math.max(1e-9, it.cell)));
+    const at = (x) => items.map((it, i) => Math.min(it.bytes.length - 1, Math.max(0, Math.floor(x + g[i]))));
+    const total = (lv) => lv.reduce((a, l, i) => a + items[i].bytes[l], 0);
+    let lo = -Math.max(...g) - 1, hi = -Math.min(...g) + 9;
+    const fine = at(lo);
+    if (total(fine) <= budget)
+        return fine;
+    const coarse = at(hi);
+    if (total(coarse) > budget)
+        return coarse;
+    for (let k = 0; k < 40; k++) {
+        const m = (lo + hi) / 2;
+        if (total(at(m)) <= budget)
+            hi = m;
+        else
+            lo = m;
+    }
+    return at(hi);
+}
 export function create() {
     let device = null, ctx = null, host = null, format = "rgba8unorm";
-    let shared = false, wantDepth = false, dead = false;
+    let shared = false, wantDepth = false, dead = false, budget = Infinity;
     let drawPipe = null, pickPipe = null, layout = null, sampler = null, dummyDep = null, boxVb = null;
     let pickTex = null, pickW = 0, pickH = 0;
     const layers = new Map(), probes = [];
@@ -306,19 +330,11 @@ export function create() {
         ctx = c;
     }
     const upCmap = (t, cmap) => { device.queue.writeTexture({ texture: t }, cmap, { bytesPerRow: 1024 }, [256, 1]); };
-    function build(d, q) {
-        const [nx, ny, nz] = d.dims, max3d = device.limits.maxTextureDimension3D || 256;
+    function mkLvl(b, q, region, single, comps) {
+        const [nx, ny, nz] = b.dims, max3d = device.limits.maxTextureDimension3D || 256;
         if (nx > max3d || ny > max3d || nz > max3d)
-            throw new Error(`voxel: dims ${d.dims.join("×")} exceed maxTextureDimension3D ${max3d}`);
-        const mk = (data, f) => {
-            const tex = device.createTexture({ label: `voxel-${q}${f ? `@${f}` : ""}`, size: [nx, ny, nz], dimension: "3d", format: d.comps === 4 ? "rgba8unorm" : "r8unorm", usage: TEX.TEXTURE_BINDING | TEX.COPY_DST });
-            device.queue.writeTexture({ texture: tex }, data, { bytesPerRow: nx * d.comps, rowsPerImage: ny }, [nx, ny, nz]);
-            return tex;
-        };
-        const texs = (d.frames?.length ? d.frames : [d.tex]).map(mk), frame = Math.min(texs.length - 1, Math.max(0, d.frame ?? 0)), tex = texs[frame];
-        const cmap = device.createTexture({ label: `voxel-cmap-${q}`, size: [256, 1], format: "rgba8unorm", usage: TEX.TEXTURE_BINDING | TEX.COPY_DST });
-        upCmap(cmap, d.cmap);
-        const A = Float64Array.from(d.A), Ainv = inv3(A), nrm = new Float32Array(9), gradM = new Float32Array(9);
+            throw new Error(`voxel: dims ${b.dims.join("×")} exceed maxTextureDimension3D ${max3d}`);
+        const A = Float64Array.from(b.A), Ainv = inv3(A), nrm = new Float32Array(9), gradM = new Float32Array(9);
         for (let c = 0; c < 3; c++) {
             const l = Math.hypot(A[c], A[3 + c], A[6 + c]) || 1;
             nrm[c * 3] = A[c] / l;
@@ -328,15 +344,96 @@ export function create() {
         for (let c = 0; c < 3; c++)
             for (let r = 0; r < 3; r++)
                 gradM[c * 3 + r] = Ainv[c * 3 + r];
-        const geo = d.shape === "region" ? d.geo ?? null : null;
-        const vb = geo ? mkBuf(geo.shell, BUF.VERTEX, `voxel-shell-${q}`) : boxVb, nVerts = geo ? geo.shell.length / 3 : 36;
+        const geo = region ? b.geo ?? null : null;
+        const vb = geo ? mkBuf(geo.shell, BUF.VERTEX, `voxel-shell-${q}@${b.s}`) : boxVb, nVerts = geo ? geo.shell.length / 3 : 36;
+        return { s: b.s, dims: [nx, ny, nz], off: single ? null : b.off, own: single ? null : b.own, anchor: b.anchor, A, Ainv, radius: b.radius || 0, nrm, gradM, geo, vb, nVerts, frames: b.frames?.length ? b.frames : [b.tex], tex: null, up: -1, bytes: nx * ny * nz * comps };
+    }
+    function useLevel(L, i) {
+        const V = L.levels[i];
+        L.lv = i;
+        L.s = V.s;
+        L.dims = V.dims;
+        L.off = V.off;
+        L.own = V.own;
+        L.anchor = V.anchor;
+        L.A = V.A;
+        L.Ainv = V.Ainv;
+        L.radius = V.radius;
+        L.nrm = V.nrm;
+        L.gradM = V.gradM;
+        L.geo = V.geo;
+        L.vb = V.vb;
+        L.nVerts = V.nVerts;
+        L.full = V.off ? L.fullFine : null;
+    }
+    /** いまの段の見せるこまを GPU に（無ければ作る・こまが違えば上げ直す）。tex が変われば bind group を作り直す */
+    function ensure(L) {
+        const V = L.levels[L.lv], f = Math.min(V.frames.length - 1, Math.max(0, L.frame)), [nx, ny, nz] = V.dims;
+        if (!V.tex)
+            V.tex = device.createTexture({ label: `voxel-${L.q}@${V.s}`, size: [nx, ny, nz], dimension: "3d", format: L.comps === 4 ? "rgba8unorm" : "r8unorm", usage: TEX.TEXTURE_BINDING | TEX.COPY_DST });
+        if (V.up !== f) {
+            device.queue.writeTexture({ texture: V.tex }, V.frames[f], { bytesPerRow: nx * L.comps, rowsPerImage: ny }, [nx, ny, nz]);
+            V.up = f;
+        }
+        if (L.tex !== V.tex) {
+            L.tex = V.tex;
+            L.bg = null;
+        }
+    }
+    const freeTex = (V) => { if (V.tex) {
+        V.tex.destroy();
+        V.tex = null;
+        V.up = -1;
+    } };
+    function build(d, q) {
+        const region = d.shape === "region", single = !d.off;
+        const levels = [mkLvl({ s: 1, off: d.off ?? [0, 0, 0], own: d.own ?? [0, 0, 0, ...d.dims], dims: d.dims, tex: d.tex, frames: d.frames ?? null, vmax: d.vmax ?? 255, anchor: d.anchor, A: d.A, radius: d.radius, geo: d.geo ?? null }, q, region, single, d.comps), ...(d.lods ?? []).map(b => mkLvl(b, q, region, false, d.comps))];
+        const cmap = device.createTexture({ label: `voxel-cmap-${q}`, size: [256, 1], format: "rgba8unorm", usage: TEX.TEXTURE_BINDING | TEX.COPY_DST });
+        upCmap(cmap, d.cmap);
         const ub = device.createBuffer({ label: `voxel-u-${q}`, size: UB_SIZE, usage: BUF.UNIFORM | BUF.COPY_DST });
-        return { q, shape: geo ? "region" : "box", dims: [nx, ny, nz], comps: d.comps, n: nx * ny * nz, anchor: d.anchor, alt: d.alt || "relative", A, Ainv, radius: d.radius || 0, mode: d.mode, thr: d.thr || 0, step: d.step || 0.5, opacity: d.opacity ?? 1, occl: d.depthTest !== false, clip: d.clip ?? null, properties: d.properties ?? null, range: d.range, tex, cmap, vb, nVerts, ub, bg: null, bgDep: null, on: d.on !== false, base: 0, dist: 0, nrm, gradM, geo, layer: d.layer ?? null, off: d.off ?? null, own: d.own ?? null, full: d.full ?? null, vmax: d.vmax ?? 255, texs, frame };
+        const V = levels[0], fullFine = d.full ?? d.dims;
+        const L = { q, shape: region ? "region" : "box", dims: V.dims, comps: d.comps, n: V.dims[0] * V.dims[1] * V.dims[2], anchor: V.anchor, alt: d.alt || "relative", A: V.A, Ainv: V.Ainv, radius: V.radius, mode: d.mode, thr: d.thr || 0, step: d.step || 0.5, opacity: d.opacity ?? 1, occl: d.depthTest !== false, clip: d.clip ?? null, properties: d.properties ?? null, range: d.range,
+            tex: null, cmap, vb: V.vb, nVerts: V.nVerts, ub, bg: null, bgDep: null, on: d.on !== false, base: 0, dist: 0, nrm: V.nrm, gradM: V.gradM, geo: V.geo, layer: d.layer ?? null, off: V.off, own: V.own, full: d.full ?? null, vmax: d.vmax ?? 255,
+            frame: Math.max(0, d.frame ?? 0), lv: 0, s: 1, levels, cell: cellOf(V.A, V.geo), fullFine: [fullFine[0], fullFine[1], fullFine[2]] };
+        useLevel(L, 0);
+        return L;
     }
     const drop = (L) => { if (!L)
-        return; for (const t of L.texs)
-        t.destroy(); L.cmap.destroy(); L.ub.destroy(); if (L.vb !== boxVb)
-        L.vb.destroy(); };
+        return; for (const V of L.levels) {
+        freeTex(V);
+        if (V.vb !== boxVb)
+            V.vb.destroy();
+    } L.cmap.destroy(); L.ub.destroy(); };
+    /** 描く層の段を決めてテクスチャを揃え、予算を越えれば描かない層のテクスチャを捨てる（voxel-gl.ts settle と同じ） */
+    function settle(vis) {
+        const lv = planLevels(vis.map(L => ({ dist: L.dist, cell: L.cell, bytes: L.levels.map(V => V.bytes) })), budget);
+        vis.forEach((L, i) => { if (L.lv !== lv[i])
+            useLevel(L, lv[i]); L.levels.forEach((V, k) => { if (k !== L.lv)
+            freeTex(V); }); ensure(L); });
+        let used = 0;
+        for (const L of layers.values())
+            for (const V of L.levels)
+                if (V.tex)
+                    used += V.bytes;
+        if (used <= budget)
+            return;
+        const on = new Set(vis);
+        for (const L of layers.values()) {
+            if (on.has(L))
+                continue;
+            for (const V of L.levels)
+                if (V.tex) {
+                    used -= V.bytes;
+                    freeTex(V);
+                }
+            if (L.tex && !L.levels.some(V => V.tex === L.tex)) {
+                L.tex = null;
+                L.bg = null;
+            }
+            if (used <= budget)
+                return;
+        }
+    }
     function rebase() { let b = 1; for (const L of layers.values()) {
         L.base = b;
         b += L.n;
@@ -347,6 +444,11 @@ export function create() {
             host?.requestDraw?.();
             return;
         }
+        if (d.type === "budget") {
+            budget = d.bytes > 0 ? d.bytes : Infinity;
+            host?.requestDraw?.();
+            return;
+        }
         if (!device)
             return;
         if (d.type === "layer") {
@@ -365,11 +467,8 @@ export function create() {
             const L = layers.get(d.q);
             if (!L)
                 return;
-            if (d.frame != null && L.texs.length > 1) {
-                L.frame = Math.min(L.texs.length - 1, Math.max(0, d.frame));
-                L.tex = L.texs[L.frame];
-                L.bg = null;
-            } // bind group は tex を抱える＝作り直す
+            if (d.frame != null)
+                L.frame = Math.max(0, d.frame); // 次に描く時に上げ直す（ensure）
             if (d.mode != null)
                 L.mode = d.mode;
             if (d.thr != null)
@@ -545,7 +644,7 @@ export function create() {
     function answerProbes(w, h) {
         if (!probes.length)
             return;
-        const ids = probes.splice(0), ls = [...layers].map(([q, L]) => ({ q, on: L.on, n: L.n, dims: L.dims, mode: L.mode, shape: L.shape, dist: +L.dist.toFixed(4), anchor: L.anchor, alt: L.alt }));
+        const ids = probes.splice(0), ls = [...layers].map(([q, L]) => ({ q, on: L.on, n: L.n, dims: L.dims, mode: L.mode, shape: L.shape, dist: +L.dist.toFixed(4), anchor: L.anchor, alt: L.alt, lv: L.lv, s: L.s, gpu: L.levels.reduce((a, V) => a + (V.tex ? V.bytes : 0), 0) }));
         readCanvas(w, h).then(data => { for (const id of ids)
             host?.post?.({ type: "pixels", id, w, h, flipY: false, layers: ls, data: id === ids[ids.length - 1] ? data : data.slice() }); })
             .catch(e => { for (const id of ids)
@@ -561,6 +660,7 @@ export function create() {
                 continue;
             vis.push(L);
         }
+        settle(vis);
         vis.sort((a, b) => b.dist - a.dist);
         const enc = device.createCommandEncoder({ label: "voxel-frame" });
         const pass = enc.beginRenderPass({ label: "voxel-draw", colorAttachments: [{ view: ctx.getCurrentTexture().createView(), loadOp: "clear", clearValue: { r: 0, g: 0, b: 0, a: 0 }, storeOp: "store" }] });
@@ -599,6 +699,7 @@ export function create() {
                 continue;
             vis.push(L);
         }
+        settle(vis);
         vis.sort((a, b) => b.dist - a.dist);
         const enc = device.createCommandEncoder({ label: "voxel-pick" });
         const pass = enc.beginRenderPass({ colorAttachments: [{ view: pickTex.createView(), loadOp: "clear", clearValue: { r: 0, g: 0, b: 0, a: 0 }, storeOp: "store" }] });

package/dist/voxel-reader.d.ts

@@ -1,16 +1,26 @@
 import type { PrepareOptions, PreparedGrid } from "./grid.js";
-/** main → worker */
+import type { SourceOptions, SourceMeta } from "./sources.js";
+/** main → worker（data＝VoxelData か、src＝URL／bytes と読み方） */
 export interface ReaderRequest {
     seq: number;
-    data: unknown;
+    data?: unknown;
+    src?: string | ArrayBuffer | Uint8Array;
+    source?: Omit<SourceOptions, "fetch">;
     opts: PrepareOptions;
 }
+/** ファイルから読んだ時に足す物 */
+export interface ReadExtra {
+    values?: Float32Array;
+    meta?: SourceMeta;
+}
 /** worker → main */
 export type ReaderReply = ({
     seq: number;
     ok: true;
-} & PreparedGrid) | {
+} & PreparedGrid & ReadExtra) | {
     seq: number;
     ok: false;
     error: string;
 };
+/** 読んで下ごしらえ（worker でも main でも同じ） */
+export declare function readAndPrepare(m: Omit<ReaderRequest, "seq">): Promise<PreparedGrid & ReadExtra>;

package/dist/voxel-reader.js

@@ -1,16 +1,19 @@
 // voxel の読み手の worker（@ortho-earth/voxel・2026-10-23）＝値の詰め（R8／RGBA8）・色表・置き方を main の外で。
 //   main（index.ts）が new Worker(new URL("./voxel-reader.js", import.meta.url), { type: "module" }) で一つ立て、層ごとに { seq, data, opts } を送る。
 //   答え＝{ seq, ok: true, …PreparedGrid（大きな格子はブリックに分けた物）}（配列は transfer）か { seq, ok: false, error }。worker が立たない器では index.ts が同じ prepareGrid を main で呼ぶ＝答えは同じ。
+//   ファイル（NetCDF／Zarr＝sources.ts）は src で送る＝取りに行くのも読むのも worker。答えに読んだ値（values・main が当たりの値に使う）と meta を足す。
 import { prepareGrid, transferables } from "./grid.js";
+import { readSource } from "./sources.js";
+/** 読んで下ごしらえ（worker でも main でも同じ） */
+export async function readAndPrepare(m) {
+    if (m.src == null)
+        return prepareGrid(m.data, m.opts);
+    const d = await readSource(m.src, m.source ?? {});
+    return { ...prepareGrid(d, m.opts), values: d.values, meta: d.meta };
+}
 const w = globalThis;
 if (typeof w.postMessage === "function" && typeof w.document === "undefined" && /Worker/.test(w.constructor?.name || "")) {
     w.onmessage = ({ data: m }) => {
-        try {
-            const p = prepareGrid(m.data, m.opts);
-            w.postMessage({ seq: m.seq, ok: true, ...p }, transferables(p));
-        }
-        catch (e) {
-            w.postMessage({ seq: m.seq, ok: false, error: String(e?.message || e) });
-        }
+        readAndPrepare(m).then(p => w.postMessage({ seq: m.seq, ok: true, ...p }, [...transferables(p), ...(p.values ? [p.values.buffer] : [])]), e => w.postMessage({ seq: m.seq, ok: false, error: String(e?.message || e) }));
     };
 }

package/package.json

@@ -1,6 +1,6 @@
 {
 	"name": "@ortho-earth/voxel",
-	"version": "0.4.0",
+	"version": "0.6.0",
 	"license": "GPL-3.0-or-later",
 	"type": "module",
 	"description": "Voxel grids for @ortho-earth/globe: fills the reserved map.addData kind \"voxel\" (VoxelData: a box or longitude/latitude/height region grid of scalar / RGBA cells), drawn as a same-frame overlay by ray marching (blocks: first solid cell, Minecraft-like; volume: trilinear, colormap, front-to-back; iso: the surface at a threshold) with clipping, scene depth, unified picking and a worker reader. WebGPU (shared device and depth) or WebGL2.",
@@ -13,7 +13,7 @@
 		"build": "tsc -p .",
 		"prepack": "node ../../scripts/pack-ts.mjs pre",
 		"postpack": "node ../../scripts/pack-ts.mjs post",
-		"test": "node tests/t-voxel.mjs"
+		"test": "node tests/t-voxel.mjs && node tests/t-sources.mjs"
 	},
 	"peerDependencies": {
 		"@ortho-earth/globe": "^2.2.0"