def METHOD_NAME(im1, im2, p1: ndarray, p2: ndarray) -> None:
    h1, w1, c = im1.shape
    h2, w2, c = im2.shape
    image = np.zeros((max(h1, h2), w1 + w2, 3), dtype=im1.dtype)
    image[0:h1, 0:w1, :] = im1
    image[0:h2, w1 : (w1 + w2), :] = im2
    p1 = features.denormalized_image_coordinates(p1, w1, h1)
    p2 = features.denormalized_image_coordinates(p2, w2, h2)
    pl.imshow(image)
    for a, b in zip(p1, p2):
        pl.plot([a[0], b[0] + w1], [a[1], b[1]], "c")
    pl.plot(p1[:, 0], p1[:, 1], "ob")
    pl.plot(p2[:, 0] + w1, p2[:, 1], "ob")