@holoscript/radio-astronomy-plugin 2.0.2

package/dist/fits/FITSParser.js

@@ -0,0 +1,230 @@
+/**
+ * FITSParser — Pure JavaScript parser for FITS (Flexible Image Transport System) files.
+ *
+ * FITS is the standard data format in astronomy. Structure:
+ *   - Header: 2880-byte blocks of 80-character ASCII "cards" (key=value pairs)
+ *   - Data: big-endian binary arrays, padded to 2880-byte boundary
+ *   - Extensions: additional HDUs (Header Data Units) with same structure
+ *
+ * Supports: BITPIX 8 (uint8), 16 (int16), 32 (int32), -32 (float32), -64 (float64)
+ * Handles: BSCALE/BZERO physical value scaling, WCS coordinate metadata
+ *
+ * @see https://fits.gsfc.nasa.gov/fits_standard.html
+ */
+// ── Constants ────────────────────────────────────────────────────────────────
+const BLOCK_SIZE = 2880;
+const CARD_SIZE = 80;
+const CARDS_PER_BLOCK = BLOCK_SIZE / CARD_SIZE; // 36
+// ── Parser ───────────────────────────────────────────────────────────────────
+/**
+ * Parse a FITS file from an ArrayBuffer.
+ * Returns the primary HDU (first header + data unit).
+ */
+export function parseFITS(buffer) {
+    const bytes = new Uint8Array(buffer);
+    const view = new DataView(buffer);
+    // ── Parse Header ─────────────────────────────────────────────────
+    const headers = new Map();
+    let headerEnd = 0;
+    outer: for (let block = 0; block * BLOCK_SIZE < buffer.byteLength; block++) {
+        for (let card = 0; card < CARDS_PER_BLOCK; card++) {
+            const offset = block * BLOCK_SIZE + card * CARD_SIZE;
+            if (offset + CARD_SIZE > buffer.byteLength)
+                break outer;
+            const cardStr = decodeASCII(bytes, offset, CARD_SIZE);
+            const keyword = cardStr.substring(0, 8).trim();
+            if (keyword === 'END') {
+                headerEnd = (block + 1) * BLOCK_SIZE; // data starts at next block boundary
+                break outer;
+            }
+            if (cardStr[8] === '=' && cardStr[9] === ' ') {
+                const valueStr = cardStr.substring(10).split('/')[0].trim();
+                headers.set(keyword, parseCardValue(valueStr));
+            }
+            else if (keyword === 'COMMENT' || keyword === 'HISTORY') {
+                // Skip comment/history cards
+            }
+        }
+    }
+    if (headerEnd === 0) {
+        throw new Error('FITS: No END card found in header');
+    }
+    // ── Extract Critical Keywords ────────────────────────────────────
+    const bitpix = getNum(headers, 'BITPIX');
+    const naxis = getNum(headers, 'NAXIS');
+    const shape = [];
+    for (let i = 1; i <= naxis; i++) {
+        shape.push(getNum(headers, `NAXIS${i}`));
+    }
+    const bscale = getNumOr(headers, 'BSCALE', 1.0);
+    const bzero = getNumOr(headers, 'BZERO', 0.0);
+    // ── Parse Data ───────────────────────────────────────────────────
+    const totalPixels = shape.reduce((a, b) => a * b, 1);
+    const bytesPerPixel = Math.abs(bitpix) / 8;
+    const dataOffset = headerEnd;
+    if (dataOffset + totalPixels * bytesPerPixel > buffer.byteLength) {
+        throw new Error(`FITS: Data section extends beyond buffer (need ${dataOffset + totalPixels * bytesPerPixel}, have ${buffer.byteLength})`);
+    }
+    const data = new Float32Array(totalPixels);
+    for (let i = 0; i < totalPixels; i++) {
+        const off = dataOffset + i * bytesPerPixel;
+        let raw;
+        switch (bitpix) {
+            case 8:
+                raw = bytes[off];
+                break;
+            case 16:
+                raw = view.getInt16(off, false); // big-endian
+                break;
+            case 32:
+                raw = view.getInt32(off, false);
+                break;
+            case -32:
+                raw = view.getFloat32(off, false);
+                break;
+            case -64:
+                raw = view.getFloat64(off, false);
+                break;
+            default:
+                throw new Error(`FITS: Unsupported BITPIX ${bitpix}`);
+        }
+        // Apply physical value transformation
+        data[i] = bscale * raw + bzero;
+    }
+    // ── Extract WCS ──────────────────────────────────────────────────
+    let wcs = null;
+    if (headers.has('CRVAL1')) {
+        wcs = {
+            crpix: [], crval: [], cdelt: [], ctype: [], cunit: [],
+        };
+        for (let i = 1; i <= naxis; i++) {
+            wcs.crpix.push(getNumOr(headers, `CRPIX${i}`, 1));
+            wcs.crval.push(getNumOr(headers, `CRVAL${i}`, 0));
+            wcs.cdelt.push(getNumOr(headers, `CDELT${i}`, 1));
+            wcs.ctype.push(getStrOr(headers, `CTYPE${i}`, ''));
+            wcs.cunit.push(getStrOr(headers, `CUNIT${i}`, ''));
+        }
+    }
+    return {
+        headers,
+        data,
+        shape,
+        wcs,
+        bitpix,
+        object: getStrOr(headers, 'OBJECT', ''),
+        telescope: getStrOr(headers, 'TELESCOP', ''),
+        dateObs: getStrOr(headers, 'DATE-OBS', ''),
+    };
+}
+// ── FITS Builder (for tests) ─────────────────────────────────────────────────
+/**
+ * Build a minimal FITS file as ArrayBuffer (for testing).
+ */
+export function buildFITS(opts) {
+    const cards = [];
+    cards.push(fmtCard('SIMPLE', true));
+    cards.push(fmtCard('BITPIX', opts.bitpix));
+    cards.push(fmtCard('NAXIS', opts.shape.length));
+    for (let i = 0; i < opts.shape.length; i++) {
+        cards.push(fmtCard(`NAXIS${i + 1}`, opts.shape[i]));
+    }
+    if (opts.bscale !== undefined)
+        cards.push(fmtCard('BSCALE', opts.bscale));
+    if (opts.bzero !== undefined)
+        cards.push(fmtCard('BZERO', opts.bzero));
+    if (opts.headers) {
+        for (const [key, val] of Object.entries(opts.headers)) {
+            cards.push(fmtCard(key, val));
+        }
+    }
+    cards.push('END'.padEnd(CARD_SIZE));
+    // Pad header to block boundary
+    while (cards.length % CARDS_PER_BLOCK !== 0) {
+        cards.push(' '.repeat(CARD_SIZE));
+    }
+    const headerBytes = new Uint8Array(cards.length * CARD_SIZE);
+    for (let i = 0; i < cards.length; i++) {
+        for (let j = 0; j < CARD_SIZE; j++) {
+            headerBytes[i * CARD_SIZE + j] = cards[i].charCodeAt(j);
+        }
+    }
+    // Write data
+    const bytesPerPixel = Math.abs(opts.bitpix) / 8;
+    const dataSize = opts.data.length * bytesPerPixel;
+    const paddedDataSize = Math.ceil(dataSize / BLOCK_SIZE) * BLOCK_SIZE;
+    const dataBytes = new ArrayBuffer(paddedDataSize);
+    const dataView = new DataView(dataBytes);
+    for (let i = 0; i < opts.data.length; i++) {
+        const off = i * bytesPerPixel;
+        switch (opts.bitpix) {
+            case 8:
+                new Uint8Array(dataBytes)[off] = opts.data[i];
+                break;
+            case 16:
+                dataView.setInt16(off, opts.data[i], false);
+                break;
+            case 32:
+                dataView.setInt32(off, opts.data[i], false);
+                break;
+            case -32:
+                dataView.setFloat32(off, opts.data[i], false);
+                break;
+            case -64:
+                dataView.setFloat64(off, opts.data[i], false);
+                break;
+        }
+    }
+    // Combine
+    const result = new Uint8Array(headerBytes.length + paddedDataSize);
+    result.set(headerBytes);
+    result.set(new Uint8Array(dataBytes), headerBytes.length);
+    return result.buffer;
+}
+// ── Helpers ──────────────────────────────────────────────────────────────────
+function decodeASCII(bytes, offset, length) {
+    let s = '';
+    for (let i = 0; i < length; i++) {
+        s += String.fromCharCode(bytes[offset + i]);
+    }
+    return s;
+}
+function parseCardValue(s) {
+    if (s === 'T')
+        return true;
+    if (s === 'F')
+        return false;
+    if (s.startsWith("'"))
+        return s.replace(/^'|'$/g, '').trim();
+    const n = Number(s);
+    return Number.isNaN(n) ? s : n;
+}
+function getNum(headers, key) {
+    const v = headers.get(key);
+    if (typeof v !== 'number')
+        throw new Error(`FITS: Missing or non-numeric header ${key}`);
+    return v;
+}
+function getNumOr(headers, key, def) {
+    const v = headers.get(key);
+    return typeof v === 'number' ? v : def;
+}
+function getStrOr(headers, key, def) {
+    const v = headers.get(key);
+    return typeof v === 'string' ? v : def;
+}
+function fmtCard(keyword, value) {
+    const kw = keyword.padEnd(8);
+    let valStr;
+    if (typeof value === 'boolean') {
+        valStr = value ? 'T' : 'F';
+        valStr = valStr.padStart(20);
+    }
+    else if (typeof value === 'number') {
+        valStr = String(value).padStart(20);
+    }
+    else {
+        valStr = `'${value}'`.padEnd(20);
+    }
+    return `${kw}= ${valStr}`.padEnd(CARD_SIZE);
+}
+//# sourceMappingURL=FITSParser.js.map

package/dist/fits/FITSParser.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FITSParser.js","sourceRoot":"","sources":["../../src/fits/FITSParser.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAoCH,gFAAgF;AAEhF,MAAM,UAAU,GAAG,IAAI,CAAC;AACxB,MAAM,SAAS,GAAG,EAAE,CAAC;AACrB,MAAM,eAAe,GAAG,UAAU,GAAG,SAAS,CAAC,CAAC,KAAK;AAErD,gFAAgF;AAEhF;;;GAGG;AACH,MAAM,UAAU,SAAS,CAAC,MAAmB;IAC3C,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAC;IACrC,MAAM,IAAI,GAAG,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC;IAElC,oEAAoE;IACpE,MAAM,OAAO,GAAG,IAAI,GAAG,EAAqC,CAAC;IAC7D,IAAI,SAAS,GAAG,CAAC,CAAC;IAElB,KAAK,EACL,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,UAAU,GAAG,MAAM,CAAC,UAAU,EAAE,KAAK,EAAE,EAAE,CAAC;QACpE,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,eAAe,EAAE,IAAI,EAAE,EAAE,CAAC;YAClD,MAAM,MAAM,GAAG,KAAK,GAAG,UAAU,GAAG,IAAI,GAAG,SAAS,CAAC;YACrD,IAAI,MAAM,GAAG,SAAS,GAAG,MAAM,CAAC,UAAU;gBAAE,MAAM,KAAK,CAAC;YAExD,MAAM,OAAO,GAAG,WAAW,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,CAAC,CAAC;YACtD,MAAM,OAAO,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;YAE/C,IAAI,OAAO,KAAK,KAAK,EAAE,CAAC;gBACtB,SAAS,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,qCAAqC;gBAC3E,MAAM,KAAK,CAAC;YACd,CAAC;YAED,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;gBAC7C,MAAM,QAAQ,GAAG,OAAO,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;gBAC5D,OAAO,CAAC,GAAG,CAAC,OAAO,EAAE,cAAc,CAAC,QAAQ,CAAC,CAAC,CAAC;YACjD,CAAC;iBAAM,IAAI,OAAO,KAAK,SAAS,IAAI,OAAO,KAAK,SAAS,EAAE,CAAC;gBAC1D,6BAA6B;YAC/B,CAAC;QACH,CAAC;IACH,CAAC;IAED,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;IACvD,CAAC;IAED,oEAAoE;IACpE,MAAM,MAAM,GAAG,MAAM,CAAC,OAAO,EAAE,QAAQ,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAEvC,MAAM,KAAK,GAAa,EAAE,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3C,CAAC;IAED,MAAM,MAAM,GAAG,QAAQ,CAAC,OAAO,EAAE,QAAQ,EAAE,GAAG,CAAC,CAAC;IAChD,MAAM,KAAK,GAAG,QAAQ,CAAC,OAAO,EAAE,OAAO,EAAE,GAAG,CAAC,CAAC;IAE9C,oEAAoE;IACpE,MAAM,WAAW,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IACrD,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC3C,MAAM,UAAU,GAAG,SAAS,CAAC;IAE7B,IAAI,UAAU,GAAG,WAAW,GAAG,aAAa,GAAG,MAAM,CAAC,UAAU,EAAE,CAAC;QACjE,MAAM,IAAI,KAAK,CAAC,kDAAkD,UAAU,GAAG,WAAW,GAAG,aAAa,UAAU,MAAM,CAAC,UAAU,GAAG,CAAC,CAAC;IAC5I,CAAC;IAED,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,WAAW,CAAC,CAAC;IAE3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,GAAG,GAAG,UAAU,GAAG,CAAC,GAAG,aAAa,CAAC;QAC3C,IAAI,GAAW,CAAC;QAEhB,QAAQ,MAAM,EAAE,CAAC;YACf,KAAK,CAAC;gBACJ,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;gBACjB,MAAM;YACR,KAAK,EAAE;gBACL,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,aAAa;gBAC9C,MAAM;YACR,KAAK,EAAE;gBACL,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;gBAChC,MAAM;YACR,KAAK,CAAC,EAAE;gBACN,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;gBAClC,MAAM;YACR,KAAK,CAAC,EAAE;gBACN,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;gBAClC,MAAM;YACR;gBACE,MAAM,IAAI,KAAK,CAAC,4BAA4B,MAAM,EAAE,CAAC,CAAC;QAC1D,CAAC;QAED,sCAAsC;QACtC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;IACjC,CAAC;IAED,oEAAoE;IACpE,IAAI,GAAG,GAAmB,IAAI,CAAC;IAC/B,IAAI,OAAO,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC;QAC1B,GAAG,GAAG;YACJ,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE;SACtD,CAAC;QACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;YAClD,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;YAClD,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;YAClD,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YACnD,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACrD,CAAC;IACH,CAAC;IAED,OAAO;QACL,OAAO;QACP,IAAI;QACJ,KAAK;QACL,GAAG;QACH,MAAM;QACN,MAAM,EAAE,QAAQ,CAAC,OAAO,EAAE,QAAQ,EAAE,EAAE,CAAC;QACvC,SAAS,EAAE,QAAQ,CAAC,OAAO,EAAE,UAAU,EAAE,EAAE,CAAC;QAC5C,OAAO,EAAE,QAAQ,CAAC,OAAO,EAAE,UAAU,EAAE,EAAE,CAAC;KAC3C,CAAC;AACJ,CAAC;AAED,gFAAgF;AAEhF;;GAEG;AACH,MAAM,UAAU,SAAS,CAAC,IAOzB;IACC,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC,CAAC;IACpC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;IAC3C,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;IAChD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3C,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,IAAI,IAAI,CAAC,MAAM,KAAK,SAAS;QAAE,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,QAAQ,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;IAC1E,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS;QAAE,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;IAEvE,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;QACjB,KAAK,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC;YACtD,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAChC,CAAC;IACH,CAAC;IAED,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;IAEpC,+BAA+B;IAC/B,OAAO,KAAK,CAAC,MAAM,GAAG,eAAe,KAAK,CAAC,EAAE,CAAC;QAC5C,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;IACpC,CAAC;IAED,MAAM,WAAW,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,MAAM,GAAG,SAAS,CAAC,CAAC;IAC7D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE,CAAC;YACnC,WAAW,CAAC,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QAC1D,CAAC;IACH,CAAC;IAED,aAAa;IACb,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAChD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,aAAa,CAAC;IAClD,MAAM,cAAc,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,UAAU,CAAC,GAAG,UAAU,CAAC;IACrE,MAAM,SAAS,GAAG,IAAI,WAAW,CAAC,cAAc,CAAC,CAAC;IAClD,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,SAAS,CAAC,CAAC;IAEzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;QAC9B,QAAQ,IAAI,CAAC,MAAM,EAAE,CAAC;YACpB,KAAK,CAAC;gBAAE,IAAI,UAAU,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBAAC,MAAM;YAC7D,KAAK,EAAE;gBAAE,QAAQ,CAAC,QAAQ,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;gBAAC,MAAM;YAC5D,KAAK,EAAE;gBAAE,QAAQ,CAAC,QAAQ,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;gBAAC,MAAM;YAC5D,KAAK,CAAC,EAAE;gBAAE,QAAQ,CAAC,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;gBAAC,MAAM;YAC/D,KAAK,CAAC,EAAE;gBAAE,QAAQ,CAAC,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;gBAAC,MAAM;QACjE,CAAC;IACH,CAAC;IAED,UAAU;IACV,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,WAAW,CAAC,MAAM,GAAG,cAAc,CAAC,CAAC;IACnE,MAAM,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC;IACxB,MAAM,CAAC,GAAG,CAAC,IAAI,UAAU,CAAC,SAAS,CAAC,EAAE,WAAW,CAAC,MAAM,CAAC,CAAC;IAC1D,OAAO,MAAM,CAAC,MAAM,CAAC;AACvB,CAAC;AAED,gFAAgF;AAEhF,SAAS,WAAW,CAAC,KAAiB,EAAE,MAAc,EAAE,MAAc;IACpE,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,CAAC,IAAI,MAAM,CAAC,YAAY,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IAC9C,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,cAAc,CAAC,CAAS;IAC/B,IAAI,CAAC,KAAK,GAAG;QAAE,OAAO,IAAI,CAAC;IAC3B,IAAI,CAAC,KAAK,GAAG;QAAE,OAAO,KAAK,CAAC;IAC5B,IAAI,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC;QAAE,OAAO,CAAC,CAAC,OAAO,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC;IAC7D,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACpB,OAAO,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACjC,CAAC;AAED,SAAS,MAAM,CAAC,OAA+C,EAAE,GAAW;IAC1E,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3B,IAAI,OAAO,CAAC,KAAK,QAAQ;QAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,GAAG,EAAE,CAAC,CAAC;IACzF,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,QAAQ,CAAC,OAA+C,EAAE,GAAW,EAAE,GAAW;IACzF,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3B,OAAO,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;AACzC,CAAC;AAED,SAAS,QAAQ,CAAC,OAA+C,EAAE,GAAW,EAAE,GAAW;IACzF,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3B,OAAO,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;AACzC,CAAC;AAED,SAAS,OAAO,CAAC,OAAe,EAAE,KAAgC;IAChE,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC7B,IAAI,MAAc,CAAC;IACnB,IAAI,OAAO,KAAK,KAAK,SAAS,EAAE,CAAC;QAC/B,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAC3B,MAAM,GAAG,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IAC/B,CAAC;SAAM,IAAI,OAAO,KAAK,KAAK,QAAQ,EAAE,CAAC;QACrC,MAAM,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IACtC,CAAC;SAAM,CAAC;QACN,MAAM,GAAG,IAAI,KAAK,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;IACnC,CAAC;IACD,OAAO,GAAG,EAAE,KAAK,MAAM,EAAE,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC;AAC9C,CAAC"}

package/dist/fits/FITSToGrid.d.ts

@@ -0,0 +1,27 @@
+/**
+ * FITSToGrid — Convert parsed FITS data to HoloScript grid structures.
+ *
+ * Maps FITS spectral cubes (RA × Dec × Freq) and 2D images to
+ * RegularGrid3D for visualization via ScalarFieldOverlay or SimResultsMesh.
+ */
+import { RegularGrid3D } from '@holoscript/engine/simulation';
+import type { FITSFile } from './FITSParser';
+/**
+ * Convert a parsed FITS file to a RegularGrid3D.
+ *
+ * - 3D cubes (NAXIS=3): maps directly to grid
+ * - 2D images (NAXIS=2): creates a 1-deep grid (nx × ny × 1)
+ * - 1D spectra (NAXIS=1): creates a 1×1×n grid
+ */
+export declare function fitsToGrid3D(fits: FITSFile): RegularGrid3D;
+/**
+ * Extract a single frequency channel (z-slice) from a 3D FITS cube.
+ * Returns a 2D Float32Array (nx × ny) for use as a ScalarFieldOverlay.
+ */
+export declare function extractChannel(fits: FITSFile, channel: number): Float32Array;
+/**
+ * Get data range (min/max) for the FITS data.
+ * Useful for colormap normalization.
+ */
+export declare function fitsDataRange(fits: FITSFile): [number, number];
+//# sourceMappingURL=FITSToGrid.d.ts.map

package/dist/fits/FITSToGrid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FITSToGrid.d.ts","sourceRoot":"","sources":["../../src/fits/FITSToGrid.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,aAAa,EAAE,MAAM,+BAA+B,CAAC;AAC9D,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,cAAc,CAAC;AAE7C;;;;;;GAMG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,QAAQ,GAAG,aAAa,CAmC1D;AAED;;;GAGG;AACH,wBAAgB,cAAc,CAAC,IAAI,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,GAAG,YAAY,CAsB5E;AAED;;;GAGG;AACH,wBAAgB,aAAa,CAAC,IAAI,EAAE,QAAQ,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAQ9D"}

package/dist/fits/FITSToGrid.js

@@ -0,0 +1,85 @@
+/**
+ * FITSToGrid — Convert parsed FITS data to HoloScript grid structures.
+ *
+ * Maps FITS spectral cubes (RA × Dec × Freq) and 2D images to
+ * RegularGrid3D for visualization via ScalarFieldOverlay or SimResultsMesh.
+ */
+import { RegularGrid3D } from '@holoscript/engine/simulation';
+/**
+ * Convert a parsed FITS file to a RegularGrid3D.
+ *
+ * - 3D cubes (NAXIS=3): maps directly to grid
+ * - 2D images (NAXIS=2): creates a 1-deep grid (nx × ny × 1)
+ * - 1D spectra (NAXIS=1): creates a 1×1×n grid
+ */
+export function fitsToGrid3D(fits) {
+    const shape = fits.shape;
+    let nx, ny, nz;
+    if (shape.length >= 3) {
+        nx = shape[0];
+        ny = shape[1];
+        nz = shape[2];
+    }
+    else if (shape.length === 2) {
+        nx = shape[0];
+        ny = shape[1];
+        nz = 1;
+    }
+    else if (shape.length === 1) {
+        nx = shape[0];
+        ny = 1;
+        nz = 1;
+    }
+    else {
+        throw new Error(`FITS: Cannot convert ${shape.length}D data to grid`);
+    }
+    const grid = new RegularGrid3D([nx, ny, nz], [nx, ny, nz]);
+    // FITS stores data in FORTRAN order (column-major: NAXIS1 varies fastest)
+    // RegularGrid3D uses row-major (x varies fastest in our convention)
+    // Since NAXIS1 = x and it varies fastest in both, direct copy works
+    const data = fits.data;
+    const gridData = grid.data;
+    const len = Math.min(data.length, gridData.length);
+    for (let i = 0; i < len; i++) {
+        gridData[i] = data[i];
+    }
+    return grid;
+}
+/**
+ * Extract a single frequency channel (z-slice) from a 3D FITS cube.
+ * Returns a 2D Float32Array (nx × ny) for use as a ScalarFieldOverlay.
+ */
+export function extractChannel(fits, channel) {
+    if (fits.shape.length < 3) {
+        // 2D image — return the whole thing
+        return new Float32Array(fits.data);
+    }
+    const nx = fits.shape[0];
+    const ny = fits.shape[1];
+    const nz = fits.shape[2];
+    if (channel < 0 || channel >= nz) {
+        throw new Error(`Channel ${channel} out of range [0, ${nz - 1}]`);
+    }
+    const slice = new Float32Array(nx * ny);
+    const offset = channel * nx * ny;
+    for (let i = 0; i < nx * ny; i++) {
+        slice[i] = fits.data[offset + i];
+    }
+    return slice;
+}
+/**
+ * Get data range (min/max) for the FITS data.
+ * Useful for colormap normalization.
+ */
+export function fitsDataRange(fits) {
+    let min = Infinity, max = -Infinity;
+    for (let i = 0; i < fits.data.length; i++) {
+        const v = fits.data[i];
+        if (v < min)
+            min = v;
+        if (v > max)
+            max = v;
+    }
+    return [min, max];
+}
+//# sourceMappingURL=FITSToGrid.js.map

package/dist/fits/FITSToGrid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FITSToGrid.js","sourceRoot":"","sources":["../../src/fits/FITSToGrid.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,aAAa,EAAE,MAAM,+BAA+B,CAAC;AAG9D;;;;;;GAMG;AACH,MAAM,UAAU,YAAY,CAAC,IAAc;IACzC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;IAEzB,IAAI,EAAU,EAAE,EAAU,EAAE,EAAU,CAAC;IAEvC,IAAI,KAAK,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACtB,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAChB,CAAC;SAAM,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC9B,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,CAAC,CAAC;IACT,CAAC;SAAM,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC9B,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,CAAC,CAAC;QACP,EAAE,GAAG,CAAC,CAAC;IACT,CAAC;SAAM,CAAC;QACN,MAAM,IAAI,KAAK,CAAC,wBAAwB,KAAK,CAAC,MAAM,gBAAgB,CAAC,CAAC;IACxE,CAAC;IAED,MAAM,IAAI,GAAG,IAAI,aAAa,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAE3D,0EAA0E;IAC1E,oEAAoE;IACpE,oEAAoE;IACpE,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;IACvB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC;IAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC;IAEnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,cAAc,CAAC,IAAc,EAAE,OAAe;IAC5D,IAAI,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC1B,oCAAoC;QACpC,OAAO,IAAI,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACrC,CAAC;IAED,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAEzB,IAAI,OAAO,GAAG,CAAC,IAAI,OAAO,IAAI,EAAE,EAAE,CAAC;QACjC,MAAM,IAAI,KAAK,CAAC,WAAW,OAAO,qBAAqB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACpE,CAAC;IAED,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IACxC,MAAM,MAAM,GAAG,OAAO,GAAG,EAAE,GAAG,EAAE,CAAC;IAEjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,aAAa,CAAC,IAAc;IAC1C,IAAI,GAAG,GAAG,QAAQ,EAAE,GAAG,GAAG,CAAC,QAAQ,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,GAAG;YAAE,GAAG,GAAG,CAAC,CAAC;QACrB,IAAI,CAAC,GAAG,GAAG;YAAE,GAAG,GAAG,CAAC,CAAC;IACvB,CAAC;IACD,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AACpB,CAAC"}

package/dist/index.d.ts

@@ -0,0 +1,25 @@
+import { RADIO_ASTRONOMY_TRAITS, RadioAstronomyTraitName } from './constants/astronomy-traits';
+import { PythonAstropyBridge, AstropyResult } from './bridge/python-runner';
+/**
+ * @holoscript/radio-astronomy-plugin
+ *
+ * Domain plugin bridging Radio Astrophysics simulation concepts into the HoloScript Universal pipeline.
+ * Extends standard traits without bloating core. Provides an astropy python bridge for logic evaluation.
+ */
+export { RADIO_ASTRONOMY_TRAITS, type RadioAstronomyTraitName };
+export { PythonAstropyBridge, type AstropyResult };
+export { parseFITS, buildFITS, type FITSFile, type WCSInfo } from './fits/FITSParser';
+export { fitsToGrid3D, extractChannel, fitsDataRange } from './fits/FITSToGrid';
+export { SpectralCubeViewer, FITSViewerPanel, type SpectralCubeViewerProps, type FITSViewerPanelProps } from './components/SpectralCubeViewer';
+/**
+ * Metadata exposing domain capabilities to the Studio / Schema Mapper.
+ */
+export declare const DOMAIN_MANIFEST: {
+    id: string;
+    name: string;
+    version: string;
+    description: string;
+    keywords: string[];
+    traits: readonly ["radio_emitter", "synchrotron", "interferometer", "em_wave", "pulsar_timing", "spectral_line"];
+};
+//# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,sBAAsB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAC/F,OAAO,EAAE,mBAAmB,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAE5E;;;;;GAKG;AAGH,OAAO,EAAE,sBAAsB,EAAE,KAAK,uBAAuB,EAAE,CAAC;AAGhE,OAAO,EAAE,mBAAmB,EAAE,KAAK,aAAa,EAAE,CAAC;AAGnD,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,KAAK,QAAQ,EAAE,KAAK,OAAO,EAAE,MAAM,mBAAmB,CAAC;AACtF,OAAO,EAAE,YAAY,EAAE,cAAc,EAAE,aAAa,EAAE,MAAM,mBAAmB,CAAC;AAChF,OAAO,EAAE,kBAAkB,EAAE,eAAe,EAAE,KAAK,uBAAuB,EAAE,KAAK,oBAAoB,EAAE,MAAM,iCAAiC,CAAC;AAE/I;;GAEG;AACH,eAAO,MAAM,eAAe;;;;;;;CAO3B,CAAC"}

package/dist/index.js

@@ -0,0 +1,28 @@
+import { RADIO_ASTRONOMY_TRAITS } from './constants/astronomy-traits';
+import { PythonAstropyBridge } from './bridge/python-runner';
+/**
+ * @holoscript/radio-astronomy-plugin
+ *
+ * Domain plugin bridging Radio Astrophysics simulation concepts into the HoloScript Universal pipeline.
+ * Extends standard traits without bloating core. Provides an astropy python bridge for logic evaluation.
+ */
+// Export vocabulary
+export { RADIO_ASTRONOMY_TRAITS };
+// Export Bridges
+export { PythonAstropyBridge };
+// Export FITS parsing and visualization
+export { parseFITS, buildFITS } from './fits/FITSParser';
+export { fitsToGrid3D, extractChannel, fitsDataRange } from './fits/FITSToGrid';
+export { SpectralCubeViewer, FITSViewerPanel } from './components/SpectralCubeViewer';
+/**
+ * Metadata exposing domain capabilities to the Studio / Schema Mapper.
+ */
+export const DOMAIN_MANIFEST = {
+    id: 'domain.science.astronomy.radio',
+    name: 'Radio Astronomy Plugin',
+    version: '1.0.0',
+    description: 'Extends HoloScript spatial environments with radio astrophysics primitives.',
+    keywords: ['interferometer', 'radio emitting', 'synchrotron radiation', 'pulsar'],
+    traits: RADIO_ASTRONOMY_TRAITS,
+};
+//# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,sBAAsB,EAA2B,MAAM,8BAA8B,CAAC;AAC/F,OAAO,EAAE,mBAAmB,EAAiB,MAAM,wBAAwB,CAAC;AAE5E;;;;;GAKG;AAEH,oBAAoB;AACpB,OAAO,EAAE,sBAAsB,EAAgC,CAAC;AAEhE,iBAAiB;AACjB,OAAO,EAAE,mBAAmB,EAAsB,CAAC;AAEnD,wCAAwC;AACxC,OAAO,EAAE,SAAS,EAAE,SAAS,EAA+B,MAAM,mBAAmB,CAAC;AACtF,OAAO,EAAE,YAAY,EAAE,cAAc,EAAE,aAAa,EAAE,MAAM,mBAAmB,CAAC;AAChF,OAAO,EAAE,kBAAkB,EAAE,eAAe,EAA2D,MAAM,iCAAiC,CAAC;AAE/I;;GAEG;AACH,MAAM,CAAC,MAAM,eAAe,GAAG;IAC7B,EAAE,EAAE,gCAAgC;IACpC,IAAI,EAAE,wBAAwB;IAC9B,OAAO,EAAE,OAAO;IAChB,WAAW,EAAE,6EAA6E;IAC1F,QAAQ,EAAE,CAAC,gBAAgB,EAAE,gBAAgB,EAAE,uBAAuB,EAAE,QAAQ,CAAC;IACjF,MAAM,EAAE,sBAAsB;CAC/B,CAAC"}

package/package.json

@@ -0,0 +1,33 @@
+{
+  "name": "@holoscript/radio-astronomy-plugin",
+  "version": "2.0.2",
+  "description": "Radio Astrophysics plugin for HoloScript spatial environments. Adds primitives for interferometers, pulsar analysis, and synchrotron emissions.",
+  "main": "src/index.ts",
+  "types": "src/index.ts",
+  "keywords": [
+    "holoscript",
+    "plugin",
+    "astrophysics",
+    "radio-astronomy",
+    "science"
+  ],
+  "author": "brianonbased-dev",
+  "license": "MIT",
+  "dependencies": {
+    "@holoscript/engine": "^6.1.3"
+  },
+  "peerDependencies": {
+    "@holoscript/core": "8.0.6"
+  },
+  "devDependencies": {
+    "typescript": "^5.0.0",
+    "vitest": "^4.1.5"
+  },
+  "scripts": {
+    "build": "tsc",
+    "test": "vitest run",
+    "lint": "eslint src --ext .ts",
+    "clean": "rimraf dist",
+    "test:coverage": "vitest run --coverage"
+  }
+}

package/python/astropy_bridge.py

@@ -0,0 +1,45 @@
+import sys
+import json
+import math
+
+def calculate_synchrotron(params):
+    """
+    Mock calculation for synchrotron radiation flux based on magnetic field strength.
+    In a fully operational deployment, this would use `astropy.units` and `astropy.constants`.
+    """
+    b_field = params.get('magnetic_field_gauss', 1e-4) # default interstellar field
+    frequency = params.get('frequency_hz', 1.4e9)      # default 1.4 GHz (HI line proximity)
+    
+    # Placeholder formula for demonstration (proportionality mock)
+    # Spectral index alpha ~ 0.75
+    flux = (b_field ** 2) * (frequency ** -0.75) * 1e12 
+    
+    return {
+        "flux_density_jy": flux,
+        "frequency_hz": frequency,
+        "wavelength_meters": 3e8 / frequency
+    }
+
+def main():
+    if len(sys.argv) < 3:
+        print(json.dumps({"error": "Insufficient arguments. Usage: script.py <command> <json_params>"}))
+        sys.exit(1)
+
+    command = sys.argv[1]
+    
+    try:
+        params = json.loads(sys.argv[2])
+    except Exception as e:
+        print(json.dumps({"error": f"Invalid JSON payload: {str(e)}"}))
+        sys.exit(1)
+
+    if command == "calc_synchrotron":
+        result = calculate_synchrotron(params)
+    else:
+        result = {"error": f"Unknown command: {command}"}
+
+    # Print pure JSON to stdout for the JS bridge to parse
+    print(json.dumps(result))
+
+if __name__ == "__main__":
+    main()

package/src/bridge/python-runner.ts

@@ -0,0 +1,62 @@
+import { spawn } from 'child_process';
+import * as path from 'path';
+
+/**
+ * Interface representing the result of a radio astronomy physics calculation.
+ */
+export interface AstropyResult {
+  wavelength_meters?: number;
+  frequency_hz?: number;
+  flux_density_jy?: number;
+  error?: string;
+}
+
+/**
+ * Runner that bridges HoloScript to the underlying astropy Python toolkit.
+ */
+export class PythonAstropyBridge {
+  private scriptPath: string;
+
+  constructor() {
+    this.scriptPath = path.resolve(__dirname, '../../python/astropy_bridge.py');
+  }
+
+  /**
+   * Dispatches a calculation to the Python environment.
+   * @param command The astronomical command to invoke (e.g. 'calc_synchrotron')
+   * @param params JSON payload of parameters
+   * @returns A promise resolving to the astronomical result payload
+   */
+  public async executeCommand(
+    command: string,
+    params: Record<string, any>
+  ): Promise<AstropyResult> {
+    return new Promise((resolve, reject) => {
+      const process = spawn('python', [this.scriptPath, command, JSON.stringify(params)]);
+
+      let output = '';
+      let errorOutput = '';
+
+      process.stdout.on('data', (data) => {
+        output += data.toString();
+      });
+
+      process.stderr.on('data', (data) => {
+        errorOutput += data.toString();
+      });
+
+      process.on('close', (code) => {
+        if (code !== 0) {
+          return reject(new Error(`Astropy bridge exited with code ${code}: ${errorOutput}`));
+        }
+        try {
+          // Attempt to parse JSON line output from python script
+          const parsed = JSON.parse(output.trim());
+          resolve(parsed as AstropyResult);
+        } catch (_e) {
+          reject(new Error(`Failed to parse bridge output: ${output}`));
+        }
+      });
+    });
+  }
+}