java空气质量指数AQI算法

我就是我 2022-11-12 09:52 413阅读 0赞

项目中要对空气质量指数AQI进行计算,空气质量指数AQI浓度限值找到的资料如下:

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2xvdmUxMDI2OTk5MDQ1_size_16_color_FFFFFF_t_70][]

计算方法如下:

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2xvdmUxMDI2OTk5MDQ1_size_16_color_FFFFFF_t_70 1][]

如何用java实现这个算法了,下面贴上计算代码:

package com.rate.rate_pz_server.util;
     
     
    import com.rate.rate_pz_server.entity.Aqi;
     
    import java.math.BigDecimal;
    import java.util.ArrayList;
    import java.util.Collections;
    import java.util.Comparator;
    import java.util.List;
     
    /** 
    * @ClassName: AqiUtil 
    * @Description: 计算aqi的值
    * @author jiangya
    * @date 2018年12月14日 上午10:29:35 
    *  
    */
    public class AqiUtil {
    	/**
    	 *  
    	 * 
    	 * @Title: getPollutionDegree 
    	 * @Description: 计算aqi值对应的等级
    	 * @param @param
    	 *            aqi
    	 * @param @return
    	 * @return int   @throws 
    	 */
    	public static int getPollutionDegree(Double aqi) {
    		int pollutionDegree = 1;
    		if (aqi <= 50) {
    			pollutionDegree = 1;
    		} else if (aqi > 50 && aqi <= 100) {
    			pollutionDegree = 2;
    		} else if (aqi > 100 && aqi <= 150) {
    			pollutionDegree = 3;
    		} else if (aqi > 150 && aqi <= 200) {
    			pollutionDegree = 4;
    		} else if (aqi > 200 && aqi <= 300) {
    			pollutionDegree = 5;
    		} else if (aqi > 300) {
    			pollutionDegree = 6;
    		}
    		return pollutionDegree;
    	}
     
    	/**
    	 *  
    	 * 
    	 * @Title: getDegree 
    	 * @Description: 计算aqi值对应的等级
    	 * @param @param
    	 *            pollutionDegree
    	 * @param @return
    	 * @return String   @throws 
    	 */
    	public static String getDegree(int pollutionDegree) {
    		if (pollutionDegree == 1) {
    			return "优";
    		} else if (pollutionDegree == 2) {
    			return "良";
    		} else if (pollutionDegree == 3) {
    			return "轻度污染";
    		} else if (pollutionDegree == 4) {
    			return "中度污染";
    		} else if (pollutionDegree == 5) {
    			return "重度污染";
    		} else if (pollutionDegree == 6) {
    			return "严重污染";
    		}
    		return "数据错误";
    	}
     
    	/**
    	 * 计算每种污染物项目 P的空气质量分指数
    	 * 
    	 * @param cp
    	 *            污染物项目P的质量浓度
    	 * @param r
    	 *            污染物项目P所在数组中的行号
    	 * @return
    	 */
    	public static double countPerIaqi(double cp, int r) {
    		double bph = 0; // 与 cp相近的污染物浓度限值的高位值
    		double bpl = 0; // 与 cp相近的污染物浓度限值的低位值
    		double iaqih = 0; // 与 bph对应的空气质量分指数
    		double iaqil = 0; // 与 bpl对应的空气质量分指数
    		double iaqip = 0; // 当前污染物项目P的空气质量分指数
    		// 空气质量分指数及对应的污染物项目浓度限值
    		int[][] aqiArr = { { 0, 50, 100, 150, 200, 300, 400, 500 }, { 0, 35, 75, 115, 150, 250, 350, 500 },
    				{ 0, 50, 150, 250, 350, 420, 500, 600 }, { 0, 2, 4, 14, 24, 36, 48, 60 },
    				{ 0, 40, 80, 180, 280, 565, 750, 940 }, { 0, 160, 200, 300, 400, 800, 1000, 1200 },
    				{ 0, 50, 150, 475, 800, 1600, 2100, 2620 }, {0,100,160,215,265,800}};
     
    		double min = aqiArr[r][0];
    		int index= aqiArr[r].length-1;
    		double max = aqiArr[r][index];
    		if (cp<=min || cp>=max){
    			return 0;
    		}else {
    			// 对每种污染物的bph、bpl、iaqih、iaqil进行赋值
    			for (int i = r; i < r + 1; i++) {
    				for (int j = 0; j < aqiArr[0].length; j++) {
    					if (cp < aqiArr[i][j]) {
    						bph = aqiArr[i][j];
    						bpl = aqiArr[i][j - 1];
    						iaqih = aqiArr[0][j];
    						iaqil = aqiArr[0][j - 1];
    						break;
    					}
    				}
    			}
    			// 计算污染物项目 P的空气质量分指数
    			iaqip = (iaqih - iaqil) / (bph - bpl) * (cp - bpl) + iaqil;
    			BigDecimal bg = new BigDecimal(Math.ceil(iaqip));
    			double f1 = bg.setScale(0, BigDecimal.ROUND_HALF_UP).doubleValue();
    			return f1;
    		}
    	}
     
    	/**
    	 * 根据提供污染物的各项指标，对AQI进行计算
    	 * 
    	 * @param pmtw
    	 *            PM2.5
    	 * @param pmte
    	 *            PM10
    	 * @param co
    	 *            一氧化碳浓度
    	 * @param no2
    	 *            二氧化氮浓度
    	 * @param o3
    	 *            臭氧浓度
    	 * @param so2
    	 *            二氧化硫浓度
    	 * @return
    	 */
    	public static Aqi CountAqi(double pmtw, double pmte, double co, double no2, double o3, double so2 ) {
    		double pmtwIaqi = getPm25IAQI(pmtw);
    		double pmteIaqi = getPm10IAQI(pmte);
    		double coIaqi = getCoIAQI(co);
    		double no2Iaqi = getNo2IAQI(no2);
    		double o3Iaqi = getO3OneHourIAQI(o3);
    		double so2Iaqi = getSo2IAQI(so2);
    		// 初始化对象数组
    		List<Aqi> aList = new ArrayList<Aqi>();
    		if(pmtwIaqi != 0){
    			aList.add(new Aqi("PM2.5", pmtwIaqi));
    		}
    		if(pmteIaqi != 0){
    			aList.add(new Aqi("PM10", pmteIaqi));
    		}
    		if(coIaqi != 0){
    			aList.add(new Aqi("CO", coIaqi));
    		}
    		if(no2Iaqi != 0){
    			aList.add(new Aqi("NO2", no2Iaqi));
    		}
     
    		if(o3Iaqi != 0){
    			aList.add(new Aqi("O3", o3Iaqi));
    		}
    		if(so2Iaqi != 0){
    			aList.add(new Aqi("SO2", so2Iaqi));
    		}
    		Collections.sort(aList, new AqiComparator());
    		Aqi aqi = null;
    		if(aList.size() > 0){
    			aqi = aList.get(aList.size() - 1);
    		}
    		return aqi;
    	}
    	public static Aqi CountAqi(double pmtw, double pmte, double co, double no2, double o3, double so2,double o3Eight) {
    		double pmtwIaqi = getPm25IAQI(pmtw);
    		double pmteIaqi = getPm10IAQI(pmte);
    		double coIaqi = getCoIAQI(co);
    		double no2Iaqi = getNo2IAQI(no2);
    		double o3Iaqi = getO3OneHourIAQI(o3);
    		double so2Iaqi = getSo2IAQI(so2);
    		double o3EightHourIAQI = getO3EightHourIAQI(o3Eight);
    		// 初始化对象数组
    		List<Aqi> aList = new ArrayList<Aqi>();
    		if(pmtwIaqi != 0){
    			aList.add(new Aqi("PM2.5", pmtwIaqi));
    		}
    		if(pmteIaqi != 0){
    			aList.add(new Aqi("PM10", pmteIaqi));
    		}
    		if(coIaqi != 0){
    			aList.add(new Aqi("CO", coIaqi));
    		}
    		if(no2Iaqi != 0){
    			aList.add(new Aqi("NO2", no2Iaqi));
    		}
     
    		if(o3Iaqi != 0){
    			aList.add(new Aqi("O3", o3Iaqi));
    		}
    		if(o3EightHourIAQI != 0){
    			aList.add(new Aqi("o3EightHour", o3EightHourIAQI));
    		}
    		if(so2Iaqi != 0){
    			aList.add(new Aqi("SO2", so2Iaqi));
    		}
    		Collections.sort(aList, new AqiComparator());
    		Aqi aqi = null;
    		if(aList.size() > 0){
    			aqi = aList.get(aList.size() - 1);
    		}
    		return aqi;
    	}
     
    	public static double getCoIAQI(double co){
    		if(co > 0){
    			return  countPerIaqi(co, 3);
    		}
             return 0;
        }
     
        public static double getNo2IAQI(double no2){
    		if(no2 > 0){
    			long round = Math.round(no2);
    			return countPerIaqi(round,  4);
    		}
            return 0;
        }
     
        public static double  getO3OneHourIAQI(double o3One){
    		if(o3One > 0){
    			long round = Math.round(o3One);
    			return countPerIaqi(round,  5);
    		}
    		return 0;
        }
     
        public static double getO3EightHourIAQI(double o3Eight){
    	    if(o3Eight > 0 && o3Eight <= 800){
    			long round = Math.round(o3Eight);
                return countPerIaqi(round, 7);
            }
    	    return 0;
        }
     
        public static double getPm10IAQI(double pmte){
    		if(pmte > 0){
    			long round = Math.round(pmte);
    			return countPerIaqi(round, 2);
    		}
    		return 0;
        }
     
        public static double getPm25IAQI(double pmtw){
    		if(pmtw > 0){
    			long round = Math.round(pmtw);
    			return countPerIaqi(round, 1);
    		}
    		return 0;
     
        }
     
        public static double getSo2IAQI(double so2){
    		if(so2 > 0){
    			long round = Math.round(so2);
    			return countPerIaqi(round, 6);
    		}
    		return 0;
     
        }
    }
     
    /**
     * 构造分类器类，对AQI对象链表进行排序
     */
    class AqiComparator implements Comparator {
    	public int compare(Object o1, Object o2) {
    		Aqi a1 = (Aqi) o1;
    		Aqi a2 = (Aqi) o2;
    		double result = a1.getAqi() - a2.getAqi();
    		return (int) result;
    	}
    }

实体类如下:

package com.rate.rate_pz_server.entity;
     
    import java.io.Serializable;
     
    /**
     * @ClassName: Aqi
     * @Description: TODO
     * @author jiangya
     * @date 2019年8月14日 下午4:37:25
     *
     */
    public class Aqi implements Serializable {
    	private static final long serialVersionUID = 1L;
    	private String name;
     
    	private Double aqi;
     
    	public Aqi(String name, Double aqi) {
    		super();
    		this.name = name;
    		this.aqi = aqi;
    	}
     
    	public String getName() {
    		return name;
    	}
     
    	public void setName(String name) {
    		this.name = name;
    	}
     
    	public Double getAqi() {
    		return aqi;
    	}
     
    	public void setAqi(Double aqi) {
    		this.aqi = aqi;
    	}
     
    }

原文链接：https://blog.csdn.net/love1026999045/article/details/99582775

# [基于C\#WPF的空气质量预报预警系统设计与实现][C_WPF] #

AQI百度百科：[https://baike.baidu.com/item/AQI/463245?fr=aladdin][https_baike.baidu.com_item_AQI_463245_fr_aladdin]

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2xvdmUxMDI2OTk5MDQ1_size_16_color_FFFFFF_t_70]: /images/20221022/482c17717be149e29eb7b07aed42925c.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2xvdmUxMDI2OTk5MDQ1_size_16_color_FFFFFF_t_70 1]: /images/20221022/ca998e2fe07d4e389cccb43afd81cd3b.png
[C_WPF]: https://www.doc88.com/p-3387324091042.html
[https_baike.baidu.com_item_AQI_463245_fr_aladdin]: https://baike.baidu.com/item/AQI/463245?fr=aladdin