运行GradCAM遇报错：Sequential层未被调用无输出，求助解决

运行GradCAM时触发如下报错：

ValueError: The layer sequential_22 has never been called and thus has no defined output

代码中已明确训练并调用过Sequential模型，尝试过卸载TensorFlow、重启内核、更换运行环境均无法解决，此前未出现该问题。

报错核心是GradCAM构建的子模型无法识别原模型的输出定义，具体修复点如下：

• 删除重复的TensorFlow导入：代码在模型训练完成后重复执行import tensorflow as tf，会重置模型上下文，导致已训练的模型状态丢失。直接移除该行即可。
• 修正GradCAM的维度计算：原代码针对2D卷积设计维度平均，你的模型使用1D卷积，需将pooled_grads = tf.reduce_mean(grads, axis=(0, 1, 2))改为pooled_grads = tf.reduce_mean(grads, axis=(0,1))，避免无效维度的平均操作。
• 确保模型前向传播状态有效：删除重复导入后，原模型通过fit和evaluate建立的前向传播记录会被正确保留，GradCAM子模型能正常获取层输出。

# Build the Model
import tensorflow as tf
import numpy as np
import matplotlib.pyplot as plt
import random

model = tf.keras.Sequential()
model.add(Conv1D(filters=64, kernel_size=3, activation='relu', input_shape=(signal_length, 1), name = 'first_conv'))
model.add(MaxPooling1D(pool_size=2))
model.add(Dropout(0.2))

model.add(Conv1D(filters = 128, kernel_size = 5, activation = 'relu', name = 'second_conv'))
model.add(BatchNormalization())
model.add(MaxPooling1D(pool_size = 2))
model.add(Dropout(0.2))

model.add(Conv1D(filters = 256, kernel_size = 10, activation = 'relu', name='last_conv'))
model.add(BatchNormalization())
model.add(MaxPooling1D(pool_size = 2))
model.add(Dropout(0.2))

model.add(Flatten())

model.add(Dense(units=128, activation='relu'))
model.add(Dropout(0.2))

model.add(Dense(units=5, activation='softmax')) # 5 classes for N, S, V, F, Q

# Compile the Model
model.compile(loss='categorical_crossentropy', optimizer='adam', metrics=['accuracy'])

# Print Model.output
model.summary()

# Train the Model
history = model.fit(X_train, y_train, epochs=5, batch_size=32, validation_data=(X_test, y_test))
loss, accuracy = model.evaluate(X_test, y_test)
print("Loss: ", loss)
print("Accuracy: ", accuracy)

plt.plot(history.history['accuracy'], label='train accuracy')
plt.plot(history.history['val_accuracy'], label='val accuracy')
plt.xlabel('Epoch')
plt.ylabel('Accuracy')
plt.legend()
plt.show()

def grad_cam(model, layer_name, data, pred_index):
    # Define a model that outputs the layer's activations and the final prediction
    grad_model = tf.keras.models.Model(
        inputs=[model.inputs], 
        outputs=[model.get_layer(layer_name).output, model.output]
    )

    with tf.GradientTape() as tape:
        conv_outputs, predictions = grad_model(data)
        loss = predictions[:, pred_index]

    grads = tape.gradient(loss, conv_outputs)
    # 针对1D卷积调整维度平均轴
    pooled_grads = tf.reduce_mean(grads, axis=(0, 1))  

    conv_outputs = conv_outputs[0]
    heatmap = tf.reduce_mean(tf.multiply(pooled_grads, conv_outputs), axis=-1)
    heatmap = np.maximum(heatmap, 0) / (np.max(heatmap) + tf.keras.backend.epsilon())  # Avoid division by zero
    return heatmap

def plot_gradcam(ecg_data, cam, prediction, class_names, title='GradCAM'):
    fig, ax = plt.subplots(figsize=(10, 5))
    ax.plot(ecg_data.flatten(), alpha=0.8, label='ECG Signal')

    cax = ax.imshow(np.expand_dims(cam, axis=0), cmap='jet', aspect='auto', alpha=0.5, extent=[0, len(ecg_data.flatten()), ax.get_ylim()[0], ax.get_ylim()[1]])
    fig.colorbar(cax, ax=ax, orientation='horizontal', label='Activation Intensity')
    ax.set_title(title)
    ax.set_xlabel('Time (ms)')
    ax.set_ylabel('Activation')
    ax.legend()

    text = f'Classification: {class_names[np.argmax(prediction)]}\nPrediction: {class_names[np.argmax(prediction)]}\nConfidence: {np.max(prediction):.2f}'
    print(text)
    plt.show()

# Example usage
sample_index = random.randint(0, X_test.shape[0] - 1)
sample_data = X_test[sample_index:sample_index + 1]
class_index = np.argmax(y_test[sample_index])

cam = grad_cam(model, 'last_conv', sample_data, class_index)
prediction = model.predict(sample_data)
class_names = list(id_to_label.values())

for i in range(10):
    sample_index = random.randint(0, X_test.shape[0] - 1)
    sample_data = X_test[sample_index:sample_index + 1]
    class_index = np.argmax(y_test[sample_index])
    cam = grad_cam(model, 'last_conv', sample_data, class_index)
    prediction = model.predict(sample_data)
    plot_gradcam(sample_data[0], cam, prediction, class_names, title=f'GradCAM for Sample {i + 1}')

内容的提问来源于stack exchange，提问作者Noah Gallego