主頁 > 後端開發 > 一、使用paddleseg套件對遙感影像預測(基礎)

一、使用paddleseg套件對遙感影像預測(基礎)

2021-02-11 12:13:09 後端開發

一、使用paddleseg套件對遙感影像預測(基礎)


目前paddleseg套件中的predict.py代碼檔案還不支持直接對遙感影像(大圖)做預測,或者說把遙感大圖直接丟進predict.py,它的預測效果非常差,

基于以上問題,本文結合paddleseg中predict.py原始碼和這篇博文代碼(遙感語意分割切圖預測之后再拼接)重新寫了predict.py代碼,希望可以幫助到使用飛槳框架做遙感影像語意分割的朋友,所以這里需要你會使用paddleseg套件或者對paddleseg原始碼有所了解,這里有位博主寫了一系列有關paddleseg原始碼的文章,值得參考學習(人工智能研習社),

重新寫的predict.py代碼主要分為四個部分
讀取和裁剪遙感大圖 網路模型推理預測小圖塊 拼接小圖塊預測結果 拼接結果寫入檔案

1、讀取待預測遙感大圖,將遙感大圖裁剪成小圖塊,這里裁剪的小圖塊相鄰之間不設定重疊度,小圖塊大小為256x256,
本部分代碼如下:

#讀取需要預測的遙感大圖img_lists[local_rank][local_rank] = /home/aistudio/data/data70483/img.png
ori_image=cv2.imread(img_lists[local_rank][local_rank])
h_step = ori_image.shape[0] // 256  #高度步數
w_step = ori_image.shape[1] // 256  #寬度步數

h_rest = -(ori_image.shape[0] - 256 * h_step)  #剩余行數
w_rest = -(ori_image.shape[1] - 256 * w_step)  #剩余列數

seg_list = []  #小圖塊的串列
predict_list = []#預測小圖塊結果的串列
# 回圈切圖
for h in range(h_step):
    for w in range(w_step):
        # 劃窗采樣
        image_sample = ori_image[(h * 256):(h * 256 + 256),
                       (w * 256):(w * 256 + 256), :]
        seg_list.append(image_sample)
    seg_list.append(ori_image[(h * 256):(h * 256 + 256), -256:, :])
for w in range(w_step - 1):
    seg_list.append(ori_image[-256:, (w * 256):(w * 256 + 256), :])
seg_list.append(ori_image[-256:, -256:, :])

2、利用網路模型推理預測小圖塊,這里的代碼改動不多,但是需要將img_lists[local_rank]引數改成存盤小圖塊的串列seg_list,其他引數的設定根據需要而定,在這里本文只對小圖塊做最普通的推理預測,既不做多尺度預測、也不做滑窗預測(多尺度和滑窗預測是原本predict.py的功能,當然在這里我們也可以用),
本部分代碼如下:

progbar_pred = progbar.Progbar(target=len(seg_list), verbose=1)
with paddle.no_grad():
    for i, im in enumerate(seg_list):
        ori_shape = im.shape[:2]  #原始圖片形狀(h,w)
        im, _ = transforms(im)    #im.shape(3, 256, 256)    _為None
        im = im[np.newaxis, ...]   #im.shape(1,3,256,256)
        im = paddle.to_tensor(im)

        if False:
            pred = infer.aug_inference(
                model,
                im,
                ori_shape=ori_shape,
                transforms=transforms.transforms,
                scales=scales,
                flip_horizontal=flip_horizontal,
                flip_vertical=flip_vertical,
                is_slide=is_slide,
                stride=None,
                crop_size=None)
        else:
            pred = infer.inference(
                model,
                im,
                ori_shape=ori_shape,
                transforms=transforms.transforms,
                is_slide=False,
                stride=None,
                crop_size=None)
        pred = paddle.squeeze(pred)  #該OP會洗掉輸入Tensor的Shape中尺寸為1的維度,查看pred的形狀 應該剩下[h,w]
        pred = pred.numpy().astype('uint8')
        predict_list.append(pred)
        progbar_pred.update(i + 1)

3、將小圖塊的預測結果進行拼接,這里的拼接思想很簡單,就是按照裁剪的順序進行拼接,
本部分代碼如下:

count_temp = 0
tmp = np.ones([ori_image.shape[0], ori_image.shape[1]])
for h in range(h_step):
    for w in range(w_step):
        tmp[
        h * 256:(h + 1) * 256,
        w * 256:(w + 1) * 256
        ] = predict_list[count_temp]
        count_temp += 1
    tmp[h * 256:(h + 1) * 256, w_rest:] = predict_list[count_temp][:, w_rest:]
    count_temp += 1
for w in range(w_step - 1):
    tmp[h_rest:, (w * 256):(w * 256 + 256)] = predict_list[count_temp][h_rest:, :]
    count_temp += 1
tmp[-257:-1, -257:-1] = predict_list[count_temp][:, :]

4、將拼接結果 tmp 寫入影像檔案中,這里使用了原先predict.py的寫入函式,只是將函式中pred引數改成了tmp,需要注意的是一定要將tmp變數提前轉換為uint8型別,不然程式會報錯,
本部分代碼如下:

tmp = tmp.astype('uint8')
# save added image
added_image = utils.visualize.visualize(args.image_path,tmp, weight=0.6)
added_image_path = os.path.join(added_saved_dir, im_file)
mkdir(added_image_path)
cv2.imwrite(added_image_path, added_image)
# save pseudo color prediction
pred_mask = utils.visualize.get_pseudo_color_map(tmp)
pred_saved_path = os.path.join(pred_saved_dir,
                                im_file.rsplit(".")[0] + ".png")
mkdir(pred_saved_path)
pred_mask.save(pred_saved_path)

到這里代碼的主體部分基本上搞定了,值得注意的是paddleseg套件中predict.py會從paddleseg.core 呼叫predict.py,而本文為了方便移植代碼,就將兩個predict.py寫成了一個predict.py,

當時寫的第一個版本predict.py將裁剪的小圖塊尺寸設定為了256,同時將inference有些引數都設定死了,所以不推薦直接copy使用,僅作為參考學習,第一個版本predict.py完整代碼如下:

import sys
import argparse
import os
import paddle
from paddleseg.cvlibs import manager, Config
from paddleseg.utils import get_sys_env, logger
import math
import cv2
import numpy as np
from paddleseg import utils
from paddleseg.core import infer
from paddleseg.utils import progbar

def mkdir(path):
    sub_dir = os.path.dirname(path)  #去掉檔案名,回傳目錄
    if not os.path.exists(sub_dir):
        os.makedirs(sub_dir)

def partition_list(arr, m):
    """split the list 'arr' into m pieces"""
    n = int(math.ceil(len(arr) / float(m)))
    return [arr[i:i + n] for i in range(0, len(arr), n)]

def parse_args():
    parser = argparse.ArgumentParser(description='Model prediction')

    # params of prediction
    parser.add_argument(
        "--config", dest="cfg", help="The config file.", default=None, type=str)
    parser.add_argument(
        '--model_path',
        dest='model_path',
        help='The path of model for prediction',
        type=str,
        default=None)
    parser.add_argument(
        '--image_path',
        dest='image_path',
        help=
        'The path of image, it can be a file or a directory including images',
        type=str,
        default=None)
    parser.add_argument(
        '--save_dir',
        dest='save_dir',
        help='The directory for saving the predicted results',
        type=str,
        default='./output/result')

    # augment for prediction
    parser.add_argument(
        '--aug_pred',
        dest='aug_pred',
        help='Whether to use mulit-scales and flip augment for prediction',
        action='store_true')
    parser.add_argument(
        '--scales',
        dest='scales',
        nargs='+',
        help='Scales for augment',
        type=float,
        default=1.0)
    parser.add_argument(
        '--flip_horizontal',
        dest='flip_horizontal',
        help='Whether to use flip horizontally augment',
        action='store_true')
    parser.add_argument(
        '--flip_vertical',
        dest='flip_vertical',
        help='Whether to use flip vertically augment',
        action='store_true')

    # sliding window prediction
    parser.add_argument(
        '--is_slide',
        dest='is_slide',
        help='Whether to prediction by sliding window',
        action='store_true')
    parser.add_argument(
        '--crop_size',
        dest='crop_size',
        nargs=2,
        help=
        'The crop size of sliding window, the first is width and the second is height.',
        type=int,
        default=None)
    parser.add_argument(
        '--stride',
        dest='stride',
        nargs=2,
        help=
        'The stride of sliding window, the first is width and the second is height.',
        type=int,
        default=None)

    return parser.parse_args()

def get_image_list(image_path):
    """Get image list"""
    valid_suffix = [
        '.JPEG', '.jpeg', '.JPG', '.jpg', '.BMP', '.bmp', '.PNG', '.png' ,'.tif'
    ]
    image_list = []
    image_dir = None
    if os.path.isfile(image_path):
        if os.path.splitext(image_path)[-1] in valid_suffix:
            image_list.append(image_path)
    elif os.path.isdir(image_path):
        image_dir = image_path
        for root, dirs, files in os.walk(image_path):  #root=image_path
            for f in files:
                if os.path.splitext(f)[-1] in valid_suffix:
                    image_list.append(os.path.join(root, f))
    else:
        raise FileNotFoundError(
            '`--image_path` is not found. it should be an image file or a directory including images'
        )

    if len(image_list) == 0:
        raise RuntimeError('There are not image file in `--image_path`')

    return image_list, image_dir   #回傳測驗檔案串列


def main(args):
    env_info = get_sys_env()
    place = 'gpu' if env_info['Paddle compiled with cuda'] and env_info[
        'GPUs used'] else 'cpu'

    paddle.set_device(place)
    if not args.cfg:
        raise RuntimeError('No configuration file specified.')

    cfg = Config(args.cfg)
    val_dataset = cfg.val_dataset   
    if not val_dataset:
        raise RuntimeError(
            'The verification dataset is not specified in the configuration file.'
        )

    msg = '\n---------------Config Information---------------\n'
    msg += str(cfg)
    msg += '------------------------------------------------'
    logger.info(msg)

    model = cfg.model
    transforms = val_dataset.transforms
    #image_list, image_dir = get_image_list('data/UAV_seg/images')   
    image_list, image_dir = get_image_list(args.image_path)#需要傳入args.image_path引數   這個引數可以是測驗圖片的路徑,也可以是單張圖片的路徑
    
    model_path=args.model_path,  #傳入訓練模型的路徑
    save_dir=args.save_dir,
    aug_pred=False,
    scales=1.0,
    flip_horizontal=True,
    flip_vertical=False,
    is_slide=False,
    stride=None,
    crop_size=None

    para_state_dict = paddle.load(model_path[0])
    model.set_dict(para_state_dict)
    model.eval()
    nranks = paddle.distributed.get_world_size()
    local_rank = paddle.distributed.get_rank()
    if nranks > 1:
        img_lists = partition_list(image_list, nranks)
    else:
        img_lists = [image_list]  #串列的串列   img_lists[0]  ->串列

    added_saved_dir = os.path.join(save_dir[0], 'added_prediction')  #偽彩色和原圖疊加
    pred_saved_dir = os.path.join(save_dir[0], 'pseudo_color_prediction')  #偽彩色預測結果

    logger.info("Start to predict...")
    

##############################  1、裁剪遙感大圖  ########################
    #讀取需要預測的遙感大圖img_lists[local_rank][local_rank] = /home/aistudio/data/data70483/img.png
    ori_image=cv2.imread(img_lists[local_rank][local_rank])
    h_step = ori_image.shape[0] // 256  #高度步數
    w_step = ori_image.shape[1] // 256  #寬度步數

    h_rest = -(ori_image.shape[0] - 256 * h_step)  #剩余行數
    w_rest = -(ori_image.shape[1] - 256 * w_step)  #剩余列數

    seg_list = []  #由遙感大圖裁剪成小圖塊的串列
    predict_list = []#預測小圖塊結果的串列
    # 回圈切圖
    for h in range(h_step):
        for w in range(w_step):
            # 劃窗采樣
            image_sample = ori_image[(h * 256):(h * 256 + 256),
                           (w * 256):(w * 256 + 256), :]
            seg_list.append(image_sample)
        seg_list.append(ori_image[(h * 256):(h * 256 + 256), -256:, :])
    for w in range(w_step - 1):
        seg_list.append(ori_image[-256:, (w * 256):(w * 256 + 256), :])
    seg_list.append(ori_image[-256:, -256:, :])
##############################裁剪結束########################

##############################  2、利用網路模型推理小圖塊  ########################
    progbar_pred = progbar.Progbar(target=len(seg_list), verbose=1)
    with paddle.no_grad():
        for i, im in enumerate(seg_list):
            ori_shape = im.shape[:2]  #原始圖片形狀(h,w)
            im, _ = transforms(im)    #im.shape(3, 256, 256)    _為None
            im = im[np.newaxis, ...]   #im.shape(1,3,256,256)
            im = paddle.to_tensor(im)

            if False:
                pred = infer.aug_inference(
                    model,
                    im,
                    ori_shape=ori_shape,
                    transforms=transforms.transforms,
                    scales=scales,
                    flip_horizontal=flip_horizontal,
                    flip_vertical=flip_vertical,
                    is_slide=is_slide,
                    stride=None,
                    crop_size=None)
            else:
                pred = infer.inference(
                    model,
                    im,
                    ori_shape=ori_shape,
                    transforms=transforms.transforms,
                    is_slide=False,
                    stride=None,
                    crop_size=None)
            pred = paddle.squeeze(pred)  #該OP會洗掉輸入Tensor的Shape中尺寸為1的維度,查看pred的形狀 應該剩下[h,w]
            pred = pred.numpy().astype('uint8')
            predict_list.append(pred)
            progbar_pred.update(i + 1)
##############################推理結束########################

############# 3、將預測后的影像塊再拼接起來 ########################
    count_temp = 0
    tmp = np.ones([ori_image.shape[0], ori_image.shape[1]])
    for h in range(h_step):
        for w in range(w_step):
            tmp[
            h * 256:(h + 1) * 256,
            w * 256:(w + 1) * 256
            ] = predict_list[count_temp]
            count_temp += 1
        tmp[h * 256:(h + 1) * 256, w_rest:] = predict_list[count_temp][:, w_rest:]
        count_temp += 1
    for w in range(w_step - 1):
        tmp[h_rest:, (w * 256):(w * 256 + 256)] = predict_list[count_temp][h_rest:, :]
        count_temp += 1
    tmp[-257:-1, -257:-1] = predict_list[count_temp][:, :]
##################拼接結束########################

    #獲取需要保存的圖片名稱,去掉前面的路徑
    # get the saved name
    if image_dir is not None:
        pass
        #im_file = im_path.replace(image_dir, '')  #例:將PaddleSeg/data/optic_disc_seg/JPEGImages/P0011.jpg替換為/P0011.jpg
    else:
        im_file = os.path.basename(img_lists[local_rank][local_rank])  #帶后綴名
    if im_file[0] == '/':  #去掉/
        im_file = im_file[1:]

#############
    tmp = tmp.astype('uint8')
    # save added image
    added_image = utils.visualize.visualize(args.image_path,tmp, weight=0.6)
    added_image_path = os.path.join(added_saved_dir, im_file)
    mkdir(added_image_path)
    cv2.imwrite(added_image_path, added_image)
    
    # save pseudo color prediction
    pred_mask = utils.visualize.get_pseudo_color_map(tmp)
    pred_saved_path = os.path.join(pred_saved_dir,
                                    im_file.rsplit(".")[0] + ".png")
    mkdir(pred_saved_path)
    pred_mask.save(pred_saved_path)

    # pred_im = utils.visualize(im_path, pred, weight=0.0)
    # pred_saved_path = os.path.join(pred_saved_dir, im_file)
    # mkdir(pred_saved_path)
    # cv2.imwrite(pred_saved_path, pred_im)
 
    #progbar_pred.update(i + 1)

if __name__ == '__main__':
    args = parse_args()
    main(args)

第二個版本的predict.py將“裁剪遙感大圖”和“拼接小圖塊的預測結果”封裝成了函式,分別為CropBigImage(ImagePath,CropScale)PinJie(predict_list , CropScale , ori_image , h_step , w_step , h_rest , w_rest) ,CropBigImage函式可以將遙感大圖裁剪成任意尺寸的小圖塊,第二個版本的predict.py完整代碼如下:

import sys
import argparse
import os
import paddle
from paddleseg.cvlibs import manager, Config
from paddleseg.utils import get_sys_env, logger
import math
import cv2
import numpy as np
from paddleseg import utils
from paddleseg.core import infer
from paddleseg.utils import progbar

def mkdir(path):
    sub_dir = os.path.dirname(path)  #去掉檔案名,回傳目錄
    if not os.path.exists(sub_dir):
        os.makedirs(sub_dir)

def partition_list(arr, m):
    """split the list 'arr' into m pieces"""
    n = int(math.ceil(len(arr) / float(m)))
    return [arr[i:i + n] for i in range(0, len(arr), n)]

def parse_args():
    parser = argparse.ArgumentParser(description='Model prediction')

    # params of prediction
    parser.add_argument(
        "--config", dest="cfg", help="The config file.", default=None, type=str)
    parser.add_argument(
        '--model_path',
        dest='model_path',
        help='The path of model for prediction',
        type=str,
        default=None)
    parser.add_argument(
        '--image_path',
        dest='image_path',
        help=
        'The path of image, it can be a file or a directory including images',
        type=str,
        default=None)
    parser.add_argument(
        '--save_dir',
        dest='save_dir',
        help='The directory for saving the predicted results',
        type=str,
        default='./output/result')

    # augment for prediction
    parser.add_argument(
        '--aug_pred',
        dest='aug_pred',
        help='Whether to use mulit-scales and flip augment for prediction',
        action='store_true')
    parser.add_argument(
        '--scales',
        dest='scales',
        nargs='+',
        help='Scales for augment',
        type=float,
        default=1.0)
    parser.add_argument(
        '--flip_horizontal',
        dest='flip_horizontal',
        help='Whether to use flip horizontally augment',
        action='store_true')
    parser.add_argument(
        '--flip_vertical',
        dest='flip_vertical',
        help='Whether to use flip vertically augment',
        action='store_true')

    # sliding window prediction
    parser.add_argument(
        '--is_slide',
        dest='is_slide',
        help='Whether to prediction by sliding window',
        action='store_true')
    parser.add_argument(
        '--crop_size',
        dest='crop_size',
        nargs=2,
        help=
        'The crop size of sliding window, the first is width and the second is height.',
        type=int,
        default=None)
    parser.add_argument(
        '--stride',
        dest='stride',
        nargs=2,
        help=
        'The stride of sliding window, the first is width and the second is height.',
        type=int,
        default=None)

    return parser.parse_args()


def get_image_list(image_path):
    """Get image list"""
    valid_suffix = [
        '.JPEG', '.jpeg', '.JPG', '.jpg', '.BMP', '.bmp', '.PNG', '.png' ,'.tif'
    ]
    image_list = []
    image_dir = None
    if os.path.isfile(image_path):
        if os.path.splitext(image_path)[-1] in valid_suffix:
            image_list.append(image_path)
    elif os.path.isdir(image_path):
        image_dir = image_path
        for root, dirs, files in os.walk(image_path):  #root=image_path
            for f in files:
                if os.path.splitext(f)[-1] in valid_suffix:
                    image_list.append(os.path.join(root, f))
    else:
        raise FileNotFoundError(
            '`--image_path` is not found. it should be an image file or a directory including images'
        )

    if len(image_list) == 0:
        raise RuntimeError('There are not image file in `--image_path`')

    return image_list, image_dir   #回傳測驗檔案串列

def CropBigImage(ImagePath,CropScale):

    ImagePath = ImagePath
    CropScale = CropScale
    seg_list = []#存盤分割的圖塊
    ori_image=cv2.imread(ImagePath)##
    h_step = ori_image.shape[0] // CropScale
    w_step = ori_image.shape[1] // CropScale

    h_rest = -(ori_image.shape[0] - CropScale * h_step)
    w_rest = -(ori_image.shape[1] - CropScale * w_step)

    # 回圈切圖
    for h in range(h_step):
        for w in range(w_step):
            # 劃窗采樣
            image_sample = ori_image[(h * CropScale):(h * CropScale + CropScale),
                           (w * CropScale):(w * CropScale + CropScale), :]
            seg_list.append(image_sample)
        seg_list.append(ori_image[(h * CropScale):(h * CropScale + CropScale), -CropScale:, :])
    for w in range(w_step - 1):
        seg_list.append(ori_image[-CropScale:, (w * CropScale):(w * CropScale + CropScale), :])
    seg_list.append(ori_image[-CropScale:, -CropScale:, :])

    return seg_list , ori_image , h_step , w_step , h_rest , w_rest

def PinJie(predict_list , CropScale , ori_image , h_step , w_step , h_rest , w_rest):

    # 將預測后的影像塊再拼接起來
    count_temp = 0
    tmp = np.ones([ori_image.shape[0], ori_image.shape[1]])
    for h in range(h_step):
        for w in range(w_step):
            tmp[
            h * CropScale:(h + 1) * CropScale,
            w * CropScale:(w + 1) * CropScale
            ] = predict_list[count_temp]
            count_temp += 1
        tmp[h * CropScale:(h + 1) * CropScale, w_rest:] = predict_list[count_temp][:, w_rest:]
        count_temp += 1
    for w in range(w_step - 1):
        tmp[h_rest:, (w * CropScale):(w * CropScale + CropScale)] = predict_list[count_temp][h_rest:, :]
        count_temp += 1
    tmp[-(CropScale+1):-1, -(CropScale+1):-1] = predict_list[count_temp][:, :]
    return tmp.astype('uint8')

def main(args):
    env_info = get_sys_env()
    place = 'gpu' if env_info['Paddle compiled with cuda'] and env_info[
        'GPUs used'] else 'cpu'

    paddle.set_device(place)
    if not args.cfg:
        raise RuntimeError('No configuration file specified.')

    cfg = Config(args.cfg)
    val_dataset = cfg.val_dataset   #用val_dataset?
    if not val_dataset:
        raise RuntimeError(
            'The verification dataset is not specified in the configuration file.'
        )

    msg = '\n---------------Config Information---------------\n'
    msg += str(cfg)
    msg += '------------------------------------------------'
    logger.info(msg)

    model = cfg.model
    transforms = val_dataset.transforms
    #image_list, image_dir = get_image_list('data/UAV_seg/images')   
    image_list, image_dir = get_image_list(args.image_path)#需要傳入args.image_path引數   這個引數可以是測驗圖片的路徑,也可以是單張圖片的路徑
    
    model_path=args.model_path  #傳入訓練模型的路徑
    save_dir=args.save_dir
    aug_pred=args.aug_pred
    scales=args.scales
    flip_horizontal=args.flip_horizontal
    flip_vertical=args.flip_vertical
    is_slide=args.is_slide
    crop_size=args.crop_size
    stride=args.stride

    para_state_dict = paddle.load(model_path)
    model.set_dict(para_state_dict)
    model.eval()
    nranks = paddle.distributed.get_world_size()
    local_rank = paddle.distributed.get_rank()
    if nranks > 1:
        img_lists = partition_list(image_list, nranks)
    else:
        img_lists = [image_list]  #是串列還是串列的串列,等待測驗   img_lists[0]  ->串列的串列

    added_saved_dir = os.path.join(save_dir, 'added_prediction')  #偽彩色和原圖疊加
    pred_saved_dir = os.path.join(save_dir, 'pseudo_color_prediction')  #偽彩色預測結果
    
    #主要將遙感大圖裁剪成固定尺寸的圖塊,生成圖塊串列
    ImagePath = img_lists[local_rank][local_rank]
    CropScale = 256
    seg_list , ori_image , h_step , w_step , h_rest , w_rest = CropBigImage(ImagePath,CropScale)

    predict_list = []
    progbar_pred = progbar.Progbar(target=len(seg_list), verbose=1)
    logger.info("Start to predict...")
    with paddle.no_grad():
        for i, im in enumerate(seg_list):
            ori_shape = im.shape[:2]  #原始圖片形狀(h,w)
            im, _ = transforms(im)    #im.shape(3, 512, 512)    _為None
            im = im[np.newaxis, ...]   #im.shape(1,3,512,512)
            im = paddle.to_tensor(im)

            if aug_pred:
                pred = infer.aug_inference(
                    model,
                    im,
                    ori_shape=ori_shape,
                    transforms=transforms.transforms,
                    scales=scales,
                    flip_horizontal=flip_horizontal,
                    flip_vertical=flip_vertical,
                    is_slide=is_slide,
                    stride=stride,
                    crop_size=crop_size)
            else:
                pred = infer.inference(
                    model,
                    im,
                    ori_shape=ori_shape,
                    transforms=transforms.transforms,
                    is_slide=is_slide,
                    stride=stride,
                    crop_size=crop_size)
            pred = paddle.squeeze(pred)  #該OP會洗掉輸入Tensor的Shape中尺寸為1的維度,查看pred的形狀 應該剩下[h,w]
            pred = pred.numpy().astype('uint8')
            predict_list.append(pred)
            progbar_pred.update(i + 1)

    #主要將圖塊的預測結果拼接成大圖
    tmp = PinJie(predict_list , CropScale , ori_image , h_step , w_step , h_rest , w_rest)
    #############
    #獲取需要保存的圖片名稱,去掉前面的路徑
    # get the saved name
    if image_dir is not None:
        pass
        #im_file = im_path.replace(image_dir, '')  #例:將PaddleSeg/data/optic_disc_seg/JPEGImages/P0011.jpg替換為/P0011.jpg
    else:
        im_file = os.path.basename(img_lists[local_rank][local_rank])  #帶后綴名
    if im_file[0] == '/':  #去掉/
        im_file = im_file[1:]

    # save added image
    added_image = utils.visualize.visualize(args.image_path,tmp, weight=0.6)
    added_image_path = os.path.join(added_saved_dir, im_file)
    mkdir(added_image_path)
    cv2.imwrite(added_image_path, added_image)
    
    # save pseudo color prediction
    pred_mask = utils.visualize.get_pseudo_color_map(tmp)
    pred_saved_path = os.path.join(pred_saved_dir,
                                    im_file.rsplit(".")[0] + ".png")
    mkdir(pred_saved_path)
    pred_mask.save(pred_saved_path)
    logger.info("-"*30+"END"+"-"*30)

if __name__ == '__main__':
    args = parse_args()
    main(args)

改寫了predict.py原始碼檔案,就要測驗下它的效果,本文用了一張無人機遙感影像,目的是作物分類,如下圖:
在這里插入圖片描述
為了減輕邊緣效應和拼接痕跡,這里使用重疊度為50%的裁剪方式將原圖裁剪成7000多張256x256的資料集,利用Unet網路對資料集進行訓練,利用本文改寫的predict.py對原圖進行預測,運行predict.py代碼參考如下:!python predict.py --config unet-uav.yml --model_path output/best_model/model.pdparams --image_path /home/aistudio/data/data70483/img.png,其中 –config–model_path–image_path 都是需要傳入的引數,有這些引數但不僅限這些引數,預測結果圖如下:
在這里插入圖片描述
在這里插入圖片描述
從語意分割的結果來看,感徑訓不錯,不過這里我的訓練集和測驗集是同一個資料集,所以并不能說明網路模型的泛化能力,只能說明網路模型的擬合能力還可以,但是本文目的已經達到了,就是對遙感影像(大圖)預測,

各位小伙伴有任何問題可以在評論中留言,下一篇博文的內容依然是使用paddleseg套件對遙感影像預測,不過下篇博文的方法和以上代碼有所差別,主要是做有重疊度裁剪待預測遙感大圖和忽略相鄰圖塊重疊部分做拼接,目的是為了減輕邊緣效應和拼接痕跡,這對語意分割來說十分重要,

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/258732.html

標籤:python

上一篇:【paddle領航團基礎python課程】三歲水課—結營大作業

下一篇:【人生苦短,我學 Python】進階篇——函式與模塊(Day13)

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more