- 积分
- 224
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-5-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我从CMIP5上下载了全球海冰数据(用于练手),使用MATLAB进行读取NC数据并想将某一天的数据输出成tif格式在ARCGIS里使用,使用了如下程序:
clc;
clear all;
InPath = 'F:\paper8\';
ncdisp(strcat(InPath,'1.nc'));
source1 = strcat(InPath,'1.nc');
lon=ncread(source1,'lon');%读取经度变量
lat=ncread(source1,'lat');%读取纬度变量
time=ncread(source1,'time');%读取时间变量
sic=ncread(source1,'sic');%获取SIC数据
sic1=sic(:,:,1);
pcolor(lat,lon-180,sic1);
[x,y]=meshgrid(lon-180,lat);
phandle=pcolor(x,y,sic1');
colorbar
imwrite(sic1,'F:\paper8\1.tif','tif')
GeoRef = georasterref('RasterSize', size(sic1),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
geotiffwrite('F:\paper8\1.tif',flip(sic1'),GeoRef)
disp('finished')
但无法只能得到一个4k的tif,并且程序显示如下错误
错误: 文件:geotiffwrite.m 行:238 列:12
多个 LHS 赋值的数组不能包含 M_STRING。
跪求大佬指点,我该如何使用matlab将NC数据输出为TIF数据。
(另外:我本科是环境类的,想研究生读海洋以及气候变化的方向,现在在尝试着自己照着别人论文做,熟悉研究流程
,求前辈指点一波明路,感激不尽!)
|
|