文档章节

node-gdal实践-创建shp文件

遥想公瑾当年
 遥想公瑾当年
发布于 2017/06/01 11:12
字数 771
阅读 10
收藏 0

一 gdal.open解析

gdal.open.png 阅读api可知,gdal.open方法用于创建和打开数据集dataset,如果以‘w’模式的话,会将更改从内存提交磁盘,否则,更改将被丢弃无效。 path:创建数据集的路径,如“e:/shp”。 mode:"r","r+","w"三种模式。 drivers:数据集的驱动,字符或数组类型。 那么问题来了,如果我要创建一个shp,path和mode都很容易确定,而drivers是个什么鬼?查询api文档查不到,于是,使用以下代码,用于浏览下到底有哪些驱动:


	var gdal=require('gdal'); 
        gdal.drivers.forEach(function(drive,i){
		console.log(drive.description);
	})

输出结果如下:

VRT
GTiff
NITF
RPFTOC
ECRGTOC
HFA
SAR_CEOS
CEOS
JAXAPALSAR
GFF
ELAS
AIG
AAIGrid
GRASSASCIIGrid
SDTS
DTED
PNG
JPEG
MEM
JDEM
ESAT
XPM
BMP
DIMAP
AirSAR
RS2
PCIDSK
ILWIS
SGI
SRTMHGT
Leveller
Terragen
ISIS3
ISIS2
PDS
VICAR
TIL
ERS
L1B
FIT
RMF
RST
INGR
GSAG
GSBG
GS7BG
COSAR
TSX
COASP
R
MAP
KMLSUPEROVERLAY
PNM
DOQ1
DOQ2
GenBin
PAux
MFF
MFF2
FujiBAS
GSC
FAST
BT
LAN
CPG
IDA
NDF
EIR
DIPEx
LCP
GTX
LOSLAS
NTv2
CTable2
ACE2
SNODAS
KRO
ROI_PAC
ENVI
EHdr
ISCE
USGSDEM
NWT_GRD
NWT_GRC
ADRG
SRP
BLX
SAGA
XYZ
HF2
CTG
E00GRID
ZMap
NGSGEOID
IRIS
ESRI Shapefile
MapInfo File
UK .NTF
OGR_SDTS
S57
DGN
OGR_VRT
REC
Memory
BNA
CSV
GML
GPX
KML
GeoJSON
OGR_GMT
WAsP
OpenFileGDB
XPlane
DXF
Geoconcept
GeoRSS
GPSTrackMaker
PGDUMP
GPSBabel
SUA
OpenAir
OGR_PDS
HTF
AeronavFAA
EDIGEO
SVG
Idrisi
SEGUKOOA
SEGY
SXF
AVCBin
AVCE00

注意,输出的驱动描述有ESRI Shapefile和OpenFileGDB,分别对应esri的shp和gdb。通过描述可知,新建一个shp的构造函数应该如下: var dataset=gdal.open('b.shp',"w",'ESRI Shapefile');

二 创建shp

2.1 创建图层

dataset是由layers组成的,我们创建了一个shp的数据集,于是,在该数据集下创建layer,查看layers.create的api如下:

create( name,srs,geomType,creation_options) 
新增图层。
参数:
name String类型,图层名称。
srs gdal.SpatialReference| Null,图层投影坐标系。
geomType Integer | FunctionGeometry,图层几何类型。
creation_options String[] | Object 特殊选项,可选。
返回值:
gdal.Layer
示例:
dataset.layers.create('layername', null, gdal.Point);

创建点图层示例如下: dataset.layers.create('b',null, gdal.Point);

2.2 图层坐标系

新建的图层,坐标系是null,我们希望建立的shp图层坐标系是4326的wgs84坐标,上文创建点图层如下:

    var ref = gdal.SpatialReference.fromEPSGA(4326);
	dataset.layers.create('b',ref, gdal.Point);

2.3 新增字段

    var layer = dataset.layers.get(0);
	layer.fields.add(new gdal.FieldDefn('name', gdal.OFTString));
	var feature = new gdal.Feature(layer);
	//给字段赋值
	feature.fields.set('name', '测试');
	//构造图形并赋值
	var pt = new gdal.Point(118.5,32.1)
	
	feature.setGeometry(pt);
	layer.features.add(feature);
	layer.flush();

运行之后,发现name值“测试”是中文乱码: result.png 理所当然,我们认为是字符集的问题。 2.4 设置shp字符集

    //GBK的code page是CP936
	gdal.config.set('SHAPE_ENCODING','CP936');

本来设置属性是gbk,报错。非要设置gbk的codepage代码,涨姿势了。

三 结果

完整代码如下:

	var gdal=require('gdal'); 
	//GBK的code page是CP936
	gdal.config.set('SHAPE_ENCODING','CP936');
	var dataset=gdal.open('b.shp',"w",'ESRI Shapefile');
	var ref = gdal.SpatialReference.fromEPSGA(4326);
	dataset.layers.create('b',ref, gdal.Point);
	var layer = dataset.layers.get(0);
	layer.fields.add(new gdal.FieldDefn('name', gdal.OFTString));
	var feature = new gdal.Feature(layer);
	//给字段赋值
	feature.fields.set('name', '测试');
	//构造图形并赋值
	var pt = new gdal.Point(118.5,32.1)
	feature.setGeometry(pt);
	layer.features.add(feature);
	layer.flush();

可视化结果如下:

b.shp.png

后人看也许很简单,因为前人填了很多坑,无论坐标系还是模式还是中文乱码,是耗尽心血点滴测试出来的,且看且珍惜啊。

© 著作权归作者所有

共有 人打赏支持
遥想公瑾当年
粉丝 6
博文 11
码字总数 17519
作品 0
南京
高级程序员
私信 提问
GDAL不支持SHP文件中的中文字符

我最近使用了一个别人的用GDAL编写的软件,是国外的,软件功能是将SHP文件中的部分属性字段导入到CSV文件中,还有就是简化SHP文件边界点,并提取每一个节点坐标和质心坐标,但在使用中不支持...

dr_gis_jyx
2014/11/15
869
1
ArcGIS Runtime加载本地Shape数据和编辑Shape数据

问题的提出 ArcGIS Runtime 支持MPK,APK,TPK,GPK这些打包之后的数据,这些也都是ArcGIS 10.1 中新增的数据格式(MPK除外,MPK是在ArcGIS 10中增加的,但是并不能支持ArcGIS Runtime),我们经...

长平狐
2012/11/28
3.5K
0
用C++语言 利用GDAL读写遥感多波段的时候,只能用读取单波段的信息,多波段功能实现不了,我用的是多波段数据,想读取部分波段信息

#include "gdal_priv.h" include"gdal.h" include include typedef unsigned char BYTE; using namespace std; int main() { GDALAllRegister();//注册所有驱动 GDALDataset* poDataset1; /......

毛毛真nice
11/04
0
0
​python版本gdal中在py文件中创建的栅格影像如何正常关闭并使得其可以删除

python版本gdal中在py文件中创建的栅格影像如何正常关闭并使得其可以删除 测试py文件代码: #coding:gb2312 import os from osgeo import gdal from osgeo import ogr from osgeo.gdalconst...

funnybone
2013/11/15
587
0
简析将shp导入Oracle并利用geoserver将导入的数据发布

文章版权由作者李晓晖和博客园共有,若转载请于明显处标明出处:http://www.cnblogs.com/naaoveGIS/。 1.环境准备 1.1 软件准备 首先要安装有支持空间数据的Oracle,其次有安装版或免安装版的...

李晓晖
2015/01/02
0
0

没有更多内容

加载失败,请刷新页面

加载更多

颜色模型与颜色应用---光的特性

电磁频谱 颜色的心理学特征

中国龙-扬科
17分钟前
2
0
android音频及强噪相关

Android AudioRecord和AudioTrack介绍(Android音频收集和播放 麦克风降噪) https://blog.csdn.net/tanningzhong/article/details/72844559 音频采集(AudioRecorder) https://www.jianshu.......

whoisliang
24分钟前
2
0
redis-持久化

RDB rdb持久化是把当前进程数据生成快照保存到磁盘的过程。触发RDB持久化过程分为手动触发和自动触发。 触发机制 bgsave执行流程 rdb优缺点 AOF 记录每次写命令,重启时再重新执行aof文件中的...

grace_233
28分钟前
3
0
电话激活Windows 中文操作系统步骤

已购买微软MAK批量授权,系统又在企业内网中,无法通过连接Internet进行激活,还可以通过电话完成激活。 前期准备 请提前准备好产品密钥,Product Key格式如下:XXXXX-XXXXX-XXXXX-XXXXX-XXX...

tonyfox
30分钟前
4
0
Apache用户认证,域名跳转,访问日志

Apache用户认证 当设置了用户认证后,用户访问网站时,需要输入用户名和密码才能访问。 可以全局设置,也可以为某几个虚拟主机单独配置。 下面以全局配置进行操作示例。 编辑httpd.conf进行配...

野雪球
35分钟前
1
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部