defaultm
创建或重置地图投影结构
说明
mstruct = defaultm(mstructIn)mstructIn。该函数更新 origin、flatlimit、flonlimit、maplatlimit 和 maplonlimit 字段，使其彼此兼容，并与 mapprojection 字段兼容。当结构表示 UTM 或 UPS 坐标系时，该函数还会更新 zone 字段。
注意
要设置地图投影结构，必须两次使用 defaultm 函数。这些步骤展示了如何设置地图投影结构。
使用语法
defaultm(projid)创建具有默认值的地图投影结构。通过为
origin、maplatlim和maplonlim等字段指定数值来自定义地图投影结构。使用语法
defaultm(mstructIn)更新地图投影结构中的空字段和依赖字段。
示例
创建一个基于墨卡托投影的地图投影结构。默认情况下，许多字段都是空的。
mstruct = defaultm("mercator")mstruct = struct with fields:
mapprojection: 'mercator'
zone: []
angleunits: 'degrees'
aspect: 'normal'
falsenorthing: []
falseeasting: []
fixedorient: []
geoid: [1 0]
maplatlimit: []
maplonlimit: []
mapparallels: 0
nparallels: 1
origin: []
scalefactor: []
trimlat: [-86 86]
trimlon: [-180 180]
frame: []
ffill: 100
fedgecolor: 'default'
ffacecolor: 'none'
flatlimit: []
flinewidth: 2
flonlimit: []
grid: []
galtitude: Inf
gcolor: 'default'
glinestyle: ':'
glinewidth: 0.5000
mlineexception: []
mlinefill: 100
mlinelimit: []
mlinelocation: []
mlinevisible: 'on'
plineexception: []
plinefill: 100
plinelimit: []
plinelocation: []
plinevisible: 'on'
fontangle: 'normal'
fontcolor: 'default'
fontname: 'Helvetica'
fontsize: 10
fontunits: 'points'
fontweight: 'normal'
labelformat: 'compass'
labelrotation: 'off'
labelunits: []
meridianlabel: []
mlabellocation: []
mlabelparallel: []
mlabelround: 0
parallellabel: []
plabellocation: []
plabelmeridian: []
plabelround: 0
指定地图原点。然后，再次使用 defaultm 函数更新空字段以及依赖于地图原点的字段，例如地图经纬度边界。
mstruct.origin = [0 90 0]; mstruct = defaultm(mstruct)
mstruct = struct with fields:
mapprojection: 'mercator'
zone: []
angleunits: 'degrees'
aspect: 'normal'
falsenorthing: 0
falseeasting: 0
fixedorient: []
geoid: [1 0]
maplatlimit: [-86 86]
maplonlimit: [-90 270]
mapparallels: 0
nparallels: 1
origin: [0 90 0]
scalefactor: 1
trimlat: [-86 86]
trimlon: [-180 180]
frame: 'off'
ffill: 100
fedgecolor: 'default'
ffacecolor: 'none'
flatlimit: [-86 86]
flinewidth: 2
flonlimit: [-180 180]
grid: 'off'
galtitude: Inf
gcolor: 'default'
glinestyle: ':'
glinewidth: 0.5000
mlineexception: []
mlinefill: 100
mlinelimit: []
mlinelocation: 30
mlinevisible: 'on'
plineexception: []
plinefill: 100
plinelimit: []
plinelocation: 15
plinevisible: 'on'
fontangle: 'normal'
fontcolor: 'default'
fontname: 'Helvetica'
fontsize: 10
fontunits: 'points'
fontweight: 'normal'
labelformat: 'compass'
labelrotation: 'off'
labelunits: 'degrees'
meridianlabel: 'off'
mlabellocation: 30
mlabelparallel: 86
mlabelround: 0
parallellabel: 'off'
plabellocation: 15
plabelmeridian: -90
plabelround: 0
通过使用地图投影结构而非基于 axesm 的地图来投影坐标。
创建一个用于正弦投影的地图投影结构。
mstruct = defaultm("sinusoid");为地图投影结构指定地图边界和参考椭球体。再次使用 defaultm 函数，根据地图边界填充结构体的附加字段。
mstruct.maplonlimit = [-150 -30]; mstruct.geoid = referenceEllipsoid("grs80","kilometers"); mstruct = defaultm(mstruct);
加载海岸线数据并将其裁剪至地图边界。然后，使用 projfwd 函数和地图投影结构投影纬度和经度坐标。
load coastlines
[lat,lon] = maptriml(coastlat,coastlon,mstruct.maplatlimit,mstruct.maplonlimit);
[x,y] = projfwd(mstruct,lat,lon);在笛卡尔坐标区上显示投影坐标。
figure
plot(x,y)
axis equal
输入参数
输出参量
提示
默认情况下，地图投影结构中的角度值字段以度为单位。若必须使用弧度单位，请更新结构体的
angleunits字段，然后重置该结构体。此代码演示了如何将地图投影结构更新为使用弧度。mstruct = defaultm("mercator"); mstruct.angleunits = "radians"; mstruct = defaultm(mstruct);
一旦将结构更新为使用弧度，在修改角度值字段（如
origin、parallels、maplatlimit和maplonlimit）时必须使用弧度。您可以通过使用
gcm函数，从基于axesm的地图中获取地图投影结构。此代码演示了如何通过使用defaultm函数，创建与基于axesm的地图相同的地图投影结构。% Set longitude limits and define reference ellipsoid lonlim = [-150 -30]; ref = referenceEllipsoid("grs80","kilometers"); % Create map projection structure using defaultm mstruct1 = defaultm("sinusoid"); mstruct1.maplonlimit = lonlim; mstruct1.geoid = ref; mstruct1 = defaultm(mstruct1); % Create map projection structure from axesm-based map abm = axesm("sinusoid","maplonlimit",lonlim,"geoid",ref); mstruct2 = gcm(abm); f = gcf; close(f) % Compare map structures isequal(mstruct1,mstruct2)
