# boundaries for various plots in R compiled by Ian Taylor
# feel free to share but realize that it could all be wrong
# don't remember the source for much of this
# if anybody is aware of any mistakes in these boundaries, tell me
library(maps)
library(mapdata)

USboundary = matrix(c(
-115, 49.00224722,
-122.76, 49.00224722,
-123.3209222, 49.00224722,
-123.0071833, 48.83139167,
-123.0071861, 48.76727222,
-123.2665722, 48.69417222,
-123.2176306, 48.54888889,
-123.1585806, 48.45368889,
-123.1138472, 48.42301111,
-123.2471444, 48.2842,
-123.5399167, 48.22474167,
-123.6777583, 48.24018611,
-124.0107361, 48.29686667,
-124.7258889, 48.49365833,
-124.7258861, 48.49366389,
-124.7869444, 48.50305556,
-124.8391667, 48.50611111,
-124.9144444, 48.50388889,
-124.9872222, 48.49916667,
-125.0016667, 48.49555556,
-125.0963889, 48.46916667,
-125.1402778, 48.45277778,
-125.1533333, 48.44638889,
-125.38, 48.33777778,
-125.4994444, 48.30611111,
-125.8966667, 48.18472222,
-126.6825, 47.82083333,
-127.1994444, 47.61305556,
-127.6897222, 47.36666667,
-128.8655556, 46.70138889,
-129.1275, 46.52972222,
-129, 32.5,
-115, 32.5,
-115, 49.00224722
),34,2,byrow=T)

#Bonilla-Tatoosh Line: A line running from the western end of Cape Flattery to Tatoosh Island Lighthouse to the buoy adjacent to Duntze Rock, then in a straight line to Bonilla Point on Vancouver Island.
#
BTline = matrix(c(
-124-44/60-12/3600, 48+23/60+30/3600, #Tatoosh Lighthouse
-124-45/60-00/3600, 48+25/60+00/3600, #Duntze Rock   !!website says 48 28' but maps looks like 48 25' http://www.pcouncil.org/salmon/salpreIII05/PreIIITable2040805.pdf
-124-43/60-00/3600, 48+35/60+30/3600  #Bonilla Point
),3,2,byrow=T)

#calculating slope of line between Bonilla Pt and Duntze Rock
m1 = (48+35/60+30/3600 - (48+25/60+00/3600))/(-124-43/60-00/3600 - (-124-45/60-00/3600))
#calculating slope of line it crosses on the border
m2 = (48.49365833 - 48.29686667)/(-124.7258889 - -124.0107361)
#calculating the x coordinate of the intersection
#(y1 - y2 + m2x2 - m1x1)/(m2-m1)
x = (48+35/60+30/3600 - 48.49365833 + m2*-124.7258889 - m1*(-124-43/60-00/3600))/(m2-m1)
y = m1*(x-(-124-43/60-00/3600)) + 48+35/60+30/3600
#> x
#[1] -124.7349
#> y
#[1] 48.49613


PSboundary = matrix(c(
-122, 49.00224722,
-122.76, 49.00224722,
-123.3209222, 49.00224722,
-123.0071833, 48.83139167,
-123.0071861, 48.76727222,
-123.2665722, 48.69417222,
-123.2176306, 48.54888889,
-123.1585806, 48.45368889,
-123.1138472, 48.42301111,
-123.2471444, 48.2842,
-123.5399167, 48.22474167,
-123.6777583, 48.24018611,
-124.0107361, 48.29686667,
-124.7349, 48.49613,  #the corner
-124-45/60-00/3600, 48+25/60+00/3600, #Duntze Rock   !!website says 48 28' but maps looks like 48 25' http://www.pcouncil.org/salmon/salpreIII05/PreIIITable2040805.pdf
-124-44/60-12/3600, 48+23/60+30/3600, #Tatoosh Lighthouse
-124-43.75/60, 48+23/60, #Cape Flattery
-123, 47,
-122, 47,
-122, 49.00224722
),20,2,byrow=T)


GSboundary = matrix(c(
-122, 49.00224722,
-122.76, 49.00224722,
-123.3209222, 49.00224722,
-123.0071833, 48.83139167,
-123.0071861, 48.76727222,
-123.2665722, 48.69417222,
-123.2176306, 48.54888889,
-123.1585806, 48.45368889,
-123.1138472, 48.42301111,
-123.2471444, 48.2842,
-123.5399167, 48.22474167,
-123.6777583, 48.24018611,
-124.0107361, 48.29686667,
-124.7349, 48.49613,  #the corner
-124-43/60-00/3600, 48+35/60+30/3600,  #Bonilla Point
-124-43/60-00/3600, 49.2,
-125.5, 50.3,
-127, 50.5,
-128.09,50.86,
-127.79,51.17,
-125, 51.5,
-122, 49.00224722
),22,2,byrow=T)

CANboundary = matrix(c(
-122, 49.00224722,
-122.76, 49.00224722,
-123.3209222, 49.00224722,
-123.0071833, 48.83139167,
-123.0071861, 48.76727222,
-123.2665722, 48.69417222,
-123.2176306, 48.54888889,
-123.1585806, 48.45368889,
-123.1138472, 48.42301111,
-123.2471444, 48.2842,
-123.5399167, 48.22474167,
-123.6777583, 48.24018611,
-124.0107361, 48.29686667,
-124.7258889, 48.49365833,
-124.7258861, 48.49366389,
-124.7869444, 48.50305556,
-124.8391667, 48.50611111,
-124.9144444, 48.50388889,
-124.9872222, 48.49916667,
-125.0016667, 48.49555556,
-125.0963889, 48.46916667,
-125.1402778, 48.45277778,
-125.1533333, 48.44638889,
-125.38, 48.33777778,
-125.4994444, 48.30611111,
-125.8966667, 48.18472222,
-126.6825, 47.82083333,
-127.1994444, 47.61305556,
-127.6897222, 47.36666667,
-128.8655556, 46.70138889,
-129.1275, 46.52972222,
-135, 54+40/60,
-128, 54+40/60,
-122, 49.00224722
),34,2,byrow=T)

NBCboundary = cbind(c(-135,-120,-120,-135,-135),c(50.4,50.4,55,55,50.4))