search

durbank / PAIPR

Functions to generate probabilistic estimates of annual accumulation from ice-penetrating radar without the need for manual layer selection or correction
3 stars 0 forks source link

Improved stacking of data #22

Closed durbank closed 3 years ago

durbank commented 5 years ago

It may be advantageous to improve the stacking routine in PAIPR by taking account the local layer slope and stacking along this trendline rather than directly horizontally. In this case, the processing order would likely be

  1. Conversion to depth
  2. Stationarize the data (remove trend and set constant variance)
  3. Set depth cutoff, interpolate to core resolution, and clip to given cutoff
  4. Radon transform the data (possibly with a larger window?) and estimate layer grad field
  5. Data stacking, using the RT data to determine stacking (increase the window to 50 m?)
  6. Vertical data smoothing
  7. Peak finding
  8. Layer picking (and etc...)
durbank commented 3 years ago

Although this would indeed be helpful, it is not a top priority right now. I will come back to this as time permits (which it most likely won't).