geo.d

import std.stdio;
import std.file;
import std.algorithm.iteration;
import std.range;
import std.conv;

void main()
{
        auto bytes = read("geo_points.txt");
        auto s = cast(string)bytes; // avoid Unicode validation of readText

        double lat_tot = 0;
        double lon_tot = 0;
        ulong line_cnt = 0;

        foreach (line; s.splitter('\n'))
        {
                ++line_cnt;
                auto fields = line.splitter(',').array;
                lat_tot += fields[0].to!double;
                lon_tot += fields[1].to!double;
        }

        double centroid_lat = lat_tot / line_cnt;
        double centroid_lon = lon_tot / line_cnt;

        writefln("Total number of points: %d", line_cnt);
        writefln("Centroid { lat: %f, lon: %f }", centroid_lat, centroid_lon);
}