I have a large point data set that has irregular geometry (random sampling) that I want to bin in 2D XY space. I've attached a sketch below to illustrate what I mean. Apologies for hand drawing, hope no one minds.
So far I've been using the fishnet option in ArcGIS+QGIS and SAGA tool polygon attribute to points. However I have reached the workable file size limit for my shapefiles (>2GB). When attempting to use feature classes ArcGIS is incredibly slow at performing the required task, so I have been mostly focusing on SAGA for its speed.
It's been suggested to me by a colleague that I could use GMT's blockmean, blockmedian, blockmode programs. However I have limited experience in this but would be willing to learn IF it is the right direction
I can see the maths of all this as a matrix operation, but the problem I am having is with the large data size. My point files are >700mb in just XY coordinates (billions of points).
I am not looking for interpolation, just looking to assign the points a 'bin' reference in XY space.