买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：辽宁省地质勘查院有限责任公司

摘要：本发明公开了一种地球化学异常快速统计方法，涉及信息提取方法技术领域，包括资料准备、单元素异常信息统计、计算面积、相交空间分析、汇总统计数据分析、其余字段建立、填充新建字段、汇总数据、结果导出、组合异常信息统计，本发明能够迅速处理大量的地球化学数据，减少了数据处理的时间和人力成本，提高了工作效率；本方法适用于不同类型的地球化学数据和不同的研究区域；无论是土壤地球化学还是水洗沉积物地球化学，该方法都能够进行快速、准确的数据处理和分析。

主权项：1.一种地球化学异常快速统计方法，其特征在于，包括以下步骤：S1、资料准备准备以下数据：1.1地球化学采样点矢量数据，属性中必须包括所有化验元素的化验结果；1.2单元素异常圈定范围矢量数据，属性中必须包括单元素异常编号，每个元素一个文件；1.3组合异常圈定范围矢量数据，属性中必须包括组合异常编号；1.4元素名称及其异常下限值；S2、单元素异常信息统计2.1计算面积在所有单元素异常圈定范围矢量中新建字段，使用字段计算器工具，计算并填充每个单元素异常矢量文件中各个单元素异常区的面积；2.2相交空间分析将地球化学采样点矢量与所有单元素异常圈定范围矢量做相交分析，使单元素异常编号与面积字段值传递至地球化学采样点，对于未处于单元素异常区内的采样点，则不会被传递，以此来做统计区分，将结果命名为“caiyangdian_intersect_元素名称”，每个元素一个文件，以备后用；2.3汇总统计数据分析使用GIS的汇总统计数据分析功能，分别对“caiyangdian_intersect_元素名称”文件做分析，以单元素异常编号字段做数据透视维度，统计出对应元素化验结果的平均值、最大值和面积，将结果命名为“caiyangdian_Intersect_元素名称_Statistics”，每个元素一个文件，以备后用；2.4其余字段建立对所有“caiyangdian_Intersect_元素名称_Statistics”表建立“元素”、“衬度”、“规模”、“浓度分带”和“异常下限”字段，其中“元素”和“浓度分带”为字符串类型，其余均为双精度类型；2.5填充新建字段①元素字段填充公式：当前处理的元素名称+异常编号；②异常下限字段填充当前处理的元素的异常下限值即可；③衬度字段填充公式：元素分析结果平均值异常下限；④规模字段填充公式：单元素异常面积*衬度；⑤浓度分带填充基于最大值和异常下限字段；2.6汇总数据将所有“caiyangdian_Intersect_元素名称_Statistics”表转换为Python的DataFrame对象，并依据元素字段对结果排序；完成后，对DataFrame进行行列转置；2.7结果导出将所有DataFrame对象转化为Excel表的sheet对象，并合并成一个Excel文件后导出，地球化学单元素异常统计完成；S3、组合异常信息统计3.1面积与范围计算为组合异常圈定范围矢量文件增加面积、组合异常X最大值、组合异常X最小值、组合异常Y最大值和组合异常Y最小值字段，类型均为双精度，长度20、小数点长度6；使用字段计算器工具，计算并填充组合异常圈定范围矢量的面积，计算空间参考系应使用工作区范围内的国家2000坐标系高斯3°带投影，面积单位选择平方千米；使用字段计算器工具，计算并填充组合异常圈定范围矢量的组合异常X最大值、组合异常X最小值、组合异常Y最大值和组合异常Y最小值，计算空间参考系应使用工作区范围内的国家2000坐标系高斯3°带投影，坐标单位选择米；3.2剔除小于异常下限的数值为不修改原始采样点化验结果数据，复制一份地球化学采样点矢量数据，并命名为“temp_caiyangdian”，根据元素异常下限，将“temp_caiyangdian”中所有元素化验结果值小于异常下限者替换为None；3.3相交空间分析将组合异常区与“temp_caiyangdian”做相交分析，使组合异常编号、面积、组合异常X最大值、组合异常X最小值、组合异常Y最大值和组合异常Y最小值传递至“temp_caiyangdian”文件中，将结果命名为“temp_Intersect”以备后用；3.4字段统计数据转表使用GIS的字段统计数据转表分析功能，对“temp_Intersect”进行分析，将所有元素化验结果字段作为分析项，以组合异常编号字段做数据透视维度，统计各个元素在组合异常区内的元素名称、最高值、平均值和标准差，将结果命名为“temp_Intersect_FieldStatisticsToTable_Numeric”，以备后用；3.5其余字段建立在“temp_Intersect_FieldStatisticsToTable_Numeric”文件中建立异常下限、衬度和浓度分带字段，其中“浓度分带”为字符串类型，其余均为双精度类型；3.6填充新建字段①异常下限字段填充元素的异常下限值即可；②衬度字段填充公式：“temp_Intersect_FieldStatisticsToTable_Numeric”表中平均值异常下限；③浓度分带填充基于“temp_Intersect_FieldStatisticsToTable_Numeric”表中最大值和当前元素的异常下限字段，其计算方式与单元素异常中浓度分带计算方式一致；3.7生成组合异常数据DataFrame将“temp_Intersect_FieldStatisticsToTable_Numeric”表转为Python的DataFrame对象，其中每个组合异常单独转为一个DataFrame对象；接着先将所有DataFrame对象转置行列，再根据异常编号拼接组合异常圈定范围矢量的组合异常X最大值、组合异常X最小值、组合异常Y最大值和组合异常Y最小值，所得结果暂存，以备后用；3.8单元素异常与组合异常相交空间分析将所有单元素异常与组合异常做相交空间分析，分别获得每个元素的单元素异常在组合异常中的部分，然后建立相交面积字段，类型为双精度，长度20、小数点长度6；使用字段计算器工具，计算并填充相交面积，计算空间参考系应使用工作区范围内的国家2000坐标系高斯3°带投影，面积单位选择平方千米；将结果命名为“temp_insect_all_元素名称”，以备后用；3.9聚合分析使用GIS的聚合分析功能，将所有“temp_insect_all_元素名称”文件分别做元素面积聚合统计，并将统计结果表转为Python的DataFrame对象，并转置行列；所得结果暂存，以备后用；3.10DataFrame拼接将3.7和3.9的DataFrame对象，依据组合异常编号配对，再分别根据每对中的元素名称进行拼接得到以组合异常编号相同的一组DataFrame对象；3.11结果导出将所有DataFrame对象按照组合异常编号排序，并转化为Excel表的sheet对象，将sheet合并成一个Excel文件后导出，地球化学组合异常统计完成。

全文数据：

权利要求：

百度查询： 辽宁省地质勘查院有限责任公司 一种地球化学异常快速统计方法