public static IRasterRenderer ClassifyRenderer(IRasterLayer raslyr,IRgbColor start_clr,IRgbColor end_clr,int count)
{
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = pRClassRend as IRasterRenderer;
IRaster pRaster = raslyr.Raster;
IRasterBandCollection pRBandCol = pRaster as IRasterBandCollection;
IRasterBand pRBand = pRBandCol.Item(0);
if (pRBand.Histogram == null)
{
pRBand.ComputeStatsAndHist();
}
rasterRenderer.Raster = pRaster;
pRClassRend.ClassCount = count;
rasterRenderer.Update();
IAlgorithmicColorRamp colorRamp = new AlgorithmicColorRampClass();
colorRamp.Size = count;
colorRamp.FromColor = start_clr;
colorRamp.ToColor = end_clr;
bool createColorRamp;
colorRamp.CreateRamp(out createColorRamp);
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < pRClassRend.ClassCount; i++)
{
fillSymbol.Color = colorRamp.get_Color(i);
pRClassRend.set_Symbol(i, fillSymbol as ISymbol);
pRClassRend.set_Label(i, pRClassRend.get_Break(i).ToString("0.00"));
}
return rasterRenderer;
}
public static void classifyRender(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
try
{
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
classify = new EqualIntervalClass();
break;
case "自然断点分级":
classify = new NaturalBreaksClass();
break;
}
classify.Classify(ref numClasses);
double[] Classes = classify.ClassBreaks as double[];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassField = "Value";
rasClassifyRender.ClassCount = count;
rasRender.Update();
IColor pColor;
ISimpleFillSymbol pSym;
for (int j = 0; j < count; j++)
{
pColor = ramp.get_Color(j * (ramp.Size - 1) / (count - 1));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasClassifyRender.set_Break(j, rasClassifyRender.get_Break(j));
}
rasRender.Update();
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
catch
{
XtraMessageBox.Show("唯一值数量已达到限制（65536）");
}
}
public static IRasterRenderer CreateSlopePrecentRiseColorRamp(Raster gRaster)
{
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
IRaster raster = rasterDataset.CreateDefaultRaster();
rasterRenderer.Raster = raster;
classifyRenderer.ClassCount = 10;
rasterRenderer.Update();
classifyRenderer.Break[0] = 0;
classifyRenderer.Label[0] = "0 to 3.5%";
classifyRenderer.Break[1] = 3.5;
classifyRenderer.Label[1] = "3.5% to 8.75%";
classifyRenderer.Break[2] = 8.75;
classifyRenderer.Label[2] = "8.75% to 15%";
classifyRenderer.Break[3] = 15;
classifyRenderer.Label[3] = "15% to 25%";
classifyRenderer.Break[4] = 25;
classifyRenderer.Label[4] = "25% to 45%";
classifyRenderer.Break[5] = 45;
classifyRenderer.Label[5] = "45% to 70%";
classifyRenderer.Break[6] = 70;
classifyRenderer.Label[6] = "70% to 100%";
classifyRenderer.Break[7] = 100;
classifyRenderer.Label[7] = "100% to 175%";
classifyRenderer.Break[8] = 175;
classifyRenderer.Label[8] = "175% to 565%";
classifyRenderer.Break[9] = 565;
classifyRenderer.Label[9] = "> 565%";
List <RgbColor> lColors = CreateSlopeColorRamp();
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
return(rasterRenderer);
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
public static IRasterRenderer CreateGrainSizeStatisticColorRamp(GCDConsoleLib.Raster gRaster, UnitsNet.Units.LengthUnit eUnits)
{
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
IRaster raster = rasterDataset.CreateDefaultRaster();
rasterRenderer.Raster = raster;
classifyRenderer.ClassCount = 5;
rasterRenderer.Update();
classifyRenderer.Break[0] = 0;
classifyRenderer.Label[0] = "Fines, Sand (0 to 2 mm)";
classifyRenderer.Break[1] = UnitsNet.Length.From(2, UnitsNet.Units.LengthUnit.Millimeter).As(eUnits);
classifyRenderer.Label[1] = "Fine Gravel (2 mm to 16 mm)";
classifyRenderer.Break[2] = UnitsNet.Length.From(16, UnitsNet.Units.LengthUnit.Millimeter).As(eUnits);
classifyRenderer.Label[2] = "Coarse Gravel (16 mm to 64 mm)";
classifyRenderer.Break[3] = UnitsNet.Length.From(64, UnitsNet.Units.LengthUnit.Millimeter).As(eUnits);
classifyRenderer.Label[3] = "Cobbles (64 mm to 256 mm)";
classifyRenderer.Break[4] = UnitsNet.Length.From(256, UnitsNet.Units.LengthUnit.Millimeter).As(eUnits);
classifyRenderer.Label[4] = "Boulders (> 256 mm)";
List <RgbColor> lColors = new List <RgbColor>();
lColors.Add(CreateRGBColor(194, 82, 60));
lColors.Add(CreateRGBColor(240, 180, 17));
lColors.Add(CreateRGBColor(123, 237, 0));
lColors.Add(CreateRGBColor(27, 168, 124));
lColors.Add(CreateRGBColor(11, 44, 122));
IFillSymbol fillSymbol = new SimpleFillSymbol();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
return(rasterRenderer);
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
public static IRasterRenderer CreateClassifyRenderer(GCDConsoleLib.Raster gRaster, int iClassCount, string sColorRampName, double dMax, bool bInvert = false)
{
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
IFillSymbol fillSymbol = new SimpleFillSymbol();
IRaster raster = rasterDataset.CreateDefaultRaster();
rasterRenderer.Raster = raster;
IColorRamp pColorRamp;
IStyleGalleryItem pStyleItem = GetESRIStyleColorRamp(out pColorRamp, sColorRampName);
classifyRenderer.ClassCount = iClassCount;
rasterRenderer.Update();
CreateClassBreaks(dMax, iClassCount, classifyRenderer);
pColorRamp.Size = iClassCount;
bool bOK = false;
pColorRamp.CreateRamp(out bOK);
List <IColor> lColors = new List <IColor>();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
lColors.Add(pColorRamp.Color[i]);
}
if (bInvert)
{
lColors.Reverse();
}
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
return(rasterRenderer);
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
public static IRasterRenderer ClassifyRenderer(ESRI.ArcGIS.Geodatabase.IRasterDataset pRasterDataset)
{
try
{
//create the classify render:创建渲染器
IRasterClassifyColorRampRenderer pClassifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasterRenderer = (IRasterRenderer)pClassifyRenderer;
IRaster pRaster = pRasterDataset.CreateDefaultRaster();
pRasterRenderer.Raster = pRaster;
pClassifyRenderer.ClassCount = 10;
pRasterRenderer.Update();
//创建颜色
IRgbColor pFromColor = new RgbColorClass();
pFromColor.Red = 255;
pFromColor.Green = 0;
pFromColor.Blue = 0;
IRgbColor pToColor = new RgbColorClass();
pToColor.Red = 0;
pToColor.Green = 255;
pToColor.Blue = 255;
//set the color ramp for the symbology:生成色带
IAlgorithmicColorRamp pRamp = new AlgorithmicColorRampClass();
pRamp.Size = 10;
pRamp.FromColor = pFromColor;
pRamp.ToColor = pToColor;
bool pBoolColorRamp;
pRamp.CreateRamp(out pBoolColorRamp);
//create the symbol for the classes.:创建符号
IFillSymbol pFillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < pClassifyRenderer.ClassCount; i++)
{
pFillSymbol.Color = pRamp.get_Color(i);
pClassifyRenderer.set_Symbol(i, (ISymbol)pFillSymbol);
pClassifyRenderer.set_Label(i, Convert.ToString(i));
}
return(pRasterRenderer);
}
catch (System.Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
/// <summary>
/// 获取手动分级渲染设置
/// </summary>
/// <param name="classBreaks">分级数</param>
/// <param name="classField">分级字段</param>
/// <returns></returns>
public IRasterClassifyColorRampRenderer GetManualRasterClassifyRenderer(List <double> classBreaks,
string classField)
{
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
//字段
pRClassRend.ClassField = classField;
//类别总数
pRClassRend.ClassCount = classBreaks.Count / 2;
//设置第一个分级数
pRClassRend.Break[0] = classBreaks[0];
//每行第二个数值为分级数
for (int i = 0; i < classBreaks.Count / 2; i++)
{
pRClassRend.Break[i + 1] = classBreaks[2 * i + 1];
}
return(pRClassRend);
}
public static void funColorForRaster_Classify(IRasterLayer pRasterLayer, int ClassNum)
{
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRRend = pRClassRend as IRasterRenderer;
IRaster pRaster = pRasterLayer.Raster;
IRasterBandCollection pRBandCol = pRaster as IRasterBandCollection;
IRasterBand pRBand = pRBandCol.Item(0);
if (pRBand.Histogram == null)
{
pRBand.ComputeStatsAndHist();
}
pRRend.Raster = pRaster;
pRClassRend.ClassCount = ClassNum;
pRRend.Update();
IRgbColor pFromColor = new RgbColorClass();
pFromColor.Red = 0;//绿
pFromColor.Green = 255;
pFromColor.Blue = 0;
IRgbColor pToColor = new RgbColorClass();
pToColor.Red = 255;//红
pToColor.Green = 0;
pToColor.Blue = 0;
IAlgorithmicColorRamp colorRamp = new AlgorithmicColorRampClass();
colorRamp.Size = ClassNum;
colorRamp.FromColor = pFromColor;
colorRamp.ToColor = pToColor;
bool createColorRamp;
colorRamp.CreateRamp(out createColorRamp);
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < pRClassRend.ClassCount; i++)
{
fillSymbol.Color = colorRamp.get_Color(i);
pRClassRend.set_Symbol(i, fillSymbol as ISymbol);
pRClassRend.set_Label(i, (i + 1).ToString());
}
pRasterLayer.Renderer = pRRend;
MainFrom.m_mapControl.AddLayer(pRasterLayer);
}
public static void userRender(IRasterLayer rastlayer, IColorRamp ramp, int count, double[] values)
{
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
//rasClassifyRender.ClassField = "Value";
rasClassifyRender.ClassCount = count;
rasRender.Update();
IColor pColor;
ISimpleFillSymbol pSym;
int j;
for (j = 0; j < count; j++)
{
pColor = ramp.get_Color(j * (ramp.Size - 1) / (count - 1));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasClassifyRender.set_Break(j, values[j - 1 < 0 ? 0 : j - 1]);
if (j == 0)
{
rasClassifyRender.set_Label(j, string.Format("{0:#######0.#####}-{1:#######0.#####}", 0, values[j]));
}
else
{
rasClassifyRender.set_Label(j, string.Format("{0:#######0.#####}-{1:#######0.#####}", values[j - 1], values[j]));
}
Console.WriteLine("自定义" + values[j] + "," + rasClassifyRender.get_Break(j));
}
rasClassifyRender.set_Break(j, values[j - 1]);
rasRender.Update();
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
/// <summary>
/// 获取分级渲染设置
/// </summary>
/// <param name="classField">字段</param>
/// <param name="classCount">类别总数</param>
/// <param name="classMethod">分级算法</param>
/// <returns></returns>
public IRasterClassifyColorRampRenderer GetRasterClassifyRenderer(string classField, int classCount,
string classMethod)
{
//先计算统计值，不然无法获取分类数据
CalculateRasterStatistics(rasterLayer);
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
//分级属性接口
IRasterClassifyUIProperties classProperties = pRClassRend as IRasterClassifyUIProperties;
//设置分级算法
SetClassMethod(classProperties, classMethod);
//字段
pRClassRend.ClassField = classField;
//类别总数
pRClassRend.ClassCount = classCount;
IRasterRenderer rasterRend = (IRasterRenderer)pRClassRend;
rasterRend.Raster = rasterLayer.Raster;
rasterRend.Update();
return(pRClassRend);
}
//Ϊդ��ͼ�����ݰ��մ���Сֵ�����ֵ�ķ�����ɫ
public static IRasterLayer SetRsLayerClassifiedColor(IRaster pRaster)
{
IRasterClassifyColorRampRenderer pClassRen = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasRen = pClassRen as RasterClassifyColorRampRendererClass;
//Set raster for the render and update
pRasRen.Raster = pRaster;
pClassRen.ClassCount = 9;
pRasRen.Update();
//Create a color ramp to use
//���������յ���ɫ
IColor pFromColor = new RgbColorClass();
IRgbColor pRgbColor = pFromColor as IRgbColor;
pRgbColor.Red = 255;
pRgbColor.Green = 200;
pRgbColor.Blue = 0;
IColor pToColor = new RgbColorClass();
pRgbColor = pToColor as IRgbColor;
pRgbColor.Red = 0;
pRgbColor.Green = 0;
pRgbColor.Blue = 255;
//������ɫ�ּ�
IAlgorithmicColorRamp pRamp = new AlgorithmicColorRampClass();
pRamp.Size = 9;
pRamp.FromColor = pFromColor;
pRamp.ToColor = pToColor;
bool ok = true;
pRamp.CreateRamp(out ok);
//���դ��ͳ����ֵ
IRasterBandCollection pRsBandCol = pRaster as IRasterBandCollection;
IRasterBand pRsBand = pRsBandCol.Item(0);
pRsBand.ComputeStatsAndHist();
IRasterStatistics pRasterStatistic = pRsBand.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
//Create symbol for the classes
IFillSymbol pFSymbol = new SimpleFillSymbolClass();
double LabelValue = Convert.ToDouble((dMaxValue - dMinValue) / 9);
for (int i = 0; i <= pClassRen.ClassCount - 1; i++)
{
pFSymbol.Color = pRamp.get_Color(i);
pClassRen.set_Symbol(i, pFSymbol as ISymbol);
double h1 = (LabelValue * i) + dMinValue;
double h2 = (LabelValue * (i + 1)) + dMinValue;
pClassRen.set_Label(i, h1.ToString() + "��" + h2.ToString());
}
//Update the renderer and plug into layer
pRasRen.Update();
IRasterLayer pRLayer = new RasterLayerClass();
pRLayer.CreateFromRaster(pRaster);
pRLayer.Renderer = pClassRen as IRasterRenderer;
return pRLayer;
}
private IRasterLayer SetRsLayerClassifiedColor(IRaster pRaster)
{
IRasterClassifyColorRampRenderer pClassRen = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasRen = (IRasterRenderer)pClassRen;
//Set raster for the render and update
pRasRen.Raster = pRaster;
pClassRen.ClassCount = 9;
pRasRen.Update();
//Create a color ramp to use
//定义起点和终点颜色
IColor pFromColor = new RgbColorClass();
IRgbColor pRgbColor = (IRgbColor)pFromColor;
pRgbColor.Red = 255;
pRgbColor.Green = 200;
pRgbColor.Blue = 0;
IColor pToColor = new RgbColorClass();
pRgbColor = (IRgbColor)pToColor;
pRgbColor.Red = 0;
pRgbColor.Green = 0;
pRgbColor.Blue = 255;
//创建颜色分级
IAlgorithmicColorRamp pRamp = new AlgorithmicColorRampClass();
pRamp.Size = 9;
pRamp.FromColor = pFromColor;
pRamp.ToColor = pToColor;
bool ok = true;
pRamp.CreateRamp(out ok);
//获得栅格统计数值
IRasterBandCollection pRsBandCol = (IRasterBandCollection)pRaster;
IRasterBand pRsBand = pRsBandCol.Item(0);
pRsBand.ComputeStatsAndHist();
IRasterStatistics pRasterStatistic = pRsBand.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
//Create symbol for the classes
IFillSymbol pFSymbol = new SimpleFillSymbolClass();
double LabelValue = Convert.ToDouble((dMaxValue - dMinValue) / 9);
for (int i = 0; i < pClassRen.ClassCount; i++)
{
pFSymbol.Color = pRamp.get_Color(i);
pClassRen.set_Symbol(i, (ISymbol)pFSymbol);
double h1 = (LabelValue * i) + dMinValue;
double h2 = (LabelValue * (i + 1)) + dMinValue;
pClassRen.set_Label(i, h1.ToString() + "-" + h2.ToString());
}
//Update the renderer and plug into layer
pRasRen.Update();
IRasterLayer pRLayer = new RasterLayerClass();
pRLayer.CreateFromRaster(pRaster);
pRLayer.Renderer = (IRasterRenderer)pClassRen;
return(pRLayer);
}
public static IRasterRenderer CreateDoDClassifyRenderer(Raster gRaster, int iClassCount, double rampRange)
{
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
IRaster raster = rasterDataset.CreateDefaultRaster();
double rMin = -2.5;
double rMax = 2.5;
if (rampRange == 0)
{
gRaster.ComputeStatistics();
Dictionary <string, decimal> rasterStats = gRaster.GetStatistics();
rMin = (double)rasterStats["min"];
rMax = (double)rasterStats["max"];
}
else
{
rMin = rampRange * -1;
rMax = rampRange;
}
rasterRenderer.Raster = raster;
if ((rMin == double.MinValue & rMax == double.MaxValue) | (rMin == double.MaxValue & rMax == double.MinValue) | (rMin == float.MinValue & rMax == float.MaxValue) | (rMin == float.MaxValue & rMax == float.MinValue))
{
classifyRenderer.ClassCount = 1;
ESRI.ArcGIS.Display.IRgbColor rgbColor = CreateRGBColor(255, 255, 255);
rgbColor.Transparency = 0;
fillSymbol.Color = rgbColor;
classifyRenderer.Symbol[0] = (ISymbol)fillSymbol;
classifyRenderer.Label[0] = "No Data (no change detected)";
return(rasterRenderer);
}
classifyRenderer.ClassCount = iClassCount;
rasterRenderer.Update();
CreateDoDClassBreaks(rMax, rMin, iClassCount, ref classifyRenderer);
List <IColor> lColors = CreateDoDColorRamp();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
return(rasterRenderer);
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
public static IRasterRenderer CreateRoughnessColorRamp(GCDConsoleLib.Raster gRaster)
{
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
IRaster raster = rasterDataset.CreateDefaultRaster();
rasterRenderer.Raster = raster;
classifyRenderer.ClassCount = 10;
rasterRenderer.Update();
classifyRenderer.Break[0] = 0;
classifyRenderer.Label[0] = "0 to 0.1";
classifyRenderer.Break[1] = 0.1;
classifyRenderer.Label[1] = "0.1 to 0.25";
classifyRenderer.Break[2] = 0.25;
classifyRenderer.Label[2] = "0.25 to 0.5";
classifyRenderer.Break[3] = 0.5;
classifyRenderer.Label[3] = "0.5 to 0.75";
classifyRenderer.Break[4] = 0.75;
classifyRenderer.Label[4] = "0.75 to 1.0";
classifyRenderer.Break[5] = 1;
classifyRenderer.Label[5] = "1.0 to 1.5";
classifyRenderer.Break[6] = 1.5;
classifyRenderer.Label[6] = "1.5 to 2.0";
classifyRenderer.Break[7] = 2;
classifyRenderer.Label[7] = "2.0 to 3.0";
classifyRenderer.Break[8] = 3;
classifyRenderer.Label[8] = "3.0 to 5.0";
classifyRenderer.Break[9] = 5;
classifyRenderer.Label[9] = "> 5.0";
List <RgbColor> lColors = new List <RgbColor>();
lColors.Add(CreateRGBColor(255, 255, 179));
lColors.Add(CreateRGBColor(252, 241, 167));
lColors.Add(CreateRGBColor(252, 230, 157));
lColors.Add(CreateRGBColor(250, 218, 145));
lColors.Add(CreateRGBColor(250, 208, 135));
lColors.Add(CreateRGBColor(237, 191, 126));
lColors.Add(CreateRGBColor(219, 167, 118));
lColors.Add(CreateRGBColor(201, 147, 111));
lColors.Add(CreateRGBColor(184, 127, 106));
lColors.Add(CreateRGBColor(166, 101, 101));
IFillSymbol fillSymbol = new SimpleFillSymbol();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
return(rasterRenderer);
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return(null);
}
}
public static IRasterRenderer CreateSlopeDegreesColorRamp(GCDConsoleLib.Raster gRaster)
{
IRasterRenderer rasterRenderer = null;
try
{
// Open file raster dataset and ensure that statistics and histograms are present (absence of histograms will cause Renderer.Update() to crash)
IRasterDataset rasterDataset = ArcMapUtilities.GetRasterDataset(gRaster);
ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRastBands = (ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection)rasterDataset;
ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand enumRasterBand = (ESRI.ArcGIS.DataSourcesRaster.IEnumRasterBand)pRastBands.Bands;
rasterDataset.PrecalculateStats(0);
ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRastBand = enumRasterBand.Next();
pRastBand.ComputeStatsAndHist();
// Create the classify renderer
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
rasterRenderer = (IRasterRenderer)classifyRenderer;
IRaster raster = rasterDataset.CreateDefaultRaster();
rasterRenderer.Raster = raster;
classifyRenderer.ClassCount = 10;
rasterRenderer.Update();
classifyRenderer.Break[0] = (double)gRaster.GetStatistics()["min"];
classifyRenderer.Label[0] = "0 to 2";
classifyRenderer.Break[1] = 2;
classifyRenderer.Label[1] = "2 to 5";
classifyRenderer.Break[2] = 5;
classifyRenderer.Label[2] = "5 to 10";
classifyRenderer.Break[3] = 10;
classifyRenderer.Label[3] = "10 to 15";
classifyRenderer.Break[4] = 15;
classifyRenderer.Label[4] = "15 to 25";
classifyRenderer.Break[5] = 25;
classifyRenderer.Label[5] = "25 to 35";
classifyRenderer.Break[6] = 35;
classifyRenderer.Label[6] = "35 to 45";
classifyRenderer.Break[7] = 45;
classifyRenderer.Label[7] = "45 to 60";
classifyRenderer.Break[8] = 60;
classifyRenderer.Label[8] = "60 to 80";
classifyRenderer.Break[9] = 80;
classifyRenderer.Label[9] = "80 to 90";
classifyRenderer.Break[10] = 90;
//Create the symbol for the classes.
List <RgbColor> lColors = CreateSlopeColorRamp();
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = lColors[i];
classifyRenderer.Symbol[i] = (ISymbol)fillSymbol;
}
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
rasterRenderer = null;
}
return(rasterRenderer);
}
public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
try
{
//计算统计图
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IRasterHistogram pRasterHistogram = band.Histogram;
double[] dblValues = pRasterHistogram.Counts as double[];
int intValueCount = dblValues.GetUpperBound(0) + 1;
double[] vValues = new double[intValueCount];
IRasterStatistics pRasterStatistic = band.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
for (int i = 0; i < intValueCount; i++)
{
vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
}
long[] longvalues = new long[dblValues.Length];
for (int i = 0; i <= dblValues.Length - 1; i++)
{
longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
}
//IClassifyGEN classify = null;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
classify = new EqualIntervalClass();
break;
case "自然断点分级":
classify = new NaturalBreaksClass();
break;
}
classify.Classify(ref numClasses);
double[] Classes = classify.ClassBreaks as double[];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassCount = count;
IColor pColor;
ISimpleFillSymbol pSym;
for (int j = 0; j < count; j++)
{
pColor = ramp.get_Color(j * (ramp.Size - 1) / (count - 1));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasRender.Update();
rasClassifyRender.set_Break(j, rasClassifyRender.get_Break(j));
//rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
rasRender.Update();
}
//IRasterProps rasterProps = (IRasterProps)rastlayer.Raster;
//rasterProps.NoDataValue = 0;
//IRasterDisplayProps props = rasClassifyRender as IRasterDisplayProps;
//props.NoDataColor = GET(255, 255, 255);
//rasRender.Update();
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
catch
{
MessageBox.Show("唯一值数量已达到限制（65536）");
}
}
/// <summary>
/// ������Ⱦ
/// </summary>
/// <param name="rasterLayer">դ��ͼ��</param>
/// <param name="VlaueFieldName">ֵ�ֶ�</param>
/// <param name="NormalLizeFieldName">��һ���ֶΣ���ѡ��</param>
/// <param name="classesCount">��������</param>
/// <param name="fromColor">ɫ��ʼ��ɫ</param>
/// <param name="toColor">ɫ������ɫ</param>
/// <param name="bShowClassBreaks">Whether or Not Show class breaks using cell values</param>
/// <param name="NodataColor">��ֵ��ɫ</param>
/// <param name="bUserHillShade">Whether or Not Use hillshade effect</param>
/// <param name="zScale">IF Use hillshade effect Set the Z value</param>
/// <returns></returns>
private void ClassifyRenderer(IRasterLayer rasterLayer, string VlaueFieldName, string NormalLizeFieldName,
int classesCount,IList<Color> lstColor,IList<double> lstBreakVal,IList<string> lstBreakLbl/*IColor fromColor, IColor toColor*/, bool bShowClassBreaks, IColor nodataColor, bool bUserHillShade,Double zScale)
{
try
{
//create the classify renderer
IRasterClassifyColorRampRenderer classifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = (IRasterRenderer)classifyRenderer;
//setup renderer properties
IRaster raster = rasterLayer.Raster;
rasterRenderer.Raster = raster;
classifyRenderer.ClassCount = classesCount;
rasterRenderer.Update();
//create symbol for the classes
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < classifyRenderer.ClassCount; i++)
{
fillSymbol.Color = Convert2ESRIColor(lstColor[i]);
classifyRenderer.set_Break(i, lstBreakVal[i]);
classifyRenderer.set_Symbol(i, (ISymbol)fillSymbol);
classifyRenderer.set_Label(i, lstBreakLbl[i]);
}
if (VlaueFieldName.Length > 0 && VlaueFieldName.ToUpper() != "NONE")
classifyRenderer.ClassField = VlaueFieldName;
if (NormalLizeFieldName.Length > 0 && NormalLizeFieldName.ToUpper()!= "NONE")
classifyRenderer.NormField = NormalLizeFieldName;
(classifyRenderer as IHillShadeInfo).UseHillShade = bUserHillShade;
if(bUserHillShade==true)
{
(classifyRenderer as IHillShadeInfo).ZScale = zScale;
}
if(nodataColor!=null)
(classifyRenderer as IRasterDisplayProps).NoDataColor = nodataColor;
rasterLayer.Renderer = rasterRenderer;
rasterLayer.Draw(esriDrawPhase.esriDPGeography, m_MainMap.ActiveView.ScreenDisplay, null);
m_MainMap.ActiveView.Refresh();
m_MainTOC.Update();
}
catch (Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return ;
}
}
/// <summary>
/// 得到Renderer从数据库中
/// </summary>
/// <returns></returns>
public static object GetRendererFromBlob(string sql, IWorkspace pConfigWks)
{
Exception err = null;
SysCommon.Gis.IGisTable pGISTable = new SysCommon.Gis.SysGisTable(pConfigWks);
//获取Renderer
byte[] renderByte = pGISTable.GetFieldValue("Render", "Render", sql, out err) as byte[];
//是否得到了序列化的blob，没得到，则返回NULL
if (renderByte == null)
{
return(null);
}
IMemoryBlobStreamVariant pMemoryBlobStreamVariant = new MemoryBlobStreamClass();
pMemoryBlobStreamVariant.ImportFromVariant((object)renderByte);
IStream pStream = pMemoryBlobStreamVariant as IStream;
//从数据库中得到符号的类型 shduan 20110721 增加RasterRenderer
string strRenderType = pGISTable.GetFieldValue("Render", "RenderType", sql, out err).ToString();
string strLyrType = pGISTable.GetFieldValue("Render", "LayerType", sql, out err).ToString();
if (strLyrType == "RasterLayer")
{
IRasterRenderer pRRenderer = null;
switch (strRenderType)
{
case "RasterClassifyColorRampRenderer":
pRRenderer = new RasterClassifyColorRampRendererClass();
break;
case "RasterUniqueValueRenderer":
pRRenderer = new RasterUniqueValueRendererClass();
break;
case "RasterDiscreteColorRenderer":
pRRenderer = new RasterDiscreteColorRendererClass();
break;
case "RasterRGBRenderer":
pRRenderer = new RasterRGBRendererClass();
break;
case "RasterStretchColorRampRenderer":
pRRenderer = new RasterStretchColorRampRendererClass();
break;
}
IPersistStream pPersistStream = pRRenderer as IPersistStream;
pPersistStream.Load(pStream);
//pRRenderer = pPersistStream as IRasterRenderer;
return(pRRenderer);
}
else
{
IFeatureRenderer pFRenderer = null;
switch (strRenderType)
{
case "AnnoColor":
string sAnnoColor = Encoding.Default.GetString(renderByte);
return((object)sAnnoColor);
case "SimpleRenderer":
pFRenderer = new SimpleRendererClass();
break;
case "UniqueValueRenderer":
pFRenderer = new UniqueValueRendererClass();
break;
case "ClassBreaksRenderer":
pFRenderer = new ClassBreaksRendererClass();
break;
case "ProportionalSymbolRenderer":
pFRenderer = new ProportionalSymbolRendererClass();
break;
case "ChartRenderer":
pFRenderer = new ChartRendererClass();
break;
}
IPersistStream pPersistStream = pFRenderer as IPersistStream;
pPersistStream.Load(pStream);
pFRenderer = pPersistStream as IFeatureRenderer;
return(pFRenderer);
}
}
private void button_IDW_TIME_Click(object sender, EventArgs e)
{
object obj = null;
ILayer pLayer = new FeatureLayerClass();
IFeatureClass featureClass = GetFeatureClass(@"D:\a_gis工程设计实践课\china\墓穴地shp\grave.shp");
IGeoDataset geo = featureClass as IGeoDataset;
object extend = geo.Extent;
object o = null;
IFeatureClassDescriptor feades = new FeatureClassDescriptorClass();
feades.Create(featureClass, null, "time_inter");
IRasterRadius rasterrad = new RasterRadiusClass();
rasterrad.SetVariable(12, ref obj);
object dCell = 0.014800000;//可以根据不同的点图层进行设置
IInterpolationOp interpla = new RasterInterpolationOpClass();
IRasterAnalysisEnvironment rasanaenv = interpla as IRasterAnalysisEnvironment;
rasanaenv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCell);
rasanaenv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extend, ref o);
IGeoDataset g_GeoDS_Raster = interpla.IDW((IGeoDataset)feades, 2, rasterrad, ref obj);
IRaster pOutRsater = (IRaster)g_GeoDS_Raster;
IRasterLayer pOutRasterlayer = new RasterLayerClass();
pOutRasterlayer.CreateFromRaster(pOutRsater);
pOutRasterlayer.Name = "时间久远程度";
PixelValue2 = Class_GetPixelValue(pOutRasterlayer);
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRRend = pRClassRend as IRasterRenderer;
IRaster pRaster = pOutRasterlayer.Raster;
IRasterBandCollection pRBandCol = pRaster as IRasterBandCollection;
IRasterBand pRBand = pRBandCol.Item(0);
if (pRBand.Histogram == null)
{
pRBand.ComputeStatsAndHist();
}
pRRend.Raster = pRaster;
pRClassRend.ClassCount = 10;
pRRend.Update();
IRgbColor pFromColor = new RgbColorClass();
pFromColor.Red = 135;//天蓝色
pFromColor.Green = 206;
pFromColor.Blue = 235;
IRgbColor pToColor = new RgbColorClass();
pToColor.Red = 124;//草坪绿
pToColor.Green = 252;
pToColor.Blue = 0;
IAlgorithmicColorRamp colorRamp = new AlgorithmicColorRampClass();
colorRamp.Size = 10;
colorRamp.FromColor = pFromColor;
colorRamp.ToColor = pToColor;
bool createColorRamp;
colorRamp.CreateRamp(out createColorRamp);
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < pRClassRend.ClassCount; i++)
{
fillSymbol.Color = colorRamp.get_Color(i);
pRClassRend.set_Symbol(i, fillSymbol as ISymbol);
pRClassRend.set_Label(i, pRClassRend.get_Break(i).ToString("0.00"));
}
pOutRasterlayer.Renderer = pRRend;
this.axMapControl1.AddLayer(pOutRasterlayer);
//两种方法的结合后
}
private void btn_2e_Click(object sender, EventArgs e)
{
ILayer pLayer = new FeatureLayerClass();
IFeatureClass featureClass = GetFeatureClass(@"D:\a_gis工程设计实践课\china\墓穴地shp\grave.shp");
IGeoDataset geo = featureClass as IGeoDataset;
object extend = geo.Extent;
object o = null;
IFeatureClassDescriptor feades = new FeatureClassDescriptorClass();
feades.Create(featureClass, null, "area");
IRasterRadius rasterrad = new RasterRadiusClass();
rasterrad.SetVariable(12, ref o);
object dCell = 0.014800000;//可以根据不同的点图层进行设置
IInterpolationOp interpla = new RasterInterpolationOpClass();
IRasterAnalysisEnvironment rasanaenv = interpla as IRasterAnalysisEnvironment;
rasanaenv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCell);
rasanaenv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extend, ref o);
IGeoDataset g_GeoDS_Raster = interpla.IDW((IGeoDataset)feades, 2, rasterrad, ref o);
IRaster pOutRsater = (IRaster)g_GeoDS_Raster;
IRasterLayer pOutRasterlayer = new RasterLayerClass();
pOutRasterlayer.CreateFromRaster(pOutRsater);
pOutRasterlayer.Name = "两个因素都考虑";
//IRaster raster = pOutRasterlayer.Raster;
IRaster2 raster2 = (IRaster2)pOutRsater;
// 指定像素块大小
IPnt pntBlock = new PntClass();
pntBlock.X = 1280;
pntBlock.Y = 1280;
//创建一个光标以给定像素块大小
//定义RasterCursor初始化，参数设为null，内部自动设置PixelBlock大小
//IRasterCursor pRasterCursor = pRaster2.CreateCursorEx(null);
IRasterCursor rasterCursor = raster2.CreateCursorEx(null);
//控制像素块级别的编辑操作
IRasterEdit rasterEdit = raster2 as IRasterEdit;
if (rasterEdit.CanEdit())
{
//得到一段光栅带
IRasterBandCollection bandCollection = (IRasterBandCollection)pOutRsater;
System.Array pixels;
IPnt pnt = null;
object value;
int bandCount = bandCollection.Count;
//创建像素块
IPixelBlock3 pixelBlock3 = null;
int blockWidth = 0;
int blockHeight = 0;
double temp = 0;
do
{
pixelBlock3 = rasterCursor.PixelBlock as IPixelBlock3;
blockWidth = pixelBlock3.Width;
blockHeight = pixelBlock3.Height;
for (int k = 0; k < bandCount; k++)
{
//指定平面的像素的数组
pixels = (System.Array)pixelBlock3.get_PixelData(k);
for (int i = 0; i < blockWidth; i++)
{
for (int j = 0; j < blockHeight; j++)
{
//value = pixels.GetValue(i, j);
////if (Convert.ToInt32(value) == 0)
//{
//设置像素的颜色值
temp = PixelValue1[i, j] + PixelValue2[i, j];
pixels.SetValue((int)temp, i, j);
//}
}
}
pixelBlock3.set_PixelData(k, pixels);
}
pnt = rasterCursor.TopLeft;
rasterEdit.Write(pnt, (IPixelBlock)pixelBlock3);
}while (rasterCursor.Next());
System.Runtime.InteropServices.Marshal.ReleaseComObject(rasterEdit);
//渲染
IRasterClassifyColorRampRenderer pRClassRend = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRRend = pRClassRend as IRasterRenderer;
IRaster pRaster = pOutRasterlayer.Raster;
IRasterBandCollection pRBandCol = pRaster as IRasterBandCollection;
IRasterBand pRBand = pRBandCol.Item(0);
if (pRBand.Histogram == null)
{
pRBand.ComputeStatsAndHist();
}
pRRend.Raster = pRaster;
pRClassRend.ClassCount = 10;
pRRend.Update();
IRgbColor pFromColor = new RgbColorClass();
pFromColor.Red = 135;//天蓝色
pFromColor.Green = 206;
pFromColor.Blue = 235;
IRgbColor pToColor = new RgbColorClass();
pToColor.Red = 124;//草坪绿
pToColor.Green = 252;
pToColor.Blue = 0;
IAlgorithmicColorRamp colorRamp = new AlgorithmicColorRampClass();
colorRamp.Size = 10;
colorRamp.FromColor = pFromColor;
colorRamp.ToColor = pToColor;
bool createColorRamp;
colorRamp.CreateRamp(out createColorRamp);
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < pRClassRend.ClassCount; i++)
{
fillSymbol.Color = colorRamp.get_Color(i);
pRClassRend.set_Symbol(i, fillSymbol as ISymbol);
pRClassRend.set_Label(i, pRClassRend.get_Break(i).ToString("0.00"));
}
pOutRasterlayer.Renderer = pRRend;
this.axMapControl1.AddLayer(pOutRasterlayer);
}
}
public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IRasterHistogram pRasterHistogram = band.Histogram;
double[] dblValues = pRasterHistogram.Counts as double[];
int intValueCount = dblValues.GetUpperBound(0) + 1;
double[] vValues = new double[intValueCount];
IRasterStatistics pRasterStatistic = band.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
if (dMinValue == 0)
{
pRasterStatistic.IgnoredValues = pRasterStatistic.Minimum;
pRasterStatistic.Recalculate();
dMinValue = pRasterStatistic.Minimum;
}
double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
for (int i = 0; i < intValueCount; i++)
{
vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
}
long[] longvalues = new long[dblValues.Length];
for (int i = 0; i <= dblValues.Length - 1; i++)
{
longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
}
//IClassifyGEN classify = null;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
EqualInterval eqclassify = new EqualIntervalClass();
eqclassify.Classify(vValues, longvalues, ref numClasses);
classify = eqclassify as IClassify;
break;
case "自然断点分级":
NaturalBreaks naclassify = new NaturalBreaksClass();
naclassify.Classify(vValues, longvalues, ref numClasses);
classify = naclassify as IClassify;
break;
}
//switch (classMethod)
//{
// case "等间距分级":
// classify = new EqualIntervalClass();
// break;
// case "自然断点分级":
// classify = new NaturalBreaksClass();
// break;
//}
//classify.Classify(vValues, longvalues, ref numClasses);
double[] Classes = classify.ClassBreaks as double[];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassCount = count;
IColor pColor;
ISimpleFillSymbol pSym;
for (int j = 0; j < count; j++)
{
IRasterProps rasterProps = (IRasterProps)rastlayer.Raster;
rasterProps.NoDataValue = 0;
pColor = ramp.get_Color(j * (ramp.Size / count));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasRender.Update();
rasClassifyRender.set_Break(j, rasClassifyRender.get_Break(j));
rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
rasRender.Update();
}
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
private void btRender_Click(object sender, EventArgs e)
{
int number = 0;
if (cbColorRamp.SelectedIndex == -1)
{
if (lang == Language.Chinese)
{
MessageBox.Show("请先选择一种色带");
}
else if (lang == Language.English)
{
MessageBox.Show("Please choose a color ribbon first");
}
return;
}
try
{
number = int.Parse(txtNumber.Text);
if (number <= 1)
{
if (lang == Language.Chinese)
{
MessageBox.Show("类别数量只能为>1的正整数");
}
else if (lang == Language.English)
{
MessageBox.Show("The class number can only be a positive integer of >1");
}
}
}
catch
{
if (lang == Language.Chinese)
{
MessageBox.Show("类别数量只能为整数");
}
else if (lang == Language.English)
{
MessageBox.Show("The class number can only be an integer");
}
}
IRasterClassifyColorRampRenderer pRasterClassifyRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasterRenderer = pRasterClassifyRenderer as IRasterRenderer;
IUniqueValues uniVal = new UniqueValuesClass();
IRasterCalcUniqueValues2 calValues = new RasterCalcUniqueValuesClass();
calValues.MaxUniqueValueCount = 10000000;
calValues.AddFromRaster(rasterLayer.Raster, cbSelectedBand.SelectedIndex, uniVal);
object vValues, vFrequences;
uniVal.GetHistogram(out vValues, out vFrequences);
double[] doubleArray = (double[])vValues;
int[] longArray = (int[])vFrequences;
int order = SortTool.CheckInOrder(doubleArray, doubleArray.Length);
if (order == 0)
{
SortTool.QuickSort3(doubleArray, 0, doubleArray.Length - 1, longArray);
}
else if (order == -1)
{
SortTool.Reverse(doubleArray, doubleArray.Length, longArray);
}
IClassifyGEN classify = GetMethod(cbMethod.Text);
classify.Classify((double[])vValues, (int[])vFrequences, ref number);
double[] classes = classify.ClassBreaks as double[];
pRasterRenderer.Raster = rasterLayer.Raster;
pRasterClassifyRenderer.ClassCount = number;
for (int i = 0; i < classes.Length; i++)
{
pRasterClassifyRenderer.set_Break(i, classes[i]);
}
pRasterRenderer.Update();
IStyleGalleryItem galleryItem = symbol.GetItem(cbColorRamp.SelectedIndex);
IColorRamp colorRamp = galleryItem.Item as IColorRamp;
IFillSymbol fill = new SimpleFillSymbolClass() as IFillSymbol;
int increase = colorRamp.Size / (number - 1);
for (int i = 0; i < pRasterClassifyRenderer.ClassCount - 1; i++)
{
fill.Color = colorRamp.get_Color(i * increase);
pRasterClassifyRenderer.set_Symbol(i, fill as ISymbol);
pRasterClassifyRenderer.set_Label(i, pRasterClassifyRenderer.get_Break(i).ToString() + "-" + pRasterClassifyRenderer.get_Break(i + 1).ToString());
}
fill.Color = colorRamp.get_Color(colorRamp.Size - 1);
pRasterClassifyRenderer.set_Symbol(pRasterClassifyRenderer.ClassCount - 1, fill as ISymbol);
pRasterClassifyRenderer.set_Label(pRasterClassifyRenderer.ClassCount - 1, pRasterClassifyRenderer.get_Break(pRasterClassifyRenderer.ClassCount - 1).ToString() + "-" + pRasterClassifyRenderer.get_Break(pRasterClassifyRenderer.ClassCount).ToString());
rasterLayer.Renderer = pRasterRenderer;
succeed = true;
this.Close();
}
public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IRasterHistogram pRasterHistogram = band.Histogram;
double[] dblValues = pRasterHistogram.Counts as double[];
int intValueCount = dblValues.GetUpperBound(0) + 1;
double[] vValues = new double[intValueCount];
IRasterStatistics pRasterStatistic = band.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
if (dMinValue == 0)
{
pRasterStatistic.IgnoredValues = pRasterStatistic.Minimum;
pRasterStatistic.Recalculate();
dMinValue = pRasterStatistic.Minimum;
}
double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
for (int i = 0; i < intValueCount; i++)
{
vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
}
long[] longvalues = new long[dblValues.Length];
for (int i = 0; i <= dblValues.Length - 1; i++)
{
longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
}
//IClassifyGEN classify = null;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
EqualInterval eqclassify = new EqualIntervalClass();
eqclassify.Classify(vValues, longvalues, ref numClasses);
classify = eqclassify as IClassify;
break;
case "自然断点分级":
NaturalBreaks naclassify = new NaturalBreaksClass();
naclassify.Classify(vValues, longvalues, ref numClasses);
classify = naclassify as IClassify;
break;
}
//switch (classMethod)
//{
// case "等间距分级":
// classify = new EqualIntervalClass();
// break;
// case "自然断点分级":
// classify = new NaturalBreaksClass();
// break;
//}
//classify.Classify(vValues, longvalues, ref numClasses);
double[] Classes = classify.ClassBreaks as double[];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassCount = count;
IColor pColor;
ISimpleFillSymbol pSym;
for (int j = 0; j < count; j++)
{
IRasterProps rasterProps = (IRasterProps)rastlayer.Raster;
rasterProps.NoDataValue = 0;
pColor = ramp.get_Color(j * (ramp.Size / count));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasRender.Update();
rasClassifyRender.set_Break(j,rasClassifyRender.get_Break(j));
rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
rasRender.Update();
}
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
//public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp,double minValue)
//{
// IRasterBand band = GetBand(rastlayer);
// if (band.Histogram == null)
// {
// band.ComputeStatsAndHist();
// }
// IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
// IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
// rasRender.Raster = rastlayer.Raster;
// rasRender.Update();
// int numClasses = count;
// IRasterHistogram pRasterHistogram = band.Histogram;
// double[] dblValues = pRasterHistogram.Counts as double[];
// int intValueCount = dblValues.GetUpperBound(0) + 1;
// double[] vValues = new double[intValueCount];
// IRasterStatistics pRasterStatistic = band.Statistics;
// //double dMaxValue = pRasterStatistic.Maximum;
// double dMaxValue = minValue;
// double dMinValue = pRasterStatistic.Minimum;
// double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
// for (int i = 0; i < intValueCount; i++)
// {
// vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
// }
// long[] longvalues = new long[dblValues.Length];
// for (int i = 0; i <= dblValues.Length - 1; i++)
// {
// longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
// }
// //IClassifyGEN classify = null;
// IClassify classify = null;
// switch (classMethod)
// {
// case "等间距分级":
// EqualInterval eqclassify = new EqualIntervalClass();
// eqclassify.Classify(vValues, longvalues, ref numClasses);
// classify = eqclassify as IClassify;
// break;
// case "自然断点分级":
// NaturalBreaks naclassify = new NaturalBreaksClass();
// naclassify.Classify(vValues, longvalues, ref numClasses);
// classify = naclassify as IClassify;
// break;
// }
// //switch (classMethod)
// //{
// // case "等间距分级":
// // classify = new EqualIntervalClass();
// // break;
// // case "自然断点分级":
// // classify = new NaturalBreaksClass();
// // break;
// //}
// //classify.Classify(vValues, longvalues, ref numClasses);
// double[] Classes = classify.ClassBreaks as double[];
// UID pUid = classify.ClassID;
// IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
// rasClassifyUI.ClassificationMethod = pUid;
// rasClassifyRender.ClassCount = count;
// IColor pColor;
// ISimpleFillSymbol pSym;
// for (int j = 0; j < count; j++)
// {
// pColor = ramp.get_Color(j * (ramp.Size / count));
// pSym = new SimpleFillSymbolClass();
// pSym.Color = pColor;
// rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
// rasRender.Update();
// if (Classes[j] == 0)
// {
// rasClassifyRender.set_Label(j, Classes[j].ToString() + "-" + Classes[j + 1].ToString("0.000"));
// rasRender.Update();
// }
// else
// {
// rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
// rasRender.Update();
// }
// }
// rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
try
{ //进行唯一值判断
IUniqueValues values = new UniqueValuesClass();
IRasterCalcUniqueValues unique = new RasterCalcUniqueValuesClass();
unique.AddFromRaster(rastlayer.Raster, 0, values);
int uniquecount = values.get_UniqueCount(0);
//计算统计图
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IRasterHistogram pRasterHistogram = band.Histogram;
double[] dblValues = pRasterHistogram.Counts as double[];
int intValueCount = dblValues.GetUpperBound(0) + 1;
double[] vValues = new double[intValueCount];
IRasterStatistics pRasterStatistic = band.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
if (dMinValue == 0)
{
pRasterStatistic.IgnoredValues = pRasterStatistic.Minimum;
pRasterStatistic.Recalculate();
dMinValue = pRasterStatistic.Minimum;
}
double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
for (int i = 0; i < intValueCount; i++)
{
vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
}
long[] longvalues = new long[dblValues.Length];
for (int i = 0; i <= dblValues.Length - 1; i++)
{
longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
}
//IClassifyGEN classify = null;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
EqualInterval eqclassify = new EqualIntervalClass();
eqclassify.Classify(vValues, longvalues, ref numClasses);
classify = eqclassify as IClassify;
break;
case "自然断点分级":
NaturalBreaks naclassify = new NaturalBreaksClass();
naclassify.Classify(vValues, longvalues, ref numClasses);
classify = naclassify as IClassify;
break;
}
#region
//switch (classMethod)
//{
// case "等间距分级":
// classify = new EqualIntervalClass();
// break;
// case "自然断点分级":
// classify = new NaturalBreaksClass();
// break;
//}
//classify.Classify(vValues, longvalues, ref numClasses);
#endregion
double[] Classes = classify.ClassBreaks as double[];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassCount = count;
IColor pColor;
ISimpleFillSymbol pSym;
//排除数值
double[] exdouble = new double[2] { 0, Classes[0] };
IRasterDataExclusion ex = rasClassifyRender as IRasterDataExclusion;
ex.ExcludeValues = exdouble;
ex.ExcludeColor = GET(255, 255, 255);
for (int j = 0; j < count; j++)
{
pColor = ramp.get_Color(j * (ramp.Size / count));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasRender.Update();
rasClassifyRender.set_Break(j, rasClassifyRender.get_Break(j));
rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
rasRender.Update();
}
//IRasterProps rasterProps = (IRasterProps)rastlayer.Raster;
//rasterProps.NoDataValue = 0;
//IRasterDisplayProps props = rasClassifyRender as IRasterDisplayProps;
//props.NoDataColor = GET(255, 255, 255);
//rasRender.Update();
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
catch(Exception ex)
{
//MessageBox.Show(ex.ToString());
if (ex.ToString().Contains("唯一值过多"))
{
MessageBox.Show("唯一值数量已达到限制（65536）");
}
else
{
MessageBox.Show("还未计算唯一值，进行唯一值计算！");
}
}
//catch
//{
// MessageBox.Show("数据需要先计算统计值，请点击确定进行计算！");
//}
}
//Ϊդ��ͼ�����ݽ��з�����ɫ
public static IRasterLayer GetRLayerClassifyColor(IRaster pRaster, double dMaxDis)
{
IRasterClassifyColorRampRenderer pClassRen = new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasRen = pClassRen as RasterClassifyColorRampRendererClass;
//Set raster for the render and update
pRasRen.Raster = pRaster;
pClassRen.ClassCount = 10;
pRasRen.Update();
//Create a color ramp to use
//���������յ���ɫ
IColor pFromColor = new RgbColorClass();
IRgbColor pRgbColor = pFromColor as IRgbColor;
pRgbColor.Red = 255;
pRgbColor.Green = 0;
pRgbColor.Blue = 0;
IColor pToColor = new RgbColorClass();
pRgbColor = pToColor as IRgbColor;
pRgbColor.Red = 0;
pRgbColor.Green = 255;
pRgbColor.Blue = 0;
//������ɫ�ּ�
IAlgorithmicColorRamp pRamp = new AlgorithmicColorRampClass();
pRamp.Size = 10;
pRamp.FromColor = pFromColor;
pRamp.ToColor = pToColor;
bool ok = true;
pRamp.CreateRamp(out ok);
//Create symbol for the classes
IFillSymbol pFSymbol = new SimpleFillSymbolClass();
double LabelValue = Convert.ToDouble(dMaxDis / 6);
for (int i = 0; i <= pClassRen.ClassCount - 1; i++)
{
pFSymbol.Color = pRamp.get_Color(i);
pClassRen.set_Symbol(i, pFSymbol as ISymbol);
double h1 = LabelValue * i;
double h2 = LabelValue * (i + 1);
pClassRen.set_Label(i, h1.ToString() + "��" + h2.ToString());
}
//Update the renderer and plug into layer
pRasRen.Update();
IRasterLayer pRLayer = new RasterLayerClass();
pRLayer.CreateFromRaster(pRaster);
pRLayer.Renderer = pClassRen as IRasterRenderer;
return pRLayer;
}
public static IRasterRenderer ClassifyRenderer(ESRI.ArcGIS.Geodatabase.IRasterDataset pRasterDataset)
{
try
{
//create the classify render:创建渲染器
IRasterClassifyColorRampRenderer pClassifyRenderer=new RasterClassifyColorRampRendererClass();
IRasterRenderer pRasterRenderer=(IRasterRenderer)pClassifyRenderer;
IRaster pRaster=pRasterDataset.CreateDefaultRaster();
pRasterRenderer.Raster=pRaster;
pClassifyRenderer.ClassCount=10;
pRasterRenderer.Update();
//创建颜色
IRgbColor pFromColor=new RgbColorClass();
pFromColor.Red=255;
pFromColor.Green=0;
pFromColor.Blue=0;
IRgbColor pToColor=new RgbColorClass();
pToColor.Red=0;
pToColor.Green=255;
pToColor.Blue=255;
//set the color ramp for the symbology:生成色带
IAlgorithmicColorRamp pRamp=new AlgorithmicColorRampClass();
pRamp.Size=10;
pRamp.FromColor=pFromColor;
pRamp.ToColor=pToColor;
bool pBoolColorRamp;
pRamp.CreateRamp(out pBoolColorRamp);
//create the symbol for the classes.:创建符号
IFillSymbol pFillSymbol=new SimpleFillSymbolClass();
for (int i=0;i<pClassifyRenderer.ClassCount;i++){
pFillSymbol.Color=pRamp.get_Color(i);
pClassifyRenderer.set_Symbol(i,(ISymbol)pFillSymbol);
pClassifyRenderer.set_Label(i,Convert.ToString(i));
}
return pRasterRenderer;
}
catch (System.Exception ex)
{
System.Diagnostics.Debug.WriteLine(ex.Message);
return null;
}
}
//public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp,double minValue)
//{
// IRasterBand band = GetBand(rastlayer);
// if (band.Histogram == null)
// {
// band.ComputeStatsAndHist();
// }
// IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
// IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
// rasRender.Raster = rastlayer.Raster;
// rasRender.Update();
// int numClasses = count;
// IRasterHistogram pRasterHistogram = band.Histogram;
// double[] dblValues = pRasterHistogram.Counts as double[];
// int intValueCount = dblValues.GetUpperBound(0) + 1;
// double[] vValues = new double[intValueCount];
// IRasterStatistics pRasterStatistic = band.Statistics;
// //double dMaxValue = pRasterStatistic.Maximum;
// double dMaxValue = minValue;
// double dMinValue = pRasterStatistic.Minimum;
// double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
// for (int i = 0; i < intValueCount; i++)
// {
// vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
// }
// long[] longvalues = new long[dblValues.Length];
// for (int i = 0; i <= dblValues.Length - 1; i++)
// {
// longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
// }
// //IClassifyGEN classify = null;
// IClassify classify = null;
// switch (classMethod)
// {
// case "等间距分级":
// EqualInterval eqclassify = new EqualIntervalClass();
// eqclassify.Classify(vValues, longvalues, ref numClasses);
// classify = eqclassify as IClassify;
// break;
// case "自然断点分级":
// NaturalBreaks naclassify = new NaturalBreaksClass();
// naclassify.Classify(vValues, longvalues, ref numClasses);
// classify = naclassify as IClassify;
// break;
// }
// //switch (classMethod)
// //{
// // case "等间距分级":
// // classify = new EqualIntervalClass();
// // break;
// // case "自然断点分级":
// // classify = new NaturalBreaksClass();
// // break;
// //}
// //classify.Classify(vValues, longvalues, ref numClasses);
// double[] Classes = classify.ClassBreaks as double[];
// UID pUid = classify.ClassID;
// IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
// rasClassifyUI.ClassificationMethod = pUid;
// rasClassifyRender.ClassCount = count;
// IColor pColor;
// ISimpleFillSymbol pSym;
// for (int j = 0; j < count; j++)
// {
// pColor = ramp.get_Color(j * (ramp.Size / count));
// pSym = new SimpleFillSymbolClass();
// pSym.Color = pColor;
// rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
// rasRender.Update();
// if (Classes[j] == 0)
// {
// rasClassifyRender.set_Label(j, Classes[j].ToString() + "-" + Classes[j + 1].ToString("0.000"));
// rasRender.Update();
// }
// else
// {
// rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
// rasRender.Update();
// }
// }
// rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
public void RasterClassify(IRasterLayer rastlayer, string classMethod, int count, IColorRamp ramp)
{
try
{ //进行唯一值判断
IUniqueValues values = new UniqueValuesClass();
IRasterCalcUniqueValues unique = new RasterCalcUniqueValuesClass();
unique.AddFromRaster(rastlayer.Raster, 0, values);
int uniquecount = values.get_UniqueCount(0);
//计算统计图
IRasterBand band = GetBand(rastlayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRasterClassifyColorRampRenderer rasClassifyRender = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasRender = rasClassifyRender as IRasterRenderer;
rasRender.Raster = rastlayer.Raster;
rasRender.Update();
int numClasses = count;
IRasterHistogram pRasterHistogram = band.Histogram;
double[] dblValues = pRasterHistogram.Counts as double[];
int intValueCount = dblValues.GetUpperBound(0) + 1;
double[] vValues = new double[intValueCount];
IRasterStatistics pRasterStatistic = band.Statistics;
double dMaxValue = pRasterStatistic.Maximum;
double dMinValue = pRasterStatistic.Minimum;
if (dMinValue == 0)
{
pRasterStatistic.IgnoredValues = pRasterStatistic.Minimum;
pRasterStatistic.Recalculate();
dMinValue = pRasterStatistic.Minimum;
}
double BinInterval = Convert.ToDouble((dMaxValue - dMinValue) / intValueCount);
for (int i = 0; i < intValueCount; i++)
{
vValues[i] = i * BinInterval + pRasterStatistic.Minimum;
}
long[] longvalues = new long[dblValues.Length];
for (int i = 0; i <= dblValues.Length - 1; i++)
{
longvalues[i] = long.Parse(Convert.ToString(dblValues[i]));
}
//IClassifyGEN classify = null;
IClassify classify = null;
switch (classMethod)
{
case "等间距分级":
EqualInterval eqclassify = new EqualIntervalClass();
eqclassify.Classify(vValues, longvalues, ref numClasses);
classify = eqclassify as IClassify;
break;
case "自然断点分级":
NaturalBreaks naclassify = new NaturalBreaksClass();
naclassify.Classify(vValues, longvalues, ref numClasses);
classify = naclassify as IClassify;
break;
}
#region
//switch (classMethod)
//{
// case "等间距分级":
// classify = new EqualIntervalClass();
// break;
// case "自然断点分级":
// classify = new NaturalBreaksClass();
// break;
//}
//classify.Classify(vValues, longvalues, ref numClasses);
#endregion
double[] Classes = classify.ClassBreaks as double[];
int n = Classes.Count();
double dn = Classes[0];
UID pUid = classify.ClassID;
IRasterClassifyUIProperties rasClassifyUI = rasClassifyRender as IRasterClassifyUIProperties;
rasClassifyUI.ClassificationMethod = pUid;
rasClassifyRender.ClassCount = count;
IColor pColor;
ISimpleFillSymbol pSym;
//排除数值
double[] exdouble = new double[2] {
0, Classes[0]
};
IRasterDataExclusion ex = rasClassifyRender as IRasterDataExclusion;
ex.ExcludeValues = exdouble;
ex.ExcludeColor = GET(255, 255, 255);
for (int j = 0; j < count; j++)
{
pColor = ramp.get_Color(j * (ramp.Size / count));
pSym = new SimpleFillSymbolClass();
pSym.Color = pColor;
rasClassifyRender.set_Symbol(j, (ISymbol)pSym);
rasRender.Update();
rasClassifyRender.set_Break(j, rasClassifyRender.get_Break(j));
rasClassifyRender.set_Label(j, Classes[j].ToString("0.000") + "-" + Classes[j + 1].ToString("0.000"));
rasRender.Update();
}
//IRasterProps rasterProps = (IRasterProps)rastlayer.Raster;
//rasterProps.NoDataValue = 0;
//IRasterDisplayProps props = rasClassifyRender as IRasterDisplayProps;
//props.NoDataColor = GET(255, 255, 255);
//rasRender.Update();
rastlayer.Renderer = rasClassifyRender as IRasterRenderer;
}
catch (Exception ex)
{
//MessageBox.Show(ex.ToString());
if (ex.ToString().Contains("唯一值过多"))
{
MessageBox.Show("唯一值数量已达到限制（65536）");
}
else
{
MessageBox.Show("还未计算唯一值，进行唯一值计算！");
}
}
//catch
//{
// MessageBox.Show("数据需要先计算统计值，请点击确定进行计算！");
//}
}
private void ChangeRender2ClassifyColorRampRenderer(IRasterLayer rasterLayer, int classCount, IColorRamp pColorRamp)
{
if (rasterlayer == null)
{
return;
}
//计算统计图
IRasterBand band = GetBand(rasterLayer);
if (band.Histogram == null)
{
band.ComputeStatsAndHist();
}
IRaster raster = rasterLayer.Raster;
//创建分类
IRasterClassifyColorRampRenderer rasterClassifyColorRampRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRenderer = rasterClassifyColorRampRenderer as IRasterRenderer;
rasterRenderer.Raster = raster;
rasterRenderer.Update();
//IClassify classify = null;
//switch (classMethod)
//{
// case "等间距分级":
// classify = new EqualIntervalClass();
// break;
// case "自然断点分级":
// classify = new NaturalBreaksClass();
// break;
//}
//int numClass = classCount;
//classify.Classify(ref numClass);
//double[] Classes = classify.ClassBreaks as double[];
//UID pUid = classify.ClassID;
//IRasterClassifyUIProperties rasClassifyUI = rasterClassifyColorRampRenderer as IRasterClassifyUIProperties;
//rasClassifyUI.ClassificationMethod = pUid;
rasterClassifyColorRampRenderer.ClassCount = classCount;
//IAlgorithmicColorRamp algorithmicColorRamp = new AlgorithmicColorRampClass();
//algorithmicColorRamp.Size = 20;
//algorithmicColorRamp.FromColor = getRGBColor(255, 0, 0);
//algorithmicColorRamp.ToColor = getRGBColor(0, 255, 0);
//bool btrue;
//algorithmicColorRamp.CreateRamp(out btrue);
//IEnumColors enumColors = pColorRamp.Colors;
IEnumColors enumColors = pColorRamp.Colors;
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
for (int i = 0; i < rasterClassifyColorRampRenderer.ClassCount - 1; i++)
{
fillSymbol.Color = enumColors.Next();
rasterClassifyColorRampRenderer.set_Symbol(i, fillSymbol as ISymbol);
//rasterClassifyColorRampRenderer.set_Label(i, "Class" + i.ToString());
}
//ISimpleFillSymbol fillSymbol = new SimpleFillSymbolClass();
//IColor pColor;
//for (int i = 0; i < classCount - 1; i++)
//{
// pColor = pColorRamp.get_Color(i * (pColorRamp.Size / classCount));
// fillSymbol = new SimpleFillSymbolClass();
// fillSymbol.Color = pColor;
// rasterClassifyColorRampRenderer.set_Symbol(i, fillSymbol as ISymbol);
// //rasterRenderer.Update();
// //rasterClassifyColorRampRenderer.set_Break(i, rasterClassifyColorRampRenderer.get_Break(i));
// //rasterClassifyColorRampRenderer.set_Label(i, "Class" + i.ToString());
// rasterRenderer.Update();
//}
rasterRenderer.Update();
rasterLayer.Renderer = rasterClassifyColorRampRenderer as IRasterRenderer;
}
