c2 = magnitude[1:, 1:][pts[:-1, :-1]]