Merge pull request #621 from gsn9/loggingForSpatial

import time

import warnings

import logging

import json

from subprocess import CalledProcessError

  FROM iid GROUP BY iid.geom) AS row WHERE array_length(intersections, 1) > 1;

"""

# set up the logger file

log = logging.getLogger(__name__)

def place_points_in_overlaps(size_threshold=0.0007,

                             distribute_points_kwargs={},

                             cam_type='csm',

    geom = overlap.geom

    valid = compgeom.distribute_points_in_geom(geom, **distribute_points_kwargs, **kwargs)

    if not valid.any():

        warnings.warn('Failed to distribute points in overlap')

        log.warning('Failed to distribute points in overlap')

        return []

    print(f'Have {len(valid)} potential points to place.')

            # Extract ORB features in a sub-image around the desired point

            image_roi = roi.Roi(node.geodata, sample, line, size_x=size, size_y=size)

            if image_roi.variance == 0:

                warnings.warn(f'Failed to find interesting features in image {node.image_name}.')

                log.warning(f'Failed to find interesting features in image {node.image_name}.')

                continue

            image = image_roi.clip()

                break

        if interesting is None:

            warnings.warn('Unable to find an interesting point, falling back to the a priori pointing')

            log.warning('Unable to find an interesting point, falling back to the a priori pointing')

            newsample = sample

            newline = line

        else: