function y = impowlaw(x, gamma) c = 255/255^(gamma); y = uint8(c*double(x).^gamma);