% bzz=1500
xg=max(max(G)); G=log(G+xg/bzz); xg=log(xg+xg/bzz);
ng=min(min(G)); cg=255/(xg-ng); image(cg*(G-ng));
