#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'