I have a global tif file in epsg:4326 and a shapefile of points distributed around the global. I want to retrieve the sum of the value of pixels around 10km of each point.
One obvious thing to do is to convert tiff to points and intersect in an equal area projection. Unfortunately conversion to points takes ages so i abandoned it.
Then focal statistics with the corresponding buffer and getting value for each point is an option, but what projection do you think is best for this?
I used mollweide and it seems that when I measure the pixels in any area of the world they have the same size.