- 积分
- 7
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2022-11-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
现有一份冻土分布图的tif 灰度图,请问需要将它叠加底图绘图是要先划分冻土分布区域再绘图吗,请问有相关程序推荐吗?目前提取tif图的经纬度数据程序如下,运行失败,请问可以帮我看一下问题,或者有什么更快捷的提取程序吗?
info = geotiffinfo('1.tif');%读取tif信息info
huidu= imread('1.tif');%读取信息、数据
[A,refmat,bbox] = geotiffread('1.tif');
[row,col]=size(huidu,1,2);
imshow('1.tif')%出图
%imfinfo、imread、imshow读取信息、数据,显示图像,并显示信息
% % 获取各格点经纬度
% for i = 1 : row
% for j = 1 : col
% [x,y] = pix2map(info.RefMatrix, i, j);
% [lat(i,j),lon(i,j)] = projinv(info, x,y);
% end
% end
% 并行计算获取各格点经纬度
lat = zeros(size(fenbu,1),size(fenbu,2));
lon = zeros(size(fenbu,1),size(fenbu,2));
%
% size(fenbu,1)
% size(fenbu,2)
parfor i = 1 : size(fenbu,1)%行固定,循环计算一行的
vlat = zeros(1, size(fenbu,2));%行固定
vlon = zeros(1, size(fenbu,2));%行固定
for j = 1 : size(fenbu,2)%固定列。子循环
[x,y] = pix2map(info.RefMatrix, i, j);
[vlat(j),vlon(j)] = projinv(info, x,y);
end
lat(i,:) = vlat;%输出固定行的经纬度
lon(i,:) = vlon;
end
scatter(lat,lon);
拜托了!!!
|
|