def test_parse_full_roipac_header(self):
# Ensures "long style" original header can be parsed correctly
hdrs = roipac.parse_header(FULL_HEADER_PATH)
# check DATE/ DATE12 fields are parsed correctly
date0 = date(2006, 6, 19) # from "DATE 060619" header
date2 = date(2006, 8, 28) # from DATE12 060619-060828
self.assertEqual(hdrs[ifc.MASTER_DATE], date0)
self.assertEqual(hdrs[ifc.SLAVE_DATE], date2)
def test_parse_short_roipac_header(self):
hdrs = roipac.parse_header(SHORT_HEADER_PATH)
self.assertEqual(hdrs[ifc.PYRATE_NCOLS], 47)
self.assertEqual(hdrs[ifc.PYRATE_NROWS], 72)
self.assertAlmostEqual(hdrs[ifc.PYRATE_LONG], 150.910)
self.assertEqual(hdrs[ifc.PYRATE_X_STEP], 0.000833333)
self.assertEqual(hdrs[ifc.PYRATE_LAT], -34.170000000)
self.assertEqual(hdrs[ifc.PYRATE_Y_STEP], -0.000833333)
self.assertEqual(hdrs[ifc.PYRATE_WAVELENGTH_METRES], 0.0562356424)
def test_parse_short_header_has_timespan(self):
# Ensures TIME_SPAN_YEAR field is added during parsing
hdrs = roipac.parse_header(SHORT_HEADER_PATH)
self.assertIn(roipac.TIME_SPAN_YEAR, hdrs.keys())
# check time span calc
master = date(2006, 6, 19)
slave = date(2006, 10, 2)
diff = (slave - master).days / ifc.DAYS_PER_YEAR
self.assertEqual(diff, hdrs[roipac.TIME_SPAN_YEAR])
def test_to_geotiff_dem(self):
hdr = roipac.parse_header(SML_TEST_DEM_HDR)
self.dest = os.path.join(TEMPDIR, "tmp_roipac_dem.tif")
write_fullres_geotiff(hdr, SML_TEST_DEM_ROIPAC, self.dest, nodata=0)
exp_path = join(SML_TEST_DEM_DIR, 'roipac_test_trimmed.tif')
exp_ds = gdal.Open(exp_path)
ds = gdal.Open(self.dest)
# compare data and geographic headers
assert_array_almost_equal(exp_ds.ReadAsArray(), ds.ReadAsArray())
self.compare_rasters(ds, exp_ds)
self.assertIsNotNone(ds.GetMetadata())
def setUpClass(cls):
hdr_path = join(PREP_TEST_OBS, 'geo_060619-061002.unw.rsc')
cls.HDRS = roipac.parse_header(hdr_path)
def test_xylast(self):
# Test the X_LAST and Y_LAST header elements are calculated
hdrs = roipac.parse_header(FULL_HEADER_PATH)
self.assertAlmostEqual(hdrs[roipac.X_LAST], 151.8519444445)
self.assertAlmostEqual(hdrs[roipac.Y_LAST], -34.625)
def test_read_full_roipac_header2(self):
# Tests header from cropped original dataset is parsed correctly
hdrs = roipac.parse_header(FULL_HEADER_PATH)
self.assertTrue(len(hdrs) is not None)
