# Date: 2018-02-01
# Author: Yi-Ying Chen
# Purpose: to get the mask area from a reference raster map
# shp2raster funtion was adopted from Amy's blog
# https://amywhiteheadresearch.wordpress.com/2014/05/01/shp2raster/
# get nc file for migration plot
inputfileA <- c("D:/RScripts/migration_map/1994.nc")
inputfileB <- c("D:/RScripts/migration_map/2015.nc")
library(maptools)
## Warning: package 'maptools' was built under R version 3.4.3
## Loading required package: sp
## Warning: package 'sp' was built under R version 3.4.3
## Checking rgeos availability: FALSE
## Note: when rgeos is not available, polygon geometry computations in maptools depend on gpclib,
## which has a restricted licence. It is disabled by default;
## to enable gpclib, type gpclibPermit()
library(raster)
## Warning: package 'raster' was built under R version 3.4.3
library(rgdal)
## Warning: package 'rgdal' was built under R version 3.4.3
## rgdal: version: 1.2-16, (SVN revision 701)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 2.2.0, released 2017/04/28
## Path to GDAL shared files: C:/Users/ychen/Documents/R/win-library/3.4/rgdal/gdal
## GDAL binary built with GEOS: TRUE
## Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
## Path to PROJ.4 shared files: C:/Users/ychen/Documents/R/win-library/3.4/rgdal/proj
## Linking to sp version: 1.2-7
#set nc file path
inputfile <- c("D:/RScripts/migration_map/2015.nc")
#inputfile <- "D:/RScripts/migration_map/CWB_5km_grid_1904.nc"
# Grab the lat and lon from the data
lat <- raster(inputfile, varname="nav_lat")
## Loading required namespace: ncdf4
lon <- raster(inputfile, varname="nav_lon")
# Convert to points and match the lat and lons
# offset <- 0.05
offset <- 0.008
plat <- rasterToPoints(lat)+ offset
plon <- rasterToPoints(lon)- offset
lonlat <- cbind(plon[,3], plat[,3])
# Specify the lonlat as spatial points with projection as long/lat
lonlat <- SpatialPoints(lonlat, proj4string = CRS("+proj=longlat +datum=WGS84"))
# Need the rgdal package to project it to the original coordinate system
# library("rgdal")
# My best guess at the proj4 string from the information given
mycrs <- CRS("+proj=longlat +datum=WGS84 +no_defs")
plonlat <- spTransform(lonlat, CRSobj = mycrs)
# Take a look
plonlat
## class : SpatialPoints
## features : 672300
## extent : 119.9688, 123.6921, 21.88386, 25.96349 (xmin, xmax, ymin, ymax)
## coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
extent(plonlat)
## class : Extent
## xmin : 119.9688
## xmax : 123.6921
## ymin : 21.88386
## ymax : 25.96349
# set the coordinate to the raster variable imported from nc file
twRaster <- raster(x = inputfile, varname="LU_TYPE")
projection(twRaster) <- mycrs
extent(twRaster) <- extent(plonlat)
# set all value to zero
twRaster[!is.na(twRaster)] <- 0
# plot(twdRaster)
# extract Taipei city polygons and convert to a world raster where cells
# associated with Australia have a value of 1 and everything else has a
# value of 0.
water_mp <- readShapeSpatial("D:/RScripts/migration_map/Taiwan_part_84.shp")
## Warning: use rgdal::readOGR or sf::st_read
## Warning: use rgdal::readOGR or sf::st_read
proj4string(water_mp) <- "+proj=longlat +datum=WGS84" # specify projection
#plot(water_mp)
# show the plot
ld_go <- FALSE
if (ld_go) {
lu_plot <- as.matrix(raster(inputfile, varname="LU_TYPE"))
LU_P <- lu_plot
LU_P[LU_P==0 | LU_P >16 ] <- 0
LU_P[LU_P==1] <- 1.2
LU_P[LU_P==4] <- 2.2
LU_P[LU_P==2] <- 3.2
LU_P[LU_P==13] <- 4.2
LU_P[LU_P==8] <- 5.5
LU_P[LU_P>6 | LU_P<0] <-0
breakpoints <- c(0,1,2,3,4,5,6)
colors <- c("lightgray","#006400","#FFB90F","#A2CD5A","#0000FF","#FF1493")
layout(matrix(data=c(1),nrow=1, ncol=1, byrow=TRUE),
widths=c(1,1,1), heights=c(1,1,1))
par( oma=c(2,2,0,2), mgp=c(2,1,0), mar=c(3.0,3.0,0.0,2.0),plt = c(0.1,0.9,0.1,0.90))
r_p <-raster(LU_P,
xmn=extent(plonlat)[1], xmx=extent(plonlat)[2],
ymn=extent(plonlat)[3], ymx=extent(plonlat)[4],
crs=CRS("+proj=longlat +datum=WGS84") )
plot(r_p, breaks=breakpoints,col=colors,
lab.breaks=c("","Forest","Agri.",
"Grass","to Water","B/Soil",""),
xlab="", ylab="",
xlim=c(119.8,122),ylim=c(21.8,25.3),cex=0.8)
plot(water_mp,lwd=2.0,add=T)
mtext("Longitude (deg)",side=1,line=2.5)
mtext("Latitude (deg)",side=2,line=2.5)
p_text <- paste("year 2015" )
mtext(p_text,side=3,line=1.0,cex=2.0)
}
m_text <- c("North reiogn", "Center region","East region","South region", "All regions")
for ( ii in (1:5) ) {
# assign the region
if (ii< 5){region_name <- as.integer(ii)}
else{region_name <- seq(1,4,1)}
# got the mask
region_mask <- rasterize(x = water_mp[region_name, ],
y = twRaster,
field=1)
dir_path = c("D:/RScripts/migration_map/")
y_st <-substr(inputfileA,start=(nchar(inputfileA)-6),stop = (nchar(inputfileA)-3))
y_ed <-substr(inputfileB,start=(nchar(inputfileB)-6),stop = (nchar(inputfileB)-3))
pp_text<- m_text[ii]
pdf_fname=c(paste( dir_path, y_st,"_",y_ed,"_region_",pp_text,".pdf",sep="") )
print(paste("i",ii,m_text[ii]),sep="")
pdf(pdf_fname,width=8, height=10)
layout(matrix(data=c(1),nrow=1, ncol=1, byrow=TRUE),
widths=c(1,1,1), heights=c(1,1,1))
par( oma=c(2,2,0,2), mgp=c(2,1,0), mar=c(3.0,3.0,0.0,2.0),plt = c(0.1,0.9,0.1,0.90))
# assign data
LU_A <- as.matrix(raster(inputfileA, varname="LU_TYPE"))*as.matrix(region_mask)
LU_B <- as.matrix(raster(inputfileB, varname="LU_TYPE"))*as.matrix(region_mask)
# diff array calculation
DIF_LU <- (LU_B-LU_A)
#
# calculate the gross change matrix
# Forest(1) Agriculture(4) Grass(2) Water(13) City/BareSoil(8)
gross_MATRIX <- array(0, dim=c(5,5))
#
# Forest Agri Grass Water City
#
# F 0 -3 -1 -12 -7
# A 3 0 2 -9 -4
# G 1 -2 0 -11 -6
# W 12 9 11 0 5
# C 7 4 6 -5 0
id <- c(1, 4, 2, 13, 8)
for (i in 1:5) {
for (j in 1:5) {
# change pixcel information
dif_value <- (id[i]-id[j])
org_value <- id[j]
#print(paste("org:",org_value, "dif:",dif_value,sep=""))
gross_MATRIX[j,i] <- length(DIF_LU[which(DIF_LU == dif_value & LU_A == org_value )])
#if(i==j) gross_MATRIX[j,i]<- 0
}
}
factor <- (.5*.5)/1000
total <- length(LU_A[which(LU_A==1)])*factor
gross_MATRIX = gross_MATRIX*factor
gross_MATRIX
# convert matrix to raster
r_dif <-raster(DIF_LU,
xmn=extent(plonlat)[1], xmx=extent(plonlat)[2],
ymn=extent(plonlat)[3], ymx=extent(plonlat)[4],
crs=CRS("+proj=longlat +datum=WGS84") )
LU_P <- DIF_LU
LU_P[LU_P==0 | LU_P >=13 ] <- 0
LU_P[LU_P==-3 | LU_P==-1 | LU_P==-12 | LU_P==-7] <- 1.2
LU_P[LU_P==3 | LU_P==2 | LU_P==-9 | LU_P==-4] <- 2.2
LU_P[LU_P==1 | LU_P==-2 | LU_P==-11 | LU_P==-6] <- 3.2
LU_P[LU_P==12 | LU_P==9 | LU_P==11 | LU_P==5] <- 4.2
LU_P[LU_P==7 | LU_P==4 | LU_P==6 | LU_P==-5] <- 5.2
LU_P[LU_P< 0] <- 0
r_p <-raster(LU_P,
xmn=extent(plonlat)[1], xmx=extent(plonlat)[2],
ymn=extent(plonlat)[3], ymx=extent(plonlat)[4],
crs=CRS("+proj=longlat +datum=WGS84") )
breakpoints <- c(0,1,2,3,4,5,6)
colors <- c("lightgray","#006400","#FFB90F","#A2CD5A","#0000FF","#FF1493")
#colors <- c("white","darkgreen","yellow","green","blue","pink")
plot(r_p, breaks=breakpoints,col=colors,
lab.breaks=c("no Change","to Forest","to Agri.",
"to Grass","to Water","to Built/S","N/A"),
xlab="", ylab="",
xlim=c(119.8,122),ylim=c(21.8,25.3),cex=0.8)
plot(water_mp,lwd=2.0,add=T)
mtext("Longitude (deg)",side=1,line=2.5)
mtext("Latitude (deg)",side=2,line=2.5)
p_text <- paste("from ", y_st, " to ", y_ed,sep="" )
mtext(p_text,side=3,line=1.0,cex=2.0)
dev.off()
}
## [1] "i 1 North reiogn"
## [1] "i 2 Center region"
## [1] "i 3 East region"
## [1] "i 4 South region"
## [1] "i 5 All regions"
library("circlize")
## Warning: package 'circlize' was built under R version 3.4.3
## ========================================
## circlize version 0.4.3
## CRAN page: https://cran.r-project.org/package=circlize
## Github page: https://github.com/jokergoo/circlize
## Documentation: http://jokergoo.github.io/circlize_book/book/
##
## If you use it in published research, please cite:
## Gu, Z. circlize implements and enhances circular visualization
## in R. Bioinformatics 2014.
## ========================================
#pdf_fname=c(paste( dir_path, y_st,"_",y_ed,"_gross_change_subregion.pdf",sep="") )
#pdf(pdf_fname,width=12, height=10)
layout(matrix(data=c(1,2,4,3),nrow=2, ncol=2, byrow=TRUE),
widths=c(1,1,1), heights=c(1,1,1))
par( oma=c(2,2,2,2), mgp=c(2,1,0), mar=c(3.0,3.,3.0,3.0),plt = c(0.05,0.95,0.05,0.95))
for ( ii in 1:4 ) {
# assign the region
if (ii< 5){region_name <- as.integer(ii)}
else{region_name <- seq(1,4,1)}
# got the mask
region_mask <- rasterize(x = water_mp[region_name, ],
y = twRaster,
field=1)
# assign data
LU_A <- as.matrix(raster(inputfileA, varname="LU_TYPE"))*as.matrix(region_mask)
LU_B <- as.matrix(raster(inputfileB, varname="LU_TYPE"))*as.matrix(region_mask)
# diff array calculation
DIF_LU <- (LU_B-LU_A)
#
# calculate the gross change matrix
# Forest(1) Agriculture(4) Grass(2) Water(13) City/BareSoil(8)
gross_MATRIX <- array(0, dim=c(5,5))
#
# Forest Agri Grass Water City
#
# F 0 -3 -1 -12 -7
# A 3 0 2 -9 -4
# G 1 -2 0 -11 -6
# W 12 9 11 0 5
# C 7 4 6 -5 0
id <- c(1, 4, 2, 13, 8)
for (i in 1:5) {
for (j in 1:5) {
# change pixcel information
dif_value <- (id[i]-id[j])
org_value <- id[j]
#print(paste("org:",org_value, "dif:",dif_value,sep=""))
gross_MATRIX[j,i] <- length(DIF_LU[which(DIF_LU == dif_value & LU_A == org_value )])
#if(i==j) gross_MATRIX[j,i]<- 0
}
}
factor <- (.5*.5)/1000
total <- length(LU_A[which(LU_A==1)])*factor
gross_MATRIX = gross_MATRIX*factor
gross_MATRIX
#Gorss chamge plot
mat <- gross_MATRIX
# the col and raw names
rownames(mat) = c("Forest","Agri","Grass","W","Built/Soil")
colnames(mat) = c("Forest","Agri","Grass","W","Built/Soil")
# try to have a circle plot
#chordDiagram(mat, grid.col = 1:5, directional = 1, row.col = 1:5)
#output the gross matrix data to a data.frame and save as a csv.file
df_matrix <- data.frame()
df_matrix <- setNames(data.frame(matrix(ncol = 3, nrow = 0)), c("From", "To", "Area") )
for (i in 1:length(rownames(mat)) ) {
for (j in 1:length(colnames(mat)) ) {
#
newRow <- data.frame(From=colnames(mat)[i],To=rownames(mat)[j],Area=mat[i,j] )
df_matrix <- rbind(df_matrix, newRow)
print(df_matrix)
} # end of j
} # end of i
write.table(x=df_matrix,
file="D:/RScripts/migration_map/gross_exchange_data.csv",
sep=",",
row.names =FALSE,
col.names=TRUE)
#read in the data.frame for both gross change and color information
df0 <- read.csv(file = "D:/RScripts/migration_map/gross_exchange_data.csv", stringsAsFactors=FALSE)
df1 <- read.csv(file = "D:/RScripts/migration_map/gross_exchange_color.csv", stringsAsFactors=FALSE)
# start plot
#plot( 0,0, xlim=c(-1,1),ylim=c(-1,1), type="n", axes=FALSE,
# main=m_text[i], xlab="", ylab="")
circos.par(start.degree = 120, clock.wise = TRUE)
chordDiagram(x = df0, grid.col = df1$col, transparency = 0.25,
order = df1$region, directional = 1,
direction.type = c("arrows", "diffHeight"), diffHeight = -0.05,
annotationTrack = "grid", annotationTrackHeight = c(0.03, 0.13),
link.arr.type = "big.arrow", link.sort = TRUE, link.largest.ontop = TRUE)
# add text and axis
circos.trackPlotRegion(
track.index = 1,
bg.border = NA,
panel.fun = function(x, y) {
xlim = get.cell.meta.data("xlim")
sector.index = get.cell.meta.data("sector.index")
reg1 = df1$reg1[df1$region == sector.index]
reg2 = df1$reg2[df1$region == sector.index]
circos.text(x = mean(xlim), y = ifelse(test = nchar(reg2) == 0, yes =8.0, no =5.0),
labels = reg1, facing = "bending", cex = 1.2)
# circos.text(x = mean(xlim), y = 3.1, labels = reg2, facing = "bending", cex = 1.2)
circos.axis(h = "top", major.at = seq(from = 0, to = xlim[2], by = ifelse(test = xlim[2]>5, yes = 2, no = 1)),
minor.ticks = 1, major.tick.percentage = 0.5,labels.niceFacing = TRUE,labels.cex = 1.0)
})
circos.clear()
mtext(m_text[ii],side=3,line=0.5,cex=1.5)
}
## From To Area
## 1 Forest Forest 4.133
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## 21 Built/Soil Forest 0.30950
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## 21 Built/Soil Forest 0.30950
## 22 Built/Soil Agri 0.06750
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## 21 Built/Soil Forest 0.30950
## 22 Built/Soil Agri 0.06750
## 23 Built/Soil Grass 0.03375
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## 21 Built/Soil Forest 0.30950
## 22 Built/Soil Agri 0.06750
## 23 Built/Soil Grass 0.03375
## 24 Built/Soil W 0.00350
## From To Area
## 1 Forest Forest 4.13300
## 2 Forest Agri 0.20225
## 3 Forest Grass 0.11800
## 4 Forest W 0.01350
## 5 Forest Built/Soil 0.21975
## 6 Agri Forest 0.65200
## 7 Agri Agri 0.40900
## 8 Agri Grass 0.16775
## 9 Agri W 0.00750
## 10 Agri Built/Soil 0.28525
## 11 Grass Forest 0.04800
## 12 Grass Agri 0.04900
## 13 Grass Grass 0.01325
## 14 Grass W 0.00150
## 15 Grass Built/Soil 0.04725
## 16 W Forest 0.05150
## 17 W Agri 0.04350
## 18 W Grass 0.01575
## 19 W W 0.01100
## 20 W Built/Soil 0.03825
## 21 Built/Soil Forest 0.30950
## 22 Built/Soil Agri 0.06750
## 23 Built/Soil Grass 0.03375
## 24 Built/Soil W 0.00350
## 25 Built/Soil Built/Soil 0.30750
## Note: 1 point is out of plotting region in sector 'Forest', track '1'.
## Note: 6 points are out of plotting region in sector 'Forest', track '1'.
## Note: 1 point is out of plotting region in sector 'Agri', track '1'.
## Note: 4 points are out of plotting region in sector 'Agri', track '1'.
## Note: 1 point is out of plotting region in sector 'Grass', track '1'.
## Note: 5 points are out of plotting region in sector 'Grass', track '1'.
## Note: 1 point is out of plotting region in sector 'W', track '1'.
##
## Note: 1 point is out of plotting region in sector 'W', track '1'.
## Note: 1 point is out of plotting region in sector 'Built/Soil', track '1'.
## Note: 10 points are out of plotting region in sector 'Built/Soil', track '1'.
## From To Area
## 1 Forest Forest 5.45625
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## 21 Built/Soil Forest 0.17925
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## 21 Built/Soil Forest 0.17925
## 22 Built/Soil Agri 0.58675
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## 21 Built/Soil Forest 0.17925
## 22 Built/Soil Agri 0.58675
## 23 Built/Soil Grass 0.04125
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## 21 Built/Soil Forest 0.17925
## 22 Built/Soil Agri 0.58675
## 23 Built/Soil Grass 0.04125
## 24 Built/Soil W 0.04700
## From To Area
## 1 Forest Forest 5.45625
## 2 Forest Agri 0.41575
## 3 Forest Grass 0.26775
## 4 Forest W 0.00625
## 5 Forest Built/Soil 0.14475
## 6 Agri Forest 0.37675
## 7 Agri Agri 1.18650
## 8 Agri Grass 0.11625
## 9 Agri W 0.02000
## 10 Agri Built/Soil 0.56550
## 11 Grass Forest 0.08000
## 12 Grass Agri 0.11350
## 13 Grass Grass 0.02350
## 14 Grass W 0.00975
## 15 Grass Built/Soil 0.09850
## 16 W Forest 0.06925
## 17 W Agri 0.05875
## 18 W Grass 0.00850
## 19 W W 0.03475
## 20 W Built/Soil 0.04175
## 21 Built/Soil Forest 0.17925
## 22 Built/Soil Agri 0.58675
## 23 Built/Soil Grass 0.04125
## 24 Built/Soil W 0.04700
## 25 Built/Soil Built/Soil 0.42750
## Note: 1 point is out of plotting region in sector 'Forest', track '1'.
## Note: 6 points are out of plotting region in sector 'Forest', track '1'.
## Note: 1 point is out of plotting region in sector 'Agri', track '1'.
## Note: 4 points are out of plotting region in sector 'Agri', track '1'.
## Note: 1 point is out of plotting region in sector 'Grass', track '1'.
## Note: 5 points are out of plotting region in sector 'Grass', track '1'.
## Note: 1 point is out of plotting region in sector 'W', track '1'.
##
## Note: 1 point is out of plotting region in sector 'W', track '1'.
## Note: 1 point is out of plotting region in sector 'Built/Soil', track '1'.
## Note: 10 points are out of plotting region in sector 'Built/Soil', track '1'.
## From To Area
## 1 Forest Forest 5.19175
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## 21 Built/Soil Forest 0.21475
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## 21 Built/Soil Forest 0.21475
## 22 Built/Soil Agri 0.08650
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## 21 Built/Soil Forest 0.21475
## 22 Built/Soil Agri 0.08650
## 23 Built/Soil Grass 0.00875
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## 21 Built/Soil Forest 0.21475
## 22 Built/Soil Agri 0.08650
## 23 Built/Soil Grass 0.00875
## 24 Built/Soil W 0.02675
## From To Area
## 1 Forest Forest 5.19175
## 2 Forest Agri 0.49225
## 3 Forest Grass 0.09325
## 4 Forest W 0.01175
## 5 Forest Built/Soil 0.02475
## 6 Agri Forest 1.05675
## 7 Agri Agri 0.45125
## 8 Agri Grass 0.05200
## 9 Agri W 0.02375
## 10 Agri Built/Soil 0.03250
## 11 Grass Forest 0.01700
## 12 Grass Agri 0.00125
## 13 Grass Grass 0.00700
## 14 Grass W 0.00000
## 15 Grass Built/Soil 0.00000
## 16 W Forest 0.02975
## 17 W Agri 0.02650
## 18 W Grass 0.00150
## 19 W W 0.01400
## 20 W Built/Soil 0.00850
## 21 Built/Soil Forest 0.21475
## 22 Built/Soil Agri 0.08650
## 23 Built/Soil Grass 0.00875
## 24 Built/Soil W 0.02675
## 25 Built/Soil Built/Soil 0.01475
## Note: 1 point is out of plotting region in sector 'Forest', track '1'.
## Note: 6 points are out of plotting region in sector 'Forest', track '1'.
## Note: 1 point is out of plotting region in sector 'Agri', track '1'.
## Note: 4 points are out of plotting region in sector 'Agri', track '1'.
## Note: 1 point is out of plotting region in sector 'Grass', track '1'.
## Note: 5 points are out of plotting region in sector 'Grass', track '1'.
## Note: 1 point is out of plotting region in sector 'W', track '1'.
##
## Note: 1 point is out of plotting region in sector 'W', track '1'.
## Note: 1 point is out of plotting region in sector 'Built/Soil', track '1'.
## Note: 10 points are out of plotting region in sector 'Built/Soil', track '1'.
## From To Area
## 1 Forest Forest 4.37425
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## 21 Built/Soil Forest 0.12450
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## 21 Built/Soil Forest 0.12450
## 22 Built/Soil Agri 0.47050
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## 21 Built/Soil Forest 0.12450
## 22 Built/Soil Agri 0.47050
## 23 Built/Soil Grass 0.07875
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## 21 Built/Soil Forest 0.12450
## 22 Built/Soil Agri 0.47050
## 23 Built/Soil Grass 0.07875
## 24 Built/Soil W 0.00975
## From To Area
## 1 Forest Forest 4.37425
## 2 Forest Agri 0.54325
## 3 Forest Grass 0.26275
## 4 Forest W 0.00300
## 5 Forest Built/Soil 0.04000
## 6 Agri Forest 0.64250
## 7 Agri Agri 1.47400
## 8 Agri Grass 0.31825
## 9 Agri W 0.04450
## 10 Agri Built/Soil 0.34825
## 11 Grass Forest 0.05500
## 12 Grass Agri 0.12400
## 13 Grass Grass 0.02850
## 14 Grass W 0.00200
## 15 Grass Built/Soil 0.02700
## 16 W Forest 0.05225
## 17 W Agri 0.22875
## 18 W Grass 0.01925
## 19 W W 0.03425
## 20 W Built/Soil 0.07200
## 21 Built/Soil Forest 0.12450
## 22 Built/Soil Agri 0.47050
## 23 Built/Soil Grass 0.07875
## 24 Built/Soil W 0.00975
## 25 Built/Soil Built/Soil 0.29850
## Note: 1 point is out of plotting region in sector 'Forest', track '1'.
## Note: 6 points are out of plotting region in sector 'Forest', track '1'.
## Note: 1 point is out of plotting region in sector 'Agri', track '1'.
## Note: 4 points are out of plotting region in sector 'Agri', track '1'.
## Note: 1 point is out of plotting region in sector 'Grass', track '1'.
## Note: 5 points are out of plotting region in sector 'Grass', track '1'.
## Note: 1 point is out of plotting region in sector 'W', track '1'.
##
## Note: 1 point is out of plotting region in sector 'W', track '1'.
## Note: 1 point is out of plotting region in sector 'Built/Soil', track '1'.
## Note: 10 points are out of plotting region in sector 'Built/Soil', track '1'.

#dev.off()
# show the plot
#plot(region_mask, ylim=c(21.5,25.5))
#par(new=FALSE)