其實作為只把程式當工具的人來說,在寫之前一定要仔細查查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
中國聯塑推出防水產品 讓居者更加安心無憂
許多使用者都會被防水問題所困擾,新房裝修並沒有做好防水,隨之而來的問題也就愈演愈烈,家裡漏水成為 水簾洞 而漏水帶來的牆面潮濕 牆壁發霉等問題也屢見不鮮。人們漸漸意識到防水的重要性,但在琳琅滿目的防水產品市場,如何選擇一款好的防水產品成為更難的問題,中國聯塑秉承 為居者構築輕鬆生活 的品牌理念,針對...
全屋淨水哪個品牌好,濱特爾打造高階淨水產品
俗話說,病從口入。尤其疫情之下,水質的安全與健康更成為大家關注的焦點。女神們都希望擁有一瓶可改善 問題的 那麼家庭想要改善水質問題又該如何選擇,全屋淨水哪個品牌好呢?據國家衛健委組織專家研究發現,在部分省份報告的確診病例中,糞便標本檢測到新冠病毒核酸陽性,提示新冠肺炎存在 糞口傳播 的可能性。雖然現...