Extracting raster values from list of x,y points

This script extracts the value of the closest pixel in a raster to a list of x,y points.

The script assumes the following:

• Input file with list of latitudes (WGS-84, EPSG: 4326)

• Input file with list of longitudes (WGS-84, EPSG: 4326)  

• Image file already in WGS-84 projection

The script requires the following packages to be installed in your Python environment:

• numpy

• pandas

• geopandas

• rasterio

• glob

• shapely