1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151
| import torch import torch.nn as nn
class EyeTAG(nn.Module): """ EyeTAG: 眼动轨迹感知视线估计 核心设计: 1. 单帧编码器提取每帧眼部特征 2. 轨迹解码器建模眼动时序 3. 眼跳检测门控区分注视/眼跳 输出: 平滑且时序一致的视线轨迹 """ def __init__( self, n_channels: int = 3, feat_dim: int = 256, n_heads: int = 4, n_layers: int = 2 ): super().__init__() self.frame_encoder = nn.Sequential( nn.Conv2d(n_channels, 32, 3, padding=1), nn.ReLU(), nn.MaxPool2d(2), nn.Conv2d(32, 64, 3, padding=1), nn.ReLU(), nn.MaxPool2d(2), nn.Conv2d(64, 128, 3, padding=1), nn.ReLU(), nn.AdaptiveAvgPool2d((1, 1)), nn.Flatten(), nn.Linear(128, feat_dim) ) self.temporal_transformer = nn.TransformerEncoder( nn.TransformerEncoderLayer( d_model=feat_dim, nhead=n_heads, dim_feedforward=feat_dim * 2, dropout=0.1, batch_first=True ), num_layers=n_layers ) self.gaze_head = nn.Sequential( nn.Linear(feat_dim, 64), nn.ReLU(), nn.Linear(64, 2) ) self.saccade_head = nn.Sequential( nn.Linear(feat_dim, 32), nn.ReLU(), nn.Linear(32, 1), nn.Sigmoid() ) def forward(self, eye_sequence: torch.Tensor) -> dict: """ Args: eye_sequence: (B, T, C, H, W) 连续T帧眼部图像 Returns: gaze: (B, T, 2) 平滑视线轨迹 saccade_prob: (B, T) 眼跳概率 """ B, T, C, H, W = eye_sequence.shape frames_flat = eye_sequence.view(B * T, C, H, W) feats = self.frame_encoder(frames_flat) feats = feats.view(B, T, -1) temporal_out = self.temporal_transformer(feats) gaze = self.gaze_head(temporal_out) saccade_prob = self.saccade_head(temporal_out).squeeze(-1) return { 'gaze': gaze, 'saccade_prob': saccade_prob, 'features': temporal_out }
def trajectory_consistency_loss( pred_gaze: torch.Tensor, target_gaze: torch.Tensor, saccade_prob: torch.Tensor, alpha: float = 0.5, beta: float = 0.1 ) -> torch.Tensor: """ 轨迹一致性损失 三部分: 1. L1回归损失 2. 时序平滑损失(注视段惩罚变化) 3. 眼跳检测损失 注视段(saccade_prob < 0.5)施加强平滑约束 眼跳段(saccade_prob > 0.5)放松平滑约束 """ l1 = (pred_gaze - target_gaze).abs().mean() diff = pred_gaze[:, 1:] - pred_gaze[:, :-1] smooth_mask = (saccade_prob[:, :-1] < 0.5).float().unsqueeze(-1) smooth_loss = (diff.abs() * smooth_mask).mean() total = l1 + alpha * smooth_loss return total
if __name__ == "__main__": model = EyeTAG(n_channels=3, feat_dim=128) eye_seq = torch.randn(4, 30, 3, 64, 32) out = model(eye_seq) print(f"输入: {eye_seq.shape}") print(f"视线轨迹: {out['gaze'].shape}") print(f"眼跳概率: {out['saccade_prob'].shape}") gaze_diff = out['gaze'][:, 1:] - out['gaze'][:, :-1] jitter = gaze_diff.abs().mean().item() print(f"帧间抖动: {jitter:.4f}° (EyeTAG)") single_frame_gaze = model.gaze_head( model.frame_encoder(eye_seq.view(4*30, 3, 64, 32)) ).view(4, 30, 2) sf_diff = single_frame_gaze[:, 1:] - single_frame_gaze[:, :-1] sf_jitter = sf_diff.abs().mean().item() print(f"帧间抖动: {sf_jitter:.4f}° (单帧)") print(f"抖动降低: {(1 - jitter/sf_jitter)*100:.1f}%")
|