Uploaded image for project: 'Instrument control development'
  1. Instrument control development
  2. INSTRM-2558

Modify criteria of saturation flag

    XMLWordPrintable

    Details

    • Sprint:
      PreRun23Jun2025

      Description

      From Slack (auto-guiding channel):

      The value at which the CCD bleeds isn't totally stable – it's not like the 65535 which you get from a value larger than the maximum value you can store in an unsigned 16-bit value [from a 16-bit A/D converter]

       

      # determine saturation off the unprocessed data
          satFlag = data[result['peak_pixel_y_pix'],result['peak_pixel_x_pix']]==satValue
          result['flags'] += satFlag*4
      

      E.g. here's the centre of a saturated star in pfs_visit_id, agc_exposure_id = 125713, 699879; AG2 

      [[ 3632.,  3738.,  3820.,  3646.,  3429.,  3165.,  2996.,  2653.],{{       }}
      [ 5054.,  5299.,  5217.,  5032.,  4582.,  4245.,  3829.,  3323.],{{       }}
      [ 7811.,  8110.,  7855.,  7411.,  6788.,  5878.,  4967.,  4116.],       [13237., 29523., 29500., 12164., 10477.,  8915.,  7183.,  5435.],       [40880., 59682., 59352., 52460., 17490., 14037., 10655.,  7630.],       [59490., 59730., 59750., 59638., 56759., 22897., 16250., 10798.],       [59700., 59729., 59748., 59769., 59425., 38507., 24197., 14974.],       [59625., 59711., 59788., 59709., 59799., 59059., 34645., 20425.],       [59588., 59764., 59808., 59826., 59884., 59644., 46327., 26123.],       [59755., 59828., 59815., 59904., 59866., 59838., 54580., 30783.],       [59652., 59800., 59832., 59882., 59826., 59851., 57338., 32382.],       [59565., 59698., 59798., 59796., 59823., 59669., 52013., 29442.],       [59455., 59730., 59738., 59840., 59806., 59368., 40194., 23795.],       [36483., 59665., 59743., 59818., 59735., 58673., 27015., 17022.],       [16849., 38153., 59521., 59688., 59093., 22821., 16808., 11501.],       [10650., 13025., 17129., 38198., 30490., 12851., 10213.,  7708.],

      The right side should be like

      data[result['peak_pixel_y_pix'],result['peak_pixel_x_pix']] >= satVal

      And satVal should be different among AG camera, channel. So there should be 12 values.

       

        Attachments

        1. max.png
          61 kB
          karr
        2. v127208 AG Cam2.png
          105 kB
          Wilfred Gee

          Issue Links

            Activity

              People

              • Assignee:
                karr karr
                Reporter:
                yuki.moritani yuki.moritani
              • Votes:
                0 Vote for this issue
                Watchers:
                4 Start watching this issue

                Dates

                • Created:
                  Updated:
                  Resolved: