`vxvy2uv`文档

vxvy2uv是南极的一部分映射工具Matlab(格林et al ., 2017)。点击这里在AMT功能的完整列表。

这个函数变换极性球面向量组件地理坐标(纬向和经向)向量组件。

内容

(U, V) = vxvy2uv(纬度、经度、vx V) [U, V] = vxvy2uv (x, y, vx, V)

(U, V) = vxvy2uv(纬度、经度、vx V)变换极性球面向量组件vx,v地理位置的引用纬度和朗地理地带性和子午组件。

(U, V) = vxvy2uv (x, y, vx, V)变换极性球面向量组件vx,v极地立体位置的引用x和y地理地带性和子午组件。极地球面坐标自动确定如果前两个输入任何值vxvy2uv超过正常的地理坐标的值。

考虑一个向量场描绘从左向右运动标准极极射赤面投影图。使其x分量1.5级无处不在,我们将添加一个轻微downward-dipping y分量,只是为了让事情变得有趣。建立网格,我们psgrid创建一个6000公里宽的网格分辨率250公里,集中在南极:

(x, y) = psgrid (“南极”、6000、250、“xy”);vx = 1.5 * 1(大小(x));v = -0.3 * 1(大小(x));

在笛卡尔坐标向量场是什么样子。我们会为背景:覆盖Bedmap2接地线

箭袋(x, y, vx, v)在bedmap2 (“gl”,“xy”)

关闭

我们可以把这个向量场这样的纬向和经向组件:

[u, v] = vxvy2uv (x, y, vx, v);

让我们显示地图上的纬向和经向组件。我们必须使用ps2ll让我们xy的地理坐标网格。初始化一个地图antmap,画一个海岸线补丁对象上下文,和情节AMT的向量场函数quivermc:

(纬度、经度)= ps2ll (x, y);antmap bedmap2“patchcoast”quivermc(纬度,经度,u, v)

看起来合适。有一个有趣的向量接近南极因为小错误角度变得疯疯癫癫。

关闭

让我们使用措施的工具箱看冰冰川运动在复苏。首先定义一个500公里宽的电网集中在复苏的冰川。使网格分辨率3公里。

(纬度、经度)= psgrid (“复苏冰川”500年,3);

的measures_interp函数可以返回冰速度组件在每个纬度、经度网格点:

vx = measures_interp (“vx”纬度,经度);v = measures_interp (“v”纬度,经度);

在笛卡尔坐标系中,向量场是这样的:

(x, y) = ll2ps(纬度、经度);箭袋(x, y, vx, v)

这是有点难以解析使用Matlab的箭袋功能,但是如果你放大箭头。例如,

轴([-640322 -553225 758771 839473)

关闭

您可能希望冰运动的纬向和经向组件。做一下,显示向量上斜坡的形象。当调用quivermc下面我color-scaling箭头使用我最喜欢的一个colormap发电机,brewermap,我不想给太多的小箭头,所以我会downsample领域的20%“密度”,20选择:

(u, v) = vxvy2uv(纬度、经度、vx v);mapzoom (“复苏冰川”)斜坡quivermc(纬度,经度,u, v,“colormap”brewermap (256,“红色”),“密度”,20)

如果这个函数或其他任何部分的南极地图工具对你很有用,请引用描述了AMT的纸。

格林,c。,Gwyther, D. E., & Blankenship, D. D. Antarctic Mapping Tools for Matlab.电脑与地球科学。104 (2017)pp.151 - 157。doi: 10.1016 / j.cageo.2016.08.003。

的vxvy2uv功能和支持文档是乍得a·格万博1manbetx林写的德克萨斯大学奥斯汀分校的地球物理研究所(UTIG), 2015年9月。