HLS tutorial (R markdown) issue with data from different UTM zones
Closed this issue · 3 comments
@amfriesz When running the HLS tutorial with own polygon data (just one polygon geometry), which happens to situate on the border of two UTM zones, the solution on rows 502-509 to extract the CRS from the first asset in the search_df list does not work, as the list contains images from two different UTM zones, i.e from two different CRS.
Therefore, the code to "create a list of raster layers for each of our bands of interest (i.e., Red, NIR, and Fmask)" on rows 521-560 runs to a certain point, until it encounters data from that other zone, then failing with this error message:
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'x' in selecting a method for function 'mask': [crop] extents do not overlap
@pjhurska, thanks for pointing that out. Would you be able to provide the polygon data/file you're working with?
Actually I ended up writing own solution, following your excellent tutorial: https://github.com/sykefi/Habitat_monitoring_HLS
Any improvements and suggestions also here are most welcome!
BTW: I was just wondering, why did you remove the cloud masking from your tutorial? I find it quite important.
Hey @pjhurska, we identified an issue with the quality filtering in the original tutorial. We have added the updated quality filtering back to the HLS tutorial.