private Element getStationPeriod(CalendarDateRange stationDateRange) {
Element retval = getDocument().createElement("sos:time");
if (stationDateRange != null) {
// time
Element timePeriod = getDocument().createElement("gml:TimePeriod");
timePeriod.setAttribute("xsi:type", "gml:TimePeriodType");
// begin
Element begin = getDocument().createElement("gml:beginPosition");
begin.setTextContent(stationDateRange.getStart().toString());
timePeriod.appendChild(begin);
// end
Element end = getDocument().createElement("gml:endPosition");
end.setTextContent(stationDateRange.getEnd().toString());
timePeriod.appendChild(end);
retval.appendChild(timePeriod);
}
return retval;
}
public File writeCoverageDataToFile(
WcsRequest.Format format,
LatLonRect bboxLatLonRect,
AxisSubset vertSubset,
List<String> rangeSubset,
CalendarDateRange timeRange)
throws WcsException {
boolean zRangeDone = false;
boolean tRangeDone = false;
try {
// Get the height range.
Range zRange = vertSubset != null ? vertSubset.getRange() : null;
zRangeDone = true;
// Get the time range.
Range tRange = null;
if (timeRange != null) {
CoordinateAxis1DTime timeAxis = this.coordSys.getTimeAxis1D();
int startIndex = timeAxis.findTimeIndexFromCalendarDate(timeRange.getStart());
int endIndex = timeAxis.findTimeIndexFromCalendarDate(timeRange.getEnd());
tRange = new Range(startIndex, endIndex);
tRangeDone = true;
}
if (format == WcsRequest.Format.GeoTIFF || format == WcsRequest.Format.GeoTIFF_Float) {
if (rangeSubset.size() != 1) {
String msg =
"GeoTIFF response encoding only available for single range field selection ["
+ rangeSubset
+ "].";
log.error("writeCoverageDataToFile(): " + msg);
throw new WcsException(WcsException.Code.InvalidParameterValue, "RangeSubset", msg);
}
String reqRangeFieldName = rangeSubset.get(0);
File dir = new File(getDiskCache().getRootDirectory());
File tifFile = File.createTempFile("WCS", ".tif", dir);
if (log.isDebugEnabled())
log.debug("writeCoverageDataToFile(): tifFile=" + tifFile.getPath());
WcsRangeField rangeField = this.range.get(reqRangeFieldName);
GridDatatype subset =
rangeField.getGridDatatype().makeSubset(tRange, zRange, bboxLatLonRect, 1, 1, 1);
Array data = subset.readDataSlice(0, 0, -1, -1);
GeotiffWriter writer = new GeotiffWriter(tifFile.getPath());
writer.writeGrid(
this.dataset.getDataset(), subset, data, format == WcsRequest.Format.GeoTIFF);
writer.close();
return tifFile;
} else if (format == WcsRequest.Format.NetCDF3) {
File dir = new File(getDiskCache().getRootDirectory());
File outFile = File.createTempFile("WCS", ".nc", dir);
if (log.isDebugEnabled())
log.debug("writeCoverageDataToFile(): ncFile=" + outFile.getPath());
// WTF ?? this.coordSys.getVerticalAxis().isNumeric();
NetcdfFileWriter writer =
NetcdfFileWriter.createNew(NetcdfFileWriter.Version.netcdf3, outFile.getAbsolutePath());
CFGridWriter2.writeFile(
this.dataset.getDataset(),
rangeSubset,
bboxLatLonRect,
null,
1,
zRange,
timeRange,
1,
true,
writer);
return outFile;
} else {
log.error(
"writeCoverageDataToFile(): Unsupported response encoding format [" + format + "].");
throw new WcsException(
WcsException.Code.InvalidFormat,
"Format",
"Unsupported response encoding format [" + format + "].");
}
} catch (InvalidRangeException e) {
String msg = "Failed to subset coverage [" + this.getName();
if (!zRangeDone) msg += "] along vertical axis [" + vertSubset + "]. ";
else if (!tRangeDone) msg += "] along time axis [" + timeRange + "]. ";
else msg += "] in horizontal plane [" + bboxLatLonRect + "]. ";
log.error("writeCoverageDataToFile(): " + msg + e.getMessage());
throw new WcsException(WcsException.Code.CoverageNotDefined, "", msg);
} catch (IOException e) {
log.error(
"writeCoverageDataToFile(): Failed to write file for requested coverage <"
+ this.getName()
+ ">: "
+ e.getMessage());
throw new WcsException(
WcsException.Code.UNKNOWN, "", "Problem creating coverage [" + this.getName() + "].");
}
}
private void datasetInfoHtml(
RadarType radarType, String pathInfo, PrintWriter pw, HttpServletResponse res)
throws IOException {
pathInfo = pathInfo.replace("/dataset.html", "");
pathInfo = pathInfo.replace("/catalog.html", "");
Element root = new Element("RadarNexrad");
Document doc = new Document(root);
if (pathInfo.startsWith("/")) pathInfo = pathInfo.substring(1);
for (int i = 0; i < datasets.size(); i++) {
InvDatasetScan ds = (InvDatasetScan) datasets.get(i);
if (!(pathInfo.equals(ds.getPath()))) {
continue;
}
// at this point a valid dataset
// fix the location
root.setAttribute("location", "/thredds/radarServer/" + ds.getPath());
// spatial range
ThreddsMetadata.GeospatialCoverage gc = ds.getGeospatialCoverage();
LatLonRect bb = new LatLonRect();
gc.setBoundingBox(bb);
String north = Double.toString(gc.getLatNorth());
String south = Double.toString(gc.getLatSouth());
String east = Double.toString(gc.getLonEast());
String west = Double.toString(gc.getLonWest());
Element LatLonBox = new Element("LatLonBox");
LatLonBox.addContent(new Element("north").addContent(north));
LatLonBox.addContent(new Element("south").addContent(south));
LatLonBox.addContent(new Element("east").addContent(east));
LatLonBox.addContent(new Element("west").addContent(west));
root.addContent(LatLonBox);
// get the time range
Element timeSpan = new Element("TimeSpan");
CalendarDateRange dr = ds.getCalendarDateCoverage();
timeSpan.addContent(new Element("begin").addContent(dr.getStart().toString()));
timeSpan.addContent(new Element("end").addContent(dr.getEnd().toString()));
root.addContent(timeSpan);
ThreddsMetadata.Variables cvs = (ThreddsMetadata.Variables) ds.getVariables().get(0);
List vl = cvs.getVariableList();
for (int j = 0; j < vl.size(); j++) {
ThreddsMetadata.Variable v = (ThreddsMetadata.Variable) vl.get(j);
Element variable = new Element("variable");
variable.setAttribute("name", v.getName());
root.addContent(variable);
}
// add pointer to the station list XML
/*
Element stnList = new Element("stationList");
stnList.setAttribute("title", "Available Stations", XMLEntityResolver.xlinkNS);
stnList.setAttribute("href", "/thredds/radarServer/"+ pathInfo +"/stations.xml",
XMLEntityResolver.xlinkNS);
root.addContent(stnList);
*/
// String[] stations = rns.stationsDS( dataLocation.get( ds.getPath() ));
// rns.printStations( stations );
// add accept list
Element a = new Element("AcceptList");
a.addContent(new Element("accept").addContent("xml"));
a.addContent(new Element("accept").addContent("html"));
root.addContent(a);
}
ServerMethods sm = new ServerMethods(log);
InputStream xslt = sm.getInputStream(contentPath + getPath() + "radar.xsl", RadarServer.class);
try {
// what's wrong here xslt = getXSLT( "radar.xsl" );
XSLTransformer transformer = new XSLTransformer(xslt);
Document html = transformer.transform(doc);
XMLOutputter fmt = new XMLOutputter(Format.getPrettyFormat());
String infoString = fmt.outputString(html);
res.setContentType("text/html; charset=iso-8859-1");
pw = res.getWriter();
pw.println(infoString);
pw.flush();
} catch (Exception e) {
log.error("radarServer reading " + contentPath + getPath() + "radar.xsl");
log.error("radarServer XSLTransformer problem for web form ", e);
} finally {
if (xslt != null) {
try {
xslt.close();
} catch (IOException e) {
log.error("radarServer radar.xsl: error closing" + contentPath + getPath() + "radar.xsl");
}
}
}
return;
}
private void datasetInfoXml(RadarType radarType, String pathInfo, PrintWriter pw)
throws IOException {
try {
pw.println(
"<catalog xmlns=\"http://www.unidata.ucar.edu/namespaces/thredds/InvCatalog/v1.0\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" name=\"Radar Data\" version=\"1.0.1\">\n");
// add service
pw.println(
" <service name=\"radarServer\" base=\"/thredds/radarServer/\" serviceType=\"DQC\" />\n");
pathInfo = pathInfo.replace("/dataset.xml", "");
pathInfo = pathInfo.replace("/catalog.xml", "");
if (pathInfo.startsWith("/")) pathInfo = pathInfo.substring(1);
for (int i = 0; i < datasets.size(); i++) {
InvDatasetScan ds = (InvDatasetScan) datasets.get(i);
if (!(pathInfo.equals(ds.getPath()))) {
continue;
}
pw.println(" <dataset ID=\"" + ds.getID() + "\" serviceName=\"radarServer\">");
pw.println(" <urlpath>" + ds.getPath() + "</urlpath>");
pw.println(" <dataType>" + ds.getDataType() + "</dataType>");
pw.println(" <dataFormat>" + ds.getDataFormatType() + "</dataFormat>");
pw.println(" <serviceName>radarServer</serviceName>");
pw.println(" <metadata inherited=\"true\">");
pw.println(" <documentation type=\"summary\">" + ds.getSummary() + "</documentation>");
CalendarDateRange dr = ds.getCalendarDateCoverage();
pw.println(" <TimeSpan>");
pw.print(" <start>");
if (pathInfo.contains("IDD")) {
pw.print(rm.getStartDateTime(ds.getPath()));
} else {
pw.print(dr.getStart().toString());
}
pw.println("</start>");
pw.println(" <end>" + dr.getEnd().toString() + "</end>");
pw.println(" </TimeSpan>");
ThreddsMetadata.GeospatialCoverage gc = ds.getGeospatialCoverage();
LatLonRect bb = new LatLonRect();
gc.setBoundingBox(bb);
pw.println(" <LatLonBox>");
pw.println(" <north>" + gc.getLatNorth() + "</north>");
pw.println(" <south>" + gc.getLatSouth() + "</south>");
pw.println(" <east>" + gc.getLonEast() + "</east>");
pw.println(" <west>" + gc.getLonWest() + "</west>");
pw.println(" </LatLonBox>");
ThreddsMetadata.Variables cvs = (ThreddsMetadata.Variables) ds.getVariables().get(0);
List vl = cvs.getVariableList();
pw.println(" <Variables>");
for (int j = 0; j < vl.size(); j++) {
ThreddsMetadata.Variable v = (ThreddsMetadata.Variable) vl.get(j);
pw.println(
" <variable name=\""
+ v.getName()
+ "\" vocabulary_name=\""
+ v.getVocabularyName()
+ "\" units=\""
+ v.getUnits()
+ "\" />");
}
pw.println(" </Variables>");
String[] stations = rm.stationsDS(radarType, dataLocation.get(ds.getPath()));
rm.printStations(stations, pw, radarType);
pw.println(" </metadata>");
pw.println(" </dataset>");
}
pw.println("</catalog>");
pw.flush();
} catch (Throwable e) {
log.error("RadarServer.datasetInfoXml", e);
}
return;
}