
function yNew = qcWG_NaNcb(src,event,x,y,ha)

axes(ha);
gx = get(ha,'XLim');
gy = get(ha,'YLim');
toNaN=x>=gx(1) & x<=gx(2) & y>=gy(1) & y<=gy(2);
f = find(toNaN);
plot(x(f),y(f),'r.');

ha.UserData(:,end+1)=toNaN;
