#3865: i.spec.unmix: map writeout broken
--------------------------+-------------------------
Reporter: neteler | Owner: grass-dev@…
Type: defect | Status: new
Priority: normal | Milestone:
Component: Imagery | Version: unspecified
Keywords: i.spec.unmix | CPU: x86-64
Platform: Unspecified |
--------------------------+-------------------------
For unknow reasons i.spec.unmix is unable to write out the resulting maps:
{{{
GRASS 7.7.dev (nc_spm_08_grass7):~ >
g.region raster=lsat7_2002_10 -p
[...]
nsres: 28.5
ewres: 28.5
rows: 475
cols: 527
cells: 250325
i.spec.unmix group=lsat_2002 matrix=spectrum.txt result=lsat7_2002_unmix
error=lsat7_2002_unmix_err iter=lsat7_2002_unmix_iterations
D0/0: buff: Matrix: 5 by 7
Opening input file no. 1 [lsat7_2002_10]
Opening input file no. 2 [lsat7_2002_20]
Opening input file no. 3 [lsat7_2002_30]
Opening input file no. 4 [lsat7_2002_40]
Opening input file no. 5 [lsat7_2002_50]
Opening input file no. 6 [lsat7_2002_70]
Opening input file no. 7 [lsat7_2002_80]
Opening output file [lsat7_2002_unmix.1]
Opening output file [lsat7_2002_unmix.2]
Opening output file [lsat7_2002_unmix.3]
Opening output file [lsat7_2002_unmix.4]
Opening output file [lsat7_2002_unmix.5]
Opening error file [lsat7_2002_unmix_err]
Opening iteration file [lsat7_2002_unmix_iterations]
Checking linear dependencies (orthogonality check) of Matrix A...
Angle between row 1 and row 2: 60.9753 degree
Angle between row 1 and row 3: 56.5953 degree
Angle between row 1 and row 4: 34.7591 degree
Angle between row 1 and row 5: 43.122 degree
Angle between row 2 and row 1: 60.9753 degree
Angle between row 2 and row 3: 42.7788 degree
Angle between row 2 and row 4: 35.7377 degree
Angle between row 2 and row 5: 40.2744 degree
Angle between row 3 and row 1: 56.5953 degree
Angle between row 3 and row 2: 42.7788 degree
Angle between row 3 and row 4: 32.3279 degree
Angle between row 3 and row 5: 28.9019 degree
Angle between row 4 and row 1: 34.7591 degree
Angle between row 4 and row 2: 35.7377 degree
Angle between row 4 and row 3: 32.3279 degree
Angle between row 4 and row 5: 28.552 degree
Angle between row 5 and row 1: 43.122 degree
Angle between row 5 and row 2: 40.2744 degree
Angle between row 5 and row 3: 28.9019 degree
Angle between row 5 and row 4: 28.552 degree
Spectral matrix is o.k. Proceeding...
Calculating for 475 x 527 pixels (7 bands) = 277729 pixelvectors.
100%
ERROR: Unable to open header file for raster map
<lsat7_2002_unmix.1@PERMANENT>
}}}
Since i.spec.unmix was ported from GRASS 5 years ago, perhaps some raster
function usage is not correct?
I have attached the input matrix of endmembers for easy testing.
Ideas welcome!
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3865>
GRASS GIS <https://grass.osgeo.org>