fit polynomial function on spectral data

[المشرف العام]

I want to fit a polynomial function (max. 3rd order) on each raster cell over all my spectral bands (Landsat 1-7) creating a new raster(stack) representing the coefficients. I got my data (including NA values) in a stack with 6 Layer (Landsat Band 1-7[excluding 6]).

I guess somehow I should tell the polynomial function on which spectral wavelength the bands are located

1. Landsat7 Wavelength (micrometers)
   • Band 1 0.45-0.52
   • Band 2 0.52-0.60
   • Band 3 0.63-0.69
   • Band 4 0.77-0.90
   • Band 5 1.55-1.75
   • Band 7 2.09-2.35

so that it can fit it properly. Has anyone an idea how to do that polynomial fitting of each cell and extracting the coefficients in R? Thanks for any help!



أكثر...