shapefile格式说明及读写代码示例
2008-01-10 11:34
316 查看
Shape files数据说明
Shape files是ESRI提供的一种矢量数据格式,它没有拓扑信息,一个Shape files由一组文件组成,其中必要的基本文件包括坐标文件(.shp)、索引文件(.shx)和属性文件(.dbf)三个文件。
坐标文件的结构说明
坐标文件(.shp)用于记录空间坐标信息。它由头文件和实体信息两部分构成(如图2.1所示)。
坐标文件的文件头
坐标文件的文件头是一个长度固定(100 bytes)的记录段,一共有9个int型和7个double型数据,主要记录内容见表2.2。
[align=center][/align]
……
……
[align=center][/align]
[align=center]图2.1 坐标文件的结构[/align]
[align=center][/align]
[align=center][/align]
[align=center]表2.2shapefiles 头文件表[/align]
[align=center][/align]
注:最后4个加星号特别标示的四个数据只有当这个Shapefile文件包含Z方向坐标或者具有Measure值时才有值,否则为0.0。所谓Measure值,是用于存储需要的附加数据,可以用来记录各种数据,例如权值、道路长度等信息。
位序
细心的读者会注意到表2.2中的数值的位序有Little和big的区别,对于位序是big的数据我们在读取时要小心。通常,数据的位序都是Little,但在有些情况下可能会是big,二者的区别在于它们位序的顺序相反。一个位序为big的数据,如果我们想得到它的真实数值,需要将它的位序转换成Little即可。转换原理非常简单,就是交换字节顺序,下面是作者实现的在两者间进行转换的程序,代码如下:
[align=center][/align]
//位序转换程序
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';
else
sscanf(mass,"%d",&k);
value=value+int(k*pow(16,7-i));
}
return (value);
}
Shapefile文件支持的几何类型(ShapeType)
Shapefile文件所支持的几何类型如表2.3所示:
[align=center][/align]
[align=center][/align]
[align=center]表2.3shapefiles文件支持的几何类型[/align]
[align=center][/align]
对于一个不是记录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:
[align=center][/align]
[align=center][/align]
[align=center]表2.4点状目标的记录内容[/align]
[align=center][/align]
下面是读取点状目标的记录内容的代码:
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中的线状目标是由一系列点坐标串构成,一个线目标可能包括多个子线段,子线段之间可以是相离的,同时子线段之间也可以相交。Shapefile允许出现多个坐标完全相同的连续点,当读取文件时一定要注意这种情况,但是不允许出现某个退化的、长度为0的子线段出现。线状目标的记录内容如表2.5:
[align=center][/align]
[align=center][/align]
[align=center]表2.5线状目标的记录内容[/align]
[align=center][/align]
具体的数据结构如下:
PolyLine
{
Double[4] Box // 当前线状目标的坐标范围
Integer NumParts // 当前线目标所包含的子线段的个数
Integer NumPoints // 当前线目标所包含的顶点个数
Integer[NumParts] Parts // 每个子线段的第一个坐标点在Points的位置
Point[NumPoints] Points // 记录所有坐标点的数组
}
这些记录项的具体含义如下:
Box记录了当前的线目标的坐标范围,它是一个double型的数组,按照Xmin、 Ymin、 Xmax、Ymax的顺序记录了坐标范围;
NumParts记录了当前线目标所包含的子线段的个数;
NumPoints记录了当前线目标的坐标点总数;
Parts记录了每个子线段的第一个坐标点在坐标数组points中的位置,以便读取数据;
Points是用于存放当前线目标的X、Y坐标的数组。
下面是读取线状目标的记录内容的代码:
OnReadLineShp(CString ShpFileName)
{
//打开坐标文件
……
//读取坐标文件头的内容开始
……
//读取线状目标的实体信息
int RecordNumber;
int ContentLength;
int num =0;
while((fread(&RecordNumber, sizeof(int), 1,ShpFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1,ShpFile_fp);
RecordNumber = OnChangeByteOrder (RecordNumber);
ContentLength = OnChangeByteOrder (ContentLength);
int shapeType;
double Box[4];
int NumParts;
int NumPoints;
int *Parts;
fread(&shapeType, sizeof(int), 1,ShpFile_fp);
//读Box
for(i=0;i<4;i++)
fread(Box+i, sizeof(double),1,ShpFile_fp);
//读NumParts和NumPoints
fread(&NumParts, sizeof(int), 1,ShpFile_fp);
fread(&NumPoints, sizeof(int), 1,ShpFile_fp);
//读Parts和Points
Parts =new int[NumParts];
for(i=0;i<NumParts;i++)
fread(Parts+i, sizeof(int), 1,ShpFile_fp);
int pointNum;
for(i=0;i<NumParts;i++)
{
if(i!=NumParts-1)
pointNum =Parts[i+1]-Parts[i];
else
pointNum =NumPoints-Parts[i];
double *PointsX;
double *PointsY;
PointsX =new double[pointNum];
PointsY =new double[pointNum];
for(j=0;j<pointNum;j++)
{
fread(PointsX+j, sizeof(double),1,ShpFile_fp);
fread(PointsY+j, sizeof(double),1,ShpFile_fp);
}
delete[] PointsX;
delete[] PointsY;
}
delete[] Parts;
}
}
面状目标
shapefile中的面状目标是由多个子环构成,每个子环是由至少四个顶点构成的封闭的、无自相交现象的环。对于含有岛的多边形,构成它的环有内外环之分,每个环的顶点的排列顺序或者方向说明了这个环到底是内环还是外环。一个内环的顶点是按照逆时针顺序排列的;而对于外环,它的顶点排列顺序是顺时针方向。如果一个多边形只由一个环构成,那么它的顶点排列顺序肯定是顺时针方向。
每条多边形记录的数据结构与线目标的数据结构完全相同,
Polygon
{
Double[4] Box // 当前面状目标的坐标范围
Integer NumParts // 当前面目标所包含的子环的个数
Integer NumPoints // 构成当前面状目标的所有顶点的个数
Integer[NumParts] Parts // 每个子环的第一个坐标点在Points的位置
Point[NumPoints] Points // 记录所有坐标点的数组
}
对于一个shapefile中的多边形,它必须满足下面三个条件:
构成多边形的每个子环都必须是闭合的,即每个子环的第一个顶点跟最后一个顶点是同一个点;
每个子环在Points数组中的排列顺序并不重要,但每个子环的顶点必须按照一定的顺序连续排列;
存储在shapefile 中的多边形必须是干净的。所谓一个干净的多边形,它必须满足两点:
没有自相交现象。这就要求任何一个子环不能跟其它的子环相交,共线的现 象也将被当作相交。但是允许两个子环的顶点重合;
对于一个不含岛的多边形或者是含岛的多边形的外环,它们的顶点排列顺序必须是顺时针方向;而对于内环,它的排列顺序必须是逆时针方向。所谓的“脏多边形”就是指顶点排列顺序为顺时针的内环。
图2.2中的多边形是一个典型的例子。这个多边形包括一个岛,所有顶点的个数为8。NumParts等于2,NumPoints等于10。请注意内环(岛)的顶点的排列顺序是逆时针的(如图2.3所示)。
[align=center][/align]
[align=center]图2.2 带岛的多边形[/align]
[align=center][/align]
[align=center]图2.3 带岛的多边形的坐标记录[/align]
[align=center][/align]
面状目标的记录内容如表2.6:
[align=center][/align]
[align=center][/align]
[align=center]表2.6面状目标的记录内容[/align]
[align=center][/align]
下面是读取面状目标的记录内容的代码:
void OnReadAreaShp(CString ShpFileName)
{
//打开坐标文件
……
//读取坐标文件头的内容开始
……
//读取面状目标的实体信息
int RecordNumber;
int ContentLength;
while((fread(&RecordNumber, sizeof(int), 1,m_ShpFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1,m_ShpFile_fp);
RecordNumber = OnChangeByteOrder (RecordNumber);
ContentLength = OnChangeByteOrder (ContentLength);
int shapeType;
double Box[4];
int NumParts;
int NumPoints;
int *Parts;
fread(&shapeType, sizeof(int), 1,m_ShpFile_fp);
//读Box
for(i=0;i<4;i++)
fread(Box+i, sizeof(double),1,m_ShpFile_fp);
//读NumParts和NumPoints
fread(&NumParts, sizeof(int), 1,m_ShpFile_fp);
fread(&NumPoints, sizeof(int), 1,m_ShpFile_fp);
//读Parts和Points
Parts =new int[NumParts];
for(i=0;i<NumParts;i++)
fread(Parts+i, sizeof(int), 1,m_ShpFile_fp);
int pointNum;
int xx;
int yy;
for(i=0;i<NumParts;i++)
{
if(i!=NumParts-1)
pointNum =Parts[i+1]-Parts[i];
else
pointNum =NumPoints-Parts[i];
double *PointsX;
double *PointsY;
PointsX =new double[pointNum];
PointsY =new double[pointNum];
for(j=0;j<pointNum;j++)
{
fread(PointsX+j, sizeof(double),1,m_ShpFile_fp);
fread(PointsY+j, sizeof(double),1,m_ShpFile_fp);
}
delete[] PointsX;
delete[] PointsY;
}
delete[] Parts;
}
}
属性文件的结构说明
属性文件(.dbf)用于记录属性信息。它是一个标准的DBF文件,也是由头文件和实体信息两部分构成(如图2.4所示)。
[align=center][/align]
……
……
[align=center][/align]
[align=center]图2.4 属性文件的结构[/align]
[align=center][/align]
属性文件的文件头
其中文件头部分的长度是不定长的,它主要对DBF文件作了一些总体说明(见表2.7),其中最主要的是对这个DBF文件的记录项的信息进行了详细地描述,比如对每个记录项的名称、数据类型、长度等信息都有具体的说明。
[align=center][/align]
[align=center][/align]
[align=center]表2.7属性文件(.dbf)的文件头[/align]
[align=center][/align]
[align=center][/align]
[align=center]表2.8记录项信息描述[/align]
[align=center][/align]
[align=center][/align]
[align=center]表2.9dbf文件中的数据类型[/align]
[align=center][/align]
属性文件的实体信息
实体信息部分就是一条条属性记录,每条记录都是由若干个记录项构成,因此只要依次循环读取每条记录就可以了。
一个读取dbf文件的例子
假设要读取一个名为soil的dbf文件(存储了土地利用信息),它含有8个记录项,记录项信息如表2.10所示:
[align=center][/align]
[align=center][/align]
[align=center]表2.10dbf文件中的数据类型[/align]
[align=center][/align]
下面是读取这个dbf文件的代码:
void OnReadDbf(CString DbfFileName)
{
FILE* m_DbfFile_fp; //****Dbf文件指针
//打开dbf文件
if((m_DbfFile_fp=fopen(DbfFileName,"rb"))==NULL)
{
return;
}
int i,j;
//////****读取dbf文件的文件头 开始
BYTE version;
fread(&version, 1, 1,m_DbfFile_fp);
BYTE date[3];
for(i=0;i<3;i++)
{
fread(date+i, 1, 1,m_DbfFile_fp);
}
int RecordNum; //******
fread(&RecordNum, sizeof(int), 1,m_DbfFile_fp);
short HeaderByteNum;
fread(&HeaderByteNum, sizeof(short), 1,m_DbfFile_fp);
short RecordByteNum
fread(&RecordByteNum, sizeof(short), 1,m_DbfFile_fp);
short Reserved1;
fread(&Reserved1, sizeof(short), 1,m_DbfFile_fp);
BYTE Flag4s;
fread(&Flag4s, sizeof(BYTE), 1,m_DbfFile_fp);
BYTE EncrypteFlag;
fread(&EncrypteFlag, sizeof(BYTE), 1,m_DbfFile_fp);
for(i=0;i<3;i++)
{
fread(&Unused, sizeof(int), 1,m_DbfFile_fp);
}
BYTE MDXFlag;
fread(&MDXFlag, sizeof(BYTE), 1,m_DbfFile_fp);
BYTE LDriID;
fread(&LDriID, sizeof(BYTE), 1,m_DbfFile_fp);
short Reserved2;
fread(&Reserved2, sizeof(short), 1,m_DbfFile_fp);
BYTE name[11];
BYTE fieldType;
int Reserved3;
BYTE fieldLength;
BYTE decimalCount;
short Reserved4;
BYTE workID;
short Reserved5[5];
BYTE mDXFlag1;
int fieldscount;
fieldscount = (HeaderByteNum - 32) / 32;
//读取记录项信息-共有8个记录项
for(i=0;i< HeaderByteNum;i++)
{
//FieldName----11 bytes
fread(name, 11, 1,m_DbfFile_fp);
//FieldType----1 bytes
fread(&fieldType, sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved3----4 bytes
Reserved3 =0;
fread(&Reserved3, sizeof(int), 1,m_DbfFile_fp);
//FieldLength--1 bytes
fread(&fieldLength,sizeof(BYTE), 1,m_DbfFile_fp);
//DecimalCount-1 bytes
fread(&decimalCount,sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved4----2 bytes
Reserved4 =0;
fread(&Reserved4, sizeof(short), 1,m_DbfFile_fp);
//WorkID-------1 bytes
fread(&workID, sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved5----10 bytes
for(j=0;j<5;j++)
{
fread(Reserved5+j,sizeof(short), 1,m_DbfFile_fp);
}
//MDXFlag1-----1 bytes
fread(&mDXFlag1, sizeof(BYTE), 1,m_DbfFile_fp);
}
BYTE terminator;
fread(&terminator, sizeof(BYTE), 1,m_DbfFile_fp);
//读取dbf文件头结束
double Area,Perimeter,Centroid_y,Centroid_x;
int Soils_,Soils_id;
CString Soil_code,suit;
BYTE deleteFlag;
char media[31];
//读取dbf文件记录 开始
for(i=0;i<RecordNum;i++)
{
fread(&deleteFlag, sizeof(BYTE), 1,m_DbfFile_fp);
//读取 Area double
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Area =atof(media);
//读取 Perimeter double
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Perimeter =atof(media);
//读取 soils_ int
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<11;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soils_ =atoi(media);
//读取 Soils_id int
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<11;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soils_id =atoi(media);
//读取 soil_code string
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<3;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soil_code =media;
//读取 suit string
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<1;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
suit =media;
//读取 Centroid_y double
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Centroid_y =atof(media);
//读取 Centroid_x double
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Centroid_x =atof(media);
}
//读取dbf文件记录 结束
}
索引文件的结构说明
索引文件(.shx)主要包含坐标文件的索引信息,文件中每个记录包含对应的坐标文件记录距离坐标文件的文件头的偏移量。通过索引文件可以很方便地在坐标文件中定位到指定目标的坐标信息。
索引文件也是由头文件和实体信息两部分构成(如图2.5),其中文件头部分是一个长度固定(100 bytes)的记录段,其内容与坐标文件的文件头基本一致。它的实体信息以记录为基本单位,每一条记录包括偏移量(offset)和记录段长度(Content Length)两个记录项,它们的位序都是big,两个记录项都是int型,见表2.11。
[align=center][/align]
……
……
[align=center][/align]
[align=center]图2.5 索引文件的结构[/align]
[align=center][/align]
[align=center][/align]
[align=center]表2.11 索引文件的记录内容[/align]
[align=center][/align]
下面是一段读取索引文件的代码:
void OnReadShx(CString ShxFileName)
{
FILE* m_ShxFile_fp; //****Shx文件指针
//打开索引文件
if((m_ShxFile_fp=fopen(ShxFileName,"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_ShxFile_fp);
FileCode = OnChangeByteOrder(FileCode);
for(i=0;i<5;i++)
fread(&Unused,sizeof(int), 1,m_ShxFile_fp);
fread(&FileLength, sizeof(int), 1,m_ShxFile_fp);
FileLength = OnChangeByteOrder(FileLength);
fread(&Version, sizeof(int), 1,m_ShxFile_fp);
fread(&ShapeType, sizeof(int), 1,m_ShxFile_fp);
fread(&Xmin, sizeof(double),1,m_ShxFile_fp);
fread(&Ymin, sizeof(double),1,m_ShxFile_fp);
fread(&Xmax, sizeof(double),1,m_ShxFile_fp);
fread(&Ymax, sizeof(double),1,m_ShxFile_fp);
fread(&Zmin, sizeof(double),1,m_ShxFile_fp);
fread(&Zmax, sizeof(double),1,m_ShxFile_fp);
fread(&Mmin, sizeof(double),1,m_ShxFile_fp);
fread(&Mmax, sizeof(double),1,m_ShxFile_fp);
//读取索引文件头的内容结束
int Offset, ContentLength;
//读取实体信息
while((fread(&Offset, sizeof(int), 1, m_ShxFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1, m_ShxFile_fp);
Offset = OnChangeByteOrder(Offset);
ContentLength = OnChangeByteOrder(ContentLength);
}
}
小结
本节介绍了MapObjects支持的各种数据,并详细介绍了shapefiles的文件结构,同时给出了读取shapefiles的坐标文件(.shp)、属性文件(.dbf)和索引文件(.shx)的程序,给出这些程序的目的在于让读者通过这些例子深入掌握shapefiles文件的格式,进而具备将特定格式的数据文件转换成shapefiles文件的能力。
Shape files是ESRI提供的一种矢量数据格式,它没有拓扑信息,一个Shape files由一组文件组成,其中必要的基本文件包括坐标文件(.shp)、索引文件(.shx)和属性文件(.dbf)三个文件。
坐标文件的结构说明
坐标文件(.shp)用于记录空间坐标信息。它由头文件和实体信息两部分构成(如图2.1所示)。
坐标文件的文件头
坐标文件的文件头是一个长度固定(100 bytes)的记录段,一共有9个int型和7个double型数据,主要记录内容见表2.2。
[align=center][/align]
文件头 | |
记录头 | 记录内容 |
记录头 | 记录内容 |
记录头 | 记录内容 |
记录头 | 记录内容 |
……
记录头 | 记录内容 |
[align=center]图2.1 坐标文件的结构[/align]
[align=center][/align]
[align=center]起始位置[/align] | [align=center]名称[/align] | [align=center]数值[/align] | [align=center]类型[/align] | [align=center]位序[/align] |
[align=center]0[/align] | [align=center]File Code[/align] | [align=center]9994[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]4[/align] | [align=center]Unused[/align] | [align=center]0[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]8[/align] | [align=center]Unused[/align] | [align=center]0[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]12[/align] | [align=center]Unused[/align] | [align=center]0[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]16[/align] | [align=center]Unused[/align] | [align=center]0[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]20[/align] | [align=center]Unused[/align] | [align=center]0[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]24[/align] | [align=center]文件长度[/align] | [align=center]文件的实际长度[/align] | [align=center]Integer[/align] | [align=center]big[/align] |
[align=center]28[/align] | [align=center]版本号[/align] | [align=center]1000[/align] | [align=center]Integer[/align] | [align=center]Little[/align] |
[align=center]32[/align] | [align=center]几何类型[/align] | [align=center]表示这个Shapefile文件所记录的空间数据的几何类型[/align] | [align=center]Integer[/align] | [align=center]Little[/align] |
[align=center]36[/align] | [align=center]Xmin[/align] | [align=center]空间数据所占空间范围的X方向最小值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]44[/align] | [align=center]Ymin[/align] | [align=center]空间数据所占空间范围的Y方向最小值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]52[/align] | [align=center]Xmax[/align] | [align=center]空间数据所占空间范围的X方向最大值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]60[/align] | [align=center]Ymax[/align] | [align=center]空间数据所占空间范围的Y方向最大值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]68*[/align] | [align=center]Zmin[/align] | [align=center]空间数据所占空间范围的Z方向最小值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]76*[/align] | [align=center]Zmax[/align] | [align=center]空间数据所占空间范围的Z方向最大值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]84*[/align] | [align=center]Mmin[/align] | [align=center]最小Measure值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]92*[/align] | [align=center]Mmax[/align] | [align=center]最大Measure值[/align] | [align=center]Double[/align] | [align=center]Little[/align] |
[align=center]表2.2shapefiles 头文件表[/align]
[align=center][/align]
注:最后4个加星号特别标示的四个数据只有当这个Shapefile文件包含Z方向坐标或者具有Measure值时才有值,否则为0.0。所谓Measure值,是用于存储需要的附加数据,可以用来记录各种数据,例如权值、道路长度等信息。
位序
细心的读者会注意到表2.2中的数值的位序有Little和big的区别,对于位序是big的数据我们在读取时要小心。通常,数据的位序都是Little,但在有些情况下可能会是big,二者的区别在于它们位序的顺序相反。一个位序为big的数据,如果我们想得到它的真实数值,需要将它的位序转换成Little即可。转换原理非常简单,就是交换字节顺序,下面是作者实现的在两者间进行转换的程序,代码如下:
[align=center][/align]
//位序转换程序
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';
else
sscanf(mass,"%d",&k);
value=value+int(k*pow(16,7-i));
}
return (value);
}
Shapefile文件支持的几何类型(ShapeType)
Shapefile文件所支持的几何类型如表2.3所示:
[align=center][/align]
[align=center]编号[/align] | [align=center]几何类型[/align] |
[align=center]0[/align] | Null Shape(表示这个Shapefile文件不含坐标) |
[align=center]1[/align] | Point(表示Shapefile文件记录的是点状目标,但不是多点) |
[align=center]3[/align] | PolyLine(表示Shapefile文件记录的是线状目标) |
[align=center]5[/align] | Polygon(表示Shapefile文件记录的是面状目标) |
[align=center]8[/align] | MultiPoint(表示Shapefile文件记录的是多点,即点集合) |
[align=center]11[/align] | PointZ(表示Shapefile文件记录的是三维点状目标) |
[align=center]13[/align] | PolyLineZ(表示Shapefile文件记录的是三维线状目标) |
[align=center]15[/align] | PolygonZ(表示Shapefile文件记录的是三维面状目标) |
[align=center]18[/align] | MultiPointZ(表示Shapefile文件记录的是三维点集合目标) |
[align=center]21[/align] | PointM(表示含有Measure值的点状目标) |
[align=center]23[/align] | PolyLineM(表示含有Measure值的线状目标) |
[align=center]25[/align] | PolygonM(表示含有Measure值的面状目标) |
[align=center]28[/align] | MultiPointM(表示含有Measure值的多点目标) |
[align=center]31[/align] | MultiPatch(表示复合目标) |
[align=center]表2.3shapefiles文件支持的几何类型[/align]
[align=center][/align]
对于一个不是记录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:
[align=center][/align]
[align=center]记录项[/align] | [align=center]数值[/align] | [align=center]数据类型[/align] | [align=center]长度[/align] | [align=center]个数[/align] | [align=center]位序[/align] |
[align=center]几何类型(ShapeType)[/align] | [align=center]1(表示点状目标)[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]X方向坐标[/align] | [align=center]X方向坐标值[/align] | [align=center]double型[/align] | [align=center]8[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]Y方向坐标[/align] | [align=center]Y方向坐标值[/align] | [align=center]double型[/align] | [align=center]8[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]表2.4点状目标的记录内容[/align]
[align=center][/align]
下面是读取点状目标的记录内容的代码:
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中的线状目标是由一系列点坐标串构成,一个线目标可能包括多个子线段,子线段之间可以是相离的,同时子线段之间也可以相交。Shapefile允许出现多个坐标完全相同的连续点,当读取文件时一定要注意这种情况,但是不允许出现某个退化的、长度为0的子线段出现。线状目标的记录内容如表2.5:
[align=center][/align]
[align=center]记录项[/align] | [align=center]数值[/align] | [align=center]数据类型[/align] | [align=center]长度[/align] | [align=center]个数[/align] | [align=center]位序[/align] |
[align=center]几何类型(ShapeType)[/align] | [align=center]3(表示线状目标)[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]坐标范围(Box)[/align] | [align=center]表示当前线目标的坐标范围[/align] | [align=center]double型[/align] | [align=center]32[/align] | [align=center]4[/align] | [align=center]Little[/align] |
[align=center]子线段个数(NumParts)[/align] | [align=center]表示构成当前线目标的子线段的个数[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]坐标点数(NumPoints)[/align] | [align=center]表示构成当前线目标所包含的坐标点个数[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]Parts数组[/align] | [align=center]记录了每个子线段的坐标在Points数组中的起始位置[/align] | [align=center]int型[/align] | [align=center]4×NumParts[/align] | [align=center]NumParts[/align] | [align=center]Little[/align] |
[align=center]Points数组[/align] | [align=center]记录了所有的坐标信息[/align] | [align=center]Point型[/align] | [align=center]根据点个数来确定[/align] | [align=center]NumPoints[/align] | [align=center]Little[/align] |
[align=center]表2.5线状目标的记录内容[/align]
[align=center][/align]
具体的数据结构如下:
PolyLine
{
Double[4] Box // 当前线状目标的坐标范围
Integer NumParts // 当前线目标所包含的子线段的个数
Integer NumPoints // 当前线目标所包含的顶点个数
Integer[NumParts] Parts // 每个子线段的第一个坐标点在Points的位置
Point[NumPoints] Points // 记录所有坐标点的数组
}
这些记录项的具体含义如下:
Box记录了当前的线目标的坐标范围,它是一个double型的数组,按照Xmin、 Ymin、 Xmax、Ymax的顺序记录了坐标范围;
NumParts记录了当前线目标所包含的子线段的个数;
NumPoints记录了当前线目标的坐标点总数;
Parts记录了每个子线段的第一个坐标点在坐标数组points中的位置,以便读取数据;
Points是用于存放当前线目标的X、Y坐标的数组。
下面是读取线状目标的记录内容的代码:
OnReadLineShp(CString ShpFileName)
{
//打开坐标文件
……
//读取坐标文件头的内容开始
……
//读取线状目标的实体信息
int RecordNumber;
int ContentLength;
int num =0;
while((fread(&RecordNumber, sizeof(int), 1,ShpFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1,ShpFile_fp);
RecordNumber = OnChangeByteOrder (RecordNumber);
ContentLength = OnChangeByteOrder (ContentLength);
int shapeType;
double Box[4];
int NumParts;
int NumPoints;
int *Parts;
fread(&shapeType, sizeof(int), 1,ShpFile_fp);
//读Box
for(i=0;i<4;i++)
fread(Box+i, sizeof(double),1,ShpFile_fp);
//读NumParts和NumPoints
fread(&NumParts, sizeof(int), 1,ShpFile_fp);
fread(&NumPoints, sizeof(int), 1,ShpFile_fp);
//读Parts和Points
Parts =new int[NumParts];
for(i=0;i<NumParts;i++)
fread(Parts+i, sizeof(int), 1,ShpFile_fp);
int pointNum;
for(i=0;i<NumParts;i++)
{
if(i!=NumParts-1)
pointNum =Parts[i+1]-Parts[i];
else
pointNum =NumPoints-Parts[i];
double *PointsX;
double *PointsY;
PointsX =new double[pointNum];
PointsY =new double[pointNum];
for(j=0;j<pointNum;j++)
{
fread(PointsX+j, sizeof(double),1,ShpFile_fp);
fread(PointsY+j, sizeof(double),1,ShpFile_fp);
}
delete[] PointsX;
delete[] PointsY;
}
delete[] Parts;
}
}
面状目标
shapefile中的面状目标是由多个子环构成,每个子环是由至少四个顶点构成的封闭的、无自相交现象的环。对于含有岛的多边形,构成它的环有内外环之分,每个环的顶点的排列顺序或者方向说明了这个环到底是内环还是外环。一个内环的顶点是按照逆时针顺序排列的;而对于外环,它的顶点排列顺序是顺时针方向。如果一个多边形只由一个环构成,那么它的顶点排列顺序肯定是顺时针方向。
每条多边形记录的数据结构与线目标的数据结构完全相同,
Polygon
{
Double[4] Box // 当前面状目标的坐标范围
Integer NumParts // 当前面目标所包含的子环的个数
Integer NumPoints // 构成当前面状目标的所有顶点的个数
Integer[NumParts] Parts // 每个子环的第一个坐标点在Points的位置
Point[NumPoints] Points // 记录所有坐标点的数组
}
对于一个shapefile中的多边形,它必须满足下面三个条件:
构成多边形的每个子环都必须是闭合的,即每个子环的第一个顶点跟最后一个顶点是同一个点;
每个子环在Points数组中的排列顺序并不重要,但每个子环的顶点必须按照一定的顺序连续排列;
存储在shapefile 中的多边形必须是干净的。所谓一个干净的多边形,它必须满足两点:
没有自相交现象。这就要求任何一个子环不能跟其它的子环相交,共线的现 象也将被当作相交。但是允许两个子环的顶点重合;
对于一个不含岛的多边形或者是含岛的多边形的外环,它们的顶点排列顺序必须是顺时针方向;而对于内环,它的排列顺序必须是逆时针方向。所谓的“脏多边形”就是指顶点排列顺序为顺时针的内环。
图2.2中的多边形是一个典型的例子。这个多边形包括一个岛,所有顶点的个数为8。NumParts等于2,NumPoints等于10。请注意内环(岛)的顶点的排列顺序是逆时针的(如图2.3所示)。
[align=center][/align]
[align=center]图2.2 带岛的多边形[/align]
[align=center][/align]
[align=center]0[/align] |
[align=center]0[/align] |
[align=center]5[/align] |
[align=center]5[/align] |
[align=center]v1[/align] |
[align=center]0[/align] |
[align=center]v2[/align] |
[align=center]1[/align] |
[align=center]v3[/align] |
[align=center]2[/align] |
[align=center]v4[/align] |
[align=center]3[/align] |
[align=center]v1[/align] |
[align=center]4[/align] |
[align=center]v5[/align] |
[align=center]5[/align] |
[align=center]v8[/align] |
[align=center]6[/align] |
[align=center]v7[/align] |
[align=center]7[/align] |
[align=center]v6[/align] |
[align=center]8[/align] |
[align=center]v5[/align] |
[align=center]9[/align] |
[align=center][/align]
面状目标的记录内容如表2.6:
[align=center][/align]
[align=center]记录项[/align] | [align=center]数值[/align] | [align=center]数据类型[/align] | [align=center]长度[/align] | [align=center]个数[/align] | [align=center]位序[/align] |
[align=center]几何类型(ShapeType)[/align] | [align=center]5(表示面状目标)[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]坐标范围(Box)[/align] | [align=center]表示当前面目标的坐标范围[/align] | [align=center]double型[/align] | [align=center]32[/align] | [align=center]4[/align] | [align=center]Little[/align] |
[align=center]子线段个数(NumParts)[/align] | [align=center]表示构成当前面状目标的子环的个数[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]坐标点数(NumPoints)[/align] | [align=center]表示构成当前面状目标所包含的坐标点个数[/align] | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Little[/align] |
[align=center]Parts数组[/align] | [align=center]记录了每个子环的坐标在Points数组中的起始位置[/align] | [align=center]int型[/align] | [align=center]4×NumParts[/align] | [align=center]NumParts[/align] | [align=center]Little[/align] |
[align=center]Points数组[/align] | [align=center]记录了所有的坐标信息[/align] | [align=center]Point型[/align] | [align=center]根据点个数来确定[/align] | [align=center]NumPoints[/align] | [align=center]Little[/align] |
[align=center]表2.6面状目标的记录内容[/align]
[align=center][/align]
下面是读取面状目标的记录内容的代码:
void OnReadAreaShp(CString ShpFileName)
{
//打开坐标文件
……
//读取坐标文件头的内容开始
……
//读取面状目标的实体信息
int RecordNumber;
int ContentLength;
while((fread(&RecordNumber, sizeof(int), 1,m_ShpFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1,m_ShpFile_fp);
RecordNumber = OnChangeByteOrder (RecordNumber);
ContentLength = OnChangeByteOrder (ContentLength);
int shapeType;
double Box[4];
int NumParts;
int NumPoints;
int *Parts;
fread(&shapeType, sizeof(int), 1,m_ShpFile_fp);
//读Box
for(i=0;i<4;i++)
fread(Box+i, sizeof(double),1,m_ShpFile_fp);
//读NumParts和NumPoints
fread(&NumParts, sizeof(int), 1,m_ShpFile_fp);
fread(&NumPoints, sizeof(int), 1,m_ShpFile_fp);
//读Parts和Points
Parts =new int[NumParts];
for(i=0;i<NumParts;i++)
fread(Parts+i, sizeof(int), 1,m_ShpFile_fp);
int pointNum;
int xx;
int yy;
for(i=0;i<NumParts;i++)
{
if(i!=NumParts-1)
pointNum =Parts[i+1]-Parts[i];
else
pointNum =NumPoints-Parts[i];
double *PointsX;
double *PointsY;
PointsX =new double[pointNum];
PointsY =new double[pointNum];
for(j=0;j<pointNum;j++)
{
fread(PointsX+j, sizeof(double),1,m_ShpFile_fp);
fread(PointsY+j, sizeof(double),1,m_ShpFile_fp);
}
delete[] PointsX;
delete[] PointsY;
}
delete[] Parts;
}
}
属性文件的结构说明
属性文件(.dbf)用于记录属性信息。它是一个标准的DBF文件,也是由头文件和实体信息两部分构成(如图2.4所示)。
[align=center][/align]
文件头 | |
记录1 | |
记录2 | |
记录3 | |
记录4 |
……
记录n |
[align=center]图2.4 属性文件的结构[/align]
[align=center][/align]
属性文件的文件头
其中文件头部分的长度是不定长的,它主要对DBF文件作了一些总体说明(见表2.7),其中最主要的是对这个DBF文件的记录项的信息进行了详细地描述,比如对每个记录项的名称、数据类型、长度等信息都有具体的说明。
[align=center][/align]
[align=center]在文件中的位置[/align] | [align=center]内容[/align] | [align=center]说明[/align] |
[align=center]0[/align] | [align=center]1个字节[/align] | 表示当前的版本信息 |
[align=center]1-3[/align] | [align=center]3个字节[/align] | 表示最近的更新日期,按照YYMMDD格式。 |
[align=center]4-7[/align] | [align=center]1个32位数[/align] | 文件中的记录条数。 |
[align=center]8-9[/align] | [align=center]1个16位数[/align] | 文件头中的字节数。 |
[align=center]10-11[/align] | [align=center]1个16位数[/align] | 一条记录中的字节长度。 |
[align=center]12-13[/align] | [align=center]2个字节[/align] | 保留字节,用于以后添加新的说明性信息时使用,这里用0来填写。 |
[align=center]14[/align] | [align=center]1个字节[/align] | 表示未完成的操作。 |
[align=center]15[/align] | [align=center]1个字节[/align] | dBASE IV编密码标记。 |
[align=center]16-27[/align] | [align=center]12个字节[/align] | 保留字节,用于多用户处理时使用。 |
[align=center]28[/align] | [align=center]1个字节[/align] | DBF文件的MDX标识。在创建一个DBF 表时,如果使用了MDX 格式的索引文件,那么 DBF 表的表头中的这个字节就自动被设置了一个标志,当你下次试图重新打开这个DBF表的时候,数据引擎会自动识别这个标志,如果此标志为真,则数据引擎将试图打开相应的MDX 文件。 |
[align=center]29[/align] | [align=center]1个字节[/align] | Language driver ID. |
[align=center]30-31[/align] | [align=center]2个字节[/align] | 保留字节,用于以后添加新的说明性信息时使用,这里用0来填写。 |
[align=center]32-X[/align] | [align=center](n*32)个字节[/align] | 记录项信息描述数组。n表示记录项的个数。这个数组的结构在表2.8中有详细的解释。 |
[align=center]X+1[/align] | [align=center]1个字节[/align] | 作为记录项终止标识。 |
[align=center]表2.7属性文件(.dbf)的文件头[/align]
[align=center][/align]
[align=center]位置[/align] | [align=center]内容[/align] | [align=center]说明[/align] |
[align=center]0-10[/align] | [align=center]11个字节[/align] | 记录项名称,是ASCII码值。 |
[align=center]11[/align] | [align=center]1个字节[/align] | 记录项的数据类型,是ASCII码值。(B、C、D、G、L、M和N,具体的解释见表2.9)。 |
[align=center]12-15[/align] | [align=center]4个字节[/align] | 保留字节,用于以后添加新的说明性信息时使用,这里用0来填写。 |
[align=center]16[/align] | [align=center]1个字节[/align] | 记录项长度,二进制型。 |
[align=center]17[/align] | [align=center]1个字节[/align] | 记录项的精度,二进制型。 |
[align=center]18-19[/align] | [align=center]2个字节[/align] | 保留字节,用于以后添加新的说明性信息时使用,这里用0来填写。 |
[align=center]20[/align] | [align=center]1个字节[/align] | 工作区ID。 |
[align=center]21-30[/align] | [align=center]10个字节[/align] | 保留字节,用于以后添加新的说明性信息时使用,这里用0来填写。 |
[align=center]31[/align] | [align=center]1个字节[/align] | MDX标识。如果存在一个MDX 格式的索引文件,那么这个记录项为真,否则为空。 |
[align=center]表2.8记录项信息描述[/align]
[align=center][/align]
[align=center]代码[/align] | [align=center]数据类型[/align] | [align=center]允许输入的数据[/align] |
[align=center]B[/align] | [align=center]二进制型[/align] | 各种字符。 |
[align=center]C[/align] | [align=center]字符型[/align] | 各种字符。 |
[align=center]D[/align] | [align=center]日期型[/align] | 用于区分年、月、日的数字和一个字符,内部存储按照YYYYMMDD格式。 |
[align=center]G[/align] | [align=center](General [/align] [align=center]or OLE)[/align] | 各种字符。 |
[align=center]N[/align] | [align=center]数值型(Numeric)[/align] | - . 0 1 2 3 4 5 6 7 8 9 |
[align=center]L[/align] | [align=center]逻辑型(Logical)[/align] | ? Y y N n T t F f (? 表示没有初始化)。 |
[align=center]M[/align] | [align=center](Memo)[/align] | 各种字符。 |
[align=center]表2.9dbf文件中的数据类型[/align]
[align=center][/align]
属性文件的实体信息
实体信息部分就是一条条属性记录,每条记录都是由若干个记录项构成,因此只要依次循环读取每条记录就可以了。
一个读取dbf文件的例子
假设要读取一个名为soil的dbf文件(存储了土地利用信息),它含有8个记录项,记录项信息如表2.10所示:
[align=center][/align]
[align=center]记录项名称[/align] | [align=center]数据类型[/align] | [align=center]长度[/align] | [align=center]小数位数[/align] |
[align=center]Area [/align] | [align=center]数值型(double)[/align] | [align=center]31[/align] | [align=center]15[/align] |
[align=center]Perimeter[/align] | [align=center]数值型(double)[/align] | [align=center]31[/align] | [align=center]15[/align] |
[align=center]soils_[/align] | [align=center]数值型(int)[/align] | [align=center]11[/align] | [align=center]0[/align] |
[align=center]soils_id[/align] | [align=center]数值型(int)[/align] | [align=center]11[/align] | [align=center]0[/align] |
[align=center]soil_code[/align] | [align=center]字符型(character)[/align] | [align=center]3[/align] | [align=center]无[/align] |
[align=center]Suit[/align] | [align=center]字符型(character)[/align] | [align=center]1[/align] | [align=center]无[/align] |
[align=center]Centroid_x[/align] | [align=center]数值型(double)[/align] | [align=center]31[/align] | [align=center]15[/align] |
[align=center]Centroid_y[/align] | [align=center]数值型(double)[/align] | [align=center]31[/align] | [align=center]15[/align] |
[align=center]表2.10dbf文件中的数据类型[/align]
[align=center][/align]
下面是读取这个dbf文件的代码:
void OnReadDbf(CString DbfFileName)
{
FILE* m_DbfFile_fp; //****Dbf文件指针
//打开dbf文件
if((m_DbfFile_fp=fopen(DbfFileName,"rb"))==NULL)
{
return;
}
int i,j;
//////****读取dbf文件的文件头 开始
BYTE version;
fread(&version, 1, 1,m_DbfFile_fp);
BYTE date[3];
for(i=0;i<3;i++)
{
fread(date+i, 1, 1,m_DbfFile_fp);
}
int RecordNum; //******
fread(&RecordNum, sizeof(int), 1,m_DbfFile_fp);
short HeaderByteNum;
fread(&HeaderByteNum, sizeof(short), 1,m_DbfFile_fp);
short RecordByteNum
fread(&RecordByteNum, sizeof(short), 1,m_DbfFile_fp);
short Reserved1;
fread(&Reserved1, sizeof(short), 1,m_DbfFile_fp);
BYTE Flag4s;
fread(&Flag4s, sizeof(BYTE), 1,m_DbfFile_fp);
BYTE EncrypteFlag;
fread(&EncrypteFlag, sizeof(BYTE), 1,m_DbfFile_fp);
for(i=0;i<3;i++)
{
fread(&Unused, sizeof(int), 1,m_DbfFile_fp);
}
BYTE MDXFlag;
fread(&MDXFlag, sizeof(BYTE), 1,m_DbfFile_fp);
BYTE LDriID;
fread(&LDriID, sizeof(BYTE), 1,m_DbfFile_fp);
short Reserved2;
fread(&Reserved2, sizeof(short), 1,m_DbfFile_fp);
BYTE name[11];
BYTE fieldType;
int Reserved3;
BYTE fieldLength;
BYTE decimalCount;
short Reserved4;
BYTE workID;
short Reserved5[5];
BYTE mDXFlag1;
int fieldscount;
fieldscount = (HeaderByteNum - 32) / 32;
//读取记录项信息-共有8个记录项
for(i=0;i< HeaderByteNum;i++)
{
//FieldName----11 bytes
fread(name, 11, 1,m_DbfFile_fp);
//FieldType----1 bytes
fread(&fieldType, sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved3----4 bytes
Reserved3 =0;
fread(&Reserved3, sizeof(int), 1,m_DbfFile_fp);
//FieldLength--1 bytes
fread(&fieldLength,sizeof(BYTE), 1,m_DbfFile_fp);
//DecimalCount-1 bytes
fread(&decimalCount,sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved4----2 bytes
Reserved4 =0;
fread(&Reserved4, sizeof(short), 1,m_DbfFile_fp);
//WorkID-------1 bytes
fread(&workID, sizeof(BYTE), 1,m_DbfFile_fp);
//Reserved5----10 bytes
for(j=0;j<5;j++)
{
fread(Reserved5+j,sizeof(short), 1,m_DbfFile_fp);
}
//MDXFlag1-----1 bytes
fread(&mDXFlag1, sizeof(BYTE), 1,m_DbfFile_fp);
}
BYTE terminator;
fread(&terminator, sizeof(BYTE), 1,m_DbfFile_fp);
//读取dbf文件头结束
double Area,Perimeter,Centroid_y,Centroid_x;
int Soils_,Soils_id;
CString Soil_code,suit;
BYTE deleteFlag;
char media[31];
//读取dbf文件记录 开始
for(i=0;i<RecordNum;i++)
{
fread(&deleteFlag, sizeof(BYTE), 1,m_DbfFile_fp);
//读取 Area double
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Area =atof(media);
//读取 Perimeter double
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Perimeter =atof(media);
//读取 soils_ int
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<11;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soils_ =atoi(media);
//读取 Soils_id int
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<11;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soils_id =atoi(media);
//读取 soil_code string
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<3;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Soil_code =media;
//读取 suit string
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<1;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
suit =media;
//读取 Centroid_y double
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Centroid_y =atof(media);
//读取 Centroid_x double
for(j=0;j<31;j++)
strcpy(media+j,"");
for(j=0;j<31;j++)
fread(media+j, sizeof(char), 1,m_DbfFile_fp);
Centroid_x =atof(media);
}
//读取dbf文件记录 结束
}
索引文件的结构说明
索引文件(.shx)主要包含坐标文件的索引信息,文件中每个记录包含对应的坐标文件记录距离坐标文件的文件头的偏移量。通过索引文件可以很方便地在坐标文件中定位到指定目标的坐标信息。
索引文件也是由头文件和实体信息两部分构成(如图2.5),其中文件头部分是一个长度固定(100 bytes)的记录段,其内容与坐标文件的文件头基本一致。它的实体信息以记录为基本单位,每一条记录包括偏移量(offset)和记录段长度(Content Length)两个记录项,它们的位序都是big,两个记录项都是int型,见表2.11。
[align=center][/align]
文件头 | |
记录1 | |
记录2 | |
记录3 | |
记录4 |
……
记录n |
[align=center]图2.5 索引文件的结构[/align]
[align=center][/align]
[align=center]记录项[/align] | [align=center]数值[/align] | [align=center]数据类型[/align] | [align=center]长度[/align] | [align=center]个数[/align] | [align=center]位序[/align] |
[align=center]位移量(Offset)[/align] | 表示坐标文件中的对应记录的起始位置相对于坐标文件起始位置的位移量。 | [align=center][/align] [align=center]int型[/align] | [align=center][/align] [align=center]4[/align] | [align=center][/align] [align=center]1[/align] | [align=center][/align] [align=center]Big[/align] |
[align=center]记录长度[/align] [align=center](Content Length)[/align] | 表示坐标文件中的对应记录的长度。 | [align=center]int型[/align] | [align=center]4[/align] | [align=center]1[/align] | [align=center]Big[/align] |
[align=center]表2.11 索引文件的记录内容[/align]
[align=center][/align]
下面是一段读取索引文件的代码:
void OnReadShx(CString ShxFileName)
{
FILE* m_ShxFile_fp; //****Shx文件指针
//打开索引文件
if((m_ShxFile_fp=fopen(ShxFileName,"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_ShxFile_fp);
FileCode = OnChangeByteOrder(FileCode);
for(i=0;i<5;i++)
fread(&Unused,sizeof(int), 1,m_ShxFile_fp);
fread(&FileLength, sizeof(int), 1,m_ShxFile_fp);
FileLength = OnChangeByteOrder(FileLength);
fread(&Version, sizeof(int), 1,m_ShxFile_fp);
fread(&ShapeType, sizeof(int), 1,m_ShxFile_fp);
fread(&Xmin, sizeof(double),1,m_ShxFile_fp);
fread(&Ymin, sizeof(double),1,m_ShxFile_fp);
fread(&Xmax, sizeof(double),1,m_ShxFile_fp);
fread(&Ymax, sizeof(double),1,m_ShxFile_fp);
fread(&Zmin, sizeof(double),1,m_ShxFile_fp);
fread(&Zmax, sizeof(double),1,m_ShxFile_fp);
fread(&Mmin, sizeof(double),1,m_ShxFile_fp);
fread(&Mmax, sizeof(double),1,m_ShxFile_fp);
//读取索引文件头的内容结束
int Offset, ContentLength;
//读取实体信息
while((fread(&Offset, sizeof(int), 1, m_ShxFile_fp)!=0))
{
fread(&ContentLength,sizeof(int), 1, m_ShxFile_fp);
Offset = OnChangeByteOrder(Offset);
ContentLength = OnChangeByteOrder(ContentLength);
}
}
小结
本节介绍了MapObjects支持的各种数据,并详细介绍了shapefiles的文件结构,同时给出了读取shapefiles的坐标文件(.shp)、属性文件(.dbf)和索引文件(.shx)的程序,给出这些程序的目的在于让读者通过这些例子深入掌握shapefiles文件的格式,进而具备将特定格式的数据文件转换成shapefiles文件的能力。
相关文章推荐
- shapefile格式说明及读写代码示例(C++)(转)
- shapefile格式的说明及读写代码示例
- shapefile格式说明及读写代码示例(C++)(转)
- shapefile格式说明及读写代码示例
- shapefile格式说明及读写代码示例(VB)(转)
- shapefile格式说明及读写代码示例
- Shapefile格式说明及读写代码示例
- 读写ESRI Shapefile格式的C#代码
- java 读写Parquet格式的数据的示例代码
- shapefile格式说明及读写代码示例
- shapefile格式说明及读写代码示例 --http://www.gispower.org/article/arcgis/arcother/2008/48/0848115049GB922C13K2I22H7G06A4.html
- jquery 验证Email邮箱格式的代码示例
- ASP.Net 2.0 窗体身份验证机制-转+自己代码注释示例与更详细的说明
- Python文件读写保存操作的示例代码
- io流操作之读写示例代码(一)
- json格式的Ajax提交示例代码
- json格式的Ajax提交示例代码
- html table表数据转Json格式示例代码
- Lua读写文件代码示例