generate a gridded set of points in geojson format

geojson_grid.py

import json

def feat(x,y):
  feat = { "type": "Feature", 
          "geometry"  : { "type": "Point", "coordinates": [ x, y ] },
          "properties": { "label": "%s,%s" % (x,y) }
      }
  return feat

def grid(minx,miny,maxx,maxy,increment):
    geojson = {
      "type": "FeatureCollection",
    }

    features = []
    
    for x in range(minx,maxx,increment):
        for y in range(miny,maxy,increment):
            features.append(feat(x,y))
            features.append(feat(x,y+increment))
            features.append(feat(x+increment,y))
            features.append(feat(x+increment,y+increment))
    
    geojson["features"] = features
    
    return json.dumps(geojson, indent=4)

if __name__ == "__main__":
    
    print grid(-180,-90,180,90,10);