[{"_1":2,"_3":-5,"_4":-5},"loaderData",{"_5":6,"_7":8,"_9":10},"actionData","errors","root",{"_204":205,"_206":207},"routes/_app.x.$handle.$name",{"_97":96,"_32":45,"_105":87,"_27":106,"_13":14,"_107":-5,"_11":108,"_109":110,"_111":112,"_113":114,"_115":116,"_20":-5,"_117":-5,"_118":-5,"_119":25,"_120":121,"_122":59,"_123":-5,"_124":-5,"_125":126},"routes/_app.x.$handle.$name.branch.$branch.version.$version.diff",{"_11":12,"_13":14,"_15":16,"_17":-5,"_18":-5},"project",{"_43":44,"_32":45,"_38":46,"_47":-5,"_48":25,"_49":50,"_51":50,"_52":53,"_54":-5,"_18":55,"_56":57,"_58":59,"_60":-5,"_61":-5,"_62":63,"_64":65,"_66":-5,"_67":-5,"_68":-5,"_69":-5,"_70":71,"_72":71,"_73":74,"_75":74,"_76":74,"_77":71,"_78":71,"_79":80,"_81":82,"_83":84,"_85":25},"version",33,"changes",[19],"aiChatInfo","userId",{"_20":21,"_22":23,"_24":25,"_26":25,"_27":28,"_29":28},"file",{"_30":31,"_32":28,"_33":34,"_35":37,"_38":39,"_40":14,"_41":42,"_24":25},"previous",{"_30":31,"_32":28,"_33":34,"_35":36},"deleted",false,"deletedHere","path","colormaps.ts","latestPath","fileId","a0962790-869c-11f1-93a2-1607ee4eb77e","name","fileType","script","content","export type Stop = [value: number, hex: string];\n\nexport interface Palette {\n label: string;\n unit: string;\n stops: Stop[];\n /** valori sotto questa soglia = trasparenti (rumore/no eco) */\n floor?: number;\n /** se true, i limiti sono +/- NI letto dal file */\n nyquist?: boolean;\n}\n\nexport const PALETTES: Record = {\n DBZH: {\n label: \"Riflettivita'\",\n unit: \"dBZ\",\n floor: 5,\n stops: [\n [5, \"#04e9e7\"], [10, \"#019ff4\"], [15, \"#0300f4\"], [20, \"#02fd02\"],\n [25, \"#01c501\"], [30, \"#008e00\"], [35, \"#fdf802\"], [40, \"#e5bc00\"],\n [45, \"#fd9500\"], [50, \"#fd0000\"], [55, \"#d40000\"], [60, \"#bc0000\"],\n [65, \"#f800fd\"], [70, \"#9854c6\"],\n ],\n },\n TH: {\n label: \"Riflettivita' totale\",\n unit: \"dBZ\",\n floor: 5,\n stops: [\n [5, \"#04e9e7\"], [10, \"#019ff4\"], [15, \"#0300f4\"], [20, \"#02fd02\"],\n [25, \"#01c501\"], [30, \"#008e00\"], [35, \"#fdf802\"], [40, \"#e5bc00\"],\n [45, \"#fd9500\"], [50, \"#fd0000\"], [55, \"#d40000\"], [60, \"#bc0000\"],\n [65, \"#f800fd\"], [70, \"#9854c6\"],\n ],\n },\n RHOHV: {\n label: \"Coefficiente di correlazione\",\n unit: \"\",\n stops: [\n [0.20, \"#2b2b3d\"], [0.65, \"#3b4cc0\"], [0.80, \"#4a9bd6\"],\n [0.90, \"#6bd0d5\"], [0.95, \"#b0f27a\"], [0.97, \"#ffd24d\"],\n [0.99, \"#ff6b4a\"], [1.05, \"#d21f3c\"],\n ],\n },\n VRADH: {\n label: \"Velocita' radiale\",\n unit: \"m/s\",\n nyquist: true,\n stops: [\n [-1.0, \"#00ff88\"], [-0.60, \"#00c853\"], [-0.25, \"#1b7a3d\"],\n [-0.04, \"#123322\"], [0.04, \"#33161a\"], [0.25, \"#a33131\"],\n [0.60, \"#e53935\"], [1.0, \"#ff8a80\"],\n ],\n },\n WRADH: {\n label: \"Ampiezza spettrale\",\n unit: \"m/s\",\n stops: [\n [0, \"#101820\"], [1, \"#2c5f8a\"], [2, \"#3ba99c\"],\n [3, \"#8fd14f\"], [4, \"#ffd166\"], [6, \"#ef476f\"],\n ],\n },\n ZDR: {\n label: \"Riflettivita' differenziale\",\n unit: \"dB\",\n stops: [\n [-2, \"#6a3d9a\"], [-0.5, \"#3b4cc0\"], [0, \"#5a6b7a\"],\n [0.5, \"#4caf50\"], [1.5, \"#ffeb3b\"], [3, \"#ff9800\"], [6, \"#e53935\"],\n ],\n },\n KDP: {\n label: \"Fase differenziale specifica\",\n unit: \"deg/km\",\n stops: [\n [-1, \"#6a3d9a\"], [-0.2, \"#3b4cc0\"], [0.1, \"#4a5560\"],\n [0.5, \"#4caf50\"], [1.5, \"#ffeb3b\"], [3, \"#ff9800\"], [6, \"#e53935\"],\n ],\n },\n PHIDP: {\n label: \"Fase differenziale\",\n unit: \"deg\",\n stops: [\n [0, \"#2b2b3d\"], [60, \"#3b4cc0\"], [120, \"#3ba99c\"],\n [180, \"#8fd14f\"], [240, \"#ffd166\"], [360, \"#ef476f\"],\n ],\n },\n};\n\nconst hexToRgb = (h: string): [number, number, number] => [\n parseInt(h.slice(1, 3), 16),\n parseInt(h.slice(3, 5), 16),\n parseInt(h.slice(5, 7), 16),\n];\n\n/** LUT a 512 livelli fra min e max, con interpolazione lineare fra gli stop. */\nexport function buildLut(stops: Stop[], min: number, max: number, n = 512): Uint8Array {\n const lut = new Uint8Array(n * 3);\n const rgb = stops.map(([, h]) => hexToRgb(h));\n for (let i = 0; i < n; i++) {\n const v = min + ((max - min) * i) / (n - 1);\n let k = 0;\n while (k < stops.length - 1 && v > stops[k + 1][0]) k++;\n let c: [number, number, number];\n if (v <= stops[0][0]) c = rgb[0];\n else if (v >= stops[stops.length - 1][0]) c = rgb[rgb.length - 1];\n else {\n const [v0] = stops[k], [v1] = stops[k + 1];\n const t = (v - v0) / (v1 - v0 || 1);\n c = [\n rgb[k][0] + t * (rgb[k + 1][0] - rgb[k][0]),\n rgb[k][1] + t * (rgb[k + 1][1] - rgb[k][1]),\n rgb[k][2] + t * (rgb[k + 1][2] - rgb[k][2]),\n ];\n }\n lut[i * 3] = c[0];\n lut[i * 3 + 1] = c[1];\n lut[i * 3 + 2] = c[2];\n }\n return lut;\n}\n\n/** Estremi della scala per una variabile (NI serve solo per la velocita'). */\nexport function scaleRange(varName: string, NI: number): [number, number] {\n const p = PALETTES[varName];\n if (p?.nyquist) return [-Math.abs(NI || 10), Math.abs(NI || 10)];\n const s = p?.stops ?? [[0, \"#000000\"] as Stop];\n return [s[0][0], s[s.length - 1][0]];\n}\n","export type Stop = [value: number, hex: string];\n\nexport interface Palette {\n label: string;\n unit: string;\n /** spiegazione in linguaggio comune, mostrata sotto la mappa */\n desc: string;\n stops: Stop[];\n /** valori sotto questa soglia = trasparenti (rumore/no eco) */\n floor?: number;\n /** se true, i limiti sono +/- NI letto dal file */\n nyquist?: boolean;\n}\n\nconst DBZ_STOPS: Stop[] = [\n [5, \"#04e9e7\"], [10, \"#019ff4\"], [15, \"#0300f4\"], [20, \"#02fd02\"],\n [25, \"#01c501\"], [30, \"#008e00\"], [35, \"#fdf802\"], [40, \"#e5bc00\"],\n [45, \"#fd9500\"], [50, \"#fd0000\"], [55, \"#d40000\"], [60, \"#bc0000\"],\n [65, \"#f800fd\"], [70, \"#9854c6\"],\n];\n\nexport const PALETTES: Record = {\n DBZH: {\n label: \"Riflettività\",\n unit: \"dBZ\",\n desc:\n \"Quanta acqua c'è nel fascio: è il radar “classico”, quello dei siti meteo. Più alto il valore, \" +\n \"più intensa la precipitazione. Sotto 20 dBZ è pioviggine, 35-45 è pioggia vera, oltre 50 è un \" +\n \"rovescio forte e sopra 55-60 quasi sempre c'è grandine. Attenzione: misura quello che il radar \" +\n \"vede in quota, non necessariamente quello che tocca terra.\",\n floor: 5,\n stops: DBZ_STOPS,\n },\n TH: {\n label: \"Riflettività totale (non filtrata)\",\n unit: \"dBZ\",\n desc:\n \"La stessa cosa della riflettività, ma senza i filtri che eliminano gli echi non meteorologici. \" +\n \"Qui vedi anche montagne, edifici, stormi e insetti. Serve per capire cosa il radar stia \" +\n \"realmente ricevendo, e per riconoscere i disturbi fissi che si ripetono uguali ogni giorno.\",\n floor: 5,\n stops: DBZ_STOPS,\n },\n RHOHV: {\n label: \"Coefficiente di correlazione\",\n unit: \"\",\n desc:\n \"Quanto si somigliano fra loro le particelle dentro il fascio. Vicino a 1 significa “tutte uguali”, \" +\n \"cioè pioggia ordinata di gocce simili. Sotto 0,9 c'è un miscuglio disordinato: grandine, neve che \" +\n \"si sta sciogliendo, oppure roba che non è affatto meteo — colline, palazzi, uccelli, insetti. \" +\n \"È il modo migliore per capire se un eco è vera precipitazione o spazzatura.\",\n stops: [\n [0.20, \"#2b2b3d\"], [0.65, \"#3b4cc0\"], [0.80, \"#4a9bd6\"],\n [0.90, \"#6bd0d5\"], [0.95, \"#b0f27a\"], [0.97, \"#ffd24d\"],\n [0.99, \"#ff6b4a\"], [1.05, \"#d21f3c\"],\n ],\n },\n VRADH: {\n label: \"Velocità radiale\",\n unit: \"m/s\",\n desc:\n \"Quanto velocemente l'aria (e la pioggia dentro di essa) si avvicina o si allontana dal radar. \" +\n \"Verde = viene verso di noi, rosso = se ne va, scuro = ferma. Importante: non è il vento completo, \" +\n \"ma solo la sua componente lungo la linea di vista. Un vento che soffia di traverso al raggio \" +\n \"risulta zero anche se è fortissimo. Le due tinte accostate bruscamente indicano aria che ruota \" +\n \"o converge, il segnale che i meteorologi cercano nei temporali.\",\n nyquist: true,\n stops: [\n [-1.0, \"#00ff88\"], [-0.60, \"#00c853\"], [-0.25, \"#1b7a3d\"],\n [-0.04, \"#123322\"], [0.04, \"#33161a\"], [0.25, \"#a33131\"],\n [0.60, \"#e53935\"], [1.0, \"#ff8a80\"],\n ],\n },\n WRADH: {\n label: \"Ampiezza spettrale\",\n unit: \"m/s\",\n desc:\n \"Quanto sono diverse fra loro le velocità nello stesso punto. Valori bassi = aria che si muove \" +\n \"compatta; valori alti = turbolenza, correnti che si mescolano, aria rimescolata dentro un temporale.\",\n stops: [\n [0, \"#101820\"], [1, \"#2c5f8a\"], [2, \"#3ba99c\"],\n [3, \"#8fd14f\"], [4, \"#ffd166\"], [6, \"#ef476f\"],\n ],\n },\n ZDR: {\n label: \"Riflettività differenziale\",\n unit: \"dB\",\n desc:\n \"Quanto sono schiacciate le particelle. Una goccia grande cade appiattita come un panino e dà \" +\n \"valori alti; una gocciolina è sferica e sta vicino a zero. Il chicco di grandine ruota su se \" +\n \"stesso mentre cade, quindi in media appare tondo: valori vicini a zero insieme a riflettività \" +\n \"altissima è la firma classica della grandine, ed è così che la si distingue dalla pioggia intensa.\",\n stops: [\n [-2, \"#6a3d9a\"], [-0.5, \"#3b4cc0\"], [0, \"#5a6b7a\"],\n [0.5, \"#4caf50\"], [1.5, \"#ffeb3b\"], [3, \"#ff9800\"], [6, \"#e53935\"],\n ],\n },\n KDP: {\n label: \"Fase differenziale specifica\",\n unit: \"°/km\",\n desc:\n \"Quanto il segnale viene “rallentato” attraversando la pioggia, chilometro per chilometro. \" +\n \"Cresce con la quantità di acqua liquida ed è quasi insensibile alla grandine, alla taratura \" +\n \"dell'antenna e all'attenuazione: quando la riflettività è dubbia, questo resta il miglior \" +\n \"indicatore di pioggia davvero forte in atto.\",\n stops: [\n [-1, \"#6a3d9a\"], [-0.2, \"#3b4cc0\"], [0.1, \"#4a5560\"],\n [0.5, \"#4caf50\"], [1.5, \"#ffeb3b\"], [3, \"#ff9800\"], [6, \"#e53935\"],\n ],\n },\n PHIDP: {\n label: \"Fase differenziale\",\n unit: \"°\",\n desc:\n \"Il dato grezzo da cui si ricava il KDP: lo sfasamento accumulato lungo tutto il percorso del \" +\n \"fascio. Cresce sempre man mano che ci si allontana dal radar, quindi preso da solo dice poco — \" +\n \"conta la sua pendenza, non il suo valore.\",\n stops: [\n [0, \"#2b2b3d\"], [60, \"#3b4cc0\"], [120, \"#3ba99c\"],\n [180, \"#8fd14f\"], [240, \"#ffd166\"], [360, \"#ef476f\"],\n ],\n },\n};\n\nconst hexToRgb = (h: string): [number, number, number] => [\n parseInt(h.slice(1, 3), 16),\n parseInt(h.slice(3, 5), 16),\n parseInt(h.slice(5, 7), 16),\n];\n\n/** LUT a 512 livelli fra min e max, con interpolazione lineare fra gli stop. */\nexport function buildLut(stops: Stop[], min: number, max: number, n = 512): Uint8Array {\n const lut = new Uint8Array(n * 3);\n const rgb = stops.map(([, h]) => hexToRgb(h));\n for (let i = 0; i < n; i++) {\n const v = min + ((max - min) * i) / (n - 1);\n let k = 0;\n while (k < stops.length - 1 && v > stops[k + 1][0]) k++;\n let c: [number, number, number];\n if (v <= stops[0][0]) c = rgb[0];\n else if (v >= stops[stops.length - 1][0]) c = rgb[rgb.length - 1];\n else {\n const [v0] = stops[k], [v1] = stops[k + 1];\n const t = (v - v0) / (v1 - v0 || 1);\n c = [\n rgb[k][0] + t * (rgb[k + 1][0] - rgb[k][0]),\n rgb[k][1] + t * (rgb[k + 1][1] - rgb[k][1]),\n rgb[k][2] + t * (rgb[k + 1][2] - rgb[k][2]),\n ];\n }\n lut[i * 3] = c[0];\n lut[i * 3 + 1] = c[1];\n lut[i * 3 + 2] = c[2];\n }\n return lut;\n}\n\n/** Estremi della scala per una variabile (NI serve solo per la velocita'). */\nexport function scaleRange(varName: string, NI: number): [number, number] {\n const p = PALETTES[varName];\n if (p?.nyquist) return [-Math.abs(NI || 10), Math.abs(NI || 10)];\n const s = p?.stops ?? [[0, \"#000000\"] as Stop];\n return [s[0][0], s[s.length - 1][0]];\n}\n","createdAt",["D",1784842683410],"branchVersion","branchId","1e41c00a-8f7b-4f3d-b14e-dcbcc1b563f8","id","c26973aa-d769-4132-bb0a-d51ae82f03d0","radarlom-polar",["D",1784813505760],"deletedAt","legacyEnvVars","privacy","public","httpPrivacy","description","Visualizzatore PPI polarimetrico radar ARPA Lombardia","imageUrl","019ba95c-929e-70ea-a632-3daf3d2c353b","pinnedFileId","d372be9e-869c-11f1-9568-1607ee4eb77e","pinPreview",true,"pinMigrated","screenshotVersion","likeCount",0,"orgId","019ba95c-929e-70ea-a632-40b422c292cb","tursoName","tursoHostname","tursoMeasuredAt","pinned","userPerm","r","valPerm","blobPerm","rw","emailPerm","sqlitePerm","projectPerm","telemetryPerm","user",{"_97":96,"_104":100},"org",{"_97":96,"_98":-5,"_99":100,"_18":55,"_101":102,"_103":25},"branch",{"_13":86,"_32":87,"_43":42,"_88":44,"_89":-5,"_90":-5,"_91":-5,"_92":-5,"_93":-5,"_94":45,"_95":96},"isOwned",52,"main","projectId","forkedBranchId","forkedBranchVersion","forkedProjectPrivacy","forkedProjectId","forkedProjectOrgId","projectName","userHandle","niko12345678","handle","displayName","avatar","https://img.clerk.com/eyJ0eXBlIjoicHJveHkiLCJzcmMiOiJodHRwczovL2ltYWdlcy5jbGVyay5kZXYvb2F1dGhfZ2l0aHViL2ltZ18zODR5N2ZXRGNYMUFxNEZyZTVSbW4zNGJHbFYifQ","tier","free","canRestrictVals","profileImageUrl","branchName","","adoptOrgId",{"_43":44,"_32":45,"_38":200,"_47":-5,"_48":25,"_49":50,"_51":50,"_52":53,"_54":-5,"_18":55,"_56":57,"_58":59,"_60":-5,"_61":-5,"_62":63,"_64":65,"_66":-5,"_67":-5,"_68":-5,"_69":-5,"_70":71,"_72":71,"_73":74,"_75":74,"_76":74,"_77":71,"_78":71,"_79":201,"_81":202,"_83":203,"_85":25},"tree",[179,180,181,182,183],"files",[179,180,181,182,183],"counts",{"_169":63,"_170":167,"_171":63,"_172":173,"_174":175},"contents",[127,128,129,130,131],"upstreamProject","upstreamBranch","emptyEnvVars","flags",[],"hasHttpFiles","grants","orgs","customDomains",{},{"_43":165,"_32":28,"_30":31,"_135":-5,"_41":42,"_40":14,"_137":59,"_38":166,"_13":167,"_33":34,"_35":37,"_24":25,"_141":-5,"_27":168},{"_43":157,"_32":158,"_30":57,"_135":-5,"_41":42,"_40":159,"_137":25,"_38":160,"_13":161,"_33":162,"_35":163,"_24":25,"_141":-5,"_27":164},{"_43":149,"_32":150,"_30":151,"_135":-5,"_41":42,"_40":152,"_137":59,"_38":153,"_13":154,"_33":34,"_35":155,"_24":25,"_141":-5,"_27":156},{"_43":143,"_32":144,"_30":145,"_135":-5,"_41":42,"_40":139,"_137":25,"_38":146,"_13":63,"_33":34,"_35":147,"_24":25,"_141":-5,"_27":148},{"_43":132,"_32":133,"_30":134,"_135":-5,"_41":42,"_40":136,"_137":25,"_38":138,"_13":139,"_33":34,"_35":140,"_24":25,"_141":-5,"_27":142},"86e15d84-dd22-4657-bd04-bfec564a5e0b","render.ts","b002aa6e-869c-11f1-ba2e-1607ee4eb77e","parentId",25,"latest",["D",1784842368720],5,"import { encode } from \"npm:fast-png@6.2.0\";\nimport type { Sweep } from \"./radar.ts\";\nimport { buildLut, PALETTES } from \"./colormaps.ts\";\n\nconst KM_PER_DEG_LAT = 110.574;\nconst kmPerDegLon = (lat: number) => 111.320 * Math.cos((lat * Math.PI) / 180);\nconst merc = (lat: number) => Math.log(Math.tan(Math.PI / 4 + (lat * Math.PI) / 360));\nconst invMerc = (y: number) => ((2 * Math.atan(Math.exp(y)) - Math.PI / 2) * 180) / Math.PI;\n\nexport interface Rendered {\n png: Uint8Array;\n bounds: [[number, number], [number, number]];\n width: number;\n height: number;\n}\n\n/**\n * Proietta uno o piu' volumi polari sulla stessa griglia Web Mercator.\n * Con piu' radar la combinazione e' pixel-wise max (composito di riflettivita'):\n * ha senso solo per DBZH/TH, dove il valore piu' alto e' quello meno attenuato\n * o meno schermato dall'orografia.\n */\nexport function renderPPI(sweeps: Sweep[], width = 1000, theme: \"dark\" | \"light\" = \"dark\"): Rendered {\n const first = sweeps[0];\n const p = PALETTES[first.varName] ?? PALETTES.DBZH;\n const nyq = !!p.nyquist;\n const lutMin = nyq ? -1 : p.stops[0][0];\n const lutMax = nyq ? 1 : p.stops[p.stops.length - 1][0];\n const lut = buildLut(p.stops, lutMin, lutMax);\n const nLut = lut.length / 3;\n const NI = Math.abs(first.NI) || 10;\n\n // bounding box: unione dei cerchi di portata\n let latS = 90, latN = -90, lonW = 180, lonE = -180;\n for (const sw of sweeps) {\n const dLat = sw.rmaxKm / KM_PER_DEG_LAT;\n const dLon = sw.rmaxKm / kmPerDegLon(sw.lat);\n latS = Math.min(latS, sw.lat - dLat);\n latN = Math.max(latN, sw.lat + dLat);\n lonW = Math.min(lonW, sw.lon - dLon);\n lonE = Math.max(lonE, sw.lon + dLon);\n }\n const yS = merc(latS), yN = merc(latN);\n const mercW = ((lonE - lonW) * Math.PI) / 180;\n const mercH = yN - yS;\n const height = Math.round((width * mercH) / mercW);\n\n const rgba = new Uint8Array(width * height * 4);\n const spanLon = lonE - lonW;\n const kmPerPx = (spanLon / width) * kmPerDegLon(first.lat);\n const ringTol = kmPerPx * 0.6;\n const rings = sweeps.length > 1 ? [50, 100, 125] : [25, 50, 75, 100, 125];\n // i ring devono contrastare con la basemap sottostante\n const ringRGB = theme === \"light\" ? [38, 46, 56] : [255, 255, 255];\n\n // stato per riga, uno per radar\n const rowDy = new Float64Array(sweeps.length);\n const rowKmLon = new Float64Array(sweeps.length);\n\n for (let y = 0; y < height; y++) {\n const lat = invMerc(yN - ((y + 0.5) / height) * mercH);\n for (let s = 0; s < sweeps.length; s++) {\n rowDy[s] = (lat - sweeps[s].lat) * KM_PER_DEG_LAT;\n rowKmLon[s] = kmPerDegLon((lat + sweeps[s].lat) / 2);\n }\n\n for (let x = 0; x < width; x++) {\n const lon = lonW + ((x + 0.5) / width) * spanLon;\n let best = -Infinity;\n let onRing = false;\n let inRange = false;\n\n for (let s = 0; s < sweeps.length; s++) {\n const sw = sweeps[s];\n const dxKm = (lon - sw.lon) * rowKmLon[s];\n const dyKm = rowDy[s];\n const r = Math.hypot(dxKm, dyKm);\n if (r > sw.rmaxKm) continue;\n inRange = true;\n\n for (let k = 0; k < rings.length; k++) {\n if (Math.abs(r - rings[k]) < ringTol) { onRing = true; break; }\n }\n\n const bin = Math.floor((r - sw.rstartKm) / sw.rscaleKm);\n if (bin < 0 || bin >= sw.nbins) continue;\n let az = Math.atan2(dxKm, dyKm) * (180 / Math.PI);\n if (az < 0) az += 360;\n const ray = Math.min(sw.nrays - 1, Math.floor((az * sw.nrays) / 360));\n const v = sw.raw[ray * sw.nbins + bin];\n if (v === sw.nodata || v === sw.undetect) continue;\n const val = sw.offset + sw.gain * v;\n if (val > best) best = val;\n }\n\n if (!inRange) continue;\n const o = (y * width + x) * 4;\n\n if (best > -Infinity && !(p.floor !== undefined && best < p.floor)) {\n const u = nyq ? Math.max(-1, Math.min(1, best / NI)) : best;\n let i = Math.round(((u - lutMin) / (lutMax - lutMin)) * (nLut - 1));\n i = Math.max(0, Math.min(nLut - 1, i));\n rgba[o] = lut[i * 3];\n rgba[o + 1] = lut[i * 3 + 1];\n rgba[o + 2] = lut[i * 3 + 2];\n rgba[o + 3] = 255;\n }\n\n if (onRing) {\n const a = rgba[o + 3] / 255;\n const w = 1 - a * 0.5;\n rgba[o] = Math.round(rgba[o] * a * 0.5 + ringRGB[0] * w);\n rgba[o + 1] = Math.round(rgba[o + 1] * a * 0.5 + ringRGB[1] * w);\n rgba[o + 2] = Math.round(rgba[o + 2] * a * 0.5 + ringRGB[2] * w);\n rgba[o + 3] = Math.max(rgba[o + 3], theme === \"light\" ? 110 : 90);\n }\n }\n }\n\n return {\n png: encode({ width, height, data: rgba, channels: 4, depth: 8 }),\n bounds: [[latS, lonW], [latN, lonE]],\n width,\n height,\n };\n}\n","label",[133],"4377e3c6-0572-4128-b0cb-e74780a7e433","radar.ts","919f531a-869c-11f1-bb07-1607ee4eb77e",["D",1784814244122],"import h5wasm from \"npm:h5wasm@0.7.8\";\n\nexport const BASE = \"https://radarlive.arpalombardia.it/Volumi\";\nexport const SITES: Record = { DES: \"Desio\", FLE: \"Flero\" };\nexport const VARS = [\"DBZH\", \"TH\", \"ZDR\", \"KDP\", \"PHIDP\", \"RHOHV\", \"VRADH\", \"WRADH\"];\n\nexport interface Sweep {\n site: string;\n varName: string;\n ts: string;\n elangle: number;\n elangles: number[];\n nrays: number;\n nbins: number;\n rscaleKm: number;\n rstartKm: number;\n rmaxKm: number;\n gain: number;\n offset: number;\n nodata: number;\n undetect: number;\n NI: number;\n lat: number;\n lon: number;\n raw: Uint16Array;\n}\n\nconst num = (v: unknown) => (typeof v === \"bigint\" ? Number(v) : Number(v));\n\nfunction attrs(f: any, path: string): Record {\n const g = f.get(path);\n if (!g) return {};\n const out: Record = {};\n for (const [k, v] of Object.entries(g.attrs as any)) out[k] = (v as any).value;\n return out;\n}\n\n/** Elenco dei timestamp disponibili (ultime 24h), dal piu' vecchio al piu' recente. */\nexport async function listTimes(site: string, varName: string): Promise {\n const res = await fetch(`${BASE}/${site}/`);\n if (!res.ok) throw new Error(`listing ${site}: HTTP ${res.status}`);\n const html = await res.text();\n const re = new RegExp(`${SITES[site]}\\\\.(\\\\d{8}T\\\\d{6}Z)_${varName}\\\\.h5\\\\.gz`, \"g\");\n const ts = [...new Set([...html.matchAll(re)].map((m) => m[1]))];\n ts.sort();\n return ts;\n}\n\n/** Scarica un volume, decomprime e ne estrae una singola elevazione. */\nexport async function readSweep(\n site: string,\n varName: string,\n ts: string,\n sweepIdx: number,\n): Promise {\n const url = `${BASE}/${site}/${SITES[site]}.${ts}_${varName}.h5.gz`;\n const res = await fetch(url);\n if (!res.ok) throw new Error(`volume non disponibile (HTTP ${res.status})`);\n const stream = res.body!.pipeThrough(new DecompressionStream(\"gzip\"));\n const bytes = new Uint8Array(await new Response(stream).arrayBuffer());\n\n const { FS } = await h5wasm.ready;\n const tmp = `/${crypto.randomUUID()}.h5`;\n FS.writeFile(tmp, bytes);\n const f = new h5wasm.File(tmp, \"r\");\n try {\n const groups = f.keys().filter((k: string) => /^dataset\\d+$/.test(k))\n .sort((a: string, b: string) => +a.slice(7) - +b.slice(7));\n const elangles = groups.map((g: string) => num(attrs(f, `${g}/where`).elangle));\n const idx = Math.min(Math.max(sweepIdx, 0), groups.length - 1);\n const g = groups[idx];\n\n const where = attrs(f, `${g}/where`);\n const how = attrs(f, `${g}/how`);\n const what = attrs(f, `${g}/data1/what`);\n const root = attrs(f, \"where\");\n\n const nbins = num(where.nbins);\n const rscaleKm = num(where.rscale) / 1000;\n const rstartKm = num(where.rstart);\n\n return {\n site,\n varName,\n ts,\n elangle: num(where.elangle),\n elangles,\n nrays: num(where.nrays),\n nbins,\n rscaleKm,\n rstartKm,\n rmaxKm: rstartKm + nbins * rscaleKm,\n gain: num(what.gain),\n offset: num(what.offset),\n nodata: num(what.nodata),\n undetect: num(what.undetect),\n NI: how.NI !== undefined ? num(how.NI) : 0,\n lat: num(root.lat),\n lon: num(root.lon),\n raw: (f.get(`${g}/data1/data`) as any).value as Uint16Array,\n };\n } finally {\n f.close();\n try {\n FS.unlink(tmp);\n } catch { /* ignore */ }\n }\n}\n",[144],"078e2108-49c7-4bf2-9982-76de4f3cecb3","probe.ts","d992698e-869a-11f1-906c-1607ee4eb77e",4,["D",1784813712462],3,"import h5wasm from \"npm:h5wasm@0.7.8\";\n\nconst BASE = \"https://radarlive.arpalombardia.it/Volumi\";\n\nconst J = (o: unknown) =>\n JSON.stringify(o, (_k, v) => {\n if (typeof v === \"bigint\") return Number(v);\n if (ArrayBuffer.isView(v)) return Array.from(v as any).slice(0, 8).map(Number);\n return v;\n });\n\nasync function latestFile(site: string, varName: string) {\n const html = await (await fetch(`${BASE}/${site}/`)).text();\n const re = new RegExp(`href=\"([^\"]*_${varName}\\\\.h5\\\\.gz)\"`, \"g\");\n const names = [...html.matchAll(re)].map((m) => m[1]);\n names.sort();\n return names[names.length - 1];\n}\n\nasync function fetchVolume(site: string, varName: string) {\n const name = await latestFile(site, varName);\n const res = await fetch(`${BASE}/${site}/${name}`);\n const ds = res.body!.pipeThrough(new DecompressionStream(\"gzip\"));\n return { name, bytes: new Uint8Array(await new Response(ds).arrayBuffer()) };\n}\n\nconst { name, bytes } = await fetchVolume(\"DES\", \"RHOHV\");\nconsole.log(\"file:\", name, bytes.length, \"bytes\");\n\nconst { FS } = await h5wasm.ready;\nFS.writeFile(\"/probe.h5\", bytes);\nconst f = new h5wasm.File(\"/probe.h5\", \"r\");\n\nconst attrs = (path: string) => {\n const g = f.get(path) as any;\n if (!g) return null;\n const out: Record = {};\n for (const [k, v] of Object.entries(g.attrs as any)) out[k] = (v as any).value;\n return out;\n};\n\nconsole.log(\"root/where:\", J(attrs(\"where\")));\nconsole.log(\"dataset1 keys:\", J((f.get(\"dataset1\") as any).keys()));\nconsole.log(\"dataset1/where:\", J(attrs(\"dataset1/where\")));\nconsole.log(\"dataset1/what:\", J(attrs(\"dataset1/what\")));\nconsole.log(\"dataset1/how:\", J(attrs(\"dataset1/how\")));\nconsole.log(\"data1 keys:\", J((f.get(\"dataset1/data1\") as any).keys()));\nconsole.log(\"data1/what:\", J(attrs(\"dataset1/data1/what\")));\n\nconst arr = f.get(\"dataset1/data1/data\") as any;\nconsole.log(\"shape:\", J(arr.shape), \"dtype:\", arr.dtype);\n\nfor (const s of f.keys().filter((k: string) => /^dataset\\d+$/.test(k))) {\n const w = attrs(`${s}/where`) as any;\n console.log(s, \"el\", Number(w.elangle), \"nrays\", Number(w.nrays), \"nbins\", Number(w.nbins), \"rscale\", Number(w.rscale));\n}\nf.close();\n",[150],"21dad4b4-48cf-4f57-80dc-91047ca6d70d","main.tsx",32,["D",1784842428731],19,"http","import { listTimes, readSweep, SITES } from \"./radar.ts\";\nimport { renderPPI } from \"./render.ts\";\nimport { PALETTES } from \"./colormaps.ts\";\n\nconst MOSAIC_SITES = [\"DES\", \"FLE\"];\nconst MOSAIC_VARS = [\"DBZH\", \"TH\"];\n\n/** Timestamp disponibili; per il mosaico, solo quelli presenti su entrambi i radar. */\nasync function times(site: string, v: string): Promise {\n if (site !== \"MOS\") return listTimes(site, v);\n const [a, b] = await Promise.all(MOSAIC_SITES.map((s) => listTimes(s, v)));\n const setB = new Set(b);\n return a.filter((t) => setB.has(t));\n}\n\nexport default async function (req: Request): Promise {\n const url = new URL(req.url);\n\n if (url.pathname === \"/api/times\") {\n const site = url.searchParams.get(\"site\") ?? \"FLE\";\n const v = url.searchParams.get(\"var\") ?? \"DBZH\";\n try {\n return Response.json({ times: await times(site, v) });\n } catch (e) {\n return Response.json({ error: String(e) }, { status: 502 });\n }\n }\n\n if (url.pathname === \"/ppi.png\") {\n const site = url.searchParams.get(\"site\") ?? \"FLE\";\n const v = url.searchParams.get(\"var\") ?? \"DBZH\";\n const elRaw = parseInt(url.searchParams.get(\"el\") ?? \"0\", 10);\n const el = Number.isFinite(elRaw) ? elRaw : 0;\n let t = url.searchParams.get(\"t\") ?? \"latest\";\n\n if (site === \"MOS\" && !MOSAIC_VARS.includes(v)) {\n return new Response(\n \"il mosaico e' definito solo per la riflettivita' (DBZH/TH): i momenti polarimetrici dipendono dalla geometria di vista e non si combinano fra radar\",\n { status: 400 },\n );\n }\n\n try {\n if (t === \"latest\") {\n const list = await times(site, v);\n if (!list.length) throw new Error(\"nessun volume disponibile\");\n t = list[list.length - 1];\n }\n\n const wanted = site === \"MOS\" ? MOSAIC_SITES : [site];\n const results = await Promise.all(\n wanted.map((s) => readSweep(s, v, t, el).then((sw) => sw).catch(() => null)),\n );\n const sweeps = results.filter((s): s is NonNullable => s !== null);\n if (!sweeps.length) throw new Error(\"volume non disponibile per questa scadenza\");\n\n const theme = url.searchParams.get(\"theme\") === \"light\" ? \"light\" : \"dark\";\n const out = renderPPI(sweeps, 1000, theme);\n const meta = {\n ts: sweeps[0].ts,\n elangle: sweeps[0].elangle,\n elangles: sweeps[0].elangles,\n NI: sweeps[0].NI,\n rmaxKm: sweeps[0].rmaxKm,\n nrays: sweeps[0].nrays,\n nbins: sweeps[0].nbins,\n bounds: out.bounds,\n centers: sweeps.map((s) => ({ site: s.site, name: SITES[s.site], lat: s.lat, lon: s.lon })),\n missing: wanted.filter((s) => !sweeps.some((sw) => sw.site === s)),\n };\n return new Response(out.png, {\n headers: {\n \"content-type\": \"image/png\",\n \"x-meta\": JSON.stringify(meta),\n \"access-control-expose-headers\": \"x-meta\",\n \"cache-control\": \"public, max-age=240\",\n },\n });\n } catch (e) {\n return new Response(String(e), { status: 502 });\n }\n }\n\n return new Response(page(), { headers: { \"content-type\": \"text/html; charset=utf-8\" } });\n}\n\nfunction page() {\n return `\n\n\nPPI polarimetrico - radar ARPA Lombardia\n\n\n

PPI polarimetrico — rete radar banda X, ARPA Lombardia

volumi ODIM HDF5 da radarlive.arpalombardia.it, decodificati e proiettati in Web Mercator

\n \n \n \n opacita' \n \n

\n \n \n \n \n \n —\n

caricamento...

\n\n`;\n}\n",[158],"668fd4a1-d2e4-4c32-949e-d17393921dee",["D",1784842683410],1,[28],"pulls","branches","forks","envVars",{"_176":63,"_178":63},"connections",{"_176":63,"_177":63},"count","unfilledCount","emptyCount",{"_43":165,"_32":28,"_30":31,"_135":-5,"_41":42,"_40":14,"_137":59,"_38":197,"_13":167,"_33":34,"_35":37,"_24":25,"_141":-5,"_185":198,"_27":199},{"_43":157,"_32":158,"_30":57,"_135":-5,"_41":42,"_40":159,"_137":25,"_38":194,"_13":161,"_33":162,"_35":163,"_24":25,"_141":-5,"_185":195,"_27":196},{"_43":149,"_32":150,"_30":151,"_135":-5,"_41":42,"_40":152,"_137":59,"_38":191,"_13":154,"_33":34,"_35":155,"_24":25,"_141":-5,"_185":192,"_27":193},{"_43":143,"_32":144,"_30":145,"_135":-5,"_41":42,"_40":139,"_137":25,"_38":188,"_13":63,"_33":34,"_35":147,"_24":25,"_141":-5,"_185":189,"_27":190},{"_43":132,"_32":133,"_30":134,"_135":-5,"_41":42,"_40":136,"_137":25,"_38":184,"_13":139,"_33":34,"_35":140,"_24":25,"_141":-5,"_185":186,"_27":187},["D",1784842368720],"children",[],[133],["D",1784814244122],[],[144],["D",1784813712462],[],[150],["D",1784842428731],[],[158],["D",1784842683410],[],[28],["D",1784813505760],{"_97":96,"_104":100},{"_97":96,"_98":-5,"_99":100,"_18":55,"_101":102,"_103":25},{"_13":86,"_32":87,"_43":42,"_88":44,"_89":-5,"_90":-5,"_91":-5,"_92":-5,"_93":-5,"_94":45,"_95":96},"ctx",{"_79":-5,"_124":237,"_238":239,"_240":241,"_242":25,"_243":25,"_244":25,"_245":59,"_246":-5,"_247":-7,"_248":-7,"_249":63,"_250":25,"_251":25,"_252":-5,"_253":25,"_254":25,"_255":25,"_256":-5},"ENV",{"_208":209,"_210":211,"_212":87,"_213":214,"_215":216,"_217":218,"_219":220,"_221":222,"_223":224,"_225":226,"_227":228,"_229":230,"_231":232,"_233":234,"_235":236},"REACT_APP_SENTRY_DSN","https://586020d86f05489db63a554690acf1e9@o1398359.ingest.sentry.io/6724499","RENDER_GIT_COMMIT","e87eb04e707cf64e7e532c8e8ec395f209e17f15","RENDER_GIT_BRANCH","POSTHOG_API_TOKEN","phc_nmoyuAkStEZuUp69rYlg1DZgrIt19Gx0Y4Hwf8tGgJP","POSTHOG_API_HOST","https://wart.val.town","WEB_API_TEMPLATE","https://{val}.web.val.run","APP_URL","https://www.val.town","MODULE_URL","https://esm.town","VANITY_WEB_API_TEMPLATE","https://{label}.val.run/","DENOLS_URL","https://lsp.val.town","NODE_ENV","production","CLOUDFLARE_TURNSTILE_SITE_KEY","0x4AAAAAACtaIZJH7s9rBMZ6","GOOGLE_CLIENT_ID","875989696899-jtsqkhlbdv5jrlsk8631p7ipjg3c6o62.apps.googleusercontent.com","STRIPE_PUBLISHABLE_KEY","pk_live_51N2zvWEpZTTLipIkVEwXrUG9AT0yeOOEYxszarZT1jpsZV58Ly6wlg3SltZpHN26CUalOBmUSY7adQRQiSkCZz2500qSXNkf8Y","API_URL","https://api.val.town",[],"preferencesKeybindings","default","preferencesTheme","light","preferencesFixers","preferencesCopilot","preferencesSelectionAi","preferencesCloseBrackets","preferencesChatPanelState","preferencesDprint","preferencesDenofmt","notificationCount","flagDenols","isBusinessOrgMember","businessOrgHandle","isBusinessOrgAdmin","isHigherTownieTier","navCollapsed","currentOrgId"]

Diff · v33