analysis/5a_get-nj-noaa-images.R

In treysp/coolit.train: Train model to identify water cooling towers from aerial orthoimagery

library(stringr)
library(sf)

index <- st_read("data/source_from-nj-website/nj-images/0tileindex")

counties <- tigris::counties("nj")
counties <- st_as_sf(counties)
counties <- counties[counties$NAME == "Union",]

index <- st_transform(index, st_crs(counties))

temp_int <- st_intersects(index, counties)

nj_images <- index[sapply(temp_int, function(x) length(x) > 0),]
nj_images$file_name <- str_match(nj_images$URL, "(.*/)*(.*)\\.tif$")[, 3]
nj_images$url_base <- str_match(nj_images$URL, "(.*/)*(.*)\\.tif$")[, 2]

dest_dir <- "data/source_from-nj-website/nj-images"

for (i in seq_len(nrow(nj_images))) {
  Sys.sleep(2)

  download.file(
    paste0(nj_images$url_base[i], nj_images$file_name[i], ".tif"),
    file.path(dest_dir, paste0(nj_images$file_name[i], ".tif")),
    mode = "wb",
    quiet = TRUE
  )

  download.file(
    paste0(nj_images$url_base[i], nj_images$file_name[i], ".tfw"),
    file.path(dest_dir, paste0(nj_images$file_name[i], ".tfw")),
    quiet = TRUE
  )
}