你忘了附上“ai_yr1.tif”。
在某个时间点提取像素的值
29日视图(30天)
显示旧的评论
我有一个大的光栅层有干旱值分配给每个像素。我也有一些位置的经度和纬度数据在美国。我想知道像素值根据我的经度和纬度位置。我第一次剪光栅对我感兴趣的领域(美国)在Arcmap我转换剪光栅.tif文件然后我下面,我将这些像素转换成纬度和经度的地方但我完成的代码有问题。我不知道该做什么。还是有一种不同的方法。有什么建议吗?
我= imread (“ai_yr1.tif”);
信息= geotiffinfo (“ai_yr1.tif”);
身高= info.Height;%的整数,表示图象的像素高度
宽度= info.Width;%的整数表示图像像素的宽度
(行,关口)= meshgrid(1:高度,1:宽度);
(遥ADlat,阿德隆)= pix2latlon(信息。RefMatrix、行峡路);
谢谢你的帮助! !
接受的答案
沃尔特·罗伯森
2017年11月29日
我= imread (“AI_clip1.tif”);
信息= geotiffinfo (“AI_clip1.tif”);
Latlon = xlsread (“Latlon.xlsx”);
身高= info.Height;%的整数,表示图象的像素高度
宽度= info.Width;%的整数表示图像像素的宽度
(行,关口)= meshgrid(1:高度,1:宽度);
(遥ADlat,阿德隆)= pix2latlon(信息。RefMatrix、行峡路);
Idub =双(I);
掩码=我= = (intmin (“int32”)+ 1);
Idub(面具)=南;
遥pixelvalues = interp2 (ADlat阿德隆,Idub。”,Latlon (: 1), Latlon (:, 2));
注意:最后的结果是南,因为那些在剪地图。
…什么
是
这些数据?当我
冲浪(Idub“edgecolor”,“没有”);colormap (parula (256))
然后它类似于山庄,但当我比较的列表最高点在美国,它似乎失踪的大部分在科罗拉多洛矶山脉。数据中的最大值点似乎有点接近圣海伦火山;数据只是
可能
对应于它的高度英寸考虑地球膨胀?但发生在落基山脉或其他高大陆山峰?
2的评论
沃尔特·罗伯森
2017年11月30日
地图的分辨率为0.00833333376795053度,我图的位置大约是710米。你指示以东大约285米的现货标有没有数据,所以它更接近于“没有数据点比数据点的东部。
你可以改变
Idub(面具)=南;
来
Idub(面具)= 0;
它使用0点标志着“没有数据”。
