遙感降水產品空間解析度的轉化

2021-10-01 02:14:14 字數 2412 閱讀 7636

其實作為只把程式當工具的人來說,在寫之前一定要仔細查查matlab有沒有相關的函式可以簡化工作量。

參考這個問題,與其不同的是,我是由高精度到低精度,其實是個求和的過程。

以下是程式:

a(從低精度到高精度的插值)

z1=interp2(x0,y0,z0,x1,y1,'method')

功能介紹:根據已知的資料(x0,y0,z0),用method方法進行插值,然後計算(x1,y1)對應的值z1.

用指定的演算法method 計算二維插值:

』linear』:雙線性插值演算法(預設演算法);

』nearest』:最臨近插值;

』spline』:三次樣條插值;

』cubic』:雙三次插值。

**以下:

x=1.1:0.1:1.5;

y=1.1:0.1:1.5;

[x,y]=meshgrid(x,y);

z=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];

x1=1.1:0.25:1.5;

y1=1.1:0.25:1.5;

[x1,y1]=meshgrid(x1,y1);

z1=interp2(x,y,z,x1,y1,'cubic');

b(從高精度到低精度,例如我需要從0.1°到0.25°,不是整數倍,可參考以下程式)

%遙感產品空間精度重取樣,0.1°到0.25°

%for時間序列迴圈讀入矩陣b

%b矩陣需要切掉2行1列

ini_year=2014;            %輸入計算起始年份

end_year=2018;            %輸入計算終止年份

ncols=20;                %資料集列數

nrows=12;                %資料集行數

for year=ini_year:end_year

%確定每年的天數是365或366天

if year/4==fix(year/4)

daynum=366;

else

daynum=365;

endfor day=1:daynum

%讀入初始解析度的日資料

b=load(strcat('f:\黑河+gpm(寫**中)\gpm data\gpmtxt\gpmv6',num2str(yymmdd(year,day)),'.txt'));

b=b(3:32,2:51);

a=zeros(12,20);

for i=1:12 

for j=1:20 

a=2*i+floor(i/2)-1;%奇數行

b=2*(i-1)+floor((i-1)/2)+1;%偶數行

c=2*j+floor(j/2)-1;%奇數列

d=2*(j-1)+floor((j-1)/2)+1;%偶數列

if mod(i,2) == 0

if mod(j,2) == 0 %3 b d

a(i,j)=sum(sum(b(b+1:b+2,d+1:d+2)))+0.5*(sum(sum(b(b,d+1:d+2)))+sum(sum(b(b+1:b+2,d))))+0.25*sum(sum(b(b,d)))

else %4 b c

a(i,j)=sum(sum(b(b+1:b+2,c:c+1)))+0.5*(sum(sum(b(b,c:c+1)))+sum(sum(b(b+1:b+2,c+2))))+0.25*sum(sum(b(b,c+2)))

endelse

if mod(j,2) == 0 %2 a d

a(i,j)=sum(sum(b(a:a+1,d+1:d+2)))+0.5*(sum(sum(b(a+2,d+1:d+2)))+sum(sum(b(a:a+1,d))))+0.25*sum(sum(b(a+2,d)))

else %1 a c

a(i,j)=sum(sum(b(a:a+1,c:c+1)))+0.5*(sum(sum(b(a+2,c:c+1)))+sum(sum(b(a:a+1,c+2))))+0.25*sum(sum(b(a+2,c+2)))

endend

%儲存模型計算所需解析度的列資料

fid=fopen(strcat('f:\黑河+gpm(寫**中)\gpm data\gpm025\gpmv6',num2str(yymmdd(year,day)),'.txt'),'wt');

for i=1:nrows

for j=1:ncols

if j==ncols

fprintf(fid,'%3.2f\n',a(i,j));

else

fprintf(fid,'%3.2f\t',a(i,j));

endend

endfclose(fid);

endend

endend

中國聯塑推出防水產品 讓居者更加安心無憂

許多使用者都會被防水問題所困擾,新房裝修並沒有做好防水,隨之而來的問題也就愈演愈烈,家裡漏水成為 水簾洞 而漏水帶來的牆面潮濕 牆壁發霉等問題也屢見不鮮。人們漸漸意識到防水的重要性,但在琳琅滿目的防水產品市場,如何選擇一款好的防水產品成為更難的問題,中國聯塑秉承 為居者構築輕鬆生活 的品牌理念,針對...

全屋淨水哪個品牌好,濱特爾打造高階淨水產品

俗話說,病從口入。尤其疫情之下,水質的安全與健康更成為大家關注的焦點。女神們都希望擁有一瓶可改善 問題的 那麼家庭想要改善水質問題又該如何選擇,全屋淨水哪個品牌好呢?據國家衛健委組織專家研究發現,在部分省份報告的確診病例中,糞便標本檢測到新冠病毒核酸陽性,提示新冠肺炎存在 糞口傳播 的可能性。雖然現...