In iGISc/igisci: Data and Code for Environmental Data Science

library(sf); library(ggplot2)
CO <- st_polygon(list(rbind(c(-109,41),c(-102,41),c(-102,37),c(-109,37),c(-109,41))))
WY <- st_polygon(list(rbind(c(-111,45),c(-104,45),c(-104,41),c(-111,41),c(-111,45))))
UT <- st_polygon(list(rbind(c(-114,42),c(-111,42),c(-111,41),c(-109,41),c(-109,37),
                            c(-114,37),c(-114,42))))
AZ <- st_polygon(list(rbind(c(-114,37),c(-109,37),c(-109,31.3),c(-111,31.3),
                            c(-114.8,32.5),c(-114.6,32.7),c(-114.1,34.3),c(-114.5,35),
                            c(-114.5,36),c(-114,36),c(-114,37))))
NM <- st_polygon(list(rbind(c(-109,37),c(-103,37),c(-103,32),c(-106.6,32),c(-106.5,31.8),
                            c(-108.2,31.8),c(-108.2,31.3),c(-109,31.3),c(-109,37))))

CA <- st_polygon(list(rbind(c(-124,42),c(-120,42),c(-120,39),c(-114.5,35),
                            c(-114.1,34.3),c(-114.6,32.7),c(-117,32.5),c(-118.5,34),c(-120.5,34.5),
                            c(-122,36.5),c(-121.8,36.8),c(-122,37),c(-122.4,37.3),c(-122.5,37.8),
                            c(-123,38),c(-123.7,39),c(-124,40),c(-124.4,40.5),c(-124,41),c(-124,42))))
NV <- st_polygon(list(rbind(c(-120,42),c(-114,42),c(-114,36),c(-114.5,36),
                            c(-114.5,35),c(-120,39),c(-120,42))))

OR <- st_polygon(list(rbind(c(-124,42),c(-124.5,43),
                            c(-124,46),c(-123,46),c(-122.7,45.5),c(-119,46),c(-117,46),
                            c(-116.5,45.5),c(-117.2,44.5),c(-117,44),
                            c(-117,42),c(-120,42),c(-124,42))))
WA <- st_polygon(list(rbind(c(-124,46),c(-124.8,48.4),c(-123,48),
                            c(-123,49),c(-117,49),
                            c(-117,46),c(-119,46),c(-122.7,45.5),c(-123,46),c(-124,46))))
ID <- st_polygon(list(rbind(c(-117,49),
                            c(-116,49),c(-116,48),c(-114.4,46.5),c(-114.4,45.5),
                            c(-114,45.6),c(-113,44.5),c(-111,44.5),
                            c(-111,42),c(-114,42),c(-117,42),
                            c(-117,44),c(-117.2,44.5),c(-116.5,45.5),
                            c(-117,46),c(-117,49))))
MT <- st_polygon(list(rbind(c(-116,49),c(-104,49),
                            c(-104,45),c(-111,45),
                            c(-111,44.5),c(-113,44.5),c(-114,45.6),c(-114.4,45.5),
                            c(-114.4,46.5),c(-116,48),c(-116,49))))

iGISc/igisci documentation built on June 6, 2024, 6:57 p.m.