def process(self) -> str:
"""
在这里提取矢量数据的元数据, 将元数据文件存储在self.file_content.work_root_dir下, 固定名称为self.FileName_MetaData, 注意返回的串中有元数据的格式
:return:
"""
result = super().process()
ds_metadata = CFactory().give_me_db(
self.file_info.db_server_id).one_row(
'''
select dsometadatatext, dsometadatajson, dsometadataxml
from dm2_storage_object
where dsoid = :object_id
''', {'object_id': self.object_id})
metadata_filename = CFile.join_file(
self.file_content.work_root_dir,
'{0}.metadata'.format(self.object_name))
text_metadata = ds_metadata.value_by_name(0, 'dsometadatatext', None)
json_metadata = ds_metadata.value_by_name(0, 'dsometadatajson', None)
xml_metadata = ds_metadata.value_by_name(0, 'dsometadataxml', None)
format_metadata = self.MetaDataFormat_Text
if json_metadata is not None:
format_metadata = self.MetaDataFormat_Json
CJson.str_2_file(json_metadata, metadata_filename)
elif xml_metadata is not None:
format_metadata = self.MetaDataFormat_XML
CXml.str_2_file(xml_metadata, metadata_filename)
elif text_metadata is not None:
format_metadata = self.MetaDataFormat_XML
CFile.str_2_file(text_metadata, metadata_filename)
else:
return result
result = CResult.merge_result_info(result, self.Name_FileName,
metadata_filename)
return CResult.merge_result_info(result, self.Name_Format,
format_metadata)
def process_vector(self) -> str:
try:
# file_name_with_full_path = r'D:\test\vector_test\石嘴山市-3xq.json'
# file_main_name = CFile.file_main_name(file_name_with_full_path)
# file_path = CFile.file_path(file_name_with_full_path)
xml_obj = self.metadata.metadata_xml()
# <editor-fold desc="1.空间坐标信息">
wkt_info = 'POLYGON((${min_x} ${max_y},${max_x} ${max_y},' \
'${max_x} ${min_y},${min_x} ${min_y},${min_x} ${max_y}))'
# 四至坐标
native_max_x = CUtils.to_decimal(
xml_obj.get_element_text_by_xpath_one('/TileMetadata/MaxLon'))
native_max_y = CUtils.to_decimal(
xml_obj.get_element_text_by_xpath_one('/TileMetadata/MaxLat'))
native_min_x = CUtils.to_decimal(
xml_obj.get_element_text_by_xpath_one('/TileMetadata/MinLon'))
native_min_y = CUtils.to_decimal(
xml_obj.get_element_text_by_xpath_one('/TileMetadata/MinLat'))
if (native_max_x is None) or (native_max_y is None) \
or (native_min_x is None) or (native_min_y is None):
native_center_wkt = None
native_bbox_wkt = None
geom_native_wkt = None
else:
dict_native = {
'max_x': CUtils.any_2_str(native_max_x),
'max_y': CUtils.any_2_str(native_max_y),
'min_x': CUtils.any_2_str(native_min_x),
'min_y': CUtils.any_2_str(native_min_y)
}
# 中心坐标
center_x = (native_max_x - native_min_x) / 2 + native_min_x
center_y = (native_max_y - native_min_y) / 2 + native_min_y
native_center_wkt = 'POINT({0} {1})'.format(center_x, center_y)
# 外边框、外包框
native_bbox_wkt = CUtils.replace_placeholder(
wkt_info[:], dict_native)
geom_native_wkt = native_bbox_wkt
file_path = self.file_content.work_root_dir
file_main_name = self.object_name
native_center_filepath = CFile.join_file(
file_path, file_main_name + '_native_center.wkt')
CFile.str_2_file(native_center_wkt, native_center_filepath)
native_bbox_filepath = CFile.join_file(
file_path, file_main_name + '_native_bbox.wkt')
CFile.str_2_file(native_bbox_wkt, native_bbox_filepath)
native_geom_filepath = CFile.join_file(
file_path, file_main_name + '_native_geom.wkt')
CFile.str_2_file(geom_native_wkt, native_geom_filepath)
projection = xml_obj.get_element_text_by_xpath_one(
'/TileMetadata/SpatialReference/PRJ')
if (projection is not None) and (not CUtils.equal_ignore_case(
projection, '')):
source_projection = osr.SpatialReference(wkt=projection)
source = source_projection.GetAttrValue('GEOGCS', 0) # 坐标系名称
prosrs = osr.SpatialReference()
if prosrs.ImportFromWkt(projection) == gdal.CE_None:
proj_wkt = prosrs.ExportToWkt()
native_wkt = proj_wkt
else:
native_wkt = projection
native_proj4 = prosrs.ExportToProj4()
spatial = None
rb = (0, 0)
lu = (0, 0)
geosrs = prosrs.CloneGeogCS()
ct = osr.CreateCoordinateTransformation(prosrs, geosrs)
if ct is not None:
rb = ct.TransformPoint(native_max_x, native_max_y)
lu = ct.TransformPoint(native_min_x, native_min_y)
wgs84_min_x = lu[0]
wgs84_max_y = lu[1]
wgs84_max_x = rb[0]
wgs84_min_y = rb[1]
dict_wgs84 = {
'max_x': CUtils.any_2_str(wgs84_max_x),
'max_y': CUtils.any_2_str(wgs84_max_y),
'min_x': CUtils.any_2_str(wgs84_min_x),
'min_y': CUtils.any_2_str(wgs84_min_y)
}
# 中心坐标
center_x = (wgs84_max_x - wgs84_min_x) / 2 + wgs84_min_x
center_y = (wgs84_max_y - wgs84_min_y) / 2 + wgs84_min_y
wgs84_center_wkt = 'POINT({0} {1})'.format(
center_x, center_y)
# 外边框、外包框
wgs84_bbox_wkt = CUtils.replace_placeholder(
wkt_info[:], dict_wgs84)
wgs84_geom_wkt = wgs84_bbox_wkt
wgs84_center_filepath = CFile.join_file(
file_path, file_main_name + '_wgs84_center.wkt')
CFile.str_2_file(wgs84_center_wkt, wgs84_center_filepath)
wgs84_bbox_filepath = CFile.join_file(
file_path, file_main_name + '_wgs84_bbox.wkt')
CFile.str_2_file(wgs84_bbox_wkt, wgs84_bbox_filepath)
wgs84_geom_filepath = CFile.join_file(
file_path, file_main_name + '_wgs84_geom.wkt')
CFile.str_2_file(wgs84_geom_wkt, wgs84_geom_filepath)
# </editor-fold>
# <editor-fold desc="2.投影信息">
native_source = CResource.Prj_Source_Data
# 坐标系
native_coordinate = None
# 3度带/6度带
native_degree = None
# 投影方式
native_project = None
# 带号
native_zone = None
# 创建SpatialReference对象，导入wkt信息
spatial_ref = osr.SpatialReference(wkt=native_wkt)
if spatial_ref.IsProjected():
native_project = spatial_ref.GetAttrValue('PROJECTION')
native_coordinate = spatial_ref.GetAttrValue('GEOGCS')
# native_degree = spatial_ref.GetAttrValue('GEOGCS|UNIT', 1)
native_degree, native_zone = self.get_prj_degree_zone(
spatial_ref)
elif spatial_ref.IsGeocentric():
native_project = None
native_coordinate = spatial_ref.GetAttrValue('GEOGCS')
native_degree = None
native_zone = None
# </editor-fold>
result = CResult.merge_result(self.Success, '处理完毕!')
result = CResult.merge_result_info(result, self.Name_Prj_Wkt,
native_wkt)
result = CResult.merge_result_info(result, self.Name_Prj_Proj4,
native_proj4)
result = CResult.merge_result_info(result,
self.Name_Prj_Project,
native_project)
result = CResult.merge_result_info(result,
self.Name_Prj_Coordinate,
native_coordinate)
result = CResult.merge_result_info(result,
self.Name_Prj_Source,
native_source)
result = CResult.merge_result_info(result, self.Name_Prj_Zone,
native_zone)
result = CResult.merge_result_info(result,
self.Name_Prj_Degree,
native_degree)
return result
# return CResult.merge_result(self.Success, '处理完毕!')
else:
return CResult.merge_result(self.Success, '处理完毕!')
except Exception as error:
CLogger().warning('矢量数据的空间信息处理出现异常, 错误信息为: {0}'.format(
error.__str__()))
return CResult.merge_result(
self.Failure,
'矢量数据的空间信息处理出现异常,错误信息为：{0}!'.format(error.__str__()))
def process_raster(self) -> str:
try:
# file_name_with_full_path = r'D:\test\raster_test\石嘴山市.json'
# file_main_name = CFile.file_main_name(file_name_with_full_path)
# file_path = CFile.file_path(file_name_with_full_path)
json_obj = self.metadata.metadata_json()
# json_obj.load_file(file_name_with_full_path)
# <editor-fold desc="1.空间坐标信息">
wkt_info = 'POLYGON((${min_x} ${max_y},${max_x} ${max_y},' \
'${max_x} ${min_y},${min_x} ${min_y},${min_x} ${max_y}))'
# 四至坐标
native_max_x = CUtils.to_decimal(json_obj.xpath_one('boundingbox.right', None))
native_max_y = CUtils.to_decimal(json_obj.xpath_one('boundingbox.top', None))
native_min_x = CUtils.to_decimal(json_obj.xpath_one('boundingbox.left', None))
native_min_y = CUtils.to_decimal(json_obj.xpath_one('boundingbox.bottom', None))
if (native_max_x is None) or (native_max_y is None) or (native_min_x is None) or (native_min_y is None):
native_center_wkt = None
native_bbox_wkt = None
native_geom_wkt = None
else:
dict_native = {'max_x': CUtils.any_2_str(native_max_x),
'max_y': CUtils.any_2_str(native_max_y),
'min_x': CUtils.any_2_str(native_min_x),
'min_y': CUtils.any_2_str(native_min_y)}
# 中心点坐标
center_x = (native_max_x - native_min_x) / 2 + native_min_x
center_y = (native_max_y - native_min_y) / 2 + native_min_y
native_center_wkt = 'POINT({0} {1})'.format(center_x, center_y)
# 外边框、外包框
native_bbox_wkt = CUtils.replace_placeholder(wkt_info[:], dict_native)
native_geom_wkt = native_bbox_wkt
file_path = self.file_content.work_root_dir
file_main_name = self.object_name
native_center_filepath = CFile.join_file(file_path, file_main_name + '_native_center.wkt')
CFile.str_2_file(native_center_wkt, native_center_filepath)
native_bbox_filepath = CFile.join_file(file_path, file_main_name + '_native_bbox.wkt')
CFile.str_2_file(native_bbox_wkt, native_bbox_filepath)
native_geom_filepath = CFile.join_file(file_path, file_main_name + '_native_geom.wkt')
CFile.str_2_file(native_geom_wkt, native_geom_filepath)
# wgs84转换后的四至坐标
wgs84_max_x = CUtils.to_decimal(json_obj.xpath_one('wgs84.boundingbox.right', None))
wgs84_max_y = CUtils.to_decimal(json_obj.xpath_one('wgs84.boundingbox.top', None))
wgs84_min_x = CUtils.to_decimal(json_obj.xpath_one('wgs84.boundingbox.left', None))
wgs84_min_y = CUtils.to_decimal(json_obj.xpath_one('wgs84.boundingbox.bottom', None))
if (wgs84_max_x is None) or (wgs84_max_y is None) or (wgs84_min_x is None) or (wgs84_min_y is None):
wgs84_center_wkt = None
wgs84_bbox_wkt = None
else:
dict_wgs84 = {'max_x': CUtils.any_2_str(wgs84_max_x),
'max_y': CUtils.any_2_str(wgs84_max_y),
'min_x': CUtils.any_2_str(wgs84_min_x),
'min_y': CUtils.any_2_str(wgs84_min_y)}
# 中心点坐标（wgs84）
center_x = (wgs84_max_x - wgs84_min_x) / 2 + wgs84_min_x
center_y = (wgs84_max_y - wgs84_min_y) / 2 + wgs84_min_y
wgs84_center_wkt = 'POINT({0} {1})'.format(center_x, center_y)
# 外边框、外包框（wgs84)
wgs84_bbox_wkt = CUtils.replace_placeholder(wkt_info[:], dict_wgs84)
# wgs84_geom_wkt = wgs84_bbox_wkt
wgs84_geom_filepath = CFile.join_file(file_path, file_main_name + '_wgs84_geom.wkt')
# CFile.str_2_file(wgs84_geom_wkt, wgs84_geom_filepath)
object_file_path = self.file_info.file_name_with_full_path
process_bus(object_file_path, wgs84_geom_filepath, None)
wgs84_center_filepath = CFile.join_file(file_path, file_main_name + '_wgs84_center.wkt')
CFile.str_2_file(wgs84_center_wkt, wgs84_center_filepath)
wgs84_bbox_filepath = CFile.join_file(file_path, file_main_name + '_wgs84_bbox.wkt')
CFile.str_2_file(wgs84_bbox_wkt, wgs84_bbox_filepath)
# </editor-fold>
# <editor-fold desc="2.投影信息">
native_wkt = json_obj.xpath_one('coordinate.wkt', None)
native_proj4 = json_obj.xpath_one('coordinate.proj4', None)
native_source = CResource.Prj_Source_Data
# 坐标系
native_coordinate = None
# 3度带/6度带
native_degree = None
# 投影方式
native_project = None
# 带号
native_zone = None
# 创建SpatialReference对象，导入wkt信息
spatial_ref = osr.SpatialReference(wkt=native_wkt)
if spatial_ref.IsProjected():
native_coordinate = spatial_ref.GetAttrValue('GEOGCS')
native_project = spatial_ref.GetAttrValue('PROJECTION')
# native_degree = spatial_ref.GetAttrValue('GEOGCS|UNIT', 1)
# pixel_size = json_obj.xpath_one('pixelsize.width', None)
# native_zone = self.get_prj_zone(spatial_ref, pixel_size)
native_degree, native_zone = self.get_prj_degree_zone(spatial_ref)
elif spatial_ref.IsGeocentric():
native_project = None
native_coordinate = spatial_ref.GetAttrValue('GEOGCS')
native_degree = None
native_zone = None
# </editor-fold>
result = CResult.merge_result(self.Success, '处理完毕!')
result = CResult.merge_result_info(result, self.Name_Prj_Wkt, native_wkt)
result = CResult.merge_result_info(result, self.Name_Prj_Proj4, native_proj4)
result = CResult.merge_result_info(result, self.Name_Prj_Project, native_project)
result = CResult.merge_result_info(result, self.Name_Prj_Coordinate, native_coordinate)
result = CResult.merge_result_info(result, self.Name_Prj_Source, native_source)
result = CResult.merge_result_info(result, self.Name_Prj_Zone, native_zone)
result = CResult.merge_result_info(result, self.Name_Prj_Degree, native_degree)
return result
# return CResult.merge_result(self.Success, '处理完毕!')
except Exception as error:
CLogger().warning('影像数据的空间信息处理出现异常, 错误信息为: {0}'.format(error.__str__()))
return CResult.merge_result(self.Failure, '影像数据的空间信息处理出现异常,错误信息为：{0}!'.format(error.__str__()))
