- 积分
- 13
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 bing 于 2015-12-22 17:22 编辑
求问:如何读取'MST2-IR1.raw'中的ir,我是将下面的这段Fortran程序改的,
parameter (m=2750)
integer*1 ir(m,m)
dimension x(m,m),y(m,m)
integer*1 ir1(1100,600)
C read IR1 channel and geolocation data ---------------
open(1,file='MST2-IR1.raw',form='binary')
read(1) ir
close(1)
open(1,file='Mtsat_IJ_LL.dat',form='binary')
read(1) x,y
close(1)
m文件程序内容:
% aera
% MST2-IR1.raw from MATLAB "ir1.m"
% Mtsat_IJ_LL.dat: Longitude and Latitude data (2750,2750)
char m = [1,1:2750];
int (sym('ir(m,m)'));
char [x(m,m),y(m,n)];
int (sym('ir1(1100,600)'));
%read IR1 channel and geolocation data ---------------
fid=fopen('MST2-IR1.raw','wb');
data=rawread('MST2-IR1.raw','ir');
fclose(fid);
dat=fopen('Mtsat_IJ_LL.dat','w');
data2=dlmread('Mtsat_IJ_LL.dat','x','y');
fclose(dat);
%--------------------------------------------------------
%interpolate to 0.1x0.1 degree grids +++++++++++++++
for j1=1:m
for i1=1:m
if (x(i1,j1)>=70)&&(x(i1,j1)<=180)&&(y(i1,j1)<=60)&&(y(i1,j1)>=0)
i=(x(i1,j1)-70)*10+1;
j=(60-y(i1,j1))*10+1;
ir1(i,j)=ir(i1,j1);
end
end
end
一直没运行正确,想知道如何正确读取'MST2-IR1.raw'中的ir,求提示
|
|