| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
|
|
| from abc import ABC, abstractmethod |
| from typing import List, Optional, Tuple, Union |
| from datasets import load_dataset |
| import torch |
| import torch.nn as nn |
| from torch.nn import CrossEntropyLoss |
| import numpy as np |
| import copy |
| import os |
| import sys |
| from PIL import Image |
| import requests |
| from io import BytesIO |
|
|
| dir_path = os.path.dirname(os.path.realpath(__file__)) |
| sys.path.insert(0, dir_path) |
|
|
| from transformers import AutoConfig, AutoModelForCausalLM, AutoTokenizer, CLIPImageProcessor, LlamaConfig, LlamaModel, LlamaForCausalLM |
| from transformers.modeling_outputs import CausalLMOutputWithPast |
|
|
| from .configuration_mplug_owl2 import MPLUGOwl2Config, MplugOwlVisionConfig, MplugOwlVisualAbstractorConfig |
| from .visual_encoder import MplugOwlVisionModel, MplugOwlVisualAbstractorModel |
| from .modeling_llama2 import replace_llama_modality_adaptive |
| from .modeling_attn_mask_utils import _prepare_4d_causal_attention_mask |
|
|
| IGNORE_INDEX = -100 |
| IMAGE_TOKEN_INDEX = -200 |
| DEFAULT_IMAGE_TOKEN = "<|image|>" |
| from icecream import ic |
|
|
| def tokenizer_image_token(prompt, tokenizer, image_token_index=IMAGE_TOKEN_INDEX, return_tensors=None): |
| prompt_chunks = [tokenizer(chunk).input_ids if len(chunk) > 0 else [] for chunk in prompt.split(DEFAULT_IMAGE_TOKEN)] |
|
|
| def insert_separator(X, sep): |
| return [ele for sublist in zip(X, [sep]*len(X)) for ele in sublist][:-1] |
|
|
| input_ids = [] |
| offset = 0 |
| if len(prompt_chunks) > 0 and len(prompt_chunks[0]) > 0 and prompt_chunks[0][0] == tokenizer.bos_token_id: |
| offset = 1 |
| input_ids.append(prompt_chunks[0][0]) |
|
|
| for x in insert_separator(prompt_chunks, [image_token_index] * (offset + 1)): |
| input_ids.extend(x[offset:]) |
|
|
| if return_tensors is not None: |
| if return_tensors == 'pt': |
| return torch.tensor(input_ids, dtype=torch.long) |
| raise ValueError(f'Unsupported tensor type: {return_tensors}') |
| return input_ids |
|
|
| def expand2square(pil_img, background_color): |
| from PIL import Image |
| width, height = pil_img.size |
| if width == height: |
| return pil_img |
| elif width > height: |
| result = Image.new(pil_img.mode, (width, width), background_color) |
| result.paste(pil_img, (0, (width - height) // 2)) |
| return result |
| else: |
| result = Image.new(pil_img.mode, (height, height), background_color) |
| result.paste(pil_img, ((height - width) // 2, 0)) |
| return result |
|
|
| def norm_cdf(x): |
| return 0.5 * (1 + torch.erf(x / torch.sqrt(torch.tensor(2.0)))) |
|
|
| def optimize_score_map_pytorch_cuda(c, seed=0, original_seed=20020, num_iterations=100): |
| device = torch.device("cuda" if torch.cuda.is_available() else "cpu") |
| |
| torch.manual_seed(seed) |
| np.random.seed(seed) |
|
|
| c = torch.tensor(c, dtype=torch.float32, device=device, requires_grad=False) |
| initial_scores = torch.rand(c.shape[0], device=device, requires_grad=True) |
| |
| optimizer = torch.optim.Adam([initial_scores], lr=0.1) |
|
|
| for _ in range(num_iterations): |
| optimizer.zero_grad() |
| sum_log_diff = torch.sum(c * torch.log(torch.maximum(norm_cdf(initial_scores[:, None] - initial_scores), torch.tensor(1e-6, device=device)))) |
| sum_squares = torch.sum(initial_scores ** 2) / 2 |
|
|
| loss = -(sum_log_diff - sum_squares) |
| loss.backward() |
| optimizer.step() |
| |
| optimized_scores = initial_scores.detach().cpu().numpy() |
| min_score, max_score = np.min(optimized_scores), np.max(optimized_scores) |
| |
| |
| scaled_scores = 100 * (optimized_scores - min_score) / (max_score - min_score) |
| |
| |
| np.random.seed(original_seed) |
| return torch.tensor(scaled_scores[-1], device=device) |
|
|
| def softmax(logits): |
| |
| probs = np.exp(logits) / np.sum(np.exp(logits)) |
| return probs |
| |
|
|
| def update_matrix(anchor_matrix, scores, indices): |
| n = anchor_matrix.shape[0] |
| new_row = np.zeros((1, n)) |
| new_col = np.zeros((n + 1, 1)) |
| new_row[0, indices] = scores |
| new_col[indices, 0] = 1-scores |
| anchor_matrix = np.vstack([anchor_matrix, new_row]) |
| anchor_matrix = np.hstack([anchor_matrix, new_col]) |
|
|
| return anchor_matrix |
| |
|
|
| class MPLUGOwl2MetaModel: |
| def __init__(self, config): |
| super(MPLUGOwl2MetaModel, self).__init__(config) |
| self.vision_model = MplugOwlVisionModel( |
| MplugOwlVisionConfig(**config.visual_config["visual_model"]) |
| ) |
| self.visual_abstractor = MplugOwlVisualAbstractorModel( |
| MplugOwlVisualAbstractorConfig(**config.visual_config["visual_abstractor"]), config.hidden_size |
| ) |
| |
| def get_vision_tower(self): |
| vision_model = getattr(self, 'vision_model', None) |
| if type(vision_model) is list: |
| vision_model = vision_model[0] |
| return vision_model |
|
|
| def get_visual_abstractor(self): |
| visual_abstractor = getattr(self, 'visual_abstractor', None) |
| if type(visual_abstractor) is list: |
| visual_abstractor = visual_abstractor[0] |
| return visual_abstractor |
|
|
|
|
| class MPLUGOwl2MetaForCausalLM(ABC): |
| @abstractmethod |
| def get_model(self): |
| pass |
|
|
| def encode_images(self, images): |
| image_features = self.get_model().vision_model(images).last_hidden_state |
| image_features = self.get_model().visual_abstractor(encoder_hidden_states=image_features).last_hidden_state |
| return image_features |
|
|
| def prepare_inputs_labels_for_multimodal( |
| self, input_ids, attention_mask, past_key_values, labels, images |
| ): |
| if images is None or input_ids.shape[1] == 1: |
| if past_key_values is not None and images is not None and input_ids.shape[1] == 1: |
| attention_mask = torch.ones((attention_mask.shape[0], past_key_values[-1][-1].shape[-2] + 1), dtype=attention_mask.dtype, device=attention_mask.device) |
| multiway_indices = torch.zeros_like(input_ids).long().to(self.device) |
| return input_ids, multiway_indices, attention_mask, past_key_values, None, labels |
| |
| if type(images) is list or images.ndim == 5: |
| concat_images = torch.cat([image for image in images], dim=0) |
| image_features = self.encode_images(concat_images) |
| split_sizes = [image.shape[0] for image in images] |
| image_features = torch.split(image_features, split_sizes, dim=0) |
| image_features = [x.flatten(0, 1) for x in image_features] |
| else: |
| image_features = self.encode_images(images) |
|
|
| new_input_embeds = [] |
| new_modality_indicators = [] |
| new_labels = [] if labels is not None else None |
| cur_image_idx = 0 |
| for batch_idx, cur_input_ids in enumerate(input_ids): |
| if (cur_input_ids == IMAGE_TOKEN_INDEX).sum() == 0: |
| |
| |
| half_len = cur_input_ids.shape[0] // 2 |
| cur_image_features = image_features[cur_image_idx] |
| cur_input_embeds_1 = self.get_model().embed_tokens(cur_input_ids[:half_len]) |
| cur_input_embeds_2 = self.get_model().embed_tokens(cur_input_ids[half_len:]) |
| cur_input_embeds = torch.cat([cur_input_embeds_1, cur_image_features[0:0], cur_input_embeds_2], dim=0) |
| new_input_embeds.append(cur_input_embeds) |
| |
| cur_modality_indicators = torch.zeros(len(cur_input_embeds)).long().to(self.device) |
| new_modality_indicators.append(cur_modality_indicators) |
| if labels is not None: |
| new_labels.append(labels[batch_idx]) |
| cur_image_idx += 1 |
| continue |
| image_token_indices = torch.where(cur_input_ids == IMAGE_TOKEN_INDEX)[0] |
| cur_new_input_embeds = [] |
| cur_modality_indicators = [] |
| if labels is not None: |
| cur_labels = labels[batch_idx] |
| cur_new_labels = [] |
| assert cur_labels.shape == cur_input_ids.shape |
| while image_token_indices.numel() > 0: |
| |
| cur_image_features = image_features[cur_image_idx] |
| image_token_start = image_token_indices[0] |
| cur_new_input_embeds.append(self.get_model().embed_tokens(cur_input_ids[:image_token_start])) |
| cur_new_input_embeds.append(cur_image_features) |
| |
| |
| assert image_token_start == len(cur_input_ids[:image_token_start]) |
| cur_modality_indicators.append(torch.zeros(len(cur_input_ids[:image_token_start])).long()) |
| cur_modality_indicators.append(torch.ones(len(cur_image_features)).long()) |
| |
| if labels is not None: |
| cur_new_labels.append(cur_labels[:image_token_start]) |
| cur_new_labels.append(torch.full((cur_image_features.shape[0],), IGNORE_INDEX, device=labels.device, dtype=labels.dtype)) |
| cur_labels = cur_labels[image_token_start+1:] |
| cur_image_idx += 1 |
| cur_input_ids = cur_input_ids[image_token_start+1:] |
| image_token_indices = torch.where(cur_input_ids == IMAGE_TOKEN_INDEX)[0] |
| if cur_input_ids.numel() > 0: |
| cur_new_input_embeds.append(self.get_model().embed_tokens(cur_input_ids)) |
| cur_modality_indicators.append(torch.zeros(len(cur_input_ids)).long()) |
| if labels is not None: |
| cur_new_labels.append(cur_labels) |
| cur_new_input_embeds = [x.to(device=self.device) for x in cur_new_input_embeds] |
| cur_new_input_embeds = torch.cat(cur_new_input_embeds, dim=0) |
| new_input_embeds.append(cur_new_input_embeds) |
| |
| |
| cur_modality_indicators = [x.to(device=self.device) for x in cur_modality_indicators] |
| cur_modality_indicators = torch.cat(cur_modality_indicators, dim=0) |
| new_modality_indicators.append(cur_modality_indicators) |
| |
| |
| if labels is not None: |
| cur_new_labels = torch.cat(cur_new_labels, dim=0) |
| new_labels.append(cur_new_labels) |
|
|
| if any(x.shape != new_input_embeds[0].shape for x in new_input_embeds): |
| max_len = max(x.shape[0] for x in new_input_embeds) |
| |
| |
| new_input_embeds_align = [] |
| for cur_new_embed in new_input_embeds: |
| cur_new_embed = torch.cat((cur_new_embed, torch.zeros((max_len - cur_new_embed.shape[0], cur_new_embed.shape[1]), dtype=cur_new_embed.dtype, device=cur_new_embed.device)), dim=0) |
| new_input_embeds_align.append(cur_new_embed) |
| new_input_embeds = torch.stack(new_input_embeds_align, dim=0) |
| |
| |
| new_modality_indicators_align = [] |
| for cur_modality_indicator in new_modality_indicators: |
| cur_new_embed = torch.cat((cur_modality_indicator, torch.zeros(max_len - cur_modality_indicator.shape[0], dtype=cur_modality_indicator.dtype, device=cur_modality_indicator.device)), dim=0) |
| new_modality_indicators_align.append(cur_new_embed) |
| new_modality_indicators = torch.stack(new_modality_indicators_align, dim=0) |
| |
| |
| if labels is not None: |
| new_labels_align = [] |
| _new_labels = new_labels |
| for cur_new_label in new_labels: |
| cur_new_label = torch.cat((cur_new_label, torch.full((max_len - cur_new_label.shape[0],), IGNORE_INDEX, dtype=cur_new_label.dtype, device=cur_new_label.device)), dim=0) |
| new_labels_align.append(cur_new_label) |
| new_labels = torch.stack(new_labels_align, dim=0) |
| |
| |
| if attention_mask is not None: |
| new_attention_mask = [] |
| for cur_attention_mask, cur_new_labels, cur_new_labels_align in zip(attention_mask, _new_labels, new_labels): |
| new_attn_mask_pad_left = torch.full((cur_new_labels.shape[0] - labels.shape[1],), True, dtype=attention_mask.dtype, device=attention_mask.device) |
| new_attn_mask_pad_right = torch.full((cur_new_labels_align.shape[0] - cur_new_labels.shape[0],), False, dtype=attention_mask.dtype, device=attention_mask.device) |
| cur_new_attention_mask = torch.cat((new_attn_mask_pad_left, cur_attention_mask, new_attn_mask_pad_right), dim=0) |
| new_attention_mask.append(cur_new_attention_mask) |
| attention_mask = torch.stack(new_attention_mask, dim=0) |
| assert attention_mask.shape == new_labels.shape |
| else: |
| new_input_embeds = torch.stack(new_input_embeds, dim=0) |
| new_modality_indicators = torch.stack(new_modality_indicators, dim=0) |
| if labels is not None: |
| new_labels = torch.stack(new_labels, dim=0) |
|
|
| if attention_mask is not None: |
| new_attn_mask_pad_left = torch.full((attention_mask.shape[0], new_input_embeds.shape[1] - input_ids.shape[1]), True, dtype=attention_mask.dtype, device=attention_mask.device) |
| attention_mask = torch.cat((new_attn_mask_pad_left, attention_mask), dim=1) |
| assert attention_mask.shape == new_input_embeds.shape[:2] |
| return None, new_modality_indicators, attention_mask, past_key_values, new_input_embeds, new_labels |
|
|
|
|
|
|
| class MPLUGOwl2LlamaModel(MPLUGOwl2MetaModel, LlamaModel): |
| config_class = MPLUGOwl2Config |
|
|
| def __init__(self, config: MPLUGOwl2Config): |
| super(MPLUGOwl2LlamaModel, self).__init__(config) |
|
|
|
|
| class MPLUGOwl2LlamaForCausalLM(LlamaForCausalLM, MPLUGOwl2MetaForCausalLM): |
| config_class = MPLUGOwl2Config |
|
|
| def __init__(self, config): |
| super(LlamaForCausalLM, self).__init__(config) |
| self.model = MPLUGOwl2LlamaModel(config) |
| self.tokenizer = AutoTokenizer.from_pretrained("q-future/Compare2Score") |
| self.image_processor = CLIPImageProcessor.from_pretrained("q-future/Compare2Score") |
|
|
| self.lm_head = nn.Linear(config.hidden_size, config.vocab_size, bias=False) |
| self.preferential_ids_ = [id_[1] for id_ in self.tokenizer(["inferior", "worse", "similar", "better", "superior"])["input_ids"]] |
| self.anchor_images = load_dataset("VQA-CityU/Anchor_images") |
| |
| self.weight_tensor = np.array([0., 0.25, 0.5, 0.75, 1.], dtype=np.float16) |
| self.anchor_matrix = np.array( |
| [[5.0000000e-01, 2.5912809e-01, 3.3130276e-04, 1.6087297e-06, 1.1803027e-09], |
| [7.4087191e-01, 5.0000000e-01, 2.4985345e-01, 9.9954158e-02, 1.8675303e-08], |
| [9.9966872e-01, 7.5014657e-01, 5.0000000e-01, 4.9968880e-01, 2.4852838e-01], |
| [9.9999839e-01, 9.0004587e-01, 5.0031120e-01, 5.0000000e-01, 2.5400183e-01], |
| [1.0000000e+00, 1.0000000e+00, 7.5147164e-01, 7.4599814e-01, 5.0000000e-01]], |
| dtype=np.float32) |
| anchor_intervals = 5 |
| num_anchor_image_per_interval = 1 |
| num_anchor_image = anchor_intervals * num_anchor_image_per_interval |
| self.anchor_indices = np.arange(0,num_anchor_image) |
| |
| self.post_init() |
| |
|
|
| def get_model(self): |
| return self.model |
| |
| def download_image(self, url): |
| response = requests.get(url) |
| return Image.open(BytesIO(response.content)).convert('RGB') |
|
|
| def load_image(self, path): |
| if path.startswith('http://') or path.startswith('https://'): |
| return self.download_image(path) |
| return Image.open(path).convert('RGB') |
| |
| def score(self, image): |
| prompt = "USER: <|image|> <|image|> Compared with the first image, what is your quality rating for second image? \nASSISTANT: The quality of the second image is" |
| input_ids = tokenizer_image_token(prompt, self.tokenizer, IMAGE_TOKEN_INDEX, return_tensors='pt').unsqueeze(0).to(self.device) |
| |
| anchor_images = [item['image'] for item in self.anchor_images['train']] |
| |
| probabilities = [] |
| for index in self.anchor_indices: |
| anchor_image = anchor_images[index] |
| |
| images = [anchor_image, image] |
| images = [expand2square(img, tuple(int(x*255) for x in self.image_processor.image_mean)) for img in images] |
| image_tensor = self.image_processor.preprocess(images, return_tensors='pt')['pixel_values'].half().to(self.device) |
| |
| with torch.inference_mode(): |
| output_logits = self(input_ids, images=image_tensor)["logits"][:, -1, self.preferential_ids_] |
| output_logits = output_logits.cpu().detach().numpy() / 100 |
| probabilities.append(np.dot(softmax(output_logits), self.weight_tensor)) |
| updated_matrix = update_matrix(self.anchor_matrix, np.squeeze(np.array(probabilities)), self.anchor_indices) |
| score = optimize_score_map_pytorch_cuda(updated_matrix, seed=0, original_seed=20020, num_iterations=100) |
| return score |
|
|
| def forward( |
| self, |
| input_ids: torch.LongTensor = None, |
| |
| attention_mask: Optional[torch.Tensor] = None, |
| past_key_values: Optional[List[torch.FloatTensor]] = None, |
| inputs_embeds: Optional[torch.FloatTensor] = None, |
| labels: Optional[torch.LongTensor] = None, |
| use_cache: Optional[bool] = None, |
| output_attentions: Optional[bool] = None, |
| output_hidden_states: Optional[bool] = None, |
| images: Optional[torch.FloatTensor] = None, |
| return_dict: Optional[bool] = None, |
| ) -> Union[Tuple, CausalLMOutputWithPast]: |
| output_attentions = output_attentions if output_attentions is not None else self.config.output_attentions |
| output_hidden_states = ( |
| output_hidden_states if output_hidden_states is not None else self.config.output_hidden_states |
| ) |
| return_dict = return_dict if return_dict is not None else self.config.use_return_dict |
| input_ids, modality_indicators, attention_mask, past_key_values, inputs_embeds, labels = \ |
| self.prepare_inputs_labels_for_multimodal(input_ids, attention_mask, past_key_values, labels, images) |
|
|
| |
| outputs = self.model( |
| input_ids=input_ids, |
| modality_indicators=modality_indicators, |
| attention_mask=attention_mask, |
| past_key_values=past_key_values, |
| inputs_embeds=inputs_embeds, |
| use_cache=use_cache, |
| output_attentions=output_attentions, |
| output_hidden_states=output_hidden_states, |
| return_dict=return_dict |
| ) |
|
|
| hidden_states = outputs[0] |
| logits = self.lm_head(hidden_states) |
|
|
| loss = None |
| if labels is not None: |
| |
| shift_logits = logits[..., :-1, :].contiguous() |
| shift_labels = labels[..., 1:].contiguous() |
| |
| loss_fct = CrossEntropyLoss() |
| shift_logits = shift_logits.view(-1, self.config.vocab_size) |
| shift_labels = shift_labels.view(-1) |
| |
| shift_labels = shift_labels.to(shift_logits.device) |
| loss = loss_fct(shift_logits, shift_labels) |
|
|
| if not return_dict: |
| output = (logits,) + outputs[1:] |
| return (loss,) + output if loss is not None else output |
|
|
| return CausalLMOutputWithPast( |
| loss=loss, |
| logits=logits, |
| past_key_values=outputs.past_key_values, |
| hidden_states=outputs.hidden_states, |
| attentions=outputs.attentions, |
| ) |
|
|
| def prepare_inputs_for_generation( |
| self, input_ids, past_key_values=None, attention_mask=None, inputs_embeds=None, **kwargs |
| ): |
| if past_key_values: |
| input_ids = input_ids[:, -1:] |
|
|
| |
| if inputs_embeds is not None and past_key_values is None: |
| model_inputs = {"inputs_embeds": inputs_embeds} |
| else: |
| model_inputs = {"input_ids": input_ids} |
|
|
| model_inputs.update( |
| { |
| "past_key_values": past_key_values, |
| "use_cache": kwargs.get("use_cache"), |
| "attention_mask": attention_mask, |
| "images": kwargs.get("images", None), |
| } |
| ) |
| return model_inputs |
|
|
| AutoConfig.register("mplug_owl2", MPLUGOwl2Config) |
| AutoModelForCausalLM.register(MPLUGOwl2Config, MPLUGOwl2LlamaForCausalLM) |
|
|
| replace_llama_modality_adaptive() |
|
|
| if __name__ == "__main__": |
| |
| from icecream import ic |
| |
| |
| model = AutoModelForCausalLM.from_pretrained('VQA-CityU/Compare2Score_1', trust_remote_code=True, |
| torch_dtype=torch.float16, device_map="auto") |
| |
| model.score("/home/zhw/IQA/code/NeurIPS24/Q-Align/playground/data/TID2013/distorted_images/i01_01_5.bmp") |
| url = "https://raw.githubusercontent.com/Q-Future/Q-Align/main/fig/singapore_flyer.jpg" |
| model.score(url) |
| |