CN106568377B - 基于ArcEngine的构造大地测量成图方法及系统 - Google Patents
基于ArcEngine的构造大地测量成图方法及系统 Download PDFInfo
- Publication number
- CN106568377B CN106568377B CN201610975859.3A CN201610975859A CN106568377B CN 106568377 B CN106568377 B CN 106568377B CN 201610975859 A CN201610975859 A CN 201610975859A CN 106568377 B CN106568377 B CN 106568377B
- Authority
- CN
- China
- Prior art keywords
- arcengine
- circular arc
- deformation
- angle
- arc
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 230000007246 mechanism Effects 0.000 claims abstract description 57
- 238000012800 visualization Methods 0.000 claims abstract description 15
- 238000010276 construction Methods 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 21
- 230000009471 action Effects 0.000 claims description 16
- 101100537937 Caenorhabditis elegans arc-1 gene Proteins 0.000 claims description 12
- 230000000295 complement effect Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000005012 migration Effects 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims 2
- 238000006073 displacement reaction Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000008520 organization Effects 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 102100037922 Disco-interacting protein 2 homolog A Human genes 0.000 description 1
- 101000805876 Homo sapiens Disco-interacting protein 2 homolog A Proteins 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Navigation (AREA)
- Instructional Devices (AREA)
Abstract
本发明提供一种基于ArcEngine的构造大地测量成图方法，借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；设置成图参数后矢量输出成图。本发明基于ArcEngine实现了构造大地测量成图的可视化，避免了使用命令行反复调整参数，提高了绘图效率。
Description
技术领域
本发明属于专题制图领域，具体涉及一种基于ArcEngine的构造大地测量成图方法及系统。
背景技术
在科学研究中分析结果的表达相当重要，结果表达的正确清楚与否直接影响读者对于分析结果的理解。通常情况下，地理图是地学工作者表达结果的重要方式，通过地图表达读者可以直观了解结果并进行主观评估。在构造大地测量学中，如何正确、生动、美观地表达构造地貌、形变场和震源机制解信息，是地学学者在科研过程中不可避免且具有重要意义的工作。
地学研究中常采用GMT(Generic Mapping tools)软件绘制各类地图，GMT是一个被地理学界广泛使用的绘图工具，可以绘制海岸线、国界、河流作为底图，调用资料库中的地理信息，表达用户输入的各种参数，比如颜色配置、边框文字配置，并结合科研表达的结果，如经度、纬度、形变大小，生成ps格式的地图文件。然而，GMT是个命令行工具，无法可视化绘图过程，用户很难做到一次性成图。为了获取理想的输出成图，用户通常需要多次调整绘图参数，时间成本较高。另外，GMT作为命令行工具，其学习成本高于可视化绘图工具。
发明内容
本发明要解决的技术问题是：提供一种基于ArcEngine的构造大地测量成图方法及系统，解决了GMT软件绘制构造地貌、形变场和震源机制解时无法可视化的问题。
本发明为解决上述技术问题所采取的技术方案为：一种基于ArcEngine的构造大地测量成图方法，其特征在于：它包括以下步骤：
构造地貌：借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；
绘制GPS形变场、水准形变场和震源机制解：
根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；
根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；
设置成图参数：
在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
矢量输出成图。
按上述方法，它还包括以下步骤：
建立矢量图形库：
将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储；
所述的设置成图参数还包括：调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
按上述方法，所述的GPS形变场和水准形变场的绘制包括以下步骤：
根据测站经纬度，确定水平形变或水准形变的起点；
在起点的经度方向和纬度方向分别加上两个方向上形变量与缩放因子之积，确定为水平形变的终点；在起点纬度方向加上水准形变量与缩放因子之积，确定为水准形变的终点；
使用ArcEngine的Polyline类组合水平形变的起点和终点、以及水准形变的起点和终点，并由制图线符号绘制用于表示形变的大小和方向的箭头，完成形变部分的绘制；
根据GPS测站经度向形变量和纬度向形变量的中误差及其相关系数，以水平形变终点作为误差椭圆中心，计算误差椭圆三要素：长轴、短轴和旋转角；然后使用ArcEngine的EllipticArc类构建缩放因子作用下的椭圆弧，并根据旋转角使用ITransform接口旋转椭圆弧；最后使用AddSegment方法转换得到大量椭圆点，以Polyline类连接椭圆点，完成GPS形变场误差椭圆的绘制；误差椭圆表示水平形变误差的大小；所述的缩放因子代表图形放大或缩小的程度，由用户设置。
按上述方法，所述的震源机制解的绘制包括以下步骤：
以地震矩震级乘以0.01倍的缩放因子，确定震源机制解的绘制半径；
以震中位置为圆心，震源机制解绘制半径为半径，使用ArcEngine中的ConstructCircle方法分别构建外圆圈和内圆圈，并借助Ring类填充外圆圈内部为白色；
将所述的两个P波初动节面记为节面1和节面2；
根据节面1的走向方位角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第一个走向圆弧的起点p1和终点p2，并以震源机制解绘制半径乘以节面1倾角的余弦值，得到第一个走向圆弧的弦高；根据起点p1、终点p2和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，在内圆圈内顺时针构建第一个走向圆弧，记为StrikeArc1；
根据节面2的走向方位角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第二个走向圆弧的起点p3和终点p4；在内圈圈上顺时针截取点p2到点p3的圆弧，记为Arc1，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若是则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；此外，进一步判断节面1的滑动角是否为负，若是则继续选取已有圆弧对应的另一部分弧段，由此保证节面1滑动角为正时圆弧Arc1为优弧；
以震源机制解绘制半径乘以节面2倾角的余弦值，得到第二个走向圆弧的弦高，根据该走向圆弧的起点p3、终点p4和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，顺时针构建第二个走向圆弧，记为StrikeArc2；
在内圈圈上顺时针截取点p4到点p1的圆弧，记为Arc2，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若否则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；
使用ArcEngine中的Ring类组合第一个走向圆弧StrikeArc1、圆弧Arc1、第二个走向圆弧StrikeArc2和、圆弧Arc2，并根据断层类型在构建的封闭区域内填充相应颜色。
按上述方法，所述的断层类型根据节面1滑动角的大小确定，当滑动角的绝对值小于等于45度或大于等于135度时，设定为走滑断层；当滑动角大于-135度并且小于-45度时，设定为正断层；当滑动角大于45度并且小于135度时，设定为逆冲断层。
按上述方法，所述的设置成图参数包括增加色带，用于在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。
按上述方法，所述的设置成图参数包括设置地图显示范围，由用户根据成图需要自行设置地图左上角点和右下角点的经纬度，确定地图显示范围。
一种基于ArcEngine的构造大地测量成图系统，其特征在于：它包括：
构造地貌模块，用于借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；
GPS形变场和水准形变场绘制模块，用于根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；
震源机制解绘制模块，用于根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；
成图参数设置模块，用于在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
成图模块，用于在设置完成图参数后，矢量输出成图。
按上述系统，它还包括矢量图形库建立模块，用于将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储在矢量图形库中；
所述的参数设置模块还用于调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
按上述系统，所述的参数设置模块还用于增加色带，在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。
本发明的有益效果为：
1、基于ArcEngine实现了构造大地测量成图的可视化：用户在绘制构造地貌、GPS形变场、水准形变场和震源机制解时，可以直接在地图上调整图形，避免了使用命令行反复调整参数，提高了绘图效率。
2、基于ArcEngine实现了构造大地测量成图的图层化：通过建立矢量图形库，将GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储，方便在ArcGIS中作为图层调用，在下一次绘图时只需要调用已有的矢量图形文件并设置一些绘制数据即可，无需二次绘制，进一步提高了绘图效率。
3、本发明的构造大地测量绘图操作更为简便直观。用户只需按格式准备数据文件，在交互界面中选择绘制数据，并设置少许必要参数后，即可快速获取高质量成图。
附图说明
图1为本发明一实施例的方法流程图。
图2为本发明一实施例的构造地貌的成图。
图3为本发明一实施例的GPS形变场的成图。
图4为本发明一实施例的水准形变场的成图。
图5为本发明一实施例的震源机制解的成图。
具体实施方式
下面结合具体实例和附图对本发明做进一步说明。
本发明提供一种基于ArcEngine的构造大地测量成图方法，如图1所示，它包括以下步骤：
S1、构造地貌：借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌。
S2、绘制GPS形变场、水准形变场和震源机制解：
2.1、根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场。
GPS形变场和水准形变场的绘制包括以下步骤：
2.1.1、根据测站经纬度，确定水平形变或水准形变的起点；
2.1.2、在起点的经度方向和纬度方向分别加上两个方向上形变量与缩放因子之积，确定为水平形变的终点；在起点纬度方向加上水准形变量与缩放因子之积，确定为水准形变的终点；
2.1.3、使用ArcEngine的Polyline类组合水平形变的起点和终点、以及水准形变的起点和终点，并由制图线符号绘制用于表示形变的大小和方向的箭头，完成形变部分的绘制；
2.1.4、根据GPS测站经度向形变量和纬度向形变量的中误差及其相关系数，以水平形变终点作为误差椭圆中心，计算误差椭圆三要素：长轴、短轴和旋转角；然后使用ArcEngine的EllipticArc类构建缩放因子作用下的椭圆弧，并根据旋转角使用ITransform接口旋转椭圆弧；最后使用AddSegment方法转换得到大量椭圆点，以Polyline类连接椭圆点，完成GPS形变场误差椭圆的绘制；误差椭圆表示水平形变误差的大小；所述的缩放因子代表图形放大或缩小的程度，由用户设置。
2.2、根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解。
震源机制解的绘制包括以下步骤：
2.2.1、以地震矩震级乘以0.01倍的缩放因子，确定震源机制解的绘制半径；
2.2.2、以震中位置为圆心，震源机制解绘制半径为半径，使用ArcEngine中的ConstructCircle方法分别构建外圆圈和内圆圈，并借助Ring类填充外圆圈内部为白色；
将所述的两个P波初动节面记为节面1和节面2；
2.2.3、根据节面1的走向方位角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第一个走向圆弧的起点p1和终点p2，并以震源机制解绘制半径乘以节面1倾角的余弦值，得到第一个走向圆弧的弦高；根据起点p1、终点p2和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，在内圆圈内顺时针构建第一个走向圆弧，记为StrikeArc1；
2.2.4、根据节面2的走向方位角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第二个走向圆弧的起点p3和终点p4；在内圈圈上顺时针截取点p2到点p3的圆弧，记为Arc1，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若是则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；此外，进一步判断节面1的滑动角是否为负，若是则继续选取已有圆弧对应的另一部分弧段，由此保证节面1滑动角为正时圆弧Arc1为优弧；
2.2.5、以震源机制解绘制半径乘以节面2倾角的余弦值，得到第二个走向圆弧的弦高，根据该走向圆弧的起点p3、终点p4和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，顺时针构建第二个走向圆弧，记为StrikeArc2；
2.2.6、在内圈圈上顺时针截取点p4到点p1的圆弧，记为Arc2，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若否则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段
2.2.7、使用ArcEngine中的Ring类组合第一个走向圆弧StrikeArc1、圆弧Arc1、第二个走向圆弧StrikeArc2和、圆弧Arc2，并根据断层类型在构建的封闭区域内填充相应颜色。其中，所述的断层类型根据节面1滑动角的大小确定，当滑动角的绝对值小于等于45度或大于等于135度时，设定为走滑断层；当滑动角大于-135度并且小于-45度时，设定为正断层；当滑动角大于45度并且小于135度时，设定为逆冲断层。
S3、设置成图参数：在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
设置成图参数包括：1)增加色带，用于在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。2)设置地图显示范围，由用户根据成图需要自行设置地图左上角点和右下角点的经纬度，确定地图显示范围。
S4、矢量输出成图。
优选的，它还包括以下步骤：
建立矢量图形库：将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件(.shp)并存储；shpfile字段包括ID、测站名“Name”、时间“Datetime”和几何形状“Shape”。
所述的设置成图参数还包括：调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
一种基于ArcEngine的构造大地测量成图系统，它包括：
构造地貌模块，用于借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；
GPS形变场和水准形变场绘制模块，用于根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；
震源机制解绘制模块，用于根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；
成图参数设置模块，用于在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
成图模块，用于在设置完成图参数后，矢量输出成图。
优选的，它还包括矢量图形库建立模块，用于将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储在矢量图形库中；
所述的参数设置模块还用于调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
所述的参数设置模块还用于增加色带，在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。所述的参数设置模块还用于设置地图显示范围，由用户根据成图需要自行设置地图左上角点和右下角点的经纬度，确定地图显示范围。
下面以2015年Nepal地震为例，该地震发生在青藏高原和印度板块交界地区。
1、构造地貌的绘制成图
本绘图方法采用的坐标系统为WGS-84。在获取该区域的数字高程模型后，借助ArcGIS图层功能可以将其作为底图绘制在地图上，成图结果如图2所示。
2、GPS形变场的绘制成图
在示例区域解算GPS测站形变，将解算结果存储为.gps文件，为方便GPS形变场的绘制，自定义文件格式见表1，其中Site列为测站名，Long(deg)列为测站经度(单位为度)，Lat(deg)列为测站纬度(单位为度)，H(m)列为测站高度(单位为米)，DisE(cm)列为东方向位移(单位为厘米)，DisN(cm)列为北方向位移(单位为厘米)，DisH(cm)列为垂直方向位移(单位为厘米)，SigE(cm)列为东方向位移中误差(单位为厘米)，SigN(cm)列为北方向位移中误差(单位为厘米)，SigH(cm)列为垂直方向位移中误差(单位为厘米)，Corr列为水平方向相关系数，Datetime列为测站观测时间(年/月/日，为方便显示，该列未列入表中)。表1同时给出用于GPS形变场绘制的实施例部分数据。
表1：实施例中GPS形变场绘制过程读取的部分形变数据
用户在GPS形变文件组织完成后，只需要在绘制GPS形变场对话框中选择该文件作为输入数据，并设定图形的缩放因子保证显示效果，点击“绘制”按钮即可使得程序执行技术方案中GPS形变场绘制的步骤。带有误差椭圆的矢量GPS形变场显示在数字高程模型底图上，输出成图如图3所示。
3、水准形变场的绘制成图
在示例区域解算测站水准形变，将解算结果存储为.lev文件，为方便水准形变场的绘制，自定义文件格式见表2，其中Site列为测站名，Long(deg)列为测站经度(单位为度)，Lat(deg)列为测站纬度(单位为度)，H(m)列为测站高度(单位为米)，DisH(cm)列为垂直方向位移(单位为厘米)，SigH(cm)列为垂直方向位移中误差(单位为厘米)，Datetime列为测站观测时间(年/月/日)。表2同时给出用于水准形变场绘制的实施例部分数据。
表2：实施例中水准形变场绘制过程读取的部分形变数据
用户在水准形变文件组织完成后，只需要在绘制水准形变场对话框中选择该文件作为输入数据，并设定图形的缩放因子保证显示效果，点击“绘制”按钮即可使得程序执行技术方案中水准形变场绘制的步骤，由于水准形变场只表达于垂直方向，且不存在误差椭圆，因此其绘制过程相比GPS形变场的绘制更为简单。绘制完成后，水准形变场显示在数字高程模型底图上，输出成图如图4所示。
4、震源机制解的绘制成图
在示例区域解算震源机制解，将解算结果存储为.fms文件，为方便震源机制解的绘制，自定义文件格式见表4，其中Name列为地震名称，Long(deg)列为震中经度(单位为度)，Depth(km)列为断层破裂深度(单位为千米)，Mw列为地震矩震级，Strike1(deg)列为节面1的走向角(单位为度)，Dip1(deg)列为节面1的倾角(单位为度)，Rake1(deg)列为节面1的滑动角(单位为度)，Strike2(deg)列为节面2的走向角(单位为度)，Dip2(deg)列为节面2的倾角(单位为度)，Rake2(deg)列为节面2的滑动角(单位为度)，Datetime列为发震时间(年/月/日时：分：秒，为方便显示，该列未列入表中)。表3同时给出用于震源机制解绘制的实施例部分数据。
表3：实施例中震源机制解绘制过程读取的部分数据
用户在震源机制解文件组织完成后，需要在绘制震源机制解对话框中选择该文件作为输入数据，并设定图形的缩放因子保证显示效果，点击“绘制”按钮即可使得程序执行技术方案中震源机制解绘制的步骤。绘制过程中ICircularArc接口的Complement方法非常重要，它的功能是选取已有圆弧对应的另一部分弧段。借助该功能本方法保证了震源机制解中颜色填充的正确性，由此正确区分了压缩区和拉张区。另外，在绘图对话框中给出了用户断层类型选项：走滑断层、正断层和逆冲断层，用户可根据显示需求选择性绘制特定断层类型的震源机制解。本实施例中选择绘制所有断层类型的震源机制解，结果显示在数字高程模型底图上，输出成图如图5所示。
5、矢量图形文件的转换与存储
在绘制完成GPS形变场、水准形变场或震源机制解后，除了将其输出成图，还可以转换为矢量图形文件.shp并存储，方便后续工作的调用和即时修改。用户只需在矢量格式转换对话框中选择形变数据或震源机制解数据，并设置矢量文件输出路径，点击“转换”按钮即可完成转换和存储过程。
以上实施例仅用于说明本发明的设计思想和特点，其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施，本发明的保护范围不限于上述实施例。所以，凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰，均在本发明的保护范围之内。
Claims (10)
1.一种基于ArcEngine的构造大地测量成图方法，其特征在于：它包括以下步骤：
构造地貌：借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；
绘制GPS形变场、水准形变场和震源机制解：
根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；
根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；
设置成图参数：
在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
矢量输出成图；
所述的震源机制解的绘制包括以下步骤：
设置缩放因子，以地震矩阵级乘以缩放因子，确定震源机制解的绘制半径；
以震中位置为圆心，震源机制解绘制半径为半径，使用ArcEngine中的ConstructCircle方法分别构建外圆圈和内圆圈，并借助Ring类填充外圆圈内部为白色；
将所述的两个P波初动节面记为节面1和节面2；
根据节面1的走向角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第一个走向圆弧的起点p1和终点p2，并以震源机制解绘制半径乘以节面1倾角的余弦值，得到第一个走向圆弧的弦高；根据起点p1、终点p2和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，在内圆圈内顺时针构建第一个走向圆弧，记为StrikeArc1；
根据节面2的走向角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第二个走向圆弧的起点p3和终点p4；在内圈圈上顺时针截取点p2到点p3的圆弧，记为Arc1，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若是则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；此外，进一步判断节面1的滑动角是否为负，若是则继续选取已有圆弧对应的另一部分弧段，由此保证节面1的滑动角为正时，圆弧Arc1为优弧；
以震源机制解绘制半径乘以节面2倾角的余弦值，得到第二个走向圆弧的弦高，根据该走向圆弧的起点p3、终点p4和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，顺时针构建第二个走向圆弧，记为StrikeArc2；
在内圈圈上顺时针截取点p4到点p1的圆弧，记为Arc2，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若否则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；
使用ArcEngine中的Ring类组合第一个走向圆弧StrikeArc1、圆弧Arc1、第二个走向圆弧StrikeArc2和圆弧Arc2，并根据断层类型在构建的封闭区域内填充相应颜色。
2.根据权利要求1所述的基于ArcEngine的构造大地测量成图方法，其特征在于：它还包括以下步骤：
建立矢量图形库：
将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储；
所述的设置成图参数还包括：调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
3.根据权利要求1或2所述的基于ArcEngine的构造大地测量成图方法，其特征在于：所述的GPS形变场和水准形变场的绘制包括以下步骤：
根据测站经纬度，确定水平形变或水准形变的起点；
在起点的经度方向和纬度方向分别加上两个方向上形变量与缩放因子之积，确定为水平形变的终点；在起点纬度方向加上水准形变量与缩放因子之积，确定为水准形变的终点；
使用ArcEngine的Polyline类组合水平形变的起点和终点、以及水准形变的起点和终点，并由制图线符号绘制用于表示形变的大小和方向的箭头，完成形变部分的绘制；
根据GPS测站经度向形变量和纬度向形变量的中误差及其相关系数，以水平形变终点作为误差椭圆中心，计算误差椭圆三要素：长轴、短轴和旋转角；然后使用ArcEngine的EllipticArc类构建缩放因子作用下的椭圆弧，并根据旋转角使用ITransform接口旋转椭圆弧；最后使用AddSegment方法转换得到大量椭圆点，以Polyline类连接椭圆点，完成GPS形变场误差椭圆的绘制；误差椭圆表示水平形变误差的大小；所述的缩放因子代表图形放大或缩小的程度，由用户设置。
4.根据权利要求1或2所述的基于ArcEngine的构造大地测量成图方法，其特征在于：所述的缩放因子为以地震矩震级乘以0.01倍。
5.根据权利要求4所述的基于ArcEngine的构造大地测量成图方法，其特征在于：所述的断层类型根据节面1滑动角的大小确定，当滑动角的绝对值小于等于45度或大于等于135度时，设定为走滑断层；当滑动角大于-135度并且小于-45度时，设定为正断层；当滑动角大于45度并且小于135度时，设定为逆冲断层。
6.根据权利要求1或2所述的基于ArcEngine的构造大地测量成图方法，其特征在于：所述的设置成图参数包括增加色带，用于在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。
7.根据权利要求1或2所述的基于ArcEngine的构造大地测量成图方法，其特征在于：所述的设置成图参数包括设置地图显示范围，由用户根据成图需要自行设置地图左上角点和右下角点的经纬度，确定地图显示范围。
8.一种基于ArcEngine的构造大地测量成图系统，其特征在于：它包括：
构造地貌模块，用于借助ArcGIS图层功能，将构造所在区域的数字高程模型作为底图加载在地图上，用于表示构造地区的地貌；
GPS形变场和水准形变场绘制模块，用于根据测站位置、水平形变值、水准形变值和形变中误差信息，于测站所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在测站位置绘制GPS形变场和水准形变场；
震源机制解绘制模块，用于根据震中位置、地震矩震级、两个P波初动节面的走向角、倾角和滑动角，于震中所在区域的数字高程模型上，使用ArcEngine调用三维地理信息系统可视化表达功能，在震中位置绘制震源机制解；
成图参数设置模块，用于在绘制完成后，基于ArcGIS的布局视图，设置成图参数；
成图模块，用于在设置完成图参数后，矢量输出成图；
所述的震源机制解的绘制包括以下步骤：
设置缩放因子，以地震矩阵级乘以缩放因子，确定震源机制解的绘制半径；
以震中位置为圆心，震源机制解绘制半径为半径，使用ArcEngine中的ConstructCircle方法分别构建外圆圈和内圆圈，并借助Ring类填充外圆圈内部为白色；
将所述的两个P波初动节面记为节面1和节面2；
根据节面1的走向角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第一个走向圆弧的起点p1和终点p2，并以震源机制解绘制半径乘以节面1倾角的余弦值，得到第一个走向圆弧的弦高；根据起点p1、终点p2和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，在内圆圈内顺时针构建第一个走向圆弧，记为StrikeArc1；
根据节面2的走向角，使用ArcEngine中的QueryPoint方法在内圆圈上寻找第二个走向圆弧的起点p3和终点p4；在内圈圈上顺时针截取点p2到点p3的圆弧，记为Arc1，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若是则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；此外，进一步判断节面1的滑动角是否为负，若是则继续选取已有圆弧对应的另一部分弧段，由此保证节面1的滑动角为正时，圆弧Arc1为优弧；
以震源机制解绘制半径乘以节面2倾角的余弦值，得到第二个走向圆弧的弦高，根据该走向圆弧的起点p3、终点p4和弦高，使用ArcEngine中的ConstructEndPointChordHeight方法，顺时针构建第二个走向圆弧，记为StrikeArc2；
在内圈圈上顺时针截取点p4到点p1的圆弧，记为Arc2，借助ICircularArc接口的IsMinor方法判断该圆弧是否为劣弧，若否则使用ICircularArc接口的Complement方法选取已有圆弧对应的另一部分弧段；
使用ArcEngine中的Ring类组合第一个走向圆弧StrikeArc1、圆弧Arc1、第二个走向圆弧StrikeArc2和圆弧Arc2，并根据断层类型在构建的封闭区域内填充相应颜色。
9.根据权利要求8所述的基于ArcEngine的构造大地测量成图系统，其特征在于：它还包括矢量图形库建立模块，用于将绘制好的GPS形变场、水准形变场和震源机制解转换为矢量图形文件并存储在矢量图形库中；
所述的参数设置模块还用于调用矢量图形库中的矢量图形文件，基于ArcGIS的布局视图，设置成图参数。
10.根据权利要求8或9所述的基于ArcEngine的构造大地测量成图系统，其特征在于：所述的参数设置模块还用于增加色带，在使用Polygon类在地图特定位置绘制的色带框中，借助PolygonElement类填充由GradientFillSymbol类构建的颜色带；颜色带的类型和颜色由用户设置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610975859.3A CN106568377B (zh) | 2016-11-07 | 2016-11-07 | 基于ArcEngine的构造大地测量成图方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610975859.3A CN106568377B (zh) | 2016-11-07 | 2016-11-07 | 基于ArcEngine的构造大地测量成图方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106568377A CN106568377A (zh) | 2017-04-19 |
| CN106568377B true CN106568377B (zh) | 2019-07-23 |
Family
ID=58540061
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610975859.3A Expired - Fee Related CN106568377B (zh) | 2016-11-07 | 2016-11-07 | 基于ArcEngine的构造大地测量成图方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106568377B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108009217A (zh) * | 2017-11-20 | 2018-05-08 | 北京林业大学 | 基于数据库驱动快速批量制作中国水资源保障风险图的方法 |
| CN107844607B (zh) * | 2017-12-05 | 2020-09-01 | 天地科技股份有限公司 | 一种基于gis的地应力展示系统 |
| CN111580728B (zh) * | 2020-04-07 | 2021-07-20 | 深圳震有科技股份有限公司 | 基于状态模式的多种军标类型动态标绘方法及装置 |
| CN111580160A (zh) * | 2020-04-28 | 2020-08-25 | 中国地震台网中心 | 地震信息的显示方法及装置 |
| KR102307615B1 (ko) * | 2021-07-02 | 2021-10-01 | 주식회사 지오앤 | 지리좌표를 이용한 측지측량데이터의 합성을 통해 지상구조물의 위치를 확인할 수 있는 측지측량시스템 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1838073A (zh) * | 2006-04-25 | 2006-09-27 | 武汉大学 | 基于虚拟机的地理信息系统通用符号系统及其构建方法 |
| CN101545982A (zh) * | 2009-05-05 | 2009-09-30 | 江西省电力信息通讯有限公司 | 地震计算模型的动态拟合方法 |
| CN103412113A (zh) * | 2013-08-26 | 2013-11-27 | 成都理工大学 | 一种震后泥石流沟敏感性判别方法及其应用 |
| CN104794236A (zh) * | 2015-05-07 | 2015-07-22 | 中国测绘科学研究院 | 一种地图制图规则构建及其结构化组织方法与系统 |
| CN105005580A (zh) * | 2015-06-05 | 2015-10-28 | 北京世纪安图数码科技发展有限责任公司 | 一种用于显示水库地形的方法及其装置 |
| CN105809616A (zh) * | 2016-04-18 | 2016-07-27 | 山东科技大学 | 基于ArcGIS的1：1万国家标准地形图的制作方法 |
| CN106461805A (zh) * | 2014-04-30 | 2017-02-22 | 雪佛龙美国公司 | 创建和解释多个地震勘测的动画拼接图的方法 |
| CN104200126B (zh) * | 2014-09-23 | 2017-04-12 | 北京师范大学 | 一种地震导致滑坡所造成的人员死亡的快速评估方法 |
-
2016
- 2016-11-07 CN CN201610975859.3A patent/CN106568377B/zh not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1838073A (zh) * | 2006-04-25 | 2006-09-27 | 武汉大学 | 基于虚拟机的地理信息系统通用符号系统及其构建方法 |
| CN101545982A (zh) * | 2009-05-05 | 2009-09-30 | 江西省电力信息通讯有限公司 | 地震计算模型的动态拟合方法 |
| CN103412113A (zh) * | 2013-08-26 | 2013-11-27 | 成都理工大学 | 一种震后泥石流沟敏感性判别方法及其应用 |
| CN106461805A (zh) * | 2014-04-30 | 2017-02-22 | 雪佛龙美国公司 | 创建和解释多个地震勘测的动画拼接图的方法 |
| CN104200126B (zh) * | 2014-09-23 | 2017-04-12 | 北京师范大学 | 一种地震导致滑坡所造成的人员死亡的快速评估方法 |
| CN104794236A (zh) * | 2015-05-07 | 2015-07-22 | 中国测绘科学研究院 | 一种地图制图规则构建及其结构化组织方法与系统 |
| CN105005580A (zh) * | 2015-06-05 | 2015-10-28 | 北京世纪安图数码科技发展有限责任公司 | 一种用于显示水库地形的方法及其装置 |
| CN105809616A (zh) * | 2016-04-18 | 2016-07-27 | 山东科技大学 | 基于ArcGIS的1：1万国家标准地形图的制作方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106568377A (zh) | 2017-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106568377B (zh) | 基于ArcEngine的构造大地测量成图方法及系统 | |
| US11042278B2 (en) | Providing graphical indication of label boundaries in digital maps | |
| CN105516311A (zh) | 一种电子地图全景图获取方法及系统 | |
| JP7278173B2 (ja) | 地図表示装置、方法およびプログラム | |
| CN107729673A (zh) | 基于bim的道路桥梁实景模型分析方法、装置及其施工方法 | |
| US11373256B2 (en) | Editing and representing property values for a digital map | |
| CN102789361B (zh) | 一种矢量电子地图信息的显示风格设置方法 | |
| CN110597940A (zh) | 基于微信的自然资源外业核查方法及系统、服务器及介质 | |
| CN103971412B (zh) | 电子地图的展现方法及装置 | |
| Jenny et al. | Design principles for Swiss-style rock drawing | |
| Loeffler et al. | Mobile UX design: learning from the Flyover Country mobile app | |
| CN104199887B (zh) | 基于地理信息系统与Delft3D‑FLOW计算引擎的水环境分析方法及系统 | |
| CN112509453A (zh) | 基于移动设备的景区实景导览图电子化导览方法及系统 | |
| Brainerd et al. | Interactive map projections and distortion | |
| CN112099781A (zh) | 一种地图可视化方法、装置、存储介质及设备 | |
| JP6533075B2 (ja) | 地図画像表示装置、ナビゲーション装置、地図画像表示プログラム、及び地図画像表示方法 | |
| Dobesova et al. | Comparison of visual languages in Geographic Information Systems | |
| Ahmed et al. | Introduction to basic GIS and spatial analysis using QGIS: Applications in Bangladesh | |
| KR20200051295A (ko) | 분양 정보를 포함한 지도를 제공하는 방법 및 컴퓨터 프로그램 | |
| JP2012063469A (ja) | 地図表示装置、地図表示方法およびプログラム | |
| JP5735939B2 (ja) | 地図表示装置、地図表示方法及び地図表示プログラム | |
| CN107967286A (zh) | 一种地质上精确定位野外复杂地形路线的方法 | |
| JP2015036694A (ja) | 地図注記処理装置、地図注記処理方法および地図情報提供装置 | |
| Promneewat et al. | Offline Interactive Map from Hybrid App Development: A Case from Geologic Map App | |
| Islam et al. | Mastering Geospatial Development with QGIS 3. x: An in-depth guide to becoming proficient in spatial data analysis using QGIS 3.4 and 3.6 with Python |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190723 Termination date: 20211107 |