def compute1didx(extar,slc):
x1=np.argmin(np.abs(extar[0]-slc[0]))
x2=np.argmin(np.abs(extar[0]-slc[1]))
if len(extar) == 2:
y1=np.argmin(np.abs(extar[1]-slc[2]))
y2=np.argmin(np.abs(extar[1]-slc[3]))
if x1==x2:
IDX=np.s_[x1,y1:y2]
elif y1==y2:
IDX=np.s_[x1:x2,y1]
if len(extar) == 3:
z1=np.argmin(np.abs(extar[2]-slc[4]))
z2=np.argmin(np.abs(extar[2]-slc[5]))
if (x1==x2 and y1==y2): IDX=np.s_[x1,y1,z1:z2]
if (y1==y2 and z1==z2): IDX=np.s_[x1:x2,y1,z1]
if (x1==x2 and z1==z2): IDX=np.s_[x1,y1:y2,z1]
else:
IDX=np.s_[x1:x2]
return IDX
评论列表
文章目录