compute the grain boundary plane distribution
GBPD = calcGBND(traces,ori)
% use a specific halfwidth GBPD = calcGBPD(traces,ori,'halfwidth',10*degree)