Thank you for highlighting this situation as I have not previously looked at it for thermal lenses 
Thanks to
LesioQ and our cruel, but beautiful universe!

The game doesn't end with calculating the airy disc diameter ! That's a common mistake.
What do you want to do next?
--> sample it at more than just one pixel per airy diameter!
Because you're not really interested in sampling the airy disc with one pixel, but sampling the tiny contrast (in amplitude AND space) resulting from merging two close diffracted spots.
You can sample at a pitch much less than the airy diameter.
I found that document that appears to explains it quite well (I had a quick glimpse at it but it seems good).
First, thank you for this article, it is really mind changing. I recommend everyone to read it.
Second, you are absolutely right, we can sample at a small pitch to get maximum possible spatial resolution that current lens is capable to provide. I found quite nice picture from one astronomy forum, that shows the advantage of high accurate airy disk sampling.
That one says a pitch of 1/6 of the airy diameter !
Here is a quote from the article, explaining why we are talking about 1/6 of the airy diameter for those, who didn't find it:
Since we know from the Rayleigh criterion that Δxmin equals the radius of the Airy disc, it follows that the best
‘match’ between sensor characteristics and the optics is achieved when the detector pitch p equals one‐third of
the Airy disc radius. In other words: three pixels suffice to resolve the Airy disc radius entirely. The reasoning
then goes that above that threshold, the optics diffraction limitation kicks in. So digitising the radius of the PSF’s
central part with more than three pixels does not improve the final spatial resolution of the image because the
detector only resolves diffraction blur at that point.
BTW, I didn't now, but FLIR has a special UltraMax technology, that allows to increase the resolution. It is very silimilar to Fluke's superresolution.
The most interesting moment is at 0:51, the main idea is the same, just sample the data in gaps between pixels, as active pixel area is really small.
Full video:
https://youtu.be/8y2lV3nIhvc