> No problem! > Make a sites file from your raster file: > r.stats -1zg input=raster_filename | s.in.ascii sites=new_sitelist_name > -Malcolm Williamson Hi, I tried the above program, but only got a empty site file. What is the problem? Thanks. Houxiang Wang