Problem 667. Image Processing 01: Vignetting Correction

Created by Richard Zapor

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","relationshipId":"rId2","target":"/matlab/output.xml"}],"parts":[{"partUri":"/matlab/document.xml","relationship":[],"contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?><w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>The task is to correct image files for IR scanning sensors that due to errant tolerancing and structural instability have varying partially vignetted(V) images. Fortunately the V artifact is in a converging ray bundle and only the edges are impacted to differing amounts. Also, fortunately the sensor has a Low source flat field calibration target that can be used to determine the Nearfield(NF) for each image. The obscuration varies cross-scan but is constant for all samples. The NF varies with the structure temperature. To repair the optics in the drones is cost prohibitive. Fortuitously, the Image Processing is based on Matlab and needs only a fast function to perform an In-Scene image adjustment good enough for govt work.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Output: Fixed image(Fimg) (Array &lt;= 512x512 of type double)</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Input: Corrupted Image file(Cimg) (double); Low Signal NF_Cal (1xncol)</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Tolerance: Good enough is 1%+1 count of Truth (Timg)</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Assume Ideal LCal file where the center is a DC offset, not Vignetting.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>LCal provides V Offset. Obscuration: Max is 25% at edge, Min 8 columns wide.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Corrupted Image creation: Start with Truth (Timg)</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>CImg = Timg(1-Obscuration)+ NF_Vignetting</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>example: Obs = [.25 .1 0...0 .1 .25]; NF_V [50 20 0...0 20 50];</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>load penny.mat; Timg=double(P);</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Cimg=Timg.*repmat(1-Obs,128,1)+repmat(NF_V,128,1);</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>NF_Cal=NF_V+5;% NF_Cal [55 25 5 5...5 25 55]</w:t></w:r></w:p></w:body></w:document>"},{"partUri":"/matlab/output.xml","contentType":"text/xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?><embeddedOutputs><metaData><evaluationState>manual</evaluationState><layoutState>code</layoutState><outputStatus>ready</outputStatus></metaData><outputArray type=\"array\"/><regionArray type=\"array\"/></embeddedOutputs>"}]}

The task is to correct image files for IR scanning sensors that due to errant tolerancing and structural instability have varying partially vignetted(V) images. Fortunately the V artifact is in a converging ray bundle and only the edges are impacted to differing amounts. Also, fortunately the sensor has a Low source flat field calibration target that can be used to determine the Nearfield(NF) for each image. The obscuration varies cross-scan but is constant for all samples. The NF varies with the structure temperature. To repair the optics in the drones is cost prohibitive. Fortuitously, the Image Processing is based on Matlab and needs only a fast function to perform an In-Scene image adjustment good enough for govt work.Output: Fixed image(Fimg) (Array &lt;= 512x512 of type double)Input: Corrupted Image file(Cimg) (double); Low Signal NF_Cal (1xncol)Tolerance: Good enough is 1%+1 count of Truth (Timg)Assume Ideal LCal file where the center is a DC offset, not Vignetting.LCal provides V Offset. Obscuration: Max is 25% at edge, Min 8 columns wide.Corrupted Image creation: Start with Truth (Timg)CImg = Timg(1-Obscuration)+ NF_Vignettingexample: 
Obs = [.25 .1 0...0 .1 .25]; NF_V [50 20 0...0 20 50];load penny.mat;
Timg=double(P);Cimg=Timg.*repmat(1-Obs,128,1)+repmat(NF_V,128,1);NF_Cal=NF_V+5;% NF_Cal [55 25 5 5...5 25 55]

The task is to correct image files for IR scanning sensors that due to errant tolerancing and structural instability have varying partially vignetted(V) images. Fortunately the V artifact is in a converging ray bundle and only the edges are impacted to differing amounts. Also, fortunately the sensor has a Low source flat field calibration target that can be used to determine the Nearfield(NF) for each image. The obscuration varies cross-scan but is constant for all samples. The NF varies with the structure temperature. To repair the optics in the drones is cost prohibitive. Fortuitously, the Image Processing is based on Matlab and needs only a fast function to perform an In-Scene image adjustment good enough for govt work.

Output: Fixed image(Fimg) (Array <= 512x512 of type double)

Input: Corrupted Image file(Cimg) (double); Low Signal NF_Cal (1xncol)

Tolerance: Good enough is 1%+1 count of Truth (Timg)

Assume Ideal LCal file where the center is a DC offset, not Vignetting.

LCal provides V Offset. Obscuration: Max is 25% at edge, Min 8 columns wide.

Corrupted Image creation: Start with Truth (Timg)

CImg = Timg(1-Obscuration)+ NF_Vignetting

example: 
Obs = [.25 .1 0...0 .1 .25]; NF_V [50 20 0...0 20 50];

load penny.mat;
Timg=double(P);

Cimg=Timg.*repmat(1-Obs,128,1)+repmat(NF_V,128,1);

NF_Cal=NF_V+5;% NF_Cal [55 25 5 5...5 25 55]

Solve

Solution Stats

16 Solutions
4 Solvers

Last Solution submitted on Nov 08, 2021

Last 200 Solutions

Problem Comments

Problem Recent Solvers4

Problem Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Problem 667. Image Processing 01: Vignetting Correction

Solution Stats

Last 200 Solutions

Problem Comments

Problem Recent Solvers4

Suggested Problems

More from this Author255

Problem Tags

Community Treasure Hunt

Problem 667. Image Processing 01: Vignetting Correction

Solution Stats

Last 200 Solutions

Problem Comments

Problem Recent Solvers4

Suggested Problems

More from this Author255

Problem Tags

Community Treasure Hunt

Players