CN105955939A - 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法 - Google Patents

一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法 Download PDF

Info

Publication number
CN105955939A
CN105955939A CN201610269353.0A CN201610269353A CN105955939A CN 105955939 A CN105955939 A CN 105955939A CN 201610269353 A CN201610269353 A CN 201610269353A CN 105955939 A CN105955939 A CN 105955939A
Authority
CN
China
Prior art keywords
inputs
input
tendency
rose
coordinate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610269353.0A
Other languages
English (en)
Inventor
罗可
钟建华
倪良田
孙宁亮
刘闯
李莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201610269353.0A priority Critical patent/CN105955939A/zh
Publication of CN105955939A publication Critical patent/CN105955939A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F40/00—Handling natural language data
    • G06F40/10—Text processing
    • G06F40/166—Editing, e.g. inserting or deleting
    • G06F40/174—Form filling; Merging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明涉及一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，通过使用Microsoft Excel2010软件制作倾斜玫瑰花图绘制的Excel模板文件。其技术方案是：1.制作Excel表格。包括在表格中输入所测量的地质体产状的原始数据；处理原始数据，对倾向及倾角进行分组统计并换算为方位角坐标；刻画方位角坐标系，求得圆坐标以及方位角刻度线坐标。2.绘制倾斜玫瑰花图。将处理所得的圆坐标、方位角刻度线坐标、倾向和倾角坐标通过散点图的形式绘制成倾斜玫瑰花图。优点是：只需输入原始数据，就能自动进行数据统计和处理并且成图，方法简单、易于操作、灵活多变、适应性强、数据容量大、准确性高、图件美观，使得倾斜玫瑰花图的绘制时间大大缩减，提高了工作效率。

Description

一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法
技术领域
本发明涉及地质学野外测量数据处理领域，尤其涉及一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法。
背景技术
倾斜玫瑰花图可以定性而形象地反映某一地质体的倾向及倾角的优势分布，在地质学中应用广泛。沉积学中将交错层理的产状绘制成倾斜玫瑰花图来分析古流向及古风向等，构造学中则将节理的产状绘制成倾斜玫瑰花图来表示一个地区节理发育的方向和特点，用于构造应力分析。目前，绘制倾斜玫瑰花图主要有以下两种做法：
一种是人工制图，这种制图方法存在明显缺陷。首先制图前需要耗费大量时间进行数据统计，统计过程中容易出错，制图过程多是机械重复；其次，不同人之间的绘图标准也不尽相同，无法统一使用；而且，手工制图的图件存储和提取不方便，无法满足数字化时代的需求，使用起来也有诸多不便。
第二种是采用Grapher软件制图。该软件绘制的玫瑰花图目前应用较多，但仍存在其局限性。首先，Grapher软件绘制的玫瑰花图只能对产状数据分区间做单纯的统计分组工作，而我们地质中实际应用的玫瑰花图是由每个统计区间数据的平均值连接绘制而成；此外，倾斜玫瑰花图中倾向和倾角表征方式不一，一般用半径的长度来代表平均倾角，Grapher软件无法做到让两者统一。因此仍然不能很好地满足绘图需要。
要解决上述不足,技术人员通过有针对性的编程当然最好,但事实上,本领域内的技术人员多数不懂编程,而专业编程人员又对倾斜玫瑰花图的原理和运用了解不深,往往不能准确理解技术人员在编程方面的设想和需求。总之，目前绘制倾斜玫瑰花图的方法存在较大不足，不能适应地质数据处理的需要。
发明内容
本发明的目的在于提供一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，该方法使得倾斜玫瑰花图的绘制快速、高效、准确，且各个统计区间内对应的倾向和倾角达到高度一致。
本发明实施例是这样实现的，一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，通过使用微软系统自带的Microsoft Excel 2010软件，利用该软件强大的数据处理、统计分析以及绘图功能，结合倾斜玫瑰花图的绘制方法与原理，设定所需参数及运算公式，自动对原始数据进行统计、处理并且成图，这是倾斜玫瑰花图绘制方法的核心。其Excel模板文件的具体实施包括如下步骤：
1.新建工作表
新建工作薄中，建立工作表标签名为“倾斜玫瑰花图绘制”的工作表。工作表主要包括原始数据输入、产状数据统计与处理和方位角坐标系刻画三个部分。
2.各单元格输入内容
首先，设置对应的表头内容如“倾向°”、“倾角°”等，以及对应的信息输入区。
A1-V1合并及居中，输入标题“倾斜玫瑰花图绘制”；
A2-C2、A3-C3合并及居中，输入“原始数据输入”，A4输入“编号”，B4输入“倾向°”，C4输入“倾角°”，此为原始数据输入部分表头；
D2-N2合并及居中，输入“产状数据统计与处理”，D3-J3合并及居中，输入“倾向参数”，K3-N3合并及居中，输入“倾角参数”，D4输入“倾向范围°”，E4输入“条数”，F4输入“倾向值和°”，G4输入“倾向范围°”，H4输入“弧度1”，I1输入“x1”，J1输入“y1”，K4输入“倾角值和°”，L4输入“平均倾角°”，M4输入“x2”，N4输入“y2”，此为产状数据统计与处理部分表头；
O2-V2合并及居中，输入“方位角坐标系刻画”，O3-P3合并及居中输入“圆”，Q2输入“半径r1”，O4输入“角度2”，P4输入“弧度2”，Q4输入“x3”，R4输入“y3”；S3-T3合并及居中输入“方位角刻度线”，U3输入“半径r2”，S4输入“角度3”，T4输入“弧度3”，U4输入“x4”，V4输入“y4”，此为方位角坐标系刻画部分表头。
进一步，输入数据并对数据进行处理。
1)原始数据输入：
自A5起的A列单元格为原始数据的编号，A5输入“1”，A5、A6、A7……分别对应的编号为“1”、“2”、“3”……，列数取决于原始数据的组数；在相对应的B列和C列单元格则分别输入相对应的倾向和倾角数值。
2)产状数据统计与处理：
第一步，对倾向参数进行统计与处理，并且进行坐标换算。
D列输入倾向范围。将0°～360°在D列中按公差10°进行分组，即D5、D6、D7……D40分别输入“0～10”、“10～20”、“20～30”……“350～360”。
E列利用Excel函数和公式统计D列各组倾向范围内对应的产状条数。
E5＝COUNTIF(B5:B10000,"<10")-COUNTIF(B5:B10000,"<0")。利用相同方法完成其它各组的产状统计，把公式分别输入E6到E40单元格内。
F列对各组倾向值求和。利用Excel数组公式求解，以F5为例，需求得倾向范围在0°～10°的所有倾向值的总和，求解公式为：“{＝SUMIF($B$5:$B$10000,"<10")-SUMIF($B$5:$B$10000,"<0")}”。利用相同方法完成其它各组的求和，把公式分别输入F6到F40单元格内。
G列求取各组倾向范围内的平均倾向值。以G5为例，G5＝IF(F5＝0,"0",F5/E5)，下拉复制公式到G40，完成其它各组的平均倾向值计算。
H列将G列的平均倾向转化为弧度，I、J列为各倾向组的坐标，x1＝条数×sin(平均倾向，弧度1)，y1＝条数×cos(平均走向，弧度1)。如H26、I26、J26代表210°～220°组坐标点，H26＝RADIANS(G26)，I26＝E26*SIN(H26)，J26＝E26*COS(H26)。用上述方法，在I、J列求出各倾向组的坐标。
第二步，对倾角参数进行统计与处理，并且进行坐标换算。
K列求取各组对应的倾角值和。以K5为例，需求得倾向范围在0°～10°的所有倾角值的总和，求解公式为：“{＝SUM(IF(($B$5:$B$10000<10)*($B$5:$B$10000>＝0),$C$5:$C$10000))}”。利用相同方法完成其它各组的求和，把公式分别输入K6到K40单元格内。
L列求取各组倾向范围内的平均倾向值。以L5为例，L5＝IF(K5＝0,"0",K5/E5)，下拉复制公式到K40，完成其它各组的平均倾向值计算。
M、N列求取各组平均倾角的坐标。平均倾角的方位角与平均倾向方位角一致，倾角的大小用半径的长度来表示，因此在R3单元格求出圆的半径r1，圆半径取E列中的最大值，求解公式为：
“{＝MAX(E5:E40)}”。x2＝sin(平均倾向)×圆半径r1×平均倾角/90°，y2＝cos(平均倾向)×圆半径r1×平均倾角/90°。如M26、N26代表210°～220°组坐标点，M26＝SIN(H26)*$R$3*L26/90，N26＝COS(H26)*$R$3*L26/90(Excel中，*代表乘号，下同)。用上述方法，在M、N列求出各倾向组倾角的坐标。
3)方位角坐标系刻画：
方位角坐标系的刻画包括以最大条数为半径的圆和以10°为公差的方位角刻度线的刻画。在O列中依次输入270～350和0～270，公差为10，即O5、O6……O13、O14……O41依次输入“270”、“280”……“350”、“0”……“270”。在P列将O列角度转化为对应弧度，再在Q、R列求出其坐标。如P5、Q5、R5代表270°所对应圆的坐标，P5＝RADIANS(O5)，Q5＝$R$3*SIN(P5)，R5＝$R$3*COS(P5)。用上述方法，求出圆的坐标。
U、V列为各角度对应刻度线的坐标，每条刻度线对应两个坐标，一个是刻度线外侧点坐标，一个 是刻度线内侧点坐标。内侧点坐标设定为圆半径的0.975倍，即r2＝0.975r1，在V3单元格中输入“＝0.975*R3”。以270°为例，S5、T5、U5、V5为刻度线内侧点坐标，T5＝RADIANS(S5)，U5＝$V$3*SIN(T5)，V5＝$V$3*COS(T5)；S6、T6、U6、V6为刻度线内侧点坐标，T6＝RADIANS(S6)，U6＝$V$3*SIN(T6)，V6＝$V$3*COS(T6)。S列排序与O列一致，相邻角度之间需插入一行空白单元格，使刻度线不连续。用上述方法，求出方位角刻度线的坐标。
3.绘制倾斜玫瑰花图
通过以上的数据统计与处理，可进行倾斜玫瑰花图的绘制。在工作表空白处，插入“带平滑线的散点图”。
第一步，添加圆的数据。右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”输入空格，“X轴系列值”选取Q5：Q41，“Y轴系列值”选取R5：R41，点击“确定”。
第二步，添加刻度线的数据。再次单击“添加”，“系列名称”输入空格，“X轴系列值”选取U5：U114，“Y轴系列值”选取V5：V114，点击“确定”。
第三步，添加平均倾向数据。再次单击“添加”，“系列名称”输入倾向，“X轴系列值”选取I5：I40，“Y轴系列值”选取J5：J40。选择所绘制的倾向玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”。
第四步，添加平均倾角数据。再次单击“添加”，“系列名称”输入倾角，“X轴系列值”选取M5：M40，“Y轴系列值”选取N5：N40。选择所绘制的倾角玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”。
第五步，加入总数据量标注、半径标注以及方位角标注。在D41单元格输入“总数”，E41单元格计算总条数，E41＝SUM(E5:E40)；在D42单元格输入“半径”，E42＝R3。右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”选取D41：E41，点击“确定”，完成总数标注添加；再次单击“添加”，“系列名称”选取D42：E42，点击“确定”，完成半径标注添加；再次单击“添加”，“系列名称”输入“0°”，点击“确定”，按此方法依次添加“90°”、“180°”、“270°”，将其分别移动至合适位置，完成方位角的标注添加。
通过上述步骤，完成了对倾斜玫瑰花图的绘制。进一步，把制作好的倾斜玫瑰花图绘制文件属性设为只读，以免在使用过程中大意以致改动，保存好以供随时使用。以后调用该工作表时，只需在工作表中输入地质体产状原始数据倾向和倾角，就能立即得到倾斜玫瑰花图。以上为整个Excel模板文件制作过程。
附图说明
图1是Excel模板文件中工作表的设计以及各单元格输入内容和格式示例图。
图2是实施例所生成的倾斜玫瑰花图。
具体实施方式
下面结合附图及实施例对本发明作进一步的说明。
1.新建工作表
新建工作薄中，建立工作表标签名为“倾斜玫瑰花图绘制”的工作表。工作表主要包括原始数据输入、产状数据统计与处理和方位角坐标系刻画三个部分，详情见图1(只截取其中局部)。
2.各单元格输入内容
首先，设置对应的表头内容如“倾向°”、“倾角°”等，以及对应的信息输入区，详情见图1。
A1-V1合并及居中，输入标题“倾斜玫瑰花图绘制”；
A2-C2、A3-C3合并及居中，输入“原始数据输入”，A4输入“编号”，B4输入“倾向°”，C4输入“倾角°”，此为原始数据输入部分表头；
D2-N2合并及居中，输入“产状数据统计与处理”，D3-J3合并及居中，输入“倾向参数”，K3-N3合并及居中，输入“倾角参数”，D4输入“倾向范围°”，E4输入“条数”，F4输入“倾向值和°”，G4输入“倾向范围°”，H4输入“弧度1”，I1输入“x1”，J1输入“y1”，K4输入“倾角值和°”，L4输入“平均倾角°”，M4输入“x2”，N4输入“y2”，此为产状数据统计与处理部分表头；
O2-V2合并及居中，输入“方位角坐标系刻画”，O3-P3合并及居中输入“圆”，Q2输入“半径r1”，O4输入“角度2”，P4输入“弧度2”，Q4输入“x3”，R4输入“y3”；S3-T3合并及居中输入“方位角刻度线”，U3输入“半径r2”，S4输入“角度3”，T4输入“弧度3”，U4输入“x4”，V4输入“y4”，此为方位角坐标系刻画部分表头。
进一步，结合实测数据对数据的输入与处理作进一步说明，详情见图1。
1)原始数据输入：
自A5起的A列单元格为原始数据的编号，A5输入“1”，A5、A6、A7……分别对应的编号为“1”、“2”、“3”……，列数取决于原始数据的组数；在相对应的B列和C列单元格则分别输入相对应的倾向和倾角数值。
2)产状数据统计与处理：
第一步，对倾向参数进行统计与处理，并且进行坐标换算。
D列输入倾向范围。将0°～360°在D列中按公差10°进行分组，即D5、D6、D7……D40分别输入“0～10”、“10～20”、“20～30”……“350～360”，从而可以分区间对产状数据进行统计和处理；
E列利用Excel函数和公式统计D列各组倾向范围内对应的产状条数。
E5＝COUNTIF(B5:B10000,"<10")-COUNTIF(B5:B10000,"<0")(为了增强表格的计算能力，实现对大量数据进行统计的功能，将计算区间设定为B5到B10000单元格，可根据实际情况进行改动，下同)，各倾向范围取下限值，不取上限值。利用相同方法完成其它各组的产状统计，把公式分别输入E6到E40单元格内。
F列对各组倾向值求和。利用Excel数组公式求解，以F5为例，需求得倾向范围在0°～10°的所有倾向值的总和，求解公式为：“{＝SUMIF($B$5:$B$10000,"<10")-SUMIF($B$5:$B$10000,"<0")}”(公式两侧的大括号并非键盘输入，而是完成等号后面的公式后，按“Ctrl+Shift+Enter”组合键，由Excel自动添加的，这是Excel中数组公式与普通公式的区别，下同)。利用相同方法完成其它各组的求和，把公式分别输入F6到F40单元格内。
G列求取各组倾向范围内的平均倾向值。以G5为例，G5＝IF(F5＝0,"0",F5/E5)，下拉复制公式到G40，完成其它各组的平均倾向值计算。
H列将G列的平均倾向转化为弧度，I、J列为各倾向组的坐标，x1＝条数×sin(平均倾向，弧度1)，y1＝条数×cos(平均走向，弧度1)。如H26、I26、J26代表210°～220°组坐标点，H26＝RADIANS(G26)，I26＝E26*SIN(H26)，J26＝E26*COS(H26)。用上述方法，在I、J列求出各倾向组的坐标。
第二步，对倾角参数进行统计与处理，并且进行坐标换算。
K列求取各组对应的倾角值和。以K5为例，需求得倾向范围在0°～10°的所有倾角值的总和，求解公式为：“{＝SUM(IF(($B$5:$B$10000<10)*($B$5:$B$10000>＝0),$C$5:$C$10000))}”。利用相同方法完成其它各组的求和，把公式分别输入K6到K40单元格内。
L列求取各组倾向范围内的平均倾向值。以L5为例，L5＝IF(K5＝0,"0",K5/E5)，下拉复制公式到K40，完成其它各组的平均倾向值计算。
M、N列求取各组平均倾角的坐标。平均倾角的方位角与平均倾向方位角一致，倾角的大小用半径的长度来表示，因此在R3单元格求出圆的半径r1，圆半径取E列中的最大值，求解公式为：
“{＝MAX(E5:E40)}”。x2＝sin(平均倾向)×圆半径r1×平均倾角/90°，y2＝cos(平均倾向)×圆半径r1×平均倾角/90°。如M26、N26代表210°～220°组坐标点，M26＝SIN(H26)*$R$3*L26/90， N26＝COS(H26)*$R$3*L26/90(Excel中，*代表乘号，下同)。用上述方法，在M、N列求出各倾向组倾角的坐标。
3)方位角坐标系刻画：
方位角坐标系的刻画包括以最大条数为半径的圆和以10°为公差的方位角刻度线的刻画。在O列中按图1中的顺序和格式输入270～350和0～270，公差为10，即O5、O6……O13、O14……O41依次输入“270”、“280”……“350”、“0”……“270”。在P列将O列角度转化为对应弧度，再在Q、R列求出其坐标。如P5、Q5、R5代表270°所对应圆的坐标，P5＝RADIANS(O5)，Q5＝$R$3*SIN(P5)，R5＝$R$3*COS(P5)。用上述方法，求出圆的坐标。
U、V列为各角度对应刻度线的坐标，每条刻度线对应两个坐标，一个是刻度线外侧点坐标，一个是刻度线内侧点坐标。内侧点坐标设定为圆半径的0.975倍，即r2＝0.975r1，在V3单元格中输入“＝0.975*R3”。以270°为例，S5、T5、U5、V5为刻度线内侧点坐标，T5＝RADIANS(S5)，U5＝$V$3*SIN(T5)，V5＝$V$3*COS(T5)；S6、T6、U6、V6为刻度线内侧点坐标，T6＝RADIANS(S6)，U6＝$V$3*SIN(T6)，V6＝$V$3*COS(T6)。S列排序与O列一致，相邻角度之间需插入一行空白单元格，使刻度线不连续。用上述方法，求出方位角刻度线的坐标。
3.绘制倾斜玫瑰花图
通过以上的数据统计与处理，可进行倾斜玫瑰花图的绘制。在工作表空白处，插入“带平滑线的散点图”。
第一步，添加圆的数据。右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”输入空格，“X轴系列值”选取Q5：Q41，“Y轴系列值”选取R5：R41，点击“确定”。
第二步，添加刻度线的数据。再次单击“添加”，“系列名称”输入空格，“X轴系列值”选取U5：U114，“Y轴系列值”选取V5：V114，点击“确定”。
第三步，添加平均倾向数据。再次单击“添加”，“系列名称”输入倾向，“X轴系列值”选取I5：I40，“Y轴系列值”选取J5：J40。选择所绘制的倾向玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”。
第四步，添加平均倾角数据。再次单击“添加”，“系列名称”输入倾角，“X轴系列值”选取M5：M40，“Y轴系列值”选取N5：N40。选择所绘制的倾角玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”。
第五步，加入总数据量标注、半径标注以及方位角标注。在D41单元格输入“总数”，E41单元格计算总条数，E41＝SUM(E5:E40)；在D42单元格输入“半径”，E42＝R3。右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”选取D41：E41，点击“确定”，完成总数标注添加；再次单击“添加”，“系列名称”选取D42：E42，点击“确定”，完成半径标注添加；再次单击“添加”，“系列名称”输入“0°”，点击“确定”，按此方法依次添加“90°”、“180°”、“270°”，将其分别移动至合适位置，完成方位角的标注添加。
通过上述步骤，完成了对倾斜玫瑰花图的绘制，实施例生成的倾斜玫瑰花图如图2，根据实际需要，可对图中文字和线条的颜色、线型、粗细进行调整。进一步，把制作好的倾斜玫瑰花图绘制文件属性设为只读，以免在使用过程中大意以致改动，保存好以供随时使用。以后调用该工作表时，只需在工作表中输入地质体产状原始数据倾向和倾角，就能立即得到倾斜玫瑰花图。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，其特征是：通过使用微软系统自带的Microsoft Excel 2010软件，利用该软件强大的数据处理、统计分析以及绘图功能，结合倾斜玫瑰花图的绘制方法与原理，设定所需参数及运算公式，自动对原始数据进行统计、处理并且成图，这是倾斜玫瑰花图绘制方法的核心。
2.如权利要求1所述的基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，其特征是：Excel表格分为三部分，第一部分为原始数据输入，输入实测的倾向、倾角数据，并进行编号；第二部分为产状数据统计与处理，对倾向及倾角进行分组统计并换算为方位角坐标；第三部分为方位角坐标系刻画，求得圆坐标以及方位角刻度线坐标。
3.如权利要求1所述的基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，其特征是：绘制倾斜玫瑰花图，将表格中处理所得的圆坐标、方位角刻度线坐标、倾向和倾角坐标通过散点图的形式绘制成倾斜玫瑰花图。
4.如权利要求1所述的基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法，其特征在于绘制倾斜玫瑰花图的Excel模板文件的具体制作,包括以下步骤：
1)新建工作表
新建工作薄中，建立工作表标签名为“倾斜玫瑰花图绘制”的工作表；
2)各单元格输入内容
首先，设置对应的表头内容如“倾向°”、“倾角°”等，以及对应的信息输入区；
A1-V1合并及居中，输入标题“倾斜玫瑰花图绘制”；
A2-C2、A3-C3合并及居中，输入“原始数据输入”，A4输入“编号”，B4输入“倾向°”，C4输入“倾角°”；
D2-N2合并及居中，输入“产状数据统计与处理”，D3-J3合并及居中，输入“倾向参数”，K3-N3合并及居中，输入“倾角参数”，D4输入“倾向范围°”，E4输入“条数”，F4输入“倾向值和°”，G4输入“倾向范围°”，H4输入“弧度1”，I1输入“x1”，J1输入“y1”，K4输入“倾角值和°”，L4输入“平均倾角°”，M4输入“x2”，N4输入“y2”；
O2-V2合并及居中，输入“方位角坐标系刻画”，O3-P3合并及居中输入“圆”，Q2输入“半径r1”，O4输入“角度2”，P4输入“弧度2”，Q4输入“x3”，R4输入“y3”；S3-T3合并及居中输入“方位角刻度线”，U3输入“半径r2”，S4输入“角度3”，T4输入“弧度3”，U4输入“x4”，V4输入“y4”；
进一步，输入数据并对数据进行处理；
(1)原始数据输入：
自A5起的A列单元格为原始数据的编号，A5输入“1”，A5、A6、A7……分别对应的编号为“1”、“2”、“3”……，列数取决于原始数据的组数；在相对应的B列和C列单元格则分别输入相对应的倾向和倾角数值；
(2)产状数据统计与处理：
第一步，对倾向参数进行统计与处理，并且进行坐标换算；
D列输入倾向范围，将0°～360°在D列中按公差10°进行分组，即D5、D6、D7……D40分别输入“0～10”、“10～20”、“20～30”……“350～360”；
E列利用Excel函数和公式统计D列各组倾向范围内对应的产状条数，E5＝COUNTIF(B5:B10000,"<10")-COUNTIF(B5:B10000,"<0")，利用相同方法完成其它各组的产状统计，把公式分别输入E6到E40单元格内；
F列对各组倾向值求和，利用Excel数组公式求解，以F5为例，需求得倾向范围在0°～10°的所有倾向值的总和，求解公式为：“{＝SUMIF($B$5:$B$10000,"<10")-SUMIF($B$5:$B$10000,"<0")}”，利用相同方法完成其它各组的求和，把公式分别输入F6到F40单元格内；
G列求取各组倾向范围内的平均倾向值，以G5为例，G5＝IF(F5＝0,"0",F5/E5)，下拉复制公式到G40，完成其它各组的平均倾向值计算；
H列将G列的平均倾向转化为弧度，I、J列为各倾向组的坐标，x1＝条数×sin(平均倾向，弧度1)，y1＝条数×cos(平均走向，弧度1)，如H26、I26、J26代表210°～220°组坐标点，H26＝RADIANS(G26)，I26＝E26*SIN(H26)，J26＝E26*COS(H26)，用上述方法，在I、J列求出各倾向组的坐标；
第二步，对倾角参数进行统计与处理，并且进行坐标换算；
K列求取各组对应的倾角值和，以K5为例，需求得倾向范围在0°～10°的所有倾角值的总和，求解公式为：“{＝SUM(IF(($B$5:$B$10000<10)*($B$5:$B$10000>＝0),$C$5:$C$10000))}”，利用相同方法完成其它各组的求和，把公式分别输入K6到K40单元格内；
L列求取各组倾向范围内的平均倾向值，以L5为例，L5＝IF(K5＝0,"0",K5/E5)，下拉复制公式到K40，完成其它各组的平均倾向值计算；
M、N列求取各组平均倾角的坐标，平均倾角的方位角与平均倾向方位角一致，倾角的大小用半径的长度来表示，因此在R3单元格求出圆的半径r1，圆半径取E列中的最大值，求解公式为：“{＝MAX(E5:E40)}”，x2＝sin(平均倾向)×圆半径r1×平均倾角/90°，y2＝cos(平均倾向)×圆半径r1×平均倾角/90°，如M26、N26代表210°～220°组坐标点，M26＝SIN(H26)*$R$3*L26/90，N26＝COS(H26)*$R$3*L26/90(Excel中，*代表乘号，下同)，用上述方法，在M、N列求出各倾向组倾角的坐标；
(3)方位角坐标系刻画：
方位角坐标系的刻画包括以最大条数为半径的圆和以10°为公差的方位角刻度线的刻画，在O列中依次输入270～350和0～270，公差为10，即O5、O6……O13、O14……O41依次输入“270”、“280”……“350”、“0”……“270”，在P列将O列角度转化为对应弧度，再在Q、R列求出其坐标，如P5、Q5、R5代表270°所对应圆的坐标，P5＝RADIANS(O5)，Q5＝$R$3*SIN(P5)，R5＝$R$3*COS(P5)，用上述方法，求出圆的坐标；
U、V列为各角度对应刻度线的坐标，每条刻度线对应两个坐标，内侧点坐标设定为圆半径的0.975倍，即r2＝0.975r1，在V3单元格中输入“＝0.975*R3”，以270°为例，S5、T5、U5、V5为刻度线内侧点坐标，T5＝RADIANS(S5)，U5＝$V$3*SIN(T5)，V5＝$V$3*COS(T5)；S6、T6、U6、V6为刻度线内侧点坐标，T6＝RADIANS(S6)，U6＝$V$3*SIN(T6)，V6＝$V$3*COS(T6)，S列排序与O列一致，相邻角度之间需插入一行空白单元格，使刻度线不连续，用上述方法，求出方位角刻度线的坐标；
3)绘制倾斜玫瑰花图
通过以上的数据统计与处理，可进行倾斜玫瑰花图的绘制，在工作表空白处，插入“带平滑线的散点图”；
第一步，添加圆的数据，右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”输入空格，“X轴系列值”选取Q5：Q41，“Y轴系列值”选取R5：R41，点击“确定”；
第二步，添加刻度线的数据，再次单击“添加”，“系列名称”输入空格，“X轴系列值”选取U5：U114，“Y轴系列值”选取V5：V114，点击“确定”；
第三步，添加平均倾向数据，再次单击“添加”，“系列名称”输入倾向，“X轴系列值”选取I5：I40，“Y轴系列值”选取J5：J40，选择所绘制的倾向玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”；
第四步，添加平均倾角数据，再次单击“添加”，“系列名称”输入倾角，“X轴系列值”选取M5：M40，“Y轴系列值”选取N5：N40，选择所绘制的倾角玫瑰花图曲线，右键单击，选择“更改系列图表类型”，选择“带直线的散点图”，点击“确定”；
第五步，加入总数据量标注、半径标注以及方位角标注，在D41单元格输入“总数”，E41单元格计算总条数，E41＝SUM(E5:E40)；在D42单元格输入“半径”，E42＝R3；右键单击图表区域，选择“选择数据”，在弹出的对话框中，单击“添加”，在弹出的对话框中，“系列名称”选取D41：E41，点击“确定”，完成总数标注添加；再次单击“添加”，“系列名称”选取D42：E42，点击“确定”，完成半径标注添加；再次单击“添加”，“系列名称”输入“0°”，点击“确定”，按此方法依次添加“90°”、“180°”、“270°”，将其分别移动至合适位置，完成方位角的标注添加。
CN201610269353.0A 2016-04-27 2016-04-27 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法 Pending CN105955939A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610269353.0A CN105955939A (zh) 2016-04-27 2016-04-27 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610269353.0A CN105955939A (zh) 2016-04-27 2016-04-27 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法

Publications (1)

Publication Number Publication Date
CN105955939A true CN105955939A (zh) 2016-09-21

Family

ID=56916578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610269353.0A Pending CN105955939A (zh) 2016-04-27 2016-04-27 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法

Country Status (1)

Country Link
CN (1) CN105955939A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874251A (zh) * 2017-03-10 2017-06-20 中国电建集团成都勘测设计研究院有限公司 一种岩体rbi和rqd指标的测量方法
CN107389555A (zh) * 2017-07-21 2017-11-24 中水北方勘测设计研究有限责任公司 一种分析钻孔岩体主要节理发育方向的方法
CN107590849A (zh) * 2017-08-09 2018-01-16 山西晋城无烟煤矿业集团有限责任公司 一种钻孔自动成图方法
CN109637324A (zh) * 2018-12-19 2019-04-16 西南石油大学 一种原始纹层产状恢复的方法
CN111063035A (zh) * 2019-11-22 2020-04-24 同济大学 一种gis中od关系的三维可视化方法及装置
CN111159862A (zh) * 2019-12-17 2020-05-15 核工业二四三大队 一种铀矿钻孔设计指示书电子化设计及编制方法
CN112767515A (zh) * 2021-02-19 2021-05-07 上海市政工程设计研究总院（集团）有限公司 测试点地温曲线图的绘制方法、系统、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101796507A (zh) * 2007-08-27 2010-08-04 兰德马克绘图国际公司，哈里伯顿公司 计算变差函数模型的系统和方法
CN102360404A (zh) * 2011-11-01 2012-02-22 中国水电顾问集团华东勘测设计研究院 一种洞室地质展示图的编绘方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101796507A (zh) * 2007-08-27 2010-08-04 兰德马克绘图国际公司，哈里伯顿公司 计算变差函数模型的系统和方法
EA020532B1 (ru) * 2007-08-27 2014-11-28 Лэндмарк Грэфикс Корпорейшн, Э Хэллибертон Кампани Способ и машиночитаемый носитель для формирования вариограммной модели свойств пласта-коллектора
CN102360404A (zh) * 2011-11-01 2012-02-22 中国水电顾问集团华东勘测设计研究院 一种洞室地质展示图的编绘方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈军: "用Excel制节理走向玫瑰花图", 《内蒙古水利》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874251A (zh) * 2017-03-10 2017-06-20 中国电建集团成都勘测设计研究院有限公司 一种岩体rbi和rqd指标的测量方法
CN107389555A (zh) * 2017-07-21 2017-11-24 中水北方勘测设计研究有限责任公司 一种分析钻孔岩体主要节理发育方向的方法
CN107590849A (zh) * 2017-08-09 2018-01-16 山西晋城无烟煤矿业集团有限责任公司 一种钻孔自动成图方法
CN109637324A (zh) * 2018-12-19 2019-04-16 西南石油大学 一种原始纹层产状恢复的方法
CN111063035A (zh) * 2019-11-22 2020-04-24 同济大学 一种gis中od关系的三维可视化方法及装置
CN111063035B (zh) * 2019-11-22 2021-09-03 同济大学 一种gis中od关系的三维可视化方法及装置
CN111159862A (zh) * 2019-12-17 2020-05-15 核工业二四三大队 一种铀矿钻孔设计指示书电子化设计及编制方法
CN111159862B (zh) * 2019-12-17 2023-09-08 核工业二四三大队 一种铀矿钻孔设计指示书电子化设计及编制方法
CN112767515A (zh) * 2021-02-19 2021-05-07 上海市政工程设计研究总院（集团）有限公司 测试点地温曲线图的绘制方法、系统、设备及介质

Similar Documents

Publication Publication Date Title
CN105955939A (zh) 一种基于Excel绘制地质体及地质要素倾斜玫瑰花图的方法
CN103136371B (zh) 多源空间数据的剖分标识生成方法以及数据检索方法
CN103425772B (zh) 一种具有多维信息的海量数据查询方法
CN112001016B (zh) 一种基于元模型和建构逻辑的建筑信息模型自动建模方法
CN111611664B (zh) 一种输电线路路径选择方法及系统
CN107564087A (zh) 一种基于屏幕的三维线状符号渲染方法
CN106568377B (zh) 基于ArcEngine的构造大地测量成图方法及系统
CN112148774B (zh) 一种高分辨率空间多区域气象数据处理系统及方法
CN101882026A (zh) 笔迹书写的方法和装置、笔迹模板分类器训练方法和装置
CN103530904A (zh) 一种基于克里金方法的水下地形数字高程建立方法
CN106528942A (zh) 一种基于中子输运特征线理论的数据可视化方法
CN110473251A (zh) 基于网格空间索引的自定义范围空间数据面积统计方法
CN106777302A (zh) 空间地理坐标的转换方法及装置
CN103198185A (zh) 一种绘制塔基断面图的方法和装置
CN116956838A (zh) 工程量清单的项目确定方法及确定装置
CN107330233A (zh) 一种输电塔设计风速分析方法和装置
CN109655044B (zh) 一种输电线路勘测定位测量成品质量控制方法
CN104331389A (zh) 基于八点法的等值线追踪算法
CN113095012B (zh) 风电场流场分区数值仿真计算结果拼接融合方法
CN107464272A (zh) 基于关键点的中心扩散式气象传真图等值线的插值方法
CN113656852B (zh) 一种精细化河道地形快速生成方法
CN206249653U (zh) 一种多功能练字田宫格
CN111090714B (zh) 一种基于cass地形图的数据入库方法和装置
CN111241221A (zh) 一种受损地形坐标数据的自动化匹配及高精度修复方法
CN108090962B (zh) 基于ansys数据文件的plaxis几何模型快速生成方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160921