@Override
protected StationHelper createStationHelper() throws IOException {
TimedCollection.Dataset td = dataCollection.getPrototype();
if (td == null) throw new RuntimeException("No datasets in the collection");
Formatter errlog = new Formatter();
try (FeatureDatasetPoint openDataset =
(FeatureDatasetPoint)
FeatureDatasetFactoryManager.open(
FeatureType.STATION, td.getLocation(), null, errlog)) {
StationHelper stationHelper = new StationHelper();
List<FeatureCollection> fcList = openDataset.getPointFeatureCollectionList();
StationTimeSeriesCollectionImpl openCollection =
(StationTimeSeriesCollectionImpl) fcList.get(0);
List<Station> stns = openCollection.getStations();
for (Station s : stns) {
StationTimeSeriesFeature stnFeature = openCollection.getStationFeature(s);
stationHelper.addStation(
new CompositeStationFeature(
s, timeUnit, altUnits, stnFeature.getFeatureData(), this.dataCollection));
}
dataVariables = openDataset.getDataVariables();
globalAttributes = openDataset.getGlobalAttributes();
return stationHelper;
}
}
/**
* Open a ucar.nc2.ft.PointFeatureCollection, write out in CF point format.
*
* @param fileIn open through TypedDatasetFactory.open(FeatureType.POINT, ..)
* @param fileOut write to this netcdf-3 file
* @param inMemory if true, read file into memory for efficiency
* @return true on success
* @throws IOException on read/write error
*/
public static boolean rewritePointFeatureDataset(String fileIn, String fileOut, boolean inMemory)
throws IOException {
System.out.println(
"Rewrite2 .nc files from " + fileIn + " to " + fileOut + " inMemory= " + inMemory);
long start = System.currentTimeMillis();
// do it in memory for speed
NetcdfFile ncfile = inMemory ? NetcdfFile.openInMemory(fileIn) : NetcdfFile.open(fileIn);
NetcdfDataset ncd = new NetcdfDataset(ncfile);
Formatter errlog = new Formatter();
FeatureDataset fd = FeatureDatasetFactoryManager.wrap(FeatureType.ANY_POINT, ncd, null, errlog);
if (fd == null) return false;
if (fd instanceof FeatureDatasetPoint) {
writePointFeatureCollection((FeatureDatasetPoint) fd, fileOut);
fd.close();
long took = System.currentTimeMillis() - start;
System.out.println(" that took " + (took - start) + " msecs");
return true;
}
return false;
}
private PointFeatureIterator getNextIterator() throws IOException {
if (!iter.hasNext()) return null;
TimedCollection.Dataset td = iter.next();
Formatter errlog = new Formatter();
currentDataset =
(FeatureDatasetPoint)
FeatureDatasetFactoryManager.open(
FeatureType.STATION, td.getLocation(), null, errlog);
List<FeatureCollection> fcList = currentDataset.getPointFeatureCollectionList();
StationTimeSeriesFeatureCollection stnCollection =
(StationTimeSeriesFeatureCollection) fcList.get(0);
Station s = stnCollection.getStation(getName());
if (s == null) {
System.out.printf(
"CompositeStationFeatureIterator dataset: %s missing station %s%n",
td.getLocation(), getName());
return getNextIterator();
}
StationTimeSeriesFeature stnFeature = stnCollection.getStationFeature(s);
if (CompositeDatasetFactory.debug)
System.out.printf(
"CompositeStationFeatureIterator open dataset: %s for %s%n",
td.getLocation(), s.getName());
return stnFeature.getPointFeatureIterator(bufferSize);
}
