public Map<Band, Variable> addBands(
Product product, Variable idxVariable, List<Variable> l3ProdVars) {
final Structure binListStruc = (Structure) idxVariable;
final Map<Band, Variable> bandToVariableMap = new HashMap<Band, Variable>();
// bandToVariableMap.put(addBand(product, "bin_num", ProductData.TYPE_UINT32),
// binListStruc.select("bin_num").findVariable("bin_num"));
bandToVariableMap.put(
addBand(product, "weights", ProductData.TYPE_FLOAT32),
binListStruc.select("weights").findVariable("weights"));
bandToVariableMap.put(
addBand(product, "nobs", ProductData.TYPE_UINT16),
binListStruc.select("nobs").findVariable("nobs"));
bandToVariableMap.put(
addBand(product, "nscenes", ProductData.TYPE_UINT16),
binListStruc.select("nscenes").findVariable("nscenes"));
// ncFile.getRootGroup().findGroup("Level-3 Binned Data").findVariable("BinList");
if (ncFile.getRootGroup().findGroup("Level-3_Binned_Data").findVariable("qual_l3") != null) {
bandToVariableMap.put(
addBand(product, "qual_l3", ProductData.TYPE_UINT8),
ncFile.getRootGroup().findGroup("Level-3_Binned_Data").findVariable("qual_l3"));
}
String groupnames = "";
for (Variable l3Var : l3ProdVars) {
String varName = l3Var.getShortName();
final int dataType = ProductData.TYPE_FLOAT32;
if (!varName.contains("Bin")
&& (!varName.startsWith("qual"))
&& (!varName.equalsIgnoreCase("SEAGrid"))
&& (!varName.equalsIgnoreCase("Input_Files"))) {
final Structure binStruc = (Structure) l3Var;
if (groupnames.length() == 0) {
groupnames = varName;
} else {
groupnames = groupnames + ":" + varName;
}
List<String> vnames = binStruc.getVariableNames();
for (String bandvar : vnames) {
bandToVariableMap.put(
addBand(product, bandvar, dataType), binStruc.select(bandvar).findVariable(bandvar));
}
// Add virtual band for product mean
StringBuilder prodname = new StringBuilder(varName);
prodname.append("_mean");
String calcmean = ComputeBinMeans(varName);
Band varmean =
new VirtualBand(
prodname.toString(),
ProductData.TYPE_FLOAT32,
product.getSceneRasterWidth(),
product.getSceneRasterHeight(),
calcmean);
varmean.setNoDataValue(Double.NaN);
varmean.setNoDataValueUsed(true);
product.addBand(varmean);
// Add virtual band for product stdev
int underscore = prodname.indexOf("_mean");
prodname.delete(underscore, underscore + 5);
prodname.append("_stdev");
String calcstdev = ComputeBinVariances(varName);
Band varstdev =
new VirtualBand(
prodname.toString(),
ProductData.TYPE_FLOAT32,
product.getSceneRasterWidth(),
product.getSceneRasterHeight(),
calcstdev);
varstdev.setNoDataValue(Double.NaN);
varstdev.setNoDataValueUsed(true);
product.addBand(varstdev);
}
}
product.setAutoGrouping(groupnames);
return bandToVariableMap;
}
private String ComputeBinVariances(String prodname) {
StringBuilder bin_stdev = new StringBuilder("weights * weights <= nscenes ? 0.0 : sqrt((((");
bin_stdev.append(prodname);
bin_stdev.append("_sum_sq/weights) - (");
bin_stdev.append(prodname);
bin_stdev.append("_sum /weights)*(");
bin_stdev.append(prodname);
bin_stdev.append("_sum /weights))");
bin_stdev.append("* weights * weights) / (weights * weights - nscenes))");
return bin_stdev.toString();
}
private boolean checkReferencedRastersIncluded() {
final Set<String> notIncludedNames = new TreeSet<String>();
final List<String> includedNodeNames = Arrays.asList(productSubsetDef.getNodeNames());
for (final String nodeName : includedNodeNames) {
final RasterDataNode rasterDataNode = product.getRasterDataNode(nodeName);
if (rasterDataNode != null) {
collectNotIncludedReferences(rasterDataNode, notIncludedNames);
}
}
boolean ok = true;
if (!notIncludedNames.isEmpty()) {
StringBuilder nameListText = new StringBuilder();
for (String notIncludedName : notIncludedNames) {
nameListText.append(" '").append(notIncludedName).append("'\n");
}
final String pattern =
"The following dataset(s) are referenced but not included\n"
+ "in your current subset definition:\n"
+ "{0}\n"
+ "If you do not include these dataset(s) into your selection,\n"
+ "you might get unexpected results while working with the\n"
+ "resulting product.\n\n"
+ "Do you wish to include the referenced dataset(s) into your\n"
+ "subset definition?\n"; /*I18N*/
final MessageFormat format = new MessageFormat(pattern);
int status =
JOptionPane.showConfirmDialog(
getJDialog(),
format.format(new Object[] {nameListText.toString()}),
"Incomplete Subset Definition", /*I18N*/
JOptionPane.YES_NO_CANCEL_OPTION);
if (status == JOptionPane.YES_OPTION) {
final String[] nodenames = notIncludedNames.toArray(new String[notIncludedNames.size()]);
productSubsetDef.addNodeNames(nodenames);
ok = true;
} else if (status == JOptionPane.NO_OPTION) {
ok = true;
} else if (status == JOptionPane.CANCEL_OPTION) {
ok = false;
}
}
return ok;
}
private String ComputeBinMeans(String prodname) {
StringBuilder bin_mean = new StringBuilder(prodname);
bin_mean.append("_");
bin_mean.append("sum / weights");
return bin_mean.toString();
}
