import csv
import os
import re
import time
from datetime import datetime

from qgis import processing
from PyQt5.QtCore import QVariant, QThread
from PyQt5.QtWidgets import QFileDialog
from qgis._core import QgsFeatureRequest, QgsMapLayer, QgsWkbTypes, QgsCoordinateReferenceSystem, \
    QgsCoordinateTransform, QgsProject, QgsGeometry, QgsGeometryCollection, QgsVectorLayer, QgsFields, QgsField, \
    QgsFeature, Qgis, QgsPointXY, QgsUnitTypes, QgsDistanceArea, QgsProcessingContext, QgsSpatialIndex
from qgis.utils import iface

from qgis.PyQt import QtWidgets, uic, QtCore
from qgis.PyQt.QtCore import pyqtSignal

FORM_CLASS, _ = uic.loadUiType(os.path.join(
    os.path.dirname(__file__), 'cyanlove_calc_intersects_area_base.ui'))


class cyanlove_calc_intersects_area(QtWidgets.QDockWidget, FORM_CLASS):
    closingPlugin = pyqtSignal()

    def __init__(self, parent=None):
        super(cyanlove_calc_intersects_area, self).__init__(parent)

        self.export_thread = None
        self.iface = iface
        self.canvas = iface.mapCanvas()
        self.setupUi(self)
        self.setFloating(True)
        self.property_dict = None  # 存储搜索出来的名称，图元中心经纬度
        self.is_running = False  # 线程状态标志
        self.setWindowFlags(QtCore.Qt.Dialog)  # 使对话框始终在前
        # 设置为模态窗口
        # self.setWindowModality(QtCore.Qt.ApplicationModal)  # 或者使用 Qt.WindowModal
        self.btn_calc_intersects_area.clicked.connect(self.calc_intersects_area)

    def calc_intersects_area(self):
        options = QFileDialog.Options()  # 创建文件对话框选项
        options |= QFileDialog.ReadOnly  # 可以设置一些选项，例如只读
        csv_file_path, _ = QFileDialog.getSaveFileName(self,
                                                       "选择文件",
                                                       "",
                                                       "CSV文件 (*.csv);;所有文件 (*)",  # 可以添加其他过滤器
                                                       options=options)  # 传递 options 参数
        # 检查用户是否选择了文件
        if csv_file_path:
            current_layer1 = self.mMapLayerComboBox.currentLayer()
            current_layer2 = self.mMapLayerComboBox_2.currentLayer()
            if current_layer1 is None or current_layer1.wkbType() not in [QgsWkbTypes.Polygon,
                                                                          QgsWkbTypes.MultiPolygon]:
                self.label_show.setText("选择的图层1不是Polygon矢量图层！")
                return
            if current_layer2 is None or current_layer2.wkbType() not in [QgsWkbTypes.Polygon,
                                                                          QgsWkbTypes.MultiPolygon]:
                self.label_show.setText("选择的图层2不是Polygon矢量图层！")
                return
            # 确保使用现有的进度条
            self.progressBar.setRange(0, 100)  # 设置进度条范围
            self.progressBar.setValue(0)  # 初始化进度条值
            self.progressBar.show()  # 显示进度条 (如果之前是隐藏状态)
            # 启动导出线程
            self.export_thread = ExportThread(current_layer1, current_layer2, csv_file_path)
            self.export_thread.progress.connect(self.update_progress)  # 连接信号到槽
            self.export_thread.updatelabel.connect(self.update_label)
            self.export_thread.finished.connect(self.export_finished)  # 导出完成后连接到槽
            self.export_thread.start()  # 启动线程

    def update_label(self, value):
        self.label_show.setText(value)  # 更新值

    def update_progress(self, value):
        self.progressBar.setValue(value)  # 更新进度条的值

    def export_finished(self):
        self.progressBar.setValue(100)


class ExportThread(QThread):
    progress = pyqtSignal(int)  # 声明一个信号
    updatelabel = pyqtSignal(str)  # 声明一个更新label新增信号

    def __init__(self, selected_layer1, selected_layer2, csv_file_path):
        super().__init__()
        self.selected_layer1 = selected_layer1
        self.selected_layer2 = selected_layer2
        self.csv_file_path = csv_file_path

    def run(self):
        start_time = time.time()
        # 创建处理上下文
        field_names1 = [f"{field.name()}" for field in self.selected_layer1.fields()] + ['Area平方米']
        field_names2 = [f"{field.name()}" for field in self.selected_layer2.fields()] + ['Area平方米']
        combined_list1 = ([f"{field}_1" for field in field_names1] +
                          [f"{field}_2" for field in field_names2] +
                          ['相交面积Area平方米']
                          # +
                          # ['Polyonwkt']
                          )

        # 创建空间索引
        self.updatelabel.emit("创建索引...")
        spatial_index = QgsSpatialIndex(self.selected_layer2.getFeatures())

        total_features = self.selected_layer1.featureCount()
        self.updatelabel.emit("开始计算...")
        context = QgsDistanceArea()
        context.setEllipsoid('WGS84')
        with open(self.csv_file_path, mode='w', newline='', encoding='gbk') as csv_file:
            writer = csv.writer(csv_file)
            writer.writerow(combined_list1)

            for index, feature1 in enumerate(self.selected_layer1.getFeatures()):
                try:
                    row = []
                    geom1 = feature1.geometry()
                    if geom1:
                        for field_name in field_names1[:-1]:  # 不包括 WKT_BJ
                            value = feature1[field_name]
                            row.append(f'{value}')

                        area1 = context.convertAreaMeasurement(context.measureArea(geom1),
                                                               QgsUnitTypes.AreaSquareMeters)
                        row.append(f'{area1}')  # 添加面积

                        # 判断是否相交
                        intersecting_ids = spatial_index.intersects(geom1.boundingBox())
                        if intersecting_ids:  # 仅处理相交的要素
                            request = QgsFeatureRequest().setFilterFids(intersecting_ids)

                            for feature2 in self.selected_layer2.getFeatures(request):
                                geom2 = feature2.geometry()
                                if geom2 and geom1.intersects(geom2):  # 检查是否相交
                                    row2 = []
                                    intersection_geom = geom1.intersection(geom2)  # 获取相交部分的几何
                                    # 获取相交面积
                                    intersection_area = context.convertAreaMeasurement(
                                        context.measureArea(intersection_geom), QgsUnitTypes.AreaSquareMeters)

                                    for field_name in field_names2[:-1]:
                                        value = feature2[field_name]
                                        row2.append(f'{value}')
                                    area2 = context.convertAreaMeasurement(context.measureArea(geom2),
                                                                           QgsUnitTypes.AreaSquareMeters)
                                    row2.append(f'{area2}')  # 添加面积

                                    if row2:
                                        # 确保 row 不为空
                                        writer.writerow(row + row2 + [intersection_area])
                                       # writer.writerow(row + row2 + [intersection_area]+[intersection_geom.asWkt()])  # 写入相交面积

                        progress_value = int((index + 1) / total_features * 100)
                        self.progress.emit(progress_value)

                except Exception as e:
                    print(f"处理要素 {index} 时出错: {str(e)}")

        end_time = time.time()
        # 计算并打印运行时间
        execution_time = end_time - start_time
        self.updatelabel.emit(f"代码运行时间: {execution_time:.2f} 秒")
