def test_calculation_of_tie_point_grid_float_coords(self): # NOTE: This does not test against unequaly sized output of get_image_windows_to_match(). # overwrite gt and prj ref = GeoArray(self.ref_path) ref.to_mem() ref.filePath = None tgt = GeoArray(self.tgt_path) tgt.to_mem() tgt.filePath = None ref.gt = [ 330000.19999996503, 10.00000001, 0.0, 5862000.7999997628, 0.0, -10.00000001 ] # ref.gt = [330000.1, 10.1, 0.0, 5862000.1, 0.0, -10.1] tgt.gt = [ 335440.19999996503, 10.00000001, 0.0, 5866490.7999997628, 0.0, -10.00000001 ] # tgt.gt = [330000.1, 10.1, 0.0, 5862000.1, 0.0, -10.1] # get instance of COREG_LOCAL object CRL = COREG_LOCAL(ref, tgt, **dict(CPUs=32, **self.coreg_kwargs)) CRL.calculate_spatial_shifts()
def test_shift_calculation_noWGS84(self): """Test if shift computation properly raises a RunTimeError if the matching window is centered at a cloudy image position. """ ref = GeoArray(self.ref_path).to_mem() tgt = GeoArray(self.tgt_path).to_mem() # force to overwrite projection ref.filePath = None tgt.filePath = None ref.prj = EPSG2WKT(3035) # ETRS89_LAEA_Europe tgt.prj = EPSG2WKT(3035) # ETRS89_LAEA_Europe CR = self.run_shift_detection_correction(ref, tgt, **dict(self.coreg_kwargs)) self.assertTrue(CR.success)
def test_shift_calculation_nonquadratic_pixels(self): """Test with default parameters - should compute X/Y shifts properly and write the de-shifted target image.""" # overwrite gt and prj ref = GeoArray(self.ref_path) ref.to_mem() ref.filePath = None ref.gt = [330000.00000001, 5.8932, 0.0, 5862000.0000001, 0.0, -10.1] tgt = GeoArray(self.tgt_path) tgt.to_mem() tgt.filePath = None tgt.gt = [335440.0000001, 5.8933, 0.0, 5866490.0000001, 0.0, -10.1] CR = self.run_shift_detection_correction(ref, tgt, **dict(self.coreg_kwargs, wp=(341500.0, 5861440.0), footprint_poly_ref=None, footprint_poly_tgt=None)) self.assertTrue(CR.success)
def test_shift_calculation_with_image_coords_only(self): """Test with default parameters - should compute X/Y shifts properly and write the de-shifted target image.""" # FIXME fails when executed alone # overwrite gt and prj ref = GeoArray(self.ref_path) ref.to_mem() ref.filePath = None ref.gt = [0, 1, 0, 0, 0, -1] ref.prj = '' tgt = GeoArray(self.tgt_path) tgt.to_mem() tgt.filePath = None tgt.gt = [0, 1, 0, 0, 0, -1] tgt.prj = '' CR = self.run_shift_detection_correction(ref, tgt, **dict(self.coreg_kwargs, wp=(1200, -1600), footprint_poly_ref=None, footprint_poly_tgt=None)) self.assertTrue(CR.success)
def test_shift_calculation_different_geographic_datum(self): """Test if shift computation properly raises a RunTimeError if the matching window is centered at a cloudy image position. """ ref = GeoArray(self.ref_path).to_mem() tgt = GeoArray(self.tgt_path).to_mem() # force to overwrite projection ref.filePath = None ref.prj = EPSG2WKT(3035) # ETRS89_LAEA_Europe with self.assertRaises(RuntimeError): self.run_shift_detection_correction(ref, tgt, **dict(self.coreg_kwargs))