功能类似拉斐尔帕拉西奥斯的deg2utm.m但使用公式从“地图投影 - 一个工作手册”,由J.P.斯奈德(1987)在南半球更好的精度。
这第二个版本允许迫使UTM区和半球中输入,但是仍然可以像以前的版本,在其中函数查找相应的区域使用。
例子:
纬度= [48.866667;34.05;-36.85];
经度= [2.333056;-118.25;174.783333];
%数据区:
[X1,Y1,utmzone1,utmhemi1] = wgs2utm(纬度,经度)
%被迫区:
[X2,Y2,utmzone2,utmhemi2] = wgs2utm(纬度,经度,60, 'S')
亚历山大schimel(2020)。wgs2utm(版本2)(//www.tianjin-qmedu.com/matlabcentral/fileexchange/14804-wgs2utm-version-2),MATLAB中央文件交换。检索。
1.2.0.0 | 完成说明和示例 |
|
1.1.0.0 | 继请求,这第二个版本允许迫使UTM输入 |
约翰·麦克道尔(查看资料)
亚历山大schimel(查看资料)
喜的答案迟到抱歉。我不使用此帐户,也不支持此功能了。万博1manbetx
我怕我没有写一个伴侣功能开启横轴墨卡托坐标转换为纬度/只要我并不需要一个。通过J.P.斯奈德(1987) - 但是,公式应该在“一个工作手册地图投影”中找到。
对于所有其他问题,随意修改功能,按您的需求。交换输入订单,或修改输出UTM区域的格式应该很容易。
否则,请尝试最近开发https://au.mathworks.com/matlabcentral/fileexchange/45699-ll2utm-and-utm2ll
wave_buoys(查看资料)
你好,如果我的UTM数据点位于什么在两个不同的区域?如何含有两种不同UTM区为长LAT这组数据转换?谢谢
Mahsa khodaee(查看资料)
艾哈迈德Pishehvari(查看资料)
Uilke Stelwagen(查看资料)
两点意见:
1)在哪里伴侣功能utm2wgs?
2)utmzone是数字而来自utmzone例如deg2utm是一个数字加上一个字符?
库尔特(查看资料)
伟大的作品,并给出了在南半球比别人对文件交换更好的结果。我想我会做不同的是输入的顺序的唯一的事情 - 因为输出为X,则Y,我会作出的投入经度纬度,然后,但显然不是一个大问题。感谢您提供
帕维尔(查看资料)
对不起,真的很简单的问题。我需要这个代码的工作,但我真的在MATLAB不是很好。
你能解释一下如何把投入?
很抱歉的问题
贡萨洛(查看资料)
非常感谢你。此功能完美的我。我工作的公共汽车在伦敦的地理定位和功能有没有问题,并帮助我很多在我的研究。
埃斯特万(查看资料)
插口(查看资料)
真的感谢。
中号Sohrabinia(查看资料)
嘿再次,一直在研究这个问题,因为我写的评论,似乎Matlab的内置功能,需要更多的调整(我使用ArcGIS检查这个太),你的函数工作accuratesly与亚米级精度。请忽略我以前的评论。
中号Sohrabinia(查看资料)
你好,感谢你的代码,
我只是做了你的代码和Matlab的之间的快速内置函数基于LatLon计算UTM值,这里的结果:
MATLAB:
mstruct = defaultm( 'UTM');
mstruct.zone = utmzone(-43.5,172.5);
mstruct = defaultm(mstruct);
[X,Y] = mfwdtran(mstruct,-43.5,172.5)
结果:621272.828768645 5182478.29916538
您的代码:
[X,Y,zone1的,utmhemi1] = wgs2utm(-43.5,172.5)
结果:621267.241003134 5182565.80639284
(有在y中几乎百米差异)。
顺便说一句,你已经使用utmzone作为输出变量名与冲突Matlab的内置函数的名字,我花了一段时间,为什么我想你的功能后,越来越怪异的错误(从你的函数输出工作停止utmzone功能)。祝好运
S K'(查看资料)
亚历山大schimel(查看资料)
仅供参考。我刚才提交的wgs2utm允许迫使UTM区和半球中输入,按要求的新版本。
您仍然可以使用它像旧的,虽然,即没有强迫,让功能选择相应的区域。新功能应该与使用旧版本的任何代码兼容。检查功能的帮助上使用的更多信息。
安德鲁·莉丝(查看资料)
谢谢你的提示 - 我会用在需要时。
亚历山大schimel(查看资料)
评价关于强制输出UTM区。
两个参数更改为每个UTM区:参考经度和假北。在当前的代码,这些值是从在线63和68的输入数据计算值:
Lon0 =地板(Lon./6)* 6 + 3。以度%参考经度
FN =(纬度<0)* 10000000。%的假北距
为了强制输出UTM区,一个可以迫使这两个参数,以他们的所需的值。如果您知道自己的参考经度,只是通过这个值替换Lon0。如果您知道区号II,由6 * II-183替换Lon0
设置FN为0北半球,或千万的南半球。
由于安德鲁,我理解你的关切。我一直在努力,从数据的区域区两者之间,我的代码没有处理这个问题。
我彻底修改了我的原代码,以便能够迫使UTM区也:
- 允许转换到本地横轴墨卡托投影(NZ突起是预设的)
- 要计算网格收敛(与定向数据打交道时真北和网格北,绝对需要的信息之间的差异)
我建议如果可能的话,你移动到局部突起。如果你的数据在多个UTM区跨越,你的变形将是巨大的。
如果你直接寄给我,我可以为你提供新的功能。我不知道搬到网上来这里的那一刻。对于其他人谁是有兴趣相同的注释。
非常感谢。任何机会,以力输出到一个单一的区?与模型我的工作,无法处理多个UTM区,但有时我的数据集跨越多个区域,所以这将是很好的重新工程一切到指定的区域(认识到这将降低精度的数据点外的正确区)。很不错不管。非常感谢!
用户,你可以请写信审查,以评估在你感兴趣的领域转型的品质?谢谢