/// <summary>
/// Returns MarkerType property for legend.
/// </summary>
/// <param name="colorRepo"></param>
/// <param name="resultDir"></param>
/// <param name="fi"></param>
/// <param name="resultName"></param>
/// <returns></returns>
private static void ProcessVectorFile(IColorRepository colorRepo, DirectoryInfo resultDir, FileInfo fi, string resultName)
{
string localFileName = fi.Name;
string localResultName = resultName;
Console.WriteLine(string.Format("Processing {0} into {1}...", localFileName, localResultName));
StringBuilder bldr = new StringBuilder();
NetTopologySuite.IO.ShapefileDataReader dataReader = new NetTopologySuite.IO.ShapefileDataReader(fi.FullName, new GeometryFactory());
NetTopologySuite.Features.FeatureCollection featureCollection = new NetTopologySuite.Features.FeatureCollection();
List<string> csvHdr = dataReader.DbaseHeader.Fields.Select(a => a.Name).ToList();
csvHdr.Add(appColorNamspace);
bldr.AppendLine(string.Join(",", csvHdr)); //write csv file header
while (dataReader.Read())
{
NetTopologySuite.Features.Feature feature = new NetTopologySuite.Features.Feature();
feature.Geometry = dataReader.Geometry;
int numFields = dataReader.DbaseHeader.NumFields + 1;
string[] keys = new string[numFields];
int colorValueField = -1;
for (int i = 0; i < numFields - 1; i++)
{
keys[i] = dataReader.DbaseHeader.Fields[i].Name;
if (keys[i].Equals(colorRepo.ColorFieldForOutput(localFileName, localResultName)))
{
colorValueField = i;
}
}
//add attributes from source attribute table
feature.Attributes = new AttributesTable();
List<string> csvLine = new List<string>();
for (int i = 0; i < numFields - 1; i++)
{
object val = dataReader.GetValue(i);
feature.Attributes.AddAttribute(keys[i], val);
csvLine.Add(val.ToString());
}
if (colorRepo.MapColorsToThisResult(localFileName, localResultName))
{
//mark outline colors in a different attribute than fill colors
string colorNs = colorRepo.IsOutlinedNotFilled(localFileName, localResultName) ? appOutlineNamespace : appColorNamspace;
keys[numFields - 1] = colorNs;
//add additional attribute for color binding
string hexClr = colorRepo.SingleColorForFile(localFileName, localResultName); //only path where colorValueField, i.e. ColorMap.clrField can be unpopulated.
if (string.IsNullOrEmpty(hexClr) && colorValueField > -1)
{
if (colorRepo.IsCategoricalMap(localFileName, resultName))
{
//categorical color map
hexClr = colorRepo.ColorsOfValueInFile(localFileName, localResultName, dataReader.GetString(colorValueField)).HexColor;
}
else
{
//numerical range color map
hexClr = colorRepo.ColorsOfValueInFile(localFileName, localResultName, dataReader.GetDouble(colorValueField)).HexColor;
}
}
if (string.IsNullOrEmpty(hexClr)) // else if (string.IsNullOrEmpty(hexClr) && colorValueField < 0)
{
throw new NotSupportedException("Cannot color a file with no attributes to bind to and no single-color given.");
}
csvLine.Add(hexClr);
feature.Attributes.AddAttribute(colorNs, hexClr);
}
bldr.AppendLine(string.Join(",", csvLine));
featureCollection.Add(feature);
}
GeoJsonWriter wtr = new GeoJsonWriter();
string layerJson = wtr.Write(featureCollection);
File.WriteAllText(resultDir.FullName + localResultName, layerJson);
File.WriteAllText(resultDir.FullName + localResultName.Replace(".json", ".csv"), bldr.ToString());
}
/// <summary>
/// Processes a raster file stored as an ESRI Grid into an 8-bit RGB .tif file based on a value-color map. Pass the function the hdr.adf file.
/// </summary>
/// <param name="args"></param>
/// <param name="colorRepo"></param>
/// <param name="resultDir"></param>
/// <param name="srcDrv"></param>
/// <param name="fi">FileInfo of the hdr.adf file of the grid.</param>
/// <param name="resultName"></param>
private static void ProcessRasterGrid(string[] args, IColorRepository colorRepo, DirectoryInfo resultDir, OSGeo.GDAL.Driver srcDrv, FileInfo fi, string resultName)
{
Console.WriteLine(string.Format("Processing {0}.adf into {1}...", fi.Directory.Name, resultName));
StringBuilder bldr = new StringBuilder();
//read data from source raster
Dataset src = Gdal.Open(fi.FullName, Access.GA_ReadOnly);
Band band = src.GetRasterBand(1);
double[] r = new double[band.XSize * band.YSize];
byte[] red = new byte[band.XSize * band.YSize];
byte[] green = new byte[band.XSize * band.YSize];
byte[] blue = new byte[band.XSize * band.YSize];
band.ReadRaster(0, 0, band.XSize, band.YSize, r, band.XSize, band.YSize, 0, 0);
//assign values to rgb rasters to produce new raster with rgb bands matching color pattern
for (int cell = 0; cell < r.Length; cell++)
{
RGBColors colors = colorRepo.ColorsOfValueInFile(fi.Directory.Name, resultName, r[cell]);
red[cell] = (byte)colors.Red;
green[cell] = (byte)colors.Green;
blue[cell] = (byte)colors.Blue;
}
//write new output
using (Dataset output = srcDrv.Create(resultDir + resultName, band.XSize, band.YSize, 3, DataType.GDT_Byte, null))
{
if (output == null)
{
Console.WriteLine("Can't create " + args[0]);
System.Environment.Exit(-1);
}
//set metadata
output.SetProjection(src.GetProjection());
double[] geotransform = new double[0];
src.GetGeoTransform(geotransform);
output.SetGeoTransform(geotransform);
//prepare data for write
int[] colorData = new int[red.Length * 3];
red.CopyTo(colorData, 0);
green.CopyTo(colorData, red.Length);
blue.CopyTo(colorData, red.Length + green.Length);
//write data to disk
output.WriteRaster(0, 0, band.XSize, band.YSize, colorData, band.XSize, band.YSize, 3, null, 0, 0, 0);
output.FlushCache();
}
}
