A matrix is called totally positive (resp. totally nonnegative) if all its minors are positive (resp. nonnegative). Consider the Ising model with free boundary conditions and no external field on a planar graph G. Let a_1,dots ,a_k,b_k,dots ,b_1 be vertices placed in a counterclockwise order on the outer face of G. We show that the ktimes k matrix of the two-point spin correlation functions Mi,j=⟨σaiσbj⟩\\documentclass[12pt]{minimal}\t\t\t\t\\usepackage{amsmath}\t\t\t\t\\usepackage{wasysym}\t\t\t\t\\usepackage{amsfonts}\t\t\t\t\\usepackage{amssymb}\t\t\t\t\\usepackage{amsbsy}\t\t\t\t\\usepackage{mathrsfs}\t\t\t\t\\usepackage{upgreek}\t\t\t\t\\setlength{\\oddsidemargin}{-69pt}\t\t\t\t\\begin{document}$$\\begin{aligned} M_{i,j} = \\langle \\sigma _{a_i} \\sigma _{b_j} \\rangle \\end{aligned}$$\\end{document}is totally nonnegative. Moreover, det M > 0 if and only if there exist k pairwise vertex-disjoint paths that connect a_i with b_i. We also compute the scaling limit at criticality of the probability that there are k parallel and disjoint connections between a_i and b_i in the double random current model. Our results are based on a new distributional relation between double random currents and random alternating flows of Talaska [37].