"""

keep track of interesting variables

"""

def __init__(self,l):

self.l=l

def listvar(self,l):

"""

set the list of variables [{"name":name1,"fct":fct1,"txt":descr1},[name2,fct2,descr2],..,[nameN,fctN,descrN]]

"""

self.l=l

def report(self,filename,mode="a",title=None):

f=open(filename,mode)

if title!=None:

f.write("---%s "%title+time.asctime()+"------\n")

else:

f.write("---Stats on "+time.asctime()+"------\n")

for v in self.l:

#if v["name"] in dir()

if v.has_key("fnc"):

exec("value=str(%s(%s))"%(v["fnc"],v["name"]))

else:

if v["name"][-1]=="*":

exec("value=str(%s.__dict__)"%(v["name"][:-1]))

else:

exec("value=str(%s)"%(v["name"]))

#else

# value="Not Exists!"

if v.has_key("txt"):

f.write("%s=%s\n"%(v["txt"],value))

else:

f.write("%s=%s\n"%(v["name"],value))

newseq=[]

for l in seq:

if l!=[]:

newseq.append(l)

bb=numpy.zeros((len(trPosImages)*5,4))#as maximum 5 persons per image in average

name=[]

cnt=0

tot=0

if maxnum==-1:

tot=len(trPosImages)

else:

tot=min(maxnum,trPosImages)

for idx in range(tot):

#print trPosImages.getImageName(idx)

#img=trPosImages.getImage(idx)

rect=trPosImages[idx]["bbox"]#.getBBox(idx,usetr=usetr,usedf=usedf)

for r in rect:

bb[cnt,:]=r[:4]

name.append(trPosImages[idx]["name"])#.getImageName(idx))

cnt+=1

#img=pylab.imread("circle.png")

util.pdone(idx,tot)

ratio=(bb[:,2]-bb[:,0])/(bb[:,3]-bb[:,1])

area=(bb[:,2]-bb[:,0])*(bb[:,3]-bb[:,1])

ftrpos=[]

for iel,el in enumerate(trpos):

ftrpos.append(numpy.zeros(cumsize[-1],dtype=numpy.float32))

ftrpos[-1][cumsize[trposcl[iel]]:cumsize[trposcl[iel]+1]-1]=trpos[iel]

#bias

ftrpos[-1][cumsize[trposcl[iel]+1]-1]=bias

keeps = {}

for name, value in globals().iteritems():

if name in keep: keeps[name] = value

globals().clear()

for name, value in keeps.iteritems():

i=a[0]

imname=a[1]

bbox=a[2]

models=a[3]

cfg=a[4]

if len(a)<=5:

imageflip=False

else:

imageflip=a[5]

img=util.myimread(imname,resize=cfg.resize)

if imageflip:

img=util.myimread(imname,True,resize=cfg.resize)

if bbox!=None:

bbox = util.flipBBox(img,bbox)

if bbox!=None:

gtbbox=[{"bbox":x,"img":imname.split("/")[-1]} for x in bbox]

else:

gtbbox=None

if cfg.show:

img=util.myimread(imname,imageflip,resize=cfg.resize)

pylab.figure(10)

pylab.ioff()

pylab.clf()

pylab.axis("off")

pylab.imshow(img,interpolation="nearest",animated=True)

notsave=False

#if cfg.__dict__.has_key("test"):

# notsave=cfg.test

#f=pyrHOG2.pyrHOG(imname,interv=10,savedir=cfg.auxdir+"/hog/",notsave=not(cfg.savefeat),notload=not(cfg.loadfeat),hallucinate=cfg.hallucinate,cformat=True,flip=imageflip,resize=cfg.resize)

f=pyrHOG2.pyrHOG(img,interv=10,savedir=cfg.auxdir+"/hog/",notsave=not(cfg.savefeat),notload=not(cfg.loadfeat),hallucinate=cfg.hallucinate,cformat=True)#,flip=imageflip,resize=cfg.resize)

res=[]

for clm,m in enumerate(models):

if cfg.useRL:

res.append(pyrHOG2RL.detectflip(f,m,gtbbox,hallucinate=cfg.hallucinate,initr=cfg.initr,ratio=cfg.ratio,deform=cfg.deform,bottomup=cfg.bottomup,usemrf=cfg.usemrf,numneg=cfg.numneg,thr=cfg.thr,posovr=cfg.posovr,minnegincl=cfg.minnegincl,small=cfg.small,show=cfg.show,cl=clm,mythr=cfg.mythr,mpos=cfg.mpos,usefather=cfg.usefather,useprior=cfg.useprior,K=cfg.k))

else:

res.append(pyrHOG2.detect(f,m,gtbbox,hallucinate=cfg.hallucinate,initr=cfg.initr,ratio=cfg.ratio,deform=cfg.deform,bottomup=cfg.bottomup,usemrf=cfg.usemrf,numneg=cfg.numneg,thr=cfg.thr,posovr=cfg.posovr,minnegincl=cfg.minnegincl,small=cfg.small,show=cfg.show,cl=clm,mythr=cfg.mythr,mpos=cfg.mpos,usefather=cfg.usefather,useprior=cfg.useprior,emptybb=False,K=cfg.k))

if cfg.show:

pylab.draw()

pylab.show()

if old==[]:

return new,newcl

unew=[]

unewcl=[]

#mold=numpy.array(old)

clst=[]

for c in range(numcl):

select=numpy.arange(len(oldcl))[numpy.array(oldcl)==c]

clst.append(numpy.array([old[i] for i in select]))

for ep,e in enumerate(new):

#print ep,"/",len(new)

print ".",

#check the cluster

#selec=numpy.arange(len(oldcl))[numpy.array(oldcl)==newcl[ep]]

apr=numpy.sum(numpy.abs(e[::100]-clst[newcl[ep]][:,::100]),1)

#print apr

#raw_input()

if numpy.all(apr>0.1):

unew.append(e)

unewcl.append(newcl[ep])

else:

if numpy.all(numpy.sum(numpy.abs(e-clst[newcl[ep]][apr<=0.1,:]),1)>0.1):

unew.append(e)

unewcl.append(newcl[ep])

print ep,"/",len(new)

print "Doules",len(new)-len(unew)

#raw_input()

ls=[];lscl=[]

for el in dtrpos:

aux=(tr.descr(dtrpos[el],flip=False,usemrf=cfg.usemrf,usefather=cfg.usefather,k=cfg.k))

ls+=aux

auxcl=tr.mixture(dtrpos[el])

lscl+=auxcl

if not(useRL):

ls+=(tr.descr(dtrpos[el],flip=True,usemrf=cfg.usemrf,usefather=cfg.usefather,k=cfg.k))

lscl+=tr.mixture(dtrpos[el])

if cumsize!=None:

dns=buildense(aux,auxcl,cumsize)

#print "Det:",dtrpos[el][0]["img"],[x["scr"] for x in dtrpos[el]],tr.mixture(dtrpos[el])

print "Det:",dtrpos[el][0]["img"],[numpy.sum(x*w)-r for x in dns],tr.mixture(dtrpos[el])

#if cumsize!=None:

# raw_input()

ls=[];lscl=[]

for el in dtrpos:

aux=(tr.descr(dtrpos[el],flip=False,usemrf=cfg.usemrf,usefather=cfg.usefather,k=cfg.k))

#ls+=aux

auxcl=tr.mixture(dtrpos[el])

#lscl+=auxcl

if not(useRL):

aux2=(tr.descr(dtrpos[el],flip=True,usemrf=cfg.usemrf,usefather=cfg.usefather,k=cfg.k))

auxcl2=tr.mixture(dtrpos[el])

for l in range(len(aux)):

ls.append(aux[l])

lscl.append(auxcl[l])

if not(useRL):

ls.append(aux2[l])

lscl.append(auxcl2[l])

if cumsize!=None:

dns=buildense(aux,auxcl,cumsize)

#print "Det:",dtrpos[el][0]["img"],[x["scr"] for x in dtrpos[el]],tr.mixture(dtrpos[el])

print "Det:",dtrpos[el][0]["img"],[numpy.sum(x*w)-r for x in dns],tr.mixture(dtrpos[el])

#if cumsize!=None:

# raw_input()

dns=buildense(trpos,trposcl,cumsize)

loss=0

lscr=[]

for l in dns:

scr=numpy.sum(l*w)

print "Scr:",scr

loss+=max(0,1-scr)

lscr.append(scr)