using Rasters
using NighttimeLights
using Plots

lat, long = 21.963683,83.8912462

R1, R2 = 1.5, 10
res = 15

img = Raster("/home/nighttimelights/VIIRS_ANNUAL/v21_median/VNL_v21_npp_2014_global_vcmslcfg_c202205302300.median_masked.dat.tif", lazy = true)
bounds, m = annular_ring(img, R1, R2, lat, long, res) 
plot(img[bounds...])
plot(img[bounds...] .* m)
xlabel!("Longitude")
ylabel!("Latitude")

img = Raster("/home/nighttimelights/VIIRS_ANNUAL/v21_median/VNL_v21_npp_2018_global_vcmslcfg_c202205302300.median_masked.dat.tif", lazy = true)
bounds, m = annular_ring(img, R1, R2, lat, long, res) 
plot(img[bounds...])
plot(img[bounds...] .* m)
xlabel!("Longitude")
ylabel!("Latitude")