package bianjie;

import java.util.List;

import zhuye.Coordinate;

public class jisuanmianjie {
	/**
	 * 计算coordinates集合中所有经纬度点围成的多边形面积
	 * 使用球面多边形面积计算公式，适用于经纬度坐标
	 * 
	 * @return 多边形面积（平方米）
	 */
	public static double calculatePolygonArea() {
		List<Coordinate> coordinates =Coordinate.coordinates;
	    if (coordinates == null || coordinates.size() < 3) {
	        System.err.println("坐标点数量不足，至少需要3个点才能构成多边形");
	        return 0.0;
	    }
	    
	    int pointCount = coordinates.size();
	    double totalArea = 0.0;
	    final double EARTH_RADIUS = 6371000.0; // 地球半径（米）
	    
	    try {
	        // 将度分格式转换为十进制度
	        double[] lats = new double[pointCount];
	        double[] lons = new double[pointCount];
	        
	        for (int i = 0; i < pointCount; i++) {
	            Coordinate coord = coordinates.get(i);
	            lats[i] = parseDMMToDegree(coord.getLatitude(), coord.getLatDirection());
	            lons[i] = parseDMMToDegree(coord.getLongitude(), coord.getLonDirection());
	        }
	        
	        // 使用球面多边形面积公式
	        for (int i = 0; i < pointCount; i++) {
	            int j = (i + 1) % pointCount;
	            
	            double lat1 = Math.toRadians(lats[i]);
	            double lon1 = Math.toRadians(lons[i]);
	            double lat2 = Math.toRadians(lats[j]);
	            double lon2 = Math.toRadians(lons[j]);
	            
	            totalArea += (lon2 - lon1) * (2 + Math.sin(lat1) + Math.sin(lat2));
	        }
	        
	        totalArea = totalArea * EARTH_RADIUS * EARTH_RADIUS / 2.0;
	        
	        // 取绝对值确保面积为正
	        return Math.abs(totalArea);
	        
	    } catch (Exception e) {
	        System.err.println("计算多边形面积时发生错误: " + e.getMessage());
	        e.printStackTrace();
	        return 0.0;
	    }
	}

	/**
	 * 将度分格式转换为十进制度
	 * 格式示例：纬度"3956.1234" + 方向"N" -> 39.93539度
	 * 
	 * @param dmm 度分格式字符串
	 * @param direction 方向（N/S/E/W）
	 * @return 十进制度
	 */
	private static double parseDMMToDegree(String dmm, String direction) {
	    if (dmm == null || dmm.isEmpty()) {
	        return 0.0;
	    }
	    
	    try {
	        // 找到小数点的位置
	        int dotIndex = dmm.indexOf('.');
	        if (dotIndex == -1) {
	            dotIndex = dmm.length();
	        }
	        
	        // 确保有足够的位数
	        if (dotIndex < 2) {
	            return 0.0;
	        }
	        
	        // 提取度和分
	        int degrees = Integer.parseInt(dmm.substring(0, dotIndex - 2));
	        double minutes = Double.parseDouble(dmm.substring(dotIndex - 2));
	        
	        // 转换为十进制度
	        double decimalDegrees = degrees + minutes / 60.0;
	        
	        // 根据方向调整符号
	        if ("S".equalsIgnoreCase(direction) || "W".equalsIgnoreCase(direction)) {
	            decimalDegrees = -decimalDegrees;
	        }
	        
	        return decimalDegrees;
	        
	    } catch (NumberFormatException e) {
	        System.err.println("坐标格式转换错误: " + dmm);
	        return 0.0;
	    }
	}

	/**
	 * 计算面积并返回带单位的字符串
	 * 
	 * @return 面积字符串（平方米和亩）
	 */
	public static String getAreaString() {
	    double areaM2 = calculatePolygonArea();
	    double areaMu = areaM2 / 666.67; // 转换为亩（1亩 ≈ 666.67平方米）
	    
	    return String.format("面积: %.2f 平方米 (约 %.2f 亩)", areaM2, areaMu);
	}

	/**
	 * 计算面积并返回带单位的字符串（指定小数位数）
	 * 
	 * @param decimalPlaces 小数位数
	 * @return 面积字符串（平方米和亩）
	 */
	public static String getAreaString(int decimalPlaces) {
	    double areaM2 = calculatePolygonArea();
	    double areaMu = areaM2 / 666.67;
	    
	    String format = "面积: %." + decimalPlaces + "f 平方米 (约 %." + decimalPlaces + "f 亩)";
	    return String.format(format, areaM2, areaMu);
	}
}