Hareket eden gökyüzünü 73 yıllık arşivde yakala

Her sorguda IRSA'nın SIA servisinden o noktayı kapsayan bütün SPHEREx görüntüleri alınır; haftalık yeni veriler kendiliğinden gelir. Yanıt ~330 KB'lık bir CSV. Eskiden karakter karakter ayrıştırılıyordu ve Cloudflare'in işlem süresi sınırını aşıyordu (hata 1102); şimdi satır satır bölünüyor: 25,7 ms → 2,4 ms, sonuç birebir aynı.

async function search(u) {
  const ra = num(u.searchParams.get("ra"), 0, 360), dec = num(u.searchParams.get("dec"), -90, 90);
  const q = new URLSearchParams({ FACILITY: "SPHEREx", POS: `circle ${ra} ${dec} 0.002`, RESPONSEFORMAT: "CSV" });
  const r = await fetch(`${IRSA}/SIA?${q}`);
  if (!r.ok) throw new Error("IRSA SIA " + r.status);
  // Hızlı ayrıştırma: gereken sütunların hepsi, tırnaklı JSON içeren cloud_access sütunundan önce ve virgülsüz.
  // (Karakter karakter ayrıştırma ~330 KB'lık yanıtta Worker'ın CPU sınırını aşıp 1102 hatası veriyordu.)
  const rows = (await r.text()).split("\n").filter(l => l.length > 1).map(l => l.split(","));
  const h = rows.shift() || [], c = k => h.indexOf(k);
  const out = rows.filter(r => r[c("dataproduct_subtype")] === "science" && /\/level2\//.test(r[c("access_url")] || "")).map(r => ({
    url: r[c("access_url")],
    det: +r[c("energy_bandpassname")].replace(/\D/g, ""),
    wmin: +(+r[c("em_min")] * 1e6).toFixed(3), wmax: +(+r[c("em_max")] * 1e6).toFixed(3),
    mjd: (+r[c("t_min")] + +r[c("t_max")]) / 2,
    obs: r[c("obs_id")], coll: r[c("obs_collection")],
    ...detPos(r[c("s_region")], ra, dec),
  })).sort((a, b) => a.mjd - b.mjd);
  return json({ ra, dec, n: out.length, guncel: new Date().toISOString(), goruntuler: out });
}

Bir SPHEREx görüntüsü ~70 MB. Worker, NASA'nın herkese açık S3 deposundan HTTP Range ile yalnız başlığı ve hedefin çevresindeki satırları okur; sıkıştırılmış bayrak katmanını (RICE) kendisi çözer. Olmazsa IRSA'nın kesit servisine döner.

// SPHEREx kesitini IRSA'nın kesit servisi yerine doğrudan herkese açık S3 deposundan (nasa-irsa-spherex) üretir.
// Yalnız gereken satırlar HTTP Range ile okunur; sıkıştırılmış FLAGS katmanı (RICE_1) burada çözülür.
// Çıktı, IRSA kesitiyle aynı biçimde: IMAGE, FLAGS, VARIANCE, ZODI katmanları, kesite göre kaydırılmış WCS ile.
// Herhangi bir sorunda hata fırlatır; çağıran taraf IRSA'ya geri döner.

async function aralik(url, a, b) {
  const r = await fetch(url, { headers: { Range: `bytes=${a}-${b}` } });
  if (r.status !== 206 && r.status !== 200) throw new Error("S3 " + r.status);
  return new Uint8Array(await r.arrayBuffer());
}

Her yıldızın her SPHEREx karesindeki ağırlık merkezi bulunur ve zamana göre ağırlıklı bir doğru uydurulur. Pencere yıldızın parlaklık yarıçapına göre büyür: doymuş parlak bir yıldızın çekirdeği düz olduğundan küçük pencere hareketi göremez.

function yildizIzle(G,Z,q,sc,T){ // sc: yay saniyesi / piksel; dönüş piksel/yıl
  const psf=Math.max(2,6.2/sc), Rw=Math.min(18,Math.max(Math.ceil(psf)+2,Math.ceil(q.r+psf+2))), P=[];
  G.forEach((f,k)=>{ const c=izMerkez(Z[k],q.x,q.y,Rw,T); if(c&&Math.hypot(c[0]-q.x,c[1]-q.y)<Rw) P.push({t:f.s.epoch,x:c[0],y:c[1],w:Math.min(c[2],50)}); }); // ağırlık: tepe değeri (S/G), doymuşta 50
  if(P.length<Math.max(3,Math.ceil(0.7*G.length))) return null;
  const W=P.reduce((a,o)=>a+o.w,0), tm=P.reduce((a,o)=>a+o.w*o.t,0)/W, Stt=P.reduce((a,o)=>a+o.w*(o.t-tm)**2,0); if(Stt<=0) return null;
  const egim=k=>{ const m=P.reduce((a,o)=>a+o.w*o[k],0)/W; return [P.reduce((a,o)=>a+o.w*(o.t-tm)*(o[k]-m),0)/Stt, m]; };
  const [vx,ox]=egim("x"), [vy,oy]=egim("y");
  const chi=P.reduce((a,o)=>a+o.w*((o.x-ox-vx*(o.t-tm))**2+(o.y-oy-vy*(o.t-tm))**2),0)/Math.max(1,2*(P.length-2));
  return {vx,vy,ox,oy,tm, hata:Math.sqrt(Math.max(chi,0.02*W/P.length)/Stt), vv:Math.hypot(vx,vy), nk:P.length, Rw};
}

Girişteki tüm gökyüzü haritası Galaktik koordinatlarda 2:1 Mollweide. Seçilen yıldızın yeri aynı formülle altın halkayla işaretlenir; 2θ + sin 2θ = π sin b denklemi Newton yöntemiyle çözülür.

function mollXY(l,b){ // l,b radyan → birim elips içinde (x: -1..1, y: -1..1)
  l=((l+Math.PI)%(2*Math.PI)+2*Math.PI)%(2*Math.PI)-Math.PI;
  let t=b; if(Math.abs(Math.abs(b)-Math.PI/2)<1e-9) t=b; else for(let i=0;i<30;i++){ const d=(2*t+Math.sin(2*t)-Math.PI*Math.sin(b))/(2+2*Math.cos(2*t)); t-=d; if(Math.abs(d)<1e-9) break; }
  return [-l*Math.cos(t)/Math.PI, Math.sin(t)]; }
function