import rasterio
import pandas as pd

df = pd.read_csv("ATT_Mobility_Network_Definition_v7.csv", dtype=str)
row = df.iloc[0]

tif_path = r"c:\users\public\clipped_dems\WQSY442_WQST669.tif"

lon = float(row["A_Longitude"])
lat = float(row["A_Latitude"])
print("lon, lat:", lon, lat)

with rasterio.open(tif_path) as src:
    print("CRS:", src.crs)
    print("bounds:", src.bounds)
    print("nodata:", src.nodata)
    val = list(src.sample([(lon, lat)]))[0][0]
    print("sampled value:", val, type(val))
