public void Execute(IArray paramvalues, ITrackCancel trackcancel, IGPEnvironmentManager envMgr, IGPMessages message)
{
// Get Parameters
IGPParameter3 inputParameter = (IGPParameter3)paramvalues.get_Element(0);
IGPParameter3 polygonParameter = (IGPParameter3)paramvalues.get_Element(1);
IGPParameter3 outputParameter = (IGPParameter3)paramvalues.get_Element(2);
IGPParameter3 fieldParameter = (IGPParameter3)paramvalues.get_Element(3);
// UnPackGPValue. This ensures you get the value from either the dataelement or GpVariable (ModelBuilder)
IGPValue inputParameterValue = m_GPUtilities.UnpackGPValue(inputParameter);
IGPValue polygonParameterValue = m_GPUtilities.UnpackGPValue(polygonParameter);
IGPValue outputParameterValue = m_GPUtilities.UnpackGPValue(outputParameter);
IGPValue fieldParameterValue = m_GPUtilities.UnpackGPValue(fieldParameter);
// Decode Input Feature Layers
IFeatureClass inputFeatureClass;
IFeatureClass polygonFeatureClass;
IQueryFilter inputFeatureClassQF;
IQueryFilter polygonFeatureClassQF;
m_GPUtilities.DecodeFeatureLayer(inputParameterValue, out inputFeatureClass, out inputFeatureClassQF);
m_GPUtilities.DecodeFeatureLayer(polygonParameterValue, out polygonFeatureClass, out polygonFeatureClassQF);
if (inputFeatureClass == null)
{
message.AddError(2, "Could not open input dataset.");
return;
}
if (polygonFeatureClass == null)
{
message.AddError(2, "Could not open clipping polygon dataset.");
return;
}
if (polygonFeatureClass.FeatureCount(null) > 1)
{
message.AddWarning("Clipping polygon feature class contains more than one feature.");
}
// Create the Geoprocessor
Geoprocessor gp = new Geoprocessor();
// Create Output Polygon Feature Class
CreateFeatureclass cfc = new CreateFeatureclass();
IName name = m_GPUtilities.CreateFeatureClassName(outputParameterValue.GetAsText());
IDatasetName dsName = name as IDatasetName;
IFeatureClassName fcName = dsName as IFeatureClassName;
IFeatureDatasetName fdsName = fcName.FeatureDatasetName as IFeatureDatasetName;
// Check if output is in a FeatureDataset or not. Set the output path parameter for CreateFeatureClass tool.
if (fdsName != null)
{
cfc.out_path = fdsName;
}
else
{
cfc.out_path = dsName.WorkspaceName.PathName;
}
// Set the output Coordinate System for CreateFeatureClass tool.
// ISpatialReference3 sr = null;
IGPEnvironment env = envMgr.FindEnvironment("outputCoordinateSystem");
// Same as Input
if (env.Value.IsEmpty())
{
IGeoDataset ds = inputFeatureClass as IGeoDataset;
cfc.spatial_reference = ds.SpatialReference as ISpatialReference3;
}
// Use the environment setting
else
{
IGPCoordinateSystem cs = env.Value as IGPCoordinateSystem;
cfc.spatial_reference = cs.SpatialReference as ISpatialReference3;
}
// Remaining properties for Create Feature Class Tool
cfc.out_name = dsName.Name;
cfc.geometry_type = "POLYGON";
// Execute Geoprocessor
gp.Execute(cfc, null);
// Get Unique Field
int iField = inputFeatureClass.FindField(fieldParameterValue.GetAsText());
IField uniqueField = inputFeatureClass.Fields.get_Field(iField);
// Extract Clipping Polygon Geometry
IFeature polygonFeature = polygonFeatureClass.GetFeature(0);
IPolygon clippingPolygon = (IPolygon)polygonFeature.Shape;
// Spatial Filter
ISpatialFilter spatialFilter = new SpatialFilterClass();
spatialFilter.Geometry = polygonFeature.ShapeCopy;
spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelContains;
// Debug Message
message.AddMessage("Generating TIN...");
// Create TIN
ITinEdit tinEdit = new TinClass();
// Advanced TIN Functions
ITinAdvanced2 tinAdv = (ITinAdvanced2)tinEdit;
try
{
// Initialize New TIN
IGeoDataset gds = inputFeatureClass as IGeoDataset;
tinEdit.InitNew(gds.Extent);
// Add Mass Points to TIN
tinEdit.StartEditing();
tinEdit.AddFromFeatureClass(inputFeatureClass, spatialFilter, uniqueField, uniqueField, esriTinSurfaceType.esriTinMassPoint);
tinEdit.Refresh();
// Get TIN Nodes
ITinNodeCollection tinNodeCollection = (ITinNodeCollection)tinEdit;
// Report Node Count
message.AddMessage("Input Node Count: " + inputFeatureClass.FeatureCount(null).ToString());
message.AddMessage("TIN Node Count: " + tinNodeCollection.NodeCount.ToString());
// Open Output Feature Class
IFeatureClass outputFeatureClass = m_GPUtilities.OpenFeatureClassFromString(outputParameterValue.GetAsText());
// Debug Message
message.AddMessage("Generating Polygons...");
// Create Voronoi Polygons
tinNodeCollection.ConvertToVoronoiRegions(outputFeatureClass, null, clippingPolygon, "", "");
// Release COM Objects
tinEdit.StopEditing(false);
System.Runtime.InteropServices.Marshal.ReleaseComObject(tinNodeCollection);
System.Runtime.InteropServices.Marshal.ReleaseComObject(tinEdit);
}
catch (Exception ex)
{
message.AddError(2, ex.Message);
}
}
public static IFeatureClass CreateFeatureClass(IGPValue gpFeatureClass, IGPEnvironmentManager environment, IFields fields = null)
{
if (gpFeatureClass == null)
throw new ArgumentException("Argument can not be null", "gpFeatureClass");
if (environment == null)
throw new ArgumentException("Argument can not be null", "environment");
IGeoProcessorSettings gpSettings = (IGeoProcessorSettings)environment;
if (gpUtilities.Exists(gpFeatureClass))
{
if (gpSettings.OverwriteOutput == true)
gpUtilities.Delete(gpFeatureClass);
else
throw new Exception("Output feature class already exists: " + gpFeatureClass.GetAsText());
}
IDEFeatureClass deFeatureClass = (IDEFeatureClass)gpUtilities.DecodeDETable(gpFeatureClass);
if (deFeatureClass == null)
throw new Exception("Data Element decode return null");
IObjectClassDescription objectClassDescription = (IObjectClassDescription)new FeatureClassDescriptionClass();
UID clsid = objectClassDescription.InstanceCLSID;
UID extclsid = objectClassDescription.ClassExtensionCLSID;
IDataElement dataElement = (IDataElement)deFeatureClass;
if (dataElement.CatalogPath == null)
throw new ArgumentException("Catalog path is null", "CatalogPath");
IFeatureClassName featureClassName = (IFeatureClassName)gpUtilities.CreateFeatureClassName(dataElement.CatalogPath);
string path = dataElement.GetPath();
string name = dataElement.Name;
IDEGeoDataset geoDataElement = (IDEGeoDataset)deFeatureClass;
ISpatialReference spatialReference = geoDataElement.SpatialReference;
IDETable deTable = (IDETable)deFeatureClass;
string shapeFieldName = deFeatureClass.ShapeFieldName;
Dictionary<string,IField> fieldBuilder = new Dictionary<string,IField>();
foreach(var input_fields in new IFields[] { deTable.Fields, fields })
{
if (input_fields == null) continue;
for (int i = 0; i < input_fields.FieldCount; i++)
{
IField field = deTable.Fields.get_Field(i);
if (fieldBuilder.ContainsKey(field.Name.ToLower()))
fieldBuilder[field.Name.ToLower()] = (IField)((IClone)field).Clone();
else
fieldBuilder.Add(field.Name.ToLower(), (IField)((IClone)field).Clone());
if (field.Type == esriFieldType.esriFieldTypeGeometry)
{
shapeFieldName = field.Name;
break;
}
}
}
IFields output_fields = new FieldsClass();
IFieldsEdit fields_edit = (IFieldsEdit)output_fields;
foreach(IField field in fieldBuilder.Values)
{
fields_edit.AddField(field);
if(field.Type == esriFieldType.esriFieldTypeGeometry)
{
IGeometryDefEdit defEdit = (IGeometryDefEdit)field.GeometryDef;
defEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;
}
}
string configKeyword = ((IGPString)environment.FindEnvironment(CONFIG_KEYWORD_PROP_NAME).Value).Value;
//if (String.IsNullOrWhiteSpace(configKeyword)) configKeyword = "DEFAULTS";
IFeatureClass ret = null;
if(featureClassName.FeatureDatasetName != null)
{
IFeatureDataset featureDataset = (IFeatureDataset)((IName)featureClassName.FeatureDatasetName).Open();
try
{
ret = featureDataset.CreateFeatureClass(name, output_fields, clsid, extclsid, esriFeatureType.esriFTSimple, shapeFieldName, configKeyword);
}
finally
{
Marshal.ReleaseComObject(featureDataset);
}
}
else
{
IWorkspace workspace = (IWorkspace)((IName)((IDatasetName)featureClassName).WorkspaceName).Open();
try
{
IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspace;
ret = featureWorkspace.CreateFeatureClass(name, output_fields, clsid, extclsid, esriFeatureType.esriFTSimple, shapeFieldName, configKeyword);
}
finally
{
Marshal.ReleaseComObject(workspace);
}
}
return ret;
}
public static IFeatureClass CreateFeatureClass(IGPValue gpFeatureClass, IGPEnvironmentManager environment, IFields fields = null)
{
if (gpFeatureClass == null)
{
throw new ArgumentException("Argument can not be null", "gpFeatureClass");
}
if (environment == null)
{
throw new ArgumentException("Argument can not be null", "environment");
}
IGeoProcessorSettings gpSettings = (IGeoProcessorSettings)environment;
if (gpUtilities.Exists(gpFeatureClass))
{
if (gpSettings.OverwriteOutput == true)
{
gpUtilities.Delete(gpFeatureClass);
}
else
{
throw new Exception("Output feature class already exists: " + gpFeatureClass.GetAsText());
}
}
IDEFeatureClass deFeatureClass = (IDEFeatureClass)gpUtilities.DecodeDETable(gpFeatureClass);
if (deFeatureClass == null)
{
throw new Exception("Data Element decode return null");
}
IObjectClassDescription objectClassDescription = (IObjectClassDescription) new FeatureClassDescriptionClass();
UID clsid = objectClassDescription.InstanceCLSID;
UID extclsid = objectClassDescription.ClassExtensionCLSID;
IDataElement dataElement = (IDataElement)deFeatureClass;
if (dataElement.CatalogPath == null)
{
throw new ArgumentException("Catalog path is null", "CatalogPath");
}
IFeatureClassName featureClassName = (IFeatureClassName)gpUtilities.CreateFeatureClassName(dataElement.CatalogPath);
string path = dataElement.GetPath();
string name = dataElement.Name;
IDEGeoDataset geoDataElement = (IDEGeoDataset)deFeatureClass;
ISpatialReference spatialReference = geoDataElement.SpatialReference;
IDETable deTable = (IDETable)deFeatureClass;
string shapeFieldName = deFeatureClass.ShapeFieldName;
Dictionary <string, IField> fieldBuilder = new Dictionary <string, IField>();
foreach (var input_fields in new IFields[] { deTable.Fields, fields })
{
if (input_fields == null)
{
continue;
}
for (int i = 0; i < input_fields.FieldCount; i++)
{
IField field = deTable.Fields.get_Field(i);
if (fieldBuilder.ContainsKey(field.Name.ToLower()))
{
fieldBuilder[field.Name.ToLower()] = (IField)((IClone)field).Clone();
}
else
{
fieldBuilder.Add(field.Name.ToLower(), (IField)((IClone)field).Clone());
}
if (field.Type == esriFieldType.esriFieldTypeGeometry)
{
shapeFieldName = field.Name;
break;
}
}
}
IFields output_fields = new FieldsClass();
IFieldsEdit fields_edit = (IFieldsEdit)output_fields;
foreach (IField field in fieldBuilder.Values)
{
fields_edit.AddField(field);
if (field.Type == esriFieldType.esriFieldTypeGeometry)
{
IGeometryDefEdit defEdit = (IGeometryDefEdit)field.GeometryDef;
defEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;
}
}
string configKeyword = ((IGPString)environment.FindEnvironment(CONFIG_KEYWORD_PROP_NAME).Value).Value;
//if (String.IsNullOrWhiteSpace(configKeyword)) configKeyword = "DEFAULTS";
IFeatureClass ret = null;
if (featureClassName.FeatureDatasetName != null)
{
IFeatureDataset featureDataset = (IFeatureDataset)((IName)featureClassName.FeatureDatasetName).Open();
try
{
ret = featureDataset.CreateFeatureClass(name, output_fields, clsid, extclsid, esriFeatureType.esriFTSimple, shapeFieldName, configKeyword);
}
finally
{
Marshal.ReleaseComObject(featureDataset);
}
}
else
{
IWorkspace workspace = (IWorkspace)((IName)((IDatasetName)featureClassName).WorkspaceName).Open();
try
{
IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspace;
ret = featureWorkspace.CreateFeatureClass(name, output_fields, clsid, extclsid, esriFeatureType.esriFTSimple, shapeFieldName, configKeyword);
}
finally
{
Marshal.ReleaseComObject(workspace);
}
}
return(ret);
}
