While there is a function julian_date() that will go the other direction, we will have to do a bit of work ourselves to go the other direction by figuring out what ephem's name is for. Free demo version that allows the user to save the inversion results for 2D data sets (including topography) with up to 84 electrodes. The elevation angle is 90° minus the inclination angle. Specify the origin of the local AER system with the geodetic coordinates lat0, lon0, and h0. Point A 700 Just follow along the index contour from point A until you find a marked elevation. Using the Fig.2 we estimate the elevation of the points marked with letters A,B,C and D. Second, it appears that PyEphem has forgotten to provide an easy way to convert dates from Julian to its own format. The spatial reference includes a coordinate system for x-, y-, and z-values as well as tolerance and resolution values for x-, y-, z-, and m-values. lat,lon,h aer2geodetic (az,elev,slantRange,lat0,lon0,h0,spheroid) transforms the local azimuth-elevation-range (AER) spherical coordinates specified by az, elev, and slantRange to the geodetic coordinates specified by lat, lon, and h. Once we know how to determine the elevation of the unmarked contour lines, we should be able determine or at least estimate the elevation of any point on the map. First, you have happened to use “elevation” and “altitude” to mean the opposite of what the two terms mean in the PyEphem library - so you are calling the spot in the sky its “elevation / azimuth' position instead of its “altitude / azimuth” position.