最近在接触地图数据转换的东西,从硬件kml的数据转换其他的格式,因为从没做过着东西, 先去了解kml文件格式

http://baike.baidu.com/view/400307.htm?fr=aladdin

原来是谷歌地球的一种数据格式,先大致看了写里面的标签,kml与xml文件差多,

kml数据:

longitude:经度  latitude:纬度   altitude:高度 其他的就自己去看了。

现在开始转换,在谷歌上找到一个开源的项目 里面是对地图数据的操作,GDAL

http://www.gdal.org/ogr/ogr_formats.html 参考这里。

kml->SHP

 private string CreateShp(List<Placemark> list)
{ //注册Ogr库
string pszDriverName = DriverType.Shapefile;
OSGeo.OGR.Ogr.RegisterAll(); //为了支持中文路径,请添加下面这句代码
OSGeo.GDAL.Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");
// 为了支持shp属性表字段支持中文,请添加下面这句
OSGeo.GDAL.Gdal.SetConfigOption("SHAPE_ENCODING", " "); //调用对Shape文件读写的Driver接口
OSGeo.OGR.Driver poDriver = OSGeo.OGR.Ogr.GetDriverByName(pszDriverName);
if (poDriver == null)
return "打开驱动失败"; //用此Driver创建Shape文件
OSGeo.OGR.DataSource poDS;
poDS = poDriver.CreateDataSource(CreateSavePath(), null);
if (poDS == null)
return "创建数据源失败"; //创建层Layer
OSGeo.OGR.Layer poLayer;
poLayer = poDS.CreateLayer(CurrentDate, null, OSGeo.OGR.wkbGeometryType.wkbPoint, null);
if (poLayer == null)
return "创建地图层失败"; //创建属性
OSGeo.OGR.FieldDefn oField = new OSGeo.OGR.FieldDefn("name", OSGeo.OGR.FieldType.OFTString);
oField.SetWidth();
OSGeo.OGR.FieldDefn oField2 = new OSGeo.OGR.FieldDefn("height", OSGeo.OGR.FieldType.OFTInteger); OSGeo.OGR.FieldDefn ofield3 = new FieldDefn("lat", FieldType.OFTInteger);
OSGeo.OGR.FieldDefn ofield4 = new FieldDefn("lng", FieldType.OFTInteger); poLayer.CreateField(oField, );
poLayer.CreateField(oField2, );
poLayer.CreateField(ofield3, );
poLayer.CreateField(ofield4, ); //创建一个Feature,一个Point
OSGeo.OGR.Feature poFeature = new Feature(poLayer.GetLayerDefn());
OSGeo.OGR.Geometry pt = new Geometry(OSGeo.OGR.wkbGeometryType.wkbPoint); foreach (Placemark item in list)
{
//属性一"名称"赋值
poFeature.SetField(, item.Name);
//属性二"高度"赋值
poFeature.SetField(, item.LookAt.altitude);
poFeature.SetField(, item.LookAt.latitude);
poFeature.SetField(, item.LookAt.longitude); //添加坐标点 x y z
pt.AddPoint(item.LookAt.longitude, item.LookAt.latitude, item.LookAt.altitude); poFeature.SetGeometry(pt);
//将带有坐标及属性的Feature要素点写入Layer中
poLayer.CreateFeature(poFeature); } //关闭文件读写
poFeature.Dispose();
poDS.Dispose(); return "转换成功";
}

SHP->mapinfo tab转换出现问题,当我的经纬度像上图一样小数位数比较多的情况,会自动截断,仔细调试都没找到解决办法,截断都是没规则的,只能继续谷歌了,╮(╯▽╰)╭。

代码如下:

 public string Convert(string driverTypeName, OnAction action)
{
string msg = string.Empty;
DataSource sourceSource;
OSGeo.OGR.Driver shpDriver;
RegisterAll(driverTypeName, out sourceSource, out shpDriver);
DataSource destSource = shpDriver.CreateDataSource(Dest, new string[] {
"AUX=YES",
"STATISTICS=YES"}); //中文
int layerCount = sourceSource.GetLayerCount();
for (int i = ; i < layerCount; i++)
{
Layer layer = sourceSource.GetLayerByIndex(i);
int featureCount = layer.GetFeatureCount();
Layer destLayer = null; //深度拷贝
//Layer destLayer = destSource.CopyLayer(layer, dest, null); #region MyRegion
for (int j = ; j < featureCount; j++)
{
Feature feature = layer.GetFeature(j); if (feature != null)
{
try
{ if (destLayer == null)
{
wkbGeometryType geoType = feature.GetGeometryRef().GetGeometryType();
//创建图层
destLayer = destSource.CreateLayer(
layer.GetName(),
action(layer),
geoType,
new string[] { });
//创建字段
FeatureDefn featureDefn = layer.GetLayerDefn();
for (int k = ; k < featureDefn.GetFieldCount(); k++)
{
destLayer.CreateField(featureDefn.GetFieldDefn(k), );
}
}
//写入要素 Feature cloneFeature = feature.Clone(); // Feature newfeature = ConvetToFeature(feature); destLayer.CreateFeature(cloneFeature); // OnFeatureConvert(featureCount, EventArgs.Empty);
}
catch (Exception ex)
{
msg = "转换失败"+ex.Message;
continue;
}
}
}
#endregion //保存
destLayer.SyncToDisk();
} msg = "转换成功";
sourceSource.Dispose();
destSource.Dispose();
shpDriver.Dispose();
return msg;
}

原来是坐标系的问题,

shp的默认与tab的坐标系不一样。参考文章:http://blog.sina.com.cn/s/blog_6e51df7f0100ui7n.html

坐标系转换参考:http://wiki.woodpecker.org.cn/moin/lilin/ogr-create

终于解决。

KML,SHP TAB互转,GDAL的更多相关文章

  1. C#、C++用GDAL读shp文件(转载)

    C#.C++用GDAL读shp文件 C#用GDAL读shp文件 (2012-08-14 17:09:45) 标签: 杂谈 分类: c#方面的总结 1.目前使用开发环境为VS2008+GDAL1.81 ...

  2. Android GIS开发系列-- 入门季(13)Gdal简单写个shp文件

    Gdal是用来读写栅格与矢量数据的,在Gdal官网,可以下载相关的资源进行平台的编译.其实Arcgis底层也是用Gdal来读取shp文件的,那在Android中可以直接读写shp文件吗,是可以的.这里 ...

  3. Java 使用GDAL 读写 shapefile

    读取shp文件,并把它转化为json import org.gdal.ogr.*; import org.gdal.ogr.Driver; import org.gdal.gdal.*; public ...

  4. shp的基本操作

    本节将介绍如何利用python完成对shp的基本操作 1.读取shp四至 import shapefile sf = shapefile.Reader(r"E:\shp\1.shp" ...

  5. webapi_uploadfile_gdal_to_geojson_and_unzipfile

    using ICSharpCode.SharpZipLib.Zip; using OSGeo.GDAL; using OSGeo.OGR; using System; using System.Col ...

  6. python使用游标访问数据

    游标是一种数据访问对象,可用于在表中迭代一组行或者向表中插入新行.游标有三种形式:搜索.插入或更新.游标通常用于读取现有几何和写入新几何. 每种类型的游标均由对应的 ArcPy 函数(SearchCu ...

  7. SublimeCodeIntel代码自动补全配置

    主要使用python3,所有配置以python3为例.其他语言同理.利用sublimeCodeIntel插件可以实现自动提示python3代码.跳转追踪自定义函数.查看系统函数等.功能还是相当强大的. ...

  8. gdal集成kml库的做法

    作者:朱金灿 来源:http://blog.csdn.net/clever101 最近要读取kml文件,具体就是把kml文件当作一个矢量文件来读取.我发现gdal是支持集成kml库的.不过集成这个km ...

  9. GDAL集成对KML文件的支持

    目录 1. 正文 1.1. 编译LibKML 1.1.1. 第三方库支持 1.1.2. 编译错误 1.2. 配置GDAL 1.3. 链接问题 2. 参考 1. 正文 GDAL可以支持将KML作为矢量文 ...

随机推荐

  1. js 发送短信倒计时、秒杀倒计时实现代码

    <!doctype html> <html> <head> <meta charset="utf-8"> <meta name ...

  2. vue启动时候报错

    使用vue时,在已经安装模块完毕的情况下,依旧会报错,如: Module build failed: Error: %1 is not a valid Win32 application. 这个时候只 ...

  3. rank() partition by 排名次

    rank()排名 partition by分组与group by相比各有优势,在这里就省略100字.... 以下为案例: create table student -- 学生表(sid integer ...

  4. day20 hashlib、hmac、subprocess、configparser模块

    hashlib模块:加密 import hashlib# 基本使用cipher = hashlib.md5('需要加密的数据的二进制形式'.encode('utf-8'))print(cipher.h ...

  5. 重大变革即将来临 5G CPE会替代光纤入户吗?

    导读 从国内的新闻报道上我们可以看到,从2018年下半年开始各大重要活动.春晚直播等,都宣布已经使用5G网络.既然支持5G网络的终端都还没有正式上市,那么5G网络是如何使用的呢?答案是5G CPE设备 ...

  6. Recovering Low-Rank Matrices From Few Coefficients In Any Basis

    目录 引 主要结果 定理2,3 定理4 直观解释 Recovering Low-Rank Matrices From Few Coefficients In Any Basis-David Gross ...

  7. 微信内分享第三方H5链接无法使用内置浏览器打开的解决方案

    很多朋友在微信内想分享转发H5链接的时候都会很容易碰到H5链接在微信内无法打开或在微信内无法打开app下载页的情况.通常这种情况微信会给个提示 “已停止访问该网址” ,那么导致这个情况的因素有哪些呢, ...

  8. (七)jdk8学习心得之join方法

    七.join方法 1. 作用:将list或者数组按照连接符进行连接,返回一个字符串. 2. 使用方法 1) String.join(“连接符”,数组对象或者list对象) 2) 首先转换成stream ...

  9. 【MySQL 读书笔记】普通索引和唯一索引应该怎么选择

    通常我们在做这个选择的时候,考虑得最多的应该是如果我们需要让 Database MySQL 来帮助我们从数据库层面过滤掉对应字段的重复数据我们会选择唯一索引,如果没有前者的需求,一般都会使用普通索引. ...

  10. 南邮 base64全家桶

    这几天不想学逆向 做做crypto(菜还瞎j2做)..... 题目: 全家桶全家桶全家桶!我怎么饿了......密文(解密前删除回车):R1pDVE1NWlhHUTNETU4yQ0dZWkRNTUpY ...