将GeoTIFF信息转换为地图投影结构
mstruct = geotiff2mstruct(项目)
mstruct = geotiff2mstruct(项目)
转换GeoTIFF投影结构,项目
,到地图投影结构,mstruct
.的长度单位mstruct
投影是米。
GeoTIFF投影结构,项目
,必须参考一个投影坐标系,用值表示“ModelTypeProjected”
在ModelType
字段。如果ModelType
的值“ModelTypeGeographic”
然后转换成地图投影结构就没有意义了,就会产生错误。
验证使用GeoTIFF投影结构的非投影坐标是否与使用地图投影结构的非投影坐标得到相同的结果。
为此,首先获取图像的GeoTIFF投影结构。调用。将角地图坐标转换为纬度和经度projinv
并指定GeoTIFF投影结构。
项目= geotiffinfo (“boston.tif”);x = proj.CornerCoords.X;y = proj.CornerCoords.Y;[latProj, lonProj] = projinv(项目,x, y);
属性从GeoTIFF投影结构获取映射投影结构geotiff2mstruct
函数。地图投影结构的长度单位为米,但地图坐标以测量尺为单位。因此,将角图坐标从测量英尺转换为米。然后,通过调用projinv
并指定地图投影结构。
mstruct = geotiff2mstruct(项目);xsf = unitsratio (“米”,“科幻小说”) * x;ysf = unitsratio (“米”,“科幻小说”) * y;[latMstruct, lonMstruct] = projinv (mstruct、xsf ysf);
验证这些值是否在一个容忍范围内。
abs(lonProj - lonMstruct) <= 1e-7 . abs(lonProj - lonMstruct) <= 1e-7 . abs(lonProj - lonMstruct) <= 1e-7 . abs(lonProj - lonMstruct
Ans = 1 1 1 1 1