Embodiments of a system and method for indexing of geospatial data using three-dimensional Cartesian space are generally described herein. In an aspect, such example methods may include calculating endpoints of a segment, wherein the endpoints are specified in Cartesian coordinates and are located on a substantially spherical surface, defining a boundary of a polygon according to the segment, computing one or more normals corresponding to one or more planes, wherein each of the one or more planes contain a test point and a boundary point associated with the boundary, obtaining a boundary sine value of an angle defined by an arc subtended by the endpoints, summing each of a group of angle values derived from the boundary sine value to obtain an angle sum, wherein the group contains the boundary sine value, and determining whether the test point is inside the polygon based on the angle sum.