Based on the precipitation data from 95 regional meteorological stations in Huangshi City since 2005 and 3 national meteorological stations, a stepwise regression analysis was conducted to simulate and impute the missing precipitation data of each station from 1958 to date, with national meteorological stations data benchmarked and 12 factors adjusted including the precipitation of each station and the temperature of Huangshi station. After imputing the missing precipitation data in different time scales, the REOF method was used to analyze the temporal and spatial annual precipitation variation in Huangshi City. The zoning method of taking the eigenvalue with the largest absolute value as the site classification value was used to divide Huangshi City into hills and water body affected areas (Area I) and mountain land affected area (Area Ⅱ). The results show that (1) decadal variation of annual precipitation is obvious, as the precipitation gradually increased from 1960s to 1990s, decreased at the beginning of 21st century and recovered to the level of 1980s and 1990s from 2010 to date. (2) Distinct regional characteristics were founded in annual precipitation data with greater influence of topography. The annual precipitation in the mountainous area of southwestern Huangshi is larger than that of hilly area and large water body affected area. (3) Precipitation intensity in Zone I and Zone Ⅱ is synchronized 72% of all the years in records, while unsynchronized years only constitute 28%.