要求
参考资料:
名词解释:
臭氧8小时,也称臭氧8小时滑动平均,即以一天中最大的连续8小时臭氧浓度均值作为评价这一天臭氧污染水平的标准。通常发生在午后光照强温度高的时段。90百分位表示精确到小数点后2位数。
package com.en.utils.util; import java.util.Arrays; / * 环境空气质量综合指数计算 */ public class AreaCompositeUtil {
/ * 日均值二级标准 */ private static final double DAILYLEVEL_2_SO2=150; private static final double DAILYLEVEL_2_NO2=80; private static final double DAILYLEVEL_2_PM10=150; private static final double DAILYLEVEL_2_PM25=75; private static final double DAILYLEVEL_2_CO=4; private static final double DAILYLEVEL_2_O3=160;//8 小时均值二级标准 / * 1、计算各污染物统计量的浓度值 * 统计各城市的 SO2、NO2、PM10、PM2.5 的月均浓度,并统计一氧化碳(CO)日均值的第 95 百分位数以及臭氧(O3)日最大 8 小时值的第 90 百分位数。 */ / * * @param SO2MonthlyMean 月均浓度 * @param NO2MonthlyMean 月均浓度 * @param PM10MonthlyMean 月均浓度 * @param PM25MonthlyMean 月均浓度 * @param CODailyMean CO 日值的集合 * @param O3Hour_8 臭氧(O3)日最大 8 小时值集合 */ public static double getPotency(double SO2MonthlyMean, double NO2MonthlyMean, double PM10MonthlyMean, double PM25MonthlyMean, double[] CODailyMean,double[] O3Hour_8){
Arrays.sort(CODailyMean); double Percentile_CO=CODailyMean[new Double(CODailyMean.length/100*1.0*95).intValue()]; double O3Acount_8=0.0; for(Double O3:O3Hour_8){
O3Acount_8+=O3; } O3Acount_8=Double.parseDouble(String.format("%.2f",O3Acount_8/8*1.0)); return getSingleIndex(SO2MonthlyMean,NO2MonthlyMean,PM10MonthlyMean,PM25MonthlyMean,Percentile_CO,O3Acount_8); } / *(2)计算各污染物的单项指数 * Ii=Ci/Si * 式中:Ci ——污染物 i 的浓度值,当 i 为 SO2、NO2、PM10 及 PM2.5 时,Ci 为月均值,当 i 为 CO 和 O3 时,Ci 为(1)中相应的百分位数浓度值; * Si ——污染物 i 的日均值二级标准(当 i 为 O3 时,为 8 小时均值二级标准)。 */ private static double getSingleIndex(double SO2MonthlyMean, double NO2MonthlyMean, double PM10MonthlyMean, double PM25MonthlyMean, double Percentile_CO,double O3Acount_8){
double Ii_SO2=SO2MonthlyMean/DAILYLEVEL_2_SO2; double Ii_NO2=NO2MonthlyMean/DAILYLEVEL_2_NO2; double Ii_PM10=PM10MonthlyMean/DAILYLEVEL_2_PM10; double Ii_PM25=PM25MonthlyMean/DAILYLEVEL_2_PM25; double Ii_CO=Percentile_CO/DAILYLEVEL_2_CO; double Ii_O3=O3Acount_8/DAILYLEVEL_2_O3; return getCompositeIndex(Ii_SO2,Ii_NO2,Ii_PM10,Ii_PM25,Ii_CO,Ii_O3); } / *3、计算空气质量综合指数 Isum *空气质量综合指数的计算需涵盖全部六项污染物 */ private static double getCompositeIndex(double Ii_SO2,double Ii_NO2,double Ii_PM10,double Ii_PM25,double Ii_CO,double Ii_O3){
return Double.parseDouble(String.format("%.2f",(Ii_SO2+Ii_NO2+Ii_PM10+Ii_PM25+Ii_CO+Ii_O3))); } }