Skip to content
GitLab
Closed
Created by initze@initze

Computing footprint error

Hi Daniel,

I am running into the following error, which I am not able to fix (see error message further down)

Data overlap. Same projection, type, source, nodata, both are overlapping

  • ref image doesn't have super great quality, this was more a technical test. Same error in other places

Data download

image

I tried to play around with no data (with, without, image size, etc.)

Platform:

  • Windows 10
  • Anaconda environment with python 3.7.X
  • arosics (1.4.3) - same error with higher version

Input

REFERENCE = r'reference_image\reference.tif'
INFILE = r'input_images\target.tif'    
outfile = r'output_files\out.tif'  
# detect and correct global spatial shift 
CR = COREG(REFERENCE, INFILE,
           max_shift=200, path_out=outfile, 
           fmt_out='GTIFF',
           out_crea_options=['COMPRESS=DEFLATE'],
           r_b4match=1, 
           s_b4match=1)

Output

Calculating footprint polygon and actual data corner coordinates for reference image...

Polygonize progress |==================================================| 100.0% Complete => 0:00:01

AssertionError Traceback (most recent call last) in 1 # detect and correct global spatial shift ----> 2 CR = COREG(REFERENCE, infile, 3 max_shift=200, path_out=outfile, 4 fmt_out='GTIFF', 5 out_crea_options=['COMPRESS=DEFLATE'],

~\AppData\Local\Continuum\anaconda3\envs\arosics2\lib\site-packages\arosics\CoReg.py in init(self, im_ref, im_tgt, path_out, fmt_out, out_crea_options, r_b4match, s_b4match, wp, ws, max_iter, max_shift, align_grids, match_gsd, out_gsd, target_xyGrid, resamp_alg_deshift, resamp_alg_calc, footprint_poly_ref, footprint_poly_tgt, data_corners_ref, data_corners_tgt, nodata, calc_corners, binary_ws, mask_baddata_ref, mask_baddata_tgt, CPUs, force_quadratic_win, progress, v, path_verbose_out, q, ignore_errors) 398 # try: 399 gdal.AllRegister() --> 400 self._get_image_params() 401 self._set_outpathes(im_ref, im_tgt) 402 self.grid2use = 'ref' if self.shift.xgsd <= self.ref.xgsd else 'shift'

~\AppData\Local\Continuum\anaconda3\envs\arosics2\lib\site-packages\arosics\CoReg.py in _get_image_params(self) 532 533 def _get_image_params(self) -> None: --> 534 self.ref = GeoArray_CoReg(self.params, 'ref') 535 self.shift = GeoArray_CoReg(self.params, 'shift') 536

~\AppData\Local\Continuum\anaconda3\envs\arosics2\lib\site-packages\arosics\CoReg.py in init(self, CoReg_params, imID) 116 117 with warnings.catch_warnings(record=True) as w: --> 118 _ = self.footprint_poly # execute getter 119 120 if len(w) > 0 and 'disjunct polygone(s) outside' in str(w[-1].message):

~\AppData\Local\Continuum\anaconda3\envs\arosics2\lib\site-packages\geoarray\baseclasses.py in footprint_poly(self) 499 500 # validation --> 501 assert not polyVertices_outside_poly(self._footprint_poly, self.box.mapPoly, tolerance=1e-5),
502 "Computing footprint polygon for %s '%s' failed. The resulting polygon is partly or completely "
503 "outside of the image bounds." % (self.class.name, self.basename)

AssertionError: Computing footprint polygon for GeoArray_CoReg 'reference' failed. The resulting polygon is partly or completely outside of the image bounds.