利用VB 6.0实现地理国情普查分幅图元数据生产
0 引 言
为全面掌握我国地理国情现状,满足经济社会发展和生态文明建设的需要,国务院决定开展全国地理国情普查工作。多尺度数字高程模型(Digital Elevation Model,DEM)是地理国情普查中基础数据的重要组成部分,由数字高程模型派生形成的坡度、坡向等数据是地形地貌统计分析的基础。数字高程模型制作完成后,需对其元数据进行编辑,用以说明该数字高程模型数据的内容、质量、状况和其他有关特征的背景信息。其意义在于,元数据可以帮助数据生产单位有效地维护和管理数据,提供数据生产单位的各种有关信息供用户查询,帮助用户了解数据,以便用户处理和转换数据。根据对GDPJ 08—2013《多尺度数字高程模型生产技术规定》中分幅元数据表的分析可知,一幅数字高程模型数据对应一个元数据文件。传统的元数据编辑方式完全采用人工判读填写模式,不仅作业效率低下,而且容易产生编辑错误;为了有效解决地理国情普查数字高程模型数据元数据快速生产问题,本文通过对分幅元数据表内容的分析,利用Microsoft Visual Basic 6.0 ,提出一种根据图名图号关系表快速生产分幅元数据表的方法。
1 元数据表分析
DEM元数据由数据源情况、数据基本情况、生产过程情况、数据质量情况和数据分发信息5部分组成。
1)数据源情况为生产该DEM所用的数据源信息,其中,主要数据源中央子午线、主要数据源分带方式、主要数据源图名、主要数据源图号与DEM数据是相一致的。而主要数据源描述、主要数据源所采用的大地基准、主要数据源投影、主要数据源坐标单位、主要数据源高程系统名、主要数据源高程基准、主要数据源生产单位等可统一填写。
目前,B公司的经销商窜货现象较多,无论对市场还是对公司正常运营都造成了一定程度的影响。一方面,营销费用支持的不合理导致邻近市场经销商进货价格差别较大,从而为窜货的出现埋下了可能性。B公司在促销费用支持时,会针对一部分产品支持进货激励,即在经销商进货价格基础上给予一定的折扣,而不同市场折扣可能会相差较大。在业绩压力较大的情况下,也为了谋求更多的利润,一部分经销商会主动发起窜货,将产品低价出货至邻近市场的门店,这时,当地市场的经销商为抵制窜货便降低了出货价、损失了利润,这些经销商若发现长期无法保证利润时,便会自动放弃对产品的经营,给企业发展造成影响。
2)数据基本情况包括图名、图号以及图廓角点经度范围、图廓角点纬度范围,西南、西北、东南、东北图廓角点XY坐标,中央子午线,高斯-克吕格投影带号等。图廓角点经纬范围以及中央子午线、投影带号可以通过图号计算得到,而图廓角点XY坐标可以通过对计算得到的经纬度进行高斯正算得到。数据基本情况还包括如数据名称、数据所有者单位名称、数据的来源项目、数据生产单位名称、椭球长半径、椭球扁率、所采用大地基准、地形图投影、坐标单位、高程系统名、高程基准、数据生产方式、高程内插方法等,这些信息在DEM中都是统一定义的,可以进行统一编辑。
3)生产过程情况包括格网排列数,西图幅、北图幅、东图幅、南图幅接边情况,这些信息需进行人工判读后才能填写,但接边情况在一个测区范围内部都是已接边,只需挑选出测区外部未接边图幅进行人工修改编辑。
4)数据质量情况包括作业员自查结果及主要问题,队级、院级、局级检查出的主要问题及处理意见,队级、院级、局级检查员,队级、院级、局级检查时间,院级检查结论、局级验收结论,修改情况及遗留问题,局级验收单位以及接边质量评价、数据质量评价、质量评价单位及时间等信息。这些信息需要根据每幅DEM在经过队级、院级、局级检查时进行人工编辑。
5)数据分发信息包括密级、数据的版本、分发介质、分发格式名称、分发单位、分发者单位名称及联系方式、元数据创建单位、时间等。这些信息一般在相同测区内是相同的,可以进行统一编辑。
综上所述,地理国情普查DEM分幅元数据表中的信息编辑方式可分为3种类型:第一种是在人工判读后进行编辑,如生产过程情况、数据质量情况;第二种是进行统一批量编辑,如数据基本情况中的数据名称、大地基准,数据源基本情况中的主要数据源高程基准、主要数据源所采用的大地基准,以及数据分发信息。第三种是DEM经纬度范围、图廓坐标、中央子午线、投影带带号等通过计算获得的信息。其中;DEM数据经纬度范围可通过图符号计算获得;根据得到的西南图廓经纬度,由高斯投影带计算可以得出中央子午线及投影带带号;根据2000国家大地坐标系高斯正算公式,可求出图廓角点X、Y坐标。
2 程序设计基本思路
通过对DEM分幅元数据表数据项的分析,本程序设计的基本思路是:1)读取图名图号对应关系Excel表,并将图名和图号保存到内存中;2)检查图号填写是否正确;3)读取元数据示例模板文件;4)根据图号计算西南图廓点经纬度及图幅经纬度范围;5)根据图号计算中央子午线及投影带号;6)根据2000国家大地坐标系高斯正算公式计算西南图廓点经纬度对应的XY坐标,并得到西南、西北、东南、东北图廓点X、Y坐标。
在土建工程中,如果混凝土结构发生变化,荷载会产生二次应力,从而导致混凝土产生裂缝,因此,设计时应更加合理。另外,很多因素都能够影响到混凝土的施工,通过不断的进行实践可以总结出,振捣浇筑能够提高混凝土自身的性能。在进行混凝土施工时,必须注意混凝土浇筑和振动应同时进行,以有效防止裂缝,一般情况下,通常使用分层浇铸。另外,在进行振捣时,为了防止因振捣不到位而造成施工裂缝的发生,必须要保障振捣的密实度。除此之外,不合格的模板也会导致裂缝的发生,在进行浇筑施工之前应检查模板的质量,同时还要在施工过程中有效的控制模板,进而可以有效的预防因模板而引起的混凝土开裂的现象。
程序设计流程如图1所示:
3 程序的编写
Visual Basic 6.0是一种面向对象的可视化程序设计语言,采用了事件驱动的编程机制,开发周期短。因此,本程序选择使用Visual Basic 6.0进行编程。
3.1 程序界面设计
本程序的操作顺序应为读取图名图号关系表、检查图号正确性、读取元数据表模版、开始元数据生成4个步骤。因此,本程序界面设计如图2所示。
3.2 图号图名文件的读取与检查
元数据表中的经纬度范围、图廓角点坐标等数据都是基于图号信息计算得出的,因此本程序使用通用对话框读取图名图号关系表,并将其存入动态数组tuminginput(),tuhaoinput()中。图名图号关系表是人工录入的,存在图号录入错误的可能性,因此还需要进行检查。方法是:判断图号字符长度是否为10,其第2至3位及最后6位是否为数值形式,如果不是,则将该单元格颜色更改为红色。
3.3 图廓经纬度范围计算
图幅编号录入无误后,可根据地图分幅与编号理论,计算其西南图廓角点经纬度。公式为:
λ=(b-31)*6°+(d - 1) * Δλ
(1)
φ=(a-1)*4°+(4/Δφ- c)*Δφ
(2)
其中,a是1∶1 000 000图幅所在纬度带的字符所对应的数字码,b是1∶1 000 000图幅所在经度带的数字码;c、 d分别是该比例尺地形图在1∶1 000 000地形图编号后的行号、列号;Δφ、Δλ分别是该比例尺地形图分幅的纬差、经差。本文涉及的DEM数据都是1:10 000标准分幅的,其纬差Δφ=2′30″、经差Δλ=3′45″。a是常数6,b是图号第2至3位,c是图号第5至7位,d是图号8至10位。可建立一个分解图号的程序过程,并使用MId函数对其进行提取,再利用公式(1)(2)进行计算。主要代码如下:
Call de(jd) ‘将经度转换度、分、秒;其变量名分别为d1、d3、d4
充分地内部调研和外部调研,最终目的是为系统 《需求方案》的合理性、充分性和完整性负责,这个过程是允许循环迭代的,充分地前期调研论证是为保证后期项目有序推进的重要基础;立项、招标过程是将《需求方案》转化为《需求实现方案》的重要手段,务必要谨慎、字斟句酌,确保最终的服务提供者和服务方案能够将需求转化为现实;以此来保证系统最终的成功运行和顺利验收。
Dim st As String
st=tuhaoinput(i)
Call fjtuhaost
人才问题是制约乡村振兴的重要制约因素,城乡发展不平衡不协调时长期困扰我国发展的重要问题,由于城市与乡村资源分配的不均衡导致优质的要素流向城市,大量优秀人才流向城市,乡村人才匮乏。实施乡村振兴战略,需要各领域、各方面的人才支持。塑造乡士人才应注重对当地人才进行发掘培养,加强农村专业人才队伍的建设,注重人才引进,纳四方贤才。发挥科技人才支撑作用,鼓励社会各界投入乡村建设。
jd=(ib-31)*6+(id-1)*0.062 5′0.062 5=3.75/60′jd代表经度,单位是度
据沈崇德介绍,澳新风险管理标准是澳大利亚/新西兰联合制定的世界上第一个国家风险管理标准,已被澳大利亚政府和世界上很多企业广泛采用,成为争相效仿的标杆,近年来,被应用到医疗安全管理方面,成效斐然。据了解,无锡市人民医院是率先将这一模式引入后勤患者安全体系的医院之一。
wd=((ia-1)*4)+(4/(150/3 600)-ic)*(150/3 600)′wd代表纬度,单位是度
End Sub
与其他脏腑的关系:督脉、任脉和冲脉,三者“一源三歧”,均起于胞中,下出于会阴。会阴,即外生殖器后方与肛门前方的部位。足太阳膀胱经经别入于肛,所以针灸上,常常通过足太阳经和任督冲脉的穴位治疗肛门疾病。在脏腑中,魄门的开合由心神主宰,与前阴同为肾之窍。饮食糟粕的排泄不仅关乎于肾,而且与脾之运化、肺之肃降,以及肝之疏泄均有密切关系,所以中医又称:“魄门亦为五脏使。”
西南图廓角点经纬度计算得出后,对经度增加3′45″得到经度范围,对纬度增加2′30″得到纬度范围。需要注意的是,经纬度为60单位进制,在程序中,应从秒位开始相加,并判断是否>60,是则向分位进一,否则不进位;再累加分位差值,判断是否>60,是则向度位进一,否则不进位。以图廓角点经度范围为例,主要程序代码如下:
Public Function jdfw(jd As Double) As String
Public Sub creatjwd(i)
d4new=Val(d4)+45
If d4new> 60 Then d4new=d4new-60∶d3new=d3+3+1
If d3new > 60 Then d3new=d3new-60∶d1new=d1+1
Else
d1new=d1
End If
综上,笔者认为社会福利政策质量是评价社会福利政策合法性、科学性及目标性的合理程度以及政策执行效果的标准。高质量的社会福利政策具有政策过程完整、执行效率和目标群体满意度高的特点,能够提高广大社会成员的物质和精神生活水平,符合大多数公民的根本利益。
Else
小学全科师范生党员学习党的时代性思想加上教师本身自主的学习教师新时代新思想,两者相互促进、相互影响,更好的践行社会主义核心价值观。
If d3new > 60 Then d3new=d3new-60∶d1new=d1+1
Else
d1new=d1
End If
End If
…….
3.4 图廓角点X、Y坐标计算
全国地理国情普查DEM数据采用2000国家大地坐标系。根据国家测绘局《现有测绘成果转换到2000国家大地坐标系统技术指南》的规定,2000国家大地坐标系的平面坐标仍采用高斯-克吕格投影,投影方法、投影分带、中央子午线等均与原有的坐标系统一致,但椭球参数及子午线长计算公式各系数精度均有所变动。
本文在获得西南、西北、东南、东北图廓角点经纬度后,可利用2000国家大地坐标系高斯正算公式进行X、Y坐标计算。公式如下:
4.开展美体教育。蔡元培所说的“美术”是广义的,包括唱歌、手工、图画、文字趣味等。“注意美术,足以生美感,既生美感,自不致苦脑力。且美术更有足重者。”也就是要通过唱歌、手工、图画、美文等来对师范生进行美感教育,起到开发智力、陶冶情操、协调身心的作用。蔡元培一直重视学生的身体素质,他告诫学生,受教育,“该养成普遍而健全的能力”;多开展体育活动,“使身体完全发达”。开展美育和体育活动,以此为桥梁,促进师范学生素质的全面发展。这是蔡元培“五育并举”教育思想最富有特色的部分。
(3)
其中,子午弧长计算公式为:
X=a(1-e2)(A′arcB-B′sin2B+C′sin4B-D′sin6B+E′sin8B-F′sin10B+G′sin12B)
(4)
据业内人士介绍,水利基础设施建设投资步伐的加快将催生一批适合水利施工的专用机械,大型水利施工中的成套设备、水下抢险施工中的遥控机器人等是集机械、电子和自动化控制等一系列技术于一体的集成系统,将促进工程机械和专用水利施工机械智能化程度的快速提高。与此同时,在大型水利项目建设和农田水利改造过程中,挖掘、铲运和平土等工作将促进挖掘机、推土机、各类钻机和装载机等设备的需求增加,同时对其衍生品的专业化程度提出更高要求,也将带动工程机械配件行业的发展。
为使计算过程中的中央子午线、高斯投影分带带号以及计算出的X、Y坐标在程序任何模块中可读取,需将相应变量定义为公共变量。主要代码如下:
Public Sub XY(bl As Double, ll As Double)
(4)式中:
Dim a1 As Double, b1 As Double, c1 As Double, e1 As Double, f1 As Double, g1 As Double
Dim l As Double, l0 As Double, xx As Double, n As Double, t As Double, yy As Double
Const a=6378137, e=0.0818191910428, e2=0.0066943800229
Const e12=0.00673949677548, pi=3.1415926535, p=206264.80625992
a1=1+3*e2/4+45*e^4/64+175*e^6/256+11025*e^8/16384+43659*e^10/
65536+693693*e^12/1048576
b1=3*e2/8+15*e^4/32+525*e^6/1024+2205*e^8/4096+72765*e^10/
131072+297297*e^12/524288
c1=15*e^4/256+105*e^6/1024+2205*e^8/16384+10395*e^10/65536+
1486485*e^12/8388608
d3new=d3+3
紧凑型核聚变反应堆未来有望应用在飞机、火箭、舰艇、车辆等装备以及发电、海水淡化等领域,为军事装备、能源和交通运输等领域带来革命性变化。在军事上的应用望催生出超长续航、上天入地、无所不能的新型运载平台,模糊海陆空天潜的界限,甚至使得跨星系航行也成为轻而易举的事。
e1=315*e^8/131072+3465*e^10/524288+99099*e^12/8388608
f1=693*e^10/1310720+9009*e^12/5242880
g1=1001*e^12/8388608
d1=35*e^6/3072+105*e^8/4096+10395*e^10/262144+55055*e^12/1048576
b=bl*pi/180 :dh=Round(ll / 3)
l0=dh*3 ‘l0为中央子午线:l=(ll-l0)*3600/p
n=a/Sqr(1-e2*(Sin(b))^2) :t=Tan(b)
yy=e12*(Cos(b)^2)
xx=a*(1-e2)*(a1*b-b1*Sin(2*b)+c1*Sin(4*b)-d1*Sin(6*b)+e1*Sin
(8*b)- f1*Sin(10*b)+g1*Sin(12*b))
x=xx+n*t*(Cos(b)^2)*l^2/2+n*t*(5-t^2+9*yy+4*yy^2)*(Cos(b)^4)*
l^4/24+n*t*(61-58*t^2+t^4)*Cos(b)^6*l^6/720
y=n*Cos(b)*l+n*(l-t^2+yy)*(Cos(b)*3)*l^3/6+n*(5-18*t^2+t^4+14
*yy-58*yy*t^2)*(Cos(b)^5)*l^5/120
End Sub
无通信互联线储能系统的直流母线协调控制策略//黎阳,王林,黄辉,张海龙,姚为正,凤勇//(1):118
最后,编写程序读取元数据模版,将计算出的经纬度范围、中央子午线、高斯投影带带号、图廓角点XY坐标,以及其他归属于统一批量编辑的元数据属性内容写入以图号命名的Excel文件,即完成了DEM元数据的批量化生产。
5.根据对假设五检验得出的结论,笔者认为,除了相关准则要完善对企业无形资产研发支出资本化的规定,会计师事务所在对企业进行审计时,也要更加注意该方面内容的审计,重视报表附注对开发支出的披露,尤其是对报表附注披露不规范的情况实施相应的审计程序,从第三方角度进行监督管理,以避免企业利用准则的漏洞做出一些违规舞弊行为,扰乱经济秩序,误导报表使用者做出错误判断。
4 结 语
本文方法应用于墨江市地理国情普查DEM分幅元数据生产项目,10 min总计批量生产出512份元数据成果。通过验证发现央子午线、高斯投影带带号、图廓经纬度范围、图廓角点XY坐标等信息正确无误,从而证明了2000国家大地坐标系高斯正算公式及程序的正确性和精度的可靠性。本文所述方法提高了地理国情普查DEM分幅元数据表的生产效率。
参考文献:
[1] 赵英志,刘永涛,郑玉军.利用VB6.0实现2000国家大地坐标系高斯正反算程序的编写[J].测绘通报, 2010(5):38-41.
[2] 刘平,罗伊,范海生.浅谈地球坐标系及2000国家大地坐标系[J].地理空间信息,2010(3):106-108.
[3] 宋华伟,周莎,马德涛.浅议地理信息元数据[J].地理空间信息,2008, 6 (1) :74-76.
[4] 孔祥元,郭际明,刘宗泉.大地测量学基础[M].武汉:武汉大学出版社,2001.