sid

Measure spectral similarity using spectral information divergence

Since R2020a

Syntax

score = sid(inputData,refSpectra)

score = sid(testSpectra,refSpectra)

Description

score = sid(inputData,refSpectra) measures the spectral similarity between the spectra of each pixel in the hyperspectral data inputData and the specified reference spectra refSpectra by using the spectral information divergence (SID) technique. Use this syntax to identify different regions or materials in a hyperspectral data cube.

score = sid(testSpectra,refSpectra) measures the spectral similarity between the specified test spectra testSpectra and reference spectra refSpectra by using the SID method. Use this syntax to compare the spectral signature of an unknown material against the reference spectra or to compute spectral variability between two spectral signatures.

Note

This function requires the Hyperspectral Imaging Library for Image Processing Toolbox™. You can install the Hyperspectral Imaging Library for Image Processing Toolbox from Add-On Explorer. For more information about installing add-ons, see Get and Manage Add-Ons.

The Hyperspectral Imaging Library for Image Processing Toolbox requires desktop MATLAB^®, as MATLAB Online™ or MATLAB Mobile™ do not support the library.

Examples

Distinguish Hyperspectral Regions Using Spectral Information Divergence

Distinguish different regions in a hyperspectral data cube by computing the spectral information divergence (SID) between each pixel spectrum and the endmember spectrum of the data cube.

Read hyperspectral data into the workspace.

 hcube = hypercube('jasperRidge2_R198.hdr');

Specify the number of spectrally distinct bands to identify in the data cube.

numEndmembers = 7;

Extract endmember spectral signatures from the data cube by using the NFINDR algorithm.

endmembers = nfindr(hcube,numEndmembers);

Plot the spectral signatures of the endmembers.

figure  
plot(endmembers)   
xlabel('Band Number')
ylabel('Data Value')
legend('Location','Bestoutside')

Compute the spectral information divergence between each endmember and the spectrum of each pixel in the data cube.

score = zeros(size(hcube.DataCube,1),size(hcube.DataCube,2),numEndmembers);
for i= 1:numEndmembers
    score(:,:,i) = sid(hcube,endmembers(:,i));
end

Compute the minimum score value from the distance scores obtained for each pixel spectrum with respect to all the endmembers. The index of each minimum score identifies the endmember spectrum to which a pixel spectrum exhibits maximum similarity. An index value, n, at the spatial location (x, y) in the score matrix indicates that the spectral signature of the pixel at spatial location (x, y) in the data cube best matches the spectral signature of the nth endmember.

[~,matchingIndx] = min(score,[],3);

Estimate an RGB image of the hyperspectral data cube by using the colorize function. Display both the RGB image and the matrix of matched index values.

rgbImg = colorize(hcube,'Method','RGB');
figure('Position',[0 0 1100 500])
subplot('Position',[0 0.2 0.4 0.7])
imagesc(rgbImg)
axis off
colormap default
title('RGB Image of Hyperspectral Data')
subplot('Position',[0.5 0.2 0.4 0.7])
imagesc(matchingIndx);
axis off
title('Indices of Matching Endmembers')
colorbar

Determine Similarity of Endmember Spectra Using SID

Read hyperspectral data into the workspace.

hcube = hypercube("jasperRidge2_R198.hdr");

Find 10 endmembers of the hyperspectral data cube by using the N-FINDR method.

numEndmembers = 10;
endmembers = nfindr(hcube,numEndmembers);

Consider the first endmember as the reference spectrum and the rest of the endmembers as the test spectrum. Compute the SID score between the reference and test spectra.

score = zeros(1,numEndmembers-1);
refSpectrum = endmembers(:,1);
for i = 2:numEndmembers
    testSpectrum = endmembers(:,i);
    score(i-1) = sid(testSpectrum,refSpectrum);
end

Find the test spectrum that exhibit maximum similarity (minimum distance) to the reference spectrum. Then find the test spectrum that exhibit minimum similarity (maximum distance) to the reference spectrum.

[minval,minidx] = min(score);
maxMatch = endmembers(:,minidx);
[maxval,maxidx] = max(score);
minMatch = endmembers(:,maxidx);

Plot the reference spectrum, maximum similarity test spectrum, and the minimum similarity test spectrum. The test spectrum with the minimum score has the highest similarity to the reference endmember. On the other hand, the test spectrum with maximum score has the highest spectral variability and characterises the spectral behavior of two different materials.

figure
plot(refSpectrum)
hold on
plot(maxMatch,"k")
plot(minMatch,"r")
legend("Reference spectrum","Minimum match test spectrum","Maximum match test spectrum", ...
       Location="southoutside");
title("Similarity Between Spectra")
annotation("textarrow",[0.3 0.3],[0.4 0.52],String="Min score: "+minval)
annotation("textarrow",[0.6 0.6],[0.4 0.55],String="Max score: "+maxval)
xlabel("Band Number")
ylabel("Data Values")

Input Arguments

inputData — Input hyperspectral data

hypercube object | 3-D numeric array

Input hyperspectral data, specified as a hypercube object or a 3-D numeric array containing the data cube. If the input is a hypercube object, the data is read from the DataCube property of the object.

testSpectra — Test spectra

C-element vector

Test spectra, specified as a C-element vector. The test spectra is the spectral signature of an unknown region or material.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

refSpectra — Reference spectra

C-element vector

Reference spectra, specified as a C-element vector. The reference spectra is the spectral signature of a known region or material. The function matches the test spectra against these values.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Output Arguments

score — SID score

scalar | matrix

SID score, returned as a scalar or matrix. The output is a

• scalar — If you specify the testSpectra input argument. The function matches the test spectral signature against the reference spectral signature and returns a scalar value. Both the test and the reference spectra must be vectors of same length.

• matrix — If you specify the inputData input argument. The function matches the spectral signature of each pixel in the data cube against the reference spectral signature and returns a matrix. If the data cube is of size M-by-N-by-C and the reference spectra is a vector of length C, the output matrix is of size M-by-N.

A smaller SAM score indicates a strong match between the test signature and the reference signature.

Data Types: single | double

Limitations

This function does not support parfor loops, as its performance is already optimized. (since R2023a)

Algorithms

The sid function normalizes the reference spectra refSpectra and test spectra testSpectra and computes the SID value using this formula.

$$SID={\displaystyle \sum _{i=1}^C{p}_i\log \left(\frac{p_i}{q_i}\right)+}{\displaystyle \sum _{i=1}^C{q}_i\log \left(\frac{q_i}{p_i}\right)}.$$

q and p are the vectors of normalized reference and test spectra, respectively. q_i and p_i are the i^th elements of the vectors q and p, respectively. C is the length of vectors q and p.

References

[1] Chein-I Chang. “An Information-Theoretic Approach to Spectral Variability, Similarity, and Discrimination for Hyperspectral Image Analysis.” IEEE Transactions on Information Theory 46, no. 5 (August 2000): 1927–32. https://doi.org/10.1109/18.857802.

Version History

Introduced in R2020a