private static void DirectRender(ILayer pLayer)
{
IGeoFeatureLayer MyGeoFeatureLayer = pLayer as IGeoFeatureLayer;
IFeatureCursor MyFeatureCursor = MyGeoFeatureLayer.Search(null, false);
IFeature MyFeature = MyFeatureCursor.NextFeature();
Render renderer = new Render();
IUniqueValueRenderer pUniqueValueR = new UniqueValueRendererClass();
pUniqueValueR.FieldCount = 1;//单值渲染
string FieldName = MyFeature.Fields.get_Field(0).Name;
pUniqueValueR.set_Field(0, FieldName);//渲染字段
while (null != MyFeature)
{
//MyFeature.Fields.FindField("OBJECTID");
string value = Convert.ToString((int)MyFeature.get_Value(0));
pUniqueValueR.AddValue(value, FieldName, renderer.GetSymbol(MyFeature));
MyFeature = MyFeatureCursor.NextFeature();
}
//pUniqueValueR.AddValue((Convert.ToString(tem)), "FID", pSymbol);
MyGeoFeatureLayer.Renderer = pUniqueValueR as IFeatureRenderer;
}
public static List <IFeature> GetAllFeaturesFromPolygonInGeoFeatureLayer(IPolygon polygon, IGeoFeatureLayer geoFeatureLayer, IActiveView activeView)
{
List <IFeature> list = new List <IFeature>();
if (polygon == null || geoFeatureLayer == null || activeView == null)
{
return(list);
}
IMap pMap = activeView.FocusMap;
ISpatialFilter pSpatialFilter = new SpatialFilterClass();
pSpatialFilter.Geometry = polygon;
pSpatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;
pSpatialFilter.set_OutputSpatialReference(geoFeatureLayer.FeatureClass.ShapeFieldName, pMap.SpatialReference);
pSpatialFilter.GeometryField = geoFeatureLayer.FeatureClass.ShapeFieldName;
IFeatureCursor pFeatureCursor = geoFeatureLayer.Search(pSpatialFilter, false);
IFeature pFeature;
while ((pFeature = pFeatureCursor.NextFeature()) != null)
{
list.Add(pFeature);
}
Marshal.ReleaseComObject(pFeatureCursor);
return(list);
}
private static Dictionary<string, int> GetAllValues(string fieldName, IGeoFeatureLayer layer)
{
var results = new Dictionary<string, int>();
IQueryFilter query = new QueryFilter { SubFields = fieldName };
using (var comReleaser = new ComReleaser())
{
IFeatureCursor cursor = layer.Search(query, true);
comReleaser.ManageLifetime(cursor);
IFeature feature;
int fieldIndex = cursor.FindField(fieldName);
while ((feature = cursor.NextFeature()) != null)
{
results[feature.Value[fieldIndex].ToString()] = 1;
}
}
return results;
}
/// <summary>
/// Occurs when this command is clicked
/// </summary>
public override void OnClick()
{
try
{
//get access to the document and the active view
pMxDocument = (IMxDocument)clsAgrcArcMapExtension.m_application.Document;
pMap = pMxDocument.FocusMap;
pActiveView = pMxDocument.ActiveView; //pActiveView = (IActiveView)pMap;
//get the map's graphics layer
ICompositeGraphicsLayer2 pComGraphicsLayer = pMap.BasicGraphicsLayer as ICompositeGraphicsLayer2;
ICompositeLayer pCompositeLayer = pComGraphicsLayer as ICompositeLayer;
ILayer pLayer;
//loop through all graphic layers in the map and check for the 'PolyVertices' layer, if found, delete it, in order to start fresh
for (int i = 0; i < pCompositeLayer.Count; i++)
{
pLayer = pCompositeLayer.get_Layer(i);
if (pLayer.Name == "PolyVertices")
{
pComGraphicsLayer.DeleteLayer("PolyVertices");
break;
}
}
if (bolVerticesOn == false)
{
IGraphicsLayer pGraphicsLayer = pComGraphicsLayer.AddLayer("PolyVertices", null);
pMap.ActiveGraphicsLayer = (ILayer)pGraphicsLayer;
IGraphicsContainer pGraphicsContainer = pComGraphicsLayer.FindLayer("PolyVertices") as IGraphicsContainer;
//make sure the user has selected a polygon or polyline layer
if (pMxDocument.SelectedLayer == null)
{
MessageBox.Show("Please select a layer.", "Select Layer", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (!(pMxDocument.SelectedLayer is IFeatureLayer))
{
MessageBox.Show("Please select a polygon or line layer.", "Polygon or Line", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
//cast the selected layer as a feature layer
IGeoFeatureLayer pGFlayer = (IGeoFeatureLayer)pMxDocument.SelectedLayer;
//check if the feaure layer is a polygon or line layer
if (pGFlayer.FeatureClass.ShapeType != esriGeometryType.esriGeometryPolygon & pGFlayer.FeatureClass.ShapeType != esriGeometryType.esriGeometryPolyline)
{
MessageBox.Show("Please select a polygon or line layer.", "Polygon or Line", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
//setup marker symbol
ISimpleMarkerSymbol pSimpleMarker = new SimpleMarkerSymbol();
ISymbol pSymbolMarker = (ISymbol)pSimpleMarker;
IRgbColor pRgbColor = new ESRI.ArcGIS.Display.RgbColorClass();
pRgbColor.Red = 223;
pRgbColor.Green = 155;
pRgbColor.Blue = 255;
pSimpleMarker.Color = pRgbColor;
pSimpleMarker.Style = esriSimpleMarkerStyle.esriSMSDiamond;
pSimpleMarker.Size = 8;
//setup line symbol
ISimpleLineSymbol pSimpleLineSymbol = new SimpleLineSymbol();
ISymbol pSymbolLine = (ISymbol)pSimpleLineSymbol;
pRgbColor = new ESRI.ArcGIS.Display.RgbColor();
pRgbColor.Red = 0;
pRgbColor.Green = 255;
pRgbColor.Blue = 0;
pSimpleLineSymbol.Color = pRgbColor;
pSimpleLineSymbol.Style = esriSimpleLineStyle.esriSLSSolid;
pSimpleLineSymbol.Width = 1;
//setup simplefill symbol
ISimpleFillSymbol pSimpleFillSymbol = new SimpleFillSymbol();
ISymbol pSymbolPolygon = (ISymbol)pSimpleFillSymbol;
pRgbColor = new ESRI.ArcGIS.Display.RgbColor();
pRgbColor.Red = 0;
pRgbColor.Green = 0;
pRgbColor.Blue = 255;
pSimpleFillSymbol.Color = pRgbColor;
pSimpleFillSymbol.Style = esriSimpleFillStyle.esriSFSSolid;
//get all the polygons in the current map extent ina cursor
IEnvelope pMapExtent = pActiveView.Extent;
ISpatialFilter pQFilter = new SpatialFilter();
pQFilter.GeometryField = "SHAPE";
pQFilter.Geometry = pMapExtent;
pQFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;
IFeatureCursor pFCursor = pGFlayer.Search(pQFilter, true);
//draw each polygon and then each polygon's point collection
IFeature pFeature = pFCursor.NextFeature();
IGeometry pGeometry;
while (pFeature != null)
{
pGeometry = pFeature.Shape;
//draw the polygon
//draw each vertex on the polygon
IPointCollection pPointCollection = pGeometry as IPointCollection;
for (int i = 0; i < pPointCollection.PointCount; i++)
{
IGeometry pPtGeom = pPointCollection.get_Point(i);
IElement pElement = new MarkerElement();
pElement.Geometry = pPtGeom;
IMarkerElement pMarkerElement = pElement as IMarkerElement;
pMarkerElement.Symbol = pSimpleMarker;
pGraphicsContainer.AddElement(pElement, 0);
}
pFeature = pFCursor.NextFeature();
}
bolVerticesOn = true;
}
else //if (bolVerticesOn == true)
{
bolVerticesOn = false;
}
//refresh the map
pActiveView.Refresh();
pActiveView.Refresh();
}
catch (Exception ex)
{
MessageBox.Show("Error Message: " + Environment.NewLine + ex.Message + Environment.NewLine + Environment.NewLine +
"Error Source: " + Environment.NewLine + ex.Source + Environment.NewLine + Environment.NewLine +
"Error Location:" + Environment.NewLine + ex.StackTrace,
"AGRC Custom Tools ArcMap Error!", MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
}
}
public override void OnMouseDown(int Button, int Shift, int X, int Y)
{
if (pDrawDenplot.cboSourceLayer.Text != "" && pDrawDenplot.cboUField.Text != "" && pDrawDenplot.cboValueField.Text != "")
{
IGeoFeatureLayer pGeoFeatureLayer = pDrawDenplot.pGeofeatureLayer;
mForm = System.Windows.Forms.Application.OpenForms["MainForm"] as MainForm;
IActiveView pActiveView = mForm.axMapControl1.ActiveView;
string strTargetLayerName = pDrawDenplot.cboSourceLayer.Text;
string strValueField = pDrawDenplot.cboValueField.Text;
string strUncerField = pDrawDenplot.cboUField.Text;
clsSnippet pSnippet = new clsSnippet();
//Using IRubberband
IEnvelope pEnvelop = pSnippet.DrawRectangle(pActiveView);
ISpatialFilter pSpatialFilter = new SpatialFilterClass();
pSpatialFilter.Geometry = pEnvelop;
pSpatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;
int intTLayerIdx = pSnippet.GetIndexNumberFromLayerName(pActiveView, strTargetLayerName);
ILayer pLayer = mForm.axMapControl1.get_Layer(intTLayerIdx);
IFeatureLayer pFLayer = (IFeatureLayer)pLayer;
string ShapeFieldName = pFLayer.FeatureClass.ShapeFieldName;
pSpatialFilter.GeometryField = pFLayer.FeatureClass.ShapeFieldName;
IFeatureClass pFeatureClass = pFLayer.FeatureClass;
int intValueFldIdx = pFeatureClass.FindField(strValueField);
int intUncerFldIdx = pFeatureClass.FindField(strUncerField);
IFeatureCursor pFCursor = pGeoFeatureLayer.Search(pSpatialFilter, true);
StackedChartRenderer pStackedChartRenderer = new StackedChartRenderer();
pStackedChartRenderer.dblError = pDrawDenplot.dblError;
pStackedChartRenderer.intValueFldIdx = intValueFldIdx;
pStackedChartRenderer.intUncerFldIdx = intUncerFldIdx;
pStackedChartRenderer.dblMaxValue = pDrawDenplot.dblMaxValue;
pStackedChartRenderer.bln3Dfeature = pDrawDenplot.bln3Dfeature;
pStackedChartRenderer.m_strOriRenderField = strValueField;
pStackedChartRenderer.m_strUncernRenderField = strUncerField;
pStackedChartRenderer.dblMaxEstimate = pDrawDenplot.dblMaxEstimate;
pStackedChartRenderer.dblMaxUncern = pDrawDenplot.dblMaxUncern;
pStackedChartRenderer.m_pQueryFilter = pSpatialFilter as IQueryFilter;
//pStackedChartRenderer.m_stackedChart = pGeoFeatureLayer.Renderer;
//pStackedChartRenderer.PrepareFilter(pFeatureClass, pSpatialFilter);
pStackedChartRenderer.m_pDisplay = pActiveView.ScreenDisplay;
//pStackedChartRenderer.Draw(pFCursor, esriDrawPhase.esriDPSelection, pDisplay, null);
//pStackedChartRenderer.CreateLegend();
pGeoFeatureLayer.Renderer = pStackedChartRenderer;
//mForm.axMapControl1.ActiveView.Refresh();
pActiveView.Refresh();
}
}
private void btnApply_Click(object sender, EventArgs e)
{
//Choropleth mapping : Exactly same with the function in frmChoroplethwithOverlay HK102915
string strLayerName = cboSourceLayer.Text;
if (cboSourceLayer.Text == "" || cboValueField.Text == "" || cboUField.Text == "")
{
MessageBox.Show("Assign proper layer and field");
return;
}
int intLIndex = pSnippet.GetIndexNumberFromLayerName(pActiveView, strLayerName);
ILayer pLayer = mForm.axMapControl1.get_Layer(intLIndex);
IFeatureLayer pFLayer = pLayer as IFeatureLayer;
IFeatureClass pFClass = pFLayer.FeatureClass;
//Create Rendering of Mean Value at Target Layer
int intGCBreakeCount = Convert.ToInt32(nudGCNClasses.Value);
string strGCRenderField = cboValueField.Text;
IGeoFeatureLayer pGeofeatureLayer;
pGeofeatureLayer = (IGeoFeatureLayer)pFLayer;
ITable pTable = (ITable)pFClass;
IClassifyGEN pClassifyGEN;
switch (cboGCClassify.Text)
{
case "Equal Interval":
pClassifyGEN = new EqualIntervalClass();
break;
case "Geometrical Interval":
pClassifyGEN = new GeometricalInterval();
break;
case "Natural Breaks":
pClassifyGEN = new NaturalBreaksClass();
break;
case "Quantile":
pClassifyGEN = new QuantileClass();
break;
case "StandardDeviation":
pClassifyGEN = new StandardDeviationClass();
break;
default:
pClassifyGEN = new NaturalBreaksClass();
break;
}
//Need to be changed 1/29/15
ITableHistogram pTableHistogram = new TableHistogramClass();
pTableHistogram.Field = strGCRenderField;
pTableHistogram.Table = pTable;
IHistogram pHistogram = (IHistogram)pTableHistogram;
object xVals, frqs;
pHistogram.GetHistogram(out xVals, out frqs);
pClassifyGEN.Classify(xVals, frqs, intGCBreakeCount);
ClassBreaksRenderer pRender = new ClassBreaksRenderer();
double[] cb = (double[])pClassifyGEN.ClassBreaks;
pRender.Field = strGCRenderField;
pRender.BreakCount = intGCBreakeCount;
pRender.MinimumBreak = cb[0];
//' create our color ramp
IAlgorithmicColorRamp pColorRamp = new AlgorithmicColorRampClass();
pColorRamp.Algorithm = esriColorRampAlgorithm.esriCIELabAlgorithm;
IRgbColor pColor1 = new RgbColor();
IRgbColor pColor2 = new RgbColor();
//Can Change the color in here!
pColor1.Red = picSymolfrom.BackColor.R;
pColor1.Green = picSymolfrom.BackColor.G;
pColor1.Blue = picSymolfrom.BackColor.B;
Boolean blnOK = true;
pColor2.Red = picSymbolTo.BackColor.R;
pColor2.Green = picSymbolTo.BackColor.G;
pColor2.Blue = picSymbolTo.BackColor.B;
pColorRamp.FromColor = pColor1;
pColorRamp.ToColor = pColor2;
pColorRamp.Size = intGCBreakeCount;
pColorRamp.CreateRamp(out blnOK);
IEnumColors pEnumColors = pColorRamp.Colors;
pEnumColors.Reset();
IRgbColor pColorOutline = new RgbColor();
//Can Change the color in here!
pColorOutline.Red = picGCLineColor.BackColor.R;
pColorOutline.Green = picGCLineColor.BackColor.G;
pColorOutline.Blue = picGCLineColor.BackColor.B;
double dblGCOutlineSize = Convert.ToDouble(nudGCLinewidth.Value);
ICartographicLineSymbol pOutLines = new CartographicLineSymbol();
pOutLines.Width = dblGCOutlineSize;
pOutLines.Color = (IColor)pColorOutline;
//' use this interface to set dialog properties
IClassBreaksUIProperties pUIProperties = (IClassBreaksUIProperties)pRender;
pUIProperties.ColorRamp = "Custom";
ISimpleFillSymbol pSimpleFillSym;
//' be careful, indices are different for the diff lists
for (int j = 0; j < intGCBreakeCount; j++)
{
pRender.Break[j] = cb[j + 1];
pRender.Label[j] = Math.Round(cb[j], 2).ToString() + " - " + Math.Round(cb[j + 1], 2).ToString();
pUIProperties.LowBreak[j] = cb[j];
pSimpleFillSym = new SimpleFillSymbolClass();
pSimpleFillSym.Color = pEnumColors.Next();
pSimpleFillSym.Outline = pOutLines;
pRender.Symbol[j] = (ISymbol)pSimpleFillSym;
}
pGeofeatureLayer.Renderer = (IFeatureRenderer)pRender;
//////////////////The Robustness
int intRoundingDigits = 2;
int intUncernBreakCount = Convert.ToInt32(nudTeNClasses.Value);
string strUncerFieldName = cboUField.Text;
int intUncernIdx = pFClass.FindField(strUncerFieldName);
int intValueIdx = pFClass.FindField(strGCRenderField);
//Calculate Robustness
//Add fld
int intTempfldIdx = 0;
string strTempfldName = txtFldName.Text;
if (chkRobustness.Checked)
{
if (pFClass.FindField(strTempfldName) == -1)
{
AddField(pFClass, strTempfldName);
}
intTempfldIdx = pFClass.FindField(strTempfldName);
}
Chart pChart = new Chart();
IFeature pFeat = null;
//IFeatureCursor pFCursor = pFClass.Search(null, false);
IFeatureCursor pFCursor = null;
if (chkRobustness.Checked)
{
pFCursor = pFClass.Update(null, false);
}
else
{
pFCursor = pFClass.Search(null, false);
}
pFeat = pFCursor.NextFeature();
double[] arrRobustness = new double[pFClass.FeatureCount(null)];
int i = 0;
while (pFeat != null)
{
for (int j = 0; j < (cb.Length - 1); j++)
{
double dblValue = Convert.ToDouble(pFeat.get_Value(intValueIdx));
double dblStd = Convert.ToDouble(pFeat.get_Value(intUncernIdx));
if (j == 0)
{
if (dblValue >= cb[j] && dblValue <= cb[j + 1])
{
double dblUpperZvalue = (cb[j + 1] - dblValue) / dblStd;
double dblUpperConfLev = pChart.DataManipulator.Statistics.NormalDistribution(dblUpperZvalue);
double dblLowerZvalue = (cb[j] - dblValue) / dblStd;
double dblLowerConfLev = pChart.DataManipulator.Statistics.NormalDistribution(dblLowerZvalue);
arrRobustness[i] = dblUpperConfLev - dblLowerConfLev;
if (chkRobustness.Checked)
{
pFeat.set_Value(intTempfldIdx, arrRobustness[i]);
}
}
}
else
{
if (dblValue > cb[j] && dblValue <= cb[j + 1])
{
double dblUpperZvalue = (cb[j + 1] - dblValue) / dblStd;
double dblUpperConfLev = pChart.DataManipulator.Statistics.NormalDistribution(dblUpperZvalue);
double dblLowerZvalue = (cb[j] - dblValue) / dblStd;
double dblLowerConfLev = pChart.DataManipulator.Statistics.NormalDistribution(dblLowerZvalue);
arrRobustness[i] = dblUpperConfLev - dblLowerConfLev;
if (chkRobustness.Checked)
{
pFeat.set_Value(intTempfldIdx, arrRobustness[i]);
}
}
}
}
if (chkRobustness.Checked)
{
pFCursor.UpdateFeature(pFeat);
}
i++;
pFeat = pFCursor.NextFeature();
}
//Define the intervals (the last class is fixed to 1)
if (intUncernBreakCount == 1)
{
return;
}
double[] arrRobustBrks = new double[intUncernBreakCount + 1];
double dblRBrksIntervals = Math.Round(1 / Convert.ToDouble(intUncernBreakCount - 1), intRoundingDigits);
arrRobustBrks[0] = 0;
for (int j = 1; j < intUncernBreakCount; j++)
{
arrRobustBrks[j] = dblRBrksIntervals * j;
}
arrRobustBrks[intUncernBreakCount] = 1;
IFeatureLayer pflUncern = new FeatureLayerClass();
pflUncern.FeatureClass = pFClass;
pflUncern.Name = "Robustness";
pflUncern.Visible = true;
IGeoFeatureLayer pGFLUncern = (IGeoFeatureLayer)pflUncern;
pFCursor = pGFLUncern.Search(null, true);
RobustnessRenderer pRobustnessRenderer = new RobustnessRenderer();
pRobustnessRenderer.arrRobustBrks = arrRobustBrks;
pRobustnessRenderer.arrRobustness = arrRobustness;
pRobustnessRenderer.dblAngle = Convert.ToDouble(nudAngleFrom.Value);
pRobustnessRenderer.dblFromSep = Convert.ToDouble(nudSeperationFrom.Value);
pRobustnessRenderer.dblLinewidth = Convert.ToDouble(nudTeLinewidth.Value);
pRobustnessRenderer.dblToSep = Convert.ToDouble(nudSeperationTo.Value);
pRobustnessRenderer.intUncernBreakCount = intUncernBreakCount;
pRobustnessRenderer.pLineColor = pSnippet.getRGB(picTeLineColor.BackColor.R, picTeLineColor.BackColor.G, picTeLineColor.BackColor.B);
IQueryFilter pQFilter = new QueryFilterClass();
pRobustnessRenderer.PrepareFilter(pFClass, pQFilter);
pRobustnessRenderer.Draw(pFCursor, esriDrawPhase.esriDPSelection, pActiveView.ScreenDisplay, null);
pRobustnessRenderer.CreateLegend();
pGFLUncern.Renderer = pRobustnessRenderer;
pActiveView.FocusMap.AddLayer(pGFLUncern as ILayer);
mForm.axMapControl1.ActiveView.Refresh();
mForm.axTOCControl1.Update();
}
