mbjoseph

mbjoseph

library(rgdal)
library(rgeos)
library(raster)
library(gdalUtils)


# gdalinfo("GFED40_MQ_199506_BA.hdf")
# sds <- get_subdatasets("GFED40_MQ_199506_BA.hdf")
# sds

mbjoseph

# Importing data from USDA into large data file ---------------------------
library(stringr)
library(tidyverse)

convert_to_df <- function(file) {
  # takes a text file as input and parses the file to generate a data.frame
  lines <- readLines(file)

  # identify blank lines which denote separate files

mbjoseph

library(raster)
library(tidyverse)
library(lubridate)
library(mgcv)

# Boulder temperature prediction ------------------------------------------
# Goal: use historical temperature data to predict the weather a week from today

# Fetch raw data from NOAA
data_url <- "https://www.esrl.noaa.gov/psd/boulder/data/boulderdaily.complete"

mbjoseph

print("Hello, world!")

mbjoseph

library(RCurl)

# Fetching burn rasters ---------------------------------------------------
# get and read the file listing from the usgs server
prefix  <- 'https://rmgsc.cr.usgs.gov/outgoing/baecv/BAECV_CONUS_v1_2017/'
download.file(prefix, destfile = 'listing.txt')
out <- readLines('listing.txt')

# now create paths to each tar.gz file by processing the html

mbjoseph

Keybase proof

I hereby claim:

• I am mbjoseph on github.
• I am mbjoseph (https://keybase.io/mbjoseph) on keybase.
• I have a public key ASABiVndWu_-6CJIwyKMRAlsmF3sEsWB33XMDPoZ2GAbvAo

To claim this, I am signing this object:

mbjoseph

## Simple bivariate Poisson model in stan
## following parameterization in Karlis and Ntzoufras 2003
## Simulate data
n <- 50

# indpendent Poisson components
theta <- c(2, 3, 1)

X_i <- array(dim=c(n, 3))
for (i in 1:3){

mbjoseph

---
title: 'Mapping the 2018 East Africa floods from space with smapr'
author: "Max Joseph"
date: "2018-09-14"
output:
  html_document:
    keep_md: yes
---

```{r setup, include=FALSE}

mbjoseph

library(tidyverse)
library(raster)
library(sf)
library(rmapshaper)
library(elevatr)
library(USAboundaries)

dir.create('out')

# RGB satellite imagery over Rocky Mountain National Park -------------