R8206/myCrop.py
81261d3c4989master

myCrop.py
View Options

			import numpy as np

			def myCrop(img,inds,fill=0):
			outsize=[]
			dstinds=[]
			srcinds=[]
			for k in range(len(inds)):
			srcbeg=max(0,inds[k].start)
			srcend=min(img.shape[k],inds[k].stop)
			dstbeg=max(0,-inds[k].start)
			dstend=dstbeg+srcend-srcbeg
			dstinds.append(slice(dstbeg,dstend))
			srcinds.append(slice(srcbeg,srcend))
			outsize.append(inds[k].stop-inds[k].start)
			crop=np.full(tuple(outsize),fill,img.dtype)
			crop[tuple(dstinds)]=img[tuple(srcinds)]
			return crop, dstinds