GEE：Landsat-T1级别影像去云和辐射校正

文章目录

• • • Landsat-9
    • • Landsat-9 去云和云阴影
      • Landsat-9 辐射校正
    • Landsat-8
    • • Landsat-8 去云
      • Landsat-8 辐射校正
    • Landsat-7
    • • Landsat-7 去云和云阴影
      • Landsat-7 辐射校正
    • Landsat-5
    • • Landsat-5 去云
      • Landsat-5 辐射校正

这几天在做指数运算的时候，发现landsatT1级别数据是需要做辐射校正的，看到很多人使用LandsatT1数据，所以记录一下T1数据做辐射校正的方法和使用的波段。使用的数据集和波段如下表所示，

Landsat-9

// 未辐射校正数据
var L5 = ee.ImageCollection("LANDSAT/LC09/C02/T1").filterBounds(roi).filterDate('2022-1-1','2022-10-30').filterMetadata('CLOUD_COVER',"less_than",15).map(radiance).max().clip(roi);

Landsat-9 去云和云阴影

function remove_cloud(image){//去云和云阴影var cloudShadowBitMask = (1 << 3);  var cloudsBitMask = (1 << 5);  var qa = image.select("QA_PIXEL");  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0).and(qa.bitwiseAnd(cloudsBitMask).eq(0));  return image.updateMask(mask);  
}

Landsat-9 辐射校正

function radiance(image){var ThermalRadiance=ee.Algorithms.Landsat.calibratedRadiance(image.select('B10')).rename('B10_radiance');return image.addBands(ThermalRadiance)
}

Landsat-8

var L8 = ee.ImageCollection("LANDSAT/LC08/C01/T1").filterBounds(roi).filterDate('2020-1-1','2020-12-31').filterMetadata('CLOUD_COVER',"less_than",15).map(remove_cloud) .map(radiance).max().clip(roi)

Landsat-8 去云

function remove_cloud(image){var image_qa = image.select('BQA');var maskComposite = image_qa.bitwiseAnd(1 << 4).eq(0)return image.updateMask(maskComposite) 
}

Landsat-8 辐射校正

function radiance(image){var ThermalRadiance=ee.Algorithms.Landsat.calibratedRadiance(image.select('B10')).rename('B10_radiance');return image.addBands(ThermalRadiance)
}

Landsat-7

// 未辐射校正数据
var L7 = ee.ImageCollection("LANDSAT/LE07/C02/T1").filterBounds(roi).filterDate('2013-1-1','2013-12-31').filterMetadata('CLOUD_COVER',"less_than",15).map(radiance).max().clip(roi);

Landsat-7 去云和云阴影

function remove_cloud(image){//去云和云阴影var cloudShadowBitMask = (1 << 4);  var cloudsBitMask = (1 << 3);  var qa = image.select("QA_PIXEL");  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0).and(qa.bitwiseAnd(cloudsBitMask).eq(0));  return image.updateMask(mask);  
}

Landsat-7 辐射校正

function radiance(image){var ThermalRadiance=ee.Algorithms.Landsat.calibratedRadiance(image.select('B6_VCID_1')).rename('B6_radiance');return image.addBands(ThermalRadiance)
}

Landsat-5

// 未辐射校正数据
var L5 = ee.ImageCollection("LANDSAT/LT05/C01/T1").filterBounds(roi).filterDate('2005-1-1','2005-12-31').filterMetadata('CLOUD_COVER',"less_than",15).map(radiance).max().clip(roi);

Landsat-5 去云

function remove_cloud(image){var image_qa = image.select('BQA');var maskComposite = image_qa.bitwiseAnd(1 << 4).eq(0)return image.updateMask(maskComposite)
}

Landsat-5 辐射校正

function radiance(image){var ThermalRadiance=ee.Algorithms.Landsat.calibratedRadiance(image.select('B6')).rename('B6_radiance');return image.addBands(ThermalRadiance)
}

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！