/** Used to execute this plugin tool. */
@Override
public void run() {
/*
* This transformation has been taken from:
* http://ij.ms3d.de/pdf/ihs_transforms.pdf
* in reference to Haydn, Dalke, and Henkel (1982)
* Note: 0 <= I <= 3, 0 <= H <= 3, 0 <= S <= 1
*/
amIActive = true;
String redHeader, greenHeader, blueHeader, intensityHeader, saturationHeader, hueHeader;
if (args.length <= 0) {
showFeedback("Plugin parameters have not been set.");
return;
}
redHeader = args[0];
greenHeader = args[1];
blueHeader = args[2];
intensityHeader = args[3];
hueHeader = args[4];
saturationHeader = args[5];
// check to see that the inputHeader and outputHeader are not null.
if (redHeader.isEmpty()
|| greenHeader.isEmpty()
|| blueHeader == null
|| intensityHeader.isEmpty()
|| hueHeader.isEmpty()
|| saturationHeader.isEmpty()) {
showFeedback("One or more of the input parameters have not been set properly.");
return;
}
try {
int row, col;
double redVal, greenVal, blueVal;
// double redRange, greenRange, blueRange;
// double redMin, greenMin, blueMin;
double r, g, b;
double i, s, h, m;
float progress;
WhiteboxRasterInfo red = new WhiteboxRasterInfo(redHeader);
int rows = red.getNumberRows();
int cols = red.getNumberColumns();
WhiteboxRasterInfo green = new WhiteboxRasterInfo(greenHeader);
if (green.getNumberRows() != rows || green.getNumberColumns() != cols) {
showFeedback("All input images must have the same dimensions.");
return;
}
WhiteboxRasterInfo blue = new WhiteboxRasterInfo(blueHeader);
if (blue.getNumberRows() != rows || blue.getNumberColumns() != cols) {
showFeedback("All input images must have the same dimensions.");
return;
}
double redNoData = red.getNoDataValue();
double greenNoData = green.getNoDataValue();
double blueNoData = blue.getNoDataValue();
WhiteboxRaster intensity =
new WhiteboxRaster(
intensityHeader, "rw", redHeader, WhiteboxRaster.DataType.FLOAT, redNoData);
WhiteboxRaster hue =
new WhiteboxRaster(hueHeader, "rw", redHeader, WhiteboxRaster.DataType.FLOAT, redNoData);
WhiteboxRaster saturation =
new WhiteboxRaster(
saturationHeader, "rw", redHeader, WhiteboxRaster.DataType.FLOAT, redNoData);
// redMin = red.getDisplayMinimum();
// greenMin = green.getDisplayMinimum();
// blueMin = blue.getDisplayMinimum();
//
// redRange = red.getDisplayMaximum() - redMin;
// greenRange = green.getDisplayMaximum() - greenMin;
// blueRange = blue.getDisplayMaximum() - blueMin;
double overallMin =
Math.min(
Math.min(red.getDisplayMinimum(), green.getDisplayMinimum()),
blue.getDisplayMinimum());
double overallMax =
Math.max(
Math.max(red.getDisplayMaximum(), green.getDisplayMaximum()),
blue.getDisplayMaximum());
double range = overallMax - overallMin;
double[] dataRed, dataGreen, dataBlue;
for (row = 0; row < rows; row++) {
dataRed = red.getRowValues(row);
dataGreen = green.getRowValues(row);
dataBlue = blue.getRowValues(row);
for (col = 0; col < cols; col++) {
redVal = dataRed[col];
greenVal = dataGreen[col];
blueVal = dataBlue[col];
if ((redVal != redNoData) && (greenVal != greenNoData) && (blueVal != blueNoData)) {
r = (redVal - overallMin) / range;
if (r < 0) {
r = 0;
}
if (r > 1) {
r = 1;
}
g = (greenVal - overallMin) / range;
if (g < 0) {
g = 0;
}
if (g > 1) {
g = 1;
}
b = (blueVal - overallMin) / range;
if (b < 0) {
b = 0;
}
if (b > 1) {
b = 1;
}
m = Math.min(Math.min(r, g), b);
i = r + g + b;
if (i == 3) {
h = 0;
} else if (m == b) {
h = (g - b) / (i - 3 * b);
} else if (m == r) {
h = (b - r) / (i - 3 * r) + 1;
} else { // m == g
h = (r - g) / (i - 3 * g) + 2;
}
if (h <= 1) {
s = (i - 3 * b) / i;
} else if (h <= 2) {
s = (i - 3 * r) / i;
} else { // H <= 3
s = (i - 3 * g) / i;
}
intensity.setValue(row, col, i);
hue.setValue(row, col, h);
saturation.setValue(row, col, s);
} else {
intensity.setValue(row, col, redNoData);
hue.setValue(row, col, redNoData);
saturation.setValue(row, col, redNoData);
}
}
if (cancelOp) {
cancelOperation();
return;
}
progress = (float) (100f * row / (rows - 1));
updateProgress((int) progress);
}
intensity.addMetadataEntry("Created by the " + getDescriptiveName() + " tool.");
intensity.addMetadataEntry("Created on " + new Date());
intensity.close();
hue.addMetadataEntry("Created by the " + getDescriptiveName() + " tool.");
hue.addMetadataEntry("Created on " + new Date());
hue.close();
saturation.addMetadataEntry("Created by the " + getDescriptiveName() + " tool.");
saturation.addMetadataEntry("Created on " + new Date());
saturation.close();
red.close();
green.close();
blue.close();
// returning a header file string displays the image.
returnData(intensityHeader);
returnData(hueHeader);
returnData(saturationHeader);
} catch (OutOfMemoryError oe) {
myHost.showFeedback("An out-of-memory error has occurred during operation.");
} catch (Exception e) {
myHost.showFeedback("An error has occurred during operation. See log file for details.");
myHost.logException("Error in " + getDescriptiveName(), e);
} finally {
updateProgress("Progress: ", 0);
// tells the main application that this process is completed.
amIActive = false;
myHost.pluginComplete();
}
}
