https://camera-optical-zoom.blogspot.com/2022/09/speed-up-numpy-nested-loop.html