post_process.py 6.2 KB
Newer Older
F
FDInSky 已提交
1 2 3 4 5 6 7 8
import six
import os
import numpy as np
from numba import jit
from .bbox import delta2bbox, clip_bbox, expand_bbox, nms


def bbox_post_process(bboxes,
9
                      bbox_nums,
F
FDInSky 已提交
10 11 12 13 14 15 16 17
                      bbox_probs,
                      bbox_deltas,
                      im_info,
                      keep_top_k=100,
                      score_thresh=0.05,
                      nms_thresh=0.5,
                      class_nums=81,
                      bbox_reg_weights=[0.1, 0.1, 0.2, 0.2]):
18 19

    new_bboxes = [[] for _ in range(len(bbox_nums))]
F
FDInSky 已提交
20
    new_bbox_nums = [0]
21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
    st_num = 0
    end_num = 0
    for i in range(len(bbox_nums)):
        bbox_num = bbox_nums[i]
        end_num += bbox_num

        bbox = bboxes[st_num:end_num, :]  # bbox 
        bbox = bbox / im_info[i][2]  # scale
        bbox_delta = bbox_deltas[st_num:end_num, :]  # bbox delta 

        # step1: decode 
        bbox = delta2bbox(bbox_delta, bbox, bbox_reg_weights)

        # step2: clip 
        bbox = clip_bbox(bbox, im_info[i][:2] / im_info[i][2])

        # step3: nms 
F
FDInSky 已提交
38
        cls_boxes = [[] for _ in range(class_nums)]
39
        scores_n = bbox_probs[st_num:end_num, :]
F
FDInSky 已提交
40 41 42
        for j in range(1, class_nums):
            inds = np.where(scores_n[:, j] > score_thresh)[0]
            scores_j = scores_n[inds, j]
43
            rois_j = bbox[inds, j * 4:(j + 1) * 4]
F
FDInSky 已提交
44 45 46 47 48 49 50 51 52 53
            dets_j = np.hstack((scores_j[:, np.newaxis], rois_j)).astype(
                np.float32, copy=False)
            keep = nms(dets_j, nms_thresh)
            nms_dets = dets_j[keep, :]
            #add labels
            label = np.array([j for _ in range(len(keep))])
            nms_dets = np.hstack((label[:, np.newaxis], nms_dets)).astype(
                np.float32, copy=False)
            cls_boxes[j] = nms_dets

54 55
        st_num += bbox_num

F
FDInSky 已提交
56 57 58 59 60 61 62 63 64 65
        # Limit to max_per_image detections **over all classes**
        image_scores = np.hstack(
            [cls_boxes[j][:, 1] for j in range(1, class_nums)])
        if len(image_scores) > keep_top_k:
            image_thresh = np.sort(image_scores)[-keep_top_k]
            for j in range(1, class_nums):
                keep = np.where(cls_boxes[j][:, 1] >= image_thresh)[0]
                cls_boxes[j] = cls_boxes[j][keep, :]
        new_bboxes_n = np.vstack([cls_boxes[j] for j in range(1, class_nums)])
        new_bboxes[i] = new_bboxes_n
66
        new_bbox_nums.append(len(new_bboxes_n))
F
FDInSky 已提交
67 68 69 70 71 72 73 74 75
        labels = new_bboxes_n[:, 0]
        scores = new_bboxes_n[:, 1]
        boxes = new_bboxes_n[:, 2:]
    new_bboxes = np.vstack([new_bboxes[k] for k in range(len(bbox_nums) - 1)])
    new_bbox_nums = np.array(new_bbox_nums)
    return new_bbox_nums, new_bboxes


@jit
76 77
def mask_post_process(bboxes, bbox_nums, masks, im_info, resolution=14):
    scale = (resolution + 2.0) / resolution
F
FDInSky 已提交
78 79
    boxes = bboxes[:, 2:]
    labels = bboxes[:, 0]
80
    segms_results = [[] for _ in range(len(bbox_nums))]
F
FDInSky 已提交
81
    sum = 0
82 83 84 85 86 87
    st_num = 0
    end_num = 0
    for i in range(len(bbox_nums)):
        bbox_num = bbox_nums[i]
        end_num += bbox_num

F
FDInSky 已提交
88
        cls_segms = []
89 90 91 92
        boxes_n = boxes[st_num:end_num]
        labels_n = labels[st_num:end_num]
        masks_n = masks[st_num:end_num]

F
FDInSky 已提交
93 94 95 96 97
        im_h = int(round(im_info[i][0] / im_info[i][2]))
        im_w = int(round(im_info[i][1] / im_info[i][2]))
        boxes_n = expand_boxes(boxes_n, scale)
        boxes_n = boxes_n.astype(np.int32)
        padded_mask = np.zeros((M + 2, M + 2), dtype=np.float32)
98
        for j in range(len(boxes_n)):
F
FDInSky 已提交
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
            class_id = int(labels_n[j])
            padded_mask[1:-1, 1:-1] = masks_n[j, class_id, :, :]

            ref_box = boxes_n[j, :]
            w = ref_box[2] - ref_box[0] + 1
            h = ref_box[3] - ref_box[1] + 1
            w = np.maximum(w, 1)
            h = np.maximum(h, 1)

            mask = cv2.resize(padded_mask, (w, h))
            mask = np.array(mask > cfg.mrcnn_thresh_binarize, dtype=np.uint8)
            im_mask = np.zeros((im_h, im_w), dtype=np.uint8)

            x_0 = max(ref_box[0], 0)
            x_1 = min(ref_box[2] + 1, im_w)
            y_0 = max(ref_box[1], 0)
            y_1 = min(ref_box[3] + 1, im_h)
            im_mask[y_0:y_1, x_0:x_1] = mask[(y_0 - ref_box[1]):(y_1 - ref_box[
                1]), (x_0 - ref_box[0]):(x_1 - ref_box[0])]
            sum += im_mask.sum()
            rle = mask_util.encode(
                np.array(
                    im_mask[:, :, np.newaxis], order='F'))[0]
            cls_segms.append(rle)
        segms_results[i] = np.array(cls_segms)[:, np.newaxis]
124
    segms_results = np.vstack([segms_results[k] for k in range(len(bbox_nums))])
F
FDInSky 已提交
125 126 127 128 129
    bboxes = np.hstack([segms_results, bboxes])
    return bboxes[:, :3]


@jit
130 131
def get_det_res(bboxes, bbox_nums, image_id, num_id_to_cat_id_map,
                batch_size=1):
F
FDInSky 已提交
132 133
    det_res = []
    k = 0
134
    for i in range(len(bbox_nums)):
F
FDInSky 已提交
135
        image_id = int(image_id[i][0])
136 137 138 139 140 141
        image_width = int(image_shape[i][1])
        image_height = int(image_shape[i][2])

        det_nums = bbox_nums[i]
        for j in range(det_nums):
            dt = bboxes[k]
F
FDInSky 已提交
142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158
            k = k + 1
            num_id, score, xmin, ymin, xmax, ymax = dt.tolist()
            category_id = num_id_to_cat_id_map[num_id]
            w = xmax - xmin + 1
            h = ymax - ymin + 1
            bbox = [xmin, ymin, w, h]
            dt_res = {
                'image_id': image_id,
                'category_id': category_id,
                'bbox': bbox,
                'score': score
            }
            det_res.append(dt_res)
    return det_res


@jit
159
def get_seg_res(masks, mask_nums, image_id, num_id_to_cat_id_map):
F
FDInSky 已提交
160 161
    seg_res = []
    k = 0
162
    for i in range(len(mask_nums)):
F
FDInSky 已提交
163
        image_id = int(image_id[i][0])
164 165 166
        det_nums = mask_nums[i]
        for j in range(det_nums):
            dt = masks[k]
F
FDInSky 已提交
167 168 169 170 171 172 173 174 175 176 177 178 179 180
            k = k + 1
            sg, num_id, score = dt.tolist()
            cat_id = num_id_to_cat_id_map[num_id]
            if six.PY3:
                if 'counts' in sg:
                    sg['counts'] = sg['counts'].decode("utf8")
            sg_res = {
                'image_id': image_id,
                'category_id': cat_id,
                'segmentation': sg,
                'score': score
            }
            seg_res.append(sg_res)
    return seg_res