空气质量综合指数JAVA算法

要求

参考资料：

名词解释：
臭氧8小时，也称臭氧8小时滑动平均，即以一天中最大的连续8小时臭氧浓度均值作为评价这一天臭氧污染水平的标准。通常发生在午后光照强温度高的时段。90百分位表示精确到小数点后2位数。

一氧化碳日均值第95百分位数浓度：
比如20个数字组成的数列，先从小到大排列，第95百分数就是第19个数。第50百分数就是中位数。

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)));
    }






}