Rを使用して、すでに色付きのポリゴンが含まれている別のシェープファイルの上にポリゴンのシェープファイルをプロットしようとしています。maptoolsとrgdalを使用しています。下部のシェープファイル（複数のポリゴンがある）を、4つのポリゴンがある上部のシェープファイル（私は思う）で区切られた領域にクリップしたいと思います。もちろん、plot（）内でcol =コマンドを使用して、一番上のポリゴンの内側を塗りつぶす方法を知っています。しかし、私がやりたいのは、上部のポリゴンの外側を白で塗りつぶし、下部の色付きのポリゴンを上部のポリゴンの内側から見えるようにすることです。これができたら、下のポリゴンの境界を再描画します。私はR2.13.2と最新版のmaptoolsとrgdal、およびそれらのすべての依存関係をWindows7マシンで使用しています。
これらのシェープファイルに依存するコードがあります：ftp： //ftpext.usgs.gov/pub/er/wi/la.crosse/McKann/bcr%20arcinfo%20files/
library(maptools)
sm=readShapeSpatial('states.shp')
pm=readShapeSpatial('province.shp')
library(rgdal)
rgm=readOGR(dsn=wd,layer='gwwa_new_range_5Jul11')
rgm2=spTransform(rgm,CRS(proj4string(sm2)))
plot(c(-120,-61),c(35,55),type='n',axes=F,bty='n',xlab='',ylab='', main=sp)
plot(sm,col=c('red','green','yellow'),add=T)##us states
plot(pm,col=c('red','green','yellow'),add=T)##canadian provinces
plot(rgm2,lwd=2,add=T,border='blue')##gives the range of this bird
範囲マップの内側を色にすることができます
plot(rgm2,add=T,col='blue')
しかし、私がやりたいのは、外側を白にして（後で州と州の境界線を描画できるようにするため）、内側を透明にして、範囲内の州と州の色を確認できるようにすることです。
さらに詳しい情報が必要な場合はお知らせください。サポートが必要な場合は事前にお知らせください。たぶんもっと簡単な方法があります...