ArcMap与ArcScene实战：利用影像和DEM数据快速构建三维场景的完整指南

1次阅读

没有评论

共计 2080 个字符，预计需要花费 6 分钟才能阅读完成。

刚接触 GIS 三维开发的朋友，常常会在构建场景时遇到几个典型问题：

数据格式混乱 ：下载的影像和 DEM 数据可能来自不同来源，格式五花八门（如 GeoTIFF、IMG、ASCII 等），直接加载报错
坐标对不齐 ：当影像叠加到 DEM 上时出现偏移，常见于未统一地理坐标系（如 WGS84 与 CGCS2000 混用）或缺少投影定义
性能卡顿 ：加载全市域高分辨率数据时，ArcScene 卡成幻灯片，缺乏有效的金字塔管理和显示优化

ArcMap 主攻二维预处理
更适合数据格式转换、坐标系定义、影像裁剪等前期工作
通过【数据管理工具】完成 DEM 重采样（Resample）和影像校正（Warp）
ArcScene 专注三维表达
优势在于实时三维浏览和光照模拟
支持 TIN 表面生成和 3D 符号化，但大规模数据处理建议先在 ArcMap 预处理

影像校正
使用【栅格校正】工具消除几何畸变
关键参数：控制点至少选取 4 个，RMS 误差控制在 0.5 个像元内
DEM 重采样
通过【重采样】工具匹配影像分辨率
地形起伏大的区域建议用【双线性插值】，平原区可用【最邻近法】

在 ArcCatalog 中右键数据集→属性→坐标系
使用【投影】工具统一为相同坐标系（推荐 UTM 带号与中央经线匹配）
注意：DEM 高程单位需与水平坐标单位一致（如都是米）

ArcScene 中加载预处理后的 DEM
右键图层→属性→基本高度→选择【浮动在自定义表面】，指定 DEM 数据
叠加影像：直接将校正后的栅格拖入，设置相同高程参考
垂直夸大：地形起伏不明显时，在场景属性中设置 Z 值增强（通常 2 - 3 倍）

# -*- coding: utf-8 -*-
import arcpy
from arcpy.sa import *

# 环境设置
arcpy.env.workspace = "./input_data"
arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(32650)  # WGS84 UTM 50N

# 1. DEM 重采样
orig_dem = "raw_dem.tif"
resampled_dem = "dem_resampled.tif"
arcpy.Resample_management(orig_dem, resampled_dem, "10 10", "BILINEAR")

# 2. 影像校正
# 假设已有控制点文件 gcps.txt
corrected_img = "image_rectified.tif"
arcpy.Warp_management("raw_image.jpg", "gcps.txt", corrected_img, "POLYORDER1")

# 3. 创建三维场景
arcpy.CheckOutExtension("3D")
scene = arcpy.mapping.CreateNewScene("my_scene.sxd")

# 添加 DEM 并设置基本高度
dem_layer = arcpy.mapping.Layer(resampled_dem)
arcpy.mapping.AddLayer(scene, dem_layer)
dem_layer_obj = arcpy.mapping.ListLayers(scene, "dem_resampled")[0]
dem_layer_obj.setBasicElevation("FLOAT_ON_SURFACE", resampled_dem)

# 叠加影像
img_layer = arcpy.mapping.Layer(corrected_img)
arcpy.mapping.AddLayer(scene, img_layer)

# 保存场景
arcpy.mapping.SaveToLayerFile(scene, "final_scene.lyr")