Uploaded image for project: 'Data Management'
  1. Data Management
  2. DM-37357

Update masking in parallel overscan

    XMLWordPrintable

    Details

    • Type: Improvement
    • Status: Done
    • Resolution: Done
    • Fix Version/s: None
    • Component/s: None
    • Labels:
      None
    • Story Points:
      8
    • Epic Link:
    • Team:
      Data Release Production

      Description

      Parallel overscan masking is still preventing potentially correctable amplifiers from being processed.  Rework the masking to check per column (or per group of columns).

        Attachments

          Issue Links

            Activity

            Hide
            czw Christopher Waters added a comment -

            It will need to be next week, but I'll make sure the MEDIAN_PER_ROW is getting the right thing (that's what I was testing with).  I'll also add additional debug plots that demonstrate the results.

            Show
            czw Christopher Waters added a comment - It will need to be next week, but I'll make sure the MEDIAN_PER_ROW is getting the right thing (that's what I was testing with).  I'll also add additional debug plots that demonstrate the results.
            Hide
            czw Christopher Waters added a comment -

            I've uploaded an update: Craig Lage was correct that I had incorrectly put the masked-value "interpolation" fix in a function that wasn't being called by the MEDIAN_PER_ROW correction (but confusingly was run by the debug code).  This resulted in the columns with bad bleeds continuing to have NaN values, causing the masking code to mask and interpolate these pixels (doing a poor job as there is also oversubtraction that pushed the interpolated values below the background level).  I've fixed that, and although it's better (the crosstalk ghosts do not appear to have significant oversubtraction), the source bleed is still oversubtracting (~65 ADU on the shown image).  Growing the input mask helps somewhat, with 2-pixel growth reducing the oversubtraction to ~25 ADU.  I'm going to look at some bleeds, and see if there's a larger grow value that yields better results.

            I've pushed the current version in case others want to look at it, and I'm hoping to have a cleaned up version ready tomorrow that can hopefully go into review.

            Show
            czw Christopher Waters added a comment - I've uploaded an update: Craig Lage  was correct that I had incorrectly put the masked-value "interpolation" fix in a function that wasn't being called by the MEDIAN_PER_ROW correction (but confusingly was run by the debug code).  This resulted in the columns with bad bleeds continuing to have NaN values, causing the masking code to mask and interpolate these pixels (doing a poor job as there is also oversubtraction that pushed the interpolated values below the background level).  I've fixed that, and although it's better (the crosstalk ghosts do not appear to have significant oversubtraction), the source bleed is still oversubtracting (~65 ADU on the shown image).  Growing the input mask helps somewhat, with 2-pixel growth reducing the oversubtraction to ~25 ADU.  I'm going to look at some bleeds, and see if there's a larger grow value that yields better results. I've pushed the current version in case others want to look at it, and I'm hoping to have a cleaned up version ready tomorrow that can hopefully go into review.
            Hide
            cslage Craig Lage added a comment - - edited

            I looked in more detail at that dark column in 2022121200752. The attached image Dark_Column... shows plots through the various regions. I concluded that the negative values in the PostISR image are not caused by a failure of the parallel overscan, but by the negative values that are already present in the raw values in the imaging region. I'm not sure what causes those.

            Ignore this. I once again fell into the trap that the exposures get modified in place. The attached Dark_Column... is corrected and it doesn't really show anything we didn't know. It looks like we need to go +/- 3-4 columns to mask out that hot column.

            Show
            cslage Craig Lage added a comment - - edited I looked in more detail at that dark column in 2022121200752. The attached image Dark_Column... shows plots through the various regions. I concluded that the negative values in the PostISR image are not caused by a failure of the parallel overscan, but by the negative values that are already present in the raw values in the imaging region. I'm not sure what causes those. Ignore this. I once again fell into the trap that the exposures get modified in place. The attached Dark_Column... is corrected and it doesn't really show anything we didn't know. It looks like we need to go +/- 3-4 columns to mask out that hot column.
            Hide
            czw Christopher Waters added a comment -

            I played with the grow parameter last night, and eventually settled on a 7 pixel growth radius.  This still isn't perfect, but the oversubtraction is down to ~5ADU, which would be lost in the noise if the low pixels weren't correlated.

            I'm going to try to get this into review, and have launched a Jenkins run to ensure there aren't any hidden errors.

            https://ci.lsst.codes/blue/organizations/jenkins/stack-os-matrix/detail/stack-os-matrix/37949/pipeline

            Show
            czw Christopher Waters added a comment - I played with the grow parameter last night, and eventually settled on a 7 pixel growth radius.  This still isn't perfect, but the oversubtraction is down to ~5ADU, which would be lost in the noise if the low pixels weren't correlated. I'm going to try to get this into review, and have launched a Jenkins run to ensure there aren't any hidden errors. https://ci.lsst.codes/blue/organizations/jenkins/stack-os-matrix/detail/stack-os-matrix/37949/pipeline
            Hide
            erykoff Eli Rykoff added a comment -

            Looks good pending final jenkins with ci_cpp.

            Show
            erykoff Eli Rykoff added a comment - Looks good pending final jenkins with ci_cpp.

              People

              Assignee:
              czw Christopher Waters
              Reporter:
              czw Christopher Waters
              Reviewers:
              Eli Rykoff
              Watchers:
              Aaron Roodman, Christopher Waters, Craig Lage, Eli Rykoff, Stuart Marshall
              Votes:
              0 Vote for this issue
              Watchers:
              5 Start watching this issue

                Dates

                Created:
                Updated:
                Resolved:

                  Jenkins

                  No builds found.