matlab大地坐标与经纬度转换,如何把经纬度转化为大地坐标
txt文件数据修改成我的数据后,替换了文件名,程序运行后出现了两个问题:
一:弹出图像框 但没有内容 图像框里是空白的,这事为什么?
二:输出的coordinate.txt文件里数据怎样处理精度问题,就像您贴出来的转换后的坐标似地能保留到小数点后很多位?请看附录二
附录一:我的程序
L=dlmread('hangji01.txt');
axesm utm
Z=utmzone(L(1,1),L(1,2));
setm(gca,'zone',Z)
h = getm(gca)
R=zeros(size(L));
for i=1:length(L)
[x,y]= mfwdtran(h,L(i,1),L(i,2));
R(i,:)=[x;y];
end
dlmwrite('coordinate.txt',R)
h =
mapprojection: 'utm'
zone: '50S'
angleunits: 'degrees'
aspect: 'normal'
falsenorthing: 0
falseeasting: 500000
fixedorient: []
geoid: [6.3774e+006 0.0817]
maplatlimit: [32 40]
maplonlimit: [114 120]
mapparallels: []
nparallels: 0
origin: [0 117 0]
scalefactor: 0.9996
trimlat: [-80 84]
trimlon: [-180 180]
frame: 'off'
ffill: 100
fedgecolor: [0 0 0]
ffacecolor: 'none'
flatlimit: [32 40]
flinewidth: 2
flonlimit: [-3 3]
grid: 'off'
galtitude: Inf
gcolor: [0 0 0]
glinestyle: ':'
glinewidth: 0.5000
mlineexception: []
mlinefill: 100
mlinelimit: []
mlinelocation: 1
mlinevisible: 'on'
plineexception: []
plinefill: 100
plinelimit: []
plinelocation: 1
plinevisible: 'on'
fontangle: 'normal'
fontcolor: [0 0 0]
fontname: 'helvetica'
fontsize: 9
fontunits: 'points'
fontweight: 'normal'
labelformat: 'compass'
labelrotation: 'off'
labelunits: 'degrees'
meridianlabel: 'off'
mlabellocation: 1
mlabelparallel: 40
mlabelround: 0
parallellabel: 'off'
plabellocation: 1
plabelmeridian: 114
plabelround: 0
附录二:coordinate.txt文件里部分数据
7.0444e+005,3.5555e+006
7.0444e+005,3.5555e+006
7.0443e+005,3.5555e+006
7.0443e+005,3.5555e+006
7.0443e+005,3.5555e+006
