function [It, Ix, Iy] = py_partial_deriv(I1,I2,wl,deriv_filter)
    addpath('code/external/utils');
    [It,Ix,Iy] = partial_deriv(cat(3,I1,I2), wl, 'bi-cubic', deriv_filter);
end