mapgis6.7二次开发vc6_demo 之四(VCT导出之 version=1,topo=0)
本文所述内容为国土行业的VCT版本为1,拓扑=0的结构
即为2002年国家标准
本文主要包括有(
vct文件头段输出功能;
vct要素类段输出功能;
vct字段结构输出功能;
vct线要素几何图形输出功能;
vct属性结构输出功能;)
其中vct点要素几何图形输出功能;
vct面要素几何图形输出功能;
vct注记几何图形输出功能;
将在后续的文章中描述;
//Vct_Export.cpp文件内容如下所示：
// // //VCT 导出功能 version=1, topo=0 // // #include "stdafx.h" #include "resource.h" #include <math.h> #include "win_area.h" #include "win_oper.h" #include "map_out.h" #include "demo.h" // extern char* AppTitle; //应用程序名称 extern AREA_HINST AHInst; extern HWND MainWnd; //主窗口句柄 extern HWND GrpWnd; //图形窗口句柄 //long _GetDiskSpace(unsigned char dsk) //参数：入口 dsk 磁盘号, 1=A,2=B... // 出口 无 //返回值 返回磁盘剩余空间的字节数 char* GetDiskSpace_lc(unsigned char dsk) { long sy=_GetDiskSpace(dsk); // long M=sy/1024; char* atitle; int mLen=sizeof(M); atitle=new char[mLen]; ltoa(M,atitle,mLen); return atitle; } char* l2cp(long n) { int nLen=sizeof(n); char* atitle=new char[nLen]; ltoa(n,atitle,nLen); return atitle; } char* i2cp(int n) { int nLen=sizeof(n); char* atitle=new char[nLen]; itoa(n,atitle,nLen); return atitle; } char* d2cp(double value) { /* char *tt; int dec, sign; int ndig = 10; tt = ecvt(value, ndig, &dec, &sign); //tt = fcvt(value, ndig, &dec, &sign); return tt;*/ char* s=new char[sizeof(value)]; sprintf(s,"%8.6f",value); char* ret_s=new char[strlen(s)]; strcpy(ret_s,s); return ret_s; } char* short2cp(short value) { char* s=new char[sizeof(value)]; sprintf(s,"%d",value); char* ret_s=new char[strlen(s)]; strcpy(ret_s,s); return ret_s; } char* a2u(const char* A,const char* B) { return strcat_lc(A,B); } char* a3u(const char* A,const char* B,const char* C) { char* t=strcat_lc(A,B); t=strcat_lc(t,C); return t; } char* a4u(const char* A,const char* B,const char* C,const char* D) { char* t=strcat_lc(A,B); t=strcat_lc(t,C); t=strcat_lc(t,D); return t; } char* a5u(const char* A,const char* B,const char* C,const char* D,const char* E) { char* t=strcat_lc(A,B); t=strcat_lc(t,C); t=strcat_lc(t,D); t=strcat_lc(t,E); return t; } char* a6u(const char* A,const char* B,const char* C,const char* D,const char* E,const char* F) { char* t=strcat_lc(A,B); t=strcat_lc(t,C); t=strcat_lc(t,D); t=strcat_lc(t,E); t=strcat_lc(t,F); return t; } char* a7u(const char* a1,const char* a2,const char* a3,const char* a4,const char* a5,const char* a6,const char* a7) { char* t=strcat_lc(a1,a2); t=strcat_lc(t,a3); t=strcat_lc(t,a4); t=strcat_lc(t,a5); t=strcat_lc(t,a6); t=strcat_lc(t,a7); return t; } char* substr_lc(const char* cs,char* fd) { return strstr(cs,fd); } char* getdir_lc(char* s,const char* d) { char* s_s=new char[strlen(s)]; strcpy(s_s,s); char* s_rev=strrev(s_s); char* t_rev=strstr(s_rev,d); char* t=strrev(t_rev); return t; } char* getfilename_lc(char* s,const char* d) { char* t_dir=getdir_lc(s,d); // char* s_s=new char[strlen(s)]; strcpy(s_s,s); int dir_len=strlen(t_dir); int fn_len=strlen(s_s)-dir_len; char* fn_ext=new char[fn_len]; memset(fn_ext,0,fn_len); // for(int i=0;i<fn_len;i++) { fn_ext[i]=s_s[dir_len+i]; } char* ret_s=new char[strlen(fn_ext)]; strcpy(ret_s,fn_ext); return ret_s; } void ShowMsg(char* msg) { ::MessageBox(NULL,msg,"提示",0); //_MapGisAbout(msg); } void TestGetDiskSpace() { char* atitle=GetDiskSpace_lc(4); //char* atitle=AppTitle; ShowMsg(atitle); //MessageBox(hWnd,(LPCSTR)m,(LPCSTR)m,NULL); //_DispAbout(); } void D2A_test() { //查看属性 _ViewAtt(MainWnd,AHInst,0); //查看mapgis图层字段结构api _ViewStruct(MainWnd,AHInst); char* s; double value; int dec,sign; int ndig=2; //字符串长度 value=35951234.876533324234; s=ecvt(value,ndig,&dec,&sign); ShowMsg(a5u(s,",dec=",i2cp(dec),",sign=",i2cp(sign))); sprintf(s,"%8.6f",value); //OK ShowMsg(s); // value=-123.45; ndig=15; s=ecvt(value,ndig,&dec,&sign); ShowMsg(a5u(s,",dec=",i2cp(dec),",sign=",i2cp(sign))); // value=0.6789e5; ndig=5; s=ecvt(value,ndig,&dec,&sign); ShowMsg(a5u(s,",dec=",i2cp(dec),",sign=",i2cp(sign))); return; } //输出VCT文件头 OK bool OutVctHeadBeginEnd() { bool rbc = false; //ShowMsg(this, "正在输出Vct文件头信息"); //输出VCT头标志 double L0 = 105; double tmp_a = 6378140; //默认为西安80椭球体参数 double tmp_b = 6356755.288157528574350308626453; //6356755; char* Spheroid = "Krassovesky Spheroid"; char* tmpDate ="20120412";// System.DateTime.Now.ToString("yyyyMMdd"); char* Parameters=a3u(d2cp(tmp_a), ",",d2cp(tmp_b)); // f_WriteLine("HeadBegin"); f_WriteLine("DataMark: CNSDTF-VCT"); //(国家标准 地球空间数据交换格式2002年) f_WriteLine("Version: 1.0"); f_WriteLine("Unit: M"); f_WriteLine("Dim: 2"); f_WriteLine("Topo: 0"); f_WriteLine("Coordinate: M"); f_WriteLine("Projection: Gauss-Kruger Projection"); f_WriteLine(a2u("Spheroid: ",Spheroid)); //Spheroid: Krassovesky Spheroid f_WriteLine(a2u("Parameters: ",Parameters)); //Parameters: 6378140,6356755 f_WriteLine(a2u("Meridinan: ", d2cp(L0))); //Meridinan: 105 f_WriteLine("MinX: 0.000000"); f_WriteLine("MinY: 0.000000"); f_WriteLine("MaxX: 99999999.999999"); f_WriteLine("MaxY: 9999999.999999"); f_WriteLine("ScaleM:10000"); f_WriteLine(a2u("Date:" , tmpDate)); f_WriteLine("Separator: ,"); f_WriteLine("HeadEnd"); f_WriteLine(""); //ShowMsg(this, "完成输出Vct文件头信息"); rbc = true; return rbc; } //输出vct要素代码 OK bool OutVctFeatureCodeBeginEnd(short* ptALst) { bool rbc = false; char* GeoVctType = ""; char* FeatClassName = ""; char* FeatureCode = ""; char* FeatureName = ""; short ai,LayerType; char* Lfilename; //this.ShowMsg(this, "正在输出要素代码信息"); //#region 输出要素代码 int indexLayer = 0; f_WriteLine("FeatureCodeBegin"); for(int k=1;k<=ptALst[0];k++) { //获取一个工作区 ai=ptALst[k]; //获取工作区文件名 Lfilename=_GetAreaFileName(ai); //获取图层类型 LayerType=_GetAreaType(ptALst[k]); GeoVctType="UnKown"; switch(LayerType) { case REG: GeoVctType="Polygon"; break; case LIN: GeoVctType="Line"; break; case PNT: GeoVctType="Point"; break; default: break; } FeatClassName =Lfilename; FeatureCode =a3u(GeoVctType,"_",i2cp(k)); FeatureName =FeatureCode; //this.ShowMsg(this, "正在输出要素代码" + zhfeatclass.GetAliasName + "[" + FeatClassName + "]"); //this.SetValue(++indexLayer); f_WriteLine(a7u(FeatureCode , "," , FeatureName , "," , GeoVctType , ",0," , FeatClassName)); //Polygon/Line/Point //f_WriteLine("B11,地类图斑,Polygon,0,DLTB"); //f_WriteLine("B12,线状地物,Line,0,XZDW"); //f_WriteLine("C10,土地用途分区,Polygon,0,YTFQ"); //f_WriteLine("C20,基本农田保护,Polygon,0,JBNT"); //f_WriteLine("C30,土地整理,Polygon,0,TDZL"); //f_WriteLine("C40,土地复垦,Polygon,0,TDFK"); //f_WriteLine("C50,土地开发,Polygon,0,TDKF"); //f_WriteLine("C71,面状建设项目,Polygon,0,MZXM"); //f_WriteLine("C72,线状建设项目,Line,0,XZXM"); //f_WriteLine("E11,面状用地,Polygon,0,MZYD"); } f_WriteLine("FeatureCodeEnd"); f_WriteLine(""); //this.ShowMsg(this, "完成输出要素代码信息"); rbc = true; return rbc; } //输出vct表结构 OK /* */ bool OutVctTableStructureBeginEnd(short* ptALst) { bool rbc = false; char* GeoVctType; char* FeatClassName; char* FeatureCode; char* FeatureName; short ai,LayerType; char* Lfilename; //this.ShowMsg(this, "正在输出表结构信息"); //输出表结构 f_WriteLine("TableStructureBegin"); for(int k=1;k<=ptALst[0];k++) { //获取一个工作区 ai=ptALst[k]; //获取工作区文件名 Lfilename=_GetAreaFileName(ai); //获取图层类型 LayerType=_GetAreaType(ptALst[k]); GeoVctType="UnKown"; switch(LayerType) { case REG: GeoVctType="Polygon"; break; case LIN: GeoVctType="Line"; break; case PNT: GeoVctType="Point"; break; default: break; } FeatClassName =Lfilename; FeatureCode =a3u(GeoVctType,"_",i2cp(k)); FeatureName =FeatureCode; //this.ShowMsg(this, "正在输出表结构" + zhfeatclass.GetAliasName + "[" + FeatClassName + "]"); //获取字段集合 CATT_STRU *stru; _GetAttStru(ai,LayerType,&stru); //弹出 [线属性查询]窗体 /* char* expStr=new char[256]; if(_InputExpresion(MainWnd,stru,expStr,256,"线属性查询",NULL)) { ShowMsg(expStr); }*/ // CINFO_HEAD cinfo=stru->hd; short numbfield=cinfo.numbfield; //字段个数 f_WriteLine(a3u(FeatClassName , "," , short2cp(numbfield))); //DLTB,25 // char* FieldTypeStr = "Char"; short FieldLength = 50; short FieldJd = -1; CFIELD_HEAD *pfd=stru[0].fldEntry; for(int k=0;k<numbfield;k++) { // char* fdName =pfd[k].fieldname; //字段名称 char fdtype =pfd[k].fieldtype; //字段类型 short fdlength =pfd[k].msk_leng; //字段长度 short fd_digits=pfd[k].point_leng; //字段小数位数 switch (fdtype) { case TEXT_TYPE: case STR_TYPE: FieldTypeStr = "Char"; FieldLength = fdlength; FieldJd = -1; break; case LONG_TYPE: case SHORT_TYPE: FieldTypeStr = "Integer"; FieldLength = fdlength; FieldJd = -1; break; case FLOAT_TYPE: case DOUBLE_TYPE: FieldTypeStr = "Float"; FieldLength = fdlength; FieldJd = fd_digits; break; case DATE_TYPE: FieldTypeStr = "Date"; FieldLength = fdlength; FieldJd = -1; break; case BIN_DATA_TYPE: FieldTypeStr = "Varbin"; FieldLength = fdlength; FieldJd = -1; break; default: FieldTypeStr = "Char"; FieldLength = fdlength; FieldJd = -1; break; } //end switch if (FieldJd != -1) { f_WriteLine(a7u(fdName , "," , FieldTypeStr , "," , short2cp(FieldLength) , "," , short2cp(FieldJd))); } else { f_WriteLine(a5u(fdName , "," , FieldTypeStr , "," , short2cp(FieldLength))); } }//end for //f_WriteLine("BSM,Integer,12"); //Integer/Char/Float/Date //f_WriteLine("TBYBH,Char,9"); //f_WriteLine("TBBH,Char,9"); }//end for f_WriteLine("TableStructureEnd"); f_WriteLine(""); //this.ShowMsg(this, "完成输出表结构信息"); rbc = true; return rbc; } //输出线坐标 OK bool OutVctLineBeginEnd(short ai) { bool rbc = false; char* GeoVctType; char* FeatClassName; char* FeatureCode; char* FeatureName; //this.ShowMsg(this, "正在输出线坐标信息"); //输出线坐标 long i,j,ln; D_DOT *lxy; long len; double x,y; char* tmp; short LayerType; char* Lfilename; f_WriteLine("LineBegin"); //获取线层总个数 _GetLinNum(ai,&i,&ln); //获取工作区文件名 Lfilename=_GetAreaFileName(ai); //获取图层类型 LayerType=_GetAreaType(ai); GeoVctType="UnKown"; switch(LayerType) { case REG: GeoVctType="Polygon"; break; case LIN: GeoVctType="Line"; break; case PNT: GeoVctType="Point"; break; default: break; } FeatClassName =Lfilename; FeatureCode =a3u(GeoVctType,"_",i2cp(1)); FeatureName =FeatureCode; if (GeoVctType == "Line") { long lineNum = 0; for(i=1;i<ln;i++) //线要素总个数 { if(_GetLin(ai,i,&lxy,&len,NULL)<=0) { continue; } lineNum += 1; f_WriteLine(l2cp(lineNum)); //第一行序号(目标标识码) f_WriteLine(FeatureCode); //第二行要素代码 f_WriteLine(FeatureName); //第三行要素名称 f_WriteLine("1"); //第四行几何要素段类型(1=折线,2=圆弧,3=圆,4=椭圆,5=光滑曲线,6=B样条曲线,100=间接坐标) if(len>0) { f_WriteLine(l2cp(len)); //第五行点数 for(j=0;j<len;j++,lxy++) //线节点 { //获取地理坐标 x=lxy->x; y=lxy->y; //写入到文本中 char* x_s=d2cp(x); char* y_s=d2cp(y); tmp=a2u(x_s,","); tmp=a2u(tmp,y_s); // f_WriteLine(tmp); //第五行点坐标 } } else { //几何对象为空时 f_WriteLine("0"); //第五行点数 } if (lineNum % 50 == 0) { f_FlushFile(); } } //end for }//end if f_WriteLine("LineEnd"); f_WriteLine(""); f_FlushFile(); //this.ShowMsg(this, "完成输出线坐标信息"); rbc = true; return rbc; } //输出属性值 ?? bool OutVctAttributeBeginEnd(short ai) { bool rbc = false; char* GeoVctType = ""; char* FeatClassName = ""; char* FeatureCode = ""; char* FeatureName = ""; long i,ln; //this.ShowMsg(this, "正在输出属性值信息"); //输出属性值 f_WriteLine("AttributeBegin"); //获取线层总个数 _GetLinNum(ai,&i,&ln); //获取工作区文件名 char* Lfilename=_GetAreaFileName(ai); //获取图层类型 short LayerType=_GetAreaType(ai); GeoVctType="UnKown"; switch(LayerType) { case REG: GeoVctType="Polygon"; break; case LIN: GeoVctType="Line"; break; case PNT: GeoVctType="Point"; break; default: break; } FeatClassName =Lfilename; FeatureCode =a3u(GeoVctType,"_",i2cp(1)); FeatureName =FeatureCode; f_WriteLine(FeatClassName); if (ln > 0) { //有要素个数>0 char* AttrLine =NULL; char* fdval = NULL; long lineNum = 0; CATT_STRU *stru=NULL; char *att=NULL; short numbfield=0; short ret_gt; // for(i=1;i<ln;i++) //线要素总个数 { lineNum += 1; AttrLine = l2cp(lineNum); //获取属性值 ret_gt=_GetAtt(ai,LayerType,i,&stru,&att); short numbfield=stru->hd.numbfield; //字段个数 CFIELD_HEAD *pfd=stru[0].fldEntry; // for(int k=0;k<numbfield;k++) { fdval=&att[k]; //?? //_CvtFldToString(stru,att,k,fdval,pfd[k].msk_leng,NULL,NULL); if (strlen(AttrLine) <= 0) { AttrLine = fdval; } else { AttrLine =a3u(AttrLine,"," , fdval); } } f_WriteLine(AttrLine); if (lineNum % 50 == 0) { f_FlushFile(); } } } f_WriteLine("TableEnd"); f_WriteLine(""); f_WriteLine("AttributeEnd"); //this.ShowMsg(this, "完成输出属性值信息"); rbc = true; return rbc; } //获取线层总个数 /* _GetLinNum(ai,&i,&ln); f_WriteLine(UnionStr("线总个数:",i2cp(ln))); for(i=1;i<ln;i++) { if(_GetLin(ai,i,&lxy,&len,NULL)<=0) { continue; } WriteToLog(i2cp(i)); //线序号 f_WriteLine(i2cp(i)); for(j=0;j<len;j++,lxy++) //线节点 { //获取窗口坐标 //x=(int)((lxy->x-wxy->x0)*wxy->r); //y=(int)((lxy->y-wxy->y0)*wxy->r); //获取地理坐标 x=lxy->x; y=lxy->y; //写入到文本中 char* x_s=d2cp(x); char* y_s=d2cp(y); tmp=a2u(x_s,","); tmp=a2u(tmp,y_s); WriteToLog(tmp); f_WriteLine(tmp); } }*/ //删除线图层的所有要素 //_GetLin(lai,&logN,&n); //for(i=1;i<n;i++) /*编者注：图元实体号应从1开始*/ // _DelLin(lai,i); /*应该是for(i=1;i<=n;i++)*/ //走钟光标 //void WINAPI _StartClockCursor(void); //void WINAPI _GoingClockCursor(void); //void WINAPI _EndClockCursor(void); //void WINAPI _StartWait(void); //开始显示等待光标 //void WINAPI _EndWait(void); //结束显示等待光标 void Vct_Export() { _StartClockCursor(); _GoingClockCursor(); _StartWait(); // short n; WriteToLog_f_delete(); //TestGetDiskSpace(); //取打开的工作区数 n=_GetAreaNumber(AHInst); char* n_str=l2cp(n); ShowMsg(n_str); //获取当前可执行文件目录 char* curdir=GetCurDir(); ShowMsg(curdir); //WriteToLog(curdir); //获取显示工作区表集合 short *ptALst=NULL; if((ptALst=_GetDispAreaListAddr(GrpWnd))!=NULL) { if(ptALst==NULL) return; ShowMsg("工作区表不为空"); WIN_ORG tmp_wxy; short ai,LayerType; char* Lfilename; if(!_GetWinDispParam(GrpWnd,&tmp_wxy)) { //获取窗口显示参数 return; } WIN_ORG *wxy=&tmp_wxy; char* total_str=i2cp(ptALst[0]); char* tt_str=a2u("工作区表总个数:",total_str); ShowMsg(tt_str); //vct输出 char* t_fn="dgx.vct";//getfilename_lc(Lfilename,"\\"); WriteToLog(t_fn); char* tmp_line_path=a2u("c:\\",t_fn); WriteToLog(tmp_line_path); // f_SetFilePath(tmp_line_path); f_OpenFile_w(); OutVctHeadBeginEnd(); //1输出VCT文件头段 OutVctFeatureCodeBeginEnd(ptALst); //2输出VCT要素代码段 OutVctTableStructureBeginEnd(ptALst); //3输出VCT表结构段 //ptALst[0]表示总个数 for(int k=1;k<=ptALst[0];k++) { //获取一个工作区 ai=ptALst[k]; //获取工作区文件名 Lfilename=_GetAreaFileName(ai); //------------------------------- //获原文件路径的文件名 char* t_dir=getdir_lc(Lfilename,"\\"); ::MessageBox(NULL,"执行完毕!","",0); //------------------------------- //获取图层类型 LayerType=_GetAreaType(ptALst[k]); //输出VCT点,线，面图层的几何图形 switch(LayerType) { case LIN: //线图层 { WriteToLog("线图层"); WriteToLog(Lfilename); // OutVctLineBeginEnd(ai); //输出vct线要素几何图形 // } break; case REG: //面图层 { WriteToLog("面图层"); } break; case PNT: //点图层 { } break; default: break; } //end switch //输出属性表 OutVctAttributeBeginEnd(ai); }//end for f_CloseFile(); } _EndWait(); _EndClockCursor(); ShowMsg("生成完毕!"); }
--the---end----
--vp:hsg
--create date:2012-04-18---
110
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
