matlab 关于冻土分布灰度图转变彩色底图处理问题

冬青

现有一份冻土分布图的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);




拜托了！！！