Visualizing presorted point cloud data


I have a big point cloud data that cannot be loaded into memory directly. I want to achieve that I can read the data within a bounding box without loading all of the data. One solution can be 1.) sort the point cloud data according to their distance to origin in the hillbert curve. 2.) when reading the data, convert the bounding box to the range of the hillbert distance and load the data within this range. I believe spatialpandas have this function, but I do not find documentation and examples related to it.

Do anybody have met the same questions or you just happen to know the answer? I will appreciate it if you tell me. I have tried to read the source code of spatialpandas but it is too hard for me to understand. Thanks!