The aim of this study is to build two mathematical models of canine ionic currents specific to right atria and left atria. The canine left atria mathematical model was firstly modified from the Ramirez-Nattel-Courtemanche (RNC) model using the recently available experimental data of ionic currents and was further developed based on our own experimental data. A model of right atria was then built by considering the differences between right atria and left atria. The two developed models well reproduced the experimental data on action potential morphology, the rate dependence, and action potential duration restitution. They are useful for investigating the mechanisms underlying the heterogeneity of canine regional action potentials and would help the simulation of whole heart excitation propagation and cardiac arrhythmia in the near future.