Arcgis基础操作
第一章:栅格数据行列号不一致问题
前言
需要跑CASA模型,根据矢量边界从GEE下载相关数据,结果和原有的LULC的行列号不匹配,这会导致模型报错,因此需要先统一行列号。
一、数据下载
用GEE按照矢量边界下载所需的aet、pet、pr、temp、srad、ndvi数据(只想大概跑一下模型看看,所以基础数据全部用的产品),顺便放个GEE月均数据下载代码(terraclimate的数据本就是月均,我这代码本来是下mod11数据的)~
代码如下:
Map.addLayer(table,{},'table')
var dataset = ee.ImageCollection("IDAHO_EPSCOR/TERRACLIMATE")
.filterDate('2019-01-31', '2019-12-31');
for (var year=2019;year<2020;year++){
for (var i=1;i<13;i++){
var startDate = ee.Date.fromYMD(ee.Number(year),ee.Number(i),1);
var endDate = ee.Date.fromYMD(ee.Number(year),ee.Number(i),28);
var AET=dataset.filterDate(startDate,endDate).select('aet').mean();
Export.image.toDrive({
image: AET,
description: 'AET'+year+i,
folder: 'AET',
region: table,
scale: 30,
crs: "EPSG:4326",
maxPixels: 1e13});
}
}
二、统一行列号
下好数据用arcgis处理,步骤如下:
1.栅格转投影(度转米~)→project raster(注:记得也检查土地利用的投影坐标,UTM)2.裁剪(在此之前我还用栅格计算器乘了个系数)→Extract by Mask
注意:裁剪的时候设定环境!!不要default!把裁剪边界选上!snap raster也选裁剪边界!
其实我一般用模型构建器批量处理但也一定一定注意设定环境边界哦
总结
这部分很基础,但本人很菜ಥ_ಥ刚发现行列号不一致时还走了点弯路,所以想分享出来,也当作记个笔记(●ˇ∀ˇ●)