2021-11MatchLocate模板匹配

最新推荐文章于 2024-06-06 23:54:12 发布
最新推荐文章于 2024-06-06 23:54:12 发布
阅读量898 收藏 2
点赞数
文章标签： perl 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44259058/article/details/121150081
版权

目录

前言

有关于微震检测定位的算法——匹配定位算法Match&Locate的学习，github主页包括有MatchLocate2张淼老师的的程序包，GPU-MatchLocate1.0 :GPU-based Match and Locate是刘敏老师的程序包，笔者正好需要用到相应程序做测试，使用的是Ubuntu系统，以此记录环境配置过程进行学习。

1、MatchLocate2.0环境配置

一、软件编译

cd src
make

二、测试例子example

1. data preparation and processing
#Put your templates and continuous seismic data in "Template" and "Trace", respectively.
#Seismic data in this example could be downloaded from Hi-net (http://www.hinet.bosai.go.jp/).
#Time window is from 2012/09/02 03:00:00.0 to 2012/09/02 03:59.59.99 in "Trace/20120902". [prepared already]
cd Template #make sure you have templates in SAC format with correct headers
unzip waveforms-of-templates.zip #unzip all templates
perl marktaup.pl #mark t1 and t2 for P and S waves, calculate slowness
#You need TauP and a local velocity model (if not, use prem or iasp91)

cd ../Trace #continuous data
#make sure both templates and continuous have the same station, sampling rate, and frequency filtering.  
perl SACH_O.pl 20120902 #All time will be relative to ZERO of the day.


2. run Match&Locate
perl RunprocAll.pl catalog.dat 
#change corresponding parameters in RunproAll.pl and PROC_MatchLocate.pl

3. select events (in case one event is detected by multiple templates)
cd MultipleTemplate
perl SelectFinal.pl 2012 09 02 #merge all events and select final events, change threshold as needed
#DetectedFinal.dat would be your final catalog
perl PlotEventWaveform.pl DetectedFinal.dat #plot waveform comparison

三、测试例子example_CCdistribution

#template event: 2012/09/02 03:24:13.12 37.788 140.001 8.2 M 3.0
#slave event: 2012/09/02 03:33:51.61 37.791 140.005 7.8 M 2.3

1. copy ../example/Template/20120902032413.12 to Template directory [skip it, prepared]
2. copy ../example/Template/20120902033351.61 to Trace directory [skip it, prepared]
3. resample both sampling interval from 0.01 to 0.001 [skip it, prepared]
（cd Trace; unzip 20120902.zip）
4. run MatchLocate to get CCdistribution.out (perl RunproAll.pl catalog.dat)
5. go to plot directory and plot the CC distribution (perl plot_ccdistribution.pl)

问题

1、PlotEventWaveform.pl代码看不懂

2、GPU-MatchLocate1.0环境配置

一、Requirements
ObsPy (https://github.com/obspy/obspy, for data downloading, processing and picking only)
SAC (http://ds.iris.edu/ds/nodes/dmc/software/downloads/sac, for data processing only)
TauP(http://www.seis.sc.edu/taup/index.html, for P/S wave arrival time marking only)
GMT (https://www.soest.hawaii.edu/gmt, for figure plotting only)
CUDA Toolkit(https://developer.nvidia.com/cuda-downloads, for compiling cuda code)
Pssac(http://gmt.soest.hawaii.edu/doc/latest/supplements/meca/pssac.html, for figure plotting only)

说明：之前已经安装SAC、GMT6.2（pssac的相关功能也在其中，古暂时不考虑安装pssac）
因此需要安装的软件包括：obspy、taup、cuda tookit

二、软件安装
🕐obspy安装

git clone https://github.com/obspy/obspy.git

若网络不行直接下载安装包👇
解压unzip obspy-master.zip
进入目录cd obspy-master
安装sudo python setup.py develop
（出错，提示先安装numpy）故sudo apt-get install python3-numpy
注意python-numpy和python3-numpy的区别，若是前者则会有报错ImportError: No module named numpy. Please install numpy first, it is needed before installing ObsPy，
（还有其他相关库，参考）
sudo apt-get install python3-matplotlib
sudo apt-get install python3-future
sudo apt-get install python3-scipy
sudo apt-get install python3-pathlib
检测安装numpy是否成功进入python，然后import numpy，没有报错即成功。
然而再次测试sudo apt-get install python3-numpy有新报错error: invalid command ‘develop’，这个问题没有找到好的解决方法，在github看都同样的问题讨论但还是找不到合适的解决方法，放弃源码安装。直接上pip源，sudo apt install python3-pip，pip3 install obspy
前面还是顺利的，然而错误又来了！
Could not find a version that satisfies the requirement setuptools (from obspy) (from versions: )No matching distribution found for setuptools (from obspy)
解决先升级pippip3 install --upgrade pip，然后再试一次成功了！
pip3 install obspy -i http://pypi.douban.com/simple/ --trusted-host pypi.douban.com
检测是否成功，进入python，然后import obspy，没有报错即成功
总结安装obspy库代码如下（还是建议在conda下安装吧，不然问题好多…）

sudo apt-get install python3-matplotlib
sudo apt-get install python3-future
sudo apt-get install python3-scipy
sudo apt-get install python3-pathlib
pip3 install --upgrade pip
pip3 install obspy

🕑taup安装

wget http://www.seis.sc.edu/downloads/TauP/TauP-2.5.0.tgz
tar zxf TauP-2.5.0.tgz
sudo mv TauP-2.5.0 /usr/local/share
添加环境变量如下：(gedit ~/.bashrc)
export TAUP_HOME=/usr/local/share/TauP-2.5.0
export PATH=${TAUP_HOME}/bin:$PATH

🕒cuda tookit
一、检查机器上显卡信息
lspci|grep -i nvidia

03:00.0 VGA compatible controller: NVIDIA Corporation GP107 [GeForce GTX 1050 Ti] (rev a1)
03:00.1 Audio device: NVIDIA Corporation GP107GL High Definition Audio Controller (rev a1)

uname -m && cat /etc/*release有信息输出说明机器支持cuda

三、测试例子

5.1 compile source code
	$cd ./src
	$make
	$cd sacCC
	$make
5.2 download continuous data (User-specified)
	$cd ../../Demo
	$python data_download.py
	
5.3 download templates and create routine catalog (User-specified)
	$python template_download.py
	
5.4 set common origin time
	$cd ./Trace
	$perl SACH_O.pl 20190704
	
5.5 mark P/S wave arrival time and calculate slowness
	$cd ../Template
	$perl marktaup_p.pl
	$perl marktaup_s.pl

5.6 calculate weighting factor for each trace
	$bash cal_weights.sh

5.7 detect and locate events
	$cd ..
	$perl RunprocAll.pl

5.8 generate new catalog
	$cd ./MultipleTemplate
	$perl MergeEvents.pl 20190704
	$perl SelectFinal.pl 2019 07 04 Allevents

5.9 compare waveforms between template and detection
	$perl PlotEventWaveform.pl DetectedFinal.dat 1

问题

1、src目录下的编译文件Makefile运行make编译报错：
/usr/bin/ld:找不到 -lsac
/usr/bin/ld:找不到 -lsacio
解决主要是关于sac的库不在/opt/sac/lib路径，而是/usr/local/sac/lib,故此需要修改相应路径！
2、出错
🕐GPU_MatchLocate编译失败

nvcc -arch=sm_61  -Xcompiler "-fopenmp -fPIC" -Xlinker -lgomp -o ../bin/GPU_MatchLocate GPU_MatchLocate.cu sacio.o -L/usr/local/sac/lib -lsac -lsacio -lm

🕑执行bash cal_weights.sh报错sh: 1: …/…/…/bin/sac_e: not found，原因是sac_e编译失败

cc -L/usr/local/sac/lib -lsac -lsacio -lm    -o ../../bin/sac_e sac_e.o -L/usr/local/sac/lib -lsac -lsacio -lm	
sac_e.o: file not recognized: file format not recognized
collect2: error: ld returned 1 exit status
make: *** [makefile:7：sac_e] 错误 1

👉perl marktaup_p.pl
👉perl marktaup_s.pl
👉bash cal_weights.sh

(base) c@c-M:~/GPU-MatchLocate1.0-master/Demo/Template$ perl marktaup_p.pl 
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

CI.CWC.HHE CI.CWC.HHN CI.CWC.HHZ CI.EDW2.HNE CI.EDW2.HNN CI.EDW2.HNZ CI.FUR.HHE CI.FUR.HHN CI.FUR.HHZ CI.GSC.HNE CI.GSC.HNN CI.GSC.HNZ CI.ISA.HNE CI.ISA.HNN CI.ISA.HNZ CI.LRL.HHE CI.LRL.HHN CI.LRL.HHZ CI.SLA.HHE CI.SLA.HHN CI.SLA.HHZ
./20190704220235.06/CI.CWC.HHE
./20190704220235.06/CI.CWC.HHN
./20190704220235.06/CI.CWC.HHZ
./20190704220235.06/CI.EDW2.HNE
./20190704220235.06/CI.EDW2.HNN
./20190704220235.06/CI.EDW2.HNZ
./20190704220235.06/CI.FUR.HHE
./20190704220235.06/CI.FUR.HHN
./20190704220235.06/CI.FUR.HHZ
./20190704220235.06/CI.GSC.HNE
./20190704220235.06/CI.GSC.HNN
./20190704220235.06/CI.GSC.HNZ
./20190704220235.06/CI.ISA.HNE
./20190704220235.06/CI.ISA.HNN
./20190704220235.06/CI.ISA.HNZ
./20190704220235.06/CI.LRL.HHE
./20190704220235.06/CI.LRL.HHN
./20190704220235.06/CI.LRL.HHZ
./20190704220235.06/CI.SLA.HHE
./20190704220235.06/CI.SLA.HHN
./20190704220235.06/CI.SLA.HHZ
./20190704220503.21/CI.CWC.HHE
./20190704220503.21/CI.CWC.HHN
./20190704220503.21/CI.EDW2.HNE
./20190704220503.21/CI.EDW2.HNN
./20190704220503.21/CI.EDW2.HNZ
./20190704220503.21/CI.FUR.HHE
./20190704220503.21/CI.FUR.HHN
./20190704220503.21/CI.FUR.HHZ
./20190704220503.21/CI.GSC.HNE
./20190704220503.21/CI.GSC.HNN
./20190704220503.21/CI.GSC.HNZ
./20190704220503.21/CI.ISA.HNE
./20190704220503.21/CI.ISA.HNN
./20190704220503.21/CI.ISA.HNZ
./20190704220503.21/CI.LRL.HHE
./20190704220503.21/CI.LRL.HHN
./20190704220503.21/CI.LRL.HHZ
./20190704220503.21/CI.SLA.HHE
./20190704220503.21/CI.SLA.HHN
./20190704220503.21/CI.SLA.HHZ
./20190704220522.27/CI.CWC.HHE
./20190704220522.27/CI.CWC.HHN
./20190704220522.27/CI.CWC.HHZ
./20190704220522.27/CI.EDW2.HNE
./20190704220522.27/CI.EDW2.HNN
./20190704220522.27/CI.EDW2.HNZ
./20190704220522.27/CI.FUR.HHE
./20190704220522.27/CI.FUR.HHN
./20190704220522.27/CI.FUR.HHZ
./20190704220522.27/CI.GSC.HNE
./20190704220522.27/CI.GSC.HNN
./20190704220522.27/CI.GSC.HNZ
./20190704220522.27/CI.ISA.HNE
./20190704220522.27/CI.ISA.HNN
./20190704220522.27/CI.ISA.HNZ
./20190704220522.27/CI.LRL.HHE
./20190704220522.27/CI.LRL.HHN
./20190704220522.27/CI.LRL.HHZ
./20190704220522.27/CI.SLA.HHE
./20190704220522.27/CI.SLA.HHN
./20190704220522.27/CI.SLA.HHZ
./20190704221208.49/CI.CWC.HHE
./20190704221208.49/CI.CWC.HHN
./20190704221208.49/CI.CWC.HHZ
./20190704221208.49/CI.EDW2.HNE
./20190704221208.49/CI.EDW2.HNN
./20190704221208.49/CI.EDW2.HNZ
./20190704221208.49/CI.FUR.HHE
./20190704221208.49/CI.FUR.HHN
./20190704221208.49/CI.FUR.HHZ
./20190704221208.49/CI.GSC.HNE
./20190704221208.49/CI.GSC.HNN
./20190704221208.49/CI.GSC.HNZ
./20190704221208.49/CI.ISA.HNE
./20190704221208.49/CI.ISA.HNN
./20190704221208.49/CI.ISA.HNZ
./20190704221208.49/CI.LRL.HHE
./20190704221208.49/CI.LRL.HHN
./20190704221208.49/CI.LRL.HHZ
./20190704221208.49/CI.SLA.HHE
./20190704221208.49/CI.SLA.HHN
./20190704221208.49/CI.SLA.HHZ
./20190704222906.57/CI.CWC.HHE
./20190704222906.57/CI.CWC.HHN
./20190704222906.57/CI.CWC.HHZ
./20190704222906.57/CI.EDW2.HNE
./20190704222906.57/CI.EDW2.HNN
./20190704222906.57/CI.EDW2.HNZ
./20190704222906.57/CI.FUR.HHE
./20190704222906.57/CI.FUR.HHN
./20190704222906.57/CI.FUR.HHZ
./20190704222906.57/CI.GSC.HNE
./20190704222906.57/CI.GSC.HNN
./20190704222906.57/CI.GSC.HNZ
./20190704222906.57/CI.ISA.HNE
./20190704222906.57/CI.ISA.HNN
./20190704222906.57/CI.ISA.HNZ
./20190704222906.57/CI.LRL.HHE
./20190704222906.57/CI.LRL.HHN
./20190704222906.57/CI.LRL.HHZ
./20190704222906.57/CI.SLA.HHE
./20190704222906.57/CI.SLA.HHN
./20190704222906.57/CI.SLA.HHZ
./20190704223007.80/CI.CWC.HHE
./20190704223007.80/CI.CWC.HHN
./20190704223007.80/CI.CWC.HHZ
./20190704223007.80/CI.EDW2.HNE
./20190704223007.80/CI.EDW2.HNN
./20190704223007.80/CI.EDW2.HNZ
./20190704223007.80/CI.FUR.HHE
./20190704223007.80/CI.FUR.HHN
./20190704223007.80/CI.FUR.HHZ
./20190704223007.80/CI.GSC.HNE
./20190704223007.80/CI.GSC.HNN
./20190704223007.80/CI.GSC.HNZ
./20190704223007.80/CI.ISA.HNE
./20190704223007.80/CI.ISA.HNN
./20190704223007.80/CI.ISA.HNZ
./20190704223007.80/CI.LRL.HHE
./20190704223007.80/CI.LRL.HHN
./20190704223007.80/CI.LRL.HHZ
./20190704223007.80/CI.SLA.HHE
./20190704223007.80/CI.SLA.HHN
./20190704223007.80/CI.SLA.HHZ
./20190704224221.78/CI.CWC.HHE
./20190704224221.78/CI.CWC.HHN
./20190704224221.78/CI.CWC.HHZ
./20190704224221.78/CI.EDW2.HNE
./20190704224221.78/CI.EDW2.HNN
./20190704224221.78/CI.EDW2.HNZ
./20190704224221.78/CI.FUR.HHE
./20190704224221.78/CI.FUR.HHN
./20190704224221.78/CI.FUR.HHZ
./20190704224221.78/CI.GSC.HNE
./20190704224221.78/CI.GSC.HNN
./20190704224221.78/CI.GSC.HNZ
./20190704224221.78/CI.ISA.HNE
./20190704224221.78/CI.ISA.HNN
./20190704224221.78/CI.ISA.HNZ
./20190704224221.78/CI.LRL.HHE
./20190704224221.78/CI.LRL.HHN
./20190704224221.78/CI.LRL.HHZ
./20190704224221.78/CI.SLA.HHE
./20190704224221.78/CI.SLA.HHN
./20190704224221.78/CI.SLA.HHZ
(base) c@c-M:~/GPU-MatchLocate1.0-master/Demo/Template$ perl marktaup_s.pl 
 SEISMIC ANALYSIS CODE [11/11/2013 (Version 101.6a)]
 Copyright 1995 Regents of the University of California

./20190704220235.06/CI.CWC.HHE
./20190704220235.06/CI.CWC.HHN
./20190704220235.06/CI.CWC.HHZ
./20190704220235.06/CI.EDW2.HNE
./20190704220235.06/CI.EDW2.HNN
./20190704220235.06/CI.EDW2.HNZ
./20190704220235.06/CI.FUR.HHE
./20190704220235.06/CI.FUR.HHN
./20190704220235.06/CI.FUR.HHZ
./20190704220235.06/CI.GSC.HNE
./20190704220235.06/CI.GSC.HNN
./20190704220235.06/CI.GSC.HNZ
./20190704220235.06/CI.ISA.HNE
./20190704220235.06/CI.ISA.HNN
./20190704220235.06/CI.ISA.HNZ
./20190704220235.06/CI.LRL.HHE
./20190704220235.06/CI.LRL.HHN
./20190704220235.06/CI.LRL.HHZ
./20190704220235.06/CI.SLA.HHE
./20190704220235.06/CI.SLA.HHN
./20190704220235.06/CI.SLA.HHZ
./20190704220503.21/CI.CWC.HHE
./20190704220503.21/CI.CWC.HHN
./20190704220503.21/CI.EDW2.HNE
./20190704220503.21/CI.EDW2.HNN
./20190704220503.21/CI.EDW2.HNZ
./20190704220503.21/CI.FUR.HHE
./20190704220503.21/CI.FUR.HHN
./20190704220503.21/CI.FUR.HHZ
./20190704220503.21/CI.GSC.HNE
./20190704220503.21/CI.GSC.HNN
./20190704220503.21/CI.GSC.HNZ
./20190704220503.21/CI.ISA.HNE
./20190704220503.21/CI.ISA.HNN
./20190704220503.21/CI.ISA.HNZ
./20190704220503.21/CI.LRL.HHE
./20190704220503.21/CI.LRL.HHN
./20190704220503.21/CI.LRL.HHZ
./20190704220503.21/CI.SLA.HHE
./20190704220503.21/CI.SLA.HHN
./20190704220503.21/CI.SLA.HHZ
./20190704220522.27/CI.CWC.HHE
./20190704220522.27/CI.CWC.HHN
./20190704220522.27/CI.CWC.HHZ
./20190704220522.27/CI.EDW2.HNE
./20190704220522.27/CI.EDW2.HNN
./20190704220522.27/CI.EDW2.HNZ
./20190704220522.27/CI.FUR.HHE
./20190704220522.27/CI.FUR.HHN
./20190704220522.27/CI.FUR.HHZ
./20190704220522.27/CI.GSC.HNE
./20190704220522.27/CI.GSC.HNN
./20190704220522.27/CI.GSC.HNZ
./20190704220522.27/CI.ISA.HNE
./20190704220522.27/CI.ISA.HNN
./20190704220522.27/CI.ISA.HNZ
./20190704220522.27/CI.LRL.HHE
./20190704220522.27/CI.LRL.HHN
./20190704220522.27/CI.LRL.HHZ
./20190704220522.27/CI.SLA.HHE
./20190704220522.27/CI.SLA.HHN
./20190704220522.27/CI.SLA.HHZ
./20190704221208.49/CI.CWC.HHE
./20190704221208.49/CI.CWC.HHN
./20190704221208.49/CI.CWC.HHZ
./20190704221208.49/CI.EDW2.HNE
./20190704221208.49/CI.EDW2.HNN
./20190704221208.49/CI.EDW2.HNZ
./20190704221208.49/CI.FUR.HHE
./20190704221208.49/CI.FUR.HHN
./20190704221208.49/CI.FUR.HHZ
./20190704221208.49/CI.GSC.HNE
./20190704221208.49/CI.GSC.HNN
./20190704221208.49/CI.GSC.HNZ
./20190704221208.49/CI.ISA.HNE
./20190704221208.49/CI.ISA.HNN
./20190704221208.49/CI.ISA.HNZ
./20190704221208.49/CI.LRL.HHE
./20190704221208.49/CI.LRL.HHN
./20190704221208.49/CI.LRL.HHZ
./20190704221208.49/CI.SLA.HHE
./20190704221208.49/CI.SLA.HHN
./20190704221208.49/CI.SLA.HHZ
./20190704222906.57/CI.CWC.HHE
./20190704222906.57/CI.CWC.HHN
./20190704222906.57/CI.CWC.HHZ
./20190704222906.57/CI.EDW2.HNE
./20190704222906.57/CI.EDW2.HNN
./20190704222906.57/CI.EDW2.HNZ
./20190704222906.57/CI.FUR.HHE
./20190704222906.57/CI.FUR.HHN
./20190704222906.57/CI.FUR.HHZ
./20190704222906.57/CI.GSC.HNE
./20190704222906.57/CI.GSC.HNN
./20190704222906.57/CI.GSC.HNZ
./20190704222906.57/CI.ISA.HNE
./20190704222906.57/CI.ISA.HNN
./20190704222906.57/CI.ISA.HNZ
./20190704222906.57/CI.LRL.HHE
./20190704222906.57/CI.LRL.HHN
./20190704222906.57/CI.LRL.HHZ
./20190704222906.57/CI.SLA.HHE
./20190704222906.57/CI.SLA.HHN
./20190704222906.57/CI.SLA.HHZ
./20190704223007.80/CI.CWC.HHE
./20190704223007.80/CI.CWC.HHN
./20190704223007.80/CI.CWC.HHZ
./20190704223007.80/CI.EDW2.HNE
./20190704223007.80/CI.EDW2.HNN
./20190704223007.80/CI.EDW2.HNZ
./20190704223007.80/CI.FUR.HHE
./20190704223007.80/CI.FUR.HHN
./20190704223007.80/CI.FUR.HHZ
./20190704223007.80/CI.GSC.HNE
./20190704223007.80/CI.GSC.HNN
./20190704223007.80/CI.GSC.HNZ
./20190704223007.80/CI.ISA.HNE
./20190704223007.80/CI.ISA.HNN
./20190704223007.80/CI.ISA.HNZ
./20190704223007.80/CI.LRL.HHE
./20190704223007.80/CI.LRL.HHN
./20190704223007.80/CI.LRL.HHZ
./20190704223007.80/CI.SLA.HHE
./20190704223007.80/CI.SLA.HHN
./20190704223007.80/CI.SLA.HHZ
./20190704224221.78/CI.CWC.HHE
./20190704224221.78/CI.CWC.HHN
./20190704224221.78/CI.CWC.HHZ
./20190704224221.78/CI.EDW2.HNE
./20190704224221.78/CI.EDW2.HNN
./20190704224221.78/CI.EDW2.HNZ
./20190704224221.78/CI.FUR.HHE
./20190704224221.78/CI.FUR.HHN
./20190704224221.78/CI.FUR.HHZ
./20190704224221.78/CI.GSC.HNE
./20190704224221.78/CI.GSC.HNN
./20190704224221.78/CI.GSC.HNZ
./20190704224221.78/CI.ISA.HNE
./20190704224221.78/CI.ISA.HNN
./20190704224221.78/CI.ISA.HNZ
./20190704224221.78/CI.LRL.HHE
./20190704224221.78/CI.LRL.HHN
./20190704224221.78/CI.LRL.HHZ
./20190704224221.78/CI.SLA.HHE
./20190704224221.78/CI.SLA.HHN
./20190704224221.78/CI.SLA.HHZ

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
👉perl RunprocAll.pl

(base) c@c-M:~/GPU-MatchLocate1.0/Demo$ ./RunprocAll.pl 
20190704
******Start Running GPU_Match&Locate on 2022年 01月 15日 星期六 18:17:17 CST
There are 1 potential locations!
The shape of data is wrong, Please check the number of stations
******End Running GPU_Match&Locate on 2022年 01月 15日 星期六 18:17:17 CST
mv: 无法获取'EventCase.out' 的文件状态(stat): 没有那个文件或目录

找到报错The shape of data is wrong, Please check the number of stations所在源码GPU_MatchLocate.cu，对代码中变量ntrace添加printf显示，重新编译，运行perl RunprocAll.pl发现提示ntrace=20，少了第一个台站HHZ的数据。(发现INPUT.in文件是在执行perl RunprocAll.pl时生成的。而20190704220503.21恰好却了HHE的数据，暂时先编一个跟HHN一样。然而还是出错，于是直接删掉台站CWI在INPUT七个事件中的记录，再次运行)
在这里插入图片描述
在这里插入图片描述
👉perl MergeEvents.pl 20190704
生成文件Allevents
👉perl SelectFinal.pl 2019 07 04 Allevents
生成文件DetectFinnal.dat
👉perl PlotEventWaveform.pl DetectedFinal.dat 1
(现在剩下gmt4画图的问题了)

(base) c@c-M:~/GPU-MatchLocate1.0/Demo/MultipleTemplate$ perl PlotEventWaveform.pl DetectedFinal.dat 1
Unquoted string "jk" may clash with future reserved word at PlotEventWaveform.pl line 93.
Name "main::tt" used only once: possible typo at PlotEventWaveform.pl line 138.
Name "main::coefav" used only once: possible typo at PlotEventWaveform.pl line 145.
20190704161343.428
sh: 1: pssac: not found
sh: 1: pssac: not found
sh: 1: psxy: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: pstext: not found
sh: 1: psxy: not found
yep吖
关注
CUDA图像处理 | 模板匹配
【黑键】
09-10 2943
模板匹配 OpenCV中的模板匹配 该部分转载自 【OpenCV3】模板匹配——cv::matchTemplate()详解 匹配方法 cv::TM_SQDIFF：该方法使用平方差进行匹配，因此最佳的匹配结果在结果为0处，值越大匹配结果越差。 cv::TM_CCORR：相关性匹配方法，该方法使用源图像与模板图像的卷积结果进行匹配，因此，最佳匹配位置在值最大处，值越小匹配结...
CUDA-GPT:CUDA和GPT图像模板匹配的算法加速
05-10
CUDA中的加速GPT 用法 编译并执行 nvcc *.cpp include/*.c* -o cu.out -arch=sm_60 --use_fast_math 从表的#define MAKETEMP 0到#define MAKETEMP 1修改include / parameter.h第142行。 编译并执行。 然后改回#define MAKETEMP 0 。 如果以上步骤未完成，则Cannot open the file! 错误将会出现。 在Tokyo Tech TSUBAME 3.0上测试 英伟达Tesla P100 Cuda编译工具，版本8.0，V8.0.61 已知的问题 CPU和GPU之间的微小结果差异 stdGpt.cpp中的bilinear_normal_projection函数中使用的multipliVect3x3中可能是CPU精度错误 参考 @inpr
图像处理(8) : 模板匹配
QtHalcon
06-09 3万+
​模板匹配指的是通过模板图像与测试图像之间的比较，找到测试图像上与模板图像相似的部分，这是通过计算模板图像与测试图像中目标的相似度来实现的，可以快速地在测试图像中定位出预定义的目标。匹配的主要思路是使用一个目标原型，根据它创建一个模板，在测试图像中搜索与该模板图像最相似的目标，并寻找与该模板的均值或方差最接近的区域。 通过模板匹配可以得到目标的相似度，旋转角度，行列坐标，缩放大小等。 针对不同的图像特征和检测环境，有多种模板匹配算法。如何选择合适的模板匹配算法，取决于具体的图像数...
opencv(C++)GPU、CPU 模板匹配
qq_40622955的博客
07-15 5373
摘要： 本文主要关注opencv常规版和cuda版的模板匹配算法，网上cuda版的资料不多，这里做个记录，以后用到也好有个参考。 文章目录摘要：1. opencv cuda版配置2. 源码3. 结果4.总结
图像处理常用算法GPU实现三：基于模板匹配的边缘检测
u010607947的博客
08-14 3080
上一篇文章给出了基于微分方法的边缘检测，该篇给出一个基于Prewitt算子的模板匹配检测算法。除了水平和垂直方向外，此处还给出了其他方向的检测，一共8个方向，算子如下： 代码如下： [cpp] view plain copy /********************************  *  Author: rabbit729 
图像拼接（二）：柱面投影+模板匹配+渐入渐出融合（GPU版）
czl389的专栏
01-18 8066
本篇博客的代码修改自图像拼接（一）：柱面投影+模板匹配+渐入渐出融合。新的代码充分利用了OpenCV库函数的GPU版本。在一些方法中使用了GpuMat这一数据结构，它是Mat的替代。Mat运行在CPU上，GpuMat运行在GPU上，使用这一数据结构，将一些费时的操作交由GPU处理，会提高处理速度。 文末附有完整代码清单，这里对代码做些说明： 1、使用OpenCV的GPU版本，需要重新编译Open
GPU驱动、CUDA和cuDNN之间的版本匹配与下载
eastking0530的博客
08-28 1万+
GPU驱动、CUDA和cuDNN之间的版本匹配与下载
模板匹配——基于模板与基于特征的两种模板匹配对比
知来者逆的博客
04-04 3628
在图像匹配算法可分为两大类别，基于模板的匹配和基于特征的匹配各具特色。 1.基于模板的方法，也被称作区域匹配法，适用于模板图像在图像中缺乏显著特征的情况。这种方法直接处理像素值，通过计算图像与模板的强度值来衡量匹配程度。 2.基于特征的方法则适用于源图像和模板图像在特征和控制点上有较高一致性的情况。此时，特征可以是点、曲线或表面模型，用于执行模板匹配。
模板匹配(Match Template) 修改版
Qingan Yan
06-16 1万+
作者：王先荣 http://www.cnblogs.com/xrwang/archive/2010/02/05/MatchTemplate.html修改者：MadTurtle------对于王先生的版本进行了修改，因为有些地方不对。前言    模板匹配是在图像中寻找目标的方法之一。Come On, Boy.我们一起来看看模板匹配到底是怎么回事。模板匹配的工作方式    模板匹配的工作方式跟直方图的反向投影基本一样，大致过程是这样的：通过在输入图像上滑动图像块对实际的图像块和输入图像进行匹配。    假设我们
基于边缘的模板匹配
10-12
Edge Based Template Matching 基于边界的模板匹配的原理及算法实现 vc++
opencv中常用cuda函数总结
Ring__Rain的博客
06-01 5638
bitwise_and() bitwise_not() bitwise_or() bitwise_xor() compare() divide()：除 exp() log() max() min() multiply() threshold()：二值化，但要指定设定阈值 blendLinear()：两幅图片的线形混合 calcHist() createBoxFilter ()：创建一个规范化的2D...
C++开发基础之初探CUDA计算环境搭建
最新发布
dotnet研习社
06-06 1670
项目中有使用到CUDA计算的相关内容。但是在早期CUDA计算环境搭建的过程中，并不是非常顺利，编写此篇文章记录下。对于刚刚开始研究的你可能会有一定的帮助。在这个项目中，我们主要体会框架的用法，以及CUDA计算环境搭建的。通过编码实践，构建项目成功实验了CUDA计算环境搭建，为接下来的工作准备好环境。
图像处理中的模板匹配c++实现
wu_lian_nan的博客
08-16 1万+
一、理论基础 基于相关的模板匹配技术可直接用于在一幅图像中，寻找某种子图像模式。对于大小为MxN的图像f(x,y)和大小为JxK的子图像模式w(x,y),f与w的相关可表示为：其种，x=0,1,2,…N-K,y=0,1,2,…M-J。此处 的目的是寻找匹配而不是对f(x,y)进行滤波操作，因此w的原点被设置在子图像的左上角， 并且式(11-3) 给出的形式也完全适用于J和K 为偶数的情况。 计算
基于模板匹配的运动目标跟踪
启功社区
07-07 5205
一、实验目的  通过该实验，初步了解模板匹配算法，体会并学习如何用该算法实现视频跟踪。  二、算法概要  模板就是一幅已知的小图像。模板匹配就是在一幅大图像中搜寻目标。已知在该图中有要寻找的目标，且该目标与模板有相同的尺寸、方向和图像，通过一定的算法可以在图中找到目标，确定其坐标位置。 如图1 所示，设模板T 叠放在搜索图S 上平移，模板覆盖下的那块搜索图叫做子图， 这块子图的左上角像素点在
【图像配准】基于灰度的模板匹配算法（一）：MAD、SAD、SSD、MSD、NCC、SSDA、SATD算法
热门推荐
hujingshuang
08-18 10万+
本文主要介绍几种基于灰度的图像匹配算法：平均绝对差算法（MAD）、绝对误差和算法（SAD）、误差平方和算法（SSD）、平均误差平方和算法（MSD）、归一化积相关算法（NCC）、序贯相似性算法（SSDA）。
模板匹配，特征点匹配-全
jacke121的专栏
07-09 7438
https://blog.csdn.net/zhuisui_woxin/article/details/84400439 模板匹配与特征匹配 模板匹配：模板匹配是一种最原始、最基本的模式识别方法，研究某一特定对象物的图案位于图像的什么地方，进而识别对象物，这就是一个匹配问题。它是图像处理中最基本、最常用的匹配方法。模板匹配具有自身的局限性，主要表现在它只能进行平行移动，若原图像中的匹配目标...
基于模板匹配的跟踪算法代码示例
Gone_HuiLin的博客
11-19 1527
#include using namespace cv; using namespace std; // Global variables static Rect box; static bool drawing_box = false; static bool gotBB = false; // bounding box mouse callback static void mouseH
GPU（CUDA）学习日记（十一）------ 深入理解CUDA线程层次以及关于设置线程数的思考
MySniper11的专栏
12-07 1万+
GPU线程以网格（grid）的方式组织，而每个网格中又包含若干个线程块，在G80/GT200系列中，每一个线程块最多可包含512个线程，Fermi架构中每个线程块支持高达1536个线程。同一线程块中的众多线程拥有相同的指令地址，不仅能够并行执行，而且能够通过共享存储器（Shared memory）和栅栏（barrier）实现块内通信。这样，同一网格内的不同块之间存在不需要通信的粗粒度并行，而一个块
模板匹配&&频域相位相关
靈飛
06-28 3074
最近在做图像拼接项目，原打算用大众化的sift点匹配然后拼接。发现sift算法运算量巨大，不能满足拼接的实时要求，同时特征点法需要不同程度的干预（比如，需要随机抽样来删选优质的匹配点对），不满足自动匹准的要求。最后，只得尝试着用基于区域匹配的方法。           基于区域的方法其实就是模板匹配了，opencv已经为我们提供函数: void cv::matchTemplate( I
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值