|
|
##### rinpan_romou.r ## Cドライブに rgis というフォルダをつくり関連するファイルを保存 setwd("C:/rgis") ### ワーキング・ディレクトリの設定 source("GIS_source.R") ### GIS関連の関数の読み込み rinp<-trk("演習林林班.trk") ### ベクターデータの読み込み rind<-trk("演習林林道.trk") sag<-trk("演習林作業道.trk") hod<-trk("演習林歩道.trk") ### 平面直角座標系範囲 X1<-3200 ; X2<-6300 ; Y1<--2600 ; Y2<-2200 ### 白地図の作成 plot(NA, xlim=c(X1,X2),ylim=c(Y1,Y2),xlab="x (m)",ylab="y (m)", main="位山演習林の林班および路網") ### 林班と路網を描く関数map() map<-function(col1="blue",col2="black"){ lapply(rinp,lines,col=col1) lapply(rind,lines,col=col2,lwd=3) lapply(sag ,lines,col=col2,lwd=2,lty=1) lapply(hod ,lines,col=col2,lwd=2,lty=2) } map() ### 林班と路網を描画 |