According to the geomagnetic observation data in China for 2003 from 135 geomagnetic stations and 35 geomagnetic observatories, and the geomagnetic data of 38 IGRF calculated points in China's adjacent regions, the distribution of the geomagnetic anomalous field in China was calculated. Two kinds of geomagnetic models are selected as the geomagnetic normal field: one is the spherical harmonic model of the international geomagnetic reference field, and the other is the Taylor polynomial model of the geomagnetic field in China. Based on the geomagnetic anomaly values (the observed value minus the calculated value of the model) at various stations, the spherical cap harmonic model of the geomagnetic anomalous field was calculated by using the spherical cap harmonic analysis method. The geomagnetic anomaly charts (∆D, ∆I, ∆F, ∆X, ∆Y and ∆Z) in China for 2003 are drawn. The geomagnetic anomalous field in China is analyzed and discussed.