gps数据导入到sde,shp导入SDE

http://www.cnblogs.com/bpcbcm/archive/2009/05/16/bpc.html

以配电网gis为例，将采集好的电网数据gps数据导入到sde数据库， 以杆塔线路为例，杆塔数据绘制好之后自动联线，同时保持网络关系。

 /// <summary>
  /// 获取文件所有要添加的杆塔的所有的线路名称
  /// /// </summary>
  /// 
  public void Getlinenames()
  {
  //gpsEdit.CurrentLayer=
  IFeatureWorkspace pFeatWS;
  IWorkspaceFactory pWorkspaceFactory;
  pWorkspaceFactory = new ShapefileWorkspaceFactoryClass();
   
  if (Fileworkspacestring == null) return;
  pFeatWS = pWorkspaceFactory.OpenFromFile(Fileworkspacestring, 0) as IFeatureWorkspace;  

  foreach (String file in listnames)
  {
  // Create a PictureBox.
  if (file == "Export_Output.shp")
  {
  try
  {
  IQueryFilter queryFilter = new QueryFilterClass();
  queryFilter.SubFields = "xlmc,zlmc";
  IFeatureClass pFeatureClass = pFeatWS.OpenFeatureClass(file);
  ISelectionSet selectionSet = pFeatureClass.Select(null, esriSelectionType.esriSelectionTypeIDSet, esriSelectionOption.esriSelectionOptionNormal, null);

  IEnumIDs enumIDs = selectionSet.IDs;
  int fieldIndex = pFeatureClass.FindField("xlmc");
  int fieldIndex1 = pFeatureClass.FindField("zlmc");

  IFeature feature;
  int iD = enumIDs.Next();
  xlname xl = new xlname();
  int isfirst = 0;
  while (iD != -1) //-1 is reutned after the last valid ID has been reached  
  {
   
  if (isfirst == 0)
  {
  feature = pFeatureClass.GetFeature(iD);
  xl.xlmc = feature.get_Value(fieldIndex).ToString();
  xl.zlmc = feature.get_Value(fieldIndex1).ToString();
  xlnames.Add(xl);
  isfirst = 1;
  }
  else
  {
  feature = pFeatureClass.GetFeature(iD);
  xl.xlmc = feature.get_Value(fieldIndex).ToString();
  xl.zlmc = feature.get_Value(fieldIndex1).ToString();
  bool iscz = true;
  foreach (xlname xlmc in xlnames)
  {

  if ((xlmc.xlmc == xl.xlmc) && (xlmc.zlmc == xl.zlmc))
  {
  iscz = false;//是否存在

  }
  }
  if (iscz == true) //不存在就添加
  {
  xlnames.Add(xl);
  }
  }

  iD = enumIDs.Next();
  }
  int dd = xlnames.Count;

  }

  catch (Exception e)
  {
  Console.WriteLine(e.Message.ToString());
  }

  }
  }  
  }
  /// <summary>
  /// 绘制点状设备
  /// </summary>
  /// <param name="iflayer"></param>
  /// <param name="pPoint"></param>
  private void Drawpoint(ILayer iflayer,IPoint pPoint)
  {
  try
  {
  if (iflayer == null) return;
  if (!(iflayer is IGeoFeatureLayer)) return;
  IFeatureLayer pFeatureLayer = (IFeatureLayer)iflayer;
  IFeatureClass mfeatureClass = pFeatureLayer.FeatureClass;
  IDataset pDataset = (IDataset)pFeatureLayer.FeatureClass;
  if (pDataset == null) return;
   
  // 开始编辑,并设置Undo/Redo 为可用
  IWorkspaceEdit pWorkspaceEdit = (IWorkspaceEdit)pDataset.Workspace;
  if (!pWorkspaceEdit.IsBeingEdited())
  {
  pWorkspaceEdit.StartEditing(true);
  pWorkspaceEdit.EnableUndoRedo();
  IFeature myfeature;

  IFeatureClassWrite fr = mfeatureClass as IFeatureClassWrite;//绘制杆塔
  pWorkspaceEdit.StartEditOperation();
  myfeature = mfeatureClass.CreateFeature();  

  myfeature.Shape = pPoint;  

  myfeature.Store();
  fr.WriteFeature(myfeature);

  pWorkspaceEdit.StopEditOperation();
  pWorkspaceEdit.StopEditing (true );
  m_pMap.ClearSelection();
  IActiveView pActiveView = (IActiveView)m_pMap;
  pActiveView.Refresh();


  }
  }
  catch (Exception e)
  {
  Console.WriteLine(e.Message.ToString());
  }

 
  }
  /// <summary>
  /// 画线路
  /// </summary>
  public void Drawline(ILayer iflayer, IPointCollection pointCollection)
  {
  try
  {
  if (iflayer == null) return;
  if (!(iflayer is IGeoFeatureLayer)) return;
  IFeatureLayer pFeatureLayer = (IFeatureLayer)iflayer;
  IFeatureClass mfeatureClass = pFeatureLayer.FeatureClass;
  IDataset pDataset = (IDataset)pFeatureLayer.FeatureClass;
  if (pDataset == null) return;

  // 开始编辑,并设置Undo/Redo 为可用
  IWorkspaceEdit pWorkspaceEdit = (IWorkspaceEdit)pDataset.Workspace;
  if (!pWorkspaceEdit.IsBeingEdited())
  {
  pWorkspaceEdit.StartEditing(true);
  pWorkspaceEdit.EnableUndoRedo();
  IFeature myfeature;

  IFeatureClassWrite fr = mfeatureClass as IFeatureClassWrite;//绘制杆塔
  pWorkspaceEdit.StartEditOperation();
  myfeature = mfeatureClass.CreateFeature();  
  pWorkspaceEdit.StartEditOperation();
  myfeature = mfeatureClass.CreateFeature();
  //IGeometry gGeometry = (IGeometry)pointCollection ;
  IPolyline pPolyline = pointCollection as IPolyline;

  myfeature.Shape = pPolyline;  

  myfeature.Store();
  fr.WriteFeature(myfeature);
  IComplexEdgeFeature Ccomplexedte = myfeature as IComplexEdgeFeature;//建立连通性
  Ccomplexedte.ConnectAtIntermediateVertices();

  pWorkspaceEdit.StopEditOperation();
  pWorkspaceEdit.StopEditing(true);



   
   

  m_pMap.ClearSelection();
  IActiveView pActiveView = (IActiveView)m_pMap;
  pActiveView.Refresh();
  }
  }
  catch (Exception e)
  {
  Console.WriteLine(e.Message.ToString());
  }
  }  
   

/// <summary>
  /// 绘制gps数据，杆塔，线路
  /// </summary>

public void Draw() 

{

  IFeatureWorkspace pFeatWS;
  IWorkspaceFactory pWorkspaceFactory;
   
  pWorkspaceFactory = new ShapefileWorkspaceFactoryClass();  
  if (Fileworkspacestring == null) return;
  pFeatWS = pWorkspaceFactory.OpenFromFile(Fileworkspacestring, 0) as IFeatureWorkspace;  
  foreach (String file in listnames)
  {  
  if (file == "Export_Output.shp")//布绘杆塔  
  {
  foreach (xlname xlmc in xlnames)
  {
  IFeature featurepoint;
   
  int idline = GetLayerid(m_pMap, "GIS.lines");//获取线图层
  ILayer ilayerline = m_pMap.get_Layer(idline);

  int idpoint = GetLayerid(m_pMap, "gt");//获取点图层
  ILayer iflayer = m_pMap.get_Layer(idpoint);

  IPolyline mypolylin = new PolylineClass();
  IPointCollection pointCollection = mypolylin as IPointCollection;
  object missingVal = System.Reflection.Missing.Value;

  IQueryFilter queryFilter = new QueryFilterClass();
  queryFilter.WhereClause = "xlmc='"+xlmc .xlmc+"'"+" and zlmc='" + xlmc.zlmc+"'";
  IFeatureClass pFeatureClass = pFeatWS.OpenFeatureClass(file);
  //IQueryFilterDefinition queryFilterDefinition = (IQueryFilterDefinition)queryFilter; 
  //queryFilterDefinition.PostfixClause = "ORDER BY gh";
  IFields fField = pFeatureClass.Fields;
  for (int fieldindex = 0; fieldindex < fField.FieldCount;fieldindex++ )
  {
  string dd = fField.get_Field(fieldindex).Name;
  string ddd = pFeatureClass.ShapeFieldName;
  }
   
  ISelectionSet selectionSet = pFeatureClass.Select(queryFilter, esriSelectionType.esriSelectionTypeIDSet, esriSelectionOption.esriSelectionOptionNormal, null);
  int fieldIndexgh = pFeatureClass.FindField("gh");
  int fieldIndexzjgh = pFeatureClass.FindField("zjgh");
  int fieldindexxlmc = pFeatureClass.FindField("xlmc");
  int fieldindexzlmc = pFeatureClass.FindField("zlmc");
  int fieldindextjxl = pFeatureClass.FindField("tjxl");
  int fieldindextjgh = pFeatureClass.FindField("tjgh");
  IEnumIDs enumIDs = selectionSet.IDs;  

  for (int j = 1; j <= selectionSet.Count;j++ )
  {
   
  bool istjgh = false;
  enumIDs.Reset();
  int iD = enumIDs.Next();
  while (iD != -1) //-1 is reutned after the last valid ID has been reached  
  {
  featurepoint = pFeatureClass.GetFeature(iD);
  IGeometry geometrygt = featurepoint.Shape;//获取shp文件中的对象
  IPoint pointgt = geometrygt as IPoint;
  gtzb gtxx = new gtzb();  
  gtxx.gh = featurepoint.get_Value(fieldIndexgh).ToString();
  gtxx.xlmc = featurepoint.get_Value(fieldindexxlmc).ToString();
  gtxx.zlmc = featurepoint.get_Value(fieldindexzlmc).ToString();
  gtxx.zjgh = featurepoint.get_Value(fieldIndexzjgh).ToString();
  if (gtxx.gh == j.ToString() && gtxx.zjgh.ToString().Trim() =="")
  {  
  Drawpoint(iflayer, pointgt); //将shp文件中的对象绘制到数据库  
  pointCollection.AddPoint(pointgt, ref missingVal, ref missingVal);
  }
  else if (gtxx.gh == j.ToString() && gtxx.zjgh.ToString().Trim() != "")
  {
   
  istjgh = true;  
   
  }
  iD = enumIDs.Next();
  }
  if (istjgh)
  {
  IQueryFilter queryFiltersub = new QueryFilterClass();
  queryFiltersub.WhereClause = "xlmc='" + xlmc.xlmc + "'" + " and zlmc='" + xlmc.zlmc + "'and zjgh='" + j + "'";
  ISelectionSet selectionSetsub = pFeatureClass.Select(queryFiltersub, esriSelectionType.esriSelectionTypeIDSet, esriSelectionOption.esriSelectionOptionNormal, null);
   
  for (int subj = 1; subj <= selectionSetsub.Count; subj++)
  {
  IEnumIDs enumsubIDs = selectionSetsub.IDs; 
  enumsubIDs.Reset();
  int iDsub = enumsubIDs.Next();
  while (iDsub != -1) //-1 is reutned after the last valid ID has been reached  
  {
  featurepoint = pFeatureClass.GetFeature(iDsub);
  IGeometry geometrygt = featurepoint.Shape;
  IPoint pointgt = geometrygt as IPoint;
  gtzb gtxx = new gtzb();
  gtxx.longitude = pointgt.X;
  gtxx.latitude = pointgt.Y;
  gtxx.elevation = pointgt.Z;
  gtxx.gh = featurepoint.get_Value(fieldIndexgh).ToString();
  gtxx.zjgh = featurepoint.get_Value(fieldIndexzjgh).ToString();
  gtxx.xlmc = featurepoint.get_Value(fieldindexxlmc).ToString();
  gtxx.zlmc = featurepoint.get_Value(fieldindexzlmc).ToString();
  if (gtxx.zjgh == subj.ToString())
  {
  Drawpoint(iflayer, pointgt);  
  pointCollection.AddPoint(pointgt, ref missingVal, ref missingVal);
  }
  iDsub = enumsubIDs.Next();
  }
  } 
  }
  }

  //绘制线路
  Drawline(ilayerline, pointCollection);
  }

  }

}