upload: example_usage.py
Browse files- example_usage.py +43 -0
example_usage.py
ADDED
|
@@ -0,0 +1,43 @@
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 1 |
+
"""
|
| 2 |
+
example_usage.py — Minimal inference with BeMAE-Hα.
|
| 3 |
+
|
| 4 |
+
This script loads a checkpoint from HuggingFace Hub, retrieves a spectrum
|
| 5 |
+
de test depuis le dataset BeSS-foundation, et calcule son embedding z_halpha.
|
| 6 |
+
"""
|
| 7 |
+
import json
|
| 8 |
+
import sys
|
| 9 |
+
|
| 10 |
+
import numpy as np
|
| 11 |
+
import torch
|
| 12 |
+
from huggingface_hub import snapshot_download
|
| 13 |
+
|
| 14 |
+
MODEL_ID = "anonym-submit-26/bemae-halpha-v1"
|
| 15 |
+
|
| 16 |
+
# 1. Download the snapshot
|
| 17 |
+
ckpt_dir = snapshot_download(MODEL_ID)
|
| 18 |
+
sys.path.insert(0, ckpt_dir)
|
| 19 |
+
|
| 20 |
+
from model import SpectralEncoderHalpha, ModelConfig # noqa: E402
|
| 21 |
+
|
| 22 |
+
# 2. Charger la config et les poids
|
| 23 |
+
with open(f"{ckpt_dir}/config.json") as f:
|
| 24 |
+
meta = json.load(f)
|
| 25 |
+
cfg = ModelConfig(**meta["model_config"])
|
| 26 |
+
|
| 27 |
+
encoder = SpectralEncoderHalpha(cfg)
|
| 28 |
+
state = torch.load(f"{ckpt_dir}/pytorch_model.bin", map_location="cpu")
|
| 29 |
+
encoder.load_state_dict(state)
|
| 30 |
+
encoder.eval()
|
| 31 |
+
|
| 32 |
+
# 3. Prepare a dummy spectrum (Hα-centred, 128 bins, pseudo-continuum at 1.0)
|
| 33 |
+
B = 4
|
| 34 |
+
flux = torch.ones(B, 128) + 0.3 * torch.randn(B, 128) * 0.01
|
| 35 |
+
wavelengths = torch.linspace(6512.8, 6612.8, 128).unsqueeze(0).expand(B, -1)
|
| 36 |
+
validity = torch.ones(B, 128)
|
| 37 |
+
|
| 38 |
+
# 4. Inference
|
| 39 |
+
with torch.no_grad():
|
| 40 |
+
z_halpha, *_ = encoder(flux, wavelengths, validity, mask=None)
|
| 41 |
+
|
| 42 |
+
print(f"z_halpha shape : {tuple(z_halpha.shape)}")
|
| 43 |
+
print(f"z_halpha[0, :8] : {z_halpha[0, :8].numpy()}")
|