/**
* sustituye el esquema de JUMP por un esquema automático de Geopista. Para una integración
* rápida. SOLO PARA PRUEBAS
*
* @param feature
*/
public static GeopistaFeature vampiriceSchema(Feature feature) {
FeatureSchema fSchema = feature.getSchema();
GeopistaSchema gSchema = new GeopistaSchema();
// crea el atributo
Table defTable = new Table("DatosCapa", "Tabla de la BBDD para esta capa.");
for (int i = 0; i < fSchema.getAttributeCount(); i++) {
// define Dominio, Columna y añade al esquema
Domain defDomain;
if (fSchema.getAttributeType(i) == AttributeType.DOUBLE
|| fSchema.getAttributeType(i) == AttributeType.FLOAT) {
defDomain = new NumberDomain("?[-INF:INF]", "Default Double Number.");
} else if (fSchema.getAttributeType(i) == AttributeType.INTEGER
|| fSchema.getAttributeType(i) == AttributeType.LONG) {
// Definir patron para excluir los decimales.
defDomain = new NumberDomain("?[-INF:INF]", "Default Integer Number");
} else {
defDomain = new StringDomain("?[.*]", "Dominio por defecto.");
}
Column col =
new Column(fSchema.getAttributeName(i), "Columna automática.", defTable, defDomain);
gSchema.addAttribute(
fSchema.getAttributeName(i), fSchema.getAttributeType(i), col, READ_WRITE);
}
GeopistaFeature gFeature = new GeopistaFeature(gSchema);
for (int i = 0; i < fSchema.getAttributeCount(); i++) {
gFeature.setAttribute(fSchema.getAttributeName(i), feature.getAttribute(i));
}
return gFeature;
}
protected void copySelectedAttributes() {
this.attributes = new FeatureSchema();
for (int i = 0; i < tempAttributes.getAttributeCount(); i++) {
for (int i2 = 0; i2 < this.onlyTypes.length; i2++) {
if (this.tempAttributes.getAttributeType(i) == this.onlyTypes[i2]) {
this.attributes.addAttribute(
tempAttributes.getAttributeName(i), tempAttributes.getAttributeType(i));
}
}
}
}
/**
* Write a dbf file with the information from the featureCollection.
*
* @param featureCollection column data from collection
* @param fname name of the dbf file to write to
*/
void writeDbf(FeatureCollection featureCollection, String fname) throws Exception {
DbfFileWriter dbf;
FeatureSchema fs;
int t;
int f;
int u;
int num;
fs = featureCollection.getFeatureSchema();
// -1 because one of the columns is geometry
DbfFieldDef[] fields = new DbfFieldDef[fs.getAttributeCount() - 1];
// dbf column type and size
f = 0;
for (t = 0; t < fs.getAttributeCount(); t++) {
AttributeType columnType = fs.getAttributeType(t);
String columnName = fs.getAttributeName(t);
if (columnType == AttributeType.INTEGER) {
fields[f] = new DbfFieldDef(columnName, 'N', 16, 0);
f++;
} else if (columnType == AttributeType.DOUBLE) {
fields[f] = new DbfFieldDef(columnName, 'N', 33, 16);
f++;
} else if (columnType == AttributeType.STRING) {
int maxlength = findMaxStringLength(featureCollection, t);
if (maxlength > 255) {
throw new Exception(
"ShapefileWriter does not support strings longer than 255 characters");
}
fields[f] = new DbfFieldDef(columnName, 'C', maxlength, 0);
f++;
} else if (columnType == AttributeType.DATE) {
fields[f] = new DbfFieldDef(columnName, 'D', 8, 0);
f++;
} else if (columnType == AttributeType.GEOMETRY) {
// do nothing - the .shp file handles this
} else {
throw new Exception("Shapewriter: unsupported AttributeType found in featurecollection.");
}
}
// write header
dbf = new DbfFileWriter(fname);
dbf.writeHeader(fields, featureCollection.size());
// write rows
num = featureCollection.size();
List features = featureCollection.getFeatures();
for (t = 0; t < num; t++) {
// System.out.println("dbf: record "+t);
Feature feature = (Feature) features.get(t);
Vector DBFrow = new Vector();
// make data for each column in this feature (row)
for (u = 0; u < fs.getAttributeCount(); u++) {
AttributeType columnType = fs.getAttributeType(u);
if (columnType == AttributeType.INTEGER) {
Object a = feature.getAttribute(u);
if (a == null) {
DBFrow.add(new Integer(0));
} else {
DBFrow.add((Integer) a);
}
} else if (columnType == AttributeType.DOUBLE) {
Object a = feature.getAttribute(u);
if (a == null) {
DBFrow.add(new Double(0.0));
} else {
DBFrow.add((Double) a);
}
} else if (columnType == AttributeType.DATE) {
Object a = feature.getAttribute(u);
if (a == null) {
DBFrow.add("");
} else {
DBFrow.add(DbfFile.DATE_PARSER.format((Date) a));
}
} else if (columnType == AttributeType.STRING) {
Object a = feature.getAttribute(u);
if (a == null) {
DBFrow.add(new String(""));
} else {
// MD 16 jan 03 - added some defensive programming
if (a instanceof String) {
DBFrow.add(a);
} else {
DBFrow.add(a.toString());
}
}
}
}
dbf.writeRecord(DBFrow);
}
dbf.close();
}
private FeatureDataset classifyAndCreatePlot(TaskMonitor monitor, final PlugInContext context)
throws Exception {
monitor.report(this.sCalculateBreaks);
// =============== get DATA and prepare ==============/
FeatureSchema fs = this.fc.getFeatureSchema();
AttributeType type = null;
if ((fs.getAttributeType(this.selAttribute) == AttributeType.DOUBLE)
|| (fs.getAttributeType(this.selAttribute) == AttributeType.INTEGER)) {
// -- move on
type = fs.getAttributeType(this.selAttribute);
} else {
// System.out.println("ClassifyAttributesPlugIn: wrong datatype of chosen attribute");
context.getWorkbenchFrame().warnUser(sWrongDataType);
return null;
}
int size = getFeatureCollectionSize(this.fc, this.selAttribute, this.nullAsZero);
if (size < 3) {
return null;
}
this.ranges = Math.min(this.ranges, size);
double[] data = new double[size];
double[][] plotdata = new double[2][size]; // for drawing 1-D scatter plot
int[] fID = new int[size];
int i = 0;
for (Iterator iter = fc.iterator(); iter.hasNext(); ) {
Feature f = (Feature) iter.next();
if (f.getAttribute(this.selAttribute) == null && !nullAsZero) continue;
fID[i] = f.getID();
plotdata[1][i] = 1;
Object val = f.getAttribute(this.selAttribute);
if (type == AttributeType.DOUBLE) {
if (val == null) data[i] = 0.0;
else data[i] = ((Double) val).doubleValue();
} else if (type == AttributeType.INTEGER) {
if (val == null) data[i] = 0;
else data[i] = ((Integer) val).intValue();
}
plotdata[0][i] = data[i];
i++;
}
/*
//-- some testdata
double[][] plotdata2 = new double[2][8];
double[] data2 = { -2, 4, 6, 5, 0, 10, 7, 1 };
double[] axis2 = { 1, 1, 1, 1, 1, 1, 1, 1 };
plotdata2[0] = data2;
plotdata2[1] = axis2;
*/
if (monitor.isCancelRequested()) {
return null;
}
// =============== find breaks according to chosen method ==============/
double[] limits = null;
if (this.useKmeans == false) {
if (this.selClassifier == Classifier1D.EQUAL_NUMBER) {
limits = Classifier1D.classifyEqualNumber(data, this.ranges);
} else if (this.selClassifier == Classifier1D.EQUAL_RANGE) {
limits = Classifier1D.classifyEqualRange(data, this.ranges);
} else if (this.selClassifier == Classifier1D.MEAN_STDEV) {
limits = Classifier1D.classifyMeanStandardDeviation(data, this.ranges);
} else if (this.selClassifier == Classifier1D.MAX_BREAKS) {
limits = Classifier1D.classifyMaxBreaks(data, this.ranges);
} else if (this.selClassifier == Classifier1D.JENKS_BREAKS) {
limits = Classifier1D.classifyNaturalBreaks(data, this.ranges);
}
} else {
if (this.selClassifier == Classifier1D.EQUAL_NUMBER) {
limits = Classifier1D.classifyKMeansOnExistingBreaks(data, this.ranges, 3);
} else if (this.selClassifier == Classifier1D.EQUAL_RANGE) {
limits = Classifier1D.classifyKMeansOnExistingBreaks(data, this.ranges, 2);
} else if (this.selClassifier == Classifier1D.MEAN_STDEV) {
limits = Classifier1D.classifyKMeansOnExistingBreaks(data, this.ranges, 4);
} else if (this.selClassifier == Classifier1D.MAX_BREAKS) {
limits = Classifier1D.classifyKMeansOnExistingBreaks(data, this.ranges, 1);
} else if (this.selClassifier == Classifier1D.JENKS_BREAKS) {
limits = Classifier1D.classifyKMeansOnExistingBreaks(data, this.ranges, 5);
}
}
if (monitor.isCancelRequested()) {
return null;
}
monitor.report(this.sDisplayBreaks);
// =============== plot data and class breaks ==============/
// -- do display here - in case we later want to allow interactive editing of the limits
// -- reformat limits
double[][] limits2show = new double[2][limits.length];
// -- due to bug in jmathplot add limits twice if only three classes = 2breaks are sought
if (limits.length == 2) {
limits2show = new double[2][limits.length * 2];
}
for (int j = 0; j < limits.length; j++) {
limits2show[0][j] = limits[j]; // x-axis
limits2show[1][j] =
Math.floor(
i / (4.0 * this.ranges)); // y-axis, estimate height of "bar" from number of items
// limits2show[1][j]= 1;
// -- due to bug in jmathplot add limits twice if only three classes are sought
if (limits.length == 2) {
limits2show[0][limits.length + j] = limits[j];
limits2show[1][limits.length + j] = Math.floor(i / (4.0 * this.ranges));
}
}
// =============== plot data and class breaks ==============/
// -- create plots
/*final Plot2DPanelOJ*/ plot = new Plot2DPanelOJ();
plot.addHistogramPlotOJ(
this.selAttribute, data, this.ranges * 3, context, selLayer, this.selAttribute);
plot.addScatterPlotOJ(this.sDatapoints, plotdata, fID, context, this.selLayer);
plot.addBarPlot(this.sClassbreaks, limits2show);
plot.plotToolBar.setVisible(true);
plot.setAxisLabel(0, this.selAttribute);
plot.setAxisLabel(1, this.sCount);
plot.addLegend("SOUTH");
// [mmichaud 2012-04-09] Moved in run method after the addLayer method
// to avoid the problem of the focus change
// JInternalFrame frame = new JInternalFrame(this.sHistogram);
// frame.setLayout(new BorderLayout());
// frame.add(plot, BorderLayout.CENTER);
// frame.setClosable(true);
// frame.setResizable(true);
// frame.setMaximizable(true);
// frame.setSize(450, 450);
// frame.setVisible(true);
// context.getWorkbenchFrame().addInternalFrame(frame);
// =============== classify data ==============/
if (monitor.isCancelRequested()) {
return null;
}
monitor.report(this.sClassifying);
int[] classes = Classifier1D.classifyData(data, limits);
// double[] classes = org.math.array.StatisticSample.one(data.length);
// context.getWorkbenchFrame().warnUser("classification not yet implemented");
// =============== add field ==============/
if (monitor.isCancelRequested()) {
return null;
}
monitor.report(sAddingField);
FeatureDataset fd = null;
ArrayList outData = new ArrayList();
FeatureSchema targetFSnew = null;
int count = 0;
Iterator iterp = fc.iterator();
String attname = this.selAttribute + "_" + this.selClassifier;
while (iterp.hasNext()) {
// count=count+1;
// if(monitor != null){
// monitor.report("item: " + count + " of " + size);
// }
Feature p = (Feature) iterp.next();
Object val = p.getAttribute(this.selAttribute);
if (val == null && !this.nullAsZero) continue;
else count++;
if (count == 1) {
FeatureSchema targetFs = p.getSchema();
targetFSnew = FeatureSchemaTools.copyFeatureSchema(targetFs);
if (targetFSnew.hasAttribute(attname)) {
// attribute will be overwriten
} else {
// add attribute
targetFSnew.addAttribute(attname, AttributeType.INTEGER);
}
}
// -- evaluate value for every polygon
Feature fcopy = FeatureSchemaTools.copyFeature(p, targetFSnew);
// fcopy.setAttribute(this.selClassifier, new Integer(classes[count-1]));
fcopy.setAttribute(attname, new Integer(classes[count - 1]));
outData.add(fcopy);
}
fd = new FeatureDataset(targetFSnew);
fd.addAll(outData);
return fd;
}