Exemple #1
0
 def __str__(self):
     name = self.name if self.name != '' else 'Empty'
     desc = self.desc if self.desc != '' else 'Empty'
     ufactor = user_units_per_meter(self.user_units)
     n, e, z = [d * ufactor for d in (self.y, self.x, self.z)]
     s = "{0},{1:.4f},{2:.4f},{3:.4f},{4}".format(name, n, e, z, desc)
     return s
Exemple #2
0
    for i in range(len(local_ref)):
        grid = grid_ref[i]
        print('\n{0:>6s}  ({1}): {2:12.4f}, {3:12.4f} (samples={4}) {5}'.format(
            grid.name, grid.type, grid.x, grid.y, grid.samples, grid.desc))

        local = local_ref[i]
        x, y = xfm.to_grid(local.x, local.y)
        dist = math.sqrt((x - grid.x)**2 + (y - grid.y)**2)
        print('{0:>6s} ({1}): {2:12.4f}, {3:12.4f} (error={4:.2f}) {5}'
              .format(local.name, local.type, local.x, local.y, dist, local.desc))
        print('{0:>6s} ({1}): {2:12.4f}, {3:12.4f} (error={4:.2f}) {5}'
              .format(local.name, 'XFORM', x, y, dist, local.desc))

    print()
    for p in grid_pts:
        ufactor = utils.user_units_per_meter('sft')
        x, y = xfm.to_local(p.x, p.y)
        x *= ufactor
        y *= ufactor
        z = p.z * ufactor
        print('{0:d},{1:.4f},{2:.4f},{3:.4f},RIVER'.format(int(p.name) + 1000, y, x, z))

    filename_local_points = 'data/150419_Stapp_PTS.txt'
    filename_grid_points = 'data/Waypoints.gpx'
    projection = ProjectionLC('SPCS', '0401')

    local_pts = Point.read_local_points(filename_local_points, 'sft')
    if not local_pts:
        print("No points. Nothing to do.")
        exit(1)