2019-05-01 19:11:26 +00:00
|
|
|
import glob
|
|
|
|
import os
|
|
|
|
from opendm import log
|
|
|
|
|
|
|
|
class GCPFile:
|
|
|
|
def __init__(self, gcp_path):
|
|
|
|
self.gcp_path = gcp_path
|
|
|
|
self.entries = []
|
|
|
|
self.srs = ""
|
|
|
|
self.read()
|
|
|
|
|
|
|
|
def read(self):
|
|
|
|
if self.exists():
|
|
|
|
with open(self.gcp_path, 'r') as f:
|
|
|
|
contents = f.read().strip()
|
|
|
|
|
|
|
|
lines = map(str.strip, contents.split('\n'))
|
|
|
|
if lines:
|
|
|
|
self.srs = lines[0] # SRS
|
|
|
|
|
|
|
|
for line in lines[1:]:
|
|
|
|
if line != "" and line[0] != "#":
|
|
|
|
parts = line.split()
|
|
|
|
if len(parts) >= 6:
|
|
|
|
self.entries.append(line)
|
|
|
|
else:
|
|
|
|
log.ODM_WARNING("Malformed GCP line: %s" % line)
|
|
|
|
|
|
|
|
def entries_dict(self):
|
|
|
|
for entry in self.entries:
|
|
|
|
parts = entry.split()
|
|
|
|
x, y, z, px, py, filename = parts[:6]
|
|
|
|
extras = " ".join(parts[6:])
|
|
|
|
yield {
|
|
|
|
'x': x,
|
|
|
|
'y': y,
|
|
|
|
'z': z,
|
|
|
|
'px': px,
|
|
|
|
'py': py,
|
|
|
|
'filename': filename,
|
|
|
|
'extras': extras
|
|
|
|
}
|
|
|
|
|
|
|
|
def entry_dict_to_s(self, entry):
|
|
|
|
return "{x} {y} {z} {px} {py} {filename} {extras}".format(**entry).rstrip()
|
|
|
|
|
|
|
|
def exists(self):
|
2019-05-01 20:38:55 +00:00
|
|
|
return self.gcp_path and os.path.exists(self.gcp_path)
|
2019-05-01 19:11:26 +00:00
|
|
|
|
|
|
|
def make_filtered_copy(self, gcp_file_output, images_dir, min_images=3):
|
|
|
|
"""
|
|
|
|
Creates a new GCP file from an existing GCP file includes
|
|
|
|
only the points that reference images existing in the images_dir directory.
|
|
|
|
If less than min_images images are referenced, no GCP copy is created.
|
|
|
|
:return gcp_file_output if successful, None if no output file was created.
|
|
|
|
"""
|
|
|
|
if not self.exists() or not os.path.exists(images_dir):
|
|
|
|
return None
|
|
|
|
|
|
|
|
if os.path.exists(gcp_file_output):
|
|
|
|
os.remove(gcp_file_output)
|
|
|
|
|
|
|
|
files = map(os.path.basename, glob.glob(os.path.join(images_dir, "*")))
|
|
|
|
|
|
|
|
output = [self.srs]
|
|
|
|
files_found = 0
|
|
|
|
|
|
|
|
for entry in self.entries_dict():
|
|
|
|
if entry['filename'] in files:
|
|
|
|
output.append(self.entry_dict_to_s(entry))
|
|
|
|
files_found += 1
|
|
|
|
|
|
|
|
if files_found >= min_images:
|
|
|
|
with open(gcp_file_output, 'w') as f:
|
|
|
|
f.write('\n'.join(output) + '\n')
|
|
|
|
|
|
|
|
return gcp_file_output
|