其他
GEE（谷歌地球引擎）已经开始收录Landsat 9数据
Google Earth Engine (GEE)已经开始收录Landsat 9数据，目前总共收录了1773张影像，离USGS之前开放下载的57,000张影像还有一定距离。影像收录还在持续进行中。相信在接下来几天（或者几周）一定会全部收录。大家拭目以待。
目前GEE数据库里已有的Landsat 9数据主要覆盖北美和南美，具体覆盖区域见下面这个图。Landsat 9在GEE的ID如下：
ee.ImageCollection('LANDSAT/LC09/C02/T1_L2')上面这个图由如下代码生成。
https://code.earthengine.google.com/393726e99992467b1afde026e45f99ad
var collection = ee.ImageCollection('LANDSAT/LC09/C02/T1_L2') .sort('CLOUD_COVER'); print(collection.size());var c = collection.median();// Applies scaling factors.function applyScaleFactors(image) { var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2); var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0); return image.addBands(opticalBands, null, true) .addBands(thermalBands, null, true);}
var dataset = applyScaleFactors(c);
var visualization = { bands: ['SR_B4', 'SR_B3', 'SR_B2'], min: 0.0, max: 0.3,};
Map.addLayer(dataset, visualization, 'True Color (432)');Map.centerObject(c);下面这两张图是Landsat 9数据以及主成分分析的结果。（Credits to Gena）
https://code.earthengine.google.com/cbabe960ef629ce679e4327a77b5dc44
var utils = require('users/gena/packages:utils')
var imagesL8 = ee.ImageCollection("LANDSAT/LC08/C02/T1_L2").filterBounds(Map.getCenter()) .filter(ee.Filter.lt('CLOUD_COVER', 5))
var imagesL9 = ee.ImageCollection("LANDSAT/LC09/C02/T1_L2").filterBounds(Map.getCenter()) .filter(ee.Filter.lt('CLOUD_COVER', 5))
var imageL9 = imagesL9.first().unmask().resample('bicubic').multiply(0.0000275).add(-0.2)Map.addLayer(imageL9, { bands: ['SR_B4', 'SR_B3', 'SR_B2'], min: 0.0, max: 0.2, gamma: 1.4 }, 'L9')
var imageL8 = imagesL8.first().unmask().resample('bicubic').multiply(0.0000275).add(-0.2)Map.addLayer(imageL8, { bands: ['SR_B4', 'SR_B3', 'SR_B2'], min: 0.0, max: 0.2, gamma: 1.4 }, 'L8', false)
var bounds = Map.getBounds(true)var scale = Map.getScale()
var imageL9PCA = utils.PCA(imageL9, bounds, scale)Map.addLayer(imageL9PCA.select([1,2,0]).add(imageL9PCA.select([3,4,5])), { min: -2.5, max: 2.5 }, 'PCA L9')
var imageL8PCA = utils.PCA(imageL8, bounds, scale)Map.addLayer(imageL8PCA.select([1,2,0]).add(imageL8PCA.select([3,4,5])), { min: -2.5, max: 2.5 }, 'PCA L8', false)