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)



