#kmlOverlay library(maptools) library(rgdal) kmlfile <- "c:/GISdata/kmlOverlay.kml" kmlpng <-"c:/GISdata/kmlOverlay.png" qk <- SpatialPointsDataFrame(quakes[, c(2:1)], quakes) proj4string(qk) <- CRS("+proj=longlat") SGqk <- GE_SpatialGrid(qk,maxPixels = 600)#<=ここがミソ png(file=paste(kmlpng, sep=""), width=SGqk$width, height=SGqk$height,bg="transparent") par(mar=c(0,0,0,0), xaxs="i", yaxs="i") plot(qk, xlim=SGqk$xlim, ylim=SGqk$ylim,pch="+",col="red") dev.off() kmlOverlay(SGqk, kmlfile,kmlpng) system(paste('"c:/Program Files/Google/Google Earth/GoogleEarth.exe"',kmlfile),wait=FALSE) #writeOGR library(maptools) library(rgdal) x<-readShapePoly("C:/GISdata/gunma.shp") #x<-readOGR("C:/GISdata/gunma.shp","gunma") proj4string(x)<-CRS("+proj=longlat") #proj4string(x)<-CRS("+init=epsg:4326") kmlfileV <- "C:/GISdata/kmlOverlayVector.kml" writeOGR(x,kmlfileV,"gunma","KML") #system(paste('"c:/Program Files/Google/Google Earth/GoogleEarth.exe"',kmlfileV),wait=FALSE) #EPSGコード epsg<-make_EPSG() epsg[1,1:3] #プロジェクションの変換 library(maptools) x<-readShapePoly("C:/GISdata/gunma.shp") proj4string(x)<-CRS("+init=epsg:4612") plot(x)