Как рассчитать масштабно-инвариантную энергию изгиба для бинарного листа в Matlab?

Question

Я хочу рассчитать энергию изгиба, как описано здесь: https://nph.onlinelibrary.wiley.com/doi/pdf/10.1111/j.1469-8137.2010.03266.x

Мой код Matlab такой:

[B,L,N,A] = bwboundaries(BW);


for k0 = 1:length(B)
    boundary = B{k0};
    bound2 = boundary(:,2);
    bound1 = boundary(:,1);

    arith = abs(gradient(bound2).*gradient(gradient(bound1)) - gradient(bound1).*gradient(gradient(bound2)));
    paron = (gradient(bound2).^2 + gradient(bound1).^2).^(3/2);
    Kappa = arith./paron;

    P1 = regionprops(BW,'perimeter');
    BW_per = P1.Perimeter;

    BE=nanmean(Kappa)
end

Моя проблема в том, что я хочу сделать его масштабно-инвариантным, чтобы рассчитать его для разных масштабов листовых изображений. Здесь описывается путь, но я не уверен, как мне достичь этого в Matlab:

Есть идеи? Заранее спасибо!