centroid_ard_grid.py

import json

import click
import cligj

from shapely.geometry import shape, mapping


@click.command()
@click.pass_context
@cligj.features_in_arg
@cligj.sequence_opt
def centroid(ctx, features, sequence):
    """Get Centroid."""
    for f in features:
        geom = shape(f['geometry'])
        f['geometry'] = mapping(geom.centroid)
        print(json.dumps(f))


if __name__ == '__main__':
    centroid()

ogr2ogr -f GeoJSON -t_srs EPSG:4326 conus_ard_grid_wgs84.geojson conus_ard_grid.shp
fio cat conus_ard_grid_wgs84.geojson | python centroid_ard.py | fio collect > conus_ard_grid_centroid.geojson