给定一个矩形图像img
和补丁s
。现在我想用边长的正方形补丁覆盖整个图像s
,以便img
使用最少数量的补丁,其中的每个像素都在至少一个补丁中。此外,我希望相邻的补丁尽可能少地重叠。
到目前为止:我在下面包含了我的代码并制定了一个示例。但是它还不能完美地工作。希望有人发现错误。
示例:Given is img
of shape:
这意味着我将在较长的一侧放置 21 个补丁,在宽度较小的一侧放置 14 个补丁,总共 21*14 = 294 个补丁(4616, 3016)
。s = 224
现在我试图弄清楚补丁如何分配补丁之间的重叠。我的补丁可以覆盖 size: 的图像(4704, 3136)
,因此我的高度必须覆盖 88 个重叠像素missing_h = ht * s - h
,宽度类似。
现在我试着弄清楚,如何在 21 个补丁上分配 88 个像素。88 = 4* 21 + 4 因此,我将有hso = 17
重叠 shso = 4
的hbo = 4
补丁和重叠 5 的补丁,宽度是类似的。
现在我只需遍历整个图像并跟踪我当前的位置(cur_h, cur_w)
。在我调整每个循环之后,cur_h, cur_w
. 我有s
,我当前的补丁号i, j
,它表明补丁是否有小或大的重叠。
import numpy as np
def part2(img, s):
h = len(img)
w = len(img[0])
ht = int(np.ceil(h / s))
wt = int(np.ceil(w / s))
missing_h = ht * s - h
missing_w = wt * s - w
hbo = missing_h % ht
wbo = missing_w % wt
hso = ht - hbo
wso = wt - wbo
shso = int(missing_h / ht)
swso = int(missing_w / wt)
patches = list()
cur_h = 0
for i in range(ht):
cur_w = 0
for j in range(wt):
patches.append(img[cur_h:cur_h + s, cur_w: cur_w + s])
cur_w = cur_w + s
if j < wbo:
cur_w = cur_w - swso - 1
else:
cur_w = cur_w - swso
cur_h = cur_h + s
if i < hbo:
cur_h = cur_h - shso - 1
else:
cur_h = cur_h - shso
if cur_h != h or cur_w != w:
print("expected (height, width)" + str((h, w)) + ", but got: " + str((cur_h, cur_w)))
if wt*ht != len(patches):
print("Expected number patches: " + str(wt*ht) + "but got: " + str(len(patches)) )
for patch in patches:
if patch.shape[0] != patch.shape[1] or patch.shape[0] != s:
print("expected shape " + str((s, s)) + ", but got: " + str(patch.shape))
return patches
def test1():
img = np.arange(0, 34 * 7).reshape((34, 7))
p = part2(img, 3)
print("Test1 successful")
def test2():
img = np.arange(0, 4616 * 3016).reshape((4616, 3016))
p = part2(img, 224)
print("Test2 successful")
test1()
test2()