This paper presents an improved algorithm for retrieving ocean surface currents from X-band marine radar images. The original polar current shell (PCS) method begins with a 3D fast Fourier transform (FFT) of the radar image sequence, followed by the extraction of the dispersion shell from the 3D image spectrum, which is then transformed into a PCS using polar coordinates. Building on this foundation, the improved approach is to analyze all data points corresponding to different wavenumber magnitudes in the PCS domain rather than analyzing each specific wavenumber magnitude separately. In addition, kernel density estimation (KDE) to identify high-density directions, interquartile range filtering to remove outliers, and symmetry-based filtering to further reduce noise by comparing data from opposite directions are also utilized for further improvement. Finally, a single curve fitting is applied to the filtered data rather than conducting multiple curve fittings as in the original method. The algorithm is validated using simulated data and real radar data from both the Decca radar, established in 2008, and the Koden radar, established in 2017. For the 2008 Decca radar data, the improved PCS method reduced the root-mean-square deviation (RMSD) for speed estimation by 0.06 m/s and for direction estimation by 3.8° while improving the correlation coefficients (CCs) for current speed by 0.06 and direction by 0.07 compared to the original PCS method. For the 2017 Koden radar data, the improved PCS method reduced the RMSD for speed by 0.02 m/s and for direction by 4.6°, with CCs being improved for current speed by 0.03 and direction by 0.05 compared to the original PCS method.