I apologize, I have one more question actually.
The docs
specify: "This class also defines a framework for decomposing the unit sphere into a hierarchy of 'cells'". How do I actually create the top-level cell to subdivide?
My current approach is to create six different top-level cells rather than just one using "from_face_pos_level" (note that my code is in python):
> top_level_cells = [cell = s2sphere.Cell.from_face_pos_level(i,0,0) for i in range(6)]
However, I think this is wrong, because one of these six cells contains a point in Ghana as well as a point in Guam (on opposite sides of the globe). Here is my code to repro that bug:
> p1 = s2sphere.LatLng(6.005519, 0.5326073).to_point() # Ghana
> p2 = s2sphere.LatLng(13.33591,144.7627).to_point() # Guam
> cell = s2sphere.Cell.from_face_pos_level(0,0,0)
> cell.contains(p1), cell.contains(p2)
OUTPUT: True, True
Plotting some points in "cell" (from above) yields this:
It thus seems that "s2sphere.Cell.from_face_pos_level(0,0,0)" is not the correct way to create a top-level cell. What is the correct way to create a top-level cell?
Sorry for the disturbance - I have been reading the docs and trying various things for a few hours, but to no avail. I would tremendously appreciate anyone who could help unblock me!