Aquila-Consortium/map2map
32
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Replace np.take by advanced indexing for cropping

Browse Source
This commit is contained in:
Yin Li 2020-07-14 13:36:09 -04:00
parent c61b6d6fab
commit c2edeb344e
1 changed files with 17 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
map2map/data/fields.py
View File

@ -172,14 +172,14 @@ class FieldDataset(Dataset):
if shift is not None:
anchor[d] += torch.randint(shift, (1,))
in_fields = crop(in_fields, anchor, self.crop, self.pad)
tgt_fields = crop(tgt_fields, anchor * self.scale_factor,
self.crop * self.scale_factor,
np.zeros_like(self.pad))
in_fields = [torch.from_numpy(f).to(torch.float32) for f in in_fields]
tgt_fields = [torch.from_numpy(f).to(torch.float32) for f in tgt_fields]
in_fields = crop(in_fields, anchor, self.crop, self.pad, self.size)
tgt_fields = crop(tgt_fields, anchor * self.scale_factor,
self.crop * self.scale_factor,
np.zeros_like(self.pad), self.size)
if self.in_norms is not None:
for norm, x in zip(self.in_norms, in_fields):
norm(x)
@ -208,12 +208,20 @@ class FieldDataset(Dataset):
return in_fields, tgt_fields
def crop(fields, anchor, crop, pad):
def crop(fields, anchor, crop, pad, size):
ndim = len(size)
assert all(len(x) == ndim for x in [anchor, crop, pad, size]), 'inconsistent ndim'
new_fields = []
for x in fields:
for d, (a, c, (p0, p1)) in enumerate(zip(anchor, crop, pad)):
begin, end = a - p0, a + c + p1
x = x.take(range(begin, end), axis=1 + d, mode='wrap')
ind = [slice(None)]
for d, (a, c, (p0, p1), s) in enumerate(zip(anchor, crop, pad, size)):
i = torch.arange(a - p0, a + c + p1)
i %= s
i = i.reshape((-1,) + (1,) * (ndim - d - 1))
ind.append(i)
x = x[ind]
new_fields.append(x)