高斯扩散模型，大气模型，环境规划，最后能够出图Matlab代码.rar

function arcgridwrite(fileName,X,Y,Z,varargin) if nargin < 4 error('Not enough input arguments'); end [mz,nz,zz]=size(Z); [mx]=size(X); [my]=size(Y); minX=min(X(:)); minY=min(Y(:)); dx=abs(diff(X(1,:))); dy=abs(diff(Y(:,1))); maxDiff=0.01; if any(diff(dx)>maxDiff) || any(diff(dy)>maxDiff) error('X- and Y- grid spacing should be non-varying'); else dx=dx(1); dy=dy(1); end if ischar(fileName)~=1 error('First input must be a string') elseif zz>1 error('Z must be 2-dimensional'); elseif mx~=mz error('X, Y and Z should be same size'); elseif my~=mz error('X, Y and Z should be same size'); elseif abs(dx-dy)>maxDiff; error('X- and Y- grid spacing should be equal'); end convert=0; dc=3; if isempty(varargin)~=1 [m1,n1]=size(varargin); opts={'precision';'convert'}; for i=1:n1; indi=strcmpi(varargin{i},opts); ind=find(indi==1); if isempty(ind)~=1 switch ind case 1 dc=varargin{i+1}; case 2 convert=1; end else end end end Z(isnan(Z))=-9999; Z=flipud(Z); if isnumeric(dc) dc=['%.',sprintf('%d',dc),'f']; elseif isnumeric(dc) && dc==0 dc=['%.',sprintf('%d',dc),'d']; end fid=fopen(fileName,'wt'); fprintf(fid,'%s\t','ncols'); fprintf(fid,'%d\n',nz); fprintf(fid,'%s\t','nrows'); fprintf(fid,'%d\n',mz); fprintf(fid,'%s\t','xllcorner'); fprintf(fid,[dc,'\n'],minX); fprintf(fid,'%s\t','yllcorner'); fprintf(fid,[dc,'\n'],minY); fprintf(fid,'%s\t','cellsize'); fprintf(fid,[dc,'\n'],dx); fprintf(fid,'%s\t','NODATA_value'); fprintf(fid,[dc,'\n'],-9999); [pname,fname,ext] = fileparts(fileName); if isempty(pname)==1; pname=pwd; end if strcmpi(pname(end),filesep)~=1 pname=[pname,filesep]; end h = waitbar(0,['Writing file: ',fname,ext... ', Please wait...']); for i=1:mz for j=1:nz if j==nz fprintf(fid,[dc,'\n'],Z(i,j)); else fprintf(fid,[dc,'\t'],Z(i,j)); end end waitbar(i/mz,h,['Writing file: ',[fname,ext],... sprintf(' %d%% complete...',round(i/mz*100))]) end close(h) fclose(fid); if convert==1 if ispc~=1; disp('"Convert" option only works on a pc') return end dos(['ARC ASCIIGRID ',pname,fname,ext, ' ',... pname,fname, ' FLOAT']); end