<p>Take your 2-d image of a star; mask out any pixels that may contain saturated data; fit a Gaussian (or whatever your model for the point spread function is) to the remaining pixels, with the amplitude and source position as free parameters; replace the masked pixels with values computed from the best-fit model you've just calculated.</p>
<p>This will work pretty well as long as the point spread function is not dependent on source intensity and can perhaps be derived from nearby, unsaturated stars.</p>

CC BY-SA 4.0
2023-12-19T00:06:21.977Z
2023-12-19T00:06:21.977Z
ProfRob
<p><strong>Main Question</strong></p>
<p>I have a project goal, where I have an image of light sources (stars or lamps doesn't matter). The image is split up into channels for each expected band of 'color'. There are occasional instances where the sensors become saturated at the highest intensities. I remember hearing at one point that astronomers &quot;Tail Fit&quot; a Gaussian distribution so that it strictly fits the unsaturated data, the saturated data is then ignored.</p>
<p>How are these fits computed?</p>
<p><strong>To put in more specific terms:</strong></p>
<p>You can simulate this by taking a photograph of a lit lamp. Your image will likely result an RGB value per pixel, most home cameras map to color ranges in [0, 255]. When any channel that registers as (255) the sensor is considered saturated. When all channels register as (255, 255, 255) then this is &quot;saturated white&quot;. Saturated regions are unable to provide useful information because true intensity measurements would not have a ceiling.</p>
<p>Models of intensity (per channel) will allow for identification of types of light. In order to make realistic models of intensity we have to correct the saturated data. This would encourage us to fit a curve as informed by only the unsaturated tails.</p>
<p>We want a curve that approximates the saturated section, informed by the unsaturated section. How do we &quot;tail fit&quot; this curve? Is this a strategy that exists or is it a fabrication I have remembered?</p>
<p>I currently expect this to be a Gaussian curve on each channel, however I'm happy to be suggested a different known curve.</p>

CC BY-SA 4.0
2022-01-29T11:56:36.323Z
Just to clarify, this is for a 2-dimensional image? This maybe of use in crafting an answer: [Point-Spread Function Fitting Photometry](https://ui.adsabs.harvard.edu/abs/1999ASPC..189...56H/abstract)
GrapefruitIsAwesome
CC BY-SA 4.0
2022-01-30T17:54:03.813Z
"we have to censor the saturated data" -- maybe you mean "correct the saturated data"? (I fixed the other instances of "censor" -- where you clearly meant *sensor* -- but wasn't sure about that last phrase.)
Peter Erwin
<star><identify-this-object><image-processing><distribution-function><color>
Using tail-fit Gaussian distribution to compensate for saturated sensors
probinso
# Using tail-fit Gaussian distribution to compensate for saturated sensors

(Asked by: username_0 on 2022-01-29)

**Main Question**

I have a project goal, where I have an image of light sources (stars or lamps doesn't matter). The image is split up into channels for each expected band of 'color'. There are occasional instances where the sensors become saturated at the highest intensities. I remember hearing at one point that astronomers "Tail Fit" a Gaussian distribution so that it strictly fits the unsaturated data, the saturated data is then ignored.

How are these fits computed?

**To put in more specific terms:**

You can simulate this by taking a photograph of a lit lamp. Your image will likely result an RGB value per pixel, most home cameras map to color ranges in [0, 255]. When any channel that registers as (255) the sensor is considered saturated. When all channels register as (255, 255, 255) then this is "saturated white". Saturated regions are unable to provide useful information because true intensity measurements would not have a ceiling.

Models of intensity (per channel) will allow for identification of types of light. In order to make realistic models of intensity we have to correct the saturated data. This would encourage us to fit a curve as informed by only the unsaturated tails.

We want a curve that approximates the saturated section, informed by the unsaturated section. How do we "tail fit" this curve? Is this a strategy that exists or is it a fabrication I have remembered?

I currently expect this to be a Gaussian curve on each channel, however I'm happy to be suggested a different known curve.

**Tags:** star, identify-this-object, image-processing, distribution-function, color

> **Comments:**
> **username_1**: Just to clarify, this is for a 2-dimensional image? This maybe of use in crafting an answer: [Point-Spread Function Fitting Photometry](https://ui.adsabs.harvard.edu/abs/1999ASPC..189...56H/abstract)  
> **username_2**: "we have to censor the saturated data" -- maybe you mean "correct the saturated data"? (I fixed the other instances of "censor" -- where you clearly meant *sensor* -- but wasn't sure about that last phrase.)  

---

## 1 answers

### Answer by username_3 on 2023-12-19. Score: 2 | Accepted answer

Take your 2-d image of a star; mask out any pixels that may contain saturated data; fit a Gaussian (or whatever your model for the point spread function is) to the remaining pixels, with the amplitude and source position as free parameters; replace the masked pixels with values computed from the best-fit model you've just calculated.

This will work pretty well as long as the point spread function is not dependent on source intensity and can perhaps be derived from nearby, unsaturated stars.


---

