/*
* GeoTools - The Open Source Java GIS Toolkit
* http://geotools.org
*
* (C) 2006-2010, Open Source Geospatial Foundation (OSGeo)
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation;
* version 2.1 of the License.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*/
package org.geotools.gce.grassraster;
import java.awt.Color;
import java.awt.Transparency;
import java.awt.color.ColorSpace;
import java.awt.image.BufferedImage;
import java.awt.image.ColorModel;
import java.awt.image.ComponentColorModel;
import java.awt.image.DataBuffer;
import java.awt.image.DataBufferByte;
import java.awt.image.Raster;
import java.awt.image.WritableRaster;
import java.io.BufferedWriter;
import java.io.File;
import java.io.FileWriter;
import java.io.IOException;
import java.util.List;
import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.GridEnvelope2D;
import org.geotools.coverage.grid.InvalidGridGeometryException;
import org.geotools.geometry.Envelope2D;
import org.geotools.geometry.jts.JTS;
import org.geotools.referencing.CRS;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;
import org.opengis.referencing.operation.TransformException;
import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.Envelope;
/**
* <p>
* A facade of often used methods by the JGrass engine
* </p>
*
* @author Andrea Antonello - www.hydrologis.com
* @since 1.1.0
*
* @source $URL: http://svn.osgeo.org/geotools/branches/2.7.x/build/maven/javadoc/../../../modules/unsupported/grassraster/src/main/java/org/geotools/gce/grassraster/JGrassUtilities.java $
*/
public class JGrassUtilities {
public static final String[] COLORS = {"snow", "GhostWhite", "WhiteSmoke", "gainsboro",
"FloralWhite", "OldLace", "linen", "antiquewhite", "papaya whip", "PapayaWhip",
"blanched almond", "BlanchedAlmond", "bisque", "peach puff", "PeachPuff",
"navajo white", "NavajoWhite", "moccasin", "cornsilk", "ivory", "lemon chiffon",
"LemonChiffon", "seashell", "honeydew", "mint cream", "MintCream", "azure",
"alice blue", "AliceBlue", "lavender", "lavender blush", "LavenderBlush", "misty rose",
"MistyRose", "white", "black", "dark slate gray", "DarkSlateGray", "dark slate grey",
"DarkSlateGrey", "dim gray", "DimGray", "dim grey", "DimGrey", "slate gray",
"SlateGray", "slate grey", "SlateGrey", "light slate gray", "LightSlateGray",
"light slate grey", "LightSlateGrey", "gray", "grey", "light grey", "LightGrey",
"light gray", "LightGray", "midnight blue", "MidnightBlue", "navy", "navy blue",
"NavyBlue", "cornflower blue", "CornflowerBlue", "dark slate blue", "DarkSlateBlue",
"slate blue", "SlateBlue", "medium slate blue", "MediumSlateBlue", "light slate blue",
"LightSlateBlue", "medium blue", "MediumBlue", "royal blue", "RoyalBlue", "blue",
"dodger blue", "DodgerBlue", "deep sky blue", "DeepSkyBlue", "sky blue", "SkyBlue",
"light sky blue", "LightSkyBlue", "steel blue", "SteelBlue", "light steel blue",
"LightSteelBlue", "light blue", "LightBlue", "powder blue", "PowderBlue",
"pale turquoise", "PaleTurquoise", "dark turquoise", "DarkTurquoise",
"medium turquoise", "MediumTurquoise", "turquoise", "cyan", "light cyan", "LightCyan",
"cadet blue", "CadetBlue", "medium aquamarine", "MediumAquamarine", "aquamarine",
"dark green", "DarkGreen", "dark olive green", "DarkOliveGreen", "dark sea green",
"DarkSeaGreen", "sea green", "SeaGreen", "medium sea green", "MediumSeaGreen",
"light sea green", "LightSeaGreen", "pale green", "PaleGreen", "spring green",
"SpringGreen", "lawn green", "LawnGreen", "green", "chartreuse", "medium spring green",
"MediumSpringGreen", "green yellow", "GreenYellow", "lime green", "LimeGreen",
"yellow green", "YellowGreen", "forest green", "ForestGreen", "olive drab",
"OliveDrab", "dark khaki", "DarkKhaki", "khaki", "pale goldenrod", "PaleGoldenrod",
"light goldenrod yellow", "LightGoldenrodYellow", "light yellow", "LightYellow",
"yellow", "gold", "light goldenrod", "LightGoldenrod", "goldenrod", "dark goldenrod",
"DarkGoldenrod", "rosy brown", "RosyBrown", "indian red", "IndianRed", "saddle brown",
"SaddleBrown", "sienna", "peru", "burlywood", "beige", "wheat", "sandy brown",
"SandyBrown", "tan", "chocolate", "firebrick", "brown", "dark salmon", "DarkSalmon",
"salmon", "light salmon", "LightSalmon", "orange", "dark orange", "DarkOrange",
"coral", "light coral", "LightCoral", "tomato", "orange red", "OrangeRed", "red",
"hot pink", "HotPink", "deep pink", "DeepPink", "pink", "light pink", "LightPink",
"pale violet red", "PaleVioletRed", "maroon", "medium violet red", "MediumVioletRed",
"violet red", "VioletRed", "magenta", "violet", "plum", "orchid", "medium orchid",
"MediumOrchid", "dark orchid", "DarkOrchid", "dark violet", "DarkViolet",
"blue violet", "BlueViolet", "purple", "medium purple", "MediumPurple", "thistle",
"snow1", "snow2", "snow3", "snow4", "seashell1", "seashell2", "seashell3", "seashell4",
"AntiqueWhite1", "AntiqueWhite2", "AntiqueWhite3", "AntiqueWhite4", "bisque1",
"bisque2", "bisque3", "bisque4", "PeachPuff1", "PeachPuff2", "PeachPuff3",
"PeachPuff4", "NavajoWhite1", "NavajoWhite2", "NavajoWhite3", "NavajoWhite4",
"LemonChiffon1", "LemonChiffon2", "LemonChiffon3", "LemonChiffon4", "cornsilk1",
"cornsilk2", "cornsilk3", "cornsilk4", "ivory1", "ivory2", "ivory3", "ivory4",
"honeydew1", "honeydew2", "honeydew3", "honeydew4", "LavenderBlush1", "LavenderBlush2",
"LavenderBlush3", "LavenderBlush4", "MistyRose1", "MistyRose2", "MistyRose3",
"MistyRose4", "azure1", "azure2", "azure3", "azure4", "SlateBlue1", "SlateBlue2",
"SlateBlue3", "SlateBlue4", "RoyalBlue1", "RoyalBlue2", "RoyalBlue3", "RoyalBlue4",
"blue1", "blue2", "blue3", "blue4", "DodgerBlue1", "DodgerBlue2", "DodgerBlue3",
"DodgerBlue4", "SteelBlue1", "SteelBlue2", "SteelBlue3", "SteelBlue4", "DeepSkyBlue1",
"DeepSkyBlue2", "DeepSkyBlue3", "DeepSkyBlue4", "SkyBlue1", "SkyBlue2", "SkyBlue3",
"SkyBlue4", "LightSkyBlue1", "LightSkyBlue2", "LightSkyBlue3", "LightSkyBlue4",
"SlateGray1", "SlateGray2", "SlateGray3", "SlateGray4", "LightSteelBlue1",
"LightSteelBlue2", "LightSteelBlue3", "LightSteelBlue4", "LightBlue1", "LightBlue2",
"LightBlue3", "LightBlue4", "LightCyan1", "LightCyan2", "LightCyan3", "LightCyan4",
"PaleTurquoise1", "PaleTurquoise2", "PaleTurquoise3", "PaleTurquoise4", "CadetBlue1",
"CadetBlue2", "CadetBlue3", "CadetBlue4", "turquoise1", "turquoise2", "turquoise3",
"turquoise4", "cyan1", "cyan2", "cyan3", "cyan4", "DarkSlateGray1", "DarkSlateGray2",
"DarkSlateGray3", "DarkSlateGray4", "aquamarine1", "aquamarine2", "aquamarine3",
"aquamarine4", "DarkSeaGreen1", "DarkSeaGreen2", "DarkSeaGreen3", "DarkSeaGreen4",
"SeaGreen1", "SeaGreen2", "SeaGreen3", "SeaGreen4", "PaleGreen1", "PaleGreen2",
"PaleGreen3", "PaleGreen4", "SpringGreen1", "SpringGreen2", "SpringGreen3",
"SpringGreen4", "green1", "green2", "green3", "green4", "chartreuse1", "chartreuse2",
"chartreuse3", "chartreuse4", "OliveDrab1", "OliveDrab2", "OliveDrab3", "OliveDrab4",
"DarkOliveGreen1", "DarkOliveGreen2", "DarkOliveGreen3", "DarkOliveGreen4", "khaki1",
"khaki2", "khaki3", "khaki4", "LightGoldenrod1", "LightGoldenrod2", "LightGoldenrod3",
"LightGoldenrod4", "LightYellow1", "LightYellow2", "LightYellow3", "LightYellow4",
"yellow1", "yellow2", "yellow3", "yellow4", "gold1", "gold2", "gold3", "gold4",
"goldenrod1", "goldenrod2", "goldenrod3", "goldenrod4", "DarkGoldenrod1",
"DarkGoldenrod2", "DarkGoldenrod3", "DarkGoldenrod4", "RosyBrown1", "RosyBrown2",
"RosyBrown3", "RosyBrown4", "IndianRed1", "IndianRed2", "IndianRed3", "IndianRed4",
"sienna1", "sienna2", "sienna3", "sienna4", "burlywood1", "burlywood2", "burlywood3",
"burlywood4", "wheat1", "wheat2", "wheat3", "wheat4", "tan1", "tan2", "tan3", "tan4",
"chocolate1", "chocolate2", "chocolate3", "chocolate4", "firebrick1", "firebrick2",
"firebrick3", "firebrick4", "brown1", "brown2", "brown3", "brown4", "salmon1",
"salmon2", "salmon3", "salmon4", "LightSalmon1", "LightSalmon2", "LightSalmon3",
"LightSalmon4", "orange1", "orange2", "orange3", "orange4", "DarkOrange1",
"DarkOrange2", "DarkOrange3", "DarkOrange4", "coral1", "coral2", "coral3", "coral4",
"tomato1", "tomato2", "tomato3", "tomato4", "OrangeRed1", "OrangeRed2", "OrangeRed3",
"OrangeRed4", "red1", "red2", "red3", "red4", "DeepPink1", "DeepPink2", "DeepPink3",
"DeepPink4", "HotPink1", "HotPink2", "HotPink3", "HotPink4", "pink1", "pink2", "pink3",
"pink4", "LightPink1", "LightPink2", "LightPink3", "LightPink4", "PaleVioletRed1",
"PaleVioletRed2", "PaleVioletRed3", "PaleVioletRed4", "maroon1", "maroon2", "maroon3",
"maroon4", "VioletRed1", "VioletRed2", "VioletRed3", "VioletRed4", "magenta1",
"magenta2", "magenta3", "magenta4", "orchid1", "orchid2", "orchid3", "orchid4",
"plum1", "plum2", "plum3", "plum4", "MediumOrchid1", "MediumOrchid2", "MediumOrchid3",
"MediumOrchid4", "DarkOrchid1", "DarkOrchid2", "DarkOrchid3", "DarkOrchid4", "purple1",
"purple2", "purple3", "purple4", "MediumPurple1", "MediumPurple2", "MediumPurple3",
"MediumPurple4", "thistle1", "thistle2", "thistle3", "thistle4", "gray0", "grey0",
"gray1", "grey1", "gray2", "grey2", "gray3", "grey3", "gray4", "grey4", "gray5",
"grey5", "gray6", "grey6", "gray7", "grey7", "gray8", "grey8", "gray9", "grey9",
"gray10", "grey10", "gray11", "grey11", "gray12", "grey12", "gray13", "grey13",
"gray14", "grey14", "gray15", "grey15", "gray16", "grey16", "gray17", "grey17",
"gray18", "grey18", "gray19", "grey19", "gray20", "grey20", "gray21", "grey21",
"gray22", "grey22", "gray23", "grey23", "gray24", "grey24", "gray25", "grey25",
"gray26", "grey26", "gray27", "grey27", "gray28", "grey28", "gray29", "grey29",
"gray30", "grey30", "gray31", "grey31", "gray32", "grey32", "gray33", "grey33",
"gray34", "grey34", "gray35", "grey35", "gray36", "grey36", "gray37", "grey37",
"gray38", "grey38", "gray39", "grey39", "gray40", "grey40", "gray41", "grey41",
"gray42", "grey42", "gray43", "grey43", "gray44", "grey44", "gray45", "grey45",
"gray46", "grey46", "gray47", "grey47", "gray48", "grey48", "gray49", "grey49",
"gray50", "grey50", "gray51", "grey51", "gray52", "grey52", "gray53", "grey53",
"gray54", "grey54", "gray55", "grey55", "gray56", "grey56", "gray57", "grey57",
"gray58", "grey58", "gray59", "grey59", "gray60", "grey60", "gray61", "grey61",
"gray62", "grey62", "gray63", "grey63", "gray64", "grey64", "gray65", "grey65",
"gray66", "grey66", "gray67", "grey67", "gray68", "grey68", "gray69", "grey69",
"gray70", "grey70", "gray71", "grey71", "gray72", "grey72", "gray73", "grey73",
"gray74", "grey74", "gray75", "grey75", "gray76", "grey76", "gray77", "grey77",
"gray78", "grey78", "gray79", "grey79", "gray80", "grey80", "gray81", "grey81",
"gray82", "grey82", "gray83", "grey83", "gray84", "grey84", "gray85", "grey85",
"gray86", "grey86", "gray87", "grey87", "gray88", "grey88", "gray89", "grey89",
"gray90", "grey90", "gray91", "grey91", "gray92", "grey92", "gray93", "grey93",
"gray94", "grey94", "gray95", "grey95", "gray96", "grey96", "gray97", "grey97",
"gray98", "grey98", "gray99", "grey99", "gray100", "grey100", "dark grey", "DarkGrey",
"dark gray", "DarkGray", "dark blue", "DarkBlue", "dark cyan", "DarkCyan",
"dark magenta", "DarkMagenta", "dark red", "DarkRed", "light"};
private static final int[][] COLORS_RGB = {{255, 250, 250}, {248, 248, 255}, {245, 245, 245},
{220, 220, 220}, {255, 250, 240}, {253, 245, 230}, {250, 240, 230}, {250, 235, 215},
{255, 239, 213}, {255, 239, 213}, {255, 235, 205}, {255, 235, 205}, {255, 228, 196},
{255, 218, 185}, {255, 218, 185}, {255, 222, 173}, {255, 222, 173}, {255, 228, 181},
{255, 248, 220}, {255, 255, 240}, {255, 250, 205}, {255, 250, 205}, {255, 245, 238},
{240, 255, 240}, {245, 255, 250}, {245, 255, 250}, {240, 255, 255}, {240, 248, 255},
{240, 248, 255}, {230, 230, 250}, {255, 240, 245}, {255, 240, 245}, {255, 228, 225},
{255, 228, 225}, {255, 255, 255}, {0, 0, 0}, {47, 79, 79}, {47, 79, 79}, {47, 79, 79},
{47, 79, 79}, {105, 105, 105}, {105, 105, 105}, {105, 105, 105}, {105, 105, 105},
{112, 128, 144}, {112, 128, 144}, {112, 128, 144}, {112, 128, 144}, {119, 136, 153},
{119, 136, 153}, {119, 136, 153}, {119, 136, 153}, {190, 190, 190}, {190, 190, 190},
{211, 211, 211}, {211, 211, 211}, {211, 211, 211}, {211, 211, 211}, {25, 25, 112},
{25, 25, 112}, {0, 0, 128}, {0, 0, 128}, {0, 0, 128}, {100, 149, 237}, {100, 149, 237},
{72, 61, 139}, {72, 61, 139}, {106, 90, 205}, {106, 90, 205}, {123, 104, 238},
{123, 104, 238}, {132, 112, 255}, {132, 112, 255}, {0, 0, 205}, {0, 0, 205},
{65, 105, 225}, {65, 105, 225}, {0, 0, 255}, {30, 144, 255}, {30, 144, 255},
{0, 191, 255}, {0, 191, 255}, {135, 206, 235}, {135, 206, 235}, {135, 206, 250},
{135, 206, 250}, {70, 130, 180}, {70, 130, 180}, {176, 196, 222}, {176, 196, 222},
{173, 216, 230}, {173, 216, 230}, {176, 224, 230}, {176, 224, 230}, {175, 238, 238},
{175, 238, 238}, {0, 206, 209}, {0, 206, 209}, {72, 209, 204}, {72, 209, 204},
{64, 224, 208}, {0, 255, 255}, {224, 255, 255}, {224, 255, 255}, {95, 158, 160},
{95, 158, 160}, {102, 205, 170}, {102, 205, 170}, {127, 255, 212}, {0, 100, 0},
{0, 100, 0}, {85, 107, 47}, {85, 107, 47}, {143, 188, 143}, {143, 188, 143},
{46, 139, 87}, {46, 139, 87}, {60, 179, 113}, {60, 179, 113}, {32, 178, 170},
{32, 178, 170}, {152, 251, 152}, {152, 251, 152}, {0, 255, 127}, {0, 255, 127},
{124, 252, 0}, {124, 252, 0}, {0, 255, 0}, {127, 255, 0}, {0, 250, 154}, {0, 250, 154},
{173, 255, 47}, {173, 255, 47}, {50, 205, 50}, {50, 205, 50}, {154, 205, 50},
{154, 205, 50}, {34, 139, 34}, {34, 139, 34}, {107, 142, 35}, {107, 142, 35},
{189, 183, 107}, {189, 183, 107}, {240, 230, 140}, {238, 232, 170}, {238, 232, 170},
{250, 250, 210}, {250, 250, 210}, {255, 255, 224}, {255, 255, 224}, {255, 255, 0},
{255, 215, 0}, {238, 221, 130}, {238, 221, 130}, {218, 165, 32}, {184, 134, 11},
{184, 134, 11}, {188, 143, 143}, {188, 143, 143}, {205, 92, 92}, {205, 92, 92},
{139, 69, 19}, {139, 69, 19}, {160, 82, 45}, {205, 133, 63}, {222, 184, 135},
{245, 245, 220}, {245, 222, 179}, {244, 164, 96}, {244, 164, 96}, {210, 180, 140},
{210, 105, 30}, {178, 34, 34}, {165, 42, 42}, {233, 150, 122}, {233, 150, 122},
{250, 128, 114}, {255, 160, 122}, {255, 160, 122}, {255, 165, 0}, {255, 140, 0},
{255, 140, 0}, {255, 127, 80}, {240, 128, 128}, {240, 128, 128}, {255, 99, 71},
{255, 69, 0}, {255, 69, 0}, {255, 0, 0}, {255, 105, 180}, {255, 105, 180},
{255, 20, 147}, {255, 20, 147}, {255, 192, 203}, {255, 182, 193}, {255, 182, 193},
{219, 112, 147}, {219, 112, 147}, {176, 48, 96}, {199, 21, 133}, {199, 21, 133},
{208, 32, 144}, {208, 32, 144}, {255, 0, 255}, {238, 130, 238}, {221, 160, 221},
{218, 112, 214}, {186, 85, 211}, {186, 85, 211}, {153, 50, 204}, {153, 50, 204},
{148, 0, 211}, {148, 0, 211}, {138, 43, 226}, {138, 43, 226}, {160, 32, 240},
{147, 112, 219}, {147, 112, 219}, {216, 191, 216}, {255, 250, 250}, {238, 233, 233},
{205, 201, 201}, {139, 137, 137}, {255, 245, 238}, {238, 229, 222}, {205, 197, 191},
{139, 134, 130}, {255, 239, 219}, {238, 223, 204}, {205, 192, 176}, {139, 131, 120},
{255, 228, 196}, {238, 213, 183}, {205, 183, 158}, {139, 125, 107}, {255, 218, 185},
{238, 203, 173}, {205, 175, 149}, {139, 119, 101}, {255, 222, 173}, {238, 207, 161},
{205, 179, 139}, {139, 121, 94}, {255, 250, 205}, {238, 233, 191}, {205, 201, 165},
{139, 137, 112}, {255, 248, 220}, {238, 232, 205}, {205, 200, 177}, {139, 136, 120},
{255, 255, 240}, {238, 238, 224}, {205, 205, 193}, {139, 139, 131}, {240, 255, 240},
{224, 238, 224}, {193, 205, 193}, {131, 139, 131}, {255, 240, 245}, {238, 224, 229},
{205, 193, 197}, {139, 131, 134}, {255, 228, 225}, {238, 213, 210}, {205, 183, 181},
{139, 125, 123}, {240, 255, 255}, {224, 238, 238}, {193, 205, 205}, {131, 139, 139},
{131, 111, 255}, {122, 103, 238}, {105, 89, 205}, {71, 60, 139}, {72, 118, 255},
{67, 110, 238}, {58, 95, 205}, {39, 64, 139}, {0, 0, 255}, {0, 0, 238}, {0, 0, 205},
{0, 0, 139}, {30, 144, 255}, {28, 134, 238}, {24, 116, 205}, {16, 78, 139},
{99, 184, 255}, {92, 172, 238}, {79, 148, 205}, {54, 100, 139}, {0, 191, 255},
{0, 178, 238}, {0, 154, 205}, {0, 104, 139}, {135, 206, 255}, {126, 192, 238},
{108, 166, 205}, {74, 112, 139}, {176, 226, 255}, {164, 211, 238}, {141, 182, 205},
{96, 123, 139}, {198, 226, 255}, {185, 211, 238}, {159, 182, 205}, {108, 123, 139},
{202, 225, 255}, {188, 210, 238}, {162, 181, 205}, {110, 123, 139}, {191, 239, 255},
{178, 223, 238}, {154, 192, 205}, {104, 131, 139}, {224, 255, 255}, {209, 238, 238},
{180, 205, 205}, {122, 139, 139}, {187, 255, 255}, {174, 238, 238}, {150, 205, 205},
{102, 139, 139}, {152, 245, 255}, {142, 229, 238}, {122, 197, 205}, {83, 134, 139},
{0, 245, 255}, {0, 229, 238}, {0, 197, 205}, {0, 134, 139}, {0, 255, 255},
{0, 238, 238}, {0, 205, 205}, {0, 139, 139}, {151, 255, 255}, {141, 238, 238},
{121, 205, 205}, {82, 139, 139}, {127, 255, 212}, {118, 238, 198}, {102, 205, 170},
{69, 139, 116}, {193, 255, 193}, {180, 238, 180}, {155, 205, 155}, {105, 139, 105},
{84, 255, 159}, {78, 238, 148}, {67, 205, 128}, {46, 139, 87}, {154, 255, 154},
{144, 238, 144}, {124, 205, 124}, {84, 139, 84}, {0, 255, 127}, {0, 238, 118},
{0, 205, 102}, {0, 139, 69}, {0, 255, 0}, {0, 238, 0}, {0, 205, 0}, {0, 139, 0},
{127, 255, 0}, {118, 238, 0}, {102, 205, 0}, {69, 139, 0}, {192, 255, 62},
{179, 238, 58}, {154, 205, 50}, {105, 139, 34}, {202, 255, 112}, {188, 238, 104},
{162, 205, 90}, {110, 139, 61}, {255, 246, 143}, {238, 230, 133}, {205, 198, 115},
{139, 134, 78}, {255, 236, 139}, {238, 220, 130}, {205, 190, 112}, {139, 129, 76},
{255, 255, 224}, {238, 238, 209}, {205, 205, 180}, {139, 139, 122}, {255, 255, 0},
{238, 238, 0}, {205, 205, 0}, {139, 139, 0}, {255, 215, 0}, {238, 201, 0},
{205, 173, 0}, {139, 117, 0}, {255, 193, 37}, {238, 180, 34}, {205, 155, 29},
{139, 105, 20}, {255, 185, 15}, {238, 173, 14}, {205, 149, 12}, {139, 101, 8},
{255, 193, 193}, {238, 180, 180}, {205, 155, 155}, {139, 105, 105}, {255, 106, 106},
{238, 99, 99}, {205, 85, 85}, {139, 58, 58}, {255, 130, 71}, {238, 121, 66},
{205, 104, 57}, {139, 71, 38}, {255, 211, 155}, {238, 197, 145}, {205, 170, 125},
{139, 115, 85}, {255, 231, 186}, {238, 216, 174}, {205, 186, 150}, {139, 126, 102},
{255, 165, 79}, {238, 154, 73}, {205, 133, 63}, {139, 90, 43}, {255, 127, 36},
{238, 118, 33}, {205, 102, 29}, {139, 69, 19}, {255, 48, 48}, {238, 44, 44},
{205, 38, 38}, {139, 26, 26}, {255, 64, 64}, {238, 59, 59}, {205, 51, 51},
{139, 35, 35}, {255, 140, 105}, {238, 130, 98}, {205, 112, 84}, {139, 76, 57},
{255, 160, 122}, {238, 149, 114}, {205, 129, 98}, {139, 87, 66}, {255, 165, 0},
{238, 154, 0}, {205, 133, 0}, {139, 90, 0}, {255, 127, 0}, {238, 118, 0},
{205, 102, 0}, {139, 69, 0}, {255, 114, 86}, {238, 106, 80}, {205, 91, 69},
{139, 62, 47}, {255, 99, 71}, {238, 92, 66}, {205, 79, 57}, {139, 54, 38},
{255, 69, 0}, {238, 64, 0}, {205, 55, 0}, {139, 37, 0}, {255, 0, 0}, {238, 0, 0},
{205, 0, 0}, {139, 0, 0}, {255, 20, 147}, {238, 18, 137}, {205, 16, 118},
{139, 10, 80}, {255, 110, 180}, {238, 106, 167}, {205, 96, 144}, {139, 58, 98},
{255, 181, 197}, {238, 169, 184}, {205, 145, 158}, {139, 99, 108}, {255, 174, 185},
{238, 162, 173}, {205, 140, 149}, {139, 95, 101}, {255, 130, 171}, {238, 121, 159},
{205, 104, 137}, {139, 71, 93}, {255, 52, 179}, {238, 48, 167}, {205, 41, 144},
{139, 28, 98}, {255, 62, 150}, {238, 58, 140}, {205, 50, 120}, {139, 34, 82},
{255, 0, 255}, {238, 0, 238}, {205, 0, 205}, {139, 0, 139}, {255, 131, 250},
{238, 122, 233}, {205, 105, 201}, {139, 71, 137}, {255, 187, 255}, {238, 174, 238},
{205, 150, 205}, {139, 102, 139}, {224, 102, 255}, {209, 95, 238}, {180, 82, 205},
{122, 55, 139}, {191, 62, 255}, {178, 58, 238}, {154, 50, 205}, {104, 34, 139},
{155, 48, 255}, {145, 44, 238}, {125, 38, 205}, {85, 26, 139}, {171, 130, 255},
{159, 121, 238}, {137, 104, 205}, {93, 71, 139}, {255, 225, 255}, {238, 210, 238},
{205, 181, 205}, {139, 123, 139}, {0, 0, 0}, {0, 0, 0}, {3, 3, 3}, {3, 3, 3},
{5, 5, 5}, {5, 5, 5}, {8, 8, 8}, {8, 8, 8}, {10, 10, 10}, {10, 10, 10}, {13, 13, 13},
{13, 13, 13}, {15, 15, 15}, {15, 15, 15}, {18, 18, 18}, {18, 18, 18}, {20, 20, 20},
{20, 20, 20}, {23, 23, 23}, {23, 23, 23}, {26, 26, 26}, {26, 26, 26}, {28, 28, 28},
{28, 28, 28}, {31, 31, 31}, {31, 31, 31}, {33, 33, 33}, {33, 33, 33}, {36, 36, 36},
{36, 36, 36}, {38, 38, 38}, {38, 38, 38}, {41, 41, 41}, {41, 41, 41}, {43, 43, 43},
{43, 43, 43}, {46, 46, 46}, {46, 46, 46}, {48, 48, 48}, {48, 48, 48}, {51, 51, 51},
{51, 51, 51}, {54, 54, 54}, {54, 54, 54}, {56, 56, 56}, {56, 56, 56}, {59, 59, 59},
{59, 59, 59}, {61, 61, 61}, {61, 61, 61}, {64, 64, 64}, {64, 64, 64}, {66, 66, 66},
{66, 66, 66}, {69, 69, 69}, {69, 69, 69}, {71, 71, 71}, {71, 71, 71}, {74, 74, 74},
{74, 74, 74}, {77, 77, 77}, {77, 77, 77}, {79, 79, 79}, {79, 79, 79}, {82, 82, 82},
{82, 82, 82}, {84, 84, 84}, {84, 84, 84}, {87, 87, 87}, {87, 87, 87}, {89, 89, 89},
{89, 89, 89}, {92, 92, 92}, {92, 92, 92}, {94, 94, 94}, {94, 94, 94}, {97, 97, 97},
{97, 97, 97}, {99, 99, 99}, {99, 99, 99}, {102, 102, 102}, {102, 102, 102},
{105, 105, 105}, {105, 105, 105}, {107, 107, 107}, {107, 107, 107}, {110, 110, 110},
{110, 110, 110}, {112, 112, 112}, {112, 112, 112}, {115, 115, 115}, {115, 115, 115},
{117, 117, 117}, {117, 117, 117}, {120, 120, 120}, {120, 120, 120}, {122, 122, 122},
{122, 122, 122}, {125, 125, 125}, {125, 125, 125}, {127, 127, 127}, {127, 127, 127},
{130, 130, 130}, {130, 130, 130}, {133, 133, 133}, {133, 133, 133}, {135, 135, 135},
{135, 135, 135}, {138, 138, 138}, {138, 138, 138}, {140, 140, 140}, {140, 140, 140},
{143, 143, 143}, {143, 143, 143}, {145, 145, 145}, {145, 145, 145}, {148, 148, 148},
{148, 148, 148}, {150, 150, 150}, {150, 150, 150}, {153, 153, 153}, {153, 153, 153},
{156, 156, 156}, {156, 156, 156}, {158, 158, 158}, {158, 158, 158}, {161, 161, 161},
{161, 161, 161}, {163, 163, 163}, {163, 163, 163}, {166, 166, 166}, {166, 166, 166},
{168, 168, 168}, {168, 168, 168}, {171, 171, 171}, {171, 171, 171}, {173, 173, 173},
{173, 173, 173}, {176, 176, 176}, {176, 176, 176}, {179, 179, 179}, {179, 179, 179},
{181, 181, 181}, {181, 181, 181}, {184, 184, 184}, {184, 184, 184}, {186, 186, 186},
{186, 186, 186}, {189, 189, 189}, {189, 189, 189}, {191, 191, 191}, {191, 191, 191},
{194, 194, 194}, {194, 194, 194}, {196, 196, 196}, {196, 196, 196}, {199, 199, 199},
{199, 199, 199}, {201, 201, 201}, {201, 201, 201}, {204, 204, 204}, {204, 204, 204},
{207, 207, 207}, {207, 207, 207}, {209, 209, 209}, {209, 209, 209}, {212, 212, 212},
{212, 212, 212}, {214, 214, 214}, {214, 214, 214}, {217, 217, 217}, {217, 217, 217},
{219, 219, 219}, {219, 219, 219}, {222, 222, 222}, {222, 222, 222}, {224, 224, 224},
{224, 224, 224}, {227, 227, 227}, {227, 227, 227}, {229, 229, 229}, {229, 229, 229},
{232, 232, 232}, {232, 232, 232}, {235, 235, 235}, {235, 235, 235}, {237, 237, 237},
{237, 237, 237}, {240, 240, 240}, {240, 240, 240}, {242, 242, 242}, {242, 242, 242},
{245, 245, 245}, {245, 245, 245}, {247, 247, 247}, {247, 247, 247}, {250, 250, 250},
{250, 250, 250}, {252, 252, 252}, {252, 252, 252}, {255, 255, 255}, {255, 255, 255},
{169, 169, 169}, {169, 169, 169}, {169, 169, 169}, {169, 169, 169}, {0, 0, 139},
{0, 0, 139}, {0, 139, 139}, {0, 139, 139}, {139, 0, 139}, {139, 0, 139}, {139, 0, 0},
{139, 0, 0}, {144, 238, 144}, {144, 238, 144}};
public static int numberOfAvailableColors() {
return COLORS.length;
}
public static int[] getColorTripletByIndex( int index ) {
return COLORS_RGB[index];
}
public static Color getColor( String clr, Color defclr ) {
if (clr == null)
return defclr;
if (clr.startsWith("#")) {
/* Try and decode it as a RGB color specification */
try {
return Color.decode(clr);
} catch (NumberFormatException nfe) {
return defclr;
}
} else {
for( int i = 0; i < COLORS.length; i++ ) {
if (COLORS[i].equalsIgnoreCase(clr)) {
int[] iclr = COLORS_RGB[i];
return new Color(iclr[0], iclr[1], iclr[2]);
}
}
}
return defclr;
}
/**
* Returns the list of files involved in the raster map issues. If for example a map has to be
* deleted, then all these files have to.
*
* @param mapsetPath - the path of the mapset
* @param mapname -the name of the map
* @return the array of strings containing the full path to the involved files
*/
public static boolean checkRasterMapConsistence( String mapsetPath, String mapname ) {
File file = null;
File file2 = null;
file = new File(mapsetPath + File.separator + JGrassConstants.FCELL + File.separator
+ mapname);
file2 = new File(mapsetPath + File.separator + JGrassConstants.CELL + File.separator
+ mapname);
// the map is in one of the two
if (!file.exists() && !file2.exists())
return false;
/*
* helper files
*/
file = new File(mapsetPath + File.separator + JGrassConstants.CELLHD + File.separator
+ mapname);
if (!file.exists())
return false;
// it is important that the folder cell_misc/mapname comes before the
// files in it
file = new File(mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator
+ mapname);
if (!file.exists())
return false;
return true;
}
/**
* create a buffered image from a set of color triplets
*
* @param data
* @param width
* @param height
* @return
*/
public static BufferedImage ByteBufferImage( byte[] data, int width, int height ) {
int[] bandoffsets = {0, 1, 2, 3};
DataBufferByte dbb = new DataBufferByte(data, data.length);
WritableRaster wr = Raster.createInterleavedRaster(dbb, width, height, width * 4, 4,
bandoffsets, null);
int[] bitfield = {8, 8, 8, 8};
ColorSpace cs = ColorSpace.getInstance(ColorSpace.CS_sRGB);
ColorModel cm = new ComponentColorModel(cs, bitfield, true, false,
Transparency.TRANSLUCENT, DataBuffer.TYPE_BYTE);
return new BufferedImage(cm, wr, false, null);
}
public static Envelope reprojectEnvelopeByEpsg( int srcEpsg, int destEpsg, Envelope srcEnvelope )
throws FactoryException, TransformException {
CoordinateReferenceSystem sourceCRS = CRS.decode("EPSG:" + srcEpsg); //$NON-NLS-1$
CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:" + destEpsg); //$NON-NLS-1$
MathTransform tr = CRS.findMathTransform(sourceCRS, targetCRS);
// From that point, I'm not sure which kind of object is returned by
// getLatLonBoundingBox(). But there is some convenience methods if CRS
// like:
return JTS.transform(srcEnvelope, tr);
}
/**
* return the rectangle of the cell of the active region, that surrounds the given coordinates
*
* @param activeRegion
* @param x the given easting coordinate
* @param y given northing coordinate
* @return the rectangle localizing the cell inside which the x and y stay
*/
public static JGrassRegion getRectangleAroundPoint( JGrassRegion activeRegion, double x,
double y ) {
double minx = activeRegion.getRectangle().getBounds2D().getMinX();
double ewres = activeRegion.getWEResolution();
double snapx = minx + (Math.round((x - minx) / ewres) * ewres);
double miny = activeRegion.getRectangle().getBounds2D().getMinY();
double nsres = activeRegion.getNSResolution();
double snapy = miny + (Math.round((y - miny) / nsres) * nsres);
double xmin = 0.0;
double xmax = 0.0;
double ymin = 0.0;
double ymax = 0.0;
if (x >= snapx) {
xmin = snapx;
xmax = xmin + ewres;
} else {
xmax = snapx;
xmin = xmax - ewres;
}
if (y <= snapy) {
ymax = snapy;
ymin = ymax - nsres;
} else {
ymin = snapy;
ymax = ymin + nsres;
}
// why do I have to put ymin, Rectangle requires the upper left
// corner?????!!!! Is it a BUG
// in the Rectangle2D class? or docu?
// Rectangle2D rect = new Rectangle2D.Double(xmin, ymin, ewres, nsres);
return new JGrassRegion(xmin, xmax, ymin, ymax, 1, 1);
}
// /**
// * This method gives the possibility to supply a point through coordinates
// and gain the value
// of
// * a certain supplied map in that point (a kind of d.what.rast)
// *
// * @param mapName - the map of which we want to know the value ("@mapset")
// * @param copt
// * @param rasterFormat - the raster format from where to read from
// * @param coordinates - the point in which we want to know the value
// (x=easting, y=northing)
// * @return
// */
// public static String getWhatValueFromMap( String mapName, Point2D
// coordinates ) {
//
// Rectangle2D pixelBound =
// dataMask.snapRectangleToClickedCell(coordinates.getX(),
// coordinates.getY());
//
// Window active = new Window(pixelBound.getMinX(), pixelBound.getMaxX(),
// pixelBound.getMinY(), pixelBound.getMaxY(), 1, 1);
//
// // define the map
// MapReader reader = MapIOFactory.CreateRasterMapReader(rasterFormat);
// RasterMap theMap = (RasterMap) copt.g.getLocation().getMaps(reader,
// mapName).elementAt(0); //
// get
// // the
// // rastermap
// theMap.setReader(rasterFormat); // the type of reader has to be set
// theMap.setNovalue(new java.lang.Double(FluidConstants.fluidnovalue)); //
// the
// // novalue
// // has to be
// // set
// theMap.setDataWindow(active); // the calculation region is set
// theMap.setOutputMatrixType(0); // set the matrix type
// theMap.setOutputDataObject(new double[0][0]); // data type
//
// // read elevation
// double[][] elevation = null;
// if (!theMap.openMap()) {
// JOptionPane.showMessageDialog(null, "Error! Problems opening map...");
// return "NaN";
// }
// while( theMap.hasMoreData() ) {
// elevation = (double[][]) theMap.getNextData();
// }
//
// if (elevation[0][0] == FluidConstants.fluidnovalue)
// return "null";
// return String.valueOf(elevation[0][0]);
// }
//
// /**
// * When we click on the map the chosen point will be somewhere inside the
// hit cell. It is a
// * matter of rounding the northing and easting of the point to get the
// right row and column
// that
// * contains the hit. Therefore this method returns the row and col in
// which the point falls.
// *
// * @param loc - the location
// * @param coordinates - the coordinates of the hit point
// * @return the reviewed point as row, col
// */
// public static Point putClickToCenterOfCell( Window active, Point2D
// coordinates ) {
// double eastingClick = coordinates.getX();
// double northingClick = coordinates.getY();
//
// double startnorth = active.getNorth();
// double startwest = active.getWest();
// double northdelta = active.getNSResolution();
// double westdelta = active.getWEResolution();
// int clickrow = 0;
// int clickcol = 0;
//
// for( int i = 0; i < active.getRows(); i++ ) {
// startnorth = startnorth - northdelta;
// if (northingClick > startnorth) {
// clickrow = i;
// break;
// }
// }
// for( int i = 0; i < active.getCols(); i++ ) {
// startwest = startwest + westdelta;
// if (eastingClick < startwest) {
// clickcol = i;
// break;
// }
// }
//
// return new Point(clickrow, clickcol);
// }
//
/**
* <p>
* Transforms row and column index of the active region into the regarding northing and easting
* coordinates. The center of the cell is taken.
* </p>
* <p>
* NOTE: basically the inverse of
* {@link JGrassUtilities#coordinateToNearestRowCol(JGrassRegion, Coordinate)}
* </p>
*
* @param active - the active region (can be null)
* @param row - row number of the point to transform
* @param col - column number of the point to transform
* @return the point in N/E coordinates of the supplied row and column
*/
public static Coordinate rowColToCenterCoordinates( JGrassRegion active, int row, int col ) {
double north = active.getNorth();
double west = active.getWest();
double nsres = active.getNSResolution();
double ewres = active.getWEResolution();
double northing = north - row * nsres - nsres / 2.0;
double easting = west + col * ewres + ewres / 2.0;
return new Coordinate(easting, northing);
}
/**
* <p>
* Return the row and column of the active region matrix for a give coordinate *
* </p>
* <p>
* NOTE: basically the inverse of
* {@link JGrassUtilities#rowColToCenterCoordinates(JGrassRegion, int, int)}
* </p>
*
* @param active the active region
* @param coord
* @return and int array containing row and col
*/
public static int[] coordinateToNearestRowCol( JGrassRegion active, Coordinate coord ) {
double easting = coord.x;
double northing = coord.y;
int[] rowcol = new int[2];
if (easting > active.getEast() || easting < active.getWest()
|| northing > active.getNorth() || northing < active.getSouth()) {
return null;
}
double minx = active.getWest();
double ewres = active.getWEResolution();
for( int i = 0; i < active.getCols(); i++ ) {
minx = minx + ewres;
if (easting < minx) {
rowcol[1] = i;
break;
}
}
double maxy = active.getNorth();
double nsres = active.getNSResolution();
for( int i = 0; i < active.getRows(); i++ ) {
maxy = maxy - nsres;
if (northing > maxy) {
rowcol[0] = i;
break;
}
}
return rowcol;
}
// /**
// * Given the mapsetpath and the mapname, the map is removed with all its accessor files
// *
// * @param mapsetPath
// * @param mapName
// */
// public static boolean removeGrassRasterMap( String mapsetPath, String mapName ) {
//
// // list of files to remove
// String mappaths[] = filesOfRasterMap(mapsetPath, mapName);
//
// // first delete the list above, which are just files
// for( int j = 0; j < mappaths.length; j++ ) {
// File filetoremove = new File(mappaths[j]);
// if (filetoremove.exists()) {
// if (!FileUtilities.deleteFileOrDir(filetoremove)) {
// return false;
// }
// }
// }
// return true;
// }
/**
* Returns the list of files involved in the raster map issues. If for example a map has to be
* deleted, then all these files have to.
*
* @param mapsetPath - the path of the mapset
* @param mapname -the name of the map
* @return the array of strings containing the full path to the involved files
*/
public static String[] filesOfRasterMap( String mapsetPath, String mapname ) {
String filesOfRaster[] = new String[]{
mapsetPath + File.separator + JGrassConstants.FCELL + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.CELL + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.CATS + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.HIST + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.CELLHD + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.COLR + File.separator + mapname,
// it is very important that the folder cell_misc/mapname comes
// before the files in it
mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator + mapname,
mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator + mapname
+ File.separator + JGrassConstants.CELLMISC_FORMAT,
mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator + mapname
+ File.separator + JGrassConstants.CELLMISC_QUANT,
mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator + mapname
+ File.separator + JGrassConstants.CELLMISC_RANGE,
mapsetPath + File.separator + JGrassConstants.CELL_MISC + File.separator + mapname
+ File.separator + JGrassConstants.CELLMISC_NULL};
return filesOfRaster;
}
/**
* Transforms row and column index of the active region into an array of the coordinates of the
* edgaes, i.e. n, s, e, w
*
* @param active - the active region (can be null)
* @param row - row number of the point to transform
* @param col - column number of the point to transform
* @return the array of north, south, east, west
*/
public static double[] rowColToNodeboundCoordinates( JGrassRegion active, int row, int col ) {
double anorth = active.getNorth();
double awest = active.getWest();
double nsres = active.getNSResolution();
double ewres = active.getWEResolution();
double[] nsew = new double[4];
nsew[0] = anorth - row * nsres;
nsew[1] = anorth - row * nsres - nsres;
nsew[2] = awest + col * ewres + ewres;
nsew[3] = awest + col * ewres;
return nsew;
}
public static int factorial( int n ) {
int fact = 1;
for( int i = 1; i <= n; i++ ) {
fact *= i;
}
return fact;
}
public static void makeColorRulesPersistent( File colrFile, List<String> rules,
double[] minMax, int alpha ) throws IOException {
if (!colrFile.getParentFile().exists()) {
colrFile.getParentFile().mkdir();
}
BufferedWriter bw = new BufferedWriter(new FileWriter(colrFile));
if (rules.size() == 0) {
throw new IllegalArgumentException("The list of colorrules can't be empty.");
}
String header = "% " + minMax[0] + " " + minMax[1] + " " + alpha;
bw.write(header + "\n");
for( String r : rules ) {
bw.write(r + "\n");
}
bw.close();
}
/**
* Calculates optimal tile size for the actual free memory.
*
* @param rows the rows of the complete image the tiles are calculated for.
* @param cols the cols of the complete image the tiles are calculated for.
* @return
*/
public static int[] getTilesBasedOnFreeMemory( int rows, int cols ) {
long freeMemory = Runtime.getRuntime().freeMemory();
int tileSizeY = 256;
int tileSizeX = 256;
if (freeMemory > 8L * cols) {
tileSizeX = cols;
tileSizeY = (int) (freeMemory / 8) / cols;
if (tileSizeY > rows) {
tileSizeY = rows;
}
}
return new int[]{tileSizeX, tileSizeY};
}
public static JGrassRegion getJGrassRegionFromGridCoverage( GridCoverage2D gridCoverage2D )
throws InvalidGridGeometryException, TransformException {
Envelope2D env = gridCoverage2D.getEnvelope2D();
GridEnvelope2D worldToGrid = gridCoverage2D.getGridGeometry().worldToGrid(env);
double xRes = env.getWidth() / worldToGrid.getWidth();
double yRes = env.getHeight() / worldToGrid.getHeight();
JGrassRegion region = new JGrassRegion(env.getMinX(), env.getMaxX(), env.getMinY(), env
.getMaxY(), xRes, yRes);
return region;
}
}