IDL> .run hot.pix
% Compiled module: $MAIN$.
     
**Reading FITS file: r621196pp.fits
% READFITS: Now reading 1024 by 1024 array
     
**Applying a median filter to the dark frame (width=10 pixels)
     
**Calculation of the stdev in the image difference (dark-smoothed_dark)
  -creating 256.000 images of 64.0000 pixels
  -calculating stdev in each subimage after a 3-sigma clipping around median
  -calculating stdev over the whole image with  3-sigma clipping
   around median on the stdev measurements of the individual sub-images
     
      ==>  stdev in whole image=5.72398 ADU
     
      ==>  Per quadrants:
            left-up    =6.00398 ADU
            right-up   =5.88131 ADU
            right-down =5.53623 ADU
            left-down  =4.89111 ADU
     
**Creating a hot pixel mask
  -considering hot pixels those in which signal>4*stdev in dark-smoothed_dark
  -Number of hot pixels detected: 21779
  -Percentage of hot pixels in the whole image: 2.07701%
  -writing mask into file: hpm_r621196pp.fits
     
** Number of hot pixels and percentage per quadrant:
     
  -Quadrant left-up: 3701pix, 1.41182%
  -Quadrant right-up: 3651pix, 1.39275%
  -Quadrant right-down: 10772pix, 4.10919%
  -Quadrant left-down: 3655pix, 1.39427%
     
**Dark level in each quadrant:
 
  -Quadrant left-up: 4.00000ADU 
  -Quadrant right-up: 3.00000ADU 
  -Quadrant right-down: 8.00000ADU 
  -Quadrant left-down: 8.00000ADU 
IDL>