I have some shapefiles with lat/lon coordinates and I want to extract elevation using the USGS-National Map. I have created a python script following this question. It works