Cold Jupiters play a crucial role in planet formation and dynamical evolution. Since their initial discovery around 47 UMa, they have attracted significant interest;" however, their formation mechanisms remain uncertain, underscoring the need to expand the known population. In this work, we combine radial velocity (RV) data with Gaia astrometry using Hipparcos-Gaia proper-motion anomalies over a 25 yr baseline. By jointly modeling both data sets with the Markov Chain Monte Carlo framework, we constrain planetary masses, orbital inclinations, and three-dimensional orbital architectures. This reduces RV degeneracies and improves mass determinations. Four cold Jupiters are reported: HD 68475 b and HD 100508 b are each the first confirmed planet in their systems, with orbital periods 7832(-323)(+463) days and 5681 +/- 42 days and dynamical masses of 5.16(-0.47)(+0.53) MJup and 1.2(-0.18)(+0.30) MJup, respectively. In multiplanet systems, HD 48265 c has a period of 10,418(-1400)(+2400) days and a mass of 3.71(-0.43)(+0.68) MJup, while HD 114386 c orbits at 444.00-0.88+0.93 days with a minimum mass of 0.37(-0.03)(+0.03) MJup. The two planets in the HD 48265 system may exhibit a significant mutual inclination, making it a target for testing the von Zeipel-Kozai-Lidov mechanism. HD 68475 b is a promising candidate for future direct imaging with the Extremely Large Telescope's Mid-infrared ELT Imager and Spectrograph. We identified a Jupiter analog with the longest known orbital period among planets with masses between 0.5 and 2 M-Jup, implying that a substantial population of cold Jupiters likely awaits discovery by Gaia. This study expands the sample of cold Jupiters with constrained orbits and dynamical masses, demonstrating the value of combining RV and astrometry in exoplanet research."