Specified lat and lon to increase readability (#319)

        points.extend(place_points_in_overlap(nodes, o.geom, dem=gd,

                                              iterative_phase_kwargs=iterative_phase_kwargs))

    session.add_all(points)

    session.commit()

        The geometry of the overlap region

    dem : GeoDataset

         The DEM used to compute point elevations. An elevation of 0 is is used

         The DEM used to compute point elevations. An elevation of 0 is used

         if no DEM is passed in.

    iterative_phase_kwargs : dict

    nodes.remove(source)

    source_camera = source.camera

    for v in valid:

        geom = shapely.geometry.Point(v[0], v[1])

        lon = v[0]

        lat = v[1]

        geom = shapely.geometry.Point(lon, lat)

        point = Points(geom=geom,

                       pointtype=2) # Would be 3 or 4 for ground

        if dem is None:

            height = 0

        else:

            px, py = dem.latlon_to_pixel(v[1], v[0])

            px, py = dem.latlon_to_pixel(lat, lon)

            height = dem.read_array(1, [px, py, 1, 1])[0][0]

        # Get the BCEF coordinate from the lon, lat

        x, y, z = pyproj.transform(lla, ecef, v[0], v[1], height)  # -3000 working well in elysium, need aeroid

        x, y, z = pyproj.transform(lla, ecef, lon, lat, height)

        gnd = csmapi.EcefCoord(x, y, z)

        sic = source_camera.groundToImage(gnd)