NJU1healer的个人博客分享 http://blog.sciencenet.cn/u/NJU1healer

博文

[转载]MATLAB中geotiff栅格影像保存(geotiffwrite)

已有 1134 次阅读 2021-5-23 19:06 |个人分类:matlab|系统分类:科研笔记|文章来源:转载

geotiffwrite

基础语法

geotiffwrite(filename,A,R)

geotiffwrite(filename,X,cmap,R)

geotiffwrite(…,Name,Value)


输入参数变量释义:

filename:输入文件的名称和位置,文件扩展名必须是 “.tif” 或者 “.TIF”,变量类型为字符向量。


A:要导出的影像数据,为MNP的矩阵。坐标系为地理坐标,且默认为WGS-84,指定坐标可通过’GeoKeyDirectoryTag’和’CoordRefSysCode’实现。数据类型

数据类型 double | single | uint8 | uint16 | uint32 | uint64 | int8 | int16 | int32 | int64 | logical


R:空间参考信息。

空间引用信息,指定为地理栅格引用对象、引用矩阵或引用向量。

如果您正在处理投影坐标系中的图像坐标,而R是一个地图栅格引用对象或引用矩阵,则相应地指定“GeoKeyDirectoryTag”或“CoordRefSysCode”。


X:待补充

cmap:RGB图像,为M*3的矩阵


示例1

% 影像的读取

basename = 'boston_ovr';

imagefile = [basename '.jpg'];

RGB = imread(imagefile);

worldfile = getworldfilename(imagefile);

R = worldfileread(worldfile, 'geographic', size(RGB)); %读取空间参考信息

%导出影像

filename = [basename '.tif'];

geotiffwrite(filename, RGB, R)

% 显示影像

figure

usamap(RGB, R)

geoshow(filename)


示例2 读取tif处理后保存tif

% 读取tif及空间坐标信息

[A, R] = geotiffread('F:\boston.tif');

info = geotiffinfo('F:\boston.tif');

% 保存tif并使用原始数据空间坐标信息

subimage = A;

subR = R;

filename = 'boston_subimage.tif';

geotiffwrite(filename, subImage, subR,  ...

       'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);

figure

mapshow(filename);

其他参数较多,可以从原始影像中读取,也可以自定义。


‘CoordRefSysCode’ :坐标系,默认EPSG:4326(即WGS84坐标系)

使用方法:'CoordRefSysCode',4326

小知识,EPSG:4326是web map中用来表示地图坐标的代码,每个坐标都有一个对应的EPSG代码,常用的WGS84代码为EPSG:4326,伪墨卡托投影代码为EPSG:3857


‘TiffTags’:结构体,内容有以下

image.png

使用方法:tifftags = sturct('Compression','LZW')

重点!!!

      之前跑数据时,发现超过4GB的数据无法写入保存,这是因为geotiff有两种格式,一种是classictiff,另外一种是针对大数据的bigtiff。解决此办法就是加入字段       geotiffwrite(name,out,R,....,'TiffType','bigtiff','Tifftags',tifftags)

      老版本的没有这个,如果报错,请注意版本问题。目前用的19a,19b没有问题。

【参考】

MATLAB中对tif格式栅格影像读取-保存

点滴分享,福泽你我!Add oil!



http://wap.sciencenet.cn/blog-3428464-1287946.html

上一篇:[转载]激活函数总结(一)
下一篇:A类会议和期刊汇总-CCF(人工智能方向)

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-10-19 08:25

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部