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

          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: