Landing footprint, as basis for further mission planning, is crucial in the aerospace vehicle reentry phase especially in some emergency scenario. The no-fly zone constraints make the boundary of the landing footprint complicated and bring difficulties to the solution. In this paper, a landing footprint calculation algorithm is proposed based on the classification of no-fly zones. The constant bank angle method is employed as fundamental method of solving landing footprint. First, the no-fly zones are classified through different avoidance methods. Second, after calculating the initial landing footprint tangent trajectories and tangent points are found by bi-partition. Third, the sub-areas are calculated according to the states of tangent points. Finally, the bank angle reversal method is presented to obtain limit boundaries. Mathematical simulation validates the effectiveness of the proposed algorithm.