function make_hierarchy_matrix(dconn_file, seed_voxels, target_voxels, output_filename)

data_dconn = readimgfile(dconn_file);
data_hierarchy = [];

for i_seed = 1:length(seed_voxels)
    for i_target = 1:length(target_voxels)
        data_hierarchy.correlation(i_seed,i_target) = data_dconn(seed_voxels(i_seed),target_voxels(i_target));
        data_hierarchy.seed_exemplar(i_seed,i_target) = seed_voxels(i_seed);
        data_hierarchy.target_exemplar(i_seed,i_target) = target_voxels(i_target);       
    end
end

save(output_filename, 'data_hierarchy')

end