GEE：空气污染指数 (API)计算

作者：CSDN @ _养乐多_

本文记录了基于以下公式计算的空气污染指数代码：

API = 0:7(NO2)+ 0:15(SO2)+ 0:1(O3)+ 0:05(CO) 计算空气污染指数 (API)。

一、完整代码

var roi = table;var dataset_NO2 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_NO2");
var dataset_CO = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_CO");
var dataset_SO2 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2");
var dataset_AER_AI = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_AER_AI');
var dataset_CH4 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_CH4");
var dataset_O3 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_O3");var imageNO2 = dataset_NO2.select('NO2_column_number_density').filterDate('2022-01-01', '2022-12-25').filterBounds(roi).mean().clip(roi.geometry());var band_vizNO2 = {min: 0,max: 0.0002,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};var imageCO = dataset_CO.filterBounds(roi).filterDate('2022-01-01', '2022-12-25').select('CO_column_number_density').mean().clip(roi);var band_vizCO = {min: 0,max: 0.05,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};var imageSO2 = dataset_SO2.filterBounds(roi).filterDate('2022-01-01', '2022-12-25').select('SO2_column_number_density').mean().clip(roi);var band_vizSO2 = {min: 0.0,max: 0.0005,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};var imagePM = dataset_AER_AI.filterBounds(roi).filterDate('2022-01-01', '2022-12-25').select('absorbing_aerosol_index').mean().clip(roi);var band_vizPM = {min: -1,max: 2.0,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};var imageO3 = dataset_O3.filterBounds(roi).filterDate('2022-01-01', '2022-12-25').select('O3_column_number_density').mean().clip(roi);var band_vizO3 = {min: 0.12,max: 0.15,palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};Map.addLayer(imageNO2, band_vizNO2, 'Mean NO2 Bengaluru');
Map.addLayer(imageCO, band_vizCO, 'Mean CO Bengaluru');
Map.addLayer(imageSO2, band_vizSO2, 'Mean SO2 Bengaluru');
Map.addLayer(imagePM, band_vizPM, 'Mean PM2.5 Bengaluru');
Map.addLayer(imageO3, band_vizO3, 'Mean O3 Bengaluru');var image = imageNO2.addBands(imageCO).addBands(imageSO2).addBands(imagePM).addBands(imageO3).addBands(imageO3)function calculateAirPollutionIndex(image) {var index = image.select('NO2_column_number_density').multiply(0.7).add(image.select('SO2_column_number_density').multiply(0.15)).add(image.select('O3_column_number_density').multiply(0.1)).add(image.select('CO_column_number_density').multiply(0.05));return image.addBands(index.rename('air_pollution_index'));
}print(calculateAirPollutionIndex(image))var visParams = {min: 0,max: 1,palette: ['green', 'yellow', 'orange', 'red']
};Map.addLayer(calculateAirPollutionIndex(image).select('air_pollution_index'), visParams, 'Air pollutin Index');Map.centerObject(roi, 11)

二、完整连接

https://code.earthengine.google.com/8e3d58d0bac4f9f1adab5b5e73222a44?noload=true

声明：
本人作为一名作者，非常重视自己的作品和知识产权。在此声明，本人的所有原创文章均受版权法保护，未经本人授权，任何人不得擅自公开发布。
本人的文章已经在一些知名平台进行了付费发布，希望各位读者能够尊重知识产权，不要进行侵权行为。任何未经本人授权而将付费文章免费或者付费（包含商用）发布在互联网上的行为，都将视为侵犯本人的版权，本人保留追究法律责任的权利。
谢谢各位读者对本人文章的关注和支持！