Adding a function plot_img

src/NighttimeLights.jl

@@ -1,11 +1,11 @@
 module NighttimeLights
 
-using Dates, ArchGDAL, GeoArrays, DataFrames, Shapefile, GeoInterface, ProgressMeter, StatsBase, Statistics, SmoothingSplines, CubicSplines, GLM, JLD, Distributions, HypothesisTests, Plots
+using Dates, ArchGDAL, GeoArrays, DataFrames, Shapefile, GeoInterface, ProgressMeter, StatsBase, Statistics, SmoothingSplines, CubicSplines, GLM, JLD, Distributions, HypothesisTests, Plots, Images
 
 export load_img, load_datacube, save_img, save_datacube, 
 lat_to_row, row_to_lat, lat_to_row, long_to_column, column_to_long, translate_geometry, Coordinate, CoordinateSystem, image_to_coordinate,
 polygon_mask, load_shapefile, mask_area, long_apply, cross_apply, apply_mask, view_img, aggregate, aggregate_timeseries, aggregate_dataframe, bias_correction, bias_correction_datacube, outlier_mask, outlier_ts, linear_interpolation,
-conventional_cleaning, PatnaikSTT2021, threshold_datacube, background_noise_mask, sparse_cube, mark_missing, make_datacube, load_example, radiance_datacube, clouds_datacube, mumbai_map, TILE1_COORDINATE_SYSTEM, TILE2_COORDINATE_SYSTEM, TILE3_COORDINATE_SYSTEM, TILE4_COORDINATE_SYSTEM, TILE5_COORDINATE_SYSTEM, TILE6_COORDINATE_SYSTEM, FULL_COORDINATE_SYSTEM_COORDINATE_SYSTEM, INDIA_COORDINATE_SYSTEM, MUMBAI_COORDINATE_SYSTEM, coordinate_to_image, aggregate_per_area_dataframe, replace_negative, bounding_box, plot_chloropleth
+conventional_cleaning, PatnaikSTT2021, threshold_datacube, background_noise_mask, sparse_cube, mark_missing, make_datacube, load_example, radiance_datacube, clouds_datacube, mumbai_map, TILE1_COORDINATE_SYSTEM, TILE2_COORDINATE_SYSTEM, TILE3_COORDINATE_SYSTEM, TILE4_COORDINATE_SYSTEM, TILE5_COORDINATE_SYSTEM, TILE6_COORDINATE_SYSTEM, FULL_COORDINATE_SYSTEM_COORDINATE_SYSTEM, INDIA_COORDINATE_SYSTEM, MUMBAI_COORDINATE_SYSTEM, coordinate_to_image, aggregate_per_area_dataframe, replace_negative, bounding_box, plot_chloropleth, plot_img
 
 include("coordinate_system.jl")
 include("polygons.jl")

src/plot.jl

@@ -39,4 +39,37 @@ plot_chloropleth(mumbai_map, :censuscode, :DISTRICT)
 """
 function plot_chloropleth(geoms::DataFrames.DataFrame, data::Symbol, names::Symbol)
  plot_chloropleth(geoms, geoms[!,string(data)], names)
-end
+end
+
+"""
+load_example()
+april2012=radiance_datacube[:,:,1]
+plot_img(april2012, MUMBAI_COORDINATE_SYSTEM)
+"""
+function plot_img(img, coordinate_system)
+
+ function lat(c, a, b)
+ return round(image_to_coordinate(c, [a, b]).latitude; digits=2)
+ end
+
+ function lon(c, a, b)
+ return round(image_to_coordinate(c, [a, b]).longitude; digits=2)
+ end
+
+ a=size(img)[1]/6
+ b=size(img)[2]/5
+ c=coordinate_system
+
+ xc=[b, 2*b, 3*b, 4*b]
+ yc=[a, 2*a, 3*a, 4*a, 5*a]
+ xcs=[lon(c,1,b), lon(c,1,2*b), lon(c,1,3*b), lon(c,1,4*b)]
+ ycs=[lat(c,a,1), lat(c,2*a,1), lat(c,3*a,1), lat(c,4*a,1), lat(c,5*a,1)]
+
+ Plots.plot(Gray.(img))
+ plot!(xc, seriestype="vline", xticks = (xc,xcs), label="", color=:red, linestyle=:dot)
+ plot!(yc, seriestype="hline", yticks = (yc,ycs), label="", color=:red, linestyle=:dot)
+ xlabel!("longitude")
+ ylabel!("latitude")
+end
+
+