我想要“修改”Mathematica的Interpolation[]函数(在一维中),当输入超出范围时,用常量值替换外推。换句话说,如果插值域是1,20,f1==7和f20==12,我想:f[x] = 12 for x>=20 (* interpolation w cutoff *)
(* sort array(* dit
.`), numpy.newaxis (`None`) and integer or boolean arrays are valid indices scidata0 = fits.getdata(data,0)
scidata0[0,0,:,:] = ndimage.interpolation.shift