pygplates专栏——Reconstruc features——reconstruct regular features

pygplates专栏——Reconstruc features——reconstruct regular features

  • Reconstruct regular features
    • 导出重构特征到文件
      • 示例代码
      • 详解
    • 计算重建距离
      • 示例代码
      • 详解

Reconstruct regular features

这个例子展示了几个不同的场景,涉及到重建地质时代的常见特征。

导出重构特征到文件

在这个例子中,我们重建常规特征并将结果导出到Shapefile。

示例代码

import pygplates# 加载板块运动模型
rotation_model = pygplates.RotationModel("Muller2019-Young2019-Cao2020_CombinedRotations.rot")
# 加载一些features
features = pygplates.FeatureCollection("Global_EarthByte_GPlates_PresentDay_Coastlines.gpmlz")
# 重建的地质时间
reconstruction_time = 50
# 输出的文件
export_filename = "6-Exported_reconstructed_features_to_a_file_reconstructed_{0}Ma.shp".format(reconstruction_time)
# 重建
pygplates.reconstruct(features, rotation_model, export_filename, reconstruction_time)

详解

首先加载板块运动模型(pygplates.RotationModel)

rotation_model = pygplates.RotationModel('rotations.rot')

加载需要重建的特征(pygplates.FeatureCollection)

features = pygplates.FeatureCollection('features.gpml')

设置重建时间

reconstruction_time = 50

最后完成重建结果并保存

pygplates.reconstruct(features, rotation_model, export_filename, reconstruction_time)

计算重建距离

示例代码

import pygplates# 返回几何要素(点/多点/折线/多边形)质心的函数。
def get_geometry_centroid(geometry):# 检查几何要素是否为多边形try:return geometry.get_interior_centroid()except AttributeError:# 不是多边形,继续pass# 检查几何要素是否为多线条或多点try:return geometry.get_centroid()except AttributeError:pass# 仅剩点return geometry
# 加载板块运动模型
rotation_model = pygplates.RotationModel("Muller2019-Young2019-Cao2020_CombinedRotations.rot")
# 加载一些特征
features = pygplates.FeatureCollection("2-output_points.gpml")
# 重建地质时间
reconstruction_time = 50
# 重建
reconstructed_feature_geometries = []
pygplates.reconstruct(features, rotation_model, reconstructed_feature_geometries, reconstruction_time)
# 遍历所有重建结果
for reconstructed_feature_geometry in reconstructed_feature_geometries:# 计算距离:今时今日几何要素的质点和重建后的质点distance_reconstructed = pygplates.GeometryOnSphere.distance(get_geometry_centroid(reconstructed_feature_geometry.get_present_day_geometry()),get_geometry_centroid(reconstructed_feature_geometry.get_reconstructed_geometry()))# 将弧度转换为公里distance_reconstructed_in_kms = distance_reconstructed * pygplates.Earth.mean_radius_in_kms# 输出相关特征名称和板块ID,以及重建的距离print("Feature: %s" % reconstructed_feature_geometry.get_feature().get_name())print("  plate ID: %d" % reconstructed_feature_geometry.get_feature().get_reconstruction_plate_id())print("  distance reconstructed: %f Kms" % distance_reconstructed_in_kms)

详解

为了某个功能,pygplates定义了许多相关的函数
如果,我们可以忽略掉几何要素是否为pygplates.PointOnSphere,pygplates.MultiPointOnSphere,pygplates.PolyineOnSphere或pygplates.PolygonOnSphere。每种几何要素类型需要不同的方法获取质点。
因此我们首先尝试pygplates.PolygonOnSphere.get_interior_centroid(),然后尝试get_centroid(),最后几何要素点就是它本身的质点。

def get_geometry_centroid(geometry):try:return geometry.get_interior_centroid()except AttributeError:passtry:return geometry.get_centroid()except AttributeError:passreturn geometry

使用pygplates.GeometryOnSphere.distance()方法来计算两个质点之间的最短距离,然后使用pygplates.Earth将其转换为公里数。

distance_reconstructed = pygplates.GeometryOnSphere.distance(get_geometry_centroid(reconstructed_feature_geometry.get_present_day_geometry()),get_geometry_centroid(reconstructed_feature_geometry.get_reconstructed_geometry()))
distance_reconstructed_in_kms = distance_reconstructed * pygplates.Earth.mean_radius_in_kms

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/4420.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【分布式】ceph存储

目录 一、存储基础单机存储设备单机存储的问题商业存储解决方案 二、 分布式存储 (软件定义的存储 SDS)分布式存储的类型Ceph 优势Ceph 架构Ceph 核心组件Pool中数据保存方式支持两种类型:Pool、PG 和 OSD 的关系OSD 存储后端BlueStore 的主要…

stable diffusion webui mov2mov

手把手教你用stable diffusion绘画ai插件mov2mov生成动画_哔哩哔哩_bilibili手把手教你用stable diffusion绘画ai插件mov2mov生成动画, 视频播放量 14552、弹幕量 3、点赞数 275、投硬币枚数 114、收藏人数 980、转发人数 75, 视频作者 懂你的冷兮, 作者简介 科技改变世界&…

Ubuntu新版静态IP设置

cd /etc/netplan直接编辑 sudo vi /etc/netplan/00-installer-config.yaml#network: # ethernets: # ens160: # dhcp4: true # version: 2network:version: 2ethernets:ens160:dhcp4: noaddresses: [172.17.10.23/24]optional: truegateway4: 172.17.10.1nameservers…

多个excel文件合并为一个excel

Python openpyxl库实现将同一目录下的excel文件合并为一个excel功能(包含格式,不含宏),运行程序后,输入要生成的excel文件名称即可 import os import copy import sys import openpyxl from openpyxl.utils import get_column_letter import…

java工作随笔

String s JSONObject.toJSONString(fixedAsset);logger.error("-----------------8------------------" s);CusFixedAssettDTO CusFixedAssettDTO3 JSONObject.parseObject(s, CusFixedAssettDTO.class);父类转子类 相同对象合并 import java.util.Date; Data T…

Bard:Google AI开始支持中文对话和看图说话了

说起时下火爆的生成式AI,并不是只有ChatGPT。Bard也是一个很优秀的产品,并且刚刚发布的很多有趣的新功能。文末告诉你如何访问Bard。 Google AI在最近的更新中发布了Bard,一个新的语言模型。Bard支持多种语言,包括中文&#xff0…

【Ceph集群应用】Ceph对象存储系统之RGW接口详解

Ceph对象存储系统之RGW接口详解 1.创建Ceph对象存储系统RGW接口2. 开启httphttps,更改监听端口3. 更改监听端口4.S3接口访问测试5.实验中遇到的故障案例 接上文基于ceph-deploy部署Ceph集群详解 1.创建Ceph对象存储系统RGW接口 (1)对象存储概念 对象存…

WPF嵌入外部exe应用程序-使用Winfom控件承载外部程序

使用Winform控件承载外部程序 在WPF中使用Winfom控件添加winform相关的程序集在XAML头中加入对这两个程序集命名空间的引用使用Winform控件效果:问题 在Winfom控件中嵌入exe程序准备Winfrom控件更换父窗体的句柄完整实现代码:实现效果: 问题和…

ffmpeg get_bits() 函数是如何工作的.

/* * 本测试程序用来说明 GetBitContext 是什么? * get_bits() 函数是如何工作的. * author: hjjdebug * date: 2023年 07月 13日 星期四 16:42:45 CST */ GetBitContext 是如下定义的, 没有什么特别之处,指明了数据指针buffer,size_in_bits. typedef struct GetBi…

perl输出中文乱码【win10】

perl输出中文乱码 运行的时候输出的内容变成了中文乱码,原因首先来查找一下自己的perl的模块里面是否有Encode-CN。请运行打开你的cmd并输入perldoc -l Encode::CN 如果出现了地址 则就是有,如果没有需要进行该模块的安装。 安装方式有很多种&#xff0…

MetaTown:一个可以自己构建数字资产的平台

摘要:华为云Solution as Code重磅推出《基于MetaTown构建数字资产平台》解决方案。 本文分享自华为云社区《基于MetaTown构建数字资产平台》,作者: 阿米托福。 华为云Solution as Code重磅推出《基于MetaTown构建数字资产平台》解决方案&…

“掌握更多的快速排序技巧:三路划分、双路快排和非递归的深入理解”

快速排序是一种基于分治思想的排序算法,它能够以极快的速度将一个乱序的数组重新排列成有序的序列。不仅如此,快速排序还具有简洁的实现代码和良好的可扩展性,成为最受欢迎的排序算法之一。接下来,让我带你了解一下它的魅力吧&…

Linux系统部署Nginx详细教程(图文讲解)

前言:本篇博客记录了我是如何使用Linux系统一步一步部署Nginx的完整过程,也是我学习之路上的一个笔记总结,每一行代码都进行了严格的测试,特此做一个技术分享! 目录 一、安装依赖 二、安装Nginx 三、配置Nginx 四、…

visio 图片转换到 latex 中

调整图片大小 在Visio中,设计–>页面设置–>大小–>适应绘图,这样会自动去除多余空白,保留部分空白作为边界,无需使用Word。 2. 将新的Visio文件另存为pdf格式文件 3. latex 中插入pdf 格式图片

手把手教你搭建SpringCloud项目(八)集成Ribbon负载均衡器

什么是微服务?一看就会系列! 一、手把手教你搭建SpringCloud项目(一)图文详解,傻瓜式操作 二、手把手教你搭建SpringCloud项目(二)生产者与消费者 三、手把手教你搭建SpringCloud项目&#x…

【数据结构】24王道考研笔记——图

六、图 目录 六、图定义及基本术语图的定义有向图以及无向图简单图以及多重图度顶点-顶点间关系连通图、强连通图子图连通分量强连通分量生成树生成森林边的权、带权网/图特殊形态的图 图的存储及基本操作邻接矩阵邻接表法十字链表邻接多重表分析对比图的基本操作 图的遍历广度…

区块链学习一:安装 Geth

区块链技术正在迅速发展,并在各个领域展现出巨大的潜力。如果您对区块链开发感兴趣,那么学习如何安装和配置以太坊客户端 Geth 是一个很好的起点。在本篇文章中,我们将介绍如何使用 DockerHub 上的 Geth Docker 映像来快速安装和运行 Geth 节…

vue学习笔记(一)

1.编辑器选择 是用vscode 和 webstrom 个人感觉 vscode的插件比较多,对vue3的支持比较好 webstorm的自动保存比较好 各有优劣吧 我学习的这个项目目前采用vscode 2.vue2 还是 vue3 框架学通了都是通用的,这个时间点来学肯定是学vue3 只是顾虑到团…

JavaScript XHR、Fetch

1 前端数据请求方式 2 Http协议的解析 3 XHR的基本用法 4 XHR的进阶和封装 5 Fetch的使用详解 6 前端文件上传流程 早期的页面都是后端做好,浏览器直接拿到页面展示的,用到的是jsp、asp、php等等的语言。 这个叫做服务器端渲染SSR。 这里后端向前端…

金融数据库的战场,太平洋保险和OceanBase打了场胜仗

点击关注 文丨刘雨琦 “数据库的国产替代,必须经过严格的考虑,保证不会出错,所以大多数企业的领导层选择按兵不动或者简单扩容。因为不换就不会错,选了很久如果选错,还可能会出现重大事故。” 某银行数据库技术人员…