GeoWithinQ is a very useful function - is this point within this area? However, it is quite slow. It's not unheard of to have large numbers of points - for instance, global satellite measurements at a high resolution. I have used satellite data with millions of pixels - ie, data points.
Let's say I have 100,000 points and I would like to check whether they are within the polygon of France. However, doing
points = RandomGeoPosition Select[points, GeoWithinQ[Entity["Country","France"]]]
times out. The documentation has another example that uses
Pick instead of
Pick[points, GeoWithinQ[Entity["Country","France"], points]
This too times out, but it is significantly faster for a smaller number of points (I've managed about 1,000 before hitting major performance issues).
My question is, how do I use
GeoWithinQ for 100,000 points and more?