❶ GIS中坐标系与偏移算法总结
一 大地坐标系
1.1 概念
大地坐标系是大地测量中以参考椭球面为基准面建立起来的坐标系。
大地坐标系根据其原点的位置不同，分为地心坐标系和参心坐标系。地心坐标系的原点与地球质心重合，参心坐标系的原点与某一地区或国家所采用的参考椭球中心重合，通常与地球质心不重合。
1.2 常用的参心坐标系与地心坐标系
北京54 参心坐标系（参心坐标系）
西安80 参心坐标系（参心坐标系）
cgcs2000 地心坐标系（地心坐标系）
wgs84 地心坐标系（地心坐标系）
我国先后建立的1954年北京坐标系、1980西安坐标系和新1954年北京坐标系，都是参心坐标系。这些坐标系为我国经济社会发展和国防建设作出了重要贡献。
但是，随着现代科技的发展，特别是全球卫星定位技术的发展和应用，世界上许多发达国家和中等发达国家都已在多年前就开始使用地心坐标系。
国务院批准自2008年7月1日启用我国的地心坐标系——2000国家大地坐标系（CGCS-2000），同时要求用8-10年的时间，完成现行国家大地坐标系向20000国家大地坐标系的过渡和转换。过渡期结束，将停止提供现行国家大地坐标系下的测绘成果。
参考：
2018年7月1日起全面使用2000国家大地坐标系，西安80和北京54坐标系正式退出历史舞台
wgs84是为GPS全球定位系统使用而建立的坐标系统；gps设备采集的数据均为wgs84坐标系。
1.3 不同坐标系之间的转换
arcgis 软件中计算完成，参考:[arcgis坐标转换与投影]( https://www.jianshu.com/p/5c437696be06 ）
二 坐标投影
2.1 投影后的坐标形式
原始经纬度:120.0397529296875,30.229220825195313
墨卡托投影后:13362764.171082955,3533048.2025558753
参考： arcgis js api：web墨卡托(3857)转经纬度坐标(4326)
2.2 投影的目的
方便工程测量、二维图展示便于理解。
2.3 根据不同需求使用不同的投影算法
例如： 墨卡托投影后的二维图导致了地球两极被拉宽，不适用于工程测量，但适合用作普通二维图的展示。而大比例尺的工程测量图对局部区域内精度要求高，而采用高斯克里格投影。
三 gcj02偏移算法（国家测绘局，被戏称为火星坐标）
注意 这是偏移算法，而不是单独的一个坐标系。习惯上人们将加了偏移算法的坐标称为gcj02坐标。经偏移算法处理的地图数据偏差一般为 300~500 米。
3.1 在我们国家发布的互联网地图按法律规定需要经过偏移算法加偏移。
例如高德地图、腾讯地图。
3.2 网络地图在gcj02基础上进行了二次加偏移，称为bd09坐标
3.3 天地图是否加偏移？
天地图采用cgcs2000坐标。
发布在互联网上的天地图并不一定都是加偏移，未加偏移的地图做了特殊处理，很多涉密地理信息在地图上找不到。
3.4 使用了加偏移的地图如何进行gis开发
对叠加到地图上的数据同样加偏移，实现与底图吻合。
参考： WGS84坐标与不同加密算法之间转换
更多参考：
你必须知道的地理坐标系和投影坐标系
arcgis坐标转换与投影变换
❷ 关于单片机两种查表指令偏移量的算法 偏移量怎么计算啊帮帮忙
找到每条指令的字节数
把要跳过的所有指令的字节数都加起来就是偏移量
❸ 基于非稳态相移法叠前深度偏移方法应用研究
王伟国1钱荣毅2
（1.广州海洋地质调查局广州510760；2.中国地质大学北京100083）
作者简介：王伟国，男（1982—），硕士，助理工程师，工作方向为地震偏移成像及定量地震解释。E-mail:[email protected]
摘要普通的相移偏移算法具有速度快、稳定性高、理论上无网格频散等优点，但在复杂地质构造条件下，由于速度横向的剧烈变化，偏移成像往往不能得到理想的效果。而非稳态相移算法是基于在非稳态滤波器理论，通过对普通的相移算子在横向上增加一个扰动量（称之为非稳态相移算子），从而在偏移成像时能很好地适应速度的横向变化。用该算法在普通的横向变化速度模型及Marmousi模型上进行叠前深度偏移试算，均取得了很好的应用效果。
关键词非稳态相移叠前深度偏移速度模型
1前言
相移算法Seddon[1]预测一个深度层波场的振幅和相位主要是基于在更浅层的已知波场和地下速度模型。总体来说是基于波动方程的差分解，在这方面我国的很多学者也做了很多相关的研究[2～5]。相移法波场外推有许多前提条件和一个主要的难题。好的方面，对于常速度相移算子理论上是精确的，无条件稳定的，没有网格频散，并且对所有的散射角都是精确的。主要的难题是横向速度变化怎样结合到相移方法中，这还不是直观的表现，因为空间坐标进行了傅氏变换。已知的波场通常是实际检波器接收到的地震记录或是正演模拟的地震记录，相移波场被用于计算地下地质结构的反射系数。当介质速度是常量时，相移波场的过程是一个相当稳定的过程，并且给出相移角度能达到90°的精确解（Gazdag,1978;Stoffa等，1990）。Stoffa（1990）的裂步傅里叶方法和Wu（1994）的相位屏方法都实现了横变速度的近似相移，而Gazdag和Sguazzero[6]通过另一种方法叫相移加插值（PSPI），对于参考速度集而言，通过计算一系列的常速度相移，并插值得到一个单一的横向变化的结果，从而使得相移能够被拓展到横向变化的速度，但它仅仅是在数学上的一种处理方法，并不具有物理意义。Black等（1984）给出了傅氏方法的一个解析表达式来适应横向速度的变化，但并没有证明表达式。
Margrave和Ferguson[7-10]证明了Black的方法是广义的PSPI，它用很多的参考速度来代替插值，并利用非稳态滤波器[11]推导了和PSPI相联系的方法，称为非稳态相移（non stationary phase shift,NSPS）。与PSPI相比，NSPS具有物理上的解释。
2非稳态相移算子及算法流程
非稳态滤波器理论由Margrave（1998）提出，他提出非稳态滤波器理论至少有两种不同的形式是可能存在的。称之为组合滤波器和褶积滤波器，两种滤波器在稳态极限下是等价的。当对每一个不同的速度算出一个参考波场时，我们可以得到方程（1），它其实是一个广义反傅氏积分，是一个非稳态，双域的组合滤波器例子，滤波器的非稳态性通过了事实证明，滤波器描述αv（x）（kx,x，ω），取决于波数和空间位置。
南海地质研究.2010
其中，
南海地质研究.2010
这里，φ（kx,0，ω）为初始波场，ΨNSPS（x，Δz，ω）为利用NSPS外推Δz后的波场，αv（x）（kx,x，ω）在NSPS算法中称为偏移算子。kx,kz分别是横向和纵向的波数。
对公式（1）作正傅氏变换得到：
南海地质研究.2010
对公式（3）作反傅氏变换得到外推波场为：
南海地质研究.2010
通过非稳态相移的理论，可以按照图1的算法流程来进行叠前深度偏移。
3模型试算
考虑到计算机硬件资源和算法精度验证的要求，建立了一个断层－背斜模型，断层主要是加强模型的横向速度变化，深部的背斜是验证算法成像的最大倾角和深部成像对速度的要求。由于是已知模型，本文选用对速度更为敏感的炮集记录来进行叠前深度偏移。模型及参数如图2所示。
所布置的模拟观测系统参数为：双边接收系统，总炮数为50炮，炮点间隔为30m，首炮位置位于速度剖面的最左端，100道接收，道间距为5m，最小偏移距5m，排列自左向右移动。炮集模拟用的是单程声波波动方程。炮集震源为模拟爆炸震源，主频为30Hz。图3所示为第25炮，即模型中间750m位置处。具体在偏移过程中，添加了零道来满足偏移处理中与速度剖面的维数相同。图4为NSPS偏移结果。
图4中圆圈处为模型的断层点，从单炮偏移的效果看，圆圈处同相轴的纵向分辨率还是很高，断点也比较清晰，对比模型，归位也很准确，反射波和绕射波都已经收敛，600m处的一个同相轴被拉平，模型此处是水平地层，800m处背斜的顶点也能比较清晰地看到；这是没有做过任何叠前处理的炮记录，能够达到这样的效果，至少可以说明NSPS算法对于该模型是良好适应的。
图1 NSPS算法流程图Fig.1 Flow of NSPS
通过图5的叠加剖面可以看出，断层点清晰可见，归位很准确，断层上下盘界面清晰，水平层位被很好地拉平，且深度都基本和模型位置的深度一致，基本没有重影，没有频散现象，背斜轮廓明显可见。但同时也可以看出，0～100m之间存在明显的直达波影响，主要是没有做叠前的一些常规数据处理造成的，这并不影响对于算法本身的验证；逆掩断层的断面及背斜的两个倾斜角度能量不强，没有很好地收敛，其实产生这样的结果主要是因为算法本身的假设条件造成的，单程波动方程偏移对于多次绕射波在理论上无法很好地成像。
图2 断层－背斜深度模型Fig.2 Depth model of fault-anticline
图3 添加零道后的单炮记录Fig.3 Single shot record after padding with zero
图4 单炮NSPS偏移结果Fig.4 Migration of single shot
图5 NSPS偏移叠加结果Fig.5 NSPS migration after stacking shots
图6 Marmousi模型Fig.6 Marmousi model
我们可以看一下 NSPS 算法在IFP（Institut Francais Petrole，法国石油研究院）Marmousi模型上的表现，图6为Marm ousi模型，该2D 模型包含240个炮集记录，一个震源波形和一个完整的速度和密度剖面，本文采用的是简化的速度模型，只含有纵波速度，不含有密度、横波等信息，炮集记录也是重新模拟生成，如图7所示。炮集设置为240炮，左端接收，接收道96道，道间距25m，最小偏移距200m，模拟炸药震源放炮，波函数为零相位雷克子波。图6为第120炮位置，即在5550m处，图7为第120炮的单炮模拟记录。为了应用NSPS算法，对模型数据进行了抽稀，抽成25m×25m的网格。
图7 第120炮地震记录Fig.7 The record of shot 120
图8是基于爆炸发射理论模拟的剖面，由于爆炸反射界面成像原理没有时间差的关系，可认为是零炮检距剖面，考虑的是单程波，以及速度近似地认为是实际速度的一半，因此必然会存在一些误差，如图9中所出现的一些归位不是很准确及反射波不收敛的地方，当然这也和模型的精度降低有关（抽稀为25m×25m）。但是偏移使Marmousi模型的三个大断裂都基本归位，两个背斜构造成像清楚，2500m深度处的油水接触界面成像也非常清楚，表明NSPS叠前深度偏移算法的准确性和可靠性。
图8 基于爆炸反射理论的有限差分正演剖面Fig.8 Finite difference forward section based on exploded reflecting theory
4 认识和讨论
1）NSPS深度偏移算法对于层位的归位还是很准确，虽然较之普通的相移算法成像效率要慢，但在保证精确度和准确性的前提下，效率也是可以接受的。
2）从文中两个模型最终的偏移效果来看，对于倾角比较大的地质界面（断层面、背斜的两个斜面）而言，该算法还不能使绕射波完全地收敛，当然这可能是由算法本身的假设条件引起[12]。
3）Marmousi模型是工业公认的叠前深度偏移算法的验证模型，由于其复杂程度接近实际的地质结构，因此利用该模型来验证偏移算法就显得很有必要；而本文的算法是基于规则网格（25×25）下的偏移算法，对于非规则网格Marmousi模型（网格为12.5×4）还有待进一步的研究和提高。
4）相速度、波数和空间采样间距是相关的，当初值为一般函数时，由于其含有各种波数成分，它们将以各种不同的相速度传播，所以波形会不断的散开形成重影，这是差分所引起的频散；在用波动方程作模拟和偏移处理时，一些同向轴由于相速度和群速度的不一致在传播过程中就会产生这种畸变和重影，还有一些是由波动方程本身近似和空间采样率所引起的。
图9 NSPS偏移剖面Fig.9 Migration section of NSPS
参考文献
[1]Gazdag,J.Wave equation migration with the phase shift method.Geophysics,1978,43:1342～1351
[2]程玖兵，王华忠，马在田.带误差补偿的有限差分法叠前深度偏移方法.石油地球物理勘探，2001,36（4）:408～413
[3]程玖兵，王华忠，于富文等.波动方程共炮检距道集叠前深度偏移.石油地球物理勘探，2001,36（5）:526～532
[4]杨辉，高亮，刘洪等.微机群并行实现Marmousi模型叠前深度偏移.地球物理学进展，2001,16（3）:68～75
[5]马在田.高阶方程偏移的分裂算法.地球物理学报，1983,26（4）:377～388
[6]Jeno Gazdag,Piero Sguazzero.Migration of seismic data by phase shift plusinterpolation.Geophysics,1984,49:124～131
[7]Margrave,Gary F.Theory of nonstationary linear filtering in the Fouier domain with application to time variant filtering.Geophysics,1998,63:244～259
[8]Margrave,Gary F.Wavefield extrapolation by nonstationary phase shift.Geophysics,1999,64:1067～1078
[9]Margrave,Gary F.RobertJ.Ferguson.An explicit,symmetric wavefield extrapolator for depth migration.69th Ann Internat Mtg Soc Expl-Geophys,Expanded Abstract[C],1999,1461～1464
[10]RobertJ.Ferguson,Gary F.Margrave.Prestack depth migration by symmetric nonstationary phase shift.CREWES research report,1999,11:1～17
[11]Pann,K.,Shin,Y.Aclass of convo1utional time-varying filters.Geophysics,1976,41:28～43
[12]贺振华等.反射地震资料偏移处理与反演方法.重庆大学出版社，1989
The Study on Pre-stack Depth Migration Based on Nonstationary Phase Shift
Wang Weiguo,Qian Rongyi
（Guangzhou Marine Geological Survey,Guangzhou 510760）
Abstract:In theory,the general phase shift migration algorithm is fast,high stability,and has no grid dispersion.But in complex geological structure,e to the huge diversity of lateral velocity,migration often can not get the desired results.However,the nonstationary phase shift algorithm is based on the theory of nonstationary filter,and it adds a disturbance of phase shift operator in horizontally（called nonstationary phase shift operator）,resulting in migration well adapted when the lateral velocity changes.The algorithm has achieved a very good application effect with pre-stack depth migration in common velocity model varied with lateral velocity and the Marmousi model.
Key words:Nonstationary phase shift;Pre-stack Depth Migration;Velocity model