R/read_dima.R

In mschmidty/AIMtools: Tools for working with AIM Data.

read_dima<-function(file_path){

  read_table<-function(connection, name) tibble::as_tibble(RODBC::sqlFetch(connection, name))

  if(Sys.getenv("R_ARCH") != "/i386"){

    stop('You are currently using a 64-bit version of R. Please change your R version to 32-bit in "Tools>Global Options" and in the dropdown change your R version to 32-bit. Microsoft Access Does not work with 64-bit R.')
  }
  con<-RODBC::odbcConnectAccess(file_path)

  data<-list(
    plots = read_table(con, "tblPlots"),
    lines = read_table(con, "tblLines"),
    lpi_detail = read_table(con, "tblLPIDetail")%>%
      dplyr::mutate_all(dplyr::na_if,""),
    lpi_header = read_table(con, "tblLPIHeader"),
    sp_rich_detail = read_table(con, "tblSpecRichDetail"),
    sp_rich_header = read_table(con, "tblSpecRichHeader"),
    gap_detail = read_table(con, "tblGapDetail"),
    gap_header = read_table(con, "tblGapHeader"),
    unknown_plants = read_table(con,"UnknownTracking")
  )

  RODBC::odbcClose(con)

  data[["join_table"]]<-dplyr::select(data$lines, PlotKey, LineKey, LineID)%>%
    dplyr::mutate(PlotKey = as.character(PlotKey))%>%
    dplyr::left_join(dplyr::select(data$plots, PlotKey, PlotID)%>%
                       dplyr::mutate(PlotKey = as.character(PlotKey))
                     , by="PlotKey")

  data[["plant_list"]]<-plant_list%>%
    dplyr::select(list_name:fulllist)

  data[["spatial"]]<-data$plots%>%
    dplyr::select(SiteKey, PlotKey, PlotID, Easting, Northing)%>%
    dplyr::filter(!is.na(Easting))%>%
    sf::st_as_sf(coords = c("Easting", "Northing"), crs = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")%>%
    dplyr::filter(!stringr::str_detect(SiteKey, "^888|^999"))

  data
}