search

nasa / HLS-Data-Resources

This repository provides guides, short how-tos, and tutorials to help users access and work with Harmonized Landsat Sentinel-2 (HLS) data.
Apache License 2.0
47 stars 6 forks source link

HLS tutorial (R markdown) issue with data from different UTM zones #17

Open pjhurska opened 5 months ago

pjhurska commented 5 months ago

@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
amfriesz commented 4 months ago

@pjhurska, thanks for pointing that out. Would you be able to provide the polygon data/file you're working with?

pjhurska commented 3 months ago

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.