kaax-mcp 0.1.6 → 0.1.8

package/dist/gis-reports.d.ts

@@ -0,0 +1,61 @@
+/**
+ * Reports — turn one or several layers into formats the user can actually
+ * share with their team: Excel spreadsheets for the office and a single
+ * markdown summary for the agent / chat.
+ *
+ * Excel: every layer becomes a worksheet. Each feature is one row,
+ * columns are the union of every property name seen, plus columns for
+ * geometry type / WKT (Polygon, LineString, etc.) or lon/lat (Point).
+ * A first "Summary" sheet gives the agronomic totals (counts, area in
+ * hectares, length in kilometres, distinct row_ids if present).
+ *
+ * Markdown: a print-ready human-readable digest with one section per
+ * layer + a final comparison table. Stays small enough to paste into
+ * Slack or a PDF, and the agent reading it can extract numbers without
+ * a second tool call.
+ */
+import { type AnyFeatureCollection } from "./gis-io.js";
+export interface ExcelExportOptions {
+    /** Limit the number of features written — defaults to no limit. */
+    limit?: number;
+    /**
+     * If true, add a `geometry_wkt` column with WKT for non-Point geometries.
+     * Off by default because WKT for big polygons blows up cell size.
+     */
+    includeWkt?: boolean;
+    /** Override the sheet name; default 'Features'. */
+    sheetName?: string;
+}
+/**
+ * Write one or several FeatureCollections to a single XLSX file. The
+ * first sheet is `Summary`; the following sheets carry the features of
+ * each layer, one feature per row.
+ */
+export declare function exportLayersToExcel(layers: Array<{
+    name: string;
+    fc: AnyFeatureCollection;
+}>, outputPath: string, options?: ExcelExportOptions): Promise<{
+    outputPath: string;
+    sheets: number;
+    totalRows: number;
+}>;
+export interface SummaryReportOptions {
+    /** Optional title for the report. */
+    title?: string;
+    /** Add a Kaax footer link at the bottom. */
+    footer?: boolean;
+}
+/**
+ * Generate a markdown digest summarising every input layer. Lands on
+ * disk as a single .md file; small enough to paste into Slack / a PR
+ * description or attach to a PDF for the field crew.
+ */
+export declare function generateMarkdownSummary(inputs: Array<{
+    name: string;
+    path: string;
+}>, outputPath: string, options?: SummaryReportOptions): Promise<{
+    outputPath: string;
+    layers: number;
+}>;
+/** Build a friendly default name from a path. */
+export declare function defaultLayerName(path: string): string;

package/dist/gis-reports.js

@@ -0,0 +1,293 @@
+/**
+ * Reports — turn one or several layers into formats the user can actually
+ * share with their team: Excel spreadsheets for the office and a single
+ * markdown summary for the agent / chat.
+ *
+ * Excel: every layer becomes a worksheet. Each feature is one row,
+ * columns are the union of every property name seen, plus columns for
+ * geometry type / WKT (Polygon, LineString, etc.) or lon/lat (Point).
+ * A first "Summary" sheet gives the agronomic totals (counts, area in
+ * hectares, length in kilometres, distinct row_ids if present).
+ *
+ * Markdown: a print-ready human-readable digest with one section per
+ * layer + a final comparison table. Stays small enough to paste into
+ * Slack or a PDF, and the agent reading it can extract numbers without
+ * a second tool call.
+ */
+import ExcelJS from "exceljs";
+import { writeFile } from "fs/promises";
+import { basename } from "path";
+import * as turf from "@turf/turf";
+import { readAny } from "./gis-io.js";
+import { summariseLayer } from "./gis-ops.js";
+/**
+ * Write one or several FeatureCollections to a single XLSX file. The
+ * first sheet is `Summary`; the following sheets carry the features of
+ * each layer, one feature per row.
+ */
+export async function exportLayersToExcel(layers, outputPath, options = {}) {
+    const workbook = new ExcelJS.Workbook();
+    workbook.creator = "kaax-mcp";
+    workbook.created = new Date();
+    // Summary sheet — populate after we count, so leave a placeholder.
+    // exceljs accepts Partial<Column> shapes at runtime but the TS type
+    // for `worksheet.columns` is the full `Column[]`. Cast at the boundary
+    // rather than fill in every property we don't use.
+    const summary = workbook.addWorksheet("Summary");
+    summary.columns = [
+        { header: "Layer", key: "layer", width: 28 },
+        { header: "Features", key: "count", width: 11 },
+        { header: "Geometry types", key: "kinds", width: 24 },
+        { header: "Area (ha)", key: "areaHa", width: 12 },
+        { header: "Length (km)", key: "lengthKm", width: 13 },
+        { header: "Distinct row_ids", key: "rows", width: 16 },
+    ];
+    summary.getRow(1).font = { bold: true };
+    let totalRows = 0;
+    for (const { name, fc } of layers) {
+        const sheet = workbook.addWorksheet(safeSheetName(name));
+        writeFeatureSheet(sheet, fc, options);
+        const stats = layerStats(fc);
+        summary.addRow({
+            layer: name,
+            count: stats.count,
+            kinds: Object.entries(stats.byGeometry)
+                .map(([k, v]) => `${k}:${v}`)
+                .join(", "),
+            areaHa: stats.totalAreaHa ? Number(stats.totalAreaHa.toFixed(3)) : null,
+            lengthKm: stats.totalLengthM
+                ? Number((stats.totalLengthM / 1000).toFixed(3))
+                : null,
+            rows: stats.distinctRows ?? null,
+        });
+        totalRows += stats.count;
+    }
+    // Tidy the summary header style.
+    summary.getRow(1).alignment = { vertical: "middle", horizontal: "left" };
+    await workbook.xlsx.writeFile(outputPath);
+    return { outputPath, sheets: layers.length, totalRows };
+}
+function writeFeatureSheet(sheet, fc, options) {
+    const limit = options.limit ?? fc.features.length;
+    const features = fc.features.slice(0, limit);
+    // Collect the union of property names across all features (sorted).
+    const propNames = new Set();
+    for (const f of features) {
+        for (const k of Object.keys(f.properties ?? {}))
+            propNames.add(k);
+    }
+    const sortedProps = [...propNames].sort();
+    // Columns: properties + geometry summary. Same Partial-vs-full Column
+    // shrug as the Summary sheet — cast at the boundary.
+    const columns = [
+        { header: "feature_index", key: "feature_index", width: 14 },
+        { header: "geometry_type", key: "geometry_type", width: 16 },
+        { header: "lon", key: "lon", width: 12 },
+        { header: "lat", key: "lat", width: 12 },
+        ...sortedProps.map((p) => ({ header: p, key: `p_${p}`, width: 18 })),
+    ];
+    if (options.includeWkt) {
+        columns.push({ header: "geometry_wkt", key: "wkt", width: 80 });
+    }
+    sheet.columns = columns;
+    sheet.getRow(1).font = { bold: true };
+    features.forEach((f, idx) => {
+        const row = {
+            feature_index: idx,
+            geometry_type: f.geometry?.type ?? "(none)",
+        };
+        const centroid = geometryRepresentativePoint(f);
+        if (centroid) {
+            row.lon = centroid[0];
+            row.lat = centroid[1];
+        }
+        for (const p of sortedProps) {
+            const v = f.properties?.[p];
+            // Excel can store strings / numbers / booleans / Date. Coerce
+            // everything else to JSON so the cell never carries [Object].
+            row[`p_${p}`] =
+                v === null || v === undefined
+                    ? null
+                    : typeof v === "object"
+                        ? JSON.stringify(v)
+                        : v;
+        }
+        if (options.includeWkt && f.geometry) {
+            row.wkt = geometryToWKT(f.geometry);
+        }
+        sheet.addRow(row);
+    });
+    // Freeze the header row.
+    sheet.views = [{ state: "frozen", ySplit: 1 }];
+}
+/** Produce a small representative point (lon, lat) for a feature. */
+function geometryRepresentativePoint(f) {
+    if (!f.geometry)
+        return null;
+    try {
+        if (f.geometry.type === "Point") {
+            const c = f.geometry.coordinates;
+            return [c[0], c[1]];
+        }
+        const ctr = turf.pointOnFeature(f);
+        const c = ctr.geometry.coordinates;
+        return [c[0], c[1]];
+    }
+    catch {
+        return null;
+    }
+}
+/**
+ * Minimal WKT serialiser — enough for QGIS-style import / paste-into-
+ * PostGIS. Skips Z values; Kaax data is 2-D.
+ */
+function geometryToWKT(g) {
+    const ring = (coords) => `(${coords.map((c) => `${c[0]} ${c[1]}`).join(", ")})`;
+    const poly = (rings) => `(${rings.map((r) => ring(r)).join(", ")})`;
+    switch (g.type) {
+        case "Point": {
+            const c = g.coordinates;
+            return `POINT (${c[0]} ${c[1]})`;
+        }
+        case "MultiPoint":
+            return `MULTIPOINT (${g.coordinates
+                .map((c) => `${c[0]} ${c[1]}`)
+                .join(", ")})`;
+        case "LineString":
+            return `LINESTRING ${ring(g.coordinates)}`;
+        case "MultiLineString":
+            return `MULTILINESTRING (${g.coordinates
+                .map((ls) => ring(ls))
+                .join(", ")})`;
+        case "Polygon":
+            return `POLYGON ${poly(g.coordinates)}`;
+        case "MultiPolygon":
+            return `MULTIPOLYGON (${g.coordinates
+                .map((p) => poly(p))
+                .join(", ")})`;
+        default:
+            return `(${g.type})`;
+    }
+}
+function safeSheetName(raw) {
+    // Excel sheet names: max 31 chars, no : \ / ? * [ ]
+    let s = raw.replace(/[:\\/?*[\]]/g, "_").slice(0, 31);
+    if (!s)
+        s = "Sheet";
+    return s;
+}
+function layerStats(fc) {
+    const base = summariseLayer(fc);
+    const rowIds = new Set();
+    let hasRowId = false;
+    for (const f of fc.features) {
+        const rid = f.properties?.row_id;
+        if (typeof rid === "number") {
+            hasRowId = true;
+            rowIds.add(rid);
+        }
+    }
+    let bbox;
+    try {
+        const b = turf.bbox(fc);
+        bbox = [b[0], b[1], b[2], b[3]];
+    }
+    catch {
+        /* skip — empty or bad geometry */
+    }
+    return {
+        ...base,
+        ...(hasRowId ? { distinctRows: rowIds.size } : {}),
+        ...(bbox ? { bbox } : {}),
+    };
+}
+/**
+ * Generate a markdown digest summarising every input layer. Lands on
+ * disk as a single .md file; small enough to paste into Slack / a PR
+ * description or attach to a PDF for the field crew.
+ */
+export async function generateMarkdownSummary(inputs, outputPath, options = {}) {
+    const sections = [];
+    sections.push(`# ${options.title ?? "Kaax MCP — Layer Report"}`);
+    sections.push("");
+    sections.push(`Generated ${new Date().toISOString()}`);
+    sections.push("");
+    sections.push(`Layers analysed: **${inputs.length}**`);
+    sections.push("");
+    // Per-layer detail
+    const allStats = [];
+    for (const inp of inputs) {
+        sections.push(`## ${inp.name}`);
+        sections.push("");
+        sections.push(`Source: \`${inp.path}\``);
+        try {
+            const fc = await readAny(inp.path);
+            const s = layerStats(fc);
+            allStats.push({ name: inp.name, stats: s, path: inp.path });
+            sections.push("");
+            sections.push(`- **Features**: ${s.count}`);
+            sections.push(`- **Geometry types**: ${Object.entries(s.byGeometry)
+                .map(([k, v]) => `${k} × ${v}`)
+                .join(", ") || "—"}`);
+            if (s.totalAreaHa !== undefined) {
+                sections.push(`- **Total area**: ${s.totalAreaHa.toFixed(2)} ha`);
+            }
+            if (s.totalLengthM !== undefined) {
+                sections.push(`- **Total length**: ${(s.totalLengthM / 1000).toFixed(2)} km`);
+            }
+            if (s.distinctRows !== undefined) {
+                sections.push(`- **Rows (distinct row_id)**: ${s.distinctRows}`);
+            }
+            if (s.bbox) {
+                sections.push(`- **Bounding box**: ${s.bbox
+                    .map((n) => n.toFixed(5))
+                    .join(", ")} (minLon, minLat, maxLon, maxLat)`);
+            }
+        }
+        catch (err) {
+            sections.push("");
+            sections.push(`> ⚠️ Failed to read this layer: ${err.message}`);
+        }
+        sections.push("");
+    }
+    // Comparison table if there are 2+ layers.
+    if (allStats.length >= 2) {
+        sections.push("## Comparison");
+        sections.push("");
+        sections.push("| Layer | Features | Area (ha) | Length (km) | Rows |");
+        sections.push("|---|---:|---:|---:|---:|");
+        for (const s of allStats) {
+            const cells = [
+                s.name,
+                String(s.stats.count),
+                s.stats.totalAreaHa !== undefined
+                    ? s.stats.totalAreaHa.toFixed(2)
+                    : "—",
+                s.stats.totalLengthM !== undefined
+                    ? (s.stats.totalLengthM / 1000).toFixed(2)
+                    : "—",
+                s.stats.distinctRows !== undefined
+                    ? String(s.stats.distinctRows)
+                    : "—",
+            ];
+            sections.push(`| ${cells.join(" | ")} |`);
+        }
+        sections.push("");
+    }
+    if (options.footer !== false) {
+        sections.push("---");
+        sections.push("");
+        sections.push("Generated by [kaax-mcp](https://www.npmjs.com/package/kaax-mcp).");
+    }
+    const md = sections.join("\n");
+    await writeFile(outputPath, md, "utf8");
+    return { outputPath, layers: inputs.length };
+}
+// ─────────────────────────────────────────────────────────────────────────
+// Helpers for tools.ts
+// ─────────────────────────────────────────────────────────────────────────
+/** Build a friendly default name from a path. */
+export function defaultLayerName(path) {
+    return basename(path, /\.[^.]+$/.exec(path)?.[0] ?? "");
+}
+//# sourceMappingURL=gis-reports.js.map

package/dist/gis-reports.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gis-reports.js","sourceRoot":"","sources":["../src/gis-reports.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,OAAO,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,QAAQ,EAAE,MAAM,MAAM,CAAC;AAKhC,OAAO,KAAK,IAAI,MAAM,YAAY,CAAC;AACnC,OAAO,EAAE,OAAO,EAA6B,MAAM,aAAa,CAAC;AACjE,OAAO,EAAE,cAAc,EAAE,MAAM,cAAc,CAAC;AAkB9C;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,mBAAmB,CACvC,MAAyD,EACzD,UAAkB,EAClB,UAA8B,EAAE;IAEhC,MAAM,QAAQ,GAAG,IAAI,OAAO,CAAC,QAAQ,EAAE,CAAC;IACxC,QAAQ,CAAC,OAAO,GAAG,UAAU,CAAC;IAC9B,QAAQ,CAAC,OAAO,GAAG,IAAI,IAAI,EAAE,CAAC;IAE9B,mEAAmE;IACnE,oEAAoE;IACpE,uEAAuE;IACvE,mDAAmD;IACnD,MAAM,OAAO,GAAG,QAAQ,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;IACjD,OAAO,CAAC,OAAO,GAAG;QAChB,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,KAAK,EAAE,EAAE,EAAE;QAC5C,EAAE,MAAM,EAAE,UAAU,EAAE,GAAG,EAAE,OAAO,EAAE,KAAK,EAAE,EAAE,EAAE;QAC/C,EAAE,MAAM,EAAE,gBAAgB,EAAE,GAAG,EAAE,OAAO,EAAE,KAAK,EAAE,EAAE,EAAE;QACrD,EAAE,MAAM,EAAE,WAAW,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,EAAE,EAAE;QACjD,EAAE,MAAM,EAAE,aAAa,EAAE,GAAG,EAAE,UAAU,EAAE,KAAK,EAAE,EAAE,EAAE;QACrD,EAAE,MAAM,EAAE,kBAAkB,EAAE,GAAG,EAAE,MAAM,EAAE,KAAK,EAAE,EAAE,EAAE;KACxB,CAAC;IACjC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC;IAExC,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,MAAM,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,MAAM,EAAE,CAAC;QAClC,MAAM,KAAK,GAAG,QAAQ,CAAC,YAAY,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC,CAAC;QACzD,iBAAiB,CAAC,KAAK,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC;QACtC,MAAM,KAAK,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC;QAC7B,OAAO,CAAC,MAAM,CAAC;YACb,KAAK,EAAE,IAAI;YACX,KAAK,EAAE,KAAK,CAAC,KAAK;YAClB,KAAK,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,UAAU,CAAC;iBACpC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;iBAC5B,IAAI,CAAC,IAAI,CAAC;YACb,MAAM,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC,KAAK,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI;YACvE,QAAQ,EAAE,KAAK,CAAC,YAAY;gBAC1B,CAAC,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;gBAChD,CAAC,CAAC,IAAI;YACR,IAAI,EAAE,KAAK,CAAC,YAAY,IAAI,IAAI;SACjC,CAAC,CAAC;QACH,SAAS,IAAI,KAAK,CAAC,KAAK,CAAC;IAC3B,CAAC;IAED,iCAAiC;IACjC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,EAAE,QAAQ,EAAE,QAAQ,EAAE,UAAU,EAAE,MAAM,EAAE,CAAC;IAEzE,MAAM,QAAQ,CAAC,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,CAAC;IAC1C,OAAO,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC;AAC1D,CAAC;AAED,SAAS,iBAAiB,CACxB,KAAwB,EACxB,EAAwB,EACxB,OAA2B;IAE3B,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,IAAI,EAAE,CAAC,QAAQ,CAAC,MAAM,CAAC;IAClD,MAAM,QAAQ,GAAG,EAAE,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAC7C,oEAAoE;IACpE,MAAM,SAAS,GAAG,IAAI,GAAG,EAAU,CAAC;IACpC,KAAK,MAAM,CAAC,IAAI,QAAQ,EAAE,CAAC;QACzB,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,UAAU,IAAI,EAAE,CAAC;YAAE,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpE,CAAC;IACD,MAAM,WAAW,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,IAAI,EAAE,CAAC;IAE1C,sEAAsE;IACtE,qDAAqD;IACrD,MAAM,OAAO,GAAG;QACd,EAAE,MAAM,EAAE,eAAe,EAAE,GAAG,EAAE,eAAe,EAAE,KAAK,EAAE,EAAE,EAAE;QAC5D,EAAE,MAAM,EAAE,eAAe,EAAE,GAAG,EAAE,eAAe,EAAE,KAAK,EAAE,EAAE,EAAE;QAC5D,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE;QACxC,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE;QACxC,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,EAAE,KAAK,CAAC,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC,CAAC;KACrE,CAAC;IACF,IAAI,OAAO,CAAC,UAAU,EAAE,CAAC;QACvB,OAAO,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,cAAc,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC,CAAC;IAClE,CAAC;IACD,KAAK,CAAC,OAAO,GAAG,OAAsC,CAAC;IACvD,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC;IAEtC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE;QAC1B,MAAM,GAAG,GAA4B;YACnC,aAAa,EAAE,GAAG;YAClB,aAAa,EAAE,CAAC,CAAC,QAAQ,EAAE,IAAI,IAAI,QAAQ;SAC5C,CAAC;QACF,MAAM,QAAQ,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;QAChD,IAAI,QAAQ,EAAE,CAAC;YACb,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;YACtB,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACxB,CAAC;QACD,KAAK,MAAM,CAAC,IAAI,WAAW,EAAE,CAAC;YAC5B,MAAM,CAAC,GAAI,CAAC,CAAC,UAAkD,EAAE,CAAC,CAAC,CAAC,CAAC;YACrE,8DAA8D;YAC9D,8DAA8D;YAC9D,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;gBACX,CAAC,KAAK,IAAI,IAAI,CAAC,KAAK,SAAS;oBAC3B,CAAC,CAAC,IAAI;oBACN,CAAC,CAAC,OAAO,CAAC,KAAK,QAAQ;wBACrB,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;wBACnB,CAAC,CAAE,CAA+B,CAAC;QAC3C,CAAC;QACD,IAAI,OAAO,CAAC,UAAU,IAAI,CAAC,CAAC,QAAQ,EAAE,CAAC;YACrC,GAAG,CAAC,GAAG,GAAG,aAAa,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QACtC,CAAC;QACD,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC;IAEH,yBAAyB;IACzB,KAAK,CAAC,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC;AACjD,CAAC;AAED,qEAAqE;AACrE,SAAS,2BAA2B,CAClC,CAA6C;IAE7C,IAAI,CAAC,CAAC,CAAC,QAAQ;QAAE,OAAO,IAAI,CAAC;IAC7B,IAAI,CAAC;QACH,IAAI,CAAC,CAAC,QAAQ,CAAC,IAAI,KAAK,OAAO,EAAE,CAAC;YAChC,MAAM,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,WAA+B,CAAC;YACrD,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtB,CAAC;QACD,MAAM,GAAG,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,WAA+B,CAAC;QACvD,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IAAC,MAAM,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;AACH,CAAC;AAED;;;GAGG;AACH,SAAS,aAAa,CAAC,CAAW;IAChC,MAAM,IAAI,GAAG,CAAC,MAAkB,EAAE,EAAE,CAClC,IAAI,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC;IACzD,MAAM,IAAI,GAAG,CAAC,KAAmB,EAAE,EAAE,CACnC,IAAI,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC;IAC9C,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC;QACf,KAAK,OAAO,CAAC,CAAC,CAAC;YACb,MAAM,CAAC,GAAG,CAAC,CAAC,WAAuB,CAAC;YACpC,OAAO,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACnC,CAAC;QACD,KAAK,YAAY;YACf,OAAO,eAAgB,CAAC,CAAC,WAA0B;iBAChD,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;iBAC7B,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC;QACnB,KAAK,YAAY;YACf,OAAO,cAAc,IAAI,CAAC,CAAC,CAAC,WAAyB,CAAC,EAAE,CAAC;QAC3D,KAAK,iBAAiB;YACpB,OAAO,oBAAqB,CAAC,CAAC,WAA4B;iBACvD,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;iBACrB,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC;QACnB,KAAK,SAAS;YACZ,OAAO,WAAW,IAAI,CAAC,CAAC,CAAC,WAA2B,CAAC,EAAE,CAAC;QAC1D,KAAK,cAAc;YACjB,OAAO,iBAAkB,CAAC,CAAC,WAA8B;iBACtD,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;iBACnB,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC;QACnB;YACE,OAAO,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC;IACzB,CAAC;AACH,CAAC;AAED,SAAS,aAAa,CAAC,GAAW;IAChC,oDAAoD;IACpD,IAAI,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,cAAc,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IACtD,IAAI,CAAC,CAAC;QAAE,CAAC,GAAG,OAAO,CAAC;IACpB,OAAO,CAAC,CAAC;AACX,CAAC;AAiBD,SAAS,UAAU,CAAC,EAAwB;IAC1C,MAAM,IAAI,GAAG,cAAc,CAAC,EAAE,CAAC,CAAC;IAChC,MAAM,MAAM,GAAG,IAAI,GAAG,EAAU,CAAC;IACjC,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC,QAAQ,EAAE,CAAC;QAC5B,MAAM,GAAG,GAAI,CAAC,CAAC,UAA0C,EAAE,MAAM,CAAC;QAClE,IAAI,OAAO,GAAG,KAAK,QAAQ,EAAE,CAAC;YAC5B,QAAQ,GAAG,IAAI,CAAC;YAChB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IACD,IAAI,IAAoC,CAAC;IACzC,IAAI,CAAC;QACH,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACxB,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;IACtC,CAAC;IAAC,MAAM,CAAC;QACP,kCAAkC;IACpC,CAAC;IACD,OAAO;QACL,GAAG,IAAI;QACP,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAClD,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;KAC1B,CAAC;AACJ,CAAC;AAaD;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,uBAAuB,CAC3C,MAA6C,EAC7C,UAAkB,EAClB,UAAgC,EAAE;IAElC,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,QAAQ,CAAC,IAAI,CAAC,KAAK,OAAO,CAAC,KAAK,IAAI,yBAAyB,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAClB,QAAQ,CAAC,IAAI,CAAC,aAAa,IAAI,IAAI,EAAE,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC;IACvD,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAClB,QAAQ,CAAC,IAAI,CAAC,sBAAsB,MAAM,CAAC,MAAM,IAAI,CAAC,CAAC;IACvD,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAElB,mBAAmB;IACnB,MAAM,QAAQ,GAA6D,EAAE,CAAC;IAC9E,KAAK,MAAM,GAAG,IAAI,MAAM,EAAE,CAAC;QACzB,QAAQ,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC;QAChC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAClB,QAAQ,CAAC,IAAI,CAAC,aAAa,GAAG,CAAC,IAAI,IAAI,CAAC,CAAC;QACzC,IAAI,CAAC;YACH,MAAM,EAAE,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YACnC,MAAM,CAAC,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC;YACzB,QAAQ,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC;YAC5D,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClB,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;YAC5C,QAAQ,CAAC,IAAI,CACX,yBACE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,UAAU,CAAC;iBACzB,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC;iBAC9B,IAAI,CAAC,IAAI,CAAC,IAAI,GACnB,EAAE,CACH,CAAC;YACF,IAAI,CAAC,CAAC,WAAW,KAAK,SAAS,EAAE,CAAC;gBAChC,QAAQ,CAAC,IAAI,CAAC,qBAAqB,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;YACpE,CAAC;YACD,IAAI,CAAC,CAAC,YAAY,KAAK,SAAS,EAAE,CAAC;gBACjC,QAAQ,CAAC,IAAI,CACX,uBAAuB,CAAC,CAAC,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAC/D,CAAC;YACJ,CAAC;YACD,IAAI,CAAC,CAAC,YAAY,KAAK,SAAS,EAAE,CAAC;gBACjC,QAAQ,CAAC,IAAI,CAAC,iCAAiC,CAAC,CAAC,YAAY,EAAE,CAAC,CAAC;YACnE,CAAC;YACD,IAAI,CAAC,CAAC,IAAI,EAAE,CAAC;gBACX,QAAQ,CAAC,IAAI,CACX,uBAAuB,CAAC,CAAC,IAAI;qBAC1B,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;qBACxB,IAAI,CAAC,IAAI,CAAC,mCAAmC,CACjD,CAAC;YACJ,CAAC;QACH,CAAC;QAAC,OAAO,GAAG,EAAE,CAAC;YACb,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClB,QAAQ,CAAC,IAAI,CACX,mCAAoC,GAAa,CAAC,OAAO,EAAE,CAC5D,CAAC;QACJ,CAAC;QACD,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,CAAC;IAED,2CAA2C;IAC3C,IAAI,QAAQ,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACzB,QAAQ,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC;QAC/B,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAClB,QAAQ,CAAC,IAAI,CACX,uDAAuD,CACxD,CAAC;QACF,QAAQ,CAAC,IAAI,CAAC,2BAA2B,CAAC,CAAC;QAC3C,KAAK,MAAM,CAAC,IAAI,QAAQ,EAAE,CAAC;YACzB,MAAM,KAAK,GAAG;gBACZ,CAAC,CAAC,IAAI;gBACN,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC;gBACrB,CAAC,CAAC,KAAK,CAAC,WAAW,KAAK,SAAS;oBAC/B,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC;oBAChC,CAAC,CAAC,GAAG;gBACP,CAAC,CAAC,KAAK,CAAC,YAAY,KAAK,SAAS;oBAChC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;oBAC1C,CAAC,CAAC,GAAG;gBACP,CAAC,CAAC,KAAK,CAAC,YAAY,KAAK,SAAS;oBAChC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,YAAY,CAAC;oBAC9B,CAAC,CAAC,GAAG;aACR,CAAC;YACF,QAAQ,CAAC,IAAI,CAAC,KAAK,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC5C,CAAC;QACD,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,CAAC;IAED,IAAI,OAAO,CAAC,MAAM,KAAK,KAAK,EAAE,CAAC;QAC7B,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACrB,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAClB,QAAQ,CAAC,IAAI,CAAC,kEAAkE,CAAC,CAAC;IACpF,CAAC;IAED,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC/B,MAAM,SAAS,CAAC,UAAU,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;IACxC,OAAO,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,CAAC;AAC/C,CAAC;AAED,4EAA4E;AAC5E,uBAAuB;AACvB,4EAA4E;AAE5E,iDAAiD;AACjD,MAAM,UAAU,gBAAgB,CAAC,IAAY;IAC3C,OAAO,QAAQ,CAAC,IAAI,EAAE,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;AAC1D,CAAC"}

package/dist/gis-rows.d.ts

@@ -0,0 +1,122 @@
+/**
+ * Row reconstruction from point clouds — the agricultural-specific GIS
+ * pipeline. Given a SHP of plant detections, build the surco geometry:
+ *
+ *   1. Detect the dominant orientation of the field (histogram of pairwise
+ *      bearings between nearby points — surcos are colinear, so the
+ *      bearing distribution has a clear modal value).
+ *   2. Cluster every point into a row id, by projecting onto the axis
+ *      perpendicular to the orientation and 1-D DBSCAN-ing the projection.
+ *   3. Generate one or more LineStrings per row, with **gap detection**:
+ *      consecutive points within N × median(spacing) get connected (these
+ *      are gaps caused by missing / dead plants — depopulation). Gaps
+ *      larger than that are presumed to be physical streets / roads
+ *      cutting the surco, so the line is split into two segments and the
+ *      gap can be emitted as a `street` feature for the user to inspect.
+ *
+ * Everything runs in pure JS — turf for haversine distance + bearing, a
+ * small local equirectangular projection for the 1-D math (constant
+ * scaling factor against the centroid latitude — accurate enough for any
+ * single field, much faster than running every cluster step through
+ * great-circle geometry).
+ *
+ * Limits:
+ *   - Pairwise bearing computation is O(n²). For n > 8 000 points we
+ *     sample randomly to keep wall-clock under ~3 s on a laptop. The
+ *     final clustering uses every point.
+ *   - Geometry assumes lon/lat (EPSG:4326). For larger geographies the
+ *     equirectangular approximation breaks down; in practice farm fields
+ *     are < 5 km wide so the error is sub-meter.
+ */
+import type { AnyFeatureCollection } from "./gis-io.js";
+export interface OrientationResult {
+    /** Dominant row angle in degrees, wrapped to [0, 180). */
+    angleDeg: number;
+    /** Same value in radians. */
+    angleRad: number;
+    /**
+     * Fraction of the considered point-pairs whose bearing fell in the
+     * dominant bin. Useful to decide whether the field has a clear
+     * orientation (>0.4 is strong, 0.2 is weak / spotty).
+     */
+    confidence: number;
+    /** Total pairs considered after distance filter and sampling. */
+    totalPairsConsidered: number;
+    /** Width of each histogram bin in degrees (5° by default). */
+    binSizeDeg: number;
+}
+interface DetectOptions {
+    /**
+     * Maximum distance in meters between two points to count them as a
+     * "candidate row neighbour". Default 5 m — most agricultural rows are
+     * 0.5–3 m apart with detections every 0.3–1.5 m along the row, so 5 m
+     * captures intra-row pairs but not unrelated points across the field.
+     */
+    maxNeighborDistanceM?: number;
+    /** Histogram bin size in degrees. Default 5°. */
+    binSizeDeg?: number;
+    /**
+     * If point count exceeds this, sample down before computing pairwise
+     * bearings. The dominant angle is statistically stable from a few
+     * thousand pairs — no need to compute all of O(n²).
+     */
+    maxPointsForPairwise?: number;
+}
+export declare function detectRowOrientation(fc: AnyFeatureCollection, options?: DetectOptions): OrientationResult;
+export interface ClusterOptions {
+    /** Pre-computed orientation. If absent, we detect it. */
+    angleRad?: number;
+    /**
+     * Expected row spacing in meters. If absent, we estimate it as the
+     * median of nearest perpendicular-distance gaps.
+     */
+    rowSpacingM?: number;
+    /**
+     * Multiplier applied to rowSpacingM to define the "max gap between
+     * adjacent clusters". 0.55 by default — anything more than ~half a row
+     * spacing is treated as a fresh cluster.
+     */
+    clusterGapFraction?: number;
+}
+/**
+ * Return a copy of `fc` where every Point feature carries a new `row_id`
+ * property (1-indexed integer). Non-point features are skipped.
+ */
+export declare function clusterPointsByRow(fc: AnyFeatureCollection, options?: ClusterOptions): AnyFeatureCollection;
+export interface RowLinesOptions {
+    /**
+     * Gaps larger than this multiplier × median(spacing within a row) are
+     * treated as streets / roads and split the row. Default: 3 — empirical
+     * sweet spot from sugarcane and pineapple datasets.
+     */
+    gapMultiplier?: number;
+    /** Also emit a separate FeatureCollection of detected street gaps. */
+    includeStreets?: boolean;
+    /**
+     * Optional pre-computed orientation; otherwise we detect it from the
+     * same points. Used to project for spacing computation.
+     */
+    angleRad?: number;
+}
+export interface RowLinesResult {
+    /** LineString features, one per contiguous row segment. */
+    lines: AnyFeatureCollection;
+    /** Optional point-pair features marking detected street gaps. */
+    streets?: AnyFeatureCollection;
+    /** Stats — useful to render as a tool output summary. */
+    stats: {
+        rows: number;
+        segments: number;
+        streetGapsDetected: number;
+        totalLengthM: number;
+    };
+}
+/**
+ * Build LineStrings from a clustered point set. Expects each point to
+ * carry a `row_id` property (output of `clusterPointsByRow`). Each row
+ * produces one or more LineStrings depending on whether any large gap
+ * was detected inside it — those gaps are presumed to be streets, not
+ * depopulation, and split the row at that point.
+ */
+export declare function generateRowLines(fc: AnyFeatureCollection, options?: RowLinesOptions): RowLinesResult;
+export {};