ucar.nc2.iosp.sigmet.Ray 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 (c) 1998 - 2010. 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.sigmet;
import ucar.ma2.IndexIterator;
import ucar.ma2.Range;
import ucar.unidata.io.RandomAccessFile;
import java.io.IOException;
import java.util.Comparator;
import java.util.Formatter;
/**
* @author yuanho
* @since Apr 7, 2010
*/
public class Ray {
private short bins;
int dataRead;
int offset;
int offset1;
private float range, step, az, elev;
private short time;
// private float[] val;
String varName;
int nsweep;
short datatype;
public Ray(float range, float step, float az, float elev, short bins, short time, int offset, int dataRead,
int offset1, int nsweep, String name,
short datatype) {
// this.val = new float[bins];
setRange(range);
setStep(step);
setAz(az);
setElev(elev);
setBins(bins);
setTime(time);
setOffset(offset);
setDataRead(dataRead);
setOffset1(offset1);
// setVal(val);
setName(name);
setNsweep(nsweep);
setDataType(datatype);
}
public short getDataType() {
return datatype;
}
public void setDataType(short datatype) {
this.datatype = datatype;
}
public float getRange() {
return range;
}
public void setRange(float range) {
this.range = range;
}
public float getStep() {
return step;
}
public void setStep(float step) {
this.step = step;
}
public int getNsweep() {
return nsweep;
}
public void setNsweep(int nsweep) {
this.nsweep = nsweep;
}
public float getAz() {
if (az < 0.0f & az > -361.0f) {
az = 360.0f + az;
}
return az;
}
public void setAz(float az) {
this.az = az;
}
public float getElev() {
return elev;
}
public void setElev(float elev) {
this.elev = elev;
}
public short getBins() {
return bins;
}
public void setBins(short bins) {
this.bins = bins;
}
public short getTime() {
return time;
}
public void setTime(short time) {
this.time = time;
}
public int getOffset() {
return offset;
}
public void setOffset(int offset) {
this.offset = offset;
}
public int getDataRead() {
return dataRead;
}
public void setDataRead(int dataRead) {
this.dataRead = dataRead;
}
public int getOffset1() {
return offset1;
}
public void setOffset1(int offset1) {
this.offset1 = offset1;
}
/* public float[] getVal() {
return val;
}
public void setVal(float[] val) {
System.arraycopy(val, 0, this.val, 0, bins);
}
*/
public void setName(String name) {
this.varName = name;
}
public String getName() {
return varName;
}
public boolean equals(Object o) {
if (this == o) {
return true;
} else if (o instanceof Ray) {
Ray oo = (Ray) o;
return (range == oo.range & step == oo.step & az == oo.az & elev == oo.elev & bins == oo.bins
& time == oo.time);
} else {
return false;
}
}
public int hashCode() {
return new Float(range).hashCode() + new Float(step).hashCode() + new Float(az).hashCode()
+ new Float(elev).hashCode() + new Short(bins).hashCode() + new Short(time).hashCode();
// val.hashCode();
}
public String toString() {
Formatter sb = new Formatter();
sb.format("Range=%f Step=%f", range, step);
if (az > -361.0f & az < 0.0f) {
az = 360.0f + az;
}
sb.format(" Az=%f Elev=%f Bins=%d Time=%d", az, elev, bins, time);
// for (int i=0; i 0) {
raf.seek(offset);
for (int i = 0; i < dataRead; i++) {
d = raf.readByte();
dd[i] = SigmetIOServiceProvider.calcData(SigmetIOServiceProvider.recHdr, getDataType(), d);
nb++;
}
}
// System.out.println("this is az " + getAz());
raf.seek(offset1);
int cur_len = offset1;
while (nb < (int) bins) {
// --- Check if the code=1 ("1" means an end of a ray)
a00 = raf.readShort();
cur_len = cur_len + 2;
if (a00 == (short) 1) {
for (int uk = 0; uk < (int) bins; uk++) {
dd[uk] = -999.99f;
}
break;
}
if (a00 < 0) { // -- This is data
int nwords = a00 & 0x7fff;
int dataRead1 = nwords * 2;
int pos = 0;
if (cur_len % REC_SIZE == 0) {
pos = 0;
break;
}
raf.seek(cur_len);
for (int i = 0; i < dataRead1; i++) {
d = raf.readByte();
dd[nb] = SigmetIOServiceProvider.calcData(SigmetIOServiceProvider.recHdr, getDataType(), d);
nb = nb + 1;
cur_len = cur_len + 1;
if (nb % REC_SIZE == 0) {
pos = i + 1;
break;
}
}
if (pos > 0) {
break;
}
} else if (a00 > 0 & a00 != 1) {
int num_zero = a00 * 2;
int dataRead1 = num_zero;
for (int k = 0; k < dataRead1; k++) {
dd[nb + k] = SigmetIOServiceProvider.calcData(SigmetIOServiceProvider.recHdr, getDataType(), (byte) 0);
}
nb = nb + dataRead1;
if (cur_len % REC_SIZE == 0) {
break;
}
}
} // ------ end of while for num_bins---------------------------------
for (int i = gateRange.first(); i <= gateRange.last(); i += gateRange.stride()) {
if (i >= bins)
ii.setFloatNext(Float.NaN);
else
ii.setFloatNext(dd[i]);
}
} // end of readData
} // class Ray end------------------------------------------