ucar.nc2.iosp.dorade.Doradeiosp 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.
The newest version!
/*
* 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.dorade;
import ucar.ma2.*;
import ucar.nc2.*;
import ucar.nc2.iosp.AbstractIOServiceProvider;
import ucar.atd.dorade.*;
import java.io.*;
import java.util.*;
/**
* IOServiceProvider implementation abstract base class to read/write Dorade files
*/
public class Doradeiosp extends AbstractIOServiceProvider {
protected boolean readonly;
protected Doradeheader headerParser;
public DoradeSweep mySweep = null;
boolean littleEndian;
public boolean isValidFile(ucar.unidata.io.RandomAccessFile raf) {
return (Doradeheader.isValidFile(raf));
}
public String getFileTypeId() {
return "DORADE";
}
public String getFileTypeDescription() {
return "DOppler RAdar Data Exchange Format";
}
/////////////////////////////////////////////////////////////////////////////
// reading
public void open(ucar.unidata.io.RandomAccessFile raf, ucar.nc2.NetcdfFile ncfile,
ucar.nc2.util.CancelTask cancelTask) throws IOException {
super.open(raf, ncfile, cancelTask);
try {
mySweep = new DoradeSweep(raf.getRandomAccessFile());
} catch (DoradeSweep.DoradeSweepException ex) {
ex.printStackTrace();
} catch (java.io.IOException ex) {
ex.printStackTrace();
}
if (mySweep.getScanMode(0) != ScanMode.MODE_SUR) {
// System.err.println("ScanMode is : " + mySweep.getScanMode(0).getName());
//System.exit(1);
}
try {
headerParser = new Doradeheader();
headerParser.read(mySweep, ncfile, null);
} catch (DoradeSweep.DoradeSweepException e) {
e.printStackTrace();
}
ncfile.finish();
}
public Array readData(ucar.nc2.Variable v2, Section section) throws IOException, InvalidRangeException {
Array outputData;
int nSensor = mySweep.getNSensors();
int nRays = mySweep.getNRays();
if (v2.getShortName().equals("elevation")) {
float[] elev = mySweep.getElevations();
outputData = readData1(v2, section, elev);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), elev);
} else if (v2.getShortName().equals("rays_time")) {
Date[] dd = mySweep.getTimes();
double[] d = new double[dd.length];
for (int i = 0; i < dd.length; i++)
d[i] = (double) dd[i].getTime();
outputData = readData2(v2, section, d);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), d);
} else if (v2.getShortName().equals("azimuth")) {
float[] azim = mySweep.getAzimuths();
outputData = readData1(v2, section, azim);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), azim);
} else if (v2.getShortName().startsWith("latitudes_")) {
float[] allLats = new float[nSensor * nRays];
float[] lats;
for (int i = 0; i < nSensor; i++) {
lats = mySweep.getLatitudes(i);
System.arraycopy(lats, 0, allLats, i * nRays, nRays);
}
outputData = readData1(v2, section, allLats);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), allLats);
} else if (v2.getShortName().startsWith("longitudes_")) {
float[] allLons = new float[nSensor * nRays];
float[] lons;
for (int i = 0; i < nSensor; i++) {
lons = mySweep.getLongitudes(i);
System.arraycopy(lons, 0, allLons, i * nRays, nRays);
}
outputData = readData1(v2, section, allLons);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), lons);
} else if (v2.getShortName().startsWith("altitudes_")) {
float[] allAlts = new float[nSensor * nRays];
float[] alts;
for (int i = 0; i < nSensor; i++) {
alts = mySweep.getAltitudes(i);
System.arraycopy(alts, 0, allAlts, i * nRays, nRays);
}
outputData = readData1(v2, section, allAlts);
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), alts);
} else if (v2.getShortName().startsWith("distance_")) {
float[] dist;
int j = 0;
for (int i = 0; i < nSensor; i++) {
String t = "" + i;
if (v2.getShortName().endsWith(t)) {
j = i;
break;
}
}
int nc = mySweep.getNCells(j);
Array data = Array.makeArray(DataType.FLOAT, nc,
(double) mySweep.getRangeToFirstCell(j), (double) mySweep.getCellSpacing(j));
dist = (float[]) data.get1DJavaArray(Float.class);
outputData = readData1(v2, section, dist);
} else if (v2.isScalar()) {
float d = 0.0f;
if (v2.getShortName().equals("Range_to_First_Cell")) {
d = mySweep.getRangeToFirstCell(0);
} else if (v2.getShortName().equals("Cell_Spacing")) {
d = mySweep.getCellSpacing(0);
} else if (v2.getShortName().equals("Fixed_Angle")) {
d = mySweep.getFixedAngle();
} else if (v2.getShortName().equals("Nyquist_Velocity")) {
d = mySweep.getUnambiguousVelocity(0);
} else if (v2.getShortName().equals("Unambiguous_Range")) {
d = mySweep.getunambiguousRange(0);
} else if (v2.getShortName().equals("Radar_Constant")) {
d = mySweep.getradarConstant(0);
} else if (v2.getShortName().equals("rcvr_gain")) {
d = mySweep.getrcvrGain(0);
} else if (v2.getShortName().equals("ant_gain")) {
d = mySweep.getantennaGain(0);
} else if (v2.getShortName().equals("sys_gain")) {
d = mySweep.getsystemGain(0);
} else if (v2.getShortName().equals("bm_width")) {
d = mySweep.gethBeamWidth(0);
}
float[] dd = new float[1];
dd[0] = d;
outputData = Array.factory(v2.getDataType().getPrimitiveClassType(), v2.getShape(), dd);
} else {
Range radialRange = section.getRange(0);
Range gateRange = section.getRange(1);
Array data = Array.factory(v2.getDataType().getPrimitiveClassType(), section.getShape());
IndexIterator ii = data.getIndexIterator();
DoradePARM dp = mySweep.lookupParamIgnoreCase(v2.getShortName());
int ncells = dp.getNCells();
float[] rayValues = new float[ncells];
/*
float[] allValues = new float[nRays * ncells];
for (int r = 0; r < nRays; r++) {
try {
rayValues = mySweep.getRayData(dp, r, rayValues);
} catch (DoradeSweep.DoradeSweepException ex) {
ex.printStackTrace();
}
System.arraycopy(rayValues, 0, allValues, r*ncells, ncells);
} */
for (int r = radialRange.first(); r <= radialRange.last(); r += radialRange.stride()) {
try {
rayValues = mySweep.getRayData(dp, r, rayValues);
} catch (DoradeSweep.DoradeSweepException ex) {
ex.printStackTrace();
}
for (int i = gateRange.first(); i <= gateRange.last(); i += gateRange.stride()) {
ii.setFloatNext(rayValues[i]);
}
}
return data;
//outputData = Array.factory( v2.getDataType().getPrimitiveClassType(), v2.getShape(), allValues);
}
return outputData;
}
public Array readData1(Variable v2, Section section, float[] values) {
Array data = Array.factory(v2.getDataType().getPrimitiveClassType(), section.getShape());
IndexIterator ii = data.getIndexIterator();
Range radialRange = section.getRange(0); // radial range can also be gate range
for (int r = radialRange.first(); r <= radialRange.last(); r += radialRange.stride()) {
ii.setFloatNext(values[r]);
}
return data;
}
public Array readData2(Variable v2, Section section, double[] values) {
Array data = Array.factory(v2.getDataType().getPrimitiveClassType(), section.getShape());
IndexIterator ii = data.getIndexIterator();
Range radialRange = section.getRange(0);
for (int r = radialRange.first(); r <= radialRange.last(); r += radialRange.stride()) {
ii.setDoubleNext(values[r]);
}
return data;
}
// for the compressed data read all out into a array and then parse into requested
protected boolean fill;
protected HashMap dimHash = new HashMap(50);
public static void main(String args[]) throws IOException, InstantiationException, IllegalAccessException {
String fileIn = "/home/yuanho/dorade/swp.1020511015815.SP0L.573.1.2_SUR_v1";
//String fileIn = "c:/data/image/Dorade/n0r_20041013_1852";
ucar.nc2.NetcdfFile.registerIOProvider(ucar.nc2.iosp.dorade.Doradeiosp.class);
ucar.nc2.NetcdfFile ncf = ucar.nc2.NetcdfFile.open(fileIn);
//List alist = ncf.getGlobalAttributes();
ucar.unidata.io.RandomAccessFile file = new ucar.unidata.io.RandomAccessFile(fileIn, "r");
//open1(file, null, null);
//ucar.nc2.Variable v = ncf.findVariable("BaseReflectivity");
//ncf.close();
}
}