This paper focuses on Kalman-Yakubovich-Popov lemma for multidimensional systems described by Roesser model that possibly includes both continuous and discrete dynamics. It is shown that, similarly to the standard 1-D case, this lemma can be studied through the lens of S-procedure. Furthermore, by virtue of this lemma, we will examine robust stability, bounded and positive realness of multidimensional systems.