Matlab影像像素坐标得到经纬度/经纬度转影像坐标（已知经纬度获取影像DN值）

我想要利用matlab读取geotif后，输入经纬度得到该图像的坐标以及DN值。

首先在网上能找到的是利用图像的坐标值得到经纬度坐标，就是这个博客。

那么首先是

[A,refmat,bbox] = geotiffread('H:\georef.tif');
info = geotiffinfo('H:\georef.tif');

A获取的是图像的矩阵本身，info好像获取的是各种图像的信息吧。

之后

[x,y] = pix2map(info.RefMatrix,3303, 3278);
[lat,lon] = projinv(info,x,y);

pix2map函数是将3303，3278这个行列号转化为直角坐标系的坐标，官网说的是Convert pixel coordinates to map coordinates，（话说这个to变成2也是很搞笑。。。我还以为是版本的意思）

projinv函数是直角坐标投影到地理坐标（直角坐标系坐标不是tif图像的行列号！！！）

就可以得到lat和lon了

那么已知经纬度获取DN值就是上述的逆过程，如下

[A,refmat,bbox] = geotiffread('H:\georef.tif');
info = geotiffinfo('H:\georef.tif');
[x1,y1] =projfwd(info,lat,lon);
[row, col ] =map2pix(info.RefMatrix,x1,y1);
z=A(int32(row),int32(col));

注意的是这个int32是为了将row和col是double类型转换成int类型

行和列千万别搞混，可以自己用数据验证一下!

Matlab影像像素坐标得到经纬度/经纬度转影像坐标（已知经纬度获取影像DN值）相关推荐

已知经纬度坐标求两点间距离，用python表示
已知经纬度坐标求两点间距离,用python表示已知地球上任意两点(lon1, lat1),( lon2, lat2)的经纬度坐标,求两点间的距离用haversine表示: 1.首先先将经纬度坐标的角 ...
python坐标表示_已知经纬度坐标求两点间距离，用python表示
已知经纬度坐标求两点间距离,用python表示已知地球上任意两点(lon1, lat1),( lon2, lat2)的经纬度坐标,求两点间的距离用haversine表示: 1.首先先将经纬度坐标的角 ...
ArcGIS已知经纬度范围裁剪图像的步骤
ARCGIS已知经纬度范围裁剪图像的步骤: (1)生成经纬度的面 (2)利用裁剪工具裁剪:分析工具-->提取分析-->裁剪其中第一步: (1)生成四个顶点的坐标点用xls保存(注意不能用 ...
matlab如何截取图像的中间部分_利用matlab提取并分割RGB图像中的某一个已知像素值的图像...
已知一副RGB图像中的的像素值,利用matlab将其分割出来并以二进制图像形式显示: %extract.m clear all; I=imread('new_original.png'); figur ...
ArcGis中已知经纬度坐标转任意坐标系（以点和80坐标系为讲解模板）
文章首发微信公众号:ArcGisdada 一般经纬度都是以度分秒的格式显示的,如图31-1所示. 图31-1: 我们转至WGS-84坐标需要进行十进制的转化,也就是把度分秒换算为以度为单位的一串数字. ...
wgs84 投影坐标转经纬度代码_ArcGis中已知经纬度坐标转任意坐标系（以点和80坐标系为讲解模板）...
一般经纬度都是以度分秒的格式显示的,如图31-1所示. 图31-1: 我们转至WGS-84坐标需要进行十进制的转化,也就是把度分秒换算为以度为单位的一串数字.操作步骤:打开excel>选中代表经 ...
高德地图JS - 已知经纬度来获取街道、城市、详细地址等信息
1.时候偶现你得引入谷歌API|JS var geoCoder = new AMap.Geocoder({city: ""//城市,默认:"全国" }); va ...
利用百度api判断已知经纬度的一个点是否在多边形内
见这个网址:http://api.map.baidu.com/library/GeoUtils/1.2/examples/simple.html 查看其源代码即可. 下面的代码可以判断某个点是否属于一 ...
Cesium已知经纬度获取其高度
初始数据 let arr = [[124.32441911641862,48.678295141121914],[124.32851385859364,48.67824798558877],[124. ...
c语言已知三点的坐标,求一个C语言算法：已知三角形的三个顶点的坐标,求三角形面积...
共回答了16个问题采纳率:87.5% //三个顶点坐标是6个参数. #include #include double S_triangle(double ax,double ay,double bx, ...

Matlab影像像素坐标得到经纬度/经纬度转影像坐标（已知经纬度获取影像DN值）

Matlab影像像素坐标得到经纬度/经纬度转影像坐标（已知经纬度获取影像DN值）相关推荐

最新文章

热门文章