# -*- coding: utf-8 -*-

from pylab import *
def sg(a):
    figure(); imshow(a,cmap='gray',vmin=0,vmax=1)
    
org=imread('../Images/libwork.png')
# org=imread("../Images/osubib.png")
sh=org.shape

# Conversion CMYK (Wpd, alg.1)
a=org
gr=a[:,:,1]
rd=a[:,:,0]

cp=1-a[:,:,0]
mp=1-a[:,:,1]
yp=1-a[:,:,2]
k=minimum(cp,minimum(mp,yp))

c=cp-k;  m=mp-k;  y=yp-k

"""
figure(); imshow(k,cmap='gray')
figure(); imshow(c,cmap='gray')
figure(); imshow(m,cmap='gray')
figure(); imshow(y,cmap='gray')
"""

# manip

c=rd*0.55
# Reconstruction (inverse)

rr=1-(c+k); gg=1-(m+k);  bb=1-(y+k)
nw=zeros(sh)
nw[:,:,0]=rr; nw[:,:,1]=gg;  nw[:,:,2]=bb;

# ======================

# figure(); imshow(nw,vmin=0,vmax=1)
imsave("../Images/pylcorr.png",nw)
# sg(k)