Gdal 之 dxf转geojson (附加坐标转换)

问题产生

  • 公司的webgis平台需要对接cad数据,因为很多客户的数据是使用cad产生的。
  • 绘制cad时采用的坐标系不尽相同,有的bj54 有的xian80 有的cgcs 2000 ,再修改一下北(东)向加常数、*经。
  • webgis采用cesium开发,坐标系框架为 wgs84

思路

  • gdal 读取dxf文件,循环遍历feature (dwg支持很弱,最好封装成只支持dxf)
  • gdal 坐标转换
  • gdal 将feature存入geojson

code

这里只记录一下demo,没有对坐标参数,dto等进行封装

GdalBase.ConfigureAll();
Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");//支持中文路径
Gdal.SetConfigOption("DXF_ENCODING", "UTF-8");

//创建北京54坐标系
var bj54 = new SpatialReference("");
bj54.ImportFromEPSG(4214);//4214
var clong = 121 + 35 / 60.0;
bj54.SetTM(0, clong, 2.0, 50825, -342129);
bj54.SetTOWGS84(-143.309, -238.731, -83.562, -3.49779, 4.016604, -3.591758, -4.44782521);

var wgs84 = new SpatialReference("");
wgs84.ImportFromEPSG(4326);

var wgs84Geo = wgs84.CloneGeogCS();

//创建转换
var coordinateTransformation = new CoordinateTransformation(bj54, wgs84Geo);

//读取文件转化为DataSource
var ds = Ogr.Open(filePath, 0);
if (ds == null)
{
    Console.WriteLine("打开文件失败");
    return;
}

var extfile = "test.geojson";

if (File.Exists(extfile))
    File.Delete(extfile);

//获取geojson驱动
var geojsonDriver = Ogr.GetDriverByName("GeoJSON");
var geojsonDataSource = geojsonDriver.CreateDataSource(extfile, null);
var geojsonLayer = geojsonDataSource.CreateLayer("JZD", wgs84Geo, wkbGeometryType.wkbUnknown, null);

//读取第一个图层,这个图层不是cad中的layer
var oLayer = ds.GetLayerByIndex(0);
if (oLayer == null)
{
    Console.WriteLine("layer 获取失败");
    return;
}

//设置cad layer 筛选
//oLayer.SetAttributeFilter("Layer = \'JZD\'");

Feature feature;
while ((feature = oLayer.GetNextFeature()) != null)
{
    //获取空间属性
    var geometry = feature.GetGeometryRef();

    //转换
    geometry.Transform(coordinateTransformation);

    //注意,转换为后时 纬度 经度的顺序 ,如果需要经度在前,需要交换xy
    geometry.SwapXY();
    geometry.CloseRings();
    var ret = geojsonLayer.CreateFeature(feature);
}

上一篇:MaskRCNN网络结构


下一篇:使用逻辑回归对鸢尾花进行分类