-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
c76937f
commit 2b4fb1d
Showing
8 changed files
with
132 additions
and
149 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file was deleted.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
1 change: 1 addition & 0 deletions
1
...in/java/org/restopt/RasterDataLoader.java → ...st/java/org/restopt/RasterDataLoader.java
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,120 @@ | ||
| package org.restopt; | ||
|
|
||
| import org.chocosolver.util.objects.setDataStructures.ISet; | ||
| import org.chocosolver.util.objects.setDataStructures.SetFactory; | ||
| import org.chocosolver.util.tools.ArrayUtils; | ||
| import org.geotools.coverage.grid.GridCoverage2D; | ||
| import org.geotools.coverage.grid.GridCoverageFactory; | ||
| import org.geotools.gce.geotiff.GeoTiffReader; | ||
| import org.geotools.gce.geotiff.GeoTiffWriter; | ||
| import org.restopt.exception.RestoptException; | ||
| import org.restopt.grid.regular.square.PartialRegularGroupedGrid; | ||
|
|
||
| import java.awt.*; | ||
| import java.awt.image.*; | ||
| import java.io.BufferedWriter; | ||
| import java.io.File; | ||
| import java.io.FileWriter; | ||
| import java.io.IOException; | ||
|
|
||
|
|
||
| public class SolutionExporter { | ||
|
|
||
| public String csvDest; | ||
| public String rastDest; | ||
| public String template; | ||
| public int[] sites; | ||
| public int[] completeData; | ||
| public RestoptProblem restoptProblem; | ||
| public RestoptSolution solution; | ||
|
|
||
| public SolutionExporter(RestoptSolution solution, String csvDest, String rastDest, double noDataValue) { | ||
| this.solution = solution; | ||
| this.restoptProblem = solution.getProblem(); | ||
| if (restoptProblem.getData() instanceof RasterDataLoader) { | ||
| RasterDataLoader dataLoader = (RasterDataLoader) restoptProblem.getData(); | ||
| PartialRegularGroupedGrid grid = restoptProblem.getGrid(); | ||
| sites = new int[grid.getNbUngroupedCells()]; | ||
| ISet set = SetFactory.makeConstantSet(solution.getRestorationPlanningUnits()); | ||
| for (int i = 0; i < grid.getNbUngroupedCells(); i++) { | ||
| if (grid.getGroupIndexFromPartialIndex(i) < grid.getNbGroups()) { | ||
| sites[i] = 2; | ||
| } else if (set.contains(grid.getGroupIndexFromPartialIndex(i))) { | ||
| sites[i] = 3; | ||
| } else { | ||
| sites[i] = 1; | ||
| } | ||
| } | ||
| this.template = dataLoader.getHabitatRasterPath(); | ||
| this.csvDest = csvDest; | ||
| this.rastDest = rastDest; | ||
| completeData = new int[grid.getNbRows() * grid.getNbCols()]; | ||
| for (int i = 0; i < completeData.length; i++) { | ||
| if (grid.getDiscardSet().contains(i)) { | ||
| if (restoptProblem.data.getHabitatData()[i] == 0) { | ||
| completeData[i] = 0; | ||
| } else { | ||
| completeData[i] = (int) noDataValue; | ||
| } | ||
| } else { | ||
| completeData[i] = sites[grid.getPartialIndex(i)]; | ||
| } | ||
| } | ||
| } | ||
| } | ||
|
|
||
| public void exportCharacteristics() throws IOException { | ||
| String[][] orderedCharacteristics = new String[2][]; | ||
| String[] allKeys = ArrayUtils.append(solution.KEYS, solution.getObjective().getAdditionalKeys()); | ||
| orderedCharacteristics[0] = allKeys; | ||
| orderedCharacteristics[1] = new String[allKeys.length]; | ||
| for (int i = 0; i < allKeys.length; i++) { | ||
| orderedCharacteristics[1][i] = solution.getCharacteristics().get(allKeys[i]); | ||
| } | ||
| BufferedWriter br = new BufferedWriter(new FileWriter(csvDest)); | ||
| StringBuilder sb = new StringBuilder(); | ||
| for (String[] line : orderedCharacteristics) { | ||
| int i = 0; | ||
| for (String s : line) { | ||
| i++; | ||
| sb.append(s); | ||
| if (i < line.length) { | ||
| sb.append(","); | ||
| } | ||
| } | ||
| sb.append("\n"); | ||
| } | ||
| br.write(sb.toString()); | ||
| br.close(); | ||
| } | ||
|
|
||
| public void generateRaster() throws IOException { | ||
| File file = new File(template); | ||
| GeoTiffReader reader = new GeoTiffReader(file); | ||
| GridCoverage2D grid = reader.read(null); | ||
| int height = grid.getRenderedImage().getHeight(); | ||
| int width = grid.getRenderedImage().getWidth(); | ||
|
|
||
| GeoTiffWriter writer = new GeoTiffWriter(new File(rastDest)); | ||
|
|
||
| DataBuffer buff = grid.getRenderedImage().getData().getDataBuffer(); | ||
| SampleModel sm = new BandedSampleModel(DataBuffer.TYPE_INT, width, height, 1); | ||
| WritableRaster rast = Raster.createWritableRaster(sm, buff, new Point(0, 0)); | ||
| rast.setPixels(0, 0, width, height, completeData); | ||
|
|
||
| GridCoverageFactory f = new GridCoverageFactory(); | ||
| GridCoverage2D destCov = f.create("rast", rast, grid.getEnvelope()); | ||
|
|
||
| writer.write(destCov, null); | ||
| } | ||
|
|
||
| public void export(boolean verbose) throws IOException, RestoptException { | ||
| exportCharacteristics(); | ||
| generateRaster(); | ||
| if (verbose) { | ||
| this.solution.printSolutionInfos(); | ||
| System.out.println("\nRaster exported at " + rastDest + ".tif"); | ||
| System.out.println("Solution characteristics exported at " + csvDest + ".csv\n"); | ||
| } | ||
| } | ||
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
File renamed without changes.