def station(dat): if dat==0: #Robledo longi=4+(14/60)+(53/3600) #W (-1) lati=40+(25/60)+(54/3600) #N longi= -longi*(np.pi/180) lati=lati*(np.pi/180) alti=720/1000 #Va con Km # Datos de WGS84: #===================================== RE= 6378137.0/1000#Semieje mayor f=1/298.257223563 #Achatamiento Tierra #===================================== return spi.georec(longi,lati,alti,RE,f)