flatearthpoly

矢量地理数据对象包含了极必须遇到或反面子午线。在工具箱中正确显示这样的多边形,它们会导致问题的函数与笛卡尔坐标,如多边形交集和布尔操作。当这些多边形被视为笛卡尔坐标,反面子午线交叉导致一个假的线段,和多边形显示为一个补丁没有内部正确了。你可以通过使用格式化这些多边形flatearthpoly函数。

全球海岸线负载的坐标。提取第一个多边形的坐标,代表南极洲。

负载海岸线firstnan =找到(isnan (coastlat), 1,“第一”);lat = coastlat (1: firstnan);朗= coastlon (1: firstnan);

情节构成多边形边界的坐标。注意,不是封闭的边界。

情节(经度、纬度)xlim(200[-200])轴平等的

把海岸线,以便它使用平面多边形拓扑和阴谋的结果。多边形边界满足反面子午线,下降到极点,扫过的经度,并遵循反面子午线的另一边反面子午线。

[latf, lonf] = flatearthpoly(纬度、经度);图mapshow (lonf latf,“DisplayType”,“多边形”)ylim(-60年[-100])

经度坐标不连续的反面子午线混淆可以设置多边形的操作。准备使用的地理数据polybool或补丁渲染,减少多边形的反面子午线flatearthpoly函数。的flatearthpoly函数返回一个多边形的点插入遵循反面子午线杆,遍历经度杆,并返回到反面子午线沿着其他反面子午线的边缘。

创建一个拼写的地球和情节的海岸线。

axesm昊图公司setm (gca),“起源”[170]);framem在;gridm在负载海岸线plotm (coastlat coastlon)

生成一个小圆,包括北极和颜色黄色。

[latc, lonc] = scircle1(75, 45岁,30);patchm (latc lonc,“y”)

平的小圈使用flatearthpoly函数。

[latf, lonf] = flatearthpoly (latc lonc);

情节产生的切圆,你只是作为一个红色的线。

plotm (latf lonf,“米”)

生成第二个小圈子,不包括。

[latc1, lonc1] = scircle1 (170, 30);

平圆和情节红线。注意,第二个小圆,不包括钢管,沿着反面子午线剪成两半。第一个小圆是绘制的多边形平面坐标来说明其扁平的形状。的flatearthpoly函数假设多边形的内部被夷为平地的半球,包含大部分的边缘点。因此产生的多边形flatearthpoly不包括超过一个半球。

[latf1, lonf1] = flatearthpoly (latc1 lonc1);plotm (latf1 lonf1,“r”)