Problem 843. Hyperspectral Processing: Determine Material Components given a Hyperspectral vector

Created by Richard Zapor

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","targetMode":"","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","targetMode":"","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\"?>\n<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>Given a hyperspectral data set and Reflectance Spectral Signature Library determine a pixel's component percentages.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:hyperlink w:docLocation=\"http://aviris.jpl.nasa.gov/aviris/index.html\"><w:r><w:t>NASA AVIRIS</w:t></w:r></w:hyperlink></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>A Ground Square is imaged by hundreds of pixels, each at a different wavelength. The signal on pixel 1(500nm to 505nm) is a sum of the components (Concrete/Tree/Grass...) by percentage of area covered times the material reflectance. Pixel 2 (510-515nm) is different by the Reflectance deltas between Concrete and a Tree.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Let S(i,j) be the response of Material i for band j</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>g( j )= %(Concrete)*S(1,j)+%(Tree)*S(2,j)+...%(Grass)*S(end,j);</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>A 300-Band 9-Material Spectral file is loaded. Comparison between foliage and rocks is quite significant. The materials are Bush, Calcite, Concrete, Conifer, Grass (not that type), Fir tree, Gypsum, Maple, Sage</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:rPr><w:b/></w:rPr><w:t>g=S*f</w:t></w:r><w:r><w:t> where f is the percentage of the imaged pixel covered by the material.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:rPr><w:b/></w:rPr><w:t>Input:</w:t></w:r><w:r><w:t> g spectral sum [301,1]; S spectral material response [301,9] Nine materials</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:rPr><w:b/></w:rPr><w:t>Output:</w:t></w:r><w:r><w:t> Solve for f ....( eg f=[0 .5 0 .25 .25 0 0 0 0]' )</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>( f should sum to 1, max(f) is 1 and min(f) is 0 )</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>The test Suite will round to 2 decimal places. Cases of \"other materials\" which will induce negative values are not tested.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>This is introductory and ignores atmospheric absorption.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>There is a matrix operation hint in the test suite for a method to solve for f.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:hyperlink w:docLocation=\"http://aviris.jpl.nasa.gov/data/free_data.html\"><w:r><w:t>AVARIS Free Data</w:t></w:r></w:hyperlink><w:r><w:t> These data files are large with 224 bands x 750 channels x 2000 samples</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>To expand these files may require a tar converter</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"http://aviris.jpl.nasa.gov/alt_locator/111013_AV_Download.readme\"><w:r><w:t>NASA readme</w:t></w:r></w:hyperlink><w:r><w:t> ...and... </w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"http://aviris.jpl.nasa.gov/alt_gulf/\"><w:r><w:t>NASA Tools bottom Left</w:t></w:r></w:hyperlink><w:r><w:t> There are some possible issues with the NASA tar tool. Two non-standard files can be found at</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"http://dll-files.org/7968/index.html\"><w:r><w:t>libiconv-2.dll</w:t></w:r></w:hyperlink><w:r><w:t> and</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"http://dll-files.org/7975/libintl-2.dll.html\"><w:r><w:t>libintl-2.dll</w:t></w:r></w:hyperlink></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>See the Test Suite for details on opening the AVIRIS Moffett Field file.</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>"}]}

Given a hyperspectral data set and Reflectance Spectral Signature Library determine a pixel's component percentages.<a href="http://aviris.jpl.nasa.gov/aviris/index.html">NASA AVIRIS</a>A Ground Square is imaged by hundreds of pixels, each at a different wavelength.
The signal on pixel 1(500nm to 505nm) is a sum of the components (Concrete/Tree/Grass...) by percentage of area covered times the material reflectance.
Pixel 2 (510-515nm) is different by the Reflectance deltas between Concrete and a Tree.Let S(i,j) be the response of Material i for band jg( j )= %(Concrete)*S(1,j)+%(Tree)*S(2,j)+...%(Grass)*S(end,j);A 300-Band 9-Material Spectral file is loaded. Comparison between foliage and rocks is quite significant. The materials are Bush, Calcite, Concrete, Conifer, Grass (not that type), Fir tree, Gypsum, Maple, Sageg=S*f where f is the percentage of the imaged pixel covered by the
material.Input: 
g spectral sum [301,1]; 
S spectral material response [301,9] Nine materialsOutput:
Solve for f ....( eg f=[0 .5 0 .25 .25 0 0 0 0]' )( f should sum to 1, max(f) is 1 and min(f) is 0 )The test Suite will round to 2 decimal places.
Cases of "other materials" which will induce negative values are not
tested.This is introductory and ignores atmospheric absorption.There is a matrix operation hint in the test suite for a method to solve for f.<a href="http://aviris.jpl.nasa.gov/data/free_data.html">AVARIS Free Data</a>
These data files are large with 224 bands x 750 channels x 2000 samplesTo expand these files may require a tar converter <a href="http://aviris.jpl.nasa.gov/alt_locator/111013_AV_Download.readme">NASA readme</a>
...and... <a href="http://aviris.jpl.nasa.gov/alt_gulf/">NASA Tools bottom Left</a>
There are some possible issues with the NASA tar tool. Two non-standard files can be found at <a href="http://dll-files.org/7968/index.html">libiconv-2.dll</a> and <a href="http://dll-files.org/7975/libintl-2.dll.html">libintl-2.dll</a>See the Test Suite for details on opening the AVIRIS Moffett Field file.

Given a hyperspectral data set and Reflectance Spectral Signature Library determine a pixel's component percentages.

<http://aviris.jpl.nasa.gov/aviris/index.html NASA AVIRIS>

A Ground Square is imaged by hundreds of pixels, each at a different wavelength.
The signal on pixel 1(500nm to 505nm) is a sum of the components (Concrete/Tree/Grass...) by percentage of area covered times the material reflectance.
Pixel 2 (510-515nm) is different by the Reflectance deltas between Concrete and a Tree.

Let S(i,j) be the response of Material i for band j

g( j )= %(Concrete)*S(1,j)+%(Tree)*S(2,j)+...%(Grass)*S(end,j);

A 300-Band 9-Material Spectral file is loaded. Comparison between foliage and rocks is quite significant. The materials are Bush, Calcite, Concrete, Conifer, Grass (not that type), Fir tree, Gypsum, Maple, Sage

*g=S*f*  where f is the percentage of the imaged pixel covered by the
material.

*Input:* 
g spectral sum [301,1]; 
S spectral material response [301,9]  Nine materials

*Output:*
Solve for f  ....( eg f=[0 .5 0 .25 .25 0 0 0 0]' )

( f should sum to 1, max(f) is 1 and min(f) is 0 )

The test Suite will round to 2 decimal places.
Cases of "other materials" which will induce negative values are not
tested.

This is introductory and ignores atmospheric absorption.

There is a matrix operation hint in the test suite for a method to solve for f.

<http://aviris.jpl.nasa.gov/data/free_data.html AVARIS Free Data>
These data files are large with 224 bands x 750 channels x 2000 samples

To expand these files may require a tar converter
<http://aviris.jpl.nasa.gov/alt_locator/111013_AV_Download.readme NASA readme>
...and... 
<http://aviris.jpl.nasa.gov/alt_gulf/ NASA Tools bottom Left>
There are some possible issues with the NASA tar tool. Two non-standard files can be found at <http://dll-files.org/7968/index.html libiconv-2.dll> and <http://dll-files.org/7975/libintl-2.dll.html libintl-2.dll>

See the Test Suite for details on opening the AVIRIS Moffett Field file.

Solve

Solution Stats

27 Solutions
9 Solvers

Last Solution submitted on Dec 24, 2014

Last 200 Solutions

Problem Comments

2 Comments

2 Comments

Alfonso Nieto-Castanon on 1 Feb 2013

links not working...

Alfonso Nieto-Castanon on 1 Feb 2013

also pretty impossible to hack, the testsuite produces errors before ever calling our function...

Problem Recent Solvers9

Problem Tags

application hyperspectral mldivide pseudo matrix inversion

Community Treasure Hunt

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

Start Hunting!