当前位置：文档之家 › shapefile格式说明及读写代码示例 收藏

shapefile格式说明及读写代码示例 收藏

  • 格式：doc
  • 大小：76.00 KB
  • 文档页数：17

shapefile格式说明及读写代码示例收藏Shape files数据说明Shape files是ESRI提供的一种矢量数据格式，它没有拓扑信息，一个Shape files由一组文件组成，其中必要的基本文件包括坐标文件（.shp）、索引文件（.shx）和属性文件（.dbf）三个文件。

坐标文件的结构说明坐标文件(.shp)用于记录空间坐标信息。

它由头文件和实体信息两部分构成（如图2.1所示）。

坐标文件的文件头坐标文件的文件头是一个长度固定(100 bytes)的记录段，一共有9个int型和7个double型数据，主要记录内容见表2.2。

文件头记录头记录内容记录头记录内容记录头记录内容记录头记录内容…………记录头记录内容图2.1 坐标文件的结构起始位置名称数值类型位序0 File Code 9994 Integer big4 Unused 0 Integer big8 Unused 0 Integer big12 Unused 0 Integer big16 Unused 0 Integer big20 Unused 0 Integer big24 文件长度文件的实际长度Integer big28 版本号1000 Integer Little32 几何类型表示这个Shapefile文件所记录的空间数据的几何类型Integer Little36 Xmin 空间数据所占空间范围的X方向最小值Double Little44 Ymin 空间数据所占空间范围的Y方向最小值Double Little52 Xmax 空间数据所占空间范围的X方向最大值Double Little60 Ymax 空间数据所占空间范围的Y方向最大值Double Little68* Zmin 空间数据所占空间范围的Z方向最小值Double Little76* Zmax 空间数据所占空间范围的Z方向最大值Double Little84* Mmin 最小Measure值Double Little92* Mmax 最大Measure值Double Little表2.2shapefiles 头文件表注：最后4个加星号特别标示的四个数据只有当这个Shapefile文件包含Z方向坐标或者具有Measure值时才有值，否则为0.0。

所谓Measure值，是用于存储需要的附加数据，可以用来记录各种数据，例如权值、道路长度等信息。

位序细心的读者会注意到表2.2中的数值的位序有Little和big的区别，对于位序是big的数据我们在读取时要小心。

通常，数据的位序都是Little，但在有些情况下可能会是big，二者的区别在于它们位序的顺序相反。

一个位序为big的数据，如果我们想得到它的真实数值，需要将它的位序转换成Little即可。

转换原理非常简单，就是交换字节顺序，下面是作者实现的在两者间进行转换的程序，代码如下：//位序转换程序unsigned long OnChangeByteOrder (int indata){char ss[8];char ee[8];unsigned long val = unsigned long(indata);_ultoa( val, ss, 16 );//将十六进制的数(val)转到一个字符串(ss)中int i;int length=strlen(ss);if(length!=8){for(i=0;i<8-length;i++)ee[i]='0';for(i=0;i<length;i++)ee[i+8-length]=ss[i];for(i=0;i<8;i++)ss[i]=ee[i];}////******进行倒序int t;t =ss[0];ss[0] =ss[6];ss[6] =t;t =ss[1];ss[1] =ss[7];ss[7] =t;t =ss[2];ss[2] =ss[4];ss[4] =t;t =ss[3];ss[3] =ss[5];ss[5] =t;////******//******将存有十六进制数(val)的字符串(ss)中的十六进制数转成十进制数int value=0;for(i=0;i<8;i++){int k;CString mass;mass=ss[i];if(ss[i]=='a' ||ss[i]=='b' ||ss[i]=='c' ||ss[i]=='d' ||ss[i]=='e' ||ss[i]=='f')k=10+ss[i]-'a';elsesscanf(mass,"%d",&k);value=value+int(k*pow(16,7-i));}return (value);}Shapefile文件支持的几何类型（ShapeType）Shapefile文件所支持的几何类型如表2.3所示：编号几何类型0 Null Shape（表示这个Shapefile文件不含坐标）1 Point（表示Shapefile文件记录的是点状目标，但不是多点）3 PolyLine（表示Shapefile文件记录的是线状目标）5 Polygon（表示Shapefile文件记录的是面状目标）8 MultiPoint（表示Shapefile文件记录的是多点，即点集合）11 PointZ（表示Shapefile文件记录的是三维点状目标）13 PolyLineZ（表示Shapefile文件记录的是三维线状目标）15 PolygonZ（表示Shapefile文件记录的是三维面状目标）18 MultiPointZ（表示Shapefile文件记录的是三维点集合目标）21 PointM（表示含有Measure值的点状目标）23 PolyLineM（表示含有Measure值的线状目标）25 PolygonM（表示含有Measure值的面状目标）28 MultiPointM（表示含有Measure值的多点目标）31 MultiPatch（表示复合目标）表2.3shapefiles文件支持的几何类型对于一个不是记录Null Shape类型的Shapefile文件，它所记录的空间目标的几何类型必须一致，不能在一个Shapefile文件中同时记录两种不同类型的几何目标。

读取坐标文件（.shp）的文件头的代码如下：void OnReadShp（CString ShpFileName）{FILE* m_ShpFile_fp; //****Shp文件指针//打开坐标文件if((m_ShpFile_fp=fopen(ShpFileName,"rb"))==NULL){return;//读取坐标文件头的内容开始int FileCode;int Unused;int FileLength;int Version;int ShapeType;double Xmin;double Ymin;double Xmax;double Ymax;double Zmin;double Zmax;double Mmin;double Mmax;fread(&FileCode, sizeof(int), 1,m_ShpFile_fp);FileCode = OnChangeByteOrder(FileCode);for(i=0;i<5;i++)fread(&Unused,sizeof(int), 1,m_ShpFile_fp);fread(&FileLength, sizeof(int), 1,m_ShpFile_fp);FileLength = OnChangeByteOrder(FileLength);fread(&Version, sizeof(int), 1,m_ShpFile_fp);fread(&ShapeType, sizeof(int), 1,m_ShpFile_fp);fread(&Xmin, sizeof(double),1,m_ShpFile_fp);fread(&Ymin, sizeof(double),1,m_ShpFile_fp);fread(&Xmax, sizeof(double),1,m_ShpFile_fp);fread(&Ymax, sizeof(double),1,m_ShpFile_fp);fread(&Zmin, sizeof(double),1,m_ShpFile_fp);fread(&Zmax, sizeof(double),1,m_ShpFile_fp);fread(&Mmin, sizeof(double),1,m_ShpFile_fp);fread(&Mmax, sizeof(double),1,m_ShpFile_fp);//读取坐标文件头的内容结束//根据几何类型读取实体信息……}实体信息的内容实体信息负责记录坐标信息，它以记录段为基本单位，每一个记录段记录一个地理实体目标的坐标信息，每个记录段分为记录头和记录内容两部分。

记录头的内容包括记录号（Record Number）和坐标记录长度(Content Length) 两个记录项。

它们的位序都是big。

记录号（Record Number）和坐标记录长度(Content Length) 两个记录项都是int型，并且shapefile文件中的记录号都是从1开始的。

记录内容包括目标的几何类型（ShapeType）和具体的坐标记录(X、Y) ，记录内容因要素几何类型的不同其具体的内容及格式都有所不同。

下面分别介绍点状目标（Point）、线状目标（PolyLine）和面状目标（Polygon）三种几何类型的.shp文件的记录内容：shapefile中的点状目标由一对X、Y坐标构成，坐标值为双精度型（double）。

点状目标的记录内容如表2.4：记录项数值数据类型长度个数位序几何类型（ShapeType）1（表示点状目标）int型4 1 LittleX方向坐标X方向坐标值double型8 1 LittleY方向坐标Y方向坐标值double型8 1 Little表2.4点状目标的记录内容下面是读取点状目标的记录内容的代码：OnReadPointShp(CString ShpFileName){//打开坐标文件……//读取坐标文件头的内容开始……//读取点状目标的实体信息int RecordNumber;int ContentLength;int num =0;while((fread(&RecordNumber, sizeof(int), 1,ShpFile_fp)!=0)){num++;fread(&ContentLength,sizeof(int), 1,ShpFile_fp);RecordNumber = OnChangeByteOrder(RecordNumber);ContentLength = OnChangeByteOrder(ContentLength);int shapeType;double x;double y;fread(&shapeType, sizeof(int), 1,ShpFile_fp);fread(&x, sizeof(double), 1,ShpFile_fp);fread(&y, sizeof(double), 1,ShpFile_fp);}}线状目标shapefile中的线状目标是由一系列点坐标串构成，一个线目标可能包括多个子线段，子线段之间可以是相离的，同时子线段之间也可以相交。

合集下载

结合C++和GDAL实现shapefile（shp）文件的读取

结合C++和GDAL实现shapefile（shp）文件的读取

//输出每个字段的值 cout<<poFeature->GetFieldAsString(iField)<<" "; } cout<<endl; OGRFeature::DestroyFeature( poFeature ); } GDALClose( poDS ); system("pause"); return 1; }
结果显示：
登录后才能查看或发表评论立即登录或者逛逛博客园首页
结合 C++和 GDAL实现 shapefile（ shp）文件的读取
工具：vs2012+GDAL 2.0
数据：中国省界SHP文件bou2_4p.shp 可
包含头文件：
#include ""
代码：
int main() {
GDALAllRegister(); GDALDataset *poDS; CPLSetConfigOption("SHAPE_ENCODING",""); //解决中文乱码问题 //读取shp文件 poDS = (GDALDataset*) GDALOpenEx("d:/shapefile/bou2_4p.shp", GDAL_OF_VECTOR, NULL, NULL, NULL );
if(poFeature->GetFieldAsDouble("AREA")<1) continue; //去掉面积过小的polygon i=i++; cout<<i<<" "; OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn(); int iField; int n=poFDefn->GetFieldCount(); //获得字段的数目，不包括前两个字段（FID,Shape); for( iField = 0; iField <n; iField++ ) {
Java使用GDAL读写shapefile

Java使用GDAL读写shapefile

Java使⽤GDAL读写shapefile 读取shp⽂件，并把它转化为jsonimport org.gdal.ogr.*;import org.gdal.ogr.Driver;import org.gdal.gdal.*;public class GdalShpTest {public static void main(String[] args) {// 注册所有的驱动ogr.RegisterAll();// 为了⽀持中⽂路径，请添加下⾯这句代码gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES");// 为了使属性表字段⽀持中⽂，请添加下⾯这句gdal.SetConfigOption("SHAPE_ENCODING","");String strVectorFile = "D:\\test\\NODE.shp";//打开⽂件DataSource ds = ogr.Open(strVectorFile,0);if (ds == null){System.out.println("打开⽂件失败！" );return;}System.out.println("打开⽂件成功！" );Driver dv = ogr.GetDriverByName("GeoJSON");if (dv == null){System.out.println("打开驱动失败！" );return;}System.out.println("打开驱动成功！" );dv.CopyDataSource(ds, "D:\\test\\node.json");System.out.println("转换成功！" );}}写shp⽂件import org.gdal.ogr.*;import org.gdal.gdal.*;class writeShp2 {public static void main(String[] args) {writeShp2 readshpObj = new writeShp2();readshpObj.WriteVectorFile();}static void WriteVectorFile() {String strVectorFile = "D:\\test\\test.shp";ogr.RegisterAll();gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");gdal.SetConfigOption("SHAPE_ENCODING", "CP936");String strDriverName = "ESRI Shapefile";org.gdal.ogr.Driver oDriver = ogr.GetDriverByName(strDriverName);if (oDriver == null) {System.out.println(strVectorFile + " 驱动不可⽤！\n");return;}DataSource oDS = oDriver.CreateDataSource(strVectorFile, null);if (oDS == null) {System.out.println("创建⽮量⽂件【" + strVectorFile + "】失败！\n");return;}Layer oLayer = oDS.CreateLayer("TestPolygon", null, ogr.wkbPolygon, null);if (oLayer == null) {System.out.println("图层创建失败！\n");return;}// 下⾯创建属性表// 先创建⼀个叫FieldID的整型属性FieldDefn oFieldID = new FieldDefn("FieldID", ogr.OFTInteger);oLayer.CreateField(oFieldID);// 再创建⼀个叫FeatureName的字符型属性，字符长度为50FieldDefn oFieldName = new FieldDefn("FieldName", ogr.OFTString);oFieldName.SetWidth(100);oLayer.CreateField(oFieldName);FeatureDefn oDefn = oLayer.GetLayerDefn();// 创建三⾓形要素Feature oFeatureTriangle = new Feature(oDefn);oFeatureTriangle.SetField(0, 0);oFeatureTriangle.SetField(1, "三⾓形");Geometry geomTriangle = Geometry.CreateFromWkt("POLYGON ((0 0,20 0,10 15,0 0))");oFeatureTriangle.SetGeometry(geomTriangle);oLayer.CreateFeature(oFeatureTriangle);// 创建矩形要素Feature oFeatureRectangle = new Feature(oDefn);oFeatureRectangle.SetField(0, 1);oFeatureRectangle.SetField(1, "矩形");Geometry geomRectangle = Geometry.CreateFromWkt("POLYGON ((30 0,60 0,60 30,30 30,30 0))");oFeatureRectangle.SetGeometry(geomRectangle);oLayer.CreateFeature(oFeatureRectangle);// 创建五⾓形要素Feature oFeaturePentagon = new Feature(oDefn);oFeaturePentagon.SetField(0, 2);oFeaturePentagon.SetField(1, "五⾓形");Geometry geomPentagon = Geometry.CreateFromWkt("POLYGON ((70 0,85 0,90 15,80 30,65 15,70 0))"); oFeaturePentagon.SetGeometry(geomPentagon);oLayer.CreateFeature(oFeaturePentagon);oDS.SyncToDisk();System.out.println("\n数据集创建完成！\n");}}得到test.dbf, test.shp, test.shx。

使用PYTHON SHAPEFILE LIBRARY创建和编辑SHAPEFILE文件

使用PYTHON SHAPEFILE LIBRARY创建和编辑SHAPEFILE文件

使用Python Shapefile Library创建和编辑Shapefile文件shapefile是GIS中非常重要的一种数据类型，在ArcGIS中被称为要素类(Feature Classes)，主要包括点(point)、线(polyline)和多边形(polygon)。

Python脚本是ArcGIS官方推荐的脚本语言，通过Python脚本能够很方便的调用ArcGIS中的各种工具和函数批量完成所需操作。

本文所介绍的这个库(Python Shapefile Library)是一个Python库，用于在Python脚本中对ArcGIS中的Shapefile 文件(.shp，.shx，.dbf等格式)进行读写操作。

1、Python Shapefile Library的下载与安装：Python Shapefile Library下载地址：https:///p/pyshp/Python Shapefile Library使用时无需安装，只需在Python程序中导入该模块文件即可(import shapefile，具体导入方法参考Python教程中模块的导入部分)2、Shapefile文件的读操作2.1 Python Shapefile Library提供了Reader类，通过创建Reader类的对象(如下面的sf)进行shapefile文件的读操作：sf = shapefile . Reader ('shapefile name')2.2 使用Python Shapefile Library读取shapefile文件的"几何数据"(Geometry)和"属性数据"(Attribute Record)"几何数据"一般有多个几何对象组成，比如一个"点文件"，每个点就是一个对象；对于一个多边形文件，每个对象可能包含有多个多边形，每个多边形又称为"块(parts)"，每个"块"由多个点组成。

Shapefile文件修改_Python地理空间分析指南（第2版）_[共2页]

Shapefile文件修改_Python地理空间分析指南（第2版）_[共2页]

5.5 Shapefile文件编辑 155也可能包含一个名为m的测量值。

该测量测值可能是和上述形状有关的用户自定义信息。

例如一个给定位置的温度信息，这是另外一种形状类型，它允许为每个形状添加一个m值而非z值。

这种形状数据大家一般叫它M型几何图形。

和z值类似，如果存在该类型数据，就会生成m属性，如果没有那么它们就不是上述类型的数据。

你一般不会遇到包含z值的Shapefile文件，碰到需要设置m值的数据的几率就更小了。

但是如果有机会遇到它们，那么也不必大惊小怪。

只需将dbf文件中的看到的字段和记录一视同仁即可。

如果你不喜欢将z和m值分开存放在List对象中，那么可以使用zip()方法通过点集List对它们进行组合。

zip()方法可以接受多个List对象作为参数，只需要用逗号对它们分隔即可。

5.5.4 Shapefile文件修改当你使用PyShp库创建了一个读取器对象后，它是只读的。

你可以在该读取器对象中修改记录的值，但是它们不会对原Shapefile文件中的记录产生影响。

为了创建一个Shapefile 文件，你也需要创建一个写者对象。

你可以在写入和读取对象中随意修改数据记录的值，因为它们都是Python的动态数据类型。

但是与此同时，你必须将需要修改的值从读取器对象拷贝到写入对象中。

PyShp库能自动处理所有头文件信息，例如边框和记录总数。

你唯一需要关注的是几何图形和属性。

你会发现这个方法比在之前OGR示例中使用的方法容易的多。

不过它的缺点是只适用于UTM投影。

为了演示这一概念，将会读取一个使用度做计量单位的点Shapefile文件，然后在写者对象保存它之前将其参照系统转换为UTM投影。

我们将会用到本章前面提到过的PyShp 和UTM库。

本示例中用到的Shapefile文件是之前重投影为WGS84坐标系统的纽约市博物馆Shapefile文件。

你还可以通过如下地址下载该文件的ZIP压缩格式：http://git.io/vLd8Y请看下列代码：>>> import shapefile>>> import utm>>> r = shapefile.Reader("NYC_MUSEUMS_GEO")>>> w = shapefile.Writer(r.shapeType)>>> w.fields = list(r.fields)>>> w.records.extend(r.records())>>> for s in r.iterShapes():... lon,lat = s.points[0]... y,x,zone,band = utm.from_latlon(lat,lon)... w.point(x,y)>>> w.save("NYC_MUSEUMS_UTM")。