rMin=minR,
rMax=fovHeight / 2.)
# %% add gradients as a function of distance to border
aryImaOtw = np.zeros((aryImaNtr.shape))
aryImaInw = np.zeros((aryImaNtr.shape))
for indIma in range(aryImaNtr.shape[-1]):
# get image
inputIma = aryImaNtr[..., indIma]
# add gradients as a function of distance to border for outward motion
aryImaOtw[..., indIma] = getDistIma(inputIma,
fovHeight=fovHeight,
pix=pix,
pixSlct="stimAprt",
borderType="right",
normalize=True)
# add gradients as a function of distance to border for inward motion
aryImaInw[..., indIma] = getDistIma(inputIma,
fovHeight=fovHeight,
pix=pix,
pixSlct="stimAprt",
borderType="left",
normalize=True)
for indIma in range(aryImaNtr.shape[-1]):
# get image
inputIma = aryImaNtr[..., indIma]
varFigName = os.path.join("/media/sf_D_DRIVE/MotDepPrf/Presentation",
"figures", "Figure_prfshift",
(cfg.pix, cfg.pix)).reshape(cfg.pix, cfg.pix, 1), opaPgDnMasks),
axis=2)
# %% create masks for the foreground (raised cosine) and group masks for the
# foreground together (raised cosine)
opaPgUpMasks = np.empty((cfg.pix, cfg.pix, opaPgDnMasks.shape[-1]),
dtype='float32')
for i in range(opaPgDnMasks.shape[-1]):
# get a single mask
binMask = opaPgDnMasks[..., i]
# check whether there is at least one element of 1
if np.any(binMask == 1):
print "image " + str(i)
# get its distance image
distIma = getDistIma(binMask, cfg.fovHeight, cfg.pix)
# assign raised cosine values to bixels less than 0.5 away from border
opaPgUpMasks[..., i] = assignBorderVals(binMask,
distIma,
borderRange=cfg.borderRange)
else:
# assign old contrast mask
opaPgUpMasks[..., i] = binMask
# %% save masks as npz array
# restructure arrays such that if indexed with linspace, outward motion results
ouwardInd = np.arange(len(combis) / cfg.numAprtCrcle * cfg.numRep,
dtype='int32').reshape((-1, cfg.numRep)).T
# add 1 to account for the zero image in the begging
ouwardInd = ouwardInd.flatten() + 1