粗糙度中方位角添加方式待修改

surfaceRoughness_oh2004/SoilMoistureMain.py

@@ -36,6 +36,7 @@ import datetime
 import numpy as np
 import sys
 import multiprocessing
+import xml.etree.ElementTree as ET
 cover_id_list = []
 threshold_of_ndvi_min = 0
 threshold_of_ndvi_max = 0
@@ -315,8 +316,15 @@ class MoistureMain:
         par_dict = CreateDict(image_path, self.processinfo, out_path1, out_path2).calu_nature(start)
 
         # TODO:从原始文件的xml中获取得到方位角信息 SurfaceRoughness\Temporary\preprocessing\影像名\ 文件夹下 -ortho.meta.xml
-        Azimuth=self.__processing_paras['META']
-        par_dict.update({"imageinfo_SatelliteAzimuth": "None"})
+        # 这个添加方式有问题 待修改
+        # 加载XML文件
+        # tree = ET.parse(self.__processing_paras['META'])
+        # ## 获取根节点
+        # root = tree.getroot()
+        # # # 找到SatelliteAzimuth标签并获取其值
+        # satellite_azimuth = root.find('.//SatelliteAzimuth')
+        # if satellite_azimuth is not  None:
+        #     par_dict.update({"productinfo_SatelliteAzimuth": satellite_azimuth})
 
 
         model_xml_path = os.path.join(tem_folder, "creat_standard.meta.xml")  # 输出xml路径