We have presented an atmospheric correction algorithm for HJ-1 CCD imagery over Lakes Taihu and Chaohu with highly turbid waters. The Rayleigh scattering radiance (Lr) is calculated using the hyperspectral Lr with a wavelength interval 1nm. The hyperspectral Lr is interpolated from Lr in the central wavelengths of MODIS bands, which are converted from the band response-averaged Lr calculated using the Rayleigh look up tables (LUTs) in SeaDAS6.1. The scattering radiance due to aerosol (La) is interpolated from La at MODIS band 869nm, which is derived from MODIS imagery using a shortwave infrared atmospheric correction scheme. The accuracy of the atmospheric correction algorithm is firstly evaluated by comparing the CCD measured remote sensing reflectance (Rrs) with MODIS measurements, which are validated by the in situ data. The CCD measured Rrs is further validated by the in situ data for a total of 30 observation stations within ± 1h time window of satellite overpass and field measurements. The validation shows the mean relative errors about 0.341, 0.259, 0.293 and 0.803 at blue, green, red and near infrared bands.