rm(list=ls()) library(raster) library(sf) library(rasterVis) library(tidyverse) library(haven) library(stringr) library(scales)
dat %>% ggplot(aes(nl_viirs_1, nl_viirs_2)) + geom_point() + theme_bw()