https://github.com/smroecker/mast-mojave

Soil Temperature in the Mojave Desert
https://github.com/smroecker/mast-mojave

mojave-desert soil-survey soil-temperature

Last synced: 7 days ago
JSON representation

Soil Temperature in the Mojave Desert

• Host: GitHub
• URL: https://github.com/smroecker/mast-mojave
• Owner: smroecker
• License: gpl-2.0
• Created: 2019-03-24T19:30:44.000Z (over 5 years ago)
• Default Branch: master
• Last Pushed: 2019-03-26T18:28:37.000Z (over 5 years ago)
• Last Synced: 2024-10-28T05:17:17.248Z (about 2 months ago)
• Topics: mojave-desert, soil-survey, soil-temperature
• Language: HTML
• Homepage: https://www.researchgate.net/publication/331977150_Modeling_and_Application_of_Soil_Temperature_in_the_Mojave_and_Lower_Colorado_Deserts_of_California/citations
• Size: 1.57 MB
• Stars: 0
• Watchers: 0
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.Rmd
  • License: LICENSE

Awesome Lists containing this project

README

        
---

output: 

  html_document: 

    keep_md: yes

---

# Code from the MLRA 30 and 31 - Soil Climate Study

- [Project Description](https://github.com/smroecker/mast-mojave/blob/master/projectdescription.md)

- [Analysis](https://smroecker.github.io/mast-mojave/analysis_stephen.html)

- [Poster](https://www.researchgate.net/publication/331977150_Modeling_and_Application_of_Soil_Temperature_in_the_Mojave_and_Lower_Colorado_Deserts_of_California/citations)

```{r ggmap, echo=FALSE, message=FALSE, warning=FALSE}

library(ggplot2)

library(ggmap)

library(sf)

# mast sites

sites  <- read.csv("D:/projects/soilTemperatureMonitoring/data/R/HOBO_List_2013_0923_master.csv") %>%

  st_as_sf(coords = c("utmeasting", "utmnorthing"), dim = "XY")

st_crs(sites) <- "+init=epsg:26911"

sites <- st_transform(sites, "+init=epsg:4326")

xy <- as.data.frame(st_coordinates(sites))

names(xy) <- c("lon", "lat")

# mlra

mlra <- read_sf(dsn = "D:/geodata/soils/mlra_a_mbr.shp", layer = "mlra_a_mbr")

st_crs(mlra) <- "+init=epsg:4326"

mojave <- subset(mlra, MLRARSYM %in% c("30", 31))

bb <- make_bbox(lat = xy$lat, lon = xy$lon)

gmap <- get_map(bb)

ggmap(gmap) +

  coord_sf(crs = st_crs(4326)) +

  geom_sf(data = mojave, fill = NA, inherit.aes = FALSE) +

  geom_point(data = xy, aes(x = lon, y = lat)) +

  ggtitle("Location of MAST sites")

```