量化地图扭曲点位置
的天梭和mdistort函数提供天气的视觉概述不同形式的地图投影误差。但是,有时候你需要数值误差估计在特定位置以量化或正确映射扭曲。这是有用的,例如,如果你是采样环境数据在统一的基础上在一个地图,和想知道多少面积与每个样本点,随位置和统计,将投影相关的。一旦你有了这些信息,你可以调整环境密度和其他统计收集地图投影面积的变化引起的。
映射工具箱™函数返回地址映射误差统计数据从当前或投影mstruct。的distortcalc函数计算失真数据一样mdistort提供的指定地点,但作为参数。您提供经度位置一次或向量。一般的形式是
[areascale, angdef maxscale、minscale merscale, parscale] =…distortcalc (lat, mstruct长)
然而,如果你正在评估当前映射图,省略了mstruct。你不需要指定任何返回值之后最后一个你感兴趣的。
使用distortcalc确定地图投影几何失真
下面的练习使用distortcalc计算的最大区域变形阿根廷从陆地区域的地图数据集。
读到北美和南美多边形:
美洲= shaperead (landareas。轴马力”、“UseGeoCoords”,真的,…“选择器”,{@(名字)…strcmpi(名字,{北美和南美的}),'名字'});设置空间程度(限制)映射到包含南美洲的南部,也包括接近南极地区:
mlatlim = (-72.0 - -20.0);mlonlim = (-75.0 - -50.0);[alat,阿龙]= maptriml([美洲。Lat),…(美洲。朗,mlatlim mlonlim);
创建一个墨卡托圆柱形正形投影使用这些限制,下方在指定一个格,然后情节大纲供参考:
图;axesm (“MapProjection”、“墨卡托”、“网格”,“上”,…mlatlim‘MapLatLimit’,‘MapLonLimit’, mlonlim,……“MLineLocation”5“PLineLocation”, 5) plotm (alat,阿龙,' b ')
地图是这样的:
样品每十块轮廓的点分析:
alats = alat(1:10:元素个数(alat));阿龙=阿龙(1:10:元素个数(alat));
计算面积扭曲(返回的第一个值
distortcalc)在采样点:adistort = distortcalc (alats,阿龙);
找到的范围面积变形在阿根廷(单位面积的百分比,在这种情况下,赤道):
adistortmm = [min (adistort) max (adistort)] adistortmm = 1.1790 - 2.7716
阿根廷占据中期南纬度,墨卡托被夸大了它的面积,和不同的错误明显从北到南。
删除任何
南年代的坐标数组和情节符号来表示相对扭曲成比例的圈子里,使用scatterm:nanIndex = isnan (adistort);alats (nanIndex) = [];阿龙(nanIndex) = [];adistort (nanIndex) = [];scatterm (alats,阿龙,20 * adistort,“红色”,“填满”)
结果图如下所示:
区域夸大的程度将大大大如果延伸向极远。看到大多少,得到变形面积50°S 60°S, 70°S:
= distortcalc (-50、-60) = 2.4203 = distortcalc (-60、-60) = 4 > > = distortcalc (-70、-60) = 8.5485
请注意
你只能使用
distortcalc查询在当前地图框架或位置mstruct限制。外点产生南作为一个结果。使用这种方法,您可以编写一个简单的脚本,它允许您查询地图反复确定畸变在任何期望的位置。您可以选择使用与图形光标位置
inputm。例如,[平台plon] = inputm(1)平台= -62.225 plon = -72.301 > > = distortcalc(平台、plon) = 4.6048
自然你得到的答案将取决于你的选择。使用这种方法,您可以编写一个简单的脚本,它允许您查询地图反复确定任何失真数据在任何期望的位置。
试着改变地图投影甚至方向向量投影的选择如何影响地图扭曲。为进一步的信息,请参见页面的引用distortcalc。
