363 lines
11 KiB
Python
363 lines
11 KiB
Python
import os
|
|
import json
|
|
import math
|
|
import numpy as np
|
|
import pandas as pd
|
|
import torch
|
|
import torch.nn as nn
|
|
import torch.nn.functional as F
|
|
import torch.optim as optim
|
|
|
|
from torch.utils.data import Dataset, DataLoader
|
|
from sklearn.model_selection import train_test_split
|
|
from sklearn.preprocessing import LabelEncoder, StandardScaler
|
|
from multiprocessing import Pool, cpu_count
|
|
from functools import partial
|
|
from tqdm import tqdm
|
|
|
|
|
|
def load_kaggle_asl_data(base_path):
|
|
train_df = pd.read_csv(os.path.join(base_path, "train.csv"))
|
|
with open(os.path.join(base_path, "sign_to_prediction_index_map.json")) as f:
|
|
sign_to_idx = json.load(f)
|
|
return train_df, sign_to_idx
|
|
|
|
|
|
def extract_hand_landmarks_from_parquet(path):
|
|
try:
|
|
df = pd.read_parquet(path)
|
|
# Take either left or right hand - prefer the one with more landmarks
|
|
left = df[df["type"] == "left_hand"]
|
|
right = df[df["type"] == "right_hand"]
|
|
|
|
hand = left if len(left) >= len(right) else right
|
|
|
|
if len(hand) == 0:
|
|
return None
|
|
|
|
frames = sorted(hand['frame'].unique())
|
|
landmarks_seq = []
|
|
|
|
for frame in frames:
|
|
lm_frame = hand[hand['frame'] == frame]
|
|
lm_list = []
|
|
for i in range(21):
|
|
lm = lm_frame[lm_frame['landmark_index'] == i]
|
|
if len(lm) == 0:
|
|
lm_list.append([0.0, 0.0, 0.0])
|
|
else:
|
|
lm_list.append([
|
|
float(lm['x'].iloc[0]),
|
|
float(lm['y'].iloc[0]),
|
|
float(lm['z'].iloc[0])
|
|
])
|
|
landmarks_seq.append(lm_list)
|
|
|
|
return np.array(landmarks_seq, dtype=np.float32) # (T, 21, 3)
|
|
except Exception:
|
|
return None
|
|
|
|
|
|
def get_features_sequence(landmarks_seq, max_frames=100):
|
|
if landmarks_seq is None or len(landmarks_seq) == 0:
|
|
return None
|
|
|
|
# Center on wrist
|
|
landmarks_seq = landmarks_seq - landmarks_seq[:, 0:1, :]
|
|
|
|
# Better scale: distance between index finger tip and middle finger tip
|
|
scale = np.linalg.norm(landmarks_seq[:, 8] - landmarks_seq[:, 12], axis=1, keepdims=True)
|
|
scale = np.maximum(scale, 1e-6)
|
|
landmarks_seq = landmarks_seq / scale
|
|
|
|
# Flatten to (T, 63)
|
|
seq = landmarks_seq.reshape(landmarks_seq.shape[0], -1)
|
|
|
|
# Pad or truncate
|
|
if len(seq) < max_frames:
|
|
pad = np.zeros((max_frames - len(seq), seq.shape[1]), dtype=np.float32)
|
|
seq = np.concatenate([seq, pad], axis=0)
|
|
else:
|
|
seq = seq[:max_frames]
|
|
|
|
return seq
|
|
|
|
|
|
def process_row(row, base_path, max_frames=100):
|
|
path = os.path.join(base_path, row['path'])
|
|
if not os.path.exists(path):
|
|
return None, None
|
|
|
|
try:
|
|
lm_seq = extract_hand_landmarks_from_parquet(path)
|
|
if lm_seq is None:
|
|
return None, None
|
|
|
|
feat_seq = get_features_sequence(lm_seq, max_frames)
|
|
if feat_seq is None:
|
|
return None, None
|
|
|
|
return feat_seq, row['sign']
|
|
except Exception:
|
|
return None, None
|
|
|
|
|
|
class PositionalEncoding(nn.Module):
|
|
def __init__(self, d_model, max_len=128):
|
|
super().__init__()
|
|
pe = torch.zeros(max_len, d_model)
|
|
position = torch.arange(0, max_len, dtype=torch.float).unsqueeze(1)
|
|
div_term = torch.exp(torch.arange(0, d_model, 2).float() * (-math.log(10000.0) / d_model))
|
|
pe[:, 0::2] = torch.sin(position * div_term)
|
|
pe[:, 1::2] = torch.cos(position * div_term)
|
|
self.register_buffer('pe', pe.unsqueeze(0))
|
|
|
|
def forward(self, x):
|
|
return x + self.pe[:, :x.size(1)]
|
|
|
|
|
|
class TransformerASL(nn.Module):
|
|
def __init__(self, input_dim=63, num_classes=250, d_model=192, nhead=6, num_layers=4):
|
|
super().__init__()
|
|
self.proj = nn.Linear(input_dim, d_model)
|
|
self.norm_in = nn.LayerNorm(d_model)
|
|
self.pos = PositionalEncoding(d_model)
|
|
|
|
encoder_layer = nn.TransformerEncoderLayer(
|
|
d_model=d_model,
|
|
nhead=nhead,
|
|
dim_feedforward=d_model * 4,
|
|
dropout=0.15,
|
|
activation='gelu',
|
|
batch_first=True,
|
|
norm_first=True
|
|
)
|
|
self.encoder = nn.TransformerEncoder(encoder_layer, num_layers=num_layers)
|
|
|
|
self.head = nn.Sequential(
|
|
nn.LayerNorm(d_model),
|
|
nn.Dropout(0.25),
|
|
nn.Linear(d_model, num_classes)
|
|
)
|
|
|
|
def forward(self, x, key_padding_mask=None):
|
|
x = self.proj(x)
|
|
x = self.norm_in(x)
|
|
x = self.pos(x)
|
|
x = self.encoder(x, src_key_padding_mask=key_padding_mask)
|
|
x = x.mean(dim=1) # global average pooling
|
|
return self.head(x)
|
|
|
|
|
|
def create_padding_mask(lengths, max_len):
|
|
return torch.arange(max_len, device=lengths.device)[None, :] >= lengths[:, None]
|
|
|
|
|
|
def main():
|
|
# ===============================
|
|
# DEVICE SETUP
|
|
# ===============================
|
|
device = torch.device("cuda" if torch.cuda.is_available() else "cpu")
|
|
print(f"Using device: {device}")
|
|
if device.type == "cuda":
|
|
print("GPU:", torch.cuda.get_device_name(0))
|
|
|
|
# ===============================
|
|
# PATHS & PARAMETERS
|
|
# ===============================
|
|
base_path = "asl_kaggle" # ← CHANGE THIS TO YOUR ACTUAL FOLDER
|
|
max_frames = 100
|
|
|
|
# ===============================
|
|
# DATA PROCESSING
|
|
# ===============================
|
|
print("Loading metadata...")
|
|
train_df, sign_to_idx = load_kaggle_asl_data(base_path)
|
|
|
|
print(f"Processing {len(train_df)} videos...")
|
|
rows = [row for _, row in train_df.iterrows()]
|
|
|
|
with Pool(cpu_count()) as pool:
|
|
results = list(tqdm(
|
|
pool.imap(
|
|
partial(process_row, base_path=base_path, max_frames=max_frames),
|
|
rows
|
|
),
|
|
total=len(rows),
|
|
desc="Processing"
|
|
))
|
|
|
|
X, y = [], []
|
|
for feat, sign in results:
|
|
if feat is not None:
|
|
X.append(feat)
|
|
y.append(sign)
|
|
|
|
if not X:
|
|
print("No valid sequences found!")
|
|
return
|
|
|
|
X = np.stack(X)
|
|
print(f"Loaded {len(X)} valid samples | shape: {X.shape}")
|
|
|
|
# Global normalization - very important!
|
|
print("Before global norm → mean:", X.mean(), "std:", X.std())
|
|
X = np.clip(X, -5.0, 5.0) # prevent crazy outliers
|
|
mean = X.mean(axis=(0, 1), keepdims=True)
|
|
std = X.std(axis=(0, 1), keepdims=True) + 1e-8
|
|
X = (X - mean) / std
|
|
print("After global norm → mean:", X.mean(), "std:", X.std())
|
|
|
|
# ===============================
|
|
# LABELS
|
|
# ===============================
|
|
le = LabelEncoder()
|
|
y = le.fit_transform(y)
|
|
num_classes = len(le.classes_)
|
|
print(f"Number of classes: {num_classes}")
|
|
|
|
# ===============================
|
|
# SPLIT
|
|
# ===============================
|
|
X_train, X_test, y_train, y_test = train_test_split(
|
|
X, y, test_size=0.15, stratify=y, random_state=42
|
|
)
|
|
|
|
# ===============================
|
|
# DATASET & LOADERS
|
|
# ===============================
|
|
class ASLSequenceDataset(Dataset):
|
|
def __init__(self, X, y):
|
|
self.X = torch.from_numpy(X).float()
|
|
self.y = torch.from_numpy(y).long()
|
|
|
|
def __len__(self):
|
|
return len(self.X)
|
|
|
|
def __getitem__(self, idx):
|
|
return self.X[idx], self.y[idx]
|
|
|
|
train_loader = DataLoader(
|
|
ASLSequenceDataset(X_train, y_train),
|
|
batch_size=64,
|
|
shuffle=True,
|
|
num_workers=4,
|
|
pin_memory=True
|
|
)
|
|
|
|
test_loader = DataLoader(
|
|
ASLSequenceDataset(X_test, y_test),
|
|
batch_size=96,
|
|
shuffle=False,
|
|
num_workers=4,
|
|
pin_memory=True
|
|
)
|
|
|
|
# ===============================
|
|
# MODEL
|
|
# ===============================
|
|
model = TransformerASL(
|
|
input_dim=63,
|
|
num_classes=num_classes,
|
|
d_model=192,
|
|
nhead=6,
|
|
num_layers=4
|
|
).to(device)
|
|
|
|
print(f"Model parameters: {sum(p.numel() for p in model.parameters()):,}")
|
|
|
|
# ===============================
|
|
# TRAINING SETUP
|
|
# ===============================
|
|
criterion = nn.CrossEntropyLoss(label_smoothing=0.05)
|
|
optimizer = optim.AdamW(model.parameters(), lr=5e-4, weight_decay=1e-4)
|
|
scheduler = optim.lr_scheduler.CosineAnnealingWarmRestarts(optimizer, T_0=10)
|
|
|
|
# ===============================
|
|
# TRAIN / EVAL FUNCTIONS
|
|
# ===============================
|
|
def train_epoch():
|
|
model.train()
|
|
total_loss = 0
|
|
correct = 0
|
|
total = 0
|
|
|
|
for x, y in tqdm(train_loader, desc="Training"):
|
|
x, y = x.to(device), y.to(device)
|
|
|
|
# Rough length estimation
|
|
lengths = (x.abs().sum(dim=2) > 1e-5).sum(dim=1)
|
|
mask = create_padding_mask(lengths, x.size(1))
|
|
|
|
optimizer.zero_grad(set_to_none=True)
|
|
logits = model(x, key_padding_mask=mask)
|
|
|
|
loss = criterion(logits, y)
|
|
loss.backward()
|
|
|
|
grad_norm = torch.nn.utils.clip_grad_norm_(model.parameters(), max_norm=0.8)
|
|
|
|
if torch.isnan(loss) or grad_norm > 20:
|
|
print(f"Warning - large grad or NaN! norm = {grad_norm:.2f}")
|
|
|
|
optimizer.step()
|
|
|
|
total_loss += loss.item()
|
|
correct += (logits.argmax(dim=-1) == y).sum().item()
|
|
total += y.size(0)
|
|
|
|
return total_loss / len(train_loader), correct / total * 100
|
|
|
|
@torch.no_grad()
|
|
def evaluate():
|
|
model.eval()
|
|
correct = 0
|
|
total = 0
|
|
for x, y in test_loader:
|
|
x, y = x.to(device), y.to(device)
|
|
lengths = (x.abs().sum(dim=2) > 1e-5).sum(dim=1)
|
|
mask = create_padding_mask(lengths, x.size(1))
|
|
|
|
logits = model(x, key_padding_mask=mask)
|
|
correct += (logits.argmax(dim=-1) == y).sum().item()
|
|
total += y.size(0)
|
|
return correct / total * 100 if total > 0 else 0
|
|
|
|
# ===============================
|
|
# TRAINING LOOP
|
|
# ===============================
|
|
best_acc = 0
|
|
patience = 15
|
|
wait = 0
|
|
epochs = 60
|
|
|
|
for epoch in range(epochs):
|
|
loss, train_acc = train_epoch()
|
|
test_acc = evaluate()
|
|
|
|
print(f"Epoch {epoch + 1:2d}/{epochs} | Loss: {loss:.4f} | Train: {train_acc:.2f}% | Test: {test_acc:.2f}%")
|
|
|
|
scheduler.step()
|
|
|
|
if test_acc > best_acc:
|
|
best_acc = test_acc
|
|
wait = 0
|
|
torch.save({
|
|
'model': model.state_dict(),
|
|
'optimizer': optimizer.state_dict(),
|
|
'label_encoder': le.classes_,
|
|
'epoch': epoch,
|
|
'acc': best_acc
|
|
}, "best_asl_transformer.pth")
|
|
print(" → New best model saved")
|
|
else:
|
|
wait += 1
|
|
if wait >= patience:
|
|
print("Early stopping triggered")
|
|
break
|
|
|
|
print(f"\nTraining finished. Best test accuracy: {best_acc:.2f}%")
|
|
|
|
|
|
if __name__ == '__main__':
|
|
main() |