Automatic topological relation generating algorithm influences directly the performance of such spatial operations in GIS as spatial relation establishment and query. Double neighboring point judgement method can be adopted so as to decrease computing complexity greatly in two key stages of this algorithm -determining both direction of polygon and point-in-polygon relation.