给矩阵添加地理信息

使用MATLAB实现

im = cell2mat(struct2cell(load('im.mat')));
%a=cell2mat(struct2cell(load('im.mat')))
%imshow(a,[])
[Data, R]= geotiffread('origin_image.tif');  % origin_image.tif 为带有地理信息的图像
info=geotiffinfo('origin_image.tif');
geotiffwrite('new_image.tif',im, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);   % new_image.tif 为写入地理信息后的新图像，new_image_matrix为待写入地理信息图像的像素值矩阵

需要保证原tif图和待写入的矩阵大小一样

否则会报错： 错误使用 geotiffwrite (line 257) R.RasterSize is inconsistent with rasterSize