Merge pull request #24 from Tumon2001/Plots

Adding functionality to plot datacubes

src/NighttimeLights.jl

@@ -5,7 +5,7 @@ using Dates, ArchGDAL, GeoArrays, DataFrames, Shapefile, GeoInterface, ProgressM
 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, plot_img
+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, plot_datacube
 
 include("coordinate_system.jl")
 include("polygons.jl")

src/plot.jl

@@ -72,4 +72,50 @@ function plot_img(img, coordinate_system)
  ylabel!("latitude")
 end
 
+"""
+using Dates
+load_example()
+dates = collect(Date(2012,4):Month(1):Date(2020, 02))
+plot_datacube(radiance_datacube, MUMBAI_COORDINATE_SYSTEM, string.(dates))
+"""
+function plot_datacube(datacube, coordinate_system, date)
+ map=[]
 
+ function plot_img(img, coordinate_system, i)
+
+ 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)]
+
+ p=plot(Gray.(img), title=date[i], xlabel="longitude", ylabel="latitude")
+ plot!(xc, seriestype=:vline, xticks = (xc,xcs), label="", color=:red, linestyle=:dot)
+ plot!(yc, seriestype=:hline, yticks = (yc,ycs), label="", color=:red, linestyle=:dot)
+
+ append!(map,[p])
+ end
+
+ for i in 1:size(datacube)[3]
+ img=datacube[:, :, i]
+ plot_img(img, coordinate_system, i)
+ end
+
+ Plots.GR.beginprint("mumbai.pdf")
+ gr(show=true)
+ for i in 1:size(datacube)[3]
+ plot(map[i])
+ end
+ Plots.GR.endprint()
+end