g.region rast=lsat7_2002_10 -p # Find pure endmembers, they are supposed to be in the corners of # the PCA feature space (or use spectral measurement data) i.pca -n input="lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70" out=pca_lsat7_2002 d.correlate layer1=pca_lsat7_2002.1 layer2=pca_lsat7_2002.2 # TODO: problem: how to obtain the unprojected coordinates for the corner pixels? # (in 1998 the xgobi software did this nicely, check today's ggobi) i.group group=lsat7_2002 subgroup=lsat7_2002 \ input=lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70 i.spec.unmix group=lsat7_2002 matrix=sample/spectrum.dat result=lsat7_2002_unmix \ error=lsat7_2002_unmix_err iter=lsat7_2002_unmix_iterations # todo: reclass to 0..100%
Mohammed Rashad (rashadkm gmail.com) (2012, update to GRASS 6)
Last changed: $Date$
Main index - imagery index - Full index
© 2003-2016 GRASS Development Team