ArcGis 获取地理、平面坐标系

 1 IProjectedCoordinateSystem LoadProjectedCoordinateSystem( string shpfile)
 2         {
 3              /*   IFeatureClass从中获取地理坐标、平面坐标
 4              string shpDir = System.IO.Path.GetDirectoryName(shpfile);
 5              string shpfilename = System.IO.Path.GetFileName(shpfile);
 6              IWorkspaceFactory wsFactory = new ShapefileWorkspaceFactoryClass();
 7              IWorkspace ws = wsFactory.OpenFromFile(shpDir, 0);
 8              IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)ws;
 9              IFeatureClass featureCls = featureWorkspace.OpenFeatureClass(shpfilename);
10              IGeoDataset pGeoDataset = featureCls as IGeoDataset;
11              var spatialReference = pGeoDataset.SpatialReference;
12 
13              IProjectedCoordinateSystem pProCoordSys = spatialReference as IProjectedCoordinateSystem;//平面
14              IGeographicCoordinateSystem pGeographCoordinateSys = pProCoordSys.GeographicCoordinateSystem;//地理
15               */
16 
17 
18              // 1、从文件中获取
19              ESRI.ArcGIS.Geometry.ISpatialReferenceFactory spatialRefFactory =  new ESRI.ArcGIS.Geometry.SpatialReferenceEnvironmentClass();
20             ESRI.ArcGIS.Geometry.ISpatialReference spatialReference = spatialRefFactory.CreateESRISpatialReferenceFromPRJFile(shpfile.ToLower().Replace( " .shp ",  " .prj "));
21 
22              // 2、自定义
23              IProjectedCoordinateSystem geographicCoordinateSystem = spatialRefFactory.CreateProjectedCoordinateSystem(( int)ESRI.ArcGIS.Geometry.esriSRProjCSType.esriSRProjCS_World_Mercator);
24 
25              // 3、从文件中获取
26              Type factoryType = Type.GetTypeFromProgID( " esriGeometry.SpatialReferenceEnvironment ");
27             System.Object obj = Activator.CreateInstance(factoryType);
28             ISpatialReferenceFactory spatialReferenceFactory = obj  as ISpatialReferenceFactory;
29 
30             IProjectedCoordinateSystem projectedCoordinateSystem = spatialReferenceFactory.CreateESRISpatialReferenceFromPRJFile(shpfile.ToLower().Replace( " .shp ",  " .prj "))  as IProjectedCoordinateSystem;
31              return projectedCoordinateSystem;
32 
33         }