做地表温度研究的人一定绕不开ASTER这组数据。它是NASA Terra卫星上的中分辨率热红外传感器L2级地表温度产品AST_08的空间分辨率能到90米在城市热岛、火山热异常、地热资源调查这些场景里比MODIS的1公里数据好用太多。但我的学生第一次从NASA Earthdata下载ASTER L2时几乎每步都卡账号怎么注册、产品怎么搜、几十个HDF文件怎么批量下、下载完用哪个软件打开还有最常见的——手里明明有shp文件却不知道可以直接拖进去圈选研究区。这篇文章就是一条龙保姆级教程把我这些年反复操作沉淀下来的流程、参数和踩坑记录全部放出来给准备入门的遥感同学和刚接项目需要LST数据的人直接抄作业。1. 下载前必须弄懂的ASTER L2地表温度数据1.1 先搞清楚你拿到的到底是什么ASTER是Terra卫星上的高级星载热发射反射辐射计可以同时获取可见光近红外15米分辨率、短波红外30米分辨率和热红外90米分辨率数据。其中L2级地表温度产品在NASA的Earthdata平台上的短名是AST_08全称是ASTER L2 Surface Temperature Product目前常用版本是V003。这个产品直接帮你做好了大气校正和温度反演不用自己处理大气参数下载下来就能看到地表热红外亮温反演后的物理温度对做城市热岛、地热异常、火山活动、干旱胁迫甚至考古探测的人都很实用。为什么很多人不用MODIS地表温度而不是ASTER关键就是空间分辨率。MODIS的地表温度产品分辨率大约是1公里适合做省域、全国甚至全球尺度ASTER L2能到90米虽然重访周期差不多16天但在城市街区、农田地块、矿区这种中小时空尺度研究里1公里往往看不出细节90米能看到道路、建筑、水体之间的温度差异。所以如果你要把地表温度结果和土地利用分类、地块边界叠加ASTER是比MODIS合适得多的选择。还有一个容易弄混的地方ASTER的L2产品不止温度一个比如AST_05是地表发射率、AST_07XT是地表反射率、AST_09T是地表辐亮度。标题里说的地表温度一定要认准AST_08不然搜索出来一堆AST_08之外的HDF文件处理时科任何意义。1.2 核心参数和格式提前避免选错ASTER L2地表温度产品的几个关键参数我习惯在下单前先理清楚卫星与传感器Terra卫星上的ASTER传感器。产品短名AST_08不要选成AST_05或AST_09T。产品版本V003。空间分辨率热红外波段90米。单景幅宽大约60千米×60千米。重访周期约16天。数据格式HDF-EOSHDF4变体不是普通GeoTIFF。单位产品里存的表面动力学温度以开尔文K为单位通常还带有比例因子。文件大小单景HDF从几十兆到一百多兆不等视内容而定。下载之前还要明确一个概念Earthdata Search里显示的是一景一景的数据产品每一景对应一条轨道上的某个60千米宽的条带不是整块无缝的大区域。如果你的研究区是100千米×100千米很可能会覆盖两三景影像需要分别下载再做拼接。而且ASTER重访周期16天时间上不可能做到每天都有搜索时一定要把时间范围放大到目标季节或月份否则可能一景都搜不到。另外这个产品在HDF文件内部通常不是一个直接可以显示的栅格而是包含多个科学数据集SDS其中核心数据一般是Surface Kinetic Temperature数组用处是温度值同时还有纬度、经度等辅助定位数组和质量控制信息。第一次打开如果看到一堆子数据集不要慌把它理解为一个压缩包里面装了好几个图层你只要把温度图层取出来用就行。1.3 需要提前准备好的账号和软件下载ASTER L2必须有一个NASA Earthdata账号。注册地址是urs.earthdata.nasa.gov用邮箱注册并激活后在Earthdata Search右上角登录一次后续所有相关服务基本都能通用。这一步千万别偷懒不登录能搜索但点下载时会让你一遍遍跳转登录页下载大文件更会因会话过期而失败。软件方面我建议准备三类下载工具基础的浏览器下载就能用文件多的话Earthdata Search会生成Python或shell下载脚本所以电脑上最好装一个Python 3环境。数据查看与预处理工具ENVI、QGIS、ArcMap/ArcGIS Pro、HDFView四选二即可。ENVI打开HDF-EOS很顺手ArcMap适合后续用shp裁剪和制图HDFView免费小巧适合快速检查波段名。命令行工具Windows可以装Git Bash或者直接用PowerShell主要用来跑下载脚本和写简单的Python读取代码。如果你打算用shp文件做研究区筛选还要事先确认shp能不能用。所谓能不能用最基本的是 .shp、.shx、.dbf、.prj 这几个文件要齐全起码要有前三个。网上很多人只下载了一个.shp文件拖进软件后就报错实际上shp是一组文件不是单文件格式。如果有.prj文件坐标系信息就完整后面拖进Earthdata Search或ArcMap时能少踩很多坑。2. 从NASA Earthdata检索并下载ASTER L2保姆级实操2.1 登录账号并打开Earthdata Search打开设备上的浏览器进入 search.earthdata.nasa.gov。首次使用点击右上角Log In输入Earthdata账号密码登录。这一步看着简单但我见过不少人卡在登录后还是显示未登录的情况主要原因一般是浏览器拦截了第三方Cookie或者登录页跳转没完成。处理办法是关掉广告拦截插件刷新一次页面再去登录。登录后顶部搜索框输入AST_08下拉结果里会出现“ASTER L2 Surface Temperature Product V003”点击它加入收藏也可以直接在左侧数据集中勾选它。页面上接下来有两样东西要设置时间和空间范围。时间范围在左侧的时间轴上选可以用下拉日历精确到天也可以直接输入起止时间。空间范围一般有两种方式手动画框或者导入shp。手动画框适合粗范围找数据而shp适合精确到行政边界、流域边界或地块边界。还有一个细节值得记下来搜索AST_08后地球上的区块会变成有颜色的多边形或方框表示这里有对应的数据场景。点击后会出现该场景的日期、轨道号、云量、浏览图等信息。云量这一项如果看不到也不要紧重点先看浏览图挑没有大片云遮挡的景。2.2 用shp文件一键圈定研究区这个拖拽技巧太重要很多教程只教你在地图上手动画框但如果你手里有研究区的行政边界或流域边界shp手动画框又费劲又不精确还可能把边界外不需要的场景全选进来。这里我强烈推荐直接用shp文件圈定研究区。在Earthdata Search的地图界面左侧工具里通常有一个导入空间范围的功能支持GeoJSON、KML和Shapefile。请牢记散装的shp文件单独的.shp拖上去浏览器经常不认最稳妥的是把.shp、.shx、.dbf、.prj打包成一个zip压缩包然后把这个zip文件直接拖拽到地图区域或者通过导入按钮选择zip包。上传成功后地图上会出现你的研究区边界数据结果也会自动限制在这个边界内搜索结果直接从几百景变成和你研究区真正相交的十几景甚至几景。如果上传后提示无法识别先检查压缩包里的文件是否齐全再看坐标系是否是WGS84或可被Earthdata识别的坐标系。经纬度坐标的shpGCS_WGS_1984是最不容易出问题的如果是CGCS2000或者高斯投影的平面坐标shp建议先在ArcMap里做一次转换再上传不然可能出现边界跑到海里或者完全对不上的情况。另外如果shp是多边形用它圈AOI没有问题但如果是道路shp这种线状要素我更建议你先用ArcMap把道路做缓冲区分析生成面状缓冲后再上传因为Earthdata的空间筛选是以面为基础的。2.3 从结果里挑景并加入Project设置好时间和shp范围后左侧结果列表会列出所有匹配的ASTER L2场景。每一条通常有日期和时间可以先把需要的时间段内所有景都选中不用急着逐景精确到轨道号。ASTER数据一景大约只有60千米宽如果你的研究区比较大跨了几条轨道结果里会出现同一日期、不同轨道号的多个数据块这些都是正常现象全部选上就行。把需要的场景加入右侧的Project后不要急着下载先在Project面板里做一次检查。检查什么第一数量是否合理。研究区100千米×100千米、时间跨度一个月的话理论上十几景到几十景都正常如果只有一景且面积感觉明显不够可能是shp范围没生效回头再看一眼边界。第二时间是否覆盖齐全不要漏掉关键日期的数据。第三有条件的话打开浏览图特别是夏季研究区如果云量覆盖太多温度产品里的有效像元会很少下回来也基本不能用。确认无误后点击Project面板里的Download链接。这里下载方式选择是关键一环按下面这个思路走能省下无数时长。2.4 批量下载直接下载、脚本下载和断点续传下载方式有三种可以按数据量选择只有一两景直接点Download浏览器下载单个HDF文件最简单。有几十景在Download窗口里选择Download Script建议生成Python脚本。脚本会保存成一个download.py文件打开命令行进入该目录运行python download.py只要网络稳定它会逐个下载。如果一个文件中途失败可以删掉这个不完整的文件再跑一遍脚本一般能接着下载。想用wget或curl手动下载也可以在Download窗口拿到链接列表。把这些链接保存成txt一行一个然后执行wget -i links.txt -c。这里的-c就是断点续传遇到断网或服务器中断不用从头开始。我特别想提一个在实际操作里反复出现的坑用脚本下载时经常会有某个请求跳到Earthdata登录页面而不是直接返回文件导致下载下来一个HTML错误页而不是HDF文件。这种情况大概率是认证状态过期。最简单的解决办法是把download.py重新生成一次再用新脚本下载如果是在命令行里用curl或wget也可以提前准备一个.netrc文件放在当前目录里面写清楚Earthdata账号密码让认证自动完成。千万别手动去浏览器复制下载链接再粘贴到软件里这种下载链接通常具有时效性过一段时间就会失效。另外不建议同时跑十几个下载线程。Earthdata虽然不限速但瞬时大量并发请求容易触发服务端限流让某些文件返回403。我实测下来脚本默认的串行或最多3-5个并行就足够了下载速度不慢而且稳定性好很多。如果你的网络本身就不太稳定可以每天晚上挂着脚本第二天检查一下断点续传结果比白天实时盯进度舒服得多。2.5 下载之后先做一轮完整性检查HDF文件下载完成后先别急着解压或转格式先看一眼文件大小。正常的ASTER L2 HDF文件通常有几十兆以上如果看到0字节或只有几KB基本可以判断下载中断或认证失败这种文件直接删掉重新下载。然后再用HDFView打开一个文件看一下能不能正常列出子数据集能列出来说明文件没有损坏。我自己习惯在本地按日期建文件夹比如ASTER_LST/20240526/把所有同一天下的文件放一起文件名里的时间字段可以辅助识别。这样到了后面按时间序列分析时找数据一目了然。不要一股脑全堆在Download目录里几十个HDF文件堆在一起文件名又长后面处理时光核对文件就够头疼的。3. 拿到HDF文件后怎么打开和处理3.1 用HDFView、ENVI和Python读取温度子数据集ASTER L2的HDF文件不是普通图片格式直接双击可能默认打开失败。先推荐免费工具HDFView打开后左侧是一棵树状结构找到Surface Kinetic Temperature这个科学数据集右键选择打开图像就能看到灰度温度影像。如果想在遥感软件里进一步做投影转换和裁剪用ENVI更顺手File → Open External File → EOS → HDF选择产品中对应的SDS即可把它当作普通栅格加载。ArcMap用户要稍微注意老版本ArcMap对HDF4的支持并不算好有时只能看到是栅格目录却无法正确拉伸显示。我通常的解决办法是先用Python的rasterio或GDAL把HDF转成GeoTIFF再拖进ArcMap这样后面按掩膜提取和shp裁剪都会省事很多。用Python打开也很简单可以参考下面这段代码import rasterio hdf_path AST_L2_xxx.hdf with rasterio.open(hdf_path) as src: subdatasets src.subdatasets print(subdatasets) # 第一个子数据集不一定是温度先打印清单确认如果要用GDAL直接定位温度子集并导出GeoTIFF可以这样from osgeo import gdal hdf_path AST_L2_xxx.hdf ds gdal.Open(hdf_path) subs ds.GetSubDatasets() for sub in subs: print(sub[1]) # 查看子数据集名称我的习惯是先打印所有子数据集名称找到包含“Kinetic Temperature”或“Temperature”的项再把它单独打开。提取温度数组之后记得检查属性里的单位ASTER L2的地表温度一般以开尔文为单位而我们要用摄氏度时需要减去273.15。这一步看起来基础但我见过很多新手直接拿K值的灰度图去和其他摄氏温度数据比较最后结果差了300度还浑然不觉。3.2 温度单位转换、NoData处理和投影当把温度数组读到软件或脚本里后通常需要做几步处理。第一是看比例因子。如果产品节点写明了scale_factor为0.1或0.01那么原始DN值需要先乘上这个系数才能得到真实的开尔文温度。第二是无效值。云覆盖区、水体边界和扫描条带边缘可能有一些填充值比如65535或0如果不设置NoData这些非温度值会被当成有效像元参与计算轻则温度均值偏高重则整个温度场出现条带状异常。投影方面ASTER L2原始HDF一般自带地理定位数组但未必是ArcMap默认能正确识别的标准GeoTIFF投影。为了能和shp文件准确叠加建议把温度数据输出为WGS84或UTM投影的GeoTIFF。用GDAL转换时可以指定输出参考系gdalwarp -t_srs EPSG:4326 -dstnodata 65535 input.hdf temp_wgs84.tif这里EPSG:4326就是常见的WGS84经纬度坐标系。如果你的研究区范围不大后续要和矢量边界叠加我更推荐转成UTM投影因为长度和面积误差更小。投影选择的原则其实很简单你的shp是什么坐标系或者你的研究区所在的通用坐标系是什么就把温度栅格转成一致的坐标系省得后面经常做投影变换。3.3 用shp对温度栅格做裁剪按掩膜提取数据到了这一步通常就要把温度结果限制在shp范围内了。在ArcMap里最直接的工具是“按掩膜提取”Extract by Mask位置在“ArcToolbox Spatial Analyst Tools Extraction Extract by Mask”。输入栅格选刚才导出的GeoTIFF温度数据输入掩膜数据选你的研究区shp输出路径给一个GeoTIFF文件点击确定即可。这里要再次强调坐标系一致性。如果栅格是WGS84经纬度shp是CGCS2000高斯投影直接做掩膜提取软件虽然不会报错但提取出来的范围可能偏出你预期的研究区甚至出现黑边或空洞。我的习惯是在做掩膜提取之前先用“分析工具 投影和变换 要素”把shp投影到和栅格一致的坐标系然后再执行提取。这一步多花十秒钟但能避免后面检查成果时怀疑人生的十几分钟。如果你用的是QGIS也可以用“栅格 提取 按掩膜层裁剪栅格”处理逻辑和ArcMap一致。命令行人群还可以用GDAL的裁剪功能gdalwarp -cutline study_area.shp -crop_to_cutline -dstnodata 65535 input_wgs84.tif output_clip.tif-crop_to_cutline参数会让输出范围严格贴合shp边界这对后续做统计非常有帮助。3.4 温度产品的质量意识和多景拼接问题ASTER L2温度数据虽然经过了官方反演但不代表每个像元都能放心使用。云覆盖区域、地形阴影区、未正确配准的扫边区域温度值可能偏差较大。我在做城市热岛研究时会同时下载对应日期的浏览图和QA信息或者在ArcMap里把温度图层调成50%透明度叠加在真彩色影像上看凡是看起来像云的地方都要在后续统计里剔掉。这个步骤可能让有效面积缩减不少但换来的结果是可靠的。如果研究区需要多景拼接建议先对每一景单独做温度定标、NoData设置和投影转换再用“镶嵌至新栅格”Mosaic to New Raster工具拼接。拼接时选择“MEAN”或“BLEND”方法都可以别选“FIRST”因为重叠区域可能出现明显接缝。另外不同日期的场景尽量不要强行拼接成一张温度图地表温度和太阳辐射、云覆盖状态高度相关跨月跨季拼接的成果在物理解释上会非常尴尬。4. 把Excel经纬度数据变成shp文件另一种拖拽玩法4.1 Excel表里的经纬度怎么才能变成点shp你在前面下载ASTER或者用shp做区域选择时手里可能只有一堆研究对象点的坐标比如气象站点、采样点、观测井坐标记录在Excel里。这些坐标能不能直接变成shp拖进ArcMap当然能而且流程很简单。先看Excel表的整理规则。通常至少需要三列点号、经度、纬度。经度是X纬度是Y。如果坐标是度分秒格式一定要先转换成十进制度比如107度30分要写成107.5度否则导入后点的位置会完全不对。字段名建议用英文字母比如ID、Lon、Lat避免ArcMap老版本对中文表头或中文路径出现编码问题。表准备好后把Excel另存为.xls格式或导出成.csv。这里有个经验如果你用的是CSV文件在ArcMap里通过“添加数据”导入时中文字段名和中文内容经常乱码另存为.xls后用“添加XY数据”功能导入成功率更高。我通常就是用.xls省心。4.2 在ArcMap中添加XY数据并导出成shp打开ArcMap执行以下操作菜单栏“文件 添加数据 添加XY数据”。在弹出窗口里选择刚刚那个Excel表。X字段选经度列Y字段选纬度列。点开坐标系按钮设置坐标系如果Excel里的经纬度是基于GPS采集通常是WGS84坐标系选“Geographic Coordinate Systems World WGS 1984”。点击确定。地图上会出现一个临时点图层名字类似“xx事件”。右键这个临时点图层选择“数据 导出数据”导出要素类型选“SHP”保存成一个真正的shp文件。很多新手做到第5步就以为完成了其实这里导出的只是一个临时事件层关掉ArcMap或清除地图文档后可能就没了。必须做第6步导出成shp才能作为文件长期保存。导出时生成的.shp、.shx、.dbf、.prj四个文件会出现在同一个文件夹里以后用的时候直接拖拽整个文件组即可。还有一类情况更隐蔽如果Excel里的坐标已经是平面坐标比如带有中央经线的高斯投影坐标或者UTM投影坐标那么添加XY数据时坐标系一定不能选WGS84经纬度而要选对应的投影坐标系比如UTM Zone 49N或CGCS2000 / 3-degree Gauss-Kruger。选错的结果就是点全部跑到一个看起来完全不相关的地方或者点都重叠在一个小角落这时候不要怀疑Excel数据先去检查坐标系设置。4.3 shp拖进ArcMap的正确姿势和常见错误现在你已经有了研究区shp或点shp怎么拖进ArcMap最稳我推荐的方式是先新建一个空白的ArcMap地图文档然后在Windows资源管理器里找到shp文件按住鼠标左键把它拖到ArcMap左侧的内容列表Table of Contents区域松开鼠标。如果拖到地图视图中央有时候会变成添加数据失败或者图层没进内容列表反而不太方便。拖进去之后ArcMap底部或图层上可能会有黄感叹号提示常见的有三类找不到.shx或.dbf文件组不完整。回文件夹确认三个文件都在并且名称完全一致。无法识别坐标系或缺少投影信息shp没有.prj文件。右键图层属性查看范围如果显示的和预期不符用“定义投影”工具补一个坐标系或者重新在源软件里导出一次。图层虽然能显示但位置和底图对不上大概率坐标系不一致。比如shp是CGCS2000底图是WGS84两者单看数字差不少尤其在高斯投影下会偏移几百米甚至几公里。这种情况需要用“投影”工具把shp转换到和底图或栅格相同的坐标系再进行叠加。如果你拿到的shp是网上下载的道路shp、行政区shp也要先做同样的检查。很多平台分享的shp文件坐标系各不相同有的没有.prj有的是Web墨卡托EPSG:3857有的则是CGCS2000。这些数据本身没有对错只是使用时必须有明确的坐标基准。我的建议是进入项目后给自己定一个统一坐标系比如WGS84 UTM或CGCS2000所有外部shp和遥感栅格都统一到这个坐标系下再分析能省掉后续大量对不齐的烦恼。本地上传。4.4 shp在Earthdata Search和ArcMap之间的联动把Excel转为shp后这个点shp还可以进一步联动到ASTER下载流程里。比如研究区有很多采样点每个点需要对应相应的地表温度值你可以把点shp打包成zip拖进Earthdata Search圈选空间范围筛选出覆盖所有样点的ASTER场景。下载完成后再把点shp拖进ArcMap用“多值提取至点”Extract Multi Values to Points工具把温度栅格值提取到每个点上。这一套流程连起来就是从“坐标表→shp→下载ASTER→抽样提取温度”的完整链路非常实用。如果要用“多值提取至点”注意输入点要素和栅格数据必须完全在同一个坐标系下否则提取结果会出现Null值或错位。我的做法是先把点shp投影到温度GeoTIFF的坐标系再做提取一点差错都不出。5. 常见问题与排查技巧5.1 Earthdata下载阶段的典型问题我把这几年在下载ASTER时遇到的高频问题整理成了一张表方便大家对照排查现象可能原因解决思路搜索结果为空产品没选对时间范围太小shp范围和研究区不符检查搜索框是否有AST_08放大时间范围确认shp边界下载文件是HTML页面认证失效下载链接跳转到登录页重新生成下载脚本或检查.netrc认证配置下载到一半断连网络不稳定或服务器断开使用wget -c断点续传删除不完整文件后重试请求返回403并发下载太多服务端限流降低线程数串行下载或最多3-5个并行HDF文件0字节下载中断或浏览器插件拦截删除该文件重新下载脚本运行报错本机Python环境问题或脚本路径有中文安装Python 3并确认命令可用将文件夹路径改成英文这里面最值得提醒的是“脚本运行报错”。很多人把download.py放在桌面下名为“新建文件夹”的中文路径里Windows命令行解析起来很容易出问题。我一般会专门建一个英文路径比如D:\aster_data\download所有脚本、链接文件和HDF文件都放这个路径下能避免很多莫名其妙的问题。5.2 HDF文件打开和处理阶段的常见问题HDF文件打不开不只是软件兼容性问题很多时候是选错了子数据集。我见过有人用ENVI打开HDF后默认加载了Geolocation数组看见经度和纬度两个数组当成了影像输出后全是经纬度值还以为是下载错了数据。所以一定要认准温度相关的SDS比如Surface Kinetic Temperature、Surface Temperature这类字段。如果温度影像显示出来是一片漆黑通常不是没有数据而是没有做对比度拉伸。HDF原始值范围往往在2000到4000之间DN值视比例因子而定直接按默认的0到255范围显示当然全黑。解决办法是在图层属性里做一次百分比拉伸或者把DN值转换为摄氏温度后再渲染。用ArcMap的朋友可以右键图层属性在符号系统里选择Stretch类型为Percent Clip设置2%到98%的裁剪范围温度梯度一下就出来了。投影和shp叠加对不上的问题多半也和坐标参考有关。这里有一点容易忽视Earthdata上有些栅格产品虽然写着WGS84但实际是原始条带数据可能不带标准投影信息。如果你在ArcMap里看到温度栅格范围显示的是0到几百、0到几百这种奇怪数字最好用GDAL或ENVI重新生成一个带地理定位信息的GeoTIFF再做后续操作不要将错就错。5.3 温度数值对不对一个必须做的快速验证无论下载和处理流程多顺利拿到温度产品后我都会做一个快速验证随便取一个已知地物或站点的温度看看数值是否合理。比如夏季白天城市地表温度经常在30℃到50℃之间夜间可能降到10℃到25℃如果读出来是300多K但没转摄氏度或者转完后数值高达300℃以上就要回头检查定标和单位换算。还有一个验证技巧在ArcMap里用“识别”工具点击水体中心点。水体热惰性大昼夜温差小如果水体像元温度比周边裸地高得离谱说明数据或处理流程可能出了问题。这个检查和下载阶段的质量控制配合起来可以让你在进入正式分析前就把数据问题扼杀在摇篮里而不是等统计结果出来了才发现异常。5.4 收藏一个长期可用的工作流模板最后说点工作流层面的经验。做遥感数据下载和处理最怕每次都临时想起一个命令、临时试一种投影这样效率很低。我自己在电脑里维护一个小文件夹里面有下载脚本模板、HDF转GeoTIFF的Python脚本、shp坐标检查脚本和一份简单的README记录每次处理时的注意点。每次新接ASTER任务就把这个模板复制一份改日期、改研究区直接跑。这样既保证了前后处理标准一致也避免重复踩同一个坑。如果你经常用ArcMap做掩膜提取还可以把“定义投影”“投影工具”“按掩膜提取”“多值提取至点”这些操作存成ArcToolbox里的模型ModelBuilder。下次只要输入新的HDF转好的GeoTIFF和新的shp模型会自动完成投影校验和裁剪提取。我实测下来这套模型用在批量处理几十景ASTER温度数据时能省下一整天的手动劳动。最后顺手再分享两个压箱底的小习惯一个是下载ASTER之前哪怕你很确定坐标系统也尽量把shp转成WGS84并在Earthdata Search里拖拽验证一次。这个动作看着小却能避免你把研究区边界搞错后辛辛苦苦下了一堆不在范围内的数据。另一个是处理完温度数据后先把结果存成一个副本再开始做拉伸、裁剪和拼接。HDF和GeoTIFF这种栅格文件处理不可逆原始数据一旦被覆盖或误调重新下载又得花很多时间。我给数据文件起名的规律是原始HDF加_raw后缀转投影后的GeoTIFF加_wgs84裁剪后的加_clip每一步都不覆盖上一步的结果。这样即使后面某一步参数没调好也能快速回到上一个阶段重新处理不用全部推倒重来。希望这些经验能帮你在ASTER L2地表温度数据的下载和处理上少走几步弯路。