function [icx,ice,ich,icr,irx,ire,irh,irr] = procim(image,Kc,Kr)
% [icx,ice,ich,icr,irx,ire,irh,irr] = procim(image,Kc,Kr)
%    Does complete processing of gene image
[m n] = size(image);
doing_cols_cba = 0
[icc,icb,ica] = qf2dcol(image,Kc);
doing_cols_xehr = 1
[icx,ice,ich,icr] = q2par(image,Kc,icc,icb,ica);
imaget = image';
doing_rows_bca = 2
[irct,irbt,irat] = qf2dcol(imaget,Kr);
doing_rows_xehr = 3
[irxt,iret,irht,irrt] = q2par(imaget,Kr,irct,irbt,irat);
irc = irct';
irb = irbt';
ira = irat';
irx = irxt';
ire = iret';
irh = irht';
irr = irrt';
%doing_pass1 = 4
%p1 = pass1(irx,icx,Kc,Kr);
%doing_pass2 = 5
%p2 = pass2(p1,irx,icx,Kc,Kr);
%doing_pass3 = 6
%p3 = pass3(p2,irr,icr);
