简介一份福建全省DEM原始高程数据TIFF格式面向GIS、遥感和测绘领域的学习者与工程师可用于地形地貌分析、水文模拟、灾害评估及规划等工作。压缩包共30个文件主体为按经纬度分块的25个zip子压缩包另含2个tif主数据文件、1个xml元数据、1个jpg预览图及1个ovr金字塔文件整体大小约362.01MB分块设计便于按区域加载与拼接。目前已有729人学习下载数据可直接导入ArcGIS等地理信息软件结合GeoTIFF坐标信息进行坡度、坡向、等高线提取、山脊山谷分析等操作。对需要福建省高精度地形底图的科研项目或工程应用而言这套数据提供了可追溯的原始高程网格配合配套预览和元数据可快速验证数据质量并开展后续分析。1. 福建DEM原始高程数据TIF拿到文件只是第一步打开福建某县的DEM TIF第一眼看到的可能不是连绵山脉而是一整屏黑块、花斑或是与矢量路网错开几百米的“漂移”。这不是数据损坏而是原始高程TIF自带的坐标系、NoData、分块和压缩方式还没有被正确解释。福建的原始高程数据多以TIF栅格分发常见的包括30米SRTM、12.5米ALOS和5米精度DSM每个像元都记录着海拔值但要让它们真正参与土方计算、流域分析、视域分析或三维建模必须经历读元数据、裁切、拼接、重投影、重采样这一整套流水线。这篇笔记写给正在做选址、水文和三维可视化却卡在“原始TIF”这一步的从业者。2. 读透福建DEM TIF元数据gdalinfo先于一切操作拿到任何DEM TIF我习惯先跑一遍gdalinfo而不是直接拖到ArcMap里看。原因很简单栅格文件是一个黑匣子外表都是.tif内部可能装着地理坐标、投影坐标、浮点或整型高程甚至装着土地分类代码。如果连它住在哪个坐标系、什么分辨率、哪个值是NoData都搞不清楚后面的裁剪、拼接和三维化都会连环翻车。2.1 用 gdalinfo 一眼看清坐标系、分辨率与NoData最省事的做法是使用GDAL命令行工具Windows、macOS和Linux都有相应安装方式。定位到数据目录后执行gdalinfo fujian_dem.tif输出里的关键行要逐条看Driver: GTiff/GeoTIFF Size is 24000, 18000 Coordinate System is: PROJCRS[CGCS2000 / 3-degree Gauss-Kruger CM 117E, BASEGEOGCRS[China Geodetic Coordinate System 2000, DATUM[China 2000]... Pixel Size (30.00000000, -30.00000000) NoData Value -9999 Band 1 Block256x256 TypeFloat32, ColorInterpGraySize对应行列数Pixel Size表示栅格分辨率这里30米TypeFloat32说明高程以单精度浮点保存NoData Value-9999说明无效高程被标成负值。Coordinate System段落会直接打印出坐标系名称这是后面一切空间操作的依据。如果输出里只有无坐标的提示说明这个TIF缺少投影定义需要立刻处理。为什么先说分辨率因为福建地形从沿海0米到武夷山2000多米不同来源的TIF分辨率相差很大如果后续要把多个来源的DEM叠加分析分辨率不一致会让坡度和汇水面积结果失真。先确认这些字段才能决定下一步用哪种方式裁切和重采样。这里也要养成一个习惯把gdalinfo的完整输出保存成文本文件随数据一起归档。后面项目复盘时你能清楚地知道用的是哪一套投影、哪个分辨率而不是靠记忆猜。很多团队协作翻车就是有人把WGS84地理坐标直接当成CGCS2000读取所有叠加分析全错最后只能从头核对源文件。2.2 福建常用的CGCS2000 3度带与WGS84 UTM选错就整体偏福建地理位置跨115°50′E到120°43′E在3度带投影下会横跨中央经线117E和120E两个分带在UTM分区下则主要落在50N和51N。国内生产的高程成果常用CGCS2000坐标系而很多从地图工具或开源数据源下载的TIF却写着WGS84地理坐标甚至是WGS84 UTM投影。处理福建数据时我几乎每次都要先问一句这数据是测绘部门出的还是卫星产品这个问题直接决定坐标系假设。遇到这种情况我的第一个判断标准是如果数据面是福建某个县市优先使用CGCS2000 3度带如果数据源是海外卫星产品多半是WGS84 UTM。三套常用参数的差异可以直接查表坐标系EPSG编号中央经线常见来源CGCS2000 / 3-degree Gauss-Kruger CM 117E4547117E福建中西部国土/测绘成果CGCS2000 / 3-degree Gauss-Kruger CM 120E4548120E福建东部沿海成果WGS 84 / UTM zone 50N32650117ESRTM、ALOS等公开DEM选错带的后果不是“数值不一样”那么简单。把4547的数据直接按4548加载同一位置的东西方向偏差可能达到数百米把WGS84地理坐标当成CGCS2000投影坐标叠加在边缘地区更是直接出图就偏。所以我在处理前会用一个笨办法验证用gdalinfo看坐标系的Datum字段再裁一个山顶点坐标与已知控制点对比确认无误才继续。如果原始TIF没有定义坐标gdalinfo会直接显示“Coordinate System is Unknown”此时不要急着用ArcGIS的“定义投影”去猜先查数据说明文件。2.3 像素深度与NoData为什么高程统计会出一堆负数DEM TIF的像素深度常见Float32和Int16两种。Float32能记录小数海拔适合做水文分析和坡度计算Int16则省空间但通常以厘米或毫米为单位保存直接读取时会得到夸张的数值。很多教程只告诉你“导入栅格”没提这两个类型的差别于是统计高程时出现一排-9999甚至-3.4e38。这种异常值不是数据坏了而是NoData没有被识别或者像素类型解释错了。另一个高频问题是如何处理NoData。原始TIF里的无效区域通常用-9999或-32768填充如果统计时没有排除NoData均值、最大值都会被污染。常见的做法是先用gdal_translate把NoData统一改成-9999再在ArcGIS或QGIS里设置环境掩膜。命令行如下gdal_translate -ot Float32 -a_nodata -9999 fujian_dem_raw.tif fujian_dem_f32.tif这段命令把原始TIF转为Float32并明确写入NoData值为-9999。参数-ot指输出像素类型-a_nodata强制设置无效值。要注意如果原始NoData是0而0又是沿海海平面真实值绝不能顺手设成NoData否则会把沿海陆地全部挖空。先看gdalinfo输出里NoData字段再决定。在ArcGIS中还可以在图层属性里勾选“忽略NoData”但这只影响显示和统计不会改写文件本身如果后续要交给其他软件还是要写进文件里才稳妥。3. 裁剪福建DEM TIFArcMap提取工具与GDAL命令两种路线有了元数据下一步就是把覆盖全省或全市的TIF裁到目标范围。这个环节最常见的痛点不是“不会裁”而是裁完之后的边界、像元对齐和分辨率变了导致后续分析对不上。所以我会把ArcMap图形化操作和GDAL命令行分别讲清楚并说明各自的适用场景。在福建这种多市县拼合数据的场景里裁剪方式选错后面做坡度、汇水面积、三维建模都会带出连锁误差。3.1 ArcMap里“裁剪”和“掩码提取”到底有什么区别一直有不少人在问arcmap中依靠面图层裁剪dem栅格tif文件和依靠面图层掩码提取到底有啥区别。这两个工具在界面里都能实现“按行政区范围取DEM”但机制完全不同。“裁剪”对应ArcToolbox中的DataManagementTools-栅格-栅格处理-裁剪输入面图层时默认使用面的包络矩形作为裁剪范围输出像元大小、坐标系和原始DEM保持一致。如果你不勾选“使用输入要素裁剪几何”得到的是矩形范围内的完整栅格面边界外依然保留原始值只有勾选该选项才会沿着面边界把外部栅格去掉。这个操作的结果像用刀切蛋糕切出来的块仍保持原来蛋糕的分层结构只是尺寸变了。“掩码提取”对应Spatial Analyst Tools-提取分析-按掩膜提取它的逻辑是先生成一个与面图层对应的掩膜栅格再逐像元判断“是否在掩膜内”输出栅格以掩膜栅格为准。这意味着它会根据环境设置调整像元大小和像元对齐如果掩膜栅格和DEM像元大小不一致就会触发重采样。这样做的代价是边缘更贴合面边界代价是可能改变每个像元的位置相当于每个格子的坐标都按照掩膜网格重新落了一次。这是我常用的选择规则只做显示和出图用裁剪并勾选“使用输入要素裁剪几何”要做地形分析、坡度、汇水面积计算尽量用按掩膜提取并在环境设置里把“捕捉栅格”设为原始DEM避免重采样。很多翻车现场都是在裁剪后直接把结果做坡度结果边缘出现一圈异常条带就是因为像元被悄悄对齐到掩膜栅格。简单说“裁剪”尽可能保留原始网格“掩码提取”则更强调输出网格与掩膜对齐二者没有绝对优劣只有适不适合当前分析。3.2 用GDAL把分块TIF拼成整幅福建影像官方分发的福建DEM往往是一块一块的TIF按标准图幅或行列号命名。我不建议用ArcGIS一个接一个地裁剪那样效率低且容易留下接缝。常见做法是先用GDAL的虚拟栅格技术做拼接索引再统一输出。gdalbuildvrt fujian_mosaic.vrt data/*.tif gdal_translate -co COMPRESSDEFLATE -co TILEDYES -co BIGTIFFIF_SAFER fujian_mosaic.vrt fujian_mosaic.tif第一行gdalbuildvrt不复制数据只生成一个指向所有原始TIF的虚拟文件速度极快第二行gdal_translate才真正把虚拟栅格合并成一块完整TIF。参数-co COMPRESSDEFLATE是使用无损压缩能显著减小文件体积-co TILEDYES让栅格以256×256瓦片存储读取速度比按行存储快很多-co BIGTIFFIF_SAFER让文件超过4GB时自动转成BigTIFF避免写到一半空间不足。如果后续还要重投影直接把第一行的输出换成gdalwarp会更省事因为gdalwarp可以一边拼接一边重投影。需要提醒的是gdalbuildvrt拼接时并不检查各块TIF之间的重叠也不自动做色彩或高程接缝处理。如果相邻两块获取时间不同拼出来可能出现“台阶”。此时要考虑是否把重叠区域平均掉而不是直接拼接。拼接完成后我通常会用gdalinfo再确认一次大小、坐标和NoData确保所有分块的属性一致。若发现某一块坐标系不同要先单独重投影否则拼接出来的文件会在那一块区域“位移”这种错误从直方图上很难发现只有放大到街区尺度才能看到。3.3 重采样方法与分辨率参数30米、12.5米怎么设当目标数据来自不同源但要做统一分析时必须把分辨率统一到一个基准。比如手头有30米SRTM和12.5米ALOS要把它们都转成30米网格参与计算就不能简单地“导出栅格”而是要明确指定重采样算法和目标像元大小。gdalwarp -t_srs EPSG:4547 -tr 30 30 -r bilinear -dstnodata -9999 -overwrite fujian_mosaic.vrt fujian_4547_30m.tif-t_srs指定输出坐标系为CGCS2000 3度带117E-tr 30 30表示输出像元为30米×30米-r bilinear表示双线性插值适合连续表面高程数据-dstnodata -9999统一写入无效值。如果是从12.5米重采样到30米相当于把4个像元合并成1个信息必然有损失“分辨率越高越好”在统一分析场景下反而是障碍因为不同分辨率混用的结果无法对比。另一个常被忽略的参数是内存。全省范围的12.5米TIF可能超过2GBgdalwarp默认内存配置不够时会反复读写临时文件速度极慢。我通常会在命令前加一句设置缓存比如set GDAL_CACHEMAX2048Windows或export GDAL_CACHEMAX2048Linux把缓存提到2GB。实际执行时如果发现进度条长时间停在同一个百分比先看磁盘剩余空间再考虑临时文件目录是否够用。重采样算法也要根据用途选择做坡度、汇水优先双线性做三维展示立方卷积可能更好看最邻近只适合土地分类这类离散栅格用到DEM上会产生明显台阶。4. 让高程TIF“显形”等高线、山体阴影与Blender三维模型处理完几何和分辨率DEM才能被真正“用”起来。对于不熟悉三维的人来说最快的验证方式是把高程转成等高线或山体阴影对于做可视化和汇报的人来说则希望直接进入Blender里建立3D模型。下面把两条路径都走通也会提到“DSM生成DEM”这个经常被混淆的概念。4.1 ArcGIS从DEM提取等高线shp等值线生成参数ArcGIS里从DEM提取等高线shp是标准操作工具位于Spatial Analyst Tools-表面分析-等值线。输入栅格选处理好的TIF输出要素类选目标shp最重要的是“等值线间距”和“起始等高线”两个参数。这个工具对应的就是“arcgis从dem的提分文件提取shp”这个经典需求只是“提分文件”实际说的是“提取等值线”。对于福建这种地形我一般给山区设20米间距沿海平原设5米间距否则山区等高线会密到完全糊成一片。起始等高线通常设为0米或5米这样输出的shp高程属性按整数值递增方便后续标注。Z因子用于单位换算如果DEM高程单位是米Z因子保持1如果是从英尺数据转来的要改成0.3048。这里有一个细节等高线工具只处理有效值NoData区域的边缘不会自动闭合而是留下断头。所以提取等高线前务必已经统一设置NoData。输出shp之后可以用ArcGIS的标注引擎给等高线加高程注记也可以直接导出到AutoCAD做土方图。如果后续要用于Blender建模等高线shp倒是可以作为挤出模型的轮廓线但精度远不如原始栅格。在项目交付时甲方往往只要一份等高线shp而不是原始TIF所以这一步也是“把DEM变成可交付成果”的常用出口。4.2 DSM生成DEM滤波去地物与保持地形的关系热词里“dsm生成dem”是很多人搜的出处。DSM记录的是地物表面高度包括建筑和树冠DEM则是剔除地物后的裸地表。在福建省内闽江口、福州和厦门周边的高分辨率DSM会看到明显的建筑屋顶轮廓如果直接拿去做洪水淹没分析水面将沿着屋顶走结果失真。从DSM得到DEM常见做法是点云分类和地面点插值但如果你拿到的只有一张TIF没有原始点云只能用图像处理方法近似。一种做法是使用低通滤波或开运算把树冠和建筑的“尖刺”抹平。我在QGIS里常用“焦点统计”取邻域最小值配合地形曲率修正能粗略还原裸地表代价是会损失部分真实微地形。顺带说一句如果要做高精度坡度和汇水分析不要在DSM上做宁可重新处理点云。这是很多项目踩过坑后才换回来的结论因为树冠和屋顶会让坡面角度虚高排水路径也会偏。4.3 Blender里用DEM建立3D模型灰度图位移与点云网格“blender 行政高程数据建立3d模型”通常指把某个县域或市域的DEM叠加行政区边界生成一个能旋转看的三维地形。这里我给两条精度等级的路线也顺带回应“别再到处找DEM数据了5米精度高程怎么用”的问题——下载只是开始让数据进入三维才是重点。如果只是快速展示用灰度图位移最快。很多教程会先让你用Global Mapper这类工具把TIF换成PNG但我更习惯直接用GDALgdal_translate -ot Byte -scale 0 1800 0 255 -of PNG fujian_dem.tif fujian_dem_gray.png这条命令把0-1800米的海拔线性映射到0-255灰度-ot Byte强制输出8位-scale是映射范围-of PNG指定格式。然后在Blender里给一个平面细分后加Displacement修改器用这张灰度图作为高度场就能得到一个立体的地形。这个方法的基本思想是灰度越亮对应海拔越高黑的地方是海平面。注意它无法保存真实坐标只能用于展示。如果需要精确尺寸和坐标要把DEM转成点云再生成网格。常见做法是使用gdal2xyz将栅格导出为XYZ坐标文本gdal2xyz -csv fujian_dem_4547.tif fujian_dem.xyz生成的文件每一行是X,Y,Z三个坐标其中X和Y是投影坐标单位是米Z是海拔。把它导入CloudCompare或Blender先做网格化再按行政区边界裁剪。BlenderGIS这类插件可以加载GeoTIFF直接生成带地理参考的网格但要注意插件对坐标系的支持并不完整遇到非WGS84的CGCS2000数据时容易错位我通常还是先用GDAL转成WGS84 UTM再交给插件。行政边界数据可以在裁剪DEM时用也可以在网格生成后用来切边前者效率更高因为参与三维化的点少得多。5. 福建DEM TIF高频踩坑记录范围飞了、黑边、负值、锯齿与超大文件这一章是血泪经验汇总。每一条都是我在项目中真实遇到的现象按“现象→原因→解决”的顺序记录方便你对照排查。能在本地一次跑通最好跑不通时这五条基本覆盖了80%的原始TIF翻车原因。5.1 范围对不上——栅格“飞了”是坐标系没定义现象加载DEM TIF后它和道路、行政边界不能叠加整体偏出几十米到几百米甚至在QGIS里出现在海洋中心。我见过有人把福建沿海的DEM加载进来位置却落在台湾海峡中央第一反应是“文件坏了”其实只是缺少投影定义。原因最常见是TIF缺少投影定义软件按默认WGS84地理坐标渲染或者是坐标系定义本身有误比如明明是CGCS2000 3度带元数据却写成WGS84 UTM。还有一种情况是数据来自两个不同分带比如闽西用117E闽东用120E合并时没有统一整体错位。解决先用gdalinfo查Coordinate System段落。如果显示Unknown用gdal_translate -a_srs EPSG:4547给文件补投影定义如果已知原始是WGS84 UTM要先用gdalwarp转到CGCS2000再参与分析。注意补定义和重投影是两回事前者只是告诉软件“你住在哪”后者才会真正改变坐标数值。补定义后用ArcGIS的“转到XY”工具找一个已知地物核对确认差了多远再决定是否重投影。5.2 裁剪结果黑边——NoData和背景值没处理现象用行政区面裁剪后边界出现一圈黑色锯齿山体顶部也有黑色空洞。在ArcMap里黑色区域不是“没有数据”而是被当作0值绘制成了黑色背景。原因裁剪输出保留了背景值0或255这些值被当成高程值绘成黑色或NoData值在裁剪时被替换成了默认0。更多时候是原始TIF的背景区域就没有设置NoData软件只能按0处理。沿海地区尤其容易出问题因为海平面的真实高程就是0如果0被当成NoData整个沿海平原会被挖空。解决在执行裁剪前用gdal_translate把NoData统一设置为-9999在ArcGIS按掩膜提取时把环境设置的“将NoData设置为”留空如果黑边来自0值用栅格计算器把值为0的像元赋为NoData。做完之后用“识别”工具点几个像元确认边缘值域在正常海拔范围。这个检查十秒就能完成但能防止把黑边数据交给下游。5.3 高程出现负值或异常大值现象对DEM做属性表统计时最小值是-9999最大值是3.4e38做成色带时整幅图一半黑一半白山体形态完全看不出来。原因NoData未参与排除或被当成真实海拔参与统计。另一个常见原因是像素深度混乱Int16编码的厘米级DEM被当成米级读取福建最高峰在文件里会变成几十万米。还有种情况更隐蔽文件里同时存在几种无效值比如-9999和-32768都出现但只排除了其中一种。解决读取前先看gdalinfo里的Type字段Float32直接读Int16需要除以100如果原数据单位是厘米。统计和显示时打开“忽略NoData”在ArcGIS的图层属性里将NoData前打勾。关于DEM处理我一直觉得“像素深度”是最像玄学的一环不同软件对NoData的默认处理不同同一份文件在QGIS和ArcMap中统计结果可能差很多所以标准做法是先把NoData统一再进软件不要依赖某个软件的默认行为。5.4 拼接处有明显接缝、台阶现象将相邻多块TIF拼成全省后在块边界出现明显的海拔突变明明地形是连续的拼起来却像台阶。在山区这种接缝会被山体阴影工具放大成一条亮线。原因不同TIF获取时间不同或来自不同卫星/传感器同一位置的高程存在系统差gdalbuildvrt只是简单按照覆盖顺序取值并不做重叠区融合。即便同一数据源相邻条带间的照射角度差异也会导致细微的高程不一致。解决优先使用同一传感器、同一时期的相邻TIF可以在元数据里查获取时间。若要处理已存在的接缝把重叠区做平均值处理在QGIS的“构建虚拟栅格”里选择平均重采样或使用gdalbuildvrt -r average参数。对于高程数据不要用颜色平衡方案硬修因为那会改变真实海拔值。如果必须精确需要布设地面控制点做整体平差这是一条很重的路一般项目不值得投入。先确认接缝两侧差值有多大如果普遍小于一个像元分辨率很多分析场景可以直接接受。5.5 TIF文件巨大、打开转圈甚至崩溃现象全省12.5米DEM文件超过2GB拖进ArcMap要转十几分钟缩放一次卡三秒最后甚至内存不足退出。有教程会说“别再到处找DEM数据了数据很小一下就加载完”但实际拿到手根本不是那回事因为原始TIF没有经过瓦片化和压缩优化。原因原始TIF大多是按行存储未压缩没有内嵌金字塔也没有瓦片化查看大范围时软件必须读取整块数据。分辨率越高问题越严重12.5米数据的像元数量是30米数据的近6倍文件体积自然成倍上涨。解决用第3章的压缩参数重新整理文件再建金字塔gdaladdo -r average fujian_mosaic.tif 2 4 8 16 32gdaladdo为文件生成2、4、8、16、32倍缩放的概视图-r average表示用平均值采样。以后在ArcGIS或QGIS里缩放软件只用当前比例尺对应的概视图速度会快一个数量级。建金字塔后文件会变大一些但绝对的载入速度收益远大于空间成本。这一步通常能解决“数据量不大但电脑扛不动”的错觉。如果文件确实大到无法加载再考虑先按地市分幅处理用VRT虚拟拼接按需读取而不是一次性导出全量。6. 验证你的DEM TIF剖面检查与批量高程提取两个落地技巧数据处理完成后验证比处理更重要。我基本会做两件事一是用剖面工具看地形是否符合常识二是用Python脚本批量提取点位高程和现场测量值对比。这两步通过数据才算真的可用。6.1 用QGIS剖面工具检查数据是否在“正确的地球”上在QGIS中加载处理后的TIF用“地形分析”工具栏里的剖面工具沿一条已知河流或山脊画一条剖面线。合格的DEM剖面应该平滑连续河流处海拔逐步下降山脊处起伏自然如果剖面出现锯齿状跳动多半是重采样算法不当或NoData混合进了有效值。在ArcGIS里可以使用3D Analyst的剖面工具效果一样。这个检查并不复杂但能拦住后面一连串错误。6.2 用Python批量读取点坐标对应的高程值站点选址时经常手头有一批候选点坐标需要快速拿到每个点的高程。常见做法是用rasterio按点的行列号取值import rasterio with rasterio.open(fujian_dem_4547.tif) as src: for lon, lat in [(119.0, 26.0), (118.8, 26.5)]: row, col src.index(lon, lat) elevation src.read(1, window((row, row1), (col, col1)))[0, 0] print(lon, lat, elevation)这段代码的要点是src.index会把地图坐标换算成栅格行列号再用window参数只读取一个像元避免把整个大TIF读进内存。若你的候选点是WGS84经纬度而栅格是CGCS2000投影坐标必须先做坐标转换否则索引会落在错得离谱的位置。rasterio本身不处理重投影我一般用pyproj.Transformer把经纬度转到栅格坐标系再执行上面逻辑。最后的个人习惯是拿到任何一份DEM原始高程数据TIF先复制一份原始文件留底所有裁剪、重投影都在副本上操作。第一次拿福建沿海某县的数据做微观选址时我直接把原始TIF覆盖成了重投影后的版本后来想核对源数据已经没有后悔药只能重新下载。从那以后原始文件永远设为只读。这份“先验证、留备份、再分析”的习惯帮我少走了很多弯路希望帮到你。本文还有配套的精品资源点击获取