`isice`文档

的isice函数返回逻辑适用于所有接地冰和根据措施南极冰架边界的国际极地年2007 - 2009卫星雷达,版本2。这是逻辑的isopenocean函数。

南极边界数据集描述的国家冰雪数据中心的网站在这里。工具的概述在这个AMT插件可以被发现在这里。

内容

loni tf = isice(意大利)= isice特遣部队(xi, yi)

loni tf = isice(意大利)对所有地理位置返回true(意大利,loni)确认为接地冰或Mask_Antactica_v2冰架。数据集tif(原文如此)。

tf = isice (xi, yi)返回true在所有地点确定为在Mask_Antactica_v2接地冰。tif(原文如此),输入坐标自动确定为极地球面(真71年代lat)米如果它们的值超过正常的纬度,经度范围。

这个示例使用重力对AMT插件加载一些示例数据,但我将向您展示如何跟随没有加载数据集。你是否有重力插件,首先创建一个300公里宽的网格以猪为中心,以0.25公里分辨率:

(纬度、经度)= psgrid (“松岛冰川”300年,0.25);

如果你有重力插件,插入重力值网格是这样的:

FA = gravity_interp(纬度,经度,“自由的空气”);

如果你没有重力插件可以简单地使用随机数据是这样的:

FA =兰德(1201);

使用isice找出哪些部分的数据集在冰盖:

冰= isice(纬度、经度);

屏蔽所有非接地冰(不使用波浪号)通过设置它南:

FA(~冰)=南;

因为isice返回的逻辑不是isopenocean,我们可以完成上述步骤设置isopenoceanNaN值,这样;这只是个人喜好的问题:

FA (isopenocean(纬度、经度))=南;

情节重力数据冰盖对应:

pcolorps(纬度,经度,FA)轴紧

上下文,覆盖红色接地冰线和黑海岸:

antbounds (“gl”,“r”)antbounds (“海岸”,“k”)

如果你使用这个数据集,请引用如下:

数据集:Mouginot, J。,B. Scheuchl, and E. Rignot. 2017. MEaSUREs Antarctic Boundaries for IPY 2007-2009 from Satellite Radar, Version 2. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center.http://dx.doi.org/10.5067/AXE4121732AD。

文献引用:里戈诺特E。,S. S. Jacobs, J. Mouginot, and B. Scheuchl. 2013. Ice-shelf melting around Antarctica, Science. 341. 266-270.http://dx.doi.org/10.1126/science.1235798。

南极映射工具:格林,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。

这个函数和支持文档写的金宝app乍得a。格林德克萨斯大学的地球物理研究所(UTIG), 2016年11月。2017年5月更新。