samplingdemo.m

来自「生成具有空间相关性的图像。整幅图像被看作是一个离散的随机场」· M 代码 · 共 26 行

function SamplingDemo()
field=[50 50];

[nImg,uImg] = generate_correlated_image (field,0.9,1,zeros(field));

pImg = poissinv(uImg,ones(field)*5);
%pImg = wblinv(uImg,ones(field)*3,ones(field)*1);
%pImg = norminv(uImg,ones(field)*250,ones(field)*4);

surface(pImg);
axis([1 50 1 50]);

%输出数据
[r,c] = find(pImg+1);
v = pImg(:);
po=[r,c,v];

rang=[1:2500];
psa=myrandint(1,50,rang,'noreplace');

sta=[mean(v),std(v),var(v)];

xlswrite('SIMUDATA.xls', po, 'DATA', 'A1');
xlswrite('SIMUDATA.xls', po(psa,:), 'SAM', 'A1');
xlswrite('SIMUDATA.xls', sta, 'STA', 'A1');