Closed jakobjassmann closed 2 years ago
Logan: Once you have a chance, it would be great if you could:
Yesterday I expanded the readme section. Still need to develop fully self-contained examples...
On the documentation thread, I am a little confused about what lsat_calibrate_rf is returning in the resulting data table. I think (?) there is a new column called [band].cal (e.g. ndvi.xcal) which I assumed was the 'corrected' ndvi that one would want to used in the downstream analyses. But then with this interpretation I think one would want ndvi.xcal to be used in the next steps in the readme/vignette, but it's still ndvi. So I think I am missing something on the documentation end that would clarify this, perhaps others as well?
Good eye @scelmendorf ! Thanks for pointing that out. Yes, the appended [band].xcal data are the corrected ndvi. As you noted, you'll either want to use ndvi.xcal in the subsequent functions, or replace the uncalibrated data with the calibrated data under the ndvi column. I updated the documentation to make that workflow clearer.
More documentation: I don't fully understand the consequences of assigning points or polygons to different regions. I think (?) the RF cross-validation works by site(?). And I understand how we have often used the region concept for synthesis papers. But analytically, I don't know what it does here.
@scelmendorf This is a really helpful comment! There is no meaning behind assigning points or polygons to regions. When I read your comment I was surprised that you got that impression, but I re-read the readme and I think that the second example for lsat_get_pixel_centers()
is perhaps confusing to the reader.
In the example the identifier "region" is used to provide a unique identifier to the different polygons, which can then be used to split the multi row polygon sf object into unique polygons to map lsat_get_pixel_centers()
over each of them. The column name "region" carries no meaning at all. I just chose "region" instead of "site", as Logan requested that we reserve "site" for the individual point coordinates. Perhaps calling the column "polygon_id" in the example would be less confusing?
Very happy to hear any other suggestions to improve clarity on this end! :)
more spelling minutia: "parameter" is misspelled in documentation for lsat_fit_phenological_curves
spar
Smoothing paramater passed to smooth.spline(), typically around 0.65 - 0.75 for this application.
Thanks @scelmendorf! Fixed the "parameter" typo.
Continued in #28 .
Jakob:
Logan: