r - レベルプロットが正しくプロットされず、値がグリッドに重ねられない

Question

汚染物質に関するデータをグリッド化しました。各グリッドは、左端の緯度と経度で示されます。この汚染物質のレベルがグリッド間でどのように変化するかを確認したいと思います。このために、レベルプロットを作成しました。以下に示す私のコードは問題なく動作しますが、唯一の問題は、プロットされた汚染物質が abline を使用して作成されたグリッドに重ならないことです。x と y の限界カットオフを指定しましたが、値はこれらのカットオフより先にプロットされます。

    library(lattice)
    library(maptools)
    imap = readShapeSpatial("IND_adm3")
    pdf("trial2.pdf",width=11, height=8)
    levelplot(surface_aod~longitude+latitude,data=s.data,
                  panel = function(...){
                  panel.levelplot(...)
                  panel.abline(v=(71:88),col="dark red")              
                             panel.abline(h=(17:31),col="dark red")    
                   sp.polygons(imap)         
                     }, 
   col.regions = heat.colors(100),aspect="iso", region=TRUE,scales=list(x=list(at=seq(from=71,to=88, by=1)), y=list(at=seq(from=17,to=31,   by=1)),cex=.7, alternating=3),xlim=c(70,90), ylim=c(16,32))
   dev.off()  

ここに私のデータのスニペットがあります。

 dput(head(s.data))
 structure(list(latitude = c(17L, 18L, 19L, 19L, 20L, 21L), longitude = c(80L, 
 82L, 80L, 81L, 82L, 75L), surface_aod = c(0.2652681, 0.0040855, 
 0.2032517, 0.0929442, 0.1194997, 0.3600747)), .Names = c("latitude", 
 "longitude", "surface_aod"), row.names = c(NA, 6L), class = "data.frame")

グリッドにオーバーレイする値を取得するにはどうすればよいですか? ありがとう。