By taking the urban greenland as testing objective,the GeoEye-1 fusion image as background, and using Moran’s index,the correlativity of greenland units has been measured.Then,the presampling was conducted based on the measured correlativity,the systematic sampling design has been carried out at the precision 95%and 85%.By using of Kriging interpolation method,the weight of different type of greenland has been divided into the stratifications.Through the stratified sampling,the estimation of greenland coverage has been obtained. The results show that the urban greenland coverage area is 14.81 km~2,accounting for 27.73%of all urban land areas.Compared with the common sampling method,the precision of that increased about 3%. This method can also be utilized in monitoring dynamic variation of urban greenland with the real time images. By taking the urban greenland as testing objective, the GeoEye-1 fusion image as background, and using Moran’s index, the correlativity of greenland units has been measured. Chen, the presampling was conducted based on the measured correlativity, the systematic sampling design has been By using the Kriging interpolation method, the weight of different type of greenland has been divided into the stratified sampling. The estimation of greenland coverage has been obtained. The results show that that the urban greenland coverage area is 14.81 km ~ 2, accounting for 27.73% of all urban land areas. Compared with the common sampling method, the precision of that increased about 3%. This method can also be utilized in monitoring dynamic variation of urban greenland with the real time images.