//OK Loop
public ZhFeatureClass GetZhFeatureClass(string featureClassName)
{
ZhFeatureClass ZHfeatcls = null;
IFeatureClass fc = GetFeatureClass(featureClassName);
if (fc != null)
{
switch (fc.ShapeType)
{
case esriGeometryType.esriGeometryPolygon:
ZHfeatcls = new ZhPolygonFeatureClass(fc);
break;
case esriGeometryType.esriGeometryPolyline:
case esriGeometryType.esriGeometryLine:
ZHfeatcls = new ZhPolylineFeatureClass(fc);
break;
case esriGeometryType.esriGeometryPoint:
case esriGeometryType.esriGeometryMultipoint:
ZHfeatcls = new ZhPointFeatureClass(fc);
break;
}
}
return(ZHfeatcls);
}
/// <summary>
/// getZHFeatureClass
/// </summary>
/// <returns></returns>
public ZhFeatureClass getZHFeatureClass()
{
IFeatureClass fc = this.getIFeatureClass();
ZhFeatureClass zhfc = new ZhPointFeatureClass(fc);
return(zhfc);
}
/// <summary>
/// deleteAllFeatures
/// </summary>
/// <returns></returns>
public bool deleteAllFeatures()
{
bool rbc = false;
try
{
IFeatureClass fc = this.getIFeatureClass();
if (fc != null)
{
ZhPointFeatureClass zhfc = new ZhPointFeatureClass(fc);
zhfc.DeleteAllFeatures();
TokayWorkspace.ComRelease(fc);
fc = null;
zhfc = null;
rbc = true;
}
}
catch (Exception ee)
{
System.Diagnostics.Debug.WriteLine(ee.ToString());
rbc = false;
}
return(rbc);
}
//
//输出外环面 shape
public bool OutPutPolygonExteriorRing(IFeatureClass featclass, string outShapeFilePath, bool IsGetExtRingOfHaveInnerRing)
{
bool rbc = false;
if (System.IO.File.Exists(outShapeFilePath) == true)
{ //输出Shape文件
LocalShapeFileOperator delshpOp = new LocalShapeFileOperator();
delshpOp.LocalShapePathFileName = outShapeFilePath;
delshpOp.DeleteShapeFile();
delshpOp.Dispose();
}
esriGeometryType geoType = featclass.ShapeType;
//创建新的Shape文件
LocalShapeFileOperator shpOp = new LocalShapeFileOperator();
shpOp.LocalShapePathFileName = outShapeFilePath;
string dir = shpOp.getDir(outShapeFilePath);
string name = shpOp.getFileName(outShapeFilePath);
ShapefileWorkspaceFactoryClass objwsf = new ShapefileWorkspaceFactoryClass();
IWorkspace objws = objwsf.OpenFromFile(dir, 0);
IFeatureWorkspace feaureWS = objws as IFeatureWorkspace;
//设置投影
ISpatialReference sr = (featclass as IGeoDataset).SpatialReference;
shpOp.ShapeSpatialReference = sr;
//设置shp文件的oid字段和几何字段
if (IsObjGeo == false)
{
if (featclass.HasOID == true)
{
this.OIDFieldName = featclass.OIDFieldName;
}
this.GeometryFieldName = featclass.ShapeFieldName;
}
shpOp.OIDFieldName = this.OIDFieldName;
shpOp.GeometryFieldName = this.GeometryFieldName;
//创建要素类
IFeatureClass fc = shpOp.CreateFeatureClass(feaureWS, name, esriFeatureType.esriFTSimple, geoType);
if (fc != null)
{
TokayWorkspace.ComRelease(fc);
fc = null;
rbc = true;
}
TokayWorkspace.ComRelease(objws);
feaureWS = null;
objws = null;
objwsf = null;
//拷贝字段结构
IFeatureClass Dfeatclass = shpOp.getIFeatureClass();
ZhFeatureClass Szhfeatclass = new ZhPointFeatureClass(featclass);
Szhfeatclass.CopyFieldsToObjectFeatureClass(Dfeatclass);
//拷贝几何对象和属性数据
ZhFeature Szhfeat = null;
IFeatureBuffer DfeatBuffer = null;
IFeatureCursor featCurInsert = Dfeatclass.Insert(true);
//获取总要素个数
int fIndex = 0;
int tmpFCount = Szhfeatclass.FeatureClass.FeatureCount(null);
frmProgressBar1 pb = new frmProgressBar1();
pb.Text = "正在输出shp文件...";
pb.Caption1.Text = "正在输出shp文件...";
pb.progressBar1.Maximum = tmpFCount + 1;
pb.progressBar1.Value = 0;
pb.Show(this.ParentForm);
Application.DoEvents();
object refobj = Type.Missing;
IPolygon s_p = null;
IRing[] ExtRingArray = null;
IRing[] InnerRingArray = null;
IFeatureCursor featcur = Szhfeatclass.FeatureClass.Search(null, false);
IFeature feat = featcur.NextFeature();
while (feat != null)
{
fIndex += 1;
if (fIndex % 200 == 0)
{
pb.Caption1.Text = "已输出要素个数:" + fIndex.ToString() + "/总" + tmpFCount.ToString();
pb.progressBar1.Value = fIndex;
Application.DoEvents();
featCurInsert.Flush();
}
s_p = feat.ShapeCopy as IPolygon;
//获取外环面几何对象
ExtRingArray = PolygonHelper.GetExteriorRings(s_p);
if (ExtRingArray != null && ExtRingArray.Length > 0)
{
foreach (IRing r in ExtRingArray)
{
if (IsGetExtRingOfHaveInnerRing == true)
{
InnerRingArray = PolygonHelper.GetInteriorRingsByExterior(r, s_p);
if (InnerRingArray != null && InnerRingArray.Length > 0)
{ //有内环
//外环构面并输出保存
DfeatBuffer = Dfeatclass.CreateFeatureBuffer();
PolygonClass pclass = new PolygonClass();
pclass.AddGeometry(r, ref refobj, ref refobj);
pclass.SimplifyPreserveFromTo();
pclass.SpatialReference = sr;
DfeatBuffer.Shape = pclass;
//拷贝属性数据
Szhfeat = new ZHFeaturePoint(feat);
Szhfeat.CopyField(ref DfeatBuffer);
featCurInsert.InsertFeature(DfeatBuffer);
//
}
}
else
{
//外环构面并输出保存
DfeatBuffer = Dfeatclass.CreateFeatureBuffer();
PolygonClass pclass = new PolygonClass();
pclass.AddGeometry(r, ref refobj, ref refobj);
pclass.SimplifyPreserveFromTo();
pclass.SpatialReference = sr;
DfeatBuffer.Shape = pclass;
//拷贝属性数据
Szhfeat = new ZHFeaturePoint(feat);
Szhfeat.CopyField(ref DfeatBuffer);
featCurInsert.InsertFeature(DfeatBuffer);
//
}
}
}
feat = featcur.NextFeature();
}
featCurInsert.Flush();
TokayWorkspace.ComRelease(featcur);
featcur = null;
TokayWorkspace.ComRelease(featCurInsert);
featCurInsert = null;
rbc = true;
pb.Close();
pb.Dispose();
pb = null;
return(rbc);
}
//输出单个图层中所有要素集合到shp文件中功能 OK
public bool OutPut(IFeatureClass featclass, string outShapeFilePath)
{
bool rbc = false;
if (System.IO.File.Exists(outShapeFilePath) == true)
{ //输出Shape文件
LocalShapeFileOperator delshpOp = new LocalShapeFileOperator();
delshpOp.LocalShapePathFileName = outShapeFilePath;
delshpOp.DeleteShapeFile();
delshpOp.Dispose();
}
esriGeometryType geoType = featclass.ShapeType;
//创建新的Shape文件
LocalShapeFileOperator shpOp = new LocalShapeFileOperator();
shpOp.LocalShapePathFileName = outShapeFilePath;
string dir = shpOp.getDir(outShapeFilePath);
string name = shpOp.getFileName(outShapeFilePath);
ShapefileWorkspaceFactoryClass objwsf = new ShapefileWorkspaceFactoryClass();
IWorkspace objws = objwsf.OpenFromFile(dir, 0);
IFeatureWorkspace feaureWS = objws as IFeatureWorkspace;
//设置投影
shpOp.ShapeSpatialReference = (featclass as IGeoDataset).SpatialReference;
//设置shp文件的oid字段和几何字段
if (IsObjGeo == false)
{
if (featclass.HasOID == true)
{
this.OIDFieldName = featclass.OIDFieldName;
}
this.GeometryFieldName = featclass.ShapeFieldName;
}
shpOp.OIDFieldName = this.OIDFieldName;
shpOp.GeometryFieldName = this.GeometryFieldName;
//创建要素类
IFeatureClass fc = shpOp.CreateFeatureClass(feaureWS, name, esriFeatureType.esriFTSimple, geoType);
if (fc != null)
{
TokayWorkspace.ComRelease(fc);
fc = null;
rbc = true;
}
TokayWorkspace.ComRelease(objws);
feaureWS = null;
objws = null;
objwsf = null;
//拷贝字段结构
IFeatureClass Dfeatclass = shpOp.getIFeatureClass();
ZhFeatureClass Szhfeatclass = new ZhPointFeatureClass(featclass);
Szhfeatclass.CopyFieldsToObjectFeatureClass(Dfeatclass);
//拷贝几何对象和属性数据
ZhFeature Szhfeat = null;
IFeatureBuffer DfeatBuffer = null;
ZhFeature Dzhfeat = null;
IFeatureCursor featCurInsert = Dfeatclass.Insert(true);
//获取总要素个数
int fIndex = 0;
int tmpFCount = Szhfeatclass.FeatureClass.FeatureCount(null);
frmProgressBar1 pb = new frmProgressBar1();
pb.Text = "正在输出shp文件...";
pb.Caption1.Text = "正在输出shp文件...";
pb.progressBar1.Maximum = tmpFCount + 1;
pb.progressBar1.Value = 0;
pb.Show(this.ParentForm);
Application.DoEvents();
//
IGeometry geo = null;
//
IFeatureCursor featcur = Szhfeatclass.FeatureClass.Search(null, false);
IFeature feat = featcur.NextFeature();
while (feat != null)
{
fIndex += 1;
if (fIndex % 500 == 0)
{
pb.Caption1.Text = "已输出要素个数:" + fIndex.ToString() + "/总" + tmpFCount.ToString();
pb.progressBar1.Value = fIndex;
Application.DoEvents();
featCurInsert.Flush();
}
DfeatBuffer = Dfeatclass.CreateFeatureBuffer();
//拷贝几何对象
geo = feat.ShapeCopy;
//去掉Z值
IZAware pZaware = geo as IZAware;
if (pZaware.ZAware == true)
{
pZaware.DropZs();
pZaware.ZAware = false;
}
//
DfeatBuffer.Shape = geo;
//
//拷贝属性数据
Szhfeat = new ZHFeaturePoint(feat);
Szhfeat.CopyField(ref DfeatBuffer);
featCurInsert.InsertFeature(DfeatBuffer);
feat = featcur.NextFeature();
}
featCurInsert.Flush();
TokayWorkspace.ComRelease(featcur);
featcur = null;
TokayWorkspace.ComRelease(featCurInsert);
featCurInsert = null;
rbc = true;
pb.Close();
pb.Dispose();
pb = null;
return(rbc);
}
//输出通过一个图层已选中的要素去选择输出其他图层shp文件功能 OK
public bool OutPutSelectedFeatureByOverlap(IFeatureLayer SelectedFeatLayer, IFeatureClass outFeatClass, string outShapeFilePath)
{
bool rbc = false;
if (System.IO.File.Exists(outShapeFilePath) == true)
{ //输出Shape文件
LocalShapeFileOperator delshpOp = new LocalShapeFileOperator();
delshpOp.LocalShapePathFileName = outShapeFilePath;
delshpOp.DeleteShapeFile();
delshpOp.Dispose();
}
esriGeometryType geoType = outFeatClass.ShapeType;
//创建新的Shape文件
LocalShapeFileOperator shpOp = new LocalShapeFileOperator();
shpOp.LocalShapePathFileName = outShapeFilePath;
shpOp.GeometryType = geoType;
//
string dir = shpOp.getDir(outShapeFilePath);
string name = shpOp.getFileName(outShapeFilePath);
ShapefileWorkspaceFactoryClass objwsf = new ShapefileWorkspaceFactoryClass();
IWorkspace objws = objwsf.OpenFromFile(dir, 0);
IFeatureWorkspace feaureWS = objws as IFeatureWorkspace;
//设置投影
shpOp.ShapeSpatialReference = (outFeatClass as IGeoDataset).SpatialReference;
//设置shp文件的oid字段和几何字段
if (IsObjGeo == false)
{
if (outFeatClass.HasOID == true)
{
this.OIDFieldName = outFeatClass.OIDFieldName;
}
this.GeometryFieldName = outFeatClass.ShapeFieldName;
}
shpOp.OIDFieldName = this.OIDFieldName;
shpOp.GeometryFieldName = this.GeometryFieldName;
//创建要素类
IFeatureClass fc = shpOp.CreateFeatureClass(feaureWS, name, esriFeatureType.esriFTSimple, geoType);
if (fc != null)
{
TokayWorkspace.ComRelease(fc);
fc = null;
rbc = true;
}
TokayWorkspace.ComRelease(objws);
feaureWS = null;
objws = null;
objwsf = null;
//拷贝字段结构
IFeatureClass Dfeatclass = shpOp.getIFeatureClass();
ZhFeatureClass Szhfeatclass = new ZhPointFeatureClass(outFeatClass);
Szhfeatclass.CopyFieldsToObjectFeatureClass(Dfeatclass);
//拷贝几何对象和属性数据
ZhFeature Szhfeat = null;
IFeatureBuffer DfeatBuffer = null;
ZhFeature Dzhfeat = null;
IFeatureCursor featCurInsert = Dfeatclass.Insert(true);
//获取总要素个数
int fIndex = 0;
int tmpFCount = Szhfeatclass.FeatureClass.FeatureCount(null);
frmProgressBar1 pb = new frmProgressBar1();
pb.Text = "正在输出shp文件...";
pb.Caption1.Text = "正在输出shp文件...";
pb.progressBar1.Maximum = tmpFCount + 1;
pb.progressBar1.Value = 0;
pb.Show(this.ParentForm);
Application.DoEvents();
//获取已选择要的要素集合
IFeatureSelection featSel = SelectedFeatLayer as IFeatureSelection;
ICursor SelectedCur = null;
featSel.SelectionSet.Search(null, false, out SelectedCur);
if (SelectedCur != null)
{
IFeatureCursor Selectedfeatcur = SelectedCur as IFeatureCursor;// Szhfeatclass.FeatureClass.Search(null, false);
if (Selectedfeatcur != null)
{
IFeature Selectedfeat = Selectedfeatcur.NextFeature();
while (Selectedfeat != null)
{
ISpatialFilter sFilter = new SpatialFilterClass();
sFilter.Geometry = Selectedfeat.ShapeCopy;
sFilter.GeometryField = Szhfeatclass.FeatureClass.ShapeFieldName;
sFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects; //相交
//总个数
pb.progressBar1.Maximum = Szhfeatclass.FeatureClass.FeatureCount(sFilter);
pb.progressBar1.Value = 0;
Application.DoEvents();
fIndex = 0;
//
IFeatureCursor featcur = Szhfeatclass.FeatureClass.Search(sFilter, false);
IFeature feat = featcur.NextFeature();
while (feat != null)
{
fIndex += 1;
//
pb.Caption1.Text = "已输出要素个数... 第" + fIndex.ToString() + "个/总" + pb.progressBar1.Maximum.ToString() + "个";
pb.progressBar1.Value = fIndex;
Application.DoEvents();
if (fIndex % 500 == 0)
{
featCurInsert.Flush();
}
//
DfeatBuffer = Dfeatclass.CreateFeatureBuffer();
//拷贝几何对象
DfeatBuffer.Shape = feat.ShapeCopy;
//拷贝属性数据
Szhfeat = new ZHFeaturePoint(feat);
Szhfeat.CopyField(ref DfeatBuffer);
//
featCurInsert.InsertFeature(DfeatBuffer);
//下一个要素
feat = featcur.NextFeature();
}
featCurInsert.Flush();
if (featcur != null)
{
TokayWorkspace.ComRelease(featcur);
featcur = null;
}
//下一个选中的要素
Selectedfeat = Selectedfeatcur.NextFeature();
}
if (Selectedfeatcur != null)
{
TokayWorkspace.ComRelease(Selectedfeatcur);
Selectedfeatcur = null;
}
}
if (featCurInsert != null)
{
TokayWorkspace.ComRelease(featCurInsert);
featCurInsert = null;
}
rbc = true;
}
pb.Close();
pb.Dispose();
pb = null;
//
return(rbc);
}
//输出单个图层中已选中的要素集合到shp文件中功能 OK
public bool OutPutSelectedFeature(IFeatureLayer featLayer, string outShapeFilePath)
{
IFeatureClass featclass = featLayer.FeatureClass;
bool rbc = false;
if (System.IO.File.Exists(outShapeFilePath) == true)
{ //输出Shape文件
LocalShapeFileOperator delshpOp = new LocalShapeFileOperator();
delshpOp.LocalShapePathFileName = outShapeFilePath;
delshpOp.DeleteShapeFile();
delshpOp.Dispose();
}
esriGeometryType geoType = featclass.ShapeType;
//创建新的Shape文件
LocalShapeFileOperator shpOp = new LocalShapeFileOperator();
shpOp.LocalShapePathFileName = outShapeFilePath;
string dir = shpOp.getDir(outShapeFilePath);
string name = shpOp.getFileName(outShapeFilePath);
ShapefileWorkspaceFactoryClass objwsf = new ShapefileWorkspaceFactoryClass();
IWorkspace objws = objwsf.OpenFromFile(dir, 0);
IFeatureWorkspace feaureWS = objws as IFeatureWorkspace;
//设置投影
shpOp.ShapeSpatialReference = (featclass as IGeoDataset).SpatialReference;
//设置shp文件的oid字段和几何字段
if (IsObjGeo == false)
{
if (featclass.HasOID == true)
{
this.OIDFieldName = featclass.OIDFieldName;
}
this.GeometryFieldName = featclass.ShapeFieldName;
}
shpOp.OIDFieldName = this.OIDFieldName;
shpOp.GeometryFieldName = this.GeometryFieldName;
//创建要素类
IFeatureClass fc = shpOp.CreateFeatureClass(feaureWS, name, esriFeatureType.esriFTSimple, geoType);
if (fc != null)
{
TokayWorkspace.ComRelease(fc);
fc = null;
rbc = true;
}
TokayWorkspace.ComRelease(objws);
feaureWS = null;
objws = null;
objwsf = null;
//拷贝字段结构
IFeatureClass Dfeatclass = shpOp.getIFeatureClass();
ZhFeatureClass Szhfeatclass = new ZhPointFeatureClass(featclass);
Szhfeatclass.CopyFieldsToObjectFeatureClass(Dfeatclass);
//拷贝几何对象和属性数据
ZhFeature Szhfeat = null;
IFeatureBuffer DfeatBuffer = null;
ZhFeature Dzhfeat = null;
IFeatureCursor featCurInsert = Dfeatclass.Insert(true);
//获取已选择要的要素集合
IFeatureSelection featSel = featLayer as IFeatureSelection;
ICursor cur = null;
featSel.SelectionSet.Search(null, false, out cur);
if (cur != null)
{
IFeatureCursor featcur = cur as IFeatureCursor;// Szhfeatclass.FeatureClass.Search(null, false);
if (featcur != null)
{
IFeature feat = featcur.NextFeature();
while (feat != null)
{
DfeatBuffer = Dfeatclass.CreateFeatureBuffer();
//拷贝几何对象
DfeatBuffer.Shape = feat.ShapeCopy;
//拷贝属性数据
Szhfeat = new ZHFeaturePoint(feat);
Szhfeat.CopyField(ref DfeatBuffer);
featCurInsert.InsertFeature(DfeatBuffer);
feat = featcur.NextFeature();
}
}
TokayWorkspace.ComRelease(featcur);
featcur = null;
}
featCurInsert.Flush();
TokayWorkspace.ComRelease(featCurInsert);
featCurInsert = null;
rbc = true;
return(rbc);
}
private void ShapeFileProjectionConvert(ISpatialReference SourceSpatialReference, ISpatialReference ObjectSpatialReference, IFeatureClass SourceFeatureClass, IFeatureClass ObjectFeatureClass)
{
int SourceFactoryCode = SourceSpatialReference.FactoryCode;
int ObjectFactoryCode = ObjectSpatialReference.FactoryCode;
//if (SourceFactoryCode == 0 || ObjectFactoryCode == 0)
//{
// MessageBox.Show("源文件是自定义的投影坐标系统,这里暂不进行转换", "投影转换提示");
// return;
//}
EnumProjectionDatum SourceProjection = EnumProjectionDatum.Bejing54; //默认
EnumStrip SourceStrip = EnumStrip.Strip3; //默认
EnumProjectionDatum ObjectProjection = EnumProjectionDatum.Xian80; //默认
EnumStrip ObjectStrip = EnumStrip.Strip3; //默认
//源中央子午线
double L0 = (SourceSpatialReference as IProjectedCoordinateSystem4GEN).GetCentralLongitude();
bool IsBigNumber_Source = true; //源文件默认是大数
bool IsBigNumber_Object = true; //目标默认是大数
//------------------------------------------------------------------------------------------
//高斯投影大数情况分析
int BJ54D3_25srid = (int)esriSRProjCS4Type.esriSRProjCS_Beijing1954_3_Degree_GK_Zone_25;
int BJ54D3_45srid = (int)esriSRProjCS4Type.esriSRProjCS_Beijing1954_3_Degree_GK_Zone_45;
int BJ54D6_13srid = (int)esriSRProjCSType.esriSRProjCS_Beijing1954GK_13;
int BJ54D6_23srid = (int)esriSRProjCSType.esriSRProjCS_Beijing1954GK_23;
int Xian80D3_25srid = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_3_Degree_GK_Zone_25;
int Xian80D3_45srid = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_3_Degree_GK_Zone_45;
int Xian80D6_13srid = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_GK_Zone_13;
int Xian80D6_23srid = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_GK_Zone_23;
//源文件分析
//BeiJing1954 3 degree
if (SourceFactoryCode >= BJ54D3_25srid && SourceFactoryCode <= BJ54D3_45srid)
{
SourceProjection = EnumProjectionDatum.Bejing54;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = true;
}
//BeiJing1954 6 degree
if (SourceFactoryCode >= BJ54D6_13srid && SourceFactoryCode <= BJ54D6_23srid)
{
SourceProjection = EnumProjectionDatum.Bejing54;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = true;
}
//Xian1980 3 degree
if (SourceFactoryCode >= Xian80D3_25srid && SourceFactoryCode <= Xian80D3_45srid)
{
SourceProjection = EnumProjectionDatum.Xian80;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = true;
}
//Xian1980 6 degree
if (SourceFactoryCode >= Xian80D6_13srid && SourceFactoryCode <= Xian80D6_23srid)
{
SourceProjection = EnumProjectionDatum.Xian80;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = true;
}
//目标文件分析
//BeiJing1954 3 degree
if (ObjectFactoryCode >= BJ54D3_25srid && ObjectFactoryCode <= BJ54D3_45srid)
{
ObjectProjection = EnumProjectionDatum.Bejing54;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = true;
}
//BeiJing1954 6 degree
if (ObjectFactoryCode >= BJ54D6_13srid && ObjectFactoryCode <= BJ54D6_23srid)
{
ObjectProjection = EnumProjectionDatum.Bejing54;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = true;
}
//Xian1980 3 degree
if (ObjectFactoryCode >= Xian80D3_25srid && ObjectFactoryCode <= Xian80D3_45srid)
{
ObjectProjection = EnumProjectionDatum.Xian80;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = true;
}
//Xian1980 6 degree
if (ObjectFactoryCode >= Xian80D6_13srid && ObjectFactoryCode <= Xian80D6_23srid)
{
ObjectProjection = EnumProjectionDatum.Xian80;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = true;
}
//-------------------------------------------------------------------------------------
////高斯投影小数情况分析
int BJ54D3_25srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Beijing1954_3_Degree_GK_CM_75E;
int BJ54D3_45srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Beijing1954_3_Degree_GK_CM_135E;
int BJ54D6_13srid_CM = (int)esriSRProjCSType.esriSRProjCS_Beijing1954GK_13N;
int BJ54D6_23srid_CM = (int)esriSRProjCSType.esriSRProjCS_Beijing1954GK_23N;
int Xian80D3_25srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_3_Degree_GK_CM_75E;
int Xian80D3_45srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_3_Degree_GK_CM_135E;
int Xian80D6_13srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_GK_CM_75E;
int Xian80D6_23srid_CM = (int)esriSRProjCS4Type.esriSRProjCS_Xian1980_GK_CM_135E;
//源文件分析
//BeiJing1954 3 degree
if (SourceFactoryCode >= BJ54D3_25srid_CM && SourceFactoryCode <= BJ54D3_45srid_CM)
{
SourceProjection = EnumProjectionDatum.Bejing54;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = false;
}
//BeiJing1954 6 degree
if (SourceFactoryCode >= BJ54D6_13srid_CM && SourceFactoryCode <= BJ54D6_23srid_CM)
{
SourceProjection = EnumProjectionDatum.Bejing54;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = false;
}
//Xian1980 3 degree
if (SourceFactoryCode >= Xian80D3_25srid_CM && SourceFactoryCode <= Xian80D3_45srid_CM)
{
SourceProjection = EnumProjectionDatum.Xian80;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = false;
}
//Xian1980 6 degree
if (SourceFactoryCode >= Xian80D6_13srid_CM && SourceFactoryCode <= Xian80D6_23srid_CM)
{
SourceProjection = EnumProjectionDatum.Xian80;
SourceStrip = EnumStrip.Strip3;
IsBigNumber_Source = false;
}
//目标文件分析
//BeiJing1954 3 degree
if (ObjectFactoryCode >= BJ54D3_25srid_CM && ObjectFactoryCode <= BJ54D3_45srid_CM)
{
ObjectProjection = EnumProjectionDatum.Bejing54;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = false;
}
//BeiJing1954 6 degree
if (ObjectFactoryCode >= BJ54D6_13srid_CM && ObjectFactoryCode <= BJ54D6_23srid_CM)
{
ObjectProjection = EnumProjectionDatum.Bejing54;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = false;
}
//Xian1980 3 degree
if (ObjectFactoryCode >= Xian80D3_25srid_CM && ObjectFactoryCode <= Xian80D3_45srid_CM)
{
ObjectProjection = EnumProjectionDatum.Xian80;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = false;
}
//Xian1980 6 degree
if (ObjectFactoryCode >= Xian80D6_13srid_CM && ObjectFactoryCode <= Xian80D6_23srid_CM)
{
ObjectProjection = EnumProjectionDatum.Xian80;
ObjectStrip = EnumStrip.Strip3;
IsBigNumber_Object = false;
}
//------------------------------------------------------------------------------
//KDFeatureClass sfc = new KDFeatureClass_TYPoint(SourceFeatureClass);
//KDFeatureClass ofc = new KDFeatureClass_TYPoint(ObjectFeatureClass);
ZhFeatureClass sfc = new ZhPointFeatureClass(SourceFeatureClass);
ZhFeatureClass ofc = new ZhPointFeatureClass(ObjectFeatureClass);
//KDFeature[] FeatArray = sfc.getSelectedFeaturesByQueryFilter(null);
//KDFeature sfeat = null;
//KDFeature ofeat = null;
ZhFeature[] FeatArray = sfc.getSelectedFeaturesByQueryFilter(null);
ZhFeature sfeat = null;
ZhFeature ofeat = null;
IGeometry geo = null;
IPoint p = null;
/*frmProgressBar1 pb = new frmProgressBar1();
* pb.Text = "坐标投影转换";
* pb.label1.Text = "总数：" + FeatArray.Length.ToString();
* pb.progressBar1.Maximum = FeatArray.Length + 1;
* pb.progressBar1.Value = 0;
* pb.Show();*/
for (int i = 0; i < FeatArray.Length; i++)
{
/*if (pb.progressBar1.Value < pb.progressBar1.Maximum)
* {
* pb.progressBar1.Value += 1;
* pb.label1.Text = "第[" + pb.progressBar1.Value.ToString() + "]个/总数[" + FeatArray.Length.ToString() + "]";
* pb.Refresh();
* }*/
sfeat = FeatArray[i];
//ofeat = ofc.CreateKDFeature();
ofeat = ofc.CreateFeature();
//copy 源字段的值到目标字段中
sfeat.CopyField(ref ofeat);
geo = sfeat.pFeature.ShapeCopy;
//geo.Project(ObjectSpatialReference);
////坐标投影转换
IGeometry ObjGeometry = null;
if (geo is IPoint)
{ //点
double x = 0;
double y = 0;
double B = 0;
double L = 0;
x = (geo as IPoint).X;
y = (geo as IPoint).Y;
this.GetBLFromXY(y, x, ref B, ref L, 0, SourceProjection, SourceStrip, L0, IsBigNumber_Source);
this.GetXYFromBL(B, L, ref x, ref y, 0, ObjectProjection, ObjectStrip, IsBigNumber_Object);
IPoint objP = new PointClass();
objP.X = y;
objP.Y = x;
objP.SpatialReference = ObjectSpatialReference;
//投影到目标点对象
ObjGeometry = objP as IGeometry;
}
else if (geo is IPolyline)
{ //线
IPolyline objPolyline = new PolylineClass();
IPointCollection objPc = objPolyline as IPointCollection;
double x = 0;
double y = 0;
double B = 0;
double L = 0;
IPointCollection pcol = geo as IPointCollection;
object miss = Type.Missing;
for (int j = 0; j < pcol.PointCount; j++)
{
p = pcol.get_Point(j);
x = p.X;
y = p.Y;
this.GetBLFromXY(y, x, ref B, ref L, 0, SourceProjection, SourceStrip, L0, IsBigNumber_Source);
this.GetXYFromBL(B, L, ref x, ref y, 0, ObjectProjection, ObjectStrip, IsBigNumber_Object);
IPoint objPoint = new PointClass();
objPoint.X = y;
objPoint.Y = x;
objPoint.SpatialReference = ObjectSpatialReference;
objPc.AddPoint(objPoint, ref miss, ref miss);
}
(objPc as ITopologicalOperator).Simplify();
ObjGeometry = objPc as IGeometry;
}
else if (geo is IPolygon)
{ //面
IPolygon objPolygon = new PolygonClass();
IGeometryCollection objPc = objPolygon as IGeometryCollection;
double x = 0;
double y = 0;
double B = 0;
double L = 0;
IGeometryCollection GeoCol = geo as IGeometryCollection;
object miss = Type.Missing;
IGeometry tpGeo = null;
IPointCollection pcol = null;
IRing newRing = null;
IPointCollection newRingPointColl = null;
for (int j = 0; j < GeoCol.GeometryCount; j++)
{
tpGeo = GeoCol.get_Geometry(j); //面内环ring(内/外环)
pcol = tpGeo as IPointCollection;
newRing = new RingClass();
newRingPointColl = newRing as IPointCollection;
for (int k = 0; k < pcol.PointCount; k++)
{
p = pcol.get_Point(k);
x = p.X;
y = p.Y;
this.GetBLFromXY(y, x, ref B, ref L, 0, SourceProjection, SourceStrip, L0, IsBigNumber_Source);
this.GetXYFromBL(B, L, ref x, ref y, 0, ObjectProjection, ObjectStrip, IsBigNumber_Object);
IPoint objPoint = new PointClass();
objPoint.X = y;
objPoint.Y = x;
newRingPointColl.AddPoint(objPoint, ref miss, ref miss);
}
newRing.SpatialReference = ObjectSpatialReference;
objPc.AddGeometry(newRing as IGeometry, ref miss, ref miss);
}
(objPc as IGeometry).SpatialReference = ObjectSpatialReference;
(objPc as ITopologicalOperator).Simplify();
ObjGeometry = objPc as IGeometry;
}
else
{
//注记 暂未写
}
ofeat.pFeature.Shape = ObjGeometry;
ofeat.SaveFeature();
}
/*pb.Close();
* pb.Dispose();
* pb = null;zk*/
}
