def geo_correction(date_pair, dem_header, dem, incidence_angle_or_map):
""" using geo coordinates to remove APS """
aps1 = pa.PyAPSPyRateGeo(os.path.join(ECMWF_DIR, ECMWF_PRE + date_pair[0] +
ECMWF_EXT),
dem_header=dem_header,
dem=dem,
grib=ECMWF,
verb=True)
aps2 = pa.PyAPSPyRateGeo(os.path.join(ECMWF_DIR, ECMWF_PRE + date_pair[1] +
ECMWF_EXT),
dem_header=dem_header,
dem=dem,
grib=ECMWF,
verb=True)
phs1 = np.zeros((aps1.ny, aps1.nx))
phs2 = np.zeros((aps2.ny, aps2.nx))
print('Without Lat Lon files')
aps1.getdelay_pyrate(phs1, dem, inc=incidence_angle_or_map)
aps2.getdelay_pyrate(phs2, dem, inc=incidence_angle_or_map)
aps_delay = phs2 - phs1 # delay in meters as we don't provide wavelength
return aps_delay
