<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="+1">Dear Maria</font>,<br>
I have found this note by chance, because I do not follow all the
traffic in the GRASS mailing list. <br>
I have to admit I have some difficulties to grasp the meaning of the
steps you describe, for example the abbreviations PR and RRN, will
have to read the articles you mention first. Generally, if the
process includes radiometric normalization, and correlation-based
normalization is applicable, use of my script should be possible. <br>
But for now, I have found full text of the article Baig et al.,
2014, on Researchgate:
<a class="moz-txt-link-freetext" href="https://www.researchgate.net/publication/262005316_Derivation_of_a_tasselled_cap_transformation_based_on_Landsat_8_at-_satellite_reflectance">https://www.researchgate.net/publication/262005316_Derivation_of_a_tasselled_cap_transformation_based_on_Landsat_8_at-_satellite_reflectance</a><br>
Hope it helps.<br>
<br>
Tomas B.<br>
<br>
<br>
<div class="moz-cite-prefix">Dne 4.7.2015 v 23:26 Kozlova Maria
napsal(a):<br>
</div>
<blockquote cite="mid:2111671436045176@web11m.yandex.ru" type="cite">
<pre wrap="">Dear colleagues,
Some time ago I asked if it is possible to derive custom TC coefficients by Grass and described the first algorithm for derivation of TC coefficients in attached file.
The steps needed to derive TC coefficients for multispectral image are not completely uniform.
Here is the description of second algorithm I'd like to follow for derivation of TC coefficients
The second recipe, based on TM tasseled cap transformation.
(Suggested by Huang et al., 2002, Lobster and Cohen, 2007 and Baig et al., 2014)
1.
Selection of ‘dummy target' (as it is called in a paper by Baig et al., 2014) by applying of TM TC-coefficients directly to OLI or ETM+ data Here, in Grass, I can choose sensor type TM for OLI data, so this step is seems to be possible to do in Grass by i.tasscap.
2.
PCA (it is clear for Grass and easily can be done in the images of interest by i.pca).
3.
PR rotation of PCA matrix seems to be the most problematic step. Generally it can be done by the finding of best alignment of pca matrix created in step2 with previously created (in step1) dummy target. Such a way: - more then 4 TC components should be included to i.tasscap - some kind of raster comparison and relative normalization must be done, I guess, but I don't know, what tool exactly should be used.
A little question directly to Tomas Brunclic:
Dear Tomas: what do you think about the application of your script on RRN (i.grid.correl.atcor....) in this case? I would be very appreciate for your answer.
4.
Get TC-coefficients by applying i.mapcalc to original Landsat data and resulting PR rotated data Such a way, TC coefficients arise from: PCA axes rotation coefficients PR of PCA matrix coefficients.
P.S.: TCT of Landsat MSS data can be done as I see only by Jackson's algorithm. So still hoping for any help with these two TCT algorithms.
</pre>
</blockquote>
<br>
</body>
</html>