Currently I have a Landsat land cover map with 30-m resolution and I want to calculate the fractional cover for each vegetation class (forests, grasslands and croplands etc.) matching MODIS swath and do a linear regression analysis linking land cover percentage with the corresponding MODIS LST for each 1-km resolution pixel. Therefore, the fractional changes will be calculated by computing at 30m the proportion of each vegetation class (e.g. forest) that occupy each MODIS 1km grid cell.
I hardly found any practical tutorial using R to achieve this.