We present the occurrence rates for rocky planets in the habitable zones (HZs) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define η ⊕ as the HZ occurrence of planets with radii between 0.5 and 1.5 R ⊕ orbiting stars with effective temperatures between 4800 and 6300 K. We find that η ⊕ for the conservative HZ is between 0.37-0.21+0.48 (errors reflect 68% credible intervals) and 0.60-0.36+0.90 planets per star, while the optimistic HZ occurrence is between 0.58-0.33+0.73 and 0.88-0.51+1.28 planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates between using Poisson likelihood Bayesian analysis and using Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with 95% confidence that, on average, the nearest HZ planet around G and K dwarfs is ∼6 pc away and there are ∼4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun.