USGS DEM文件格式说明

转自http://blog.sina.com.cn/learnc

 表1

数据元素	类型	格式声明			说明
		ASCII格式	起始字节	结尾字节
文件名	ALPHA	A40	1	40	在DEM文件头中不应出现除美国外的其他国家如墨西哥、加拿大等国家的名称
自由格式文字	ALPHA	A40	41	80	有关数字化过程的有用信息，例如仪器、相片好、扫描孔径等
占位符			81	109	空格
东南角地理坐标	INTEGER*2, REAL*8	2(I4,I2,F7.4)	110	135	东南角地理矩形坐标顺序为： x = Longitude = SDDDMMSS.SSSS y = Latitude = SDDDMMSS.SSSS (负号右对齐，不能以零起头，而增加正号)
处理方法代号	ALPHA	A1	136	136	1 采用自动相关和简单的双线性重采样算法 2 人工采集和GRIDEM简单双线性算法 3 DLG/或地形图CTOG 8向线性内插 4 从航测系统采集的等高线数据4向内插而来 5 DLG/地形图 LINETRACE，LT4X 复杂线性内插 6 DLG/地形图 CPS-3，ANUDEM，GRASS 复杂多项式 7 电子影像（非摄影测量）被动或主动传感器系统
占位符			137	137	空格
组合指示	ALPHA	A3	138	140	这个参数是为30分DEM准备的
生产者代码	ALPHA	A4	141	144
Level等级码	INTEGER*2	I6	145	150	1=DEM-1 2=DEM-2 3=DEM-3 4=DEM-4
DEM高程点图形定义码	INTEGER*2	I6	151	156	1=规则网格 2=随机分布
水平参考系定义代码	INTEGER*2	I6	157	162	0=地理 1=UTM投影 2=国家平面
UTM分带代码	INTEGER*2	I6	163	168	UTM 分带代码， 如果采用地理坐标这项参数为0
地图投影参数	REAL*8	15D24.15	169	528	对UTM及地理坐标系，该参数全部置0
平面参考坐标系单位代码	INTEGER*2	I6	529	534	0=弧度  1=英尺 2=米  3=弧秒 一般对于7.5分DEM，使用代码2。而对于30分、1度和Alaska DEM一般采用代码3。
高程坐标单位	INTEGER*2	I6	535	540	1=英尺 2=米 一般7.5分，30分，1度和Alaska DEM采用代码2
DEM边数	INTEGER*2	I6	541	546	通常为4
4对实数表示DEM边界	REAL*8	4(2D24.15)	547	738	坐标以(东经,北纬)为记录顺序,四个角以西南角为起点按顺时针排列。
两个元素分别表示整个数据中DEM的最大和最小值	REAL*8	2D24.15	739	786	DEM以元素9的代码表示的单位为单位，遵从元素6的分带划分。
以弧度作为记录单位，逆时针记录DEM参考系与地平面轴的夹角	REAL*8	D24.15	787	810	如果采用元素5中代码定义的坐标系统，这项参数为0
精度代码	INTEGER*2	I6	811	816	0=精度未知 1=精度信息在type C 中表达
X、Y、Z三方向分辨率，单位参考元素8，9	REAL*4	3E12.6	817	852	X、Y方向的分辨率只允许取整数。对所有的USGS DEM，除1度DEM以外，Z分辨率如采用英尺为单位，则最多允许小数点后1位，以米为单位最多允许小数点后2位。一些典型的表示如对于7.5分DEM可以使用30,30,1；和10,10,1等
DEM中包含的断面文件行列数（m,n）	INTEGER*2	2I6	853	864	当行数m=1时，n等与DEM的列数
旧格式在此结束
最大计曲线间隔	INTEGER*2	I5	865	869	仅当存在两种以上计曲线间隔时存在
源曲线间隔单位	INTEGER*1	I1	870	870	对应于源曲线最大间隔单位 0=N.A 1=英尺，2=米（仅适应于 level 2DEM）
最小计曲线间距	INTEGER*2	I5	871	875	最小或仅有的计曲线间隔
源曲线间隔单位	INTEGER*1	I1	876	876	对应于地图最小计曲线间隔单位  1=英尺 2=米
数据生产日期	INTEGER*2	I4	877	880	“YYYY”一个4字符表示的年份（例如 1975，1997等）即可以是编绘日期也可以是航摄日期
数据检查和生产日期	INTEGER*2	I4	881	884	“YYYY”一个4字符表示的年份，可以是完成日期或者修测日期
检查标志	ALPHA*1	A1	885	885	“I”表示第三部分包括的所有质量控制均已实现
数据确定性标志	INTEGER*1	I1	886	886	0=没有做过确定过程 1=利用检查点计算过RMSE，未进行DEM编辑 2=利用检查点计算过RMSE，批处理水体编辑3=利用检查点计算过RMSE，编辑包括水体编辑 4=Level 1 DEM 编辑，水体编辑，利用检查点计算过RMSE。 5= Level 2 和Level 3 DEM编辑。利用检查点计算过RMSE。
可疑和空旷标志	INTEGER*1	I2	887	888	0=无 1=可疑区域 2=空旷区域 3=可疑区域和空旷区域
高程基准	INTEGER*1	I2	889	890	1=局域平均海平面 2=国家高程测量基准 1929（NGVD 29） 3=北美垂直基准（NAVD88）
水平基准	INTEGER*1	I2	891	892	1=北美基准 1927（NAD27） 2=世界测量系统1972（WGS72） 3=WGS 84 4=NAD 83 5=旧 夏威夷基准 6=波多黎各基准
数据编辑	INTEGER*2	I4	893	896	01-99最初为DMA特别设置区（USGS一般设为01）
空旷比例	INTEGER*2	I4	897	900	如果元素25指出存在空旷区域，这个区域（右对齐）说明文件存在相应比例空旷区点（-32,767）
边匹配标志	INTEGER	4I2	901	908	边匹配标志，按照西、北、东和南顺序排列
垂直标准位移	REAL*8	F7.2	909	915	垂直标准位移，形式为SFFF.DD 取值来自于VERTCON计划的矩形四条边的平均位移。通常加此值后转化为VAVD88

 表2

数据元素	类型	格式声明			说明
		ASCII格式	起始字节	结尾字节
一个二维数组表达当前断面文件的序号	INTEGER*2	2I6	1	12	行、列的取值范围均为1,m和1,n。行数通常设为1，列通常表示当前的断面序数
一个二维数组表示DEM的断面数	INTEGER*2	2I6	13	24	一个参数对应断面文件的行数，第二个参数通常设为1
一个二维数组表征断面文件的起点坐标	REAL*8	2D24.15	25	72
断面文件的局部高程基准	REAL*8	D24.15	73	96	单位采用逻辑记录A中元素9的定义
一个二维数组表达当前断面文件中DEM的最大最小值	REAL*8	2D24.15	97	144	单位采用逻辑记录A中元素9的定义
DEM数据实体	INTEGER*4	mn(I6)			6×(146 or 170),146第一个数据块最大数据数，170后续数据块为了节约存储空间，无论断面文件有多长，DEM都以1024为单位存储，不足1024的以空格补齐。这样的一个数据块所能容纳的数据量为（1024-144）/6=146，后续数据块为所能容纳的数据量为1024/6=170。

 

 

表3

数据元素	类型	格式声明			说明
		ASCII格式	起始字节	结尾字节
表明统计存在性的代码	INTEGER*2	I6	1	6	1=有统计 0=无统计
相对于绝对平面的均方根误差RMSE	INTEGER*2	3I6	7	24	RMSE的单位跟逻辑记录A中8，9
RMSE取样数	INTEGER*2	I6	25	30	如果此值为0则说明，RMSE来自于估计而不是计算
元素5是否存在	INTEGER*2	I6	31	36	1=存在 0=不存在
相对于绝对平面的均方根误差RMSE	INTEGER*2	3I6	37	54	RMSE的单位跟逻辑记录A中8，9
元素5的取样数	INTEGER*2	I6	55	60	如果此值为0则说明，RMSE来自于估计而不是计算