org.wamblee.gpx.GpxPlotter Maven / Gradle / Ivy

Go to download
/*
 * Copyright 2006 the original author or authors.
 * 
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 * 
 *      http://www.apache.org/licenses/LICENSE-2.0
 * 
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */ 

package org.wamblee.gpx;

import java.awt.Color;
import java.awt.Image;
import java.io.File;
import java.io.FileInputStream;
import java.io.FileOutputStream;
import java.io.IOException;
import java.util.ArrayList;
import java.util.List;

import org.jfree.chart.ChartFactory;
import org.jfree.chart.ChartFrame;
import org.jfree.chart.ChartUtilities;
import org.jfree.chart.JFreeChart;
import org.jfree.chart.axis.NumberAxis;
import org.jfree.chart.plot.PlotOrientation;
import org.jfree.chart.plot.XYPlot;
import org.jfree.chart.renderer.xy.XYItemRenderer;
import org.jfree.chart.renderer.xy.XYLineAndShapeRenderer;
import org.jfree.data.xy.XYSeries;
import org.jfree.data.xy.XYSeriesCollection;
import org.wamblee.general.Pair;
import org.wamblee.gps.geometry.Plane;
import org.wamblee.gps.geometry.Point;
import org.wamblee.gps.geometry.ReferenceCoordinateSystem;
import org.wamblee.gps.track.TrackSegment;
import org.wamblee.utils.JpegUtils;

/**
 * Parses a GPX file and prints out a data file with each trackpoints distance from the start of the 
 * track and its elevation, separated 0by a space. 
 *
 * @author Erik Brakkee
 */
public class GpxPlotter {
    
    public static void main(String[] aArgs) throws Exception { 
        File file = new File(aArgs[0]);
        GpxParser parser = new GpxParser(); 
        TrackSegment track = parser.parse(file.getName(), new FileInputStream(file));
        
        TrackStatistics profile = new TrackStatistics(track);
        profile.writeHeightProfilePng(new FileOutputStream("x.png"), 600, 300);
    
        List> elevationProfile = computeElevationProfile(track);
        printTrack(elevationProfile); 
        computeTotalClimb(elevationProfile);
        plotElevationProfile(elevationProfile);
        List> trackXy = computeTrackXY(track);
        List> trackLatLon = computeTrackLatLon(track);
        plotTrack(trackLatLon);
    }
    
    private static List> computeElevationProfile(TrackSegment aTrack) {
        List> results = new ArrayList>();
        double distance = 0.0; 
        for (int i = 0; i < aTrack.size(); i++) { 
            Point point = aTrack.getPoint(i);
            results.add(new Pair(distance, point.getCoordinates().getX3()));
            if ( i+1 < aTrack.size()) { 
                Point nextPoint = aTrack.getPoint(i+1); 
                distance += ReferenceCoordinateSystem.distance(point, nextPoint);
            }
        }
        return results; 
    }
    
    private static List> computeTrackXY(TrackSegment aTrack) {
        Point reference = aTrack.getPoint(0);
        Plane plane = new Plane(reference, reference); // assume the earth is spherical.
        List> results = new ArrayList>();
        for (int i = 0; i < aTrack.size(); i++) { 
            Point point = aTrack.getPoint(i);
            Pair projection = plane.normalizedProjection(point);
            results.add(projection);
            System.out.println(point);
        }
        return results; 
    }
    
    private static List> computeTrackLatLon(TrackSegment aTrack) {
        List> results = new ArrayList>();
        for (int i = 0; i < aTrack.size(); i++) { 
            Point point = aTrack.getPoint(i);
            results.add(new Pair(point.getCoordinates().getX1(), point.getCoordinates().getX2()));
        }
        return results; 
    }
    
    
    
    private static void printTrack(List> aHeightProfile) { 
       for (Pair point: aHeightProfile) { 
           System.out.println(point.getFirst() + " " + point.getSecond());
       }
    }
    
    private static void computeTotalClimb(List> aHeightProfile) {
        double result = 0.0;
        
        double lastHeight = aHeightProfile.get(0).getSecond();
        for ( int i = 1; i < aHeightProfile.size(); i++) { 
            double height = aHeightProfile.get(i).getSecond();
            if ( height > lastHeight) { 
                result += (height-lastHeight); 
            }
            lastHeight = height; 
        }
        System.out.println("Total climb: " + result);
    }
    
    private static void plotElevationProfile(List> aHeightProfile) throws IOException {
        XYSeriesCollection dataset = createDataset(aHeightProfile, "height");
        JFreeChart chart = ChartFactory.createXYLineChart(
                "Height Profile", 
                "Distance(m)",
                "Height(m)",
                dataset,
                PlotOrientation.VERTICAL,
                true,
                true,
                false);
        ChartUtilities.writeChartAsPNG(new FileOutputStream("height.png"), chart, 600, 300);
        ChartFrame frame = new ChartFrame("test", chart);
        frame.pack();
        frame.setVisible(true);
    }
    
    private static void plotTrack(List> aPoints) throws IOException, InterruptedException {
        XYSeriesCollection dataset = createDataset(aPoints, "track");
        JFreeChart chart = createLineChart(dataset);
        
        Pair,Pair> bounds = getBounds(aPoints); 
        
        chart.getXYPlot().getDomainAxis().setLowerBound(bounds.getFirst().getFirst());
        chart.getXYPlot().getDomainAxis().setUpperBound(bounds.getFirst().getSecond());
        
        chart.getXYPlot().getRangeAxis().setLowerBound(bounds.getSecond().getFirst());
        chart.getXYPlot().getRangeAxis().setUpperBound(bounds.getSecond().getSecond());
        
        Image background = JpegUtils.loadJpegImage(new FileInputStream("/home/erik/vakantie.jpg"));
        chart.getPlot().setBackgroundImage(background);
       
        XYLineAndShapeRenderer renderer = (XYLineAndShapeRenderer)chart.getXYPlot().getRenderer();
        renderer.setShapesVisible(true);
        renderer.setShapesFilled(true);
        renderer.setPaint(Color.BLACK);
     
        ChartUtilities.writeChartAsPNG(new FileOutputStream("test.png"), chart, 1280, 800);
        ChartFrame frame = new ChartFrame("test", chart);
        frame.pack();
        frame.setVisible(true);
    }

    /**
     * @param dataset
     * @return
     */
    private static JFreeChart createLineChart(XYSeriesCollection dataset) {
        NumberAxis xAxis = new NumberAxis("S->N");
        xAxis.setAutoRangeIncludesZero(false);
        
        NumberAxis yAxis = new NumberAxis("W->E");
        yAxis.setAutoRangeIncludesZero(false);
       
        XYItemRenderer renderer = new XYLineAndShapeRenderer(true, false);
        XYPlot plot = new ZoomableBackgroundXYPlot(dataset, xAxis, yAxis, renderer);
        plot.setOrientation(PlotOrientation.HORIZONTAL);
      
        JFreeChart chart = new JFreeChart(
            "Track", JFreeChart.DEFAULT_TITLE_FONT, plot, true
        );

        return chart;
        /*
        JFreeChart chart = ChartFactory.createXYLineChart(
                "Track", 
                "S->N",
                "W->E",
                dataset,
                PlotOrientation.HORIZONTAL,
                true,
                true,
                false);
        return chart;
        */
    }

    /**
     * @param aHeightProfile
     * @return
     */
    private static XYSeriesCollection createDataset(List> aHeightProfile, String aName) {
        XYSeries series = new XYSeries(aName, false);
        for (Pair point: aHeightProfile) { 
            series.add(point.getFirst(), point.getSecond());
        }
        XYSeriesCollection dataset = new XYSeriesCollection(series);
        return dataset;
    }

    private static Pair,Pair> getBounds(List> aList) { 
        Pair first = aList.get(0); 
        double minx= first.getFirst();
        double maxx = minx;
        double miny = first.getSecond();
        double maxy = miny;
        
        for (int i = 0; i < aList.size(); i++) { 
            Pair value = aList.get(i);
            minx = Math.min(minx, value.getFirst());
            maxx = Math.max(maxx, value.getFirst());
            miny = Math.min(miny, value.getSecond());
            maxy = Math.max(maxy, value.getSecond());
        }
        if ( maxx == minx ) { 
            maxx += 1.0; // to avoid problems. 
        }
        if ( maxy == miny ) { 
            maxy += 1.0; // to avoid problems.
        }
        final double paddingFactor = 0.3; // allow some space around min and max
        return new Pair,Pair>(
                new Pair( minx - paddingFactor*(maxx-minx), 
                                         maxx + paddingFactor*(maxx-minx)),
                new Pair( miny - paddingFactor*(maxy-miny), 
                                         maxy + paddingFactor*(maxy-miny))
                );
    }
}