ucar.nc2.iosp.adde.StationDB Maven / Gradle / Ivy
Go to download
Show more of this group Show more artifacts with this name
Show all versions of netcdf Show documentation
Show all versions of netcdf Show documentation
The NetCDF-Java Library is a Java interface to NetCDF files,
as well as to many other types of scientific data formats.
/*
* 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);
}
}