3/7/24 - Meet about GEE code for S2 turbidity
- Relax/remove cloud threshold and mask clouds instead - get as many images as possible
- ROI inside bay - Can't see it.
- Time series starts in 2019 according to plot(?)
- Need to multiply B4 (and all S2 bands) vals by 0.0001 to get reflectance
- Implement algorithm below to get turbidity.
- Identify images associated with spikes in time series
- Display those images from (6) individually and look to see what is going on
- From Chowdhury et al. 2023 (https://doi.org/10.3389/fmars.2023)
- Look at B4 at 665nm and B5 at Rrs665
- If Rrs665<0.042: Turbidity=TurbidityRrs665
- If Rrs665>0.042 and Rrs704>0.062: Turbidity=TurbidityRrs704
- If Rrs665≥0.042 and Rrs704≤0.062: Turbidity = Turbidityweighted
- Turbidityweighted=α×TurbidityRrs665+β×TurbidityRrs704
- α = 0.062−Rrs665/0.062−0.042 , β=1−α
Interesting way to plot transect data:
Equations from Chowdhurry:
/srv/imars-objects/tpa_pgs/usvi/L3_A250
- L1A files ordered 4/22/24 to cover PR and USVI (17.5 to 18.9N; -67.5 to -64E)
- Files that won't process: 2020159 - 2022038