ucar.nc2.iosp.adde.StationDB Maven / Gradle / Ivy

Go to download
/*
 * Copyright 1998-2009 University Corporation for Atmospheric Research/Unidata
 *
 * Portions of this software were developed by the Unidata Program at the
 * University Corporation for Atmospheric Research.
 *
 * Access and use of this software shall impose the following obligations
 * and understandings on the user. The user is granted the right, without
 * any fee or cost, to use, copy, modify, alter, enhance and distribute
 * this software, and any derivative works thereof, and its supporting
 * documentation for any purpose whatsoever, provided that this entire
 * notice appears in all copies of the software, derivative works and
 * supporting documentation.  Further, UCAR requests that the user credit
 * UCAR/Unidata in any publications that result from the use of this
 * software or in any product that includes this software. The names UCAR
 * and/or Unidata, however, may not be used in any advertising or publicity
 * to endorse or promote any products or commercial entity unless specific
 * written permission is obtained from UCAR/Unidata. The user also
 * understands that UCAR/Unidata is not obligated to provide the user with
 * any support, consulting, training or assistance of any kind with regard
 * to the use, operation and performance of this software nor to provide
 * the user with any updates, revisions, new versions or "bug fixes."
 *
 * THIS SOFTWARE IS PROVIDED BY UCAR/UNIDATA "AS IS" AND ANY EXPRESS OR
 * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
 * DISCLAIMED. IN NO EVENT SHALL UCAR/UNIDATA BE LIABLE FOR ANY SPECIAL,
 * INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING
 * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
 * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
 * WITH THE ACCESS, USE OR PERFORMANCE OF THIS SOFTWARE.
 */
package ucar.nc2.iosp.adde;

import ucar.nc2.util.CancelTask;

import java.io.*;
import java.util.*;
import java.net.URL;

import edu.wisc.ssec.mcidas.adde.AddePointDataReader;
import edu.wisc.ssec.mcidas.adde.AddeException;
import edu.wisc.ssec.mcidas.McIDASUtil;

/**
 * Reads station table ascii output from mcidas STNLIST command
 */
public class StationDB {
  static private boolean debugOpen = false, debugCall = false, debugParse = false;

  private ArrayList stations = new ArrayList();

  /**
   * Reads station table ascii output from mcidas STNLIST command
   */
  public StationDB(String urlString) throws IOException {
    long start = System.currentTimeMillis();

    InputStream ios = null;
    if (urlString.startsWith("http:")) {
      URL url = new URL(urlString);
      ios = url.openStream();
      if (debugOpen) System.out.println("opened URL " + urlString);
    } else if (urlString.startsWith("resource:")) {
      ClassLoader cl = getClass().getClassLoader();
      ios = cl.getResourceAsStream(urlString.substring(9));
      if (debugOpen) System.out.println("opened resource " + urlString);
    } else {
      ios = new FileInputStream(urlString);
      if (debugOpen) System.out.println("opened file " + urlString);
    }

    // DataInputStream dataIS = new DataInputStream( new BufferedInputStream(ios, 20000));
    BufferedReader dataIS = new BufferedReader(new InputStreamReader(ios));

    dataIS.readLine();
    dataIS.readLine();
    int count = 0;
    while (true) {
      String line = dataIS.readLine();
      if (line == null) break;
      if (line.length() < 85) break;

      String idn = line.substring(0, 5);
      String id = line.substring(6, 10);
      String name = line.substring(13, 33);
      String type = line.substring(34, 51);
      String state = line.substring(52, 54);
      String country = line.substring(55, 57);
      String latS = line.substring(58, 68);
      String lonS = line.substring(69, 79);
      String elevS = line.substring(80, 85);
      Station s = new Station(idn, id, name, type, state, country, latS, lonS, elevS);
      stations.add(s);

      if (debugParse) System.out.println(s);
      count++;
      // if (count > 1) break;
    }

    dataIS.close();

    if (debugCall) {
      long took = System.currentTimeMillis() - start;
      System.out.println(" read " + urlString + " count=" + stations.size() + " took=" + took + " msec ");
    }
  }

  /**
   * Reads station info by making a call to ADDE server: very slow!!!
   */
  StationDB(String location, CancelTask cancel) throws IOException {
    HashMap hashStations = new HashMap(5000);

    try {
      if (debugOpen) System.out.println("Call ADDE Server " + location);
      long start = System.currentTimeMillis();
      AddePointDataReader reader = AddeStationObsDataset.callAdde(location + "&num=all¶m=ID LAT LON ZS"); // LOOK
      long took = System.currentTimeMillis() - start;
      if (debugOpen) System.out.println(" time = " + took);

      int[][] stationData = reader.getData();
      int nparams = stationData.length;
      int nstations = stationData[0].length;
      System.out.println(" nparams= " + nparams + " nstations=" + nstations);
      System.out.println(" size= " + (nparams * nstations * 4) + " bytes");

      int[] scales = reader.getScales();
      double[] scaleFactor = new double[scales.length];
      for (int i = 0; i < nparams; i++) {
        scaleFactor[i] = (scales[i] == 0) ? 1.0 : 1.0 / Math.pow(10.0, (double) scales[i]);
      }

      if ((cancel != null) && cancel.isCancel()) return;

      int last = 0;
      for (int i = 0; i < nstations; i++) {
        String stnId = McIDASUtil.intBitsToString(stationData[0][i]);
        if (!hashStations.containsKey(stnId)) {
          last = i;
          double lat = stationData[1][i] * scaleFactor[1];
          double lon = stationData[2][i] * scaleFactor[2];
          double elev = stationData[3][i] * scaleFactor[3];
          hashStations.put(stnId, new Station(stnId, "", lat, lon, elev));
        }
        if ((cancel != null) && cancel.isCancel()) return;
      }
      if (debugCall) System.out.println(" hashStations count= " + hashStations.size() + " last = " + last);
      List stationList = new ArrayList(hashStations.keySet());
      Collections.sort(stationList);
      stations = new ArrayList();
      for (int i = 0; i < stationList.size(); i++) {
        String id = (String) stationList.get(i);
        stations.add(hashStations.get(id));
        if ((cancel != null) && cancel.isCancel()) return;
      }

    } catch (AddeException e) {
      e.printStackTrace();
      throw new IOException(e.getMessage());
    }
  }

  public ArrayList getStations() {
    return stations;
  }

  public class Station extends ucar.unidata.geoloc.StationImpl {
    String idn, id, name, type, state, country, desc;
    double lat, lon, elev;

    Station(String idn, String id, String name, String type, String state, String country, String latS,
            String lonS, String elevS) {
      if (debugParse)
        System.out.println("-" + idn + "-" + id + "-" + name + "-" + type + "-" + state + "-" + country + "-" + latS + "-" + lonS + "-" + elevS);

      this.idn = idn.trim();
      this.id = id.trim();
      this.name = name.trim();
      this.type = type;
      this.state = state.trim();
      this.country = country.trim();
      this.lat = parseDegree(latS);
      this.lon = -1.0 * parseDegree(lonS); // LOOK : degrees west ?? !!

      if (this.state.length() > 0)
        desc = this.name + ", " + this.state + ", " + this.country;
      else
        desc = this.name + ", " + this.country;

      try {
        elev = Double.parseDouble(elevS);
      } catch (NumberFormatException e) {
        e.printStackTrace();
      }
    }

    private double parseDegree(String s) {
      String degS = s.substring(0, 4);
      String minS = s.substring(5, 7);
      String secS = s.substring(8, 10);

      try {
        double deg = Double.parseDouble(degS);
        double min = Double.parseDouble(minS);
        double sec = Double.parseDouble(secS);
        if (deg < 0)
          return deg - min / 60 - sec / 3600;
        else
          return deg + min / 60 + sec / 3600;
      } catch (NumberFormatException e) {
        e.printStackTrace();
      }
      return 0.0;
    }

    Station(String id, String desc, double lat, double lon, double elev) {
      this.id = id;
      this.desc = desc;
      this.lat = lat;
      this.lon = lon;
      this.elev = elev;
    }

    public String toString() {
      return idn + " " + id + " " + name + " " + type + " " + state + " " + country + " " + lat + " " + lon + " " + elev;
    }

    public String getName() {
      return id;
    }

    public String getDescription() {
      return desc;
    }

    public String getWmoId() {
      return idn; // ??
    }

    public double getLatitude() {
      return lat;
    }

    public double getLongitude() {
      return lon;
    }

    public double getAltitude() {
      return elev;
    }

    public int compareTo(Station so) {
      return name.compareTo(so.getName());
    }
  }

  static String testName = "C:/data/station/adde/STNDB.TXT";
  //static String testName = "M:/temp/STNDB.TXT";
  static String testName2 = "http://localhost:8080/test/STNDB.TXT";

  static public void main(String[] args) throws IOException {
    long start = System.currentTimeMillis();
    StationDB stnDB = new StationDB(testName);
    long took = System.currentTimeMillis() - start;
    List list = stnDB.getStations();
    double sum = 0.0;
    for (int i = 0; i < list.size(); i++) {
      Station s = (Station) list.get(i);
      sum += s.getLatitude();
    }
    System.out.println(" read " + testName + " count=" + list.size() + " took=" + took + " msec " + " sum= " + sum);
  }
}