I want to create 50km buffers around a large set of polylines that have a global extent. My objective is to calculate average raster values within each buffer zone.
The matter of confusion: I want to use an equidistant projection to minimize the buffer distortion introduced by the global extent of my data, but ultimately I wish to use an equal area projection to calculate zonal statistics of the buffer zones.
How do I reconcile these two types of distortion? I cannot find a related question that explicitly addresses distortion when calculating zonal statistics of buffer zones.