Skip to content

Instantly share code, notes, and snippets.

@adv27

adv27/merge-geojsons.py

Forked from migurski/merge-geojsons.py
Created December 3, 2019 15:17
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save adv27/34e77dc9648bb3560baaac58a18f4da0 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save adv27/34e77dc9648bb3560baaac58a18f4da0 to your computer and use it in GitHub Desktop.
Download ZIP

Revisions

  1. @migurski migurski created this gist Sep 21, 2012.
    56 changes: 56 additions & 0 deletions merge-geojsons.py
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,56 @@
    from json import load, JSONEncoder
    from optparse import OptionParser
    from re import compile

    float_pat = compile(r'^-?\d+\.\d+(e-?\d+)?$')
    charfloat_pat = compile(r'^[\[,\,]-?\d+\.\d+(e-?\d+)?$')

    parser = OptionParser(usage="""%prog [options]
    Group multiple GeoJSON files into one output file.
    Example:
    python %prog -p 2 input-1.json input-2.json output.json""")

    defaults = dict(precision=6)

    parser.set_defaults(**defaults)

    parser.add_option('-p', '--precision', dest='precision',
    type='int', help='Digits of precision, default %(precision)d.' % defaults)

    if __name__ == '__main__':

    options, args = parser.parse_args()
    infiles, outfile = args[:-1], args[-1]

    outjson = dict(type='FeatureCollection', features=[])

    for infile in infiles:
    injson = load(open(infile))

    if injson.get('type', None) != 'FeatureCollection':
    raise Exception('Sorry, "%s" does not look like GeoJSON' % infile)

    if type(injson.get('features', None)) != list:
    raise Exception('Sorry, "%s" does not look like GeoJSON' % infile)

    outjson['features'] += injson['features']

    encoder = JSONEncoder(separators=(',', ':'))
    encoded = encoder.iterencode(outjson)

    format = '%.' + str(options.precision) + 'f'
    output = open(outfile, 'w')

    for token in encoded:
    if charfloat_pat.match(token):
    # in python 2.7, we see a character followed by a float literal
    output.write(token[0] + format % float(token[1:]))

    elif float_pat.match(token):
    # in python 2.6, we see a simple float literal
    output.write(format % float(token))

    else:
    output.write(token)