in endorsed/src/org.apache.sis.storage/main/org/apache/sis/storage/csv/Store.java [340:482]
private GeneralEnvelope parseEnvelope(final List<String> elements) throws DataStoreException, FactoryException {
@SuppressWarnings("LocalVariableHidesMemberVariable")
int spatialDimensionCount = 2; // Another result of this method to be computed as a side-effect.
CoordinateReferenceSystem crs = null;
boolean isDimExplicit = false;
double[] lowerCorner = ArraysExt.EMPTY_DOUBLE;
double[] upperCorner = ArraysExt.EMPTY_DOUBLE;
Instant startTime = null;
Instant endTime = null;
Unit<Time> timeUnit = Units.SECOND;
boolean isTimeAbsolute = false;
int ordinal = -1;
for (final String element : elements) {
ordinal++;
if (!element.isEmpty()) {
switch (ordinal) {
case 0: continue; // The "@stboundedby" header.
case 1: crs = CRS.forCode(element); continue;
case 2: if (element.length() == 2 && Character.toUpperCase(element.charAt(1)) == 'D') {
isDimExplicit = true;
spatialDimensionCount = element.charAt(0) - '0';
if (spatialDimensionCount < 1 || spatialDimensionCount > 3) {
throw new DataStoreReferencingException(errors().getString(
Errors.Keys.IllegalCoordinateSystem_1, element));
}
continue;
}
/*
* According the Moving Feature specification, the [dim] element is optional.
* If we did not recognized the dimension, assume that we have the next element
* (i.e. the lower corner). Fall-through so we can process it.
*/
ordinal++; // Fall through
case 3: lowerCorner = CharSequences.parseDoubles(element, ORDINATE_SEPARATOR); continue;
case 4: upperCorner = CharSequences.parseDoubles(element, ORDINATE_SEPARATOR); continue;
case 5: startTime = LenientDateFormat.parseInstantUTC(element); continue;
case 6: endTime = LenientDateFormat.parseInstantUTC(element); continue;
case 7: switch (element.toLowerCase(Locale.US)) {
case "sec":
case "second": /* Already SECOND. */ continue;
case "minute": timeUnit = Units.MINUTE; continue;
case "hour": timeUnit = Units.HOUR; continue;
case "day": timeUnit = Units.DAY; continue;
case "absolute": isTimeAbsolute = true; continue;
default: throw new DataStoreReferencingException(errors().getString(Errors.Keys.UnknownUnit_1, element));
}
}
// If we reach this point, there is some remaining unknown elements. Ignore them.
break;
}
}
/*
* Complete the CRS by adding a vertical component if needed, then a temporal component.
* Only after the CRS has been completed we can create the envelope.
*
* Vertical component:
* Ideally, should be part of the CRS created from the authority code. But if the authority
* code is only for a two-dimensional CRS, we default to an arbitrary height component.
*
* Temporal component:
* Assumed never part of the authority code. We need to build the temporal component ourselves
* in order to set the origin to the start time.
*/
@SuppressWarnings("LocalVariableHidesMemberVariable")
final GeneralEnvelope envelope; // Value will be assigned to `this.envelope` by the caller.
if (crs != null) {
int count = 0;
final CoordinateReferenceSystem[] components = new CoordinateReferenceSystem[3];
components[count++] = crs;
/*
* If the coordinates are three-dimensional but the CRS is 2D, add a vertical axis.
* The vertical axis shall be the third one, however we do not enforce that rule
* since Apache SIS should work correctly even if the vertical axis is elsewhere.
*/
int dimension = crs.getCoordinateSystem().getDimension();
if (isDimExplicit) {
if (spatialDimensionCount > dimension) {
components[count++] = CommonCRS.Vertical.MEAN_SEA_LEVEL.crs();
dimension++;
}
if (dimension != spatialDimensionCount) {
throw new DataStoreReferencingException(errors().getString(
Errors.Keys.MismatchedDimension_3, "@stboundedby(CRS)", spatialDimensionCount, dimension));
}
}
if (dimension >= Numerics.MAXIMUM_MATRIX_SIZE) {
throw new DataStoreReferencingException(errors().getString(
Errors.Keys.ExcessiveNumberOfDimensions_1, dimension));
}
spatialDimensionCount = dimension;
/*
* Add a temporal axis if we have a start time (no need for end time).
* This block presumes that the CRS does not already have a time axis.
* If a time axis was already present, an exception will be thrown at
* builder.createCompoundCRS(…) invocation time.
*/
final GeodeticObjectBuilder builder = new GeodeticObjectBuilder();
String name = crs.getName().getCode();
if (startTime != null) {
final TemporalCRS temporal;
if (isTimeAbsolute) {
temporal = TimeEncoding.DEFAULT.crs();
timeEncoding = TimeEncoding.ABSOLUTE;
} else {
temporal = builder.createTemporalCRS(startTime, timeUnit);
timeEncoding = new TimeEncoding(PseudoDatum.of(temporal), timeUnit);
}
components[count++] = temporal;
name = name + " + " + temporal.getName().getCode();
}
crs = builder.addName(name).createCompoundCRS(ArraysExt.resize(components, count));
envelope = new GeneralEnvelope(crs);
} else {
/*
* While illegal in principle, Apache SIS accepts missing CRS.
* In such case, use only the number of dimensions.
*/
int dim = spatialDimensionCount;
if (startTime != null) dim++; // Same criterion as in above block.
envelope = new GeneralEnvelope(dim);
}
/*
* At this point we got the three- or four-dimensional spatiotemporal CRS.
* We can now set the envelope coordinate values, including temporal values.
*/
int dim;
if ((dim = lowerCorner.length) != spatialDimensionCount ||
(dim = upperCorner.length) != spatialDimensionCount)
{
throw new DataStoreReferencingException(errors().getString(
Errors.Keys.MismatchedDimension_3, "@stboundedby(BBOX)", spatialDimensionCount, dim));
}
for (int i=0; i<spatialDimensionCount; i++) {
envelope.setRange(i, lowerCorner[i], upperCorner[i]);
}
if (startTime != null) {
envelope.setRange(spatialDimensionCount, timeEncoding.toCRS(startTime),
(endTime == null) ? Double.NaN : timeEncoding.toCRS(endTime));
}
this.spatialDimensionCount = (short) spatialDimensionCount;
return envelope;
}