A previously proposed method for estimating antisolvent crystallization kinetics has been successfully applied to cooling crystallization to yield kinetic parameters for simulating the crystallization process. In particular, the primary nucleation kinetics were analyzed using the modi ed induction time data, essentially as proposed by Kubota. The number density at the detection point, (N/M) det , is an important value in the kinetic analysis but is di cult to estimate experimentally. Therefore, the numerical optimization was used to estimate the value of (N/M) det in the same manner as for antisolvent crystallization. Primary nucleation was detected for both antisolvent and cooling crystallization with a sensitive detection method, visual detection, to reduce the e ects of secondary nucleation. As a result, the (N/M) det values for each method are small and almost identical (approximately 200 #/kg-solvent). The secondary nucleation and growth rate parameters for cooling crystallization were also successfully determined by numerical optimization. All the kinetic parameters determined for cooling crystallization were evaluated with experimental data. Consequently, it is con rmed that these rate parameters can simulate trends in concentration as well as the nal number mean diameter of the product crystals with acceptable accuracy. Moreover, the crystallization rate parameters determined for both antisolvent and cooling crystallization were validated with simulations and experimental data from combined crystallization, where antisolvent crystallization is followed by cooling crystallization. The simulation and experimental results for the concentration trend and number mean diameter of the produced crystal were in good agreement. The applicability of our estimation method to both antisolvent and cooling crystallization indicates the broad utility of this method for various crystallization process in the pharmaceutical industry.