获取并修改图层投影坐标系

10 篇文章 0 订阅
1 篇文章 0 订阅
//获取已知原始坐标系
IFeatureClass yFeatureClass=null;
ISpatialReference pSpatialRef = (yFeatureClass as IGeoDataset).SpatialReference;

//修改现在要素的坐标系,继承原坐标系
IFeatureClass  xFeatureClass=null;
IGeoDataset pGeoDataset =xFeatureClass as IGeoDataset;
IGeoDatasetSchemaEdit pGeoDatasetSchemaEdit = pGeoDataset as IGeoDatasetSchemaEdit;
//空间坐标系可改变
if (pGeoDatasetSchemaEdit.CanAlterSpatialReference == true)
 {
      创建SpatialReferenceEnvironmentClass对象
      //ISpatialReferenceFactory2 pSpaRefFactory = new SpatialReferenceEnvironmentClass();
      创建地理坐标系对象
      //IGeographicCoordinateSystem pNewGeoSys = pSpaRefFactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);//4214代表Beijing1954;
      //pGeoDatasetSchemaEdit.AlterSpatialReference(pNewGeoSys);
      pGeoDatasetSchemaEdit.AlterSpatialReference(pSpatialRef);
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要使用`osgeo.osr`模块从矢量图层获取原始坐标系,并将其转换为投影坐标系,可以按照以下步骤操作: ```python from osgeo import ogr, osr # 打开矢量图层 shapefile_path = 'your_shapefile.shp' driver = ogr.GetDriverByName('ESRI Shapefile') dataSource = driver.Open(shapefile_path, 0) layer = dataSource.GetLayer() # 获取原始坐标系 sourceSpatialRef = layer.GetSpatialRef() # 定义目标坐标系(例如 EPSG:3857) targetSpatialRef = osr.SpatialReference() targetSpatialRef.ImportFromEPSG(3857) # 创建坐标转换对象 transform = osr.CoordinateTransformation(sourceSpatialRef, targetSpatialRef) # 创建投影图层 projShapefile = 'projected_shapefile.shp' projDriver = ogr.GetDriverByName('ESRI Shapefile') projDataSource = projDriver.CreateDataSource(projShapefile) projLayer = projDataSource.CreateLayer('projected_layer', targetSpatialRef, geom_type=ogr.wkbPolygon) # 遍历原始图层的要素,进行投影转换并添加到投影图层 feature = layer.GetNextFeature() while feature: geometry = feature.GetGeometryRef() geometry.Transform(transform) projFeature = ogr.Feature(projLayer.GetLayerDefn()) projFeature.SetGeometry(geometry) projLayer.CreateFeature(projFeature) feature = layer.GetNextFeature() # 释放资源 dataSource.Destroy() projDataSource.Destroy() ``` 在这个示例中,我们首先使用`ogr`模块打开矢量图层文件,并获取图层对象。然后,我们使用`GetSpatialRef`方法获取原始坐标系。 接下来,我们定义了目标坐标系(例如EPSG:3857),并使用`osr.SpatialReference`创建了一个目标坐标系对象。 然后,我们使用`osr.CoordinateTransformation`创建了一个坐标转换对象,用于将原始坐标系转换为目标坐标系。 接着,我们创建了一个新的投影图层,并使用`CreateFeature`方法遍历原始图层的要素。在循环中,我们获取要素的几何对象,并使用转换对象对其进行投影转换。然后,我们创建了一个新的投影要素,并将投影后的几何对象添加到投影图层中。 最后,我们释放了资源并关闭文件。 请确保你已经安装了GDAL库,并将代码中的`your_shapefile.shp`替换为你自己的Shapefile文件路径,并将`projected_shapefile.shp`替换为你希望保存投影结果的文件路径。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值