using CSV using Rasters using NighttimeLights using Shapefile using Plots using DataFrames capex = CSV.read("completed_2013_2021_projects_lat_long.csv", DataFrame) world_shapefile = Shapefile.Table("../centre_of_mass/world_shapefile/in_primary_countries.shp") indiapolygon = world_shapefile.geometry[78] p = plot(indiapolygon, color = "white") for i in capex.lat_long pi = parse.(Float64, split(i, ",")) if pi[2] <50 continue end scatter!(p, [pi[2]], [pi[1]],color = :red, markerstrokewidth=0, legend = false, markersize = 0.3) end p