Convert csv RGB layer data to geojson points

convert_csv_to_geo_json.py

import csv, json, pdb;
from geojson import Feature, FeatureCollection, Point

def convertBox2MidPoint(lower_lat, lower_lon, upper_lat, upper_lon):
    delta_lon = abs(lower_lon - upper_lon) / 2
    delta_lat = abs(lower_lat - upper_lat) / 2
    mid_lon = lower_lon + delta_lon
    mid_lat = lower_lat + delta_lat
    # geojson is lon, lat ordering
    return (mid_lon, mid_lat)

features = []
with open('csv/overlay_0.csv') as csvfile:
    # csv.QUOTE_NONNUMERIC -> convert non quoted inputs to floats
    reader = csv.reader(csvfile, quoting=csv.QUOTE_NONNUMERIC)
    for row in reader:
        # only iterested in the red channel...as it's the only data
        # lower_* is left, upper_* is right in the source data
        # skip null data
        intensity = row[4]
        if intensity == 0:
            continue
        # csv format is lower_lat, lower_long, upper_lat, upper_long
        result_tuple = convertBox2MidPoint(*row[0:4])
        features.append(
            Feature(
                geometry = Point(result_tuple),
                properties = {
                    'intensity': intensity
                }
            )
        )

collection = FeatureCollection(features)
with open("geo_json/overlay_0.json", "w") as file:
    file.write('%s' % collection)