]
log("NED reference location:", ref)
# local surface approximation
srtm.initialize(ref, 6000, 6000, 30)
surface.update_srtm_elevations(proj)
proj.save_images_info()
# camera calibration
K = camera.get_K()
print("K:", K)
# fire up the matcher
matcher.configure()
matcher.find_matches(proj.image_list,
K,
transform=args.filter,
sort=True,
review=False)
feature_count = 0
image_count = 0
for image in proj.image_list:
feature_count += image.num_features
image_count += 1
log("Average # of features per image found = %.0f" %
(feature_count / image_count))
# The following code is deprecated ...
do_old_match_consolodation = False
if do_old_match_consolodation:
# build a list of all 'unique' keypoints. Include an index to each
matcher_node.setFloat('ground_m', args.ground)
# save any config changes
proj.save()
# camera calibration
K = camera.get_K()
# print("K:", K)
log("Matching features")
# fire up the matcher
matcher.configure()
matcher.find_matches(proj,
K,
strategy=args.match_strategy,
transform=args.filter,
sort=True,
review=False)
feature_count = 0
image_count = 0
for image in proj.image_list:
feature_count += image.num_features
image_count += 1
log("Average # of features per image found = %.0f" %
(feature_count / image_count))
state.update("STEP3a")
matches_name = os.path.join(proj.analysis_dir, "matches_grouped")