#DRCMoutfile = '../testimgdir/DRCMsubdimg.tif'   # dark row column average subtraction
#CMoutfile = '../testimgdir/CM10xsubdimg.tif'       # column average subtraction
#CSMoutfile = '../testimgdir/CM10xsubdimg.tif'       # column sigma range average subtraction

# load
print 'loading image...'
img = imgutil.loadimg(datpic)
#fullimg = imgutil.loadimg(datpic,'full')

# do subtractions
print 'doing subtraction...'
#DRCMimg = submethods.darkcolsub(fullimg)
#CMimg = submethods.colmeansub(img)
#CSMimg = submethods.colsigsub(img)
ix,iy = 1080,1080-10
iw,ih = 4,5
window,(x,y),(cx,cy) = submethods.windowsub(img,(ix,iy),(iw,ih))


# display images
#imgutil.dispimg(DRCMimg)
print 'displaying image...'
#imgutil.dispimg(CMimg,10) # displayed with an optional view factor of 10
imgutil.dispimg(window)
#imgutil.dispimg(img[100:100+int(5*2.5)+1][:,100:100+int(4*2.5)+1])

# save images
#imgutil.saveimg(DRCMimg,DRCMoutfile)
#imgutil.saveimg(CMimg,CMoutfile,10) # saved with an optional view factor of 10

# image path setup
datpic = '../testimgdir/img_1348355543_717188_00015_00000_0.dat'
#DRCMoutfile = '../testimgdir/DRCMsubdimg.tif'   # dark row column average subtraction
#CMoutfile = '../testimgdir/CM10xsubdimg.tif'       # column average subtraction
#CSMoutfile = '../testimgdir/CM10xsubdimg.tif'       # column sigma range average subtraction

# load
print 'loading image...'
img = imgutil.loadimg(datpic)
#fullimg = imgutil.loadimg(datpic,'full')

# do subtractions
print 'doing subtraction...'
#DRCMimg = submethods.darkcolsub(fullimg)
#CMimg = submethods.colmeansub(img)
#CSMimg = submethods.colsigsub(img)
ix, iy = 1080, 1080 - 10
iw, ih = 4, 5
window, (x, y), (cx, cy) = submethods.windowsub(img, (ix, iy), (iw, ih))

# display images
#imgutil.dispimg(DRCMimg)
print 'displaying image...'
#imgutil.dispimg(CMimg,10) # displayed with an optional view factor of 10
imgutil.dispimg(window)
#imgutil.dispimg(img[100:100+int(5*2.5)+1][:,100:100+int(4*2.5)+1])

# save images
#imgutil.saveimg(DRCMimg,DRCMoutfile)
#imgutil.saveimg(CMimg,CMoutfile,10) # saved with an optional view factor of 10
# find good star
centers = centroid.findstars(img)
print centers
goodstar = []
goodstar.append(centers[10])
#goodstar = list((tuple(centers[10])) # good centroid at about 1417,275

(x, y) = goodstar[0][0]
(w, h) = goodstar[0][1]

centroid_iwc = centroid.imgcentroid(img, goodstar, method="iwc")
centroid_cog = centroid.imgcentroid(img, goodstar, method="cog")
centroid_gauss = centroid.imgcentroid(img, goodstar, method="gauss")

(frame, (xframe, yframe), frame_centroid) = sm.windowsub(img, (x, y), (w, h),
                                                         neg=True,
                                                         scale=1)

cents = []
cents.append((centroid_iwc[0][0] - xframe, centroid_iwc[0][1] - yframe))
cents.append((centroid_cog[0][0] - xframe, centroid_cog[0][1] - yframe))
cents.append((centroid_gauss[0][0] - xframe, centroid_gauss[0][1] - yframe))
print cents

pl.figure()
pl.gray()
pl.imshow(frame,
          cmap=None,
          norm=None,
          aspect=None,
          interpolation='nearest',
img = imgutil.loadimg(name,from_database=True)
# find good star
centers = centroid.findstars(img)
print centers
goodstar = []
goodstar.append(centers[10])
#goodstar = list((tuple(centers[10])) # good centroid at about 1417,275

(x,y) = goodstar[0][0]
(w,h) = goodstar[0][1]

centroid_iwc   = centroid.imgcentroid(img, goodstar, method="iwc")
centroid_cog   = centroid.imgcentroid(img, goodstar, method="cog")
centroid_gauss = centroid.imgcentroid(img, goodstar, method="gauss")
        
(frame, (xframe,yframe), frame_centroid) = sm.windowsub(img, (x,y), (w,h), neg=True, scale=1)

cents = []
cents.append((centroid_iwc[0][0]-xframe,centroid_iwc[0][1]-yframe)) 
cents.append((centroid_cog[0][0]-xframe,centroid_cog[0][1]-yframe)) 
cents.append((centroid_gauss[0][0]-xframe,centroid_gauss[0][1]-yframe)) 
print cents

pl.figure()
pl.gray()
pl.imshow(frame, cmap=None, norm=None, aspect=None,
            interpolation='nearest', vmin=0, vmax=2048, origin='upper')

radius = 1 
#iwc            
color = 'r'