A method and a system for determining a center of resistance point of a tooth for orthodontic treatment planning are provided. The method comprises: obtaining a tooth mesh from image data associated with a tooth crown of a patient; identifying a mesiodistal center of the tooth crown; determining a reference plane in the image data to extend through the mesiodistal center; determining an intersection curve based on an intersection of the reference plane and the tooth mesh, the intersection curve following a shape of the surface of the crown at the reference plane; determining a tooth axis of the tooth crown based on the intersection curve; determining a crown height of the tooth crown based on the tooth axis; and determining the center of resistance of the tooth based on the determined crown height and the determined tooth axis.