Java Code Examples for ucar.nc2.constants.AxisType#GeoZ
The following examples show how to use
ucar.nc2.constants.AxisType#GeoZ .
You can vote up the ones you like or vote down the ones you don't like,
and go to the original project or source file by following the links above each example. You may check out the related API usage on the sidebar.
Example 1
| Source File: M3IOConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 6 votes |
|
@Override
protected AxisType getAxisType(VariableDS.Builder ve) {
String vname = ve.shortName;
if (vname.equalsIgnoreCase("x"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("y"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("time"))
return AxisType.Time;
if (vname.equalsIgnoreCase("level"))
return AxisType.GeoZ;
return null;
}
Example 2
| Source File: CoordSysBuilder.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 6 votes |
|
/**
* Turn the variable into a coordinate axis, if not already. Add to the dataset, replacing variable if needed.
*
* @return variable as a coordinate axis
*/
public CoordinateAxis makeIntoCoordinateAxis() {
if (axis != null)
return axis;
// if not a CoordinateAxis, will turn into one
v = axis = ds.addCoordinateAxis((VariableDS) v);
if (axisType != null) {
axis.setAxisType(axisType);
axis.addAttribute(new Attribute(_Coordinate.AxisType, axisType.toString()));
if (((axisType == AxisType.Height) || (axisType == AxisType.Pressure) || (axisType == AxisType.GeoZ))
&& (positive != null)) {
axis.setPositive(positive);
axis.addAttribute(new Attribute(_Coordinate.ZisPositive, positive));
}
}
return axis;
}
Example 3
| Source File: CoverageCoordSys.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
public CoverageCoordAxis getZAxis() {
for (String axisName : getAxisNames()) {
CoverageCoordAxis axis = dataset.findCoordAxis(axisName);
if (axis == null)
throw new IllegalStateException("Cant find axis with name " + axisName);
if (axis.getAxisType() == AxisType.GeoZ || axis.getAxisType() == AxisType.Height
|| axis.getAxisType() == AxisType.Pressure)
return axis;
}
return null;
}
Example 4
| Source File: AWIPSConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
@Override
@Nullable
protected AxisType getAxisType(VariableDS.Builder v) {
String vname = v.shortName;
if (vname.equalsIgnoreCase("x"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("record"))
return AxisType.Time;
String dimName = v.getFirstDimensionName();
if ((dimName != null) && dimName.equalsIgnoreCase("record"))
return AxisType.Time;
String unit = v.getUnits();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
// otherwise guess
return AxisType.GeoZ;
}
Example 5
| Source File: WRFConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
@Override
@Nullable
protected AxisType getAxisType(VariableDS.Builder v) {
String vname = v.shortName;
if (vname.equalsIgnoreCase("x") || vname.equalsIgnoreCase("x_stag"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || vname.equalsIgnoreCase("y_stag"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("z") || vname.equalsIgnoreCase("z_stag"))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("Z"))
return AxisType.Height;
if (vname.equalsIgnoreCase("time") || vname.equalsIgnoreCase("times"))
return AxisType.Time;
String unit = v.getUnits();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
return null;
}
Example 6
| Source File: NUWGConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
protected AxisType getAxisType(NetcdfDataset ds, VariableEnhanced ve) {
Variable v = (Variable) ve;
String vname = v.getShortName();
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase(xaxisName))
return AxisType.GeoX;
if (vname.equalsIgnoreCase(yaxisName))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("record"))
return AxisType.Time;
Dimension dim = v.getDimension(0);
if ((dim != null) && dim.getShortName().equalsIgnoreCase("record")) { // wow thats bad!
return AxisType.Time;
}
String unit = ve.getUnitsString();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
if (SimpleUnit.isCompatible("sec", unit))
return null;
}
return AxisType.GeoZ; // AxisType.GeoZ;
}
Example 7
| Source File: NUWGConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
@Override
protected AxisType getAxisType(VariableDS.Builder v) {
String vname = v.shortName;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase(xaxisName))
return AxisType.GeoX;
if (vname.equalsIgnoreCase(yaxisName))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("record"))
return AxisType.Time;
String dimName = v.getFirstDimensionName();
if ((dimName != null) && dimName.equalsIgnoreCase("record")) { // wow thats bad!
return AxisType.Time;
}
String unit = v.getUnits();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
if (SimpleUnit.isCompatible("sec", unit))
return null;
}
return AxisType.GeoZ; // AxisType.GeoZ;
}
Example 8
| Source File: ADASConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
protected AxisType getAxisType(NetcdfDataset ds, VariableEnhanced ve) {
Variable v = (Variable) ve;
String vname = v.getShortName();
if (vname.equalsIgnoreCase("x") || vname.equalsIgnoreCase("x_stag"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || vname.equalsIgnoreCase("y_stag"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("z") || vname.equalsIgnoreCase("z_stag"))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("Z"))
return AxisType.Height;
if (vname.equalsIgnoreCase("time"))
return AxisType.Time;
String unit = ve.getUnitsString();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
return null;
}
Example 9
| Source File: GDVConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
@Override
protected AxisType getAxisType(VariableDS.Builder v) {
String vname = v.shortName;
if (vname.equalsIgnoreCase("x") || findAlias(v).equalsIgnoreCase("x"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon") || vname.equalsIgnoreCase("longitude") || findAlias(v).equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || findAlias(v).equalsIgnoreCase("y"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat") || vname.equalsIgnoreCase("latitude") || findAlias(v).equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("lev") || findAlias(v).equalsIgnoreCase("lev")
|| (vname.equalsIgnoreCase("level") || findAlias(v).equalsIgnoreCase("level")))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("z") || findAlias(v).equalsIgnoreCase("z")
|| (vname.equalsIgnoreCase("altitude") || vname.equalsIgnoreCase("depth")))
return AxisType.Height;
if (vname.equalsIgnoreCase("time") || findAlias(v).equalsIgnoreCase("time"))
return AxisType.Time;
return super.getAxisType(v);
}
Example 10
| Source File: AWIPSConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
protected AxisType getAxisType(NetcdfDataset ds, VariableEnhanced ve) {
Variable v = (Variable) ve;
String vname = v.getShortName();
if (vname.equalsIgnoreCase("x"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("record"))
return AxisType.Time;
Dimension dim = v.getDimension(0);
if ((dim != null) && dim.getShortName().equalsIgnoreCase("record"))
return AxisType.Time;
String unit = ve.getUnitsString();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
return AxisType.GeoZ;
}
Example 11
| Source File: WRFConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
protected AxisType getAxisType(NetcdfDataset ds, VariableEnhanced ve) {
Variable v = (Variable) ve;
String vname = v.getShortName();
if (vname.equalsIgnoreCase("x") || vname.equalsIgnoreCase("x_stag"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || vname.equalsIgnoreCase("y_stag"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("z") || vname.equalsIgnoreCase("z_stag"))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("Z"))
return AxisType.Height;
if (vname.equalsIgnoreCase("time") || vname.equalsIgnoreCase("times"))
return AxisType.Time;
String unit = ve.getUnitsString();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
return null;
}
Example 12
| Source File: DefaultConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
private AxisType getAxisTypeCoards(NetcdfDataset ncDataset, VariableEnhanced v) {
String unit = v.getUnitsString();
if (unit == null)
return null;
if (unit.equalsIgnoreCase("degrees_east") || unit.equalsIgnoreCase("degrees_E") || unit.equalsIgnoreCase("degreesE")
|| unit.equalsIgnoreCase("degree_east") || unit.equalsIgnoreCase("degree_E")
|| unit.equalsIgnoreCase("degreeE"))
return AxisType.Lon;
if (unit.equalsIgnoreCase("degrees_north") || unit.equalsIgnoreCase("degrees_N")
|| unit.equalsIgnoreCase("degreesN") || unit.equalsIgnoreCase("degree_north")
|| unit.equalsIgnoreCase("degree_N") || unit.equalsIgnoreCase("degreeN"))
return AxisType.Lat;
if (SimpleUnit.isDateUnit(unit)) // || SimpleUnit.isTimeUnit(unit)) removed dec 18, 2008
return AxisType.Time;
// look for other z coordinate
// if (SimpleUnit.isCompatible("m", unit))
// return AxisType.Height;
if (SimpleUnit.isCompatible("mbar", unit))
return AxisType.Pressure;
if (unit.equalsIgnoreCase("level") || unit.equalsIgnoreCase("layer") || unit.equalsIgnoreCase("sigma_level"))
return AxisType.GeoZ;
String positive = ncDataset.findAttValueIgnoreCase((Variable) v, "positive", null);
if (positive != null) {
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
else
return AxisType.GeoZ;
}
return null;
}
Example 13
| Source File: ADASConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 5 votes |
|
@Override
protected AxisType getAxisType(VariableDS.Builder vb) {
String vname = vb.shortName;
if (vname.equalsIgnoreCase("x") || vname.equalsIgnoreCase("x_stag"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || vname.equalsIgnoreCase("y_stag"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("z") || vname.equalsIgnoreCase("z_stag"))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("Z"))
return AxisType.Height;
if (vname.equalsIgnoreCase("time"))
return AxisType.Time;
String unit = vb.getUnits();
if (unit != null) {
if (SimpleUnit.isCompatible("millibar", unit))
return AxisType.Pressure;
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
}
return null;
}
Example 14
| Source File: COARDSConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
protected AxisType getAxisType(NetcdfDataset ncDataset, VariableEnhanced v) {
String unit = v.getUnitsString();
if (unit == null)
return null;
unit = unit.trim();
if (unit.equalsIgnoreCase("degrees_east") || unit.equalsIgnoreCase("degrees_E") || unit.equalsIgnoreCase("degreesE")
|| unit.equalsIgnoreCase("degree_east") || unit.equalsIgnoreCase("degree_E")
|| unit.equalsIgnoreCase("degreeE"))
return AxisType.Lon;
if (unit.equalsIgnoreCase("degrees_north") || unit.equalsIgnoreCase("degrees_N")
|| unit.equalsIgnoreCase("degreesN") || unit.equalsIgnoreCase("degree_north")
|| unit.equalsIgnoreCase("degree_N") || unit.equalsIgnoreCase("degreeN"))
return AxisType.Lat;
if (SimpleUnit.isDateUnit(unit)) {
return AxisType.Time;
}
// look for other z coordinate
if (SimpleUnit.isCompatible("mbar", unit))
return AxisType.Pressure;
if (unit.equalsIgnoreCase("level") || unit.equalsIgnoreCase("layer") || unit.equalsIgnoreCase("sigma_level"))
return AxisType.GeoZ;
String positive = ncDataset.findAttValueIgnoreCase((Variable) v, CF.POSITIVE, null);
if (positive != null) {
if (SimpleUnit.isCompatible("m", unit))
return AxisType.Height;
else
return AxisType.GeoZ;
}
// a bad idea, but CDC SST relies on it :
// :Source = "NOAA/National Climatic Data Center";
// :Contact = "Dick Reynolds, email: [email protected] & Chunying Liu, email: [email protected]";
// :netcdf_Convention = "COARDS";
// if (checkForMeter && SimpleUnit.isCompatible("m", unit))
// return AxisType.Height;
return null;
}
Example 15
| Source File: DefaultConventions.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
@Nullable
private AxisType getAxisTypeCoards(VariableDS.Builder vb) {
String unit = vb.getUnits();
if (unit == null) {
return null;
}
if (unit.equalsIgnoreCase("degrees_east") || unit.equalsIgnoreCase("degrees_E") || unit.equalsIgnoreCase("degreesE")
|| unit.equalsIgnoreCase("degree_east") || unit.equalsIgnoreCase("degree_E")
|| unit.equalsIgnoreCase("degreeE")) {
return AxisType.Lon;
}
if (unit.equalsIgnoreCase("degrees_north") || unit.equalsIgnoreCase("degrees_N")
|| unit.equalsIgnoreCase("degreesN") || unit.equalsIgnoreCase("degree_north")
|| unit.equalsIgnoreCase("degree_N") || unit.equalsIgnoreCase("degreeN")) {
return AxisType.Lat;
}
if (SimpleUnit.isDateUnit(unit)) // || SimpleUnit.isTimeUnit(unit)) removed dec 18, 2008
{
return AxisType.Time;
}
// look for other z coordinate
// if (SimpleUnit.isCompatible("m", unit))
// return AxisType.Height;
if (SimpleUnit.isCompatible("mbar", unit)) {
return AxisType.Pressure;
}
if (unit.equalsIgnoreCase("level") || unit.equalsIgnoreCase("layer") || unit.equalsIgnoreCase("sigma_level")) {
return AxisType.GeoZ;
}
String positive = vb.getAttributeContainer().findAttributeString("positive", null);
if (positive != null) {
if (SimpleUnit.isCompatible("m", unit)) {
return AxisType.Height;
} else {
return AxisType.GeoZ;
}
}
return null;
}
Example 16
| Source File: DefaultConvention.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
protected AxisType getAxisType(NetcdfDataset ds, VariableEnhanced ve) {
AxisType result = getAxisTypeCoards(ds, ve);
if (result != null)
return result;
Variable v = (Variable) ve;
String vname = v.getShortName();
if (vname == null)
return null;
String unit = v.getUnitsString();
if (unit == null)
unit = "";
String desc = v.getDescription();
if (desc == null)
desc = "";
if (vname.equalsIgnoreCase("x") || findAlias(ds, v).equalsIgnoreCase("x"))
return AxisType.GeoX;
if (vname.equalsIgnoreCase("lon") || vname.equalsIgnoreCase("longitude")
|| findAlias(ds, v).equalsIgnoreCase("lon"))
return AxisType.Lon;
if (vname.equalsIgnoreCase("y") || findAlias(ds, v).equalsIgnoreCase("y"))
return AxisType.GeoY;
if (vname.equalsIgnoreCase("lat") || vname.equalsIgnoreCase("latitude") || findAlias(ds, v).equalsIgnoreCase("lat"))
return AxisType.Lat;
if (vname.equalsIgnoreCase("lev") || findAlias(ds, v).equalsIgnoreCase("lev")
|| (vname.equalsIgnoreCase("level") || findAlias(ds, v).equalsIgnoreCase("level")))
return AxisType.GeoZ;
if (vname.equalsIgnoreCase("z") || findAlias(ds, v).equalsIgnoreCase("z") || vname.equalsIgnoreCase("altitude")
|| desc.contains("altitude") || vname.equalsIgnoreCase("depth") || vname.equalsIgnoreCase("elev")
|| vname.equalsIgnoreCase("elevation")) {
if (SimpleUnit.isCompatible("m", unit)) // units of meters
return AxisType.Height;
}
if (vname.equalsIgnoreCase("time") || findAlias(ds, v).equalsIgnoreCase("time")) {
if (SimpleUnit.isDateUnit(unit))
return AxisType.Time;
}
if (vname.equalsIgnoreCase("time") && v.getDataType() == DataType.STRING) {
try {
Array firstValue = v.read("0");
if (firstValue instanceof ArrayObject.D1) {
ArrayObject.D1 sarry = (ArrayObject.D1) firstValue;
String firstStringValue = (String) sarry.get(0);
if (CalendarDate.parseISOformat(null, firstStringValue) != null) // valid iso date string
return AxisType.Time;
}
} catch (IOException | InvalidRangeException e) {
logger.warn("time string error", e);
}
}
return null;
}
Example 17
| Source File: CoordinateSystem.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
protected CoordinateSystem(Builder<?> builder, NetcdfDataset ncd, List<CoordinateAxis> axes,
List<CoordinateTransform> allTransforms) {
this.ds = ncd;
this.isImplicit = builder.isImplicit;
// find referenced coordinate axes
List<CoordinateAxis> axesList = new ArrayList<>();
StringTokenizer stoker = new StringTokenizer(builder.coordAxesNames);
while (stoker.hasMoreTokens()) {
String vname = stoker.nextToken();
for (CoordinateAxis a : axes) {
String aname = a.getFullName();
if (aname.equals(vname)) {
axesList.add(a);
}
}
}
this.coordAxes = axesList;
// calculated
this.name = makeName(coordAxes);
for (CoordinateAxis axis : this.coordAxes) {
// look for AxisType
AxisType axisType = axis.getAxisType();
if (axisType != null) {
if (axisType == AxisType.GeoX)
xAxis = lesserRank(xAxis, axis);
if (axisType == AxisType.GeoY)
yAxis = lesserRank(yAxis, axis);
if (axisType == AxisType.GeoZ)
zAxis = lesserRank(zAxis, axis);
if (axisType == AxisType.Time)
tAxis = lesserRank(tAxis, axis);
if (axisType == AxisType.Lat)
latAxis = lesserRank(latAxis, axis);
if (axisType == AxisType.Lon)
lonAxis = lesserRank(lonAxis, axis);
if (axisType == AxisType.Height)
hAxis = lesserRank(hAxis, axis);
if (axisType == AxisType.Pressure)
pAxis = lesserRank(pAxis, axis);
if (axisType == AxisType.Ensemble)
ensAxis = lesserRank(ensAxis, axis);
if (axisType == AxisType.RadialAzimuth)
aziAxis = lesserRank(aziAxis, axis);
if (axisType == AxisType.RadialDistance)
radialAxis = lesserRank(radialAxis, axis);
if (axisType == AxisType.RadialElevation)
elevAxis = lesserRank(elevAxis, axis);
}
// collect dimensions
List<Dimension> dims = axis.getDimensionsAll();
domain.addAll(dims);
}
// Find the named coordinate transforms in allTransforms.
for (String wantTransName : builder.transNames) {
CoordinateTransform got = allTransforms.stream()
.filter(ct -> (wantTransName.equals(ct.getName())
|| ct.getAttributeContainer() != null && wantTransName.equals(ct.getAttributeContainer().getName())))
.findFirst().orElse(null);
if (got != null) {
coordTrans.add(got);
}
}
}
Example 18
| Source File: DefaultConventions.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
@Override
@Nullable
protected AxisType getAxisType(VariableDS.Builder vb) {
AxisType result = getAxisTypeCoards(vb);
if (result != null) {
return result;
}
String vname = vb.shortName;
if (vname == null) {
return null;
}
String unit = vb.getUnits();
if (unit == null) {
unit = "";
}
String desc = vb.getDescription();
if (desc == null) {
desc = "";
}
if (vname.equalsIgnoreCase("x") || findAlias(vb).equalsIgnoreCase("x")) {
return AxisType.GeoX;
}
if (vname.equalsIgnoreCase("lon") || vname.equalsIgnoreCase("longitude") || findAlias(vb).equalsIgnoreCase("lon")) {
return AxisType.Lon;
}
if (vname.equalsIgnoreCase("y") || findAlias(vb).equalsIgnoreCase("y")) {
return AxisType.GeoY;
}
if (vname.equalsIgnoreCase("lat") || vname.equalsIgnoreCase("latitude") || findAlias(vb).equalsIgnoreCase("lat")) {
return AxisType.Lat;
}
if (vname.equalsIgnoreCase("lev") || findAlias(vb).equalsIgnoreCase("lev")
|| (vname.equalsIgnoreCase("level") || findAlias(vb).equalsIgnoreCase("level"))) {
return AxisType.GeoZ;
}
if (vname.equalsIgnoreCase("z") || findAlias(vb).equalsIgnoreCase("z") || vname.equalsIgnoreCase("altitude")
|| desc.contains("altitude") || vname.equalsIgnoreCase("depth") || vname.equalsIgnoreCase("elev")
|| vname.equalsIgnoreCase("elevation")) {
if (SimpleUnit.isCompatible("m", unit)) // units of meters
{
return AxisType.Height;
}
}
if (vname.equalsIgnoreCase("time") || findAlias(vb).equalsIgnoreCase("time")) {
if (SimpleUnit.isDateUnit(unit)) {
return AxisType.Time;
}
}
if (vname.equalsIgnoreCase("time") && vb.dataType == DataType.STRING) {
if (vb.orgVar != null) {
try {
Array firstValue = vb.orgVar.read("0");
if (firstValue instanceof ArrayObject.D1) {
ArrayObject.D1 sarry = (ArrayObject.D1) firstValue;
String firstStringValue = (String) sarry.get(0);
if (CalendarDate.parseISOformat(null, firstStringValue) != null) { // valid iso date string LOOK
return AxisType.Time;
}
}
} catch (IOException | InvalidRangeException e) {
logger.warn("time string error", e);
}
}
}
return null;
}
Example 19
| Source File: CoordinateSystem.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
/**
* Constructor.
*
* @param ds the containing dataset
* @param axes Collection of type CoordinateAxis, must be at least one.
* @param coordTrans Collection of type CoordinateTransform, may be empty or null.
* @deprecated Use CoordinateSystem.builder()
*/
@Deprecated
public CoordinateSystem(NetcdfDataset ds, Collection<CoordinateAxis> axes,
Collection<CoordinateTransform> coordTrans) {
this.ds = ds;
this.coordAxes = new ArrayList<>(axes);
this.name = makeName(coordAxes);
if (coordTrans != null)
this.coordTrans = new ArrayList<>(coordTrans);
for (CoordinateAxis axis : coordAxes) {
// look for AxisType
AxisType axisType = axis.getAxisType();
if (axisType != null) {
if (axisType == AxisType.GeoX)
xAxis = lesserRank(xAxis, axis);
if (axisType == AxisType.GeoY)
yAxis = lesserRank(yAxis, axis);
if (axisType == AxisType.GeoZ)
zAxis = lesserRank(zAxis, axis);
if (axisType == AxisType.Time)
tAxis = lesserRank(tAxis, axis);
if (axisType == AxisType.Lat)
latAxis = lesserRank(latAxis, axis);
if (axisType == AxisType.Lon)
lonAxis = lesserRank(lonAxis, axis);
if (axisType == AxisType.Height)
hAxis = lesserRank(hAxis, axis);
if (axisType == AxisType.Pressure)
pAxis = lesserRank(pAxis, axis);
if (axisType == AxisType.Ensemble)
ensAxis = lesserRank(ensAxis, axis);
if (axisType == AxisType.RadialAzimuth)
aziAxis = lesserRank(aziAxis, axis);
if (axisType == AxisType.RadialDistance)
radialAxis = lesserRank(radialAxis, axis);
if (axisType == AxisType.RadialElevation)
elevAxis = lesserRank(elevAxis, axis);
}
// collect dimensions
List<Dimension> dims = axis.getDimensionsAll();
domain.addAll(dims);
}
}
Example 20
| Source File: GridCoordSys.java From netcdf-java with BSD 3-Clause "New" or "Revised" License | 4 votes |
|
/**
* Add this coordinate system to the netCDF file
*
* @param ncfile the netCDF file
* @param g the group to add to
*/
void addToNetcdfFile(NetcdfFile ncfile, Group g) {
if (dontUseVertical) {
return;
}
if (g == null) {
g = ncfile.getRootGroup();
}
String dims = "time";
if (!dontUseVertical) {
dims = dims + " " + verticalName;
}
if (hcs.isLatLon()) {
dims = dims + " lat lon";
} else {
dims = dims + " y x";
}
// Collections.sort( levels);
int nlevs = levels.size();
// ncfile.addDimension(g, new Dimension(verticalName, nlevs, true));
// coordinate axis and coordinate system Variable
Variable v = new Variable(ncfile, g, null, verticalName);
v.setDataType(DataType.DOUBLE);
v.addAttribute(new Attribute("long_name", lookup.getLevelDescription(record)));
v.addAttribute(new Attribute("units", lookup.getLevelUnit(record)));
// positive attribute needed for CF-1 Height and Pressure
if (positive != null) {
v.addAttribute(new Attribute("positive", positive));
}
if (units != null) {
AxisType axisType;
if (SimpleUnit.isCompatible("millibar", units)) {
axisType = AxisType.Pressure;
} else if (SimpleUnit.isCompatible("m", units)) {
axisType = AxisType.Height;
} else {
axisType = AxisType.GeoZ;
}
v.addAttribute(new Attribute("grid_level_type", Integer.toString(record.getLevelType1())));
v.addAttribute(new Attribute(_Coordinate.AxisType, axisType.toString()));
v.addAttribute(new Attribute(_Coordinate.Axes, dims));
if (!hcs.isLatLon()) {
v.addAttribute(new Attribute(_Coordinate.Transforms, hcs.getGridName()));
}
}
double[] data = new double[nlevs];
for (int i = 0; i < levels.size(); i++) {
Double d = levels.get(i);
data[i] = d;
}
Array dataArray = Array.factory(DataType.DOUBLE, new int[] {nlevs}, data);
v.setDimensions(verticalName);
v.setCachedData(dataArray, false);
ncfile.addVariable(g, v);
// look for vertical transforms
if (record.getLevelType1() == 109) {
findCoordinateTransform(g, "Pressure", record.getLevelType1());
}
}
