读取geotiff,对该数组进行操作之后进行存储,并且将地理信息一并存储进去
[I, R] = geotiffread('H:\20170409\test\HHjiajihua\20170409_clip2_tif.tif');
info = geotiffinfo('H:\20170409\test\HHjiajihua\20170409_clip2_tif.tif');
geotiffwrite('H:\20170409\test\HHjiajihua\20170409_clip2_tif_geo.tif', finalI, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);
另一个就是没有地理信息,只有影像四个角的经纬度,使用如下
R = georasterref('RasterSize',size(db_image),'LatitudeLimits',[min(min(lat2)),max(max(lat2))],........
'LongitudeLimits',[min(min(lon2)),max(max(lon2))]);
geotiffwrite('H:\S1_ICE_Preprocess\20170409_HV_image_geo.tif',db_image,R);