[GRASS-user] Derivation of Landsat 8 Tasseled Cap coefficients for certain scene in Grass
Kozlova Maria
kclo at yandex.ru
Sat Jul 4 14:26:16 PDT 2015
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.
--
Best regards,
Maria V. Kozlova
/------------------------------------\
Information System Lab. (ISysLab)
State Oceanographic Institute,
Kropotkinsky per., 6,
MOSCOW, 119034,
RUSSIA
------------------------------------
tel. +7 499 246 6448
fax. +7 499 246 7288
mailto: kclo at yandex.ru
\------------------------------------/
11.06.2015, 18:32, "Kozlova Maria" <kclo at ya.ru>:
> Dear colleagues,
>
> I attach the document containing description of one of algorithms for calculating TC coefficients.
> Can it be done by Grass?
>
> --
> Best regards,
>
> Maria V. Kozlova
>
> /------------------------------------\
> Information System Lab. (ISysLab)
> State Oceanographic Institute,
> Kropotkinsky per., 6,
> MOSCOW, 119034,
> RUSSIA
> ------------------------------------
> tel. +7 499 246 6448
> fax. +7 499 246 7288
> mailto: kclo at yandex.ru
> \------------------------------------/
>
> 11.06.2015, 11:43, "Kozlova Maria" <kclo at yandex.ru>:
>
>> Well...While I'm trying to type the algorithm step by step you may look through the following papers which can hopefully clarify the subject
>> And thank you again for your support!
>>
>> https://disk.yandex.ru/client/disk/%D0%94%D0%BE%D0%BA%D1%83%…
>>
>> And thank you again for your support!
>> --
>> Best regards,
>>
>> Maria V. Kozlova
>>
>> /------------------------------------\
>> Information System Lab. (ISysLab)
>> State Oceanographic Institute,
>> Kropotkinsky per., 6,
>> MOSCOW, 119034,
>> RUSSIA
>> ------------------------------------
>> tel. +7 499 246 6448
>> fax. +7 499 246 7288
>> mailto: kclo at yandex.ru
>> \------------------------------------/
>>
>> 10.06.2015, 21:58, "Kozlova Maria" <kclo at yandex.ru>:
>>> Well...
>>> While I'm trying to type the algorithm step by step you may look through the following papers which can hopefully clarify the subject
>>> And thank you again for your support!
>>>
>>> --
>>> Best regards,
>>>
>>> Maria V. Kozlova
>>>
>>> /------------------------------------\
>>> Information System Lab. (ISysLab)
>>> State Oceanographic Institute,
>>> Kropotkinsky per., 6,
>>> MOSCOW, 119034,
>>> RUSSIA
>>> ------------------------------------
>>> tel. +7 499 246 6448
>>> fax. +7 499 246 7288
>>> mailto: kclo at yandex.ru
>>> \------------------------------------/
>>>
>>> 10.06.2015, 14:58, "Moritz Lennert" <mlennert at club.worldonline.be>:
>>>
>>>> On 10/06/15 13:18, Nikos Alexandris wrote:
>>>>
>>>>> * Kozlova Maria <kclo at yandex.ru> [2015-06-10 11:48:44 +0300]:
>>>>>
>>>>>> Dear colleagues,
>>>>>>
>>>>>> Thank you for response.
>>>>>> a little bit later I will type the full procedure for derivation of TC coefficients.
>>>>>> PCA itself cannot be applied here as it is based only on statistical information of the image. TCT is physically based transformation and PCA is only the first step.
>>>>>> The main problem I see here is the 'soil line' which must be determined (manually or automatically) and PCA axes rotation to align 1st PCA component in accordance to the soil line.
>>>>>> Though I have to describe it in more details later...
>>>>>
>>>>> This introduction <http://www.sjsu.edu/faculty/watkins/tassel.htm> looks
>>>>> good for me :-)
>>>>
>>>> I don't have access to the paper used for the OLI version, but from the
>>>> abstract the authors do "just" use PCA + rotation to get a result that
>>>> is as close as possible to a direct tasselled cap transformation using
>>>> the original 1970s coefficients.
>>>>
>>>> And when you run PCA on MSS images you often get quite close results to
>>>> what Kauth and Thomas arrived at through their clustering approach.
>>>>
>>>> Moritz
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
More information about the grass-user
mailing list