DEM格式转换

在进行SAR数据处理时我们需要用到DEM数据，这里面以Aster数据为例介绍在SARSCAPE和GAMMA下所需的数据格式的转换。

1、SARSCAPE下DEM数据格式转换

SARSCAPE下可以直接利用Digital Elevation Model Extraction下载DEM数据，这里不详细说了。我下载的DEM数据是TIFF格式的SARSCAPE下需要的是_DEM格式，这就需要DEM进行格式转换。

首先，如果我们研究区的范围较大，DEM先进行拼接。我的研究区需要6景DEM，如下图所示：

选择Mosaicking--->Georeferenced，把DEM数据导进去，如下图所示，然后点击apply生成拼接的DEM，这里生成的DEM是ENVI的标准格式，

其次，我们对生成的ENVI标准格式的DEM进行格式转换，SARSCAPE菜单下选择Import Data中的ENVI original，界面如下图所示

选择输入文件，在输出文件中会自动出现输出文件路径和输出文件名称，我的输出文件名自动是mosaic_bil，把文件名改为mosaic_DEM就可以了。Data Units我们选择DEM，点击start执行。就可以生成SARSCAPE所需要的_DEM数据。

在SARSCAPE中如果我们的DEM数据是TIFF格式的，我们在选择Import Data的数据类型需要选择TIFF，界面如下图所示，这里面可以选择

有地理编码和无地理编码的，我们选择有地理编码，点击start执行

2、GAMMA下DEM数据格式转换

1. 数据格式转换

   ERDAS下选择import工具栏，如下图1所示，输出数据，数据格式为GenericBinary，点击OK后，弹出图2对话框。

图1

图2

在图2的对话框中，选择dataformat为BIL，这里选择默认。data Type这里选择Signed 16 Bit，这里是根据图像自身信息获得，其信息如图3所示。 ByteOrder选择Motorola(Big Endian)，这个是gamma所要求的格式。除了上述信息可在图中得出，我们还需要记住下图中矩形框内的信息，这在gamma的下一步处理中将要用到。

图3

2、在Gamma中进行格式转换

在ERDAS得到.dem数据需要在Gamma中进行格式转换，利用swap_bytes函数进行格式转换，在命令窗口中输入如下命令：swap_bytes xiangge_16.dem xiangge_16_swap.dem 2 回车执行转换命令。

3、创建DEM参数文件

执行Create_dem_par xiangge_16_swap.dem_par命令，出现下图4窗口，要求输入DEM projection，根据DEM自身信息我们选择EQA（等角投影），在图中的参数中我们需要特别注意的是DEM width、DEM length、latitude and longitude posting、latitude andlongitude of first DEM sample这些参数需要根据每个图像自身信息输入，我们在ERDAS中所记录的信息就在这个地方用到 ，生成的参数文件如下所示：

                               图4

利用参数文件进行显示，输入命令disdem_par xiangge_16_swap.dem xiangge_16_swap.dem_par，显示图像如下所示：

特别说明的一点：如果要转换成其他投影，我们可以利用DEM_TRANS进行转换，这里不详细说了。同时如果DEM图像中有空洞，我们需要利用replace_values函数去除空洞