## ----eval = FALSE-------------------------------------------------------------
#  data <- dsmSearch::get_dsm_30(bbox = c(-83.783557,42.241833,-83.696525,42.310420), key = "Your KEY")

## ----eval = FALSE-------------------------------------------------------------
#  # search for lidar data information using bbox
#  search_result <- dsmSearch::lidar_search(bbox = c(-83.742282,
#                                                    42.273389,
#                                                    -83.733442,
#                                                    42.278724),
#                                           preview = FALSE)
#  search_result

## ----eval = FALSE-------------------------------------------------------------
#  # try coordinates -83.741289,42.270146 (in south Michigan, USA)
#  # radius is 1000ft
#  las <- dsmSearch::get_lidar(x = -83.741289,
#                              y = 42.270146,
#                              r = 1000,
#                              epsg = 2253)
#  # download with bbox
#  las <- dsmSearch::get_lidar(bbox = c(-83.742282,42.273389,-83.733442,42.278724),
#                              epsg = 2253)

## ----eval = FALSE-------------------------------------------------------------
#  # Create DTM
#  dtm_ <- lidR::rasterize_terrain(las, res = 5, lidR::tin())
#  terra::plot(dtm_)

## ----eval = FALSE-------------------------------------------------------------
#  # Create DSM
#  dsm_ <- lidR::rasterize_canopy(las, res = 5, lidR::dsmtin())
#  raster::plot(dsm_)