[RS-ecology] rasterCVA gives a blank map as output

Benjamin Leutner benjamin.leutner at uni-wuerzburg.de
Thu Oct 6 09:59:50 CEST 2016


Natalia,

it looks like your tasseled cap transformation of 1988 is off by a 
factor of about 10000.
If I calculate it from scratch it is fine:

a <- brick("raster_data/final/p224r63_1988.grd")
b <- brick("raster_data/final/p224r63_2011.grd")

ac <- tasseledCap(dropLayer(a, "B6_bt"), sat = "Landsat5TM")
bc <- tasseledCap(dropLayer(b, "B6_bt"), sat = "Landsat5TM")

rva <- rasterCVA(ac[[1:2]], bc[[1:2]])


Best,
Benjamin


On 14.09.2016 21:16, Natalia Quinteros wrote:
> Hi List,
>
> In chapter 9, when attempting a spectral change vector analysis the 
> rasterCVA gave as an input a blank map.
>
> changeVec <- rasterCVA(tc2011[[1:2]], tc1988[[1:2]])
>
> Any ideas why?
>
> Thanks,
>
> Natalia
>
>
> _______________________________________________
> RS-ecology mailing list
> RS-ecology at lists.uni-wuerzburg.de
> https://lists.uni-wuerzburg.de/mailman/listinfo/rs-ecology

-- 
Benjamin Leutner M.Sc.
PhD candidate

Department of Remote Sensing
University of Wuerzburg
Campus Hubland Nord 86
97074 Wuerzburg, Germany

Tel: +49-(0)931-31 89594
Fax: +49-(0)931-31 89594-0
Email: benjamin.leutner at uni-wuerzburg.de
Web: http://www.fernerkundung.uni-wuerzburg.de

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.uni-wuerzburg.de/pipermail/rs-ecology/attachments/20161006/09272187/attachment.html>


More information about the RS-ecology mailing list