In the analysis of spatial interactions, their representation can significantly contribute to understanding of the discussed spatial phenomenon. In the graduation thesis we automated the procedures of some, in the literature proposed standard spatial interaction mapping methods, and proposed some of our own solutions. We have automated spatial interaction mapping methods using spider web method, net origins and net destinations and vectors of asymmetry. In addition to these, we also automated the calculations and procedures of mapping spatial interactions spider web gross and net connections, destination emphasis, origin emphasis, interpolation of vectors of asymmetry approach and the difference of vectors of asymmetry. The automation of the mapping was performed in the ArcGIS software environment, and the code is written in the Python programming language. The results of this diploma thesis is the ArcGIS software environment toolkit and instructions for mapping of spatial interactions with examples, which are directly applicable in spatial planning processes.