Solutions to NoData in Remote Sensing Images

Before processing

Insert picture description here

After treatment

Insert picture description here

Code

def NoData_kill(in_path, out_path):
    im_proj, im_geotrans, im_width, im_height, im_data = gdalTools.read_img(in_path)
    mask = np.isnan(im_data)
    c, w, h = mask.shape
    mask_list = []
    for i in range(c):
        if mask[i].__contains__(True):
            mask_list.append(mask[i])

    for m in mask_list:
        m = m + 0
        m = np.uint8(m)
        inpainted_img = cv2.inpaint(im_data, m, inpaintRadius=3, flags=cv2.INPAINT_TELEA)
        im_data = inpainted_img

    gdalTools.write_img(out_path, im_proj, im_geotrans, im_data)

Guess you like

Origin blog.csdn.net/weixin_42990464/article/details/113307642