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

Merge branch 'master' of github.com:eelregit/map2map

Browse Source
This commit is contained in:
Yin Li 2020-01-23 11:35:57 -06:00
commit 1b3f32b791
1 changed files with 7 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
map2map/data/fields.py
View File

@ -149,23 +149,22 @@ def crop(fields, start, crop, pad, scale_factor=1):
new_fields = []
for x in fields:
for d, (i, N, (p0, p1)) in enumerate(zip(start, crop, pad)):
begin, stop = i - p0, i + N + p1
begin, end = i - p0, i + N + p1
if scale_factor > 1: # add buffer for linear interpolation
begin, stop = begin - 1, stop + 1
begin, end = begin - 1, end + 1
x = x.take(range(begin, stop), axis=1 + d, mode='wrap')
x = x.take(range(begin, end), axis=1 + d, mode='wrap')
if scale_factor > 1:
x = np.expand_dims(x,axis=0)
x = torch.from_numpy(x)
x = torch.from_numpy(x).unsqueeze(0)
x = F.interpolate(x, scale_factor=scale_factor, mode='trilinear')
x = x[0].numpy()
x = x.numpy().squeeze(0)
# remove buffer
for d, (N, (p0, p1)) in enumerate(zip(crop, pad)):
begin, stop = scale_factor, N + p0 + p1 - scale_factor
x = x.take(range(begin, stop), axis=1 + d)
begin, end = scale_factor, N + p0 + p1 - scale_factor
x = x.take(range(begin, end), axis=1 + d)
new_fields.append(x)