USGSDEM——USGS ASCII DEM(和CDED)

司机简称

USGSDEM

Driver built-in by default

This driver is built-in by default

GDAL支持读取USGS ASCII-DEM文件。这是美国地质勘探局在被SDT取代之前使用的传统格式,也是加拿大CDED DEM数据产品使用的格式。应支持USGS DEM文件上最流行的变体,包括正确识别坐标系和地理参考定位。

7.5分钟(UTM网格)USGS DEM文件的边缘通常会有丢失数据的区域,这些区域会用nodata值正确标记。USGS DEM文件中的高程值可以以米或英尺为单位,这将由GDALRasterBand::GetUnitType()的返回值指示(“m”或“ft”)。

请注意,USGS DEM文件表示为一个大平铺。如果GDAL平铺缓存大小较小,这可能会导致缓存抖动问题。当读取第一个像素时,它也会导致相当大的延迟,因为整个文件将被摄取。

实现usgsdemdataset.cpp的一些代码是由Ben Discoe从VTP代码派生的。见 Virtual Terrain 有关VTP的详细信息。

驱动程序功能

Supports Create()

This driver supports the GDALDriver::Create() operation

Supports Georeferencing

This driver supports georeferencing

Supports VirtualIO

This driver supports virtual I/O operations (/vsimem/, etc.)

创建问题

GDAL支持地理(和UTM)USGS DEM和CDED数据文件的导出,包括根据加拿大联邦政府规范生成CDED 2.0 50K产品的能力。

输入数据必须已在地理或UTM坐标系中采样。默认情况下,输入文件的整个区域将被输出,但对于CDED50K产品,输出文件将在产品指定的分辨率和产品平铺边界上进行采样。

如果输入文件具有适当的坐标系信息集,则导出到特定的产品格式可以接受不同坐标系中的输入(即从Albers投影到NAD83 geographic,用于CDED50K生产)。

创建选项:

  • PRODUCT=DEFAULT/CDED50K :当指定CDED50K时,输出文件将强制遵守CDED50K产品规范。输出将始终为1201x1201,通常是15分钟乘15分钟的平铺(尽管在遥远的北方地区经度更宽)。

  • TOPLEFT=long,lat :对于CDED50K产品,这用于指定要生成的磁贴的左上角。它应该在15分钟的边界上,可以用十进制度数或度数和分钟表示(例如左上角=117d15w,52d30n)。

  • RESAMPLE=Nearest/Bilinear/Cubic/CubicSpline :设置用于将数据重采样到目标网格的重采样内核。只有在生产CDED50K等特定产品时才有效果。默认为双线性。

  • DEMLevelCode=integer DEM级别(1、2或3,如果设置)。默认为1。

  • DataSpecVersion=integer :数据和规范版本/修订版(如1020)

  • PRODUCER=text :在生成的文件的producer字段中最多输入60个字符。

  • OriginCode=text :在生成的文件(YT代表育空)的源代码字段中最多输入4个字符。

  • ProcessCode=code :要放入生成文件的处理代码字段中的一个字符(8=ANUDEM,9=FME,A=TopoGrid)。

  • TEMPLATE=filename :对于任何输出文件,都可以指定模板文件。许多字段(包括数据生成器)将从模板文件(如果提供)中复制,否则留空。

  • ZRESOLUTION=float :DEM将高程信息存储为正整数,这些整数使用“z分辨率”缩放。默认情况下,此分辨率写为1.0。但是,如果希望将整数缩放为浮点数,可以在此处指定不同的分辨率。

  • NTS=name :NTS映射表名称,用于派生左上角。只有在生产CDED50K等特定产品时才有效果。

  • INTERNALNAME=name :数据集名称写入文件头。只有在生产CDED50K等特定产品时才有效果。

示例:下面将生成单个CDED50K磁贴,从左上角为117w、60n的磁贴的更大DEM覆盖范围yk_3arcsec中提取。文件yk_template.DEM用于设置某些产品字段,包括数据生产者、过程代码和源代码字段。

gdal_translate -of USGSDEM -co PRODUCT=CDED50K -co TEMPLATE=yk_template.dem \
               -co TOPLEFT=-117w,60n yk_3arcsec 031a01_e.dem

注:执行为 gdal/frmts/usgsdem/usgsdemdataset.cpp .

GDAL中的USGS DEM读取代码是从 VTP 软件。出口能力是在育空环境部的财政支持下发展起来的。