ucar.nc2.ft.point.StationHelper Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of org.gecko.ucar.netcdf Show documentation
Re-packaged Unidata netCDF
The newest version!
/*
 * Copyright (c) 1998-2018 John Caron and University Corporation for Atmospheric Research/Unidata
 * See LICENSE for license information.
 */
package ucar.nc2.ft.point;

import ucar.unidata.geoloc.LatLonPoint;
import ucar.unidata.geoloc.LatLonRect;
import ucar.unidata.geoloc.Station;
import java.util.ArrayList;
import java.util.HashMap;
import java.util.List;
import java.util.Map;

/**
 * Helper class for Station Collections.
 * This assumes that calling getData( Station s) is cheap, ie that theres no cheaper filtering to do.
 * 
 * @author caron
 * @since Feb 5, 2008
 */
public class StationHelper {
  private List stations;
  private Map stationHash;
  private static final boolean debug = false;

  public StationHelper() {
    stations = new ArrayList<>();
    stationHash = new HashMap<>();
  }

  public void addStation(StationFeature s) {
    stations.add(s);
    stationHash.put(s.getName(), s);
  }

  public void setStations(List nstations) {
    stations = new ArrayList<>();
    stationHash = new HashMap<>();
    for (StationFeature s : nstations)
      addStation(s);
  }

  private LatLonRect rect;

  public LatLonRect getBoundingBox() {
    if (rect == null) {
      if (stations.isEmpty())
        return null;

      Station s = stations.get(0);
      LatLonPoint llpt = LatLonPoint.create(s.getLatitude(), s.getLongitude());
      rect = new LatLonRect(llpt, 0, 0);
      if (debug)
        System.out.println("start=" + s.getLatitude() + " " + s.getLongitude() + " rect= " + rect.toString2());

      for (int i = 1; i < stations.size(); i++) {
        s = stations.get(i);
        rect.extend(LatLonPoint.create(s.getLatitude(), s.getLongitude()));
        if (debug)
          System.out.println("add=" + s.getLatitude() + " " + s.getLongitude() + " rect= " + rect.toString2());
      }
    }
    if (rect.crossDateline() && rect.getWidth() > 350.0) { // call it global - less confusing
      double lat_min = rect.getLowerLeftPoint().getLatitude();
      double deltaLat = rect.getUpperLeftPoint().getLatitude() - lat_min;
      rect = new LatLonRect(LatLonPoint.create(lat_min, -180.0), deltaLat, 360.0);
    }

    // To give a little "wiggle room", we're going to slightly expand the bounding box.
    double newLowerLeftLat = rect.getLowerLeftPoint().getLatitude() - .0005;
    double newLowerLeftLon = rect.getLowerLeftPoint().getLongitude() - .0005;
    LatLonPoint newLowerLeftPoint = LatLonPoint.create(newLowerLeftLat, newLowerLeftLon);

    double newUpperRightLat = rect.getUpperRightPoint().getLatitude() + .0005;
    double newUpperRightLon = rect.getUpperRightPoint().getLongitude() + .0005;
    LatLonPoint newUpperRightPoint = LatLonPoint.create(newUpperRightLat, newUpperRightLon);

    rect.extend(newLowerLeftPoint);
    rect.extend(newUpperRightPoint);

    return rect;
  }

  public List getStations(LatLonRect boundingBox) {
    if (boundingBox == null)
      return getStations();

    List result = new ArrayList<>();
    for (StationFeature s : stations) {
      LatLonPoint latlonPt = LatLonPoint.create(s.getLatitude(), s.getLongitude());
      if (boundingBox.contains(latlonPt))
        result.add(s);
    }
    return result;
  }

  public List getStationFeatures(LatLonRect boundingBox) {
    if (boundingBox == null)
      return stations;

    List result = new ArrayList<>();
    for (StationFeature s : stations) {
      LatLonPoint latlonPt = LatLonPoint.create(s.getLatitude(), s.getLongitude());
      if (boundingBox.contains(latlonPt))
        result.add(s);
    }
    return result;
  }

  public StationFeature getStation(String name) {
    return stationHash.get(name);
  }

  public List getStationFeatures() {
    return stations;
  }

  public List getStations() {
    List result = new ArrayList<>(stations.size());
    result.addAll(stations);
    return result;
  }

  public List getStationFeaturesFromNames(List stnNames) {
    List result = new ArrayList<>(stnNames.size());
    for (String ss : stnNames) {
      StationFeature s = stationHash.get(ss);
      if (s != null)
        result.add(s);
    }
    return result;
  }

  public List getStationFeatures(List stations) {
    List result = new ArrayList<>(stations.size());
    for (Station s : stations) {
      StationFeature ss = stationHash.get(s.getName());
      if (ss != null)
        result.add(ss);
    }
    return result;
  }

  public StationFeature getStationFeature(Station stn) {
    return stationHash.get(stn.getName());
  }

  public List getStations(List stnNames) {
    List result = new ArrayList<>(stnNames.size());
    for (String ss : stnNames) {
      Station s = stationHash.get(ss);
      if (s != null)
        result.add(s);
    }
    return result;
  }

  public StationHelper subset(LatLonRect bb) {
    StationHelper result = new StationHelper();
    result.setStations(getStationFeatures(bb));
    return result;
  }

  public StationHelper subset(List stns) {
    StationHelper result = new StationHelper();
    result.setStations(stns);
    return result;
  }

}