This study aimed to identify subunits of the basal ganglia and thalamus and to investigate the functional connectivity among these anatomically segregated subdivisions and the cerebral cortex in healthy subjects. For this purpose, we introduced multilevel independent component analysis (ICA) of the resting-state functional magnetic resonance imaging (fMRI). After applying ICA to the whole brain gray matter, we applied second-level ICA restrictively to the basal ganglia and the thalamus area to identify discrete functional subunits of those regions. As a result, the basal ganglia and the thalamus were parcelled into 31 functional subdivisions according to their temporal activity patterns. The extracted parcels showed functional network connectivity between hemispheres, between subdivisions of the basal ganglia and thalamus, and between the extracted subdivisions and cerebral functional components. Grossly, these findings correspond to cortico-striato-thalamo-cortical circuits in the brain. This study also showed the utility of multilevel ICA of resting state fMRI in brain network research.