public IFeatureClass GetFeatureClass(string path, string fileName)
{
IFeatureClass ifc = null;
try
{
IWorkspaceFactory wrkSpaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace featureWorkspace = null;
featureWorkspace = (IFeatureWorkspace)wrkSpaceFactory.OpenFromFile(path, 0);
ifc = featureWorkspace.OpenFeatureClass(fileName);
}
catch (COMException COMEx)
{
eLog.Debug(COMEx);
//System.Windows.Forms.MessageBox.Show(COMEx.GetBaseException().ToString(), "COM Error: " + COMEx.ErrorCode.ToString());
}
catch (System.Exception ex)
{
eLog.Debug(ex);
//System.Windows.Forms.MessageBox.Show(ex.Source + " ");//+ ex.InnerException.ToString());
}
return(ifc);
}
//注册要素类
/// <summary>
/// 返回true说明FeatureClass存在，返回false说明不存在，重新创建
/// </summary>
/// <param name="IN_ShapePath"></param>
/// <returns></returns>
private bool PRV_AddFeatureClass(string IN_ShapePath)
{
string Temp_Direction = System.IO.Path.GetDirectoryName(IN_ShapePath); //该Shp文件的目录
string Temp_Name = System.IO.Path.GetFileNameWithoutExtension(IN_ShapePath); //该Shp文件的名称
IWorkspaceFactory Temp_ShapeWorkFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace Temp_ShapeWorkspace = Temp_ShapeWorkFactory.OpenFromFile(Temp_Direction, 0) as IFeatureWorkspace;
IWorkspaceFactory Temp_AccessWorkFactory = new AccessWorkspaceFactory();
IFeatureWorkspace Temp_Workspace = Temp_AccessWorkFactory.OpenFromFile(S_MDBFile, 0) as IFeatureWorkspace;
IFeatureClassContainer tem_FeatureClassContainer = (IFeatureClassContainer)FDS_Featuredataset;
IEnumFeatureClass pEnumFeatureClass = (IEnumFeatureClass)tem_FeatureClassContainer.Classes;
IFeatureClass tem_FeatureClass = pEnumFeatureClass.Next();
while (null != tem_FeatureClass)
{
if (Temp_Name == tem_FeatureClass.AliasName)
{// return true;
}
tem_FeatureClass = pEnumFeatureClass.Next();
}
IFeatureClass Temp_FeatureClass = Temp_ShapeWorkspace.OpenFeatureClass(Temp_Name);
FeatureClassToFeatureClass Temp_FCToFC = new FeatureClassToFeatureClass(IN_ShapePath, S_MDBFile + "\\" + FDS_Featuredataset.Name, Temp_Name);//将Shp文件导入要素数据集
GP_Progress = GP_Tool.ExecuteAsync(Temp_FCToFC);
TH_TimeSpan = new Thread(PRV_GetStatus);//开辟线程计时
TH_TimeSpan.Start();
TH_TimeSpan.Join();
return(false);
// IFeatureClassContainer ss = (FDS_Featuredataset.Workspace as IFeatureWorkspace).OpenFeatureDataset(FDS_Featuredataset.Name) as IFeatureClassContainer;
// Console.WriteLine("完成");
}
public static IFeatureLayer OpenShapefile(ref string shapefilePath)
{
OpenFileDialog dg = new OpenFileDialog();
dg.Filter = "Shapefile文件(.shp)|*.shp";
dg.Multiselect = true;
if (dg.ShowDialog() == System.Windows.Forms.DialogResult.OK)
{
string fileName = dg.FileName;
shapefilePath = fileName;
if (string.IsNullOrEmpty(fileName))
{
MessageBox.Show("请选择要打开的shapefile文件");
return(null);
}
IWorkspaceFactory pWksFactory = new ShapefileWorkspaceFactory();
IWorkspace pWks = pWksFactory.OpenFromFile(System.IO.Path.GetDirectoryName(fileName), 0);
IFeatureWorkspace pFeatureWks = (IFeatureWorkspace)pWks;
IFeatureClass pFc = pFeatureWks.OpenFeatureClass(System.IO.Path.GetFileNameWithoutExtension(fileName));
IFeatureLayer pFlayer = new FeatureLayerClass();
pFlayer.FeatureClass = pFc;
pFlayer.Name = pFc.AliasName;
return(pFlayer);
}
return(null);
}
public IFeatureClass GetFeatureClass(string inFileName)
{
string path = string.Empty;
string fileName = string.Empty;
IFeatureClass ifc = null;
try
{
this.GetPathAndFileName(inFileName, out path, out fileName);
IWorkspaceFactory wrkSpaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace featureWorkspace = null;
featureWorkspace = (IFeatureWorkspace)wrkSpaceFactory.OpenFromFile(path, 0);
ifc = featureWorkspace.OpenFeatureClass(fileName);
}
catch (COMException COMEx)
{
eLog.Debug(COMEx);
eLog.Debug("Orginal File Name was " + inFileName);
eLog.Debug("so the path is " + path);
eLog.Debug("just the file name is" + fileName);
//System.Windows.Forms.MessageBox.Show(COMEx.GetBaseException().ToString(), "COM Error: " + COMEx.ErrorCode.ToString());
}
catch (System.Exception ex)
{
eLog.Debug(ex);
//System.Windows.Forms.MessageBox.Show(ex.Source + " ");//+ ex.InnerException.ToString());
}
return(ifc);
}
//添加导入的要素类
private void button5_Click(object sender, EventArgs e)
{
int index;
string path, name;
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace2;
IFeatureWorkspace pFeatureWorkspace;
OpenFileDialog openFileDialog = new OpenFileDialog();
openFileDialog.Title = "选择要添加的要素类";
openFileDialog.Multiselect = true;
openFileDialog.Filter = "shp文件(*.shp)|*.shp";
if (openFileDialog.ShowDialog() == DialogResult.OK)
{
List <string> FileList = new List <string>();
foreach (string s in openFileDialog.FileNames)
{
listBox1.Items.Add(s);
index = s.LastIndexOf("\\");
path = s.Substring(0, index);
name = s.Substring(index + 1);
pWorkspace2 = pWorkspaceFactory.OpenFromFile(path, 0);
pFeatureWorkspace = pWorkspace2 as IFeatureWorkspace;
pFeatureClassList.Add(pFeatureWorkspace.OpenFeatureClass(name));
}
}
}
private void UpdateSpatialInfoElements(string strSPInfo)
{
string strTableName = txtFlowFile.Text;
string strSourceFolder = null;
string strSourceFile = null;
if (strTableName == "")
{
return;
}
else
{
strSourceFolder = System.IO.Path.GetDirectoryName(strTableName);
strSourceFile = System.IO.Path.GetFileName(strTableName);
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace ws = workspaceFactory.OpenFromFile(strSourceFolder, 0);
IFeatureWorkspace fws = ws as IFeatureWorkspace;
m_tblFlow = fws.OpenTable(strSourceFile);
IFields fields = m_tblFlow.Fields;
if (strSPInfo == "From ShapeFile")
{
cbo1.Items.Clear();
cbo2.Items.Clear();
cbo3.Items.Clear();
cbo4.Items.Clear();
for (int i = 0; i < m_pForm.axMapControl1.LayerCount; i++)
{
cbo1.Items.Add(m_pForm.axMapControl1.get_Layer(i).Name);
}
for (int i = 0; i < fields.FieldCount; i++)
{
cbo3.Items.Add(fields.get_Field(i).Name);
cbo4.Items.Add(fields.get_Field(i).Name);
}
}
else if (strSPInfo == "From Table")
{
cbo1.Items.Clear();
cbo2.Items.Clear();
cbo3.Items.Clear();
cbo4.Items.Clear();
for (int i = 0; i < fields.FieldCount; i++)
{
if (m_pSnippet.FindNumberFieldType(fields.get_Field(i)))
{
cbo1.Items.Add(fields.get_Field(i).Name);
cbo2.Items.Add(fields.get_Field(i).Name);
cbo3.Items.Add(fields.get_Field(i).Name);
cbo4.Items.Add(fields.get_Field(i).Name);
}
}
}
}
}
protected IFeatureClass OpenShapefile(string path, string Name, int hWnd)
{
IFeatureWorkspace pFWS;
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
pFWS = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(path, hWnd); //, Application.hWnd);
return(pFWS.OpenFeatureClass(Name));
}
public static IFeatureWorkspace CreateFeatureWorkspace(string shapefilePath)
{
IWorkspaceFactory pWksFactory = new ShapefileWorkspaceFactory();
IWorkspace pWks = pWksFactory.OpenFromFile(System.IO.Path.GetDirectoryName(shapefilePath), 0);
IFeatureWorkspace pFeatureWks = pWks as IFeatureWorkspace;
return(pFeatureWks);
}
public static string GetNameShapeFile(string gdb_Path, string file_Name)
{
IWorkspaceFactory wsf = new ShapefileWorkspaceFactory();
IWorkspace ws = wsf.OpenFromFile(gdb_Path, 0);
IEnumDatasetName edn = ws.get_DatasetNames(esriDatasetType.esriDTFeatureClass);
return(Check(edn, file_Name, string.Empty));
}
public static void Mlayer_IDW_Click()
{
// 用反距离IDW插值生成的栅格图像。如下：
IInterpolationOp pInterpolationOp = new RasterInterpolationOpClass();
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
string pPath = Application.StartupPath + @"\MakeContours\Cont.shp";
string pFolder = System.IO.Path.GetDirectoryName(pPath);
string pFileName = System.IO.Path.GetFileName(pPath);
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pFolder, 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass oFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
IFeatureClassDescriptor pFCDescriptor = new FeatureClassDescriptorClass();
pFCDescriptor.Create(oFeatureClass, null, "shape.z");
IRasterRadius pRadius = new RasterRadiusClass();
object objectMaxDistance = null;
object objectbarrier = null;
object missing = Type.Missing;
pRadius.SetVariable(12, ref objectMaxDistance);
object dCellSize = 1;
object snapRasterData = Type.Missing;
IEnvelope pExtent;
pExtent = new EnvelopeClass();
Double xmin = 27202;
Double xmax = 31550;
Double ymin = 19104;
Double ymax = 22947;
pExtent.PutCoords(xmin, ymin, xmax, ymax);
object extentProvider = pExtent;
IRasterAnalysisEnvironment pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCellSize);
pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extentProvider, ref snapRasterData);
IGeoDataset poutGeoDataset = pInterpolationOp.IDW((IGeoDataset)pFCDescriptor, 2, pRadius, ref objectbarrier);
ISurfaceOp surOp = new RasterSurfaceOpClass();
IRaster pOutRaster = poutGeoDataset as IRaster;
IRasterLayer pOutRasLayer = new RasterLayer();
pOutRasLayer.CreateFromRaster(pOutRaster);
IMap pMap = Common.DataEditCommon.g_pMap;
pMap.AddLayer(pOutRasLayer);
Common.DataEditCommon.g_axTocControl.Refresh();
Common.DataEditCommon.g_pAxMapControl.ActiveView.Refresh();
}
//编写读取shapefile数据的方法：
public IFeatureClass OpenShapeFeatureClass(string shpPath, string shpName)
{
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workspaceFactory.OpenFromFile(shpPath, 0); //调用工厂的方法；
IFeatureWorkspace shpWorkspace = (IFeatureWorkspace)workspace; //强制转换接口；
IFeatureClass featureClass = shpWorkspace.OpenFeatureClass(shpName); //声明一个要素类对象，调用工作空间的方法打开要素并赋值，传入参数为文件名
return(featureClass);
}
protected bool ShowSaveFileDialog(IGxDialog pGxDialog, out string path, out string name, out string genericName, out string datasetType)
{
path = ""; name = ""; genericName = ""; datasetType = "";
pGxDialog.Title = "Save File as:";
if (!pGxDialog.DoModalSave(0))
{
return(false);
}
datasetType = pGxDialog.ObjectFilter.Name;
// ' delete the existing shapefile if user wants to replace it.
if (pGxDialog.ReplacingObject)
{
switch (datasetType)
{
case "Shapefiles":
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFWS = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(pGxDialog.FinalLocation.FullName, 0);// mFile(pGxDialullName, 0)
IFeatureClass pFeatClass = pFWS.OpenFeatureClass(pGxDialog.Name);
if (!DeleteShapefile(pFeatClass))
{
MessageBox.Show("Please specify a different name for the output shapefile", "Error", MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return(false);
}
break;
}
}
if (!(pGxDialog.Name.IndexOf(".") == 0))
{
genericName = pGxDialog.Name;
}
else
{
genericName = pGxDialog.Name.Substring(0, pGxDialog.Name.IndexOf("."));
}
//if (sTelNo.IndexOf('/') >= 0)
//{
// sTelNo = sTelNo.Substring(0, sTelNo.IndexOf('/'));
//}
switch (datasetType)
{
case "Shapefiles":
name = genericName + ".shp";
break;
case "Raster Datasets":
name = genericName + ".img";
break;
}
path = pGxDialog.FinalLocation.FullName;
return(true);
}
private void button1_Click(object sender, EventArgs e)
{
IWorkspaceFactory pWsf = new ShapefileWorkspaceFactory();
IWorkspace pWs = pWsf.OpenFromFile("",0);
IFeatureWorkspace pFws = pWs as IFeatureWorkspace;
IFeatureClass pFc = pFws.OpenFeatureClass("");
//pWs = new WorkspaceClass();
}
private static void DeleteShpFile(string pathName, string fileName)
{
IWorkspaceFactory pWSF = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFWS = pWSF.OpenFromFile(pathName, 0) as IFeatureWorkspace;
IFeatureClass pFeatureClass = pFWS.OpenFeatureClass(fileName);
IDataset pFeatureDataset = pFeatureClass as IDataset;
pFeatureDataset.Delete();
}
public static void Mlayer_Krige_Click()
{
// 用克里金Krige插值生成的栅格图像。如下：
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
string pPath = Application.StartupPath + @"\MakeContours\Cont.shp";
string pFolder = System.IO.Path.GetDirectoryName(pPath);
string pFileName = System.IO.Path.GetFileName(pPath);
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pFolder, 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass oFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
IFeatureClassDescriptor pFCDescriptor = new FeatureClassDescriptorClass();
pFCDescriptor.Create(oFeatureClass, null, "shape.z");
IInterpolationOp pInterpolationOp = new RasterInterpolationOpClass();
IRasterAnalysisEnvironment pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
object Cellsize = 0.004;//Cell size for output raster;0.004
pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref Cellsize);
//设置输出范围
//27202 19104;27202 22947;31550 22947;31550 19104
object snapRasterData = Type.Missing;
IEnvelope pExtent;
pExtent = new EnvelopeClass();
Double xmin = 27202;
Double xmax = 31550;
Double ymin = 19104;
Double ymax = 22947;
pExtent.PutCoords(xmin, ymin, xmax, ymax);
object extentProvider = pExtent;
pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extentProvider, ref snapRasterData);
Double dSearchD = 10;
object pSearchCount = 3;
object missing = Type.Missing;
IRasterRadius pRadius = new RasterRadius();
pRadius.SetFixed(dSearchD, ref pSearchCount);
//pRadius.SetVariable((int)pSearchCount, ref dSearchD);
IGeoDataset poutGeoDataset = pInterpolationOp.Krige((IGeoDataset)pFCDescriptor, esriGeoAnalysisSemiVariogramEnum.esriGeoAnalysisGaussianSemiVariogram, pRadius, false, ref missing);
IRaster pOutRaster = poutGeoDataset as IRaster;
IRasterLayer pOutRasLayer = new RasterLayer();
pOutRasLayer.CreateFromRaster(pOutRaster);
IMap pMap = Common.DataEditCommon.g_pMap;
pMap.AddLayer(pOutRasLayer);
Common.DataEditCommon.g_axTocControl.Refresh();
Common.DataEditCommon.g_pAxMapControl.ActiveView.Refresh();
}
public void Test()
{
//定义一个几何字段，类型为点类型
ISpatialReferenceFactory2 originalSpatialReferenceFactory = new SpatialReferenceEnvironmentClass();
ISpatialReference originalSpatialReference = originalSpatialReferenceFactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
IGeometryDefEdit pGeoDef = new GeometryDefClass();
IGeometryDefEdit pGeoDefEdit = pGeoDef as IGeometryDefEdit;
pGeoDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
pGeoDefEdit.SpatialReference_2 = originalSpatialReference;
//定义一个字段集合对象
IFields pFields = new FieldsClass();
IFieldsEdit pFieldsEdit = (IFieldsEdit)pFields;
//定义shape字段
IField pField = new FieldClass();
IFieldEdit pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "SHAPE";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
pFieldsEdit.AddField(pField);
pFieldEdit.GeometryDef_2 = pGeoDef;
//定义单个的字段，并添加到字段集合中
pField = new FieldClass();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "STCD";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pFieldsEdit.AddField(pField);
//定义单个的字段，并添加到字段集合中
pField = new FieldClass();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "SLM10";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pFieldsEdit.AddField(pField);
//定义单个的字段，并添加到字段集合中
pField = new FieldClass();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "SLM20";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pFieldsEdit.AddField(pField);
//定义单个的字段，并添加到字段集合中
pField = new FieldClass();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "SLM40";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pFieldsEdit.AddField(pField);
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFWS = pWorkspaceFactory.OpenFromFile(workspaceDirectory, 0) as IFeatureWorkspace;
IFeatureClass pFtClass = pFWS.CreateFeatureClass("Test", pFields, null, null, esriFeatureType.esriFTSimple, "SHAPE", null);
}
private void btnSHP_Click(object sender, EventArgs e)
{
OpenFileDialog openFile = new OpenFileDialog();
openFile.Filter = "*.shp|*.shp";
openFile.Title = "选择shp文件";
openFile.Multiselect = true;
ESRI.ArcGIS.Geometry.ISpatialReference pSpa = null;
if (openFile.ShowDialog() == DialogResult.OK)
{
string[] strs = openFile.FileNames;
this.lstLyrFile.Items.Clear();
for (int i = 0; i < strs.GetLength(0); i++)
{
string str = strs[i];
this.lstLyrFile.Items.Add(System.IO.Path.GetFileNameWithoutExtension(str));
//获得一个源数据的空间参考
//获得空间参看 源的
if (pSpa == null)
{
IWorkspaceFactory pFac = new ShapefileWorkspaceFactory();
IWorkspace pWks = pFac.OpenFromFile(System.IO.Path.GetDirectoryName(str), 0);
if (pWks == null)
{
continue;
}
IFeatureWorkspace pFeaWks = pWks as IFeatureWorkspace;
IFeatureClass pFeaCls = pFeaWks.OpenFeatureClass(System.IO.Path.GetFileNameWithoutExtension(str));
if (pFeaCls == null)
{
continue;
}
IGeoDataset pGeoData = pFeaCls as IGeoDataset;
pSpa = pGeoData.SpatialReference;
if (pSpa == null)
{
continue;
}
m_strInPrj = ExportToESRISpatialReference(pSpa);
pFeaWks = null;
}
}
for (int i = 0; i < this.lstLyrFile.Items.Count; i++)
{
this.lstLyrFile.Items[i].Checked = true;
}
this.txtSource.Text = System.IO.Path.GetDirectoryName(strs[0]);
}
}
protected override void OnClick()
{
//
// TODO: Sample code showing how to access button host
//
var referenceFactory2 = (ISpatialReferenceFactory2) new SpatialReferenceEnvironment();
ISpatialReference WGS84 = referenceFactory2.CreateSpatialReference((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
IPointCollection pointCollection = new PolylineClass();
// ------------- Ensure that both points have negative longitude values -------------------
IPoint point = new PointClass();
point.PutCoords(-170, 10); // Equivalent to 170 degrees WEST
point.SpatialReference = WGS84;
pointCollection.AddPoint(point);
point = new PointClass();
point.PutCoords(-200, 10); // Equivalent to 160 degrees EAST
point.SpatialReference = WGS84;
pointCollection.AddPoint(point);
// -----------------------------------------------------------------------
IPolyline polyline = (IPolyline)pointCollection;
polyline.SpatialReference = WGS84;
var geometryDefEdit = (IGeometryDefEdit) new GeometryDef();
geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolyline;
geometryDefEdit.SpatialReference_2 = WGS84;
var field = (IFieldEdit) new Field();
field.Name_2 = "Shape";
field.Type_2 = esriFieldType.esriFieldTypeGeometry;
field.GeometryDef_2 = geometryDefEdit;
var fields = (IFieldsEdit) new Fields();
fields.AddField(field);
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
var featureWorkspace = (IFeatureWorkspace)workspaceFactory.OpenFromFile("C:\\Temp\\", 0);
IFeatureClass featureClass = featureWorkspace.CreateFeatureClass
("test3.shp", fields, null, null, esriFeatureType.esriFTSimple, "Shape", "");
var feature = featureClass.CreateFeature();
feature.Shape = polyline;
feature.Store();
ArcMap.Application.CurrentTool = null;
}
private IFeatureClass openShapefile(string pFilePath, string pFileName)
{
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pFilePath, 0);
IFeatureWorkspace pFeatureWorkspace = (IFeatureWorkspace)pWorkspace;
IFeatureClass pFeatureClass;
pFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
return(pFeatureClass);
}
public void PracticeExercise2()
{
IMxDocument pMxDoc;
pMxDoc = (IMxDocument)ArcMap.Application.Document;
IMap pMap;
pMap = pMxDoc.FocusMap;
IWorkspaceFactory pWSFactory;
pWSFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace;
pWorkspace = pWSFactory.OpenFromFile("c:\\WCGIS\\G5224P\\Lesson1_2_3", ArcMap.Application.hWnd);
//QI
IFeatureWorkspace pFWorkspace;
pFWorkspace = (IFeatureWorkspace)pWorkspace;
IFeatureClass pStatesFClass;
pStatesFClass = pFWorkspace.OpenFeatureClass("us_boundaries");
IFeatureLayer pStatesLayer;
pStatesLayer = new FeatureLayer();
pStatesLayer.FeatureClass = pStatesFClass;
pStatesLayer.Name = "U.S. States";
pMap.AddLayer(pStatesLayer);
IFeatureClass pCitiesFClass;
pCitiesFClass = pFWorkspace.OpenFeatureClass("us_cities");
IFeatureLayer pCitiesLayer;
pCitiesLayer = new FeatureLayer();
pCitiesLayer.FeatureClass = pCitiesFClass;
pCitiesLayer.Name = "U.S. Cities";
pMap.AddLayer(pCitiesLayer);
IActiveView pActiveView;
pActiveView = (IActiveView)pMap;
pActiveView.Refresh();
pMxDoc.UpdateContents();
}
public static void AddShapeFile()
{
//ClearAllData();
try
{
OpenFileDialog pOpenFileDialog = new OpenFileDialog
{
CheckFileExists = true,
Title = "打开Shape文件",
Filter = "Shape文件（*.shp）|*.shp"
};
pOpenFileDialog.ShowDialog();
////获取文件路径
//FileInfo pFileInfo = new FileInfo(pOpenFileDialog.FileName);
//string pPath = pOpenFileDialog.FileName.Substring(0, pOpenFileDialog.FileName.Length - pFileInfo.Name.Length);
//mainMapControl.AddShapeFile(pPath, pFileInfo.Name);
IWorkspaceFactory pWorkspaceFactory;
IFeatureWorkspace pFeatureWorkspace;
IFeatureLayer pFeatureLayer;
string pFullPath = pOpenFileDialog.FileName;
if (pFullPath == "")
{
return;
}
int pIndex = pFullPath.LastIndexOf("\\");
string pFilePath = pFullPath.Substring(0, pIndex); //文件路径
string pFileName = pFullPath.Substring(pIndex + 1); //文件名
//实例化ShapefileWorkspaceFactory工作空间，打开Shape文件
pWorkspaceFactory = new ShapefileWorkspaceFactory();
pFeatureWorkspace = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(pFilePath, 0);
//创建并实例化要素集
IFeatureClass pFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
pFeatureLayer = new FeatureLayer
{
FeatureClass = pFeatureClass
};
pFeatureLayer.Name = pFeatureLayer.FeatureClass.AliasName;
ClearAllData(); //新增删除数据
Form1.form1.axMapControl1.Map.AddLayer(pFeatureLayer);
Form1.form1.axMapControl1.ActiveView.Refresh();
//同步鹰眼
SynchronizeEagleEye();
}
catch (Exception ex)
{
MessageBox.Show("图层加载失败！" + ex.Message);
}
}
public IFeatureClass GetIFeatureClass(string FeatureType)
{
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkSpace = pWorkspaceFactory.OpenFromFile(@"res\poi_shp_jw", 0);
//初始化要素工作空间
IFeatureWorkspace pFeatWorkSpace = pWorkSpace as IFeatureWorkspace;
//IFeatureClass pFeatureClass = pFeatWorkSpace.OpenFeatureClass(FeatureType);
IFeatureClass pFeatureClass = pFeatWorkSpace.OpenFeatureClass("poi.shp");
return(pFeatureClass);
}
IFeatureClass LoadFeatureClassFromPath(string path)
{
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
string inputDir = "";
string inputName = "";
FileUtil.ExtractFileNameFromFullPath(path, ref inputDir, ref inputName);
IFeatureWorkspace pFeatureWorkspace = pWorkspaceFactory.OpenFromFile(inputDir, 0) as IFeatureWorkspace;
return(pFeatureWorkspace.OpenFeatureClass(inputName));
}
public static IWorkspace CreatOrGetSHPWorkspace(string strFilePath)
{
if (Directory.Exists(strFilePath) == false)
{
Directory.CreateDirectory(strFilePath);
}
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workspaceFactory.OpenFromFile(strFilePath, 0);
return(workspace);
}
/// <summary>
/// 生成signature文件
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
private void SC_CreateSampleFiles_ItemClick(object sender, DevExpress.XtraBars.ItemClickEventArgs e)
{
if (SC_dataGridView.Rows.Count == 0)
{
return;
}
SaveFileDialog SaveSignatureFile = new SaveFileDialog();
SaveSignatureFile.Title = "生成Signature文件";
SaveSignatureFile.Filter = "样本文件|*.gsg";
if (SaveSignatureFile.ShowDialog() == DialogResult.OK)
{
IGeoDataset inputraster = SampleLayerCombox.Tag as IGeoDataset;
//在临时文件夹生成featureclass，根据featureclass生成signature文件
//判断临时文件夹下是否有重名
int changefilename = 0;
while (System.IO.File.Exists(Application.StartupPath + "\\temp\\TempSample" + changefilename + ".shp"))
{
changefilename++;
}
//新建featureclass字段
IFields pFields = new FieldsClass();
IFieldsEdit pFieldsEdit = pFields as IFieldsEdit;
IField pField = new FieldClass();
IFieldEdit pFieldEdit = pField as IFieldEdit;
pFieldEdit.Name_2 = "Shape";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
//设置geometry definition
IGeometryDef pGeometryDef = new GeometryDefClass();
IGeometryDefEdit pGeometryDefEdit = pGeometryDef as IGeometryDefEdit;
pGeometryDefEdit.GeometryType_2 = ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPolygon;
pGeometryDefEdit.SpatialReference_2 = inputraster.SpatialReference;
pFieldEdit.GeometryDef_2 = pGeometryDef;
pFieldsEdit.AddField(pField);
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFeatureWorkspace = pWorkspaceFactory.OpenFromFile(Application.StartupPath + "\\temp", 0) as IFeatureWorkspace;
IFeatureClass featureclass = pFeatureWorkspace.CreateFeatureClass("TempSample" + changefilename + ".shp", pFields, null, null, esriFeatureType.esriFTSimple, "Shape", "");
//根据单元格tag中存储的多边形生成对应要素
for (int i = 0; i < SC_dataGridView.Rows.Count; i++)
{
IFeature feature = featureclass.CreateFeature();
feature.Shape = SC_dataGridView.Rows[i].Cells["color"].Tag as IPolygon;
feature.Store();
}
//生成signature文件
IGeoDataset Sampledataset = featureclass as IGeoDataset;
IMultivariateOp Multivariateop = new RasterMultivariateOpClass();
Multivariateop.CreateSignatures(inputraster, Sampledataset, SaveSignatureFile.FileName + ".gsg", true);
}
}
public void loadShapefile(string shpFileName)
{
// 文件格式判断
if (shpFileName.Substring(shpFileName.LastIndexOf(".")) != ".shp")
{
MessageBox.Show("请选择shp文件");
return;
}
// 文件打开
string filePath = "";
string fileName = "";
int index = shpFileName.LastIndexOf("\\"); // 14
Console.WriteLine(index);
if (index > 0)
{
filePath = shpFileName.Substring(0, index);
string temp = shpFileName.Substring(index + 1);
fileName = temp.Substring(0, temp.LastIndexOf("."));
}
IWorkspaceFactory pWorkspaceFactory;
IFeatureWorkspace pFeatureWorkspace;
IFeatureClass pFeatureClass;
IFeatureLayer pFeatureLayer;
pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(filePath, 0);
pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
// 将shp数据读入内存 --- >> pFeatureClass
pFeatureClass = pFeatureWorkspace.OpenFeatureClass(fileName); // 练习ArcMap
//pFeatureLayer = new FeatureLayer();
//pFeatureLayer.FeatureClass = pFeatureClass;
//pFeatureLayer.Name = "region";
//axMapControl1.AddLayer(pFeatureLayer);
//axMapControl1.Refresh();
// 用于图形渲染 ----> pFeatureLayer
pFeatureLayer = new FeatureLayer();
pFeatureLayer.FeatureClass = pFeatureClass;
IDataset pDataset = pFeatureClass as IDataset;
pFeatureLayer.Name = pDataset.Name; // set
//pFeatureLayer.Name = "afdsafdsafs"; // set
//string s = pFeatureLayer.Name; // get
ILayer pLayer = pFeatureLayer as ILayer;
this.mapControl.Map.AddLayer(pLayer);
this.mapControl.Refresh();
}
private void 加载图像ToolStripMenuItem_Click(object sender, EventArgs e)
{
//打开选择文件对话框用于选取图形文件
string fullFilePath; //存储打开文件的全路径
//设置OpenFileDialog的属性，使其能打开多种类型文件
OpenFileDialog openFile = new OpenFileDialog();
openFile.Multiselect = true;
openFile.Filter = "shape文件(*.shp)|*.shp";
openFile.Filter += "|栅格数据(*.jpg,*.bmp,*.tiff)|*.jpg;*.bmp;*.tiff;*.tif;*.img";
openFile.Title = "打开文件";
if (openFile.ShowDialog() == DialogResult.OK)
{
fullFilePath = openFile.FileName;
//获得文件路径
int index = fullFilePath.LastIndexOf("\\");
string filePath = fullFilePath.Substring(0, index);
//获得文件名称
string fileName = fullFilePath.Substring(index + 1);
//加载shape文件
if (openFile.FilterIndex == 1)
{
//打开工作空间工厂
IWorkspaceFactory workspcFac = new ShapefileWorkspaceFactory();
IFeatureWorkspace featureWorkspc;
IFeatureLayer featureLay = new FeatureLayerClass();
//打开路径
featureWorkspc = workspcFac.OpenFromFile(filePath, 0) as IFeatureWorkspace;
//打开类要素
featureLay.FeatureClass = featureWorkspc.OpenFeatureClass(fileName);
featureLay.Name = fileName;
//添加图层
axMapControl1.AddLayer(featureLay);
axMapControl1.Refresh();
}
//加载栅格图像
else if (openFile.FilterIndex == 2)
{
IWorkspaceFactory workspcFac = new RasterWorkspaceFactory();
IRasterWorkspace rasterWorkspc;
IRasterDataset rasterDatset = new RasterDatasetClass();
IRasterLayer rasterLay = new RasterLayerClass();
rasterWorkspc = workspcFac.OpenFromFile(filePath, 0) as IRasterWorkspace;
rasterDatset = rasterWorkspc.OpenRasterDataset(fileName);
rasterLay.CreateFromDataset(rasterDatset);
axMapControl1.ClearLayers();
rasterLay.Name = fileName;
axMapControl1.AddLayer(rasterLay);
axMapControl1.Refresh();
}
}
}
public static IFeatureClass GetShpFeatureClass(string directory, string fileName)
{
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workspaceFactory.OpenFromFile(directory, 0);
IFeatureClass featureClass = workspace.GetFeatureClass(fileName);
IWorkspaceFactoryLockControl factorylock = workspaceFactory as IWorkspaceFactoryLockControl;
if (factorylock.SchemaLockingEnabled)
{
factorylock.DisableSchemaLocking();
}
return(featureClass);
}
private void btnOpenSourceData_Click(object sender, EventArgs e)
{
OpenFileDlg.FileName = "";
OpenFileDlg.Filter = "GIS（*.shp）|*.shp";
OpenFileDlg.FilterIndex = 1;
OpenFileDlg.ShowDialog();
if (OpenFileDlg.FileName != "")
{
comboBoxInputData.Text = OpenFileDlg.FileName;
IWorkspaceFactory pWorkspaceFactory;
pWorkspaceFactory = new ShapefileWorkspaceFactory();
String wsp;
int r;
int l;
r = comboBoxInputData.Text.LastIndexOf("\\");
l = comboBoxInputData.Text.LastIndexOf(".");
wsp = comboBoxInputData.Text.Substring(0, r);
IFeatureWorkspace pFWS;
pFWS = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(wsp, 0);
try
{
String sp;
IFeatureClass pFeatureClass;
sp = comboBoxInputData.Text.Substring(r + 1, l - r - 1);
pFeatureClass = pFWS.OpenFeatureClass(sp);
comboBoxZValueField.Items.Clear();
comboBoxZValueField.Items.Add("无");
try
{
for (int j = 0; j < pFeatureClass.Fields.FieldCount; j++)
{
if (pFeatureClass.Fields.get_Field(j).Type == esriFieldType.esriFieldTypeDouble || pFeatureClass.Fields.get_Field(j).Type == esriFieldType.esriFieldTypeSingle || pFeatureClass.Fields.get_Field(j).Type == esriFieldType.esriFieldTypeSmallInteger || pFeatureClass.Fields.get_Field(j).Type == esriFieldType.esriFieldTypeInteger)
{
comboBoxZValueField.Items.Add(pFeatureClass.Fields.get_Field(j).Name);
}
}
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
catch (Exception ex)
{
MessageBox.Show(ex.ToString(), "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
bDataPath = true;
}
}
// 利用 文件路径 为 axMapControl1 添加图层
private void addLayer(string path)
{
// GetDirectoryName 方法可以解析出目录名
string directoryName = Path.GetDirectoryName(path);
// GetFileNameWithoutExtension 方法可以解析出
// 不包含扩展名的文件名称（不包含路径名）
string fileNameWithoutExtension = Path.GetFileNameWithoutExtension(path);
IWorkspaceFactory shapefileWorkspaceFactory = new ShapefileWorkspaceFactory();
// 调用 OpenFromFile 通过目录名方法打开 Shapefile 数据库
IWorkspace workspace = shapefileWorkspaceFactory.OpenFromFile(directoryName, 0);
if (workspace != null)
{
// ESRI.ArcGIS.Geodatabase.IEnumDataset
IEnumDataset enumDataset =
workspace.get_Datasets(esriDatasetType.esriDTFeatureClass);
enumDataset.Reset();
IDataset shpDataset = enumDataset.Next();
while (shpDataset != null)
{
if (shpDataset.Name == fileNameWithoutExtension)
{
// ESRI.ArcGIS.Carto.IFeatureLayer
IFeatureLayer newLayer = new FeatureLayerClass();
newLayer.FeatureClass = shpDataset as IFeatureClass;
newLayer.Name = fileNameWithoutExtension;
axMapControl1.Map.AddLayer(newLayer);
break;
}
shpDataset = enumDataset.Next();
}
}
}
protected override void OnMouseUp(MouseEventArgs arg)
{
base.OnMouseUp(arg);
IPoint pPoint;
pPoint = pMxDoc.CurrentLocation;
string strNewIdField;
strNewIdField = "NEWID";
string strRoadIdVal;
strRoadIdVal = GetFeatureVal(pPoint, pRoadIdxLayer, strNewIdField);
//** Reports road id # clicked. Remove this line if you do not want this functionality.
MessageBox.Show("Adding the road id #" + strRoadIdVal);
IFeatureLayer pFLayer;
pFLayer = new FeatureLayer();
IWorkspaceFactory pWSFactory;
pWSFactory = new ShapefileWorkspaceFactory();
//* Loads file
IWorkspace pWorkspace;
pWorkspace = pWSFactory.OpenFromFile("C:/Users/jcarpenter/Documents/PennState/Geog489/Lesson3/Data/roads", ArcMap.Application.hWnd);
IFeatureWorkspace pFWorkspace;
pFWorkspace = (IFeatureWorkspace)pWorkspace;
//* Concatenates file
IFeatureClass pFClass;
pFClass = pFWorkspace.OpenFeatureClass("roads_" + strRoadIdVal + ".shp");
pFLayer.FeatureClass = pFClass;
pFLayer.Name = "Road ID " + strRoadIdVal;
//* Adds the road layer
pMap.AddLayer(pFLayer);
//** Set style
SetStyle(pFLayer);
}
protected override void OnClick()
{
//
// TODO: Sample code showing how to access button host
//
var referenceFactory2 = (ISpatialReferenceFactory2)new SpatialReferenceEnvironment();
ISpatialReference WGS84 = referenceFactory2.CreateSpatialReference((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
IPointCollection pointCollection = new PolylineClass();
// ------------- Ensure that both points have negative longitude values -------------------
IPoint point = new PointClass();
point.PutCoords(-170, 10); // Equivalent to 170 degrees WEST
point.SpatialReference = WGS84;
pointCollection.AddPoint(point);
point = new PointClass();
point.PutCoords(-200, 10); // Equivalent to 160 degrees EAST
point.SpatialReference = WGS84;
pointCollection.AddPoint(point);
// -----------------------------------------------------------------------
IPolyline polyline = (IPolyline)pointCollection;
polyline.SpatialReference = WGS84;
var geometryDefEdit = (IGeometryDefEdit)new GeometryDef();
geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolyline;
geometryDefEdit.SpatialReference_2 = WGS84;
var field = (IFieldEdit)new Field();
field.Name_2 = "Shape";
field.Type_2 = esriFieldType.esriFieldTypeGeometry;
field.GeometryDef_2 = geometryDefEdit;
var fields = (IFieldsEdit)new Fields();
fields.AddField(field);
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
var featureWorkspace = (IFeatureWorkspace)workspaceFactory.OpenFromFile("C:\\Temp\\", 0);
IFeatureClass featureClass = featureWorkspace.CreateFeatureClass
("test3.shp", fields, null, null, esriFeatureType.esriFTSimple, "Shape", "");
var feature = featureClass.CreateFeature();
feature.Shape = polyline;
feature.Store();
ArcMap.Application.CurrentTool = null;
}
private void GetUniqueValueList()
{
//get all the unique value of the classitem
if ((m_pEditLayer.classitem == null) || (m_pEditLayer.classitem.Length == 0))
return;
//label
string strDataPath = Utility.GetAbsolutePathOfData(m_pEditLayer);
IWorkspaceFactory pWsf = new ShapefileWorkspaceFactory();
IWorkspace pWs = pWsf.OpenFromFile(strDataPath, 0);
if (pWs != null)
{
IFeatureWorkspace pFws = pWs as IFeatureWorkspace;
string strDataName = System.IO.Path.GetFileNameWithoutExtension(strDataPath);
IFeatureClass pFeatureClass = pFws.OpenFeatureClass(strDataName);
if (pFeatureClass != null)
{
ICursor pCursor = pFeatureClass.Search(null, false) as ICursor;
IDataStatistics pDataStatistics = new DataStatisticsClass();
pDataStatistics.Field = m_pEditLayer.classitem;
pDataStatistics.Cursor = pCursor;
m_pUniqueEnumerator = pDataStatistics.UniqueValues;
//pDataStatistics
}
}
}
private void Create_Click(object sender, EventArgs e)
{
try
{
IWorkspaceFactory pShpWksFact = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFeatWks;
pFeatWks = (IFeatureWorkspace)pShpWksFact.OpenFromFile(filePath, 0);
const string strShapeFieldName = "Shape";
//定义属性字段
IFields pFields = new Fields();
IFieldsEdit pFieldsEdit;
pFieldsEdit = pFields as IFieldsEdit;
IField pField = new Field();
IFieldEdit pFieldEdit = new Field() as IFieldEdit;
pFieldEdit.Name_2 = strShapeFieldName;
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
pField = pFieldEdit as IField;
//定义几何属性
IGeometryDef pGeomDef = new GeometryDef();
IGeometryDefEdit pGeomDefEdit = new GeometryDef() as IGeometryDefEdit;
pGeomDefEdit = pGeomDef as IGeometryDefEdit;
switch (shpTypeComboBox.Text)
{
case "Point":
pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
break;
case "Polyline":
pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolyline;
break;
case "Polygon":
pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;
break;
}
pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
pGeomDefEdit.SpatialReference_2 = new UnknownCoordinateSystem() as ISpatialReference;
pFieldEdit.GeometryDef_2 = pGeomDef;
pFieldsEdit.AddField(pField);
pFields = pFieldsEdit as IFields;
IFeatureClass pFeatureClass;
pFeatureClass = pFeatWks.CreateFeatureClass(fileName, pFields, null, null,
esriFeatureType.esriFTSimple, strShapeFieldName, "");
//添加属性字段
for (int i = 0; i < addFieldListBox.Items.Count; i++)
{
IField pfield = new Field();
IFieldEdit pfieldEdit = new Field() as IFieldEdit;
pfieldEdit.Name_2 = addFieldListBox.Items[i].ToString();
pfieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pfield = pfieldEdit as IField;
pFeatureClass.AddField(pfield);
}
//绘制点
for (int i = 0; i < excelDataGridViewX.Rows.Count - 1; i++)
{
DataGridViewRow dataRow = excelDataGridViewX.Rows[i];
double pointX, pointY;
pointX = double.Parse(dataRow.Cells[xComboBoxEx.Text].Value.ToString());
pointY = double.Parse(dataRow.Cells[yComboBoxEx.Text].Value.ToString());
IPoint pPoint = new ESRI.ArcGIS.Geometry.Point() as IPoint;
pPoint.PutCoords(pointX, pointY);
IFeature pFeature = pFeatureClass.CreateFeature();
pFeature.Shape = pPoint;
//为该点添加属性值
for (int j = 0; j < addFieldListBox.Items.Count; j++)
{
string fieldName = addFieldListBox.Items[j].ToString();
pFeature.set_Value(pFeature.Fields.FindField(fieldName),
dataRow.Cells[fieldName].Value.ToString());
}
pFeature.Store();
}
//添加新建的数据至Map中
axMapControl.AddShapeFile(filePath, fileName);
this.Hide();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
//�����������Ч����¥�ϵ�1/50������
//public static string GetAbsolutePath(string dir, string path)
//{
//string strDir = Directory.GetCurrentDirectory();
//Directory.SetCurrentDirectory(dir);
//string strNew = Path.GetFullPath(path);
//Directory.SetCurrentDirectory(strDir);
//return strNew;
//}
/// <summary>
/// get all fields' name of the layer
/// </summary>
/// <param name="layer"></param>
/// <returns></returns>
public static StringCollection GetFieldsOfLayer(layerObj layer)
{
StringCollection pStringCollection = new StringCollection();
string strDataPath = GetAbsolutePathOfData(layer);
IWorkspaceFactory pWsf = new ShapefileWorkspaceFactory();
IWorkspace pWs = pWsf.OpenFromFile(Path.GetDirectoryName(strDataPath), 0);
if (pWs != null)
{
IFeatureWorkspace pFws = pWs as IFeatureWorkspace;
string strDataName = Path.GetFileNameWithoutExtension(strDataPath);
IFeatureClass pFeatureClass = pFws.OpenFeatureClass(strDataName);
if (pFeatureClass != null)
{
IFields pFields = pFeatureClass.Fields;
if(pFields != null)
{
int iCount = pFields.FiledCount;
for (int i = 0; i < iCount; i++)
{
pStringCollection.Add(pFields.GetField(i).Name);
}
}
}
}
return pStringCollection;
}
public static ITable ouvrirDbf(string sCheminDbf)
{
IWorkspaceFactory pFact;
IWorkspace pWorkspace;
IFeatureWorkspace pFeatws;
ITable pTable;
pFact = new ShapefileWorkspaceFactory();
pWorkspace = pFact.OpenFromFile(System.IO.Path.GetDirectoryName(sCheminDbf), 0);
pFeatws = (IFeatureWorkspace)pWorkspace;
pTable = pFeatws.OpenTable(System.IO.Path.GetFileName(sCheminDbf));
return pTable;
}
/// <summary>
/// 连接Shp的Workspace
/// </summary>
/// <param name="Path">Shp路径</param>
public IWorkspace Workspace_SHP(string Path)
{
IWorkspaceFactory WorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace WS = null;
try
{
WS = WorkspaceFactory.OpenFromFile(Path, 0);
}
catch
{
MessageBox.Show("请检查Shp路径是否正确");
}
return WS;
}
public void DeleteShapeFile(string shapeFilePath)
{
string fcName = System.IO.Path.GetFileName(shapeFilePath);
string folderName = System.IO.Path.GetDirectoryName(shapeFilePath);
using (ComReleaser oComReleaser = new ComReleaser())
{
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workspaceFactory.OpenFromFile(folderName, 0);
IFeatureWorkspace fWorkspace = (IFeatureWorkspace)workspace;
IDataset ipDs = fWorkspace.OpenFeatureClass(fcName) as IDataset;
ipDs.Delete();
File.Delete(shapeFilePath);
System.Runtime.InteropServices.Marshal.FinalReleaseComObject(workspace);
workspace = null;
fWorkspace = null;
ipDs = null;
}
GC.Collect();
}
//
//��������
//
public static IWorkspace OpenWorkspace(string strWorkspaceName)
{
IWorkspaceFactory workspaceFactory;
workspaceFactory = new ShapefileWorkspaceFactory();
return workspaceFactory.OpenFromFile(strWorkspaceName, 0);
}
public static void addShapeFileToLayer(string sPath, String sShapeName, bool isVisible = false)
{
IWorkspaceFactory pWorkspaceFactory;
IFeatureWorkspace pFeatureWorkspace;
IFeatureLayer pFeaturelayer;
try
{
// Create a new ShapefileWorkspaceFactory object and open a shapefile folder
pWorkspaceFactory = new ShapefileWorkspaceFactory();
pFeatureWorkspace = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(sPath, 0);
// Create a new FeatureLayer and assign a shapefile to it
pFeaturelayer = new FeatureLayer();
pFeaturelayer.FeatureClass = pFeatureWorkspace.OpenFeatureClass(sShapeName);
pFeaturelayer.Name = pFeaturelayer.FeatureClass.AliasName;
pFeaturelayer.Visible = isVisible;
// Add the FeatureLayer to the focus map
SF10_extMain.m_pDoc.AddLayer(pFeaturelayer);
SF10_extMain.m_pDoc.UpdateContents();
}
catch
{
}
pWorkspaceFactory = null;
pFeatureWorkspace = null;
pFeaturelayer = null;
}
public static IFeatureClass createNewEmptyShapefilePoint(string sP, string sN, int sPj)
{
String strFolder;
String strName;
const string strShapeFieldName = "Shape";
try
{
strFolder = sP;
strName = sN;
// Open the folder to contain the shapefile as a workspace
IFeatureWorkspace pFWS;
IWorkspaceFactory pWorkspaceFactory;
pWorkspaceFactory = new ShapefileWorkspaceFactory();
pFWS = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(strFolder, 0);
// Set up a simple fields collection
IFields pFields;
IFieldsEdit pFieldsEdit;
pFields = new Fields();
pFieldsEdit = (IFieldsEdit)pFields;
IField pField;
IFieldEdit pFieldEdit;
// Make the shape field
// it will need a geometry definition, with a spatial reference
pField = new Field();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = strShapeFieldName;
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
IGeometryDef pGeomDef;
IGeometryDefEdit pGeomDefEdit;
pGeomDef = new GeometryDef();
pGeomDefEdit = (IGeometryDefEdit)pGeomDef;
pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint; //Creates point shapefile
if (sPj == 0)
{
pGeomDefEdit.SpatialReference_2 = new UnknownCoordinateSystemClass();
}
else
{
ISpatialReferenceFactory2 pSpatRefFact;
pSpatRefFact = new SpatialReferenceEnvironment() as ISpatialReferenceFactory2;
IProjectedCoordinateSystem pGeoCoordSys;
pGeoCoordSys = pSpatRefFact.CreateProjectedCoordinateSystem(sPj);
pGeomDefEdit.SpatialReference_2 = pGeoCoordSys;
}
pFieldEdit.GeometryDef_2 = pGeomDef;
pFieldsEdit.AddField(pField);
// Add another miscellaneous text field
pField = new Field();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Length_2 = 25;
pFieldEdit.Name_2 = "TextField";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
pFieldsEdit.AddField(pField);
// Create the shapefile
return pFWS.CreateFeatureClass(strName, pFields, null, null, esriFeatureType.esriFTSimple, strShapeFieldName, "");
}
catch
{ return null; }
}
public void CreateRasterFromPoints(IMap pMap, IFeatureLayer pp,String Path, String Name, String Field)
{
//1.��Shape�ļ���ȡ��FeatureClass
//2.����FeatureClass����IFeatureClassDescriptor
//3.����IRasterRaduis ����
//����Cell
//4.��ֵ�����ɱ���
object obj = null;
IWorkspaceFactory pShapSpace;
pShapSpace = new ShapefileWorkspaceFactory();
IWorkspace pW;
pW = pShapSpace.OpenFromFile(Path, 0);
IFeatureWorkspace pFeatureWork;
pFeatureWork = pW as IFeatureWorkspace;
IFeatureClass pFeatureClass = pFeatureWork.OpenFeatureClass("Name");
IGeoDataset Geo = pFeatureClass as IGeoDataset;
object extend = Geo.Extent;
object o = null;
IFeatureClassDescriptor pFeatureClassDes = new FeatureClassDescriptorClass();
pFeatureClassDes.Create(pFeatureClass, null, Field);
IRasterRadius pRasterrad = new RasterRadiusClass();
pRasterrad.SetVariable(10, ref obj);
object dCell = 0.5;//���Ը��ݲ�ͬ�ĵ�ͼ���������
IInterpolationOp Pinterpla = new RasterInterpolationOpClass();
IRasterAnalysisEnvironment pRasterAnaEn = Pinterpla as IRasterAnalysisEnvironment;
pRasterAnaEn.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCell);
pRasterAnaEn.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extend, ref o);
IGeoDataset pGRaster;
object hh = 3;
pGRaster = Pinterpla.IDW((IGeoDataset)pFeatureClassDes , 2, pRasterrad, ref hh);
ISurfaceOp pSurF;
pSurF = new RasterSurfaceOpClass();
IGeoDataset pCour;
object o1 = 0;
pCour = pSurF.Contour(pGRaster, 5, ref o1);
IFeatureLayer pLa;
pLa = new FeatureLayerClass();
IFeatureClass pClass = pCour as IFeatureClass;
pLa.FeatureClass = pClass;
pLa.Name = pClass.AliasName;
pMap.AddLayer(pLa as ILayer);
}
private void AddTable(string strPath, string strFileName)
{
try
{
IWorkspaceFactory pFact = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace = pFact.OpenFromFile(strPath, 0);
IFeatureWorkspace pFeatws = pWorkspace as IFeatureWorkspace;
ITable ptable = pFeatws.OpenTable(strFileName);
IMxDocument pmxdoc = ArcMap.Document as IMxDocument;
IMap pmap = pmxdoc.FocusMap;
IStandaloneTable pStTab = new StandaloneTableClass();
pStTab.Table = ptable;
IStandaloneTableCollection pStTabColl = pmap as IStandaloneTableCollection;
pStTabColl.AddStandaloneTable(pStTab);
pmxdoc.UpdateContents();
ITableWindow2 ptableWindow = new TableWindowClass();
ptableWindow.Application = ArcMap.Application;
ptableWindow.StandaloneTable = pStTab;
ptableWindow.Show(true);
System.Runtime.InteropServices.Marshal.ReleaseComObject(pFact);
System.Runtime.InteropServices.Marshal.ReleaseComObject(pWorkspace);
System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatws);
System.Runtime.InteropServices.Marshal.ReleaseComObject(ptable);
System.Runtime.InteropServices.Marshal.ReleaseComObject(pStTab);
System.Runtime.InteropServices.Marshal.ReleaseComObject(pStTabColl);
System.Runtime.InteropServices.Marshal.ReleaseComObject(ptableWindow);
GC.Collect();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
public static void Mlayer_Krige_Click()
{
// 用克里金Krige插值生成的栅格图像。如下：
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
string pPath = Application.StartupPath + @"\MakeContours\Cont.shp";
string pFolder = System.IO.Path.GetDirectoryName(pPath);
string pFileName = System.IO.Path.GetFileName(pPath);
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pFolder, 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass oFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
IFeatureClassDescriptor pFCDescriptor = new FeatureClassDescriptorClass();
pFCDescriptor.Create(oFeatureClass, null, "shape.z");
IInterpolationOp pInterpolationOp = new RasterInterpolationOpClass();
IRasterAnalysisEnvironment pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
object Cellsize = 0.004;//Cell size for output raster;0.004
pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref Cellsize);
//设置输出范围
//27202 19104;27202 22947;31550 22947;31550 19104
object snapRasterData = Type.Missing;
IEnvelope pExtent;
pExtent = new EnvelopeClass();
Double xmin = 27202;
Double xmax = 31550;
Double ymin = 19104;
Double ymax = 22947;
pExtent.PutCoords(xmin, ymin, xmax, ymax);
object extentProvider = pExtent;
pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extentProvider, ref snapRasterData);
Double dSearchD = 10;
object pSearchCount = 3;
object missing = Type.Missing;
IRasterRadius pRadius = new RasterRadius();
pRadius.SetFixed(dSearchD, ref pSearchCount);
//pRadius.SetVariable((int)pSearchCount, ref dSearchD);
IGeoDataset poutGeoDataset = pInterpolationOp.Krige((IGeoDataset)pFCDescriptor, esriGeoAnalysisSemiVariogramEnum.esriGeoAnalysisGaussianSemiVariogram, pRadius, false, ref missing);
IRaster pOutRaster = poutGeoDataset as IRaster;
IRasterLayer pOutRasLayer = new RasterLayer();
pOutRasLayer.CreateFromRaster(pOutRaster);
IMap pMap = Common.DataEditCommon.g_pMap;
pMap.AddLayer(pOutRasLayer);
Common.DataEditCommon.g_axTocControl.Refresh();
Common.DataEditCommon.g_pAxMapControl.ActiveView.Refresh();
}
public static void Mlayer_IDW_Click()
{
// 用反距离IDW插值生成的栅格图像。如下：
IInterpolationOp pInterpolationOp = new RasterInterpolationOpClass();
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
string pPath = Application.StartupPath + @"\MakeContours\Cont.shp";
string pFolder = System.IO.Path.GetDirectoryName(pPath);
string pFileName = System.IO.Path.GetFileName(pPath);
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(pFolder, 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass oFeatureClass = pFeatureWorkspace.OpenFeatureClass(pFileName);
IFeatureClassDescriptor pFCDescriptor = new FeatureClassDescriptorClass();
pFCDescriptor.Create(oFeatureClass, null, "shape.z");
IRasterRadius pRadius = new RasterRadiusClass();
object objectMaxDistance = null;
object objectbarrier = null;
object missing = Type.Missing;
pRadius.SetVariable(12, ref objectMaxDistance);
object dCellSize =1;
object snapRasterData = Type.Missing;
IEnvelope pExtent;
pExtent = new EnvelopeClass();
Double xmin = 27202;
Double xmax = 31550;
Double ymin = 19104;
Double ymax = 22947;
pExtent.PutCoords(xmin, ymin, xmax, ymax);
object extentProvider = pExtent;
IRasterAnalysisEnvironment pEnv = pInterpolationOp as IRasterAnalysisEnvironment;
pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCellSize);
pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, ref extentProvider, ref snapRasterData);
IGeoDataset poutGeoDataset = pInterpolationOp.IDW((IGeoDataset)pFCDescriptor, 2, pRadius, ref objectbarrier);
ISurfaceOp surOp = new RasterSurfaceOpClass();
IRaster pOutRaster = poutGeoDataset as IRaster;
IRasterLayer pOutRasLayer = new RasterLayer();
pOutRasLayer.CreateFromRaster(pOutRaster);
IMap pMap = Common.DataEditCommon.g_pMap;
pMap.AddLayer(pOutRasLayer);
Common.DataEditCommon.g_axTocControl.Refresh();
Common.DataEditCommon.g_pAxMapControl.ActiveView.Refresh();
}
/// <summary>
/// Add Shapefile layer
/// </summary>
/// <param name="path"></param>
/// <param name="fileName"></param>
private void AddShapeFileLayer(string path, string fileName)
{
try
{
// Create a new ShapefileWorkspaceFactory object and
// open a shapefile folder - the path works with standard 9.3 installation
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)
workspaceFactory.OpenFromFile(path, 0);
// Create a new FeatureLayer and assign a shapefile to it
IFeatureLayer featureLayer = new ESRI.ArcGIS.Carto.FeatureLayer();
featureLayer.FeatureClass = featureWorkspace.OpenFeatureClass(fileName);
ILayer layer = (ILayer)featureLayer;
layer.Name = featureLayer.FeatureClass.AliasName;
// Add the Layer to the focus map
IMap map = ArcMap.Document.FocusMap;
map.AddLayer(layer);
}
catch (Exception e)
{
Log.WriteLine("AddShapeFileLayer: " + e.Message);
}
}
/// <summary>
/// 通过TIN文件，建立*.shp文件
/// </summary>
/// <param name="TINDir">TIN文件的文件夹路径</param>
private void CreateContourData(string TINDir)
{
try
{
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
CopyFolder(Application.StartupPath + "\\Convert\\TemContour", Application.StartupPath + "\\Convert\\ContourResult\\" + this.Random);
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(Application.StartupPath + "\\Convert\\ContourResult\\" + this.Random, 0);
IWorkspaceEdit w = pWorkspace as IWorkspaceEdit;
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass tFeatureClass = pFeatureWorkspace.OpenFeatureClass("Contour.shp");
ITinWorkspace pTinWorkspace;
IWorkspace pWS;
IWorkspaceFactory pWSFact = new TinWorkspaceFactoryClass();
pWS = pWSFact.OpenFromFile(TINDir, 0);
ITin pTin;
pTinWorkspace = pWS as ITinWorkspace;
w.StartEditing(true);
w.StartEditOperation();
if (pTinWorkspace.get_IsTin(this.Random))
{
pTin = pTinWorkspace.OpenTin(this.Random);
ITinSurface pTinSurface = pTin as ITinSurface;
pTinSurface.Contour(0, this.Interval, tFeatureClass, "Z", 0);
w.StopEditOperation();
w.StopEditing(true);
}
}
catch (Exception)
{
}
}
/// <summary>
/// 创建等高线*.shp文件的主方法
/// </summary>
/// <param name="lux">左上投影x坐标</param>
/// <param name="luy">左上投影y坐标</param>
/// <param name="luz">左上投影z坐标</param>
/// <param name="rlx">右下投影x坐标</param>
/// <param name="rly">右下投影y坐标</param>
/// <param name="rlz">右下投影z坐标</param>
/// <param name="m_interval">采样间隔</param>
/// <returns>等高线*.shp文件的文件名</returns>
public string CreateContourShape(ISGWorld61 sgworld, double lux, double luy, double luz, double rlx, double rly, double rlz,double m_interval)
{
double width = sgworld.CoordServices.GetDistance(lux, luy, rlx, luy);
double hight = sgworld.CoordServices.GetDistance(lux, luy, lux, rly);
if (width<5||hight<10)
{
MessageBox.Show("范围过小！","SUNZ",MessageBoxButtons.OK,MessageBoxIcon.Warning);
return null;
}
if ((width/1000)*(hight/1000)>40)
{
if (MessageBox.Show("范围超出40平分公里,系统计算时间比较长，是否继续？","提示",MessageBoxButtons.OKCancel,MessageBoxIcon.Question) == DialogResult.OK)
{
}
else
{
return null;
}
}
this.Interval = m_interval;
this.Random = System.Guid.NewGuid().ToString().Substring(0, 6).ToLower();
CopyFolder(Application.StartupPath + "\\Convert\\TemPoints", Application.StartupPath + "\\Convert\\PointsResult\\" + this.Random);
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(Application.StartupPath + "\\Convert\\PointsResult\\" + this.Random, 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass pFeatureClass = pFeatureWorkspace.OpenFeatureClass(System.IO.Path.GetFileNameWithoutExtension(Application.StartupPath + "\\Convert\\PointsResult\\" + Random + "\\TemPointLayer.shp"));
IFeatureClassWrite fr = (IFeatureClassWrite)pFeatureClass;
IWorkspaceEdit w = (pFeatureClass as IDataset).Workspace as IWorkspaceEdit;
IFeature f;
w.StartEditing(true);
w.StartEditOperation();
int EastPointCount = (int)(width / 5);
int NothPointCount = (int)(hight / 5);
for (int j = 0; j < NothPointCount; j++)
{
ESRI.ArcGIS.Geometry.IPoint p;
p = new ESRI.ArcGIS.Geometry.PointClass();
ESRI.ArcGIS.Geometry.IZAware iz = p as ESRI.ArcGIS.Geometry.IZAware;
iz.ZAware = true;
p.X = lux;
p.Y = luy;
p.Z = luz;
ESRI.ArcGIS.Geometry.IGeometry peo;
peo = p;
f = pFeatureClass.CreateFeature();
f.Shape = peo;
f.set_Value(3, p.Z);
f.Store();
fr.WriteFeature(f);
for (int i = 0; i < EastPointCount; i++)
{
IPosition61 positionLU = null;
ICoord2D pPoint = null;
if ((j + 2) % 2 == 0)
{
pPoint = sgworld.CoordServices.MoveCoord(lux, luy, 5, 0);
}
else
{
pPoint = sgworld.CoordServices.MoveCoord(lux, luy, -5, 0);
}
IWorldPointInfo61 pW = sgworld.Terrain.GetGroundHeightInfo(pPoint.X, pPoint.Y, AccuracyLevel.ACCURACY_BEST_FROM_MEMORY, true);
positionLU = pW.Position;
luz = positionLU.Altitude;
lux = positionLU.X;
luy = positionLU.Y;
p = new ESRI.ArcGIS.Geometry.PointClass();
iz = p as ESRI.ArcGIS.Geometry.IZAware;
iz.ZAware = true;
p.X = lux;
p.Y = luy;
p.Z = luz;
peo = p;
f = pFeatureClass.CreateFeature();
f.Shape = peo;
f.set_Value(3, p.Z);
f.Store();
fr.WriteFeature(f);
}
ICoord2D pPointL = sgworld.CoordServices.MoveCoord(lux, luy, 0, -5);
IWorldPointInfo61 pWL = sgworld.Terrain.GetGroundHeightInfo(pPointL.X, pPointL.Y, AccuracyLevel.ACCURACY_FORCE_BEST_RENDERED, true);
luz = pWL.Position.Altitude;
lux = pWL.Position.X;
luy = pWL.Position.Y;
}
w.StopEditOperation();
w.StopEditing(true);
CreateTinFromFeature(pFeatureClass);
return this.Random;
}
private void setupmap()
{
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workspaceFactory.OpenFromFile(@"C:\Users\huty\Desktop\武汉实习数据", 0);
IFeatureWorkspace featureworkspace = (IFeatureWorkspace)workspace;
mainroad = featureworkspace.OpenFeatureClass("主干道.shp");
mainroadlayer = new FeatureLayerClass();
mainroadlayer.FeatureClass = mainroad;
mainroadlayer.Name = "mainroad";
secondroad = featureworkspace.OpenFeatureClass("次干道.shp");
secondroadlayer = new FeatureLayerClass();
secondroadlayer.FeatureClass = secondroad;
secondroadlayer.Name = "secondroad";
waterregion = featureworkspace.OpenFeatureClass("水域.shp");
waterregionlayer = new FeatureLayerClass();
waterregionlayer.FeatureClass = waterregion;
waterregionlayer.Name = "waterregion";
supermarket = featureworkspace.OpenFeatureClass("超市.shp");
supermarketlayer = new FeatureLayerClass();
supermarketlayer.FeatureClass = supermarket;
supermarketlayer.Name = "supermarket";
address = featureworkspace.OpenFeatureClass("地名.shp");
addresslayer = new FeatureLayerClass();
addresslayer.FeatureClass = address;
addresslayer.Name = "address";
hotel = featureworkspace.OpenFeatureClass("饭店.shp");
hotellayer = new FeatureLayerClass();
hotellayer.FeatureClass = hotel;
hotellayer.Name = "hotel";
trainstation = featureworkspace.OpenFeatureClass("火车站.shp");
trainstationlayer = new FeatureLayerClass();
trainstationlayer.FeatureClass = trainstation;
trainstationlayer.Name = "trainstation";
bank = featureworkspace.OpenFeatureClass("金融银行.shp");
banklayer = new FeatureLayerClass();
banklayer.FeatureClass = bank;
banklayer.Name = "bank";
sight = featureworkspace.OpenFeatureClass("景点.shp");
sightlayer = new FeatureLayerClass();
sightlayer.FeatureClass = sight;
sightlayer.Name = "sight";
busstation = featureworkspace.OpenFeatureClass("汽车站.shp");
busstationlayer = new FeatureLayerClass();
busstationlayer.FeatureClass = busstation;
busstationlayer.Name = "busstation";
government = featureworkspace.OpenFeatureClass("区政府.shp");
governmentlayer = new FeatureLayerClass();
governmentlayer.FeatureClass = government;
governmentlayer.Name = "government";
hospital = featureworkspace.OpenFeatureClass("医院.shp");
hospitallayer = new FeatureLayerClass();
hospitallayer.FeatureClass = hospital;
hospitallayer.Name = "hospital";
recreation = featureworkspace.OpenFeatureClass("娱乐场所.shp");
recreationlayer = new FeatureLayerClass();
recreationlayer.FeatureClass = recreation;
recreationlayer.Name = "recreation";
railway = featureworkspace.OpenFeatureClass("铁路.shp");
railwaylayer = new FeatureLayerClass();
railwaylayer.FeatureClass = railway;
railwaylayer.Name = "railway";
}
/// <summary>
/// 叠置分析
/// </summary>
/// <param name="axmapcontrol"></param>
/// <param name="i">图层索引</param>
/// <param name="name">输出图层名字</param>
/// <returns></returns>
public static FeatureLayer Stastic(AxMapControl axmapcontrol,int i,string name)
{
if (axmapcontrol == null || axmapcontrol.LayerCount < 2)
{
return null;
}
// pLayer = axmapcontrol.Map.Layer(0);
ILayer pLayer = axmapcontrol.ActiveView.FocusMap.get_Layer(i);
IFeatureLayer pFeatureLayer = pLayer as IFeatureLayer;
IFeatureClass pFeatureClass = pFeatureLayer.FeatureClass;
ITable table = pFeatureClass as ITable;
pLayer = axmapcontrol.ActiveView.FocusMap.get_Layer(axmapcontrol.LayerCount - 1);
IFeatureLayer ppFeatureLayer=(IFeatureLayer)pLayer;
IFeatureClass ppFeatureClass=ppFeatureLayer.FeatureClass;
ITable ptable = (ITable)ppFeatureClass;
//Intersect intersect=new Intersect()
Double tol=0.0;
IWorkspaceFactory wspfctry = new ShapefileWorkspaceFactory();
IWorkspace pWorkSpace = wspfctry.OpenFromFile(SystemSet.Base_Map + "\\处理数据库\\实验数据", 0);
IDataset pDataset = (IDataset)pWorkSpace;
IName pName=pDataset.FullName;
IDatasetName pDatasetName = new FeatureClassNameClass();
pDatasetName.Name = name;
pDatasetName.WorkspaceName = (IWorkspaceName)pName;
IFeatureClassName fcn = (IFeatureClassName)pDatasetName;
IBasicGeoprocessor pBG=new BasicGeoprocessor();
IFeatureClass pOutFeatureClass = pBG.Intersect(table, false, ptable, false, tol, fcn);
IFeatureLayer pOutFeatLayer = new FeatureLayerClass();
FeatureLayer featurelayer = (FeatureLayer)pOutFeatLayer;
featurelayer.FeatureClass = pOutFeatureClass;
featurelayer.Name = pOutFeatureClass.AliasName;
return featurelayer;
}
//���һ����ͼ�㵽mapfile
public bool addLayer(string dataPath)
{
if (_Map == null)
{
if (ErrorOccured != null)
ErrorOccured(null, new MSErrorEventArgs("û�д�mapfile"));
return false;
}
layerObj pLayer = new layerObj(_Map);
pLayer.status = 1;
string strName = Path.GetFileNameWithoutExtension(dataPath);
pLayer.name = strName;
pLayer.group = strName;
string strExt = Path.GetExtension(dataPath).ToLower();
//wanliyun:2009.05.06
//add the mapinfo tab format support
if ((strExt == ".shp") || (strExt == ".tab"))
{
IFeatureClass pFc = null;
if (strExt == ".shp")
{
IWorkspaceFactory pWsf = new ShapefileWorkspaceFactory();
IFeatureWorkspace pFws = pWsf.OpenFromFile(dataPath, 0) as IFeatureWorkspace;
pFc = pFws.OpenFeatureClass(strName);
pLayer.connectiontype = MS_CONNECTION_TYPE.MS_SHAPEFILE;
pLayer.data = dataPath;
}
else if (strExt == ".tab")
{
IWorkspaceFactory pWsf = new OgrWorkspaceFactory();
IFeatureWorkspace pFws = pWsf.OpenFromFile(dataPath, 0) as IFeatureWorkspace;
pFc = pFws.OpenFeatureClass(strName);
pLayer.connectiontype = MS_CONNECTION_TYPE.MS_OGR;
pLayer.connection = dataPath;
}
IEnvelope pEnvelope = ((IGeoDataset)pFc).Extent;
if (_Map.numlayers == 1)
{
_Map.extent.minx = pEnvelope.XMin;
_Map.extent.miny = pEnvelope.YMin;
_Map.extent.maxx = pEnvelope.XMax;
_Map.extent.maxy = pEnvelope.YMax;
}
else
{
//if the layer extents lager than the map extents
//reset the map extents
if (pEnvelope.XMin < _Map.extent.minx)
_Map.extent.minx = pEnvelope.XMin;
if (pEnvelope.YMin < _Map.extent.miny)
_Map.extent.miny = pEnvelope.YMin;
if (pEnvelope.XMax > _Map.extent.maxx)
_Map.extent.maxx = pEnvelope.XMax;
if (pEnvelope.YMax > _Map.extent.maxy)
_Map.extent.maxy = pEnvelope.YMax;
}
//rerecord the max extents
m_MapMinx = _Map.extent.minx;
m_MapMiny = _Map.extent.miny;
m_MapMaxx = _Map.extent.maxx;
m_MapMaxy = _Map.extent.maxy;
classObj myClass = new classObj(pLayer);
myClass.name = "random";
styleObj myStyle = new styleObj(myClass);
Color pRandomColor = GetRandomColor();
myStyle.color.red = pRandomColor.R;
myStyle.color.green = pRandomColor.G;
myStyle.color.blue = pRandomColor.B;
myStyle.opacity = 0;
fmapGeometryType pGeomType = pFc.ShapeType;
if ((pGeomType == fmapGeometryType.Point) ||
(pGeomType == fmapGeometryType.Point25D) ||
(pGeomType == fmapGeometryType.MultiPoint) ||
(pGeomType == fmapGeometryType.MultiPoint25D))
{
pLayer.type = MS_LAYER_TYPE.MS_LAYER_POINT;
}
else if ((pGeomType == fmapGeometryType.LineString) ||
(pGeomType == fmapGeometryType.LineString) ||
(pGeomType == fmapGeometryType.MultiLineString) ||
(pGeomType == fmapGeometryType.MultiLineString25D))
{
pLayer.type = MS_LAYER_TYPE.MS_LAYER_LINE;
}
else if ((pGeomType == fmapGeometryType.Polygon) ||
(pGeomType == fmapGeometryType.Polygon25D) ||
(pGeomType == fmapGeometryType.MultiPolygon) ||
(pGeomType == fmapGeometryType.MultiPolygon25D))
{
pRandomColor = GetRandomColor();
myStyle.outlinecolor.red = pRandomColor.R;
myStyle.outlinecolor.green = pRandomColor.G;
myStyle.outlinecolor.blue = pRandomColor.B;
pLayer.type = MS_LAYER_TYPE.MS_LAYER_POLYGON;
}
else
{
pLayer.type = MS_LAYER_TYPE.MS_LAYER_POLYGON;
}
if (pFc != null)
pFc.Dispose();
}
else if ((strExt == ".tif") || (strExt == ".img"))
{
pLayer.data = dataPath;
pLayer.connectiontype = MS_CONNECTION_TYPE.MS_RASTER;
pLayer.type = MS_LAYER_TYPE.MS_LAYER_RASTER;
}
else
return false;
if (OnLayerChanged != null)
OnLayerChanged(null, new LayerChangedEventArgs(pLayer.index,strName, true));
RecordLayersState();
//Console.WriteLine(_Map.numlayers.ToString());
return true;
}
public static void cp(IPolygon pPolygon)
{
SaveFileDialog saveFileDialog = new SaveFileDialog();
saveFileDialog.Filter = "Shape文件（*.shp）|*.shp";
saveFileDialog.Title = "新建面shp文件";
saveFileDialog.CheckFileExists = true;
saveFileDialog.RestoreDirectory = true;
DialogResult dialogResult = saveFileDialog.ShowDialog();
string shpPath;
string shpName;
IWorkspaceFactory pWorkSpaceFactory = new ShapefileWorkspaceFactory();
if (dialogResult == DialogResult.OK)
{
IFeatureClass pFeatureClass;
string date = DateTime.Now.ToLongDateString().ToString();
//saveFileDialog.FileName = date + "采集的坐标点";
string fileFullName = saveFileDialog.FileName;
int index = fileFullName.LastIndexOf("\\");
shpName = fileFullName.Substring(index + 1);
shpPath = fileFullName.Substring(0, index);
IFeatureWorkspace pFeatureWorkspace = (IFeatureWorkspace)pWorkSpaceFactory.OpenFromFile(shpPath, 0);
IWorkspaceEdit pWorkspaceEdit = (IWorkspaceEdit)pFeatureWorkspace;
IGeometryDef pGeometryDef = new GeometryDef();
IGeometryDefEdit pGeometryDefEdit = pGeometryDef as IGeometryDefEdit;
pGeometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;//点、线、面什么的
ISpatialReferenceFactory2 ipSpaRefFa = new SpatialReferenceEnvironmentClass();
IGeographicCoordinateSystem ipGeoCorSys = new GeographicCoordinateSystemClass();
ipGeoCorSys = ipSpaRefFa.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
ISpatialReference ipSpaRef = ipGeoCorSys;
pGeometryDefEdit.SpatialReference_2 = ipSpaRef;
/*if (System.IO.File.Exists(saveFileDialog.FileName))
{
if (MessageBox.Show("是否替换原文件？", "询问", MessageBoxButtons.OKCancel, MessageBoxIcon.Question) == DialogResult.OK)
{
pFeatureClass = pFeatureWorkspace.OpenFeatureClass(shpName);
IDataset pDataset = (IDataset)pFeatureClass;
pDataset.Delete();
}
else
return;
}
else
{
fileFullName = null;
return;
}*/
//设置字段
double area;
IFields pFields = new FieldsClass();
IFieldsEdit pFieldsEdit = (IFieldsEdit)pFields;
IField pField = new FieldClass();
IFieldEdit pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "SHAPE";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
pFieldEdit.GeometryDef_2 = pGeometryDef;
pFieldsEdit.AddField(pField);
pField = new FieldClass();
pFieldEdit = (IFieldEdit)pField;
pFieldEdit.Name_2 = "面积";
pFieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble;
pFieldsEdit.AddField(pField);
//创建shp
pWorkspaceEdit.StartEditOperation();//启动编辑
pFeatureClass = pFeatureWorkspace.CreateFeatureClass(shpName, pFields, null, null, esriFeatureType.esriFTSimple, "SHAPE", "");
pPolygon = new PolygonClass();
IArea s = pPolygon as IArea;
area = s.Area;
IFeature pfeature = pFeatureClass.CreateFeature();
pfeature.Shape = pPolygon;
pfeature.set_Value(pfeature.Fields.FindField("面积"), area.ToString());
pfeature.Store();
pWorkspaceEdit.StopEditOperation();
IFeatureLayer pFeaturelayer = new FeatureLayerClass();
pFeaturelayer.FeatureClass = pFeatureClass;
pFeaturelayer.Name = "采集坐标图层";
}
}
public static void Save(IFeatureClass SourceFeatureClass, string WhereClause, string SaveFilePath)
{
IWorkspaceFactory workpsaceFactory = new ShapefileWorkspaceFactory();
IWorkspace workspace = workpsaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(SaveFilePath), 0);
var FieldDict = GetFieldDict(SourceFeatureClass);
IFeatureWorkspace FeatureWorkspace = workspace as IFeatureWorkspace;
IFeatureClass NewFeatureClass = Create(FeatureWorkspace, System.IO.Path.GetFileNameWithoutExtension(SaveFilePath), esriGeometryType.esriGeometryPolyline, FieldDict);
Dictionary<string, int> NewIndexDict = GetFieldIndexDict(NewFeatureClass);
Dictionary<string, int> SourceIndexDict = GetFieldIndexDict(SourceFeatureClass);
IQueryFilter queryFilter = null;
if (!string.IsNullOrEmpty(WhereClause))
{
queryFilter = new QueryFilterClass();
queryFilter.WhereClause = WhereClause;
}
IFeatureBuffer featureBuffer;
IFeatureCursor InsertFeatureCursor=null;
IFeatureCursor SearchfeatureCursor = SourceFeatureClass.Search(queryFilter, false);
IFeature feature = SearchfeatureCursor.NextFeature();
while (feature != null)
{
featureBuffer = NewFeatureClass.CreateFeatureBuffer();
InsertFeatureCursor = NewFeatureClass.Insert(true);
featureBuffer.Shape = feature.Shape;
foreach (var key in NewIndexDict.Keys)
{
if (SourceIndexDict.ContainsKey(key))
{
object val = feature.get_Value(SourceIndexDict[key]);
if (val!=DBNull.Value)
{
featureBuffer.set_Value(NewIndexDict[key], val);
}
}
}
InsertFeatureCursor.InsertFeature(featureBuffer);
InsertFeatureCursor.Flush();
feature = SearchfeatureCursor.NextFeature();
}
System.Runtime.InteropServices.Marshal.ReleaseComObject(SearchfeatureCursor);
if (InsertFeatureCursor != null)
{
System.Runtime.InteropServices.Marshal.ReleaseComObject(InsertFeatureCursor);
}
}
/// <summary>
/// 生成等值线
/// </summary>
/// <params name="targetworkspace"></params>
/// <params name="nameOftargetFeatureClass"></params>
/// <params name="sLayerAliasName"></params>
private bool CreateContours(IWorkspace targetworkspace, string nameOftargetFeatureClass,
string sLayerAliasName, string extent, ref string failInfo)
{
double douElevation;
//设置一个最小值
progressBarControl1.Properties.Minimum = 0;
//设置一个最大值
progressBarControl1.Properties.Maximum = 6;
//设置步长，即每次增加的数
progressBarControl1.Properties.Step = 1;
//设置进度条的样式
progressBarControl1.Properties.ProgressViewStyle = ProgressViewStyle.Solid;
progressBarControl1.Position = 0;
try
{
Geoprocessor GP = new Geoprocessor();
if (withIn(m_strDataFilePath) == false)
return false;
if (TB_Interval.Text == "")
{
MessageBox.Show("请输入正确的等高距！", "提示！");
return false;
}
if (double.TryParse(TB_Interval.Text.Trim(), out douElevation))
{ }
else
{
TB_Interval.Text = null;
MessageBox.Show("请输入正确的等高距！", "提示!");
return false;
}
string featurOut = m_MakeContoursFolder + "\\Cont.shp";
int k = 0;
while (File.Exists(featurOut))
{
k++;
featurOut = m_MakeContoursFolder + "\\Cont" + k.ToString() + ".shp";
}
int countCont = Directory.GetFiles(m_MakeContoursFolder, "Cont*").Length;
if (countCont > 0)
{
featurOut = m_MakeContoursFolder + "\\Cont" + countCont.ToString() + ".shp";
}
if (DrawContours.ConvertASCIIDescretePoint2FeatureClass(GP, m_strDataFilePath, featurOut) == false)
return false;
//执行步长
this.progressBarControl1.PerformStep();
string ContourRaster = m_MakeContoursFolder + "\\Spline";
//string ContourRaster = m_MakeContoursFolder + "\\" + nameOftargetFeatureClass + "_Render";
if (extent != "")//右上左下 xmax ymax xmin ymin
GP.SetEnvironmentValue("Extent", extent);
bool suc = false;
string sRasterLayerAliasName = sLayerAliasName + "渲染图";
ILayer rsLayer = DataEditCommon.GetLayerByName(DataEditCommon.g_pMap, sRasterLayerAliasName);
if (rsLayer != null)
{
DataEditCommon.g_pMap.DeleteLayer(rsLayer);
}
if (Directory.Exists(ContourRaster))
{
try
{
Directory.Delete(ContourRaster, true);
}
catch
{
k = 1;
while (Directory.Exists(ContourRaster))
{
ContourRaster = m_MakeContoursFolder + "\\Spline" + k.ToString();
k++;
}
}
}
switch (CB_InterpolationMethod.Text)
{
case "样条函数插值法":
suc = DrawContours.Interpolate2RasterSpline(GP, featurOut, ContourRaster, CB_SplineType.Text);
break;
case "自然邻域插值法":
suc = DrawContours.Interpolate2RasterNN(GP, featurOut, ContourRaster);
break;
case "克里格插值法":
suc = DrawContours.Interpolate2RasterKriging(GP, featurOut, ContourRaster,
CB_semiVariogramProp.Text, CB_searchRadiusProp.Text);
break;
case "反距离权重插值法":
suc = DrawContours.Interpolate2RasterIDW(GP, featurOut, ContourRaster);
break;
case "趋势面插值法":
suc = DrawContours.TrendToRaster(GP, featurOut, ContourRaster);
break;
}
if (suc == false)
return false;
this.progressBarControl1.PerformStep();
GP = new Geoprocessor();
string R2Contour = m_MakeContoursFolder + "\\Contour.shp";
k = 1;
while (File.Exists(R2Contour))
{
R2Contour = m_MakeContoursFolder + "\\Contour" + k.ToString() + ".shp";
k++;
}
int countContour = Directory.GetFiles(m_MakeContoursFolder, "Contour*").Length;
if (countContour > 0)
{
R2Contour = m_MakeContoursFolder + "\\Contour" + countContour.ToString() + ".shp";
}
if (DrawContours.SplineRasterToContour(GP, ContourRaster, R2Contour, douElevation) == false)
return false;
this.progressBarControl1.PerformStep();
string EvEContour = m_MakeContoursFolder + "\\EvEContour.shp";
k = 1;
while (File.Exists(EvEContour))
{
EvEContour = m_MakeContoursFolder + "\\EvEContour" + k.ToString() + ".shp";
k++;
}
int countEvEContour = Directory.GetFiles(m_MakeContoursFolder, "EvEContour*").Length;
if (countEvEContour > 0)
{
EvEContour = m_MakeContoursFolder + "\\EvEContour" + countEvEContour.ToString() + ".shp";
}
if (DrawContours.FeaturesTo3D(GP, R2Contour, EvEContour) == false)
return false;
this.progressBarControl1.PerformStep();
//获得等值线Shp文件所在的工作空间
IWorkspaceFactory workspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace sourceworkspace = workspaceFactory.OpenFromFile(m_MakeContoursFolder, 0);
//裁切
IPolygon bianjie = GetPolygon();
//导入等值线Shp文件到数据库中
DrawSpecialCommon drawspecial = new DrawSpecialCommon();
string nameOfsourceFeatureClass = EvEContour.Substring(EvEContour.LastIndexOf("\\") + 1);
nameOfsourceFeatureClass = nameOfsourceFeatureClass.Substring(0, nameOfsourceFeatureClass.LastIndexOf("."));
bool Import = false;
List<ziduan> list = new List<ziduan>
{
new ziduan("BID", ""),
new ziduan("mingcheng", sLayerAliasName),
new ziduan("mcid", "0"),
new ziduan("date", DateTime.Now.ToString()),
new ziduan("type", CB_InterpolationMethod.Text)
};
list.Add(new ziduan("BID", ""));
list.Add(new ziduan("mingcheng", sLayerAliasName));
list.Add(new ziduan("mcid", "0"));
list.Add(new ziduan("date", DateTime.Now.ToString()));
list.Add(new ziduan("type", CB_InterpolationMethod.Text));
IFeatureLayer pFeatureLayer = DataEditCommon.GetLayerByName(DataEditCommon.g_pMap, EditLayerName) as IFeatureLayer;
string WhereClause = "mcid='" + "0" + "'";
DataEditCommon.DeleteFeatureByWhereClause(pFeatureLayer, WhereClause);
if (radioBtnKJ.Checked && bianjie != null)
Import = IntersectAll(DataEditCommon.GetFeatureClassByName(sourceworkspace, nameOfsourceFeatureClass), bianjie, list);
else
Import = drawspecial.ShapeImportGDB(sourceworkspace, targetworkspace, nameOfsourceFeatureClass, nameOftargetFeatureClass, list);
this.progressBarControl1.PerformStep();
if (Import == false)
{
MessageBox.Show(sLayerAliasName + "导入数据库失败！");
DataEditCommon.g_axTocControl.Update();
DataEditCommon.g_pAxMapControl.Refresh();
failInfo = sLayerAliasName;
return false;
}
//添加相应的渲染图
IRasterLayer newRasterLayer = new RasterLayerClass();
newRasterLayer.CreateFromFilePath(ContourRaster);
newRasterLayer = IntersectRaster(newRasterLayer, bianjie);
//newRasterLayer.CreateFromDataset(newRasterDs);
newRasterLayer.Name = sRasterLayerAliasName;
UsingRasterStretchColorRampRender(newRasterLayer);
//groupLayer.Add(newRasterLayer as ILayer);
int indexPosition = DataEditCommon.g_pAxMapControl.LayerCount;//GroupLayer序号
//判断MapControl中是否存在该图层
ILayer mLayer = DataEditCommon.GetLayerByName(DataEditCommon.g_pMap, sRasterLayerAliasName);
if (mLayer != null)
{
DataEditCommon.g_pMyMapCtrl.Map.DeleteLayer(mLayer);
indexPosition--;
}
DataEditCommon.g_pAxMapControl.AddLayer(newRasterLayer as ILayer, indexPosition);//添加到最下面
DataEditCommon.g_pAxMapControl.ActiveView.Extent = newRasterLayer.AreaOfInterest;
DataEditCommon.g_axTocControl.Update();
DataEditCommon.g_pAxMapControl.Refresh();
this.progressBarControl1.PerformStep();
return true;
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
return false;
}
}
private void ExportTable2DBF(ITable ptable, string strPath, string strFilename)
{
IDataset pDataset = ptable as IDataset;
IDatasetName pDatasetName = pDataset.FullName as IDatasetName;
IWorkspaceFactory pWF = new ShapefileWorkspaceFactory();
IWorkspace pWorkspace = pWF.OpenFromFile(strPath, 0);
IDataset pDatasetOUT = pWorkspace as IDataset;
IWorkspaceName pWorkspaceName = pDatasetOUT.FullName as IWorkspaceName;
IDatasetName pOUTDatasetName = new TableNameClass ();
pOUTDatasetName.Name = strFilename;
pOUTDatasetName.WorkspaceName = pWorkspaceName;
IExportOperation pExpOp = new ExportOperationClass();
pExpOp.ExportTable(pDatasetName, null, null, pOUTDatasetName, 0);
}
