source: tmcsimulator/trunk/webapps/GTEC/data_layers/vds_intersect.py @ 556

import json, sys
# Read 2 filenames from command line
# Output to stdout features in file2 not in file1

if len(sys.argv) == 1:
    print "Must provide 2 command line arguments"
    print 'Number of arguments:', len(sys.argv)-1
    exit
    
with open(sys.argv[1], "r") as read_file:
    j = json.load(read_file)

features1 = j["features"]

# Create a set of ID's from first file of features
alpha = set()
for f in features1:
    alpha.add(f["id"]);

with open(sys.argv[2], "r") as read_file:
    j = json.load(read_file)

features2 = j["features"]

# Lookup each feature ID from file2 in set     
bravo = []
for f in features2:
    target = f["id"]
    if target not in alpha:
        bravo.append(target);

# bravo contains items in dat2 not in dat1
for b in bravo:
    print b
    
# Note features dicts are saved for future use