def init_extractors(self, gam_type):
if gam_type == 'scoregam':
self.cam_extractors = [ScoreCAM(self.cam_model, self.conv_layer, self.input_layer)]
elif gam_type == 'gradcampp':
self.cam_extractors = [GradCAMpp(self.cam_model, self.conv_layer)]
elif gam_type == 'smoothgradcampp':
self.cam_extractors = [SmoothGradCAMpp(self.cam_model, self.conv_layer, self.input_layer)]
def main(args):
if args.device is None:
args.device = 'cuda:0' if torch.cuda.is_available() else 'cpu'
device = torch.device(args.device)
# Pretrained imagenet model
model = models.__dict__[args.model](pretrained=True).eval().to(device=device)
conv_layer = MODEL_CONFIG[args.model]['conv_layer']
input_layer = MODEL_CONFIG[args.model]['input_layer']
fc_layer = MODEL_CONFIG[args.model]['fc_layer']
# Image
if args.img.startswith('http'):
img_path = BytesIO(requests.get(args.img).content)
pil_img = Image.open(img_path, mode='r').convert('RGB')
# Preprocess image
img_tensor = normalize(to_tensor(resize(pil_img, (224, 224))),
[0.485, 0.456, 0.406], [0.229, 0.224, 0.225]).to(device=device)
# Hook the corresponding layer in the model
cam_extractors = [CAM(model, conv_layer, fc_layer), GradCAM(model, conv_layer),
GradCAMpp(model, conv_layer), SmoothGradCAMpp(model, conv_layer, input_layer),
ScoreCAM(model, conv_layer, input_layer), SSCAM(model, conv_layer, input_layer),
ISSCAM(model, conv_layer, input_layer)]
# Don't trigger all hooks
for extractor in cam_extractors:
extractor._hooks_enabled = False
fig, axes = plt.subplots(1, len(cam_extractors), figsize=(7, 2))
for idx, extractor in enumerate(cam_extractors):
extractor._hooks_enabled = True
model.zero_grad()
scores = model(img_tensor.unsqueeze(0))
# Select the class index
class_idx = scores.squeeze(0).argmax().item() if args.class_idx is None else args.class_idx
# Use the hooked data to compute activation map
activation_map = extractor(class_idx, scores).cpu()
# Clean data
extractor.clear_hooks()
extractor._hooks_enabled = False
# Convert it to PIL image
# The indexing below means first image in batch
heatmap = to_pil_image(activation_map, mode='F')
# Plot the result
result = overlay_mask(pil_img, heatmap)
axes[idx].imshow(result)
axes[idx].axis('off')
axes[idx].set_title(extractor.__class__.__name__, size=8)
plt.tight_layout()
if args.savefig:
plt.savefig(args.savefig, dpi=200, transparent=True, bbox_inches='tight', pad_inches=0)
plt.show()
def visualize_cam_on_img(img_name, model):
cam_extractor = SmoothGradCAMpp(model)
img = read_image(str(img_name))
input_tensor = normalize(
resize(img, (224, 224)) / 255., [0.485, 0.456, 0.406],
[0.229, 0.224, 0.225]).cuda()
out = model(input_tensor.unsqueeze(0))
activation_map = cam_extractor(out.squeeze(0).argmax().item(), out)
result = overlay_mask(to_pil_image(img),
to_pil_image(activation_map, mode='F'),
alpha=0.5)
plt.imshow(result)
plt.axis('off')
plt.tight_layout()
plt.show()
def for_vis(args):
transform = transforms.Compose([
transforms.Resize((args.img_size, args.img_size)),
transforms.ToTensor(),
])
# Con-text
if args.dataset == 'ConText':
train, val = MakeList(args).get_data()
dataset_val = ConText(val, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
data = iter(data_loader_val).next()
image = data["image"][0]
label = data["label"][0]
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
elif args.dataset == 'ImageNet':
train, val = MakeListImage(args).get_data()
dataset_val = ConText(val, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
iter_loader = iter(data_loader_val)
for i in range(0, 1):
data = iter_loader.next()
image = data["image"][0]
label = data["label"][0].item()
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
# MNIST
elif args.dataset == 'MNIST':
dataset_val = datasets.MNIST('./data/mnist',
train=False,
transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
image = iter(data_loader_val).next()[0][0]
label = ''
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8)[0], mode='L')
image = transform(image_orl)
transform = transforms.Compose(
[transforms.Normalize((0.1307, ), (0.3081, ))])
# CUB
elif args.dataset == 'CUB200':
dataset_val = CUB_200(args, train=False, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
data = iter(data_loader_val).next()
image = data["image"][0]
label = data["label"][0]
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
image = transform(image)
image = image.unsqueeze(0)
device = torch.device(args.device)
model = load_backbone(args)
model = model.to(device)
RESNET_CONFIG = dict(input_layer='conv1',
conv_layer='layer4',
fc_layer='fc')
MODEL_CONFIG = {**RESNET_CONFIG}
conv_layer = MODEL_CONFIG['conv_layer']
input_layer = MODEL_CONFIG['input_layer']
fc_layer = MODEL_CONFIG['fc_layer']
model.eval()
# Hook the corresponding layer in the model
cam_extractors = [
CAM(model, conv_layer, fc_layer),
GradCAM(model, conv_layer),
GradCAMpp(model, conv_layer),
SmoothGradCAMpp(model, conv_layer, input_layer),
ScoreCAM(model, conv_layer, input_layer),
SSCAM(model, conv_layer, input_layer),
# ISSCAM(model, conv_layer, input_layer),
]
cam_extractors_names = [
'CAM',
'GradCAM',
'GradCAMpp',
'SmoothGradCAMpp',
'ScoreCAM',
'SSCAM',
# 'ISSCAM',
]
for idx, extractor in enumerate(cam_extractors):
model.zero_grad()
output1 = model(image.to(device))
output = F.softmax(output1, dim=1)
prediction_score, pred_label_idx = torch.topk(output, 1)
pred_label_idx.squeeze_()
predicted_label = str(pred_label_idx.item() + 1)
print('Predicted:', predicted_label, '(',
prediction_score.squeeze().item(), ')')
make_grad(extractor, output1, image_orl, args.grad_min_level,
cam_extractors_names[idx])
extractor.clear_hooks()
def for_vis(args):
transform = transforms.Compose([
transforms.Resize((args.img_size, args.img_size)),
transforms.ToTensor(),
])
# Con-text
if args.dataset == 'ConText':
train, val = MakeList(args).get_data()
dataset_val = ConText(val, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
data = iter(data_loader_val).next()
image = data["image"][0]
label = data["label"][0]
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
elif args.dataset == 'ImageNet':
train, val = MakeListImage(args).get_data()
dataset_val = ConText(val, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
iter_loader = iter(data_loader_val)
for i in range(0, 1):
data = iter_loader.next()
image = data["image"][0]
label = data["label"][0].item()
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
# MNIST
elif args.dataset == 'MNIST':
dataset_val = datasets.MNIST('./data/mnist',
train=False,
transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
image = iter(data_loader_val).next()[0][0]
label = ''
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8)[0], mode='L')
image = transform(image_orl)
transform = transforms.Compose(
[transforms.Normalize((0.1307, ), (0.3081, ))])
# CUB
elif args.dataset == 'CUB200':
dataset_val = CUB_200(args, train=False, transform=transform)
data_loader_val = torch.utils.data.DataLoader(dataset_val,
args.batch_size,
shuffle=False,
num_workers=1,
pin_memory=True)
data = iter(data_loader_val).next()
image = data["image"][0]
label = data["label"][0]
image_orl = Image.fromarray(
(image.cpu().detach().numpy() * 255).astype(np.uint8).transpose(
(1, 2, 0)),
mode='RGB')
image = transform(image_orl)
transform = transforms.Compose([
transforms.Normalize([0.485, 0.456, 0.406], [0.229, 0.224, 0.225])
])
image = transform(image)
image = image.unsqueeze(0)
device = torch.device(args.device)
### IGOS
model = load_backbone(args)
model = model.to(device)
model.eval()
image_orl_for_blur = np.float32(image_orl) / 255.
img, blurred_img, logitori = Get_blurred_img(image_orl_for_blur,
label,
model,
resize_shape=(260, 260),
Gaussian_param=[51, 50],
Median_param=11,
blur_type='Gaussian',
use_cuda=1)
for target_index in tqdm(range(0, args.num_classes)):
mask, upsampled_mask, imgratio, curvetop, curve1, curve2, category = Integrated_Mask(
img,
blurred_img,
model,
label,
max_iterations=15,
integ_iter=20,
tv_beta=2,
l1_coeff=0.01 * 100,
tv_coeff=0.2 * 100,
size_init=8,
use_cuda=1) #
mask = upsampled_mask.cpu().detach().numpy()[0, 0]
mask = -mask + mask.max() * 2.
mask = np.maximum(mask, 0)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
mask = np.maximum(mask, args.grad_min_level)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
# mask = Image.fromarray(mask*255, mode='L').resize((args.img_size, args.img_size), Image.BILINEAR)
# mask = np.uint8(mask)
image_orl = image_orl.resize((args.img_size, args.img_size),
Image.BILINEAR)
# heatmap = np.array(heatmap)
show_cam_on_image(image_orl, mask, target_index, 'IGOS')
del model
### torchray (RISE)
model = load_backbone(args)
model = model.to(device)
model.eval()
for target_index in tqdm(range(0, args.num_classes)):
mask = rise(model, image.to(device), target_index)
mask = mask.cpu().numpy()[0, 0]
mask = np.maximum(mask, 0)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
mask = np.maximum(mask, args.grad_min_level)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
image_orl = image_orl.resize((args.img_size, args.img_size),
Image.BILINEAR)
# heatmap = np.array(heatmap)
show_cam_on_image(image_orl, mask, target_index, 'RISE')
del model
### torchray (Extremal)
model = load_backbone(args)
model = model.to(device)
model.eval()
for target_index in tqdm(range(0, args.num_classes)):
mask, _ = extremal_perturbation(model, image.to(device), target_index)
mask = mask.cpu().numpy()[0, 0]
mask = np.maximum(mask, 0)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
mask = np.maximum(mask, args.grad_min_level)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
image_orl = image_orl.resize((args.img_size, args.img_size),
Image.BILINEAR)
# heatmap = np.array(heatmap)
show_cam_on_image(image_orl, mask, target_index, 'Extremal')
del model
### IBA
model = load_backbone(args)
model = model.to(device)
model.eval()
imagenet_dir = '../../data/imagenet/ILSVRC/Data/CLS-LOC/validation'
# Add a Per-Sample Bottleneck at layer conv4_1
iba = IBA(model.layer4)
# Estimate the mean and variance of the feature map at this layer.
val_set = get_imagenet_folder(imagenet_dir)
val_loader = DataLoader(val_set,
batch_size=64,
shuffle=True,
num_workers=4)
iba.estimate(model, val_loader, n_samples=5000, progbar=True)
for target_index in tqdm(range(0, args.num_classes)):
# Closure that returns the loss for one batch
model_loss_closure = lambda x: -torch.log_softmax(
model(x.to(device)), dim=1)[:, target_index].mean()
# Explain class target for the given image
saliency_map = iba.analyze(image, model_loss_closure, beta=10)
# display result
model_loss_closure = lambda x: -torch.log_softmax(
model(x.to(device)), 1)[:, target_index].mean()
heatmap = iba.analyze(image, model_loss_closure)
mask = heatmap
mask = np.maximum(mask, 0)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
mask = np.maximum(mask, args.grad_min_level)
mask = mask - np.min(mask)
mask = mask / np.max(mask)
image_orl = image_orl.resize((args.img_size, args.img_size),
Image.BILINEAR)
# heatmap = np.array(heatmap)
show_cam_on_image(image_orl, mask, target_index, 'IBA')
# plot_saliency_map(heatmap, tensor_to_np_img(image[0]))
RESNET_CONFIG = dict(input_layer='conv1',
conv_layer='layer4',
fc_layer='fc')
MODEL_CONFIG = {**RESNET_CONFIG}
conv_layer = MODEL_CONFIG['conv_layer']
input_layer = MODEL_CONFIG['input_layer']
fc_layer = MODEL_CONFIG['fc_layer']
### torchcam
del model
model = load_backbone(args)
model = model.to(device)
model.eval()
# Hook the corresponding layer in the model
cam_extractors = [
CAM(model, conv_layer, fc_layer),
GradCAM(model, conv_layer),
GradCAMpp(model, conv_layer),
SmoothGradCAMpp(model, conv_layer, input_layer),
ScoreCAM(model, conv_layer, input_layer),
SSCAM(model, conv_layer, input_layer),
# ISSCAM(model, conv_layer, input_layer),
]
cam_extractors_names = [
'CAM',
'GradCAM',
'GradCAMpp',
'SmoothGradCAMpp',
'ScoreCAM',
'SSCAM',
# 'ISSCAM',
]
for idx, extractor in enumerate(cam_extractors):
model.zero_grad()
output1 = model(image.to(device))
output = F.softmax(output1, dim=1)
prediction_score, pred_label_idx = torch.topk(output, 1)
pred_label_idx.squeeze_()
predicted_label = str(pred_label_idx.item())
print('Predicted:', predicted_label, '(',
prediction_score.squeeze().item(), ')')
make_grad(extractor, output1, image_orl, args.grad_min_level,
cam_extractors_names[idx])
extractor.clear_hooks()