I found that this works for me. Suggestions are more than welcome. Thanks
1)self define a colormap, e.g
%Yellow [1 1 0]
YellowColormap=zeros(256,3);
for i=1:256
YellowColormap(i,1) = (i-1)/256;
YellowColormap(i,2) = (i-1)/256;
end
2)X=imread(16bit image);
imshow(X,'DisplayRange',[0 65535],'Colormap',YellowColorMap)