Dear Roger, I want to make a raster for my shapefile. I am try to apply function "rasterize". However I get a problem of the unit of plot. If I make a plot for my shapefile, then the unit is "degree". However, when I apply raster() and make a plot for it, the unit is not "degree". Please see the two plots I attached. I do not know if there is any problem... Is it due to the proj4string? And how to fix it? Many thanks for your time! Best regads Xijia -------------- next part -------------- A non-text attachment was scrubbed... Name: Rplot1.png Type: image/png Size: 7760 bytes Desc: Rplot1.png URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20140921/541e824b/attachment.png> -------------- next part -------------- A non-text attachment was scrubbed... Name: Rplot01.png Type: image/png Size: 10369 bytes Desc: Rplot01.png URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20140921/541e824b/attachment-0001.png>