功能mapexwmsanimate (useInternet datadir)

如果~ (“使用互联网”，'var'）Undernernet = false;结尾

如果~ (“datadir”，'var') datadir = fullfile(matlabroot，“例子”，“地图”，'数据'）;结尾如果~exist（datadir，“dir”mkdir (datadir)结尾

Datafile = @(filename) fullfile(datadir,filename);

卡特里娜=wmsfind(“卡特里娜”）;谁是卡特里娜州

卡特里娜(1)

卡特里娜=精炼（卡特里娜州，“卡特里娜goes-12 * *可见*近* *动画的）;卡特里娜=精炼（卡特里娜州，'svs.gsfc.nasa.gov'，'searchfield'，“服务器URL”）;卡特里娜=卡特里娜（1）;谁是卡特里娜州

cachefile =丢失(“卡特里娜，马特”）;如果/ / / / / / / / / /如果~存在(cachefile“文件”)保存(cachefile“卡特里娜”）结尾别的缓存=负载(cachefile);卡特里娜= cache.katrina;结尾

摘要= katrina.Abstract;endOfLine =找到(isstrprop(抽象,“cntrl”),1); 抽象=抽象（1：内弗林）；numSpaces=60；尽管(~ is空(抽象))k = find(isspace(抽象));n = find(k > numSpaces,1);如果~ isempty (n)流(“%s\n”，abstract(1:k(n)) = [];别的流(“%s\n”文摘)=''；结尾结尾

katrina.Details.Attributes

katrina.details.dimension.

katrina.Details.Dimension.Extent

cachefile =丢失(“katrina.tif”）;如果Userneternet [katrinamap，r] = wmsread（katrina）;如果~存在(cachefile“文件”) geotiffwrite (cachefile katrinaMap右)结尾别的[katrinaMap，R]=readgeoraster（缓存文件）；结尾

州= shaperead (“usastatehi.shp”，“UseGeoCoords”,真正的);图usamap(卡特里娜飓风。geoshow(katrinaMap,R) geoshow(states，'facecholor'，'没有任何'）标题（{katrina.layertitle，katrina.details.diminess.default}，．..“口译员”，'没有任何'）

IEMLAYERS = WMSFIND（“IEM WMS服务”，“搜索字段”，'servertitle'）;nexrad =提炼(iemLayers,“nexrad-n0r-wmst”）;

cachefile =丢失('nexrad.mat'）;如果使用因特网nexrad = wmsupdate(nexrad);如果~存在(cachefile“文件”)保存(cachefile“nexrad”）结尾别的缓存=负载(cachefile);nexrad = cache.nexrad;结尾

nexrad.details.dimension.

nexrad.details.dimension.default.

nexradTime = [katrina.Details.Dimension.Default (1: end-1)'：00z'];

流(“% s % d \ n”，“NEXRAD层的南方纬度界限:”, nexrad.Latlim(1))流(“% s % d \ n”，卡特里娜飓风的南纬极限katrina.Latlim (1))

Latlim = katrina.latlim;lonlim = katrina.lonlim;

imageHeight = katrina.Details.Attributes.FixedHeight;imageWidth = katrina.Details.Attributes.FixedWidth;

黑色= [0 0 0];cachefile =丢失(“nexrad.tif”）;如果useInternet [nexradMap,R] = wmsread(nexrad, R)．..“Latlim”，拉特林，“Lonlim”，lonlim，“时间”nexradTime,．..“背景色”黑色“ImageFormat”，“图像/ png”，．..“ImageHeight”imageHeight,“ImageWidth”，图像宽度）；如果~存在(cachefile“文件”）地理迪迪特（CacheFile，Nexradmap，R）结尾别的[nexradMap R] = readgeoraster (cachefile);结尾

图usamap(latlim,lonlim) geoshow(nexradMap,R)'facecholor'，'没有任何'，“EdgeColor”，“白色”）标题（{nexrad.layertitle，nexradtime}，“口译员”，'没有任何'）;

阈值= 0;index = any(nexradMap > threshold, 3);索引= repmat(Index，[1 1 3]);

组合=卡特里娜飓风图；组合（指数）=uint8（下一次飓风图（指数）*255）；

图usamap(latlim,lonlim) geoshow(combination,R)'facecholor'，'没有任何'） 标题（{“12幅卡特里娜飓风图像”，．..“与NEXRAD雷达合成”, nexradTime})

程度上= katrina.Details.Dimension.Extent;削减=' / '；slashIndex=strfind（extent，slash）；startTime=extent（1:slashIndex（1）-1）；endTime=extent（slashIndex（1）+1:slashIndex（2）-1）；

连字符=“- - -”；hyphenIndex = strfind(开始时间,连字符);dayIndex = [- 1， - 2];startDay = str2double(开始时间(dayIndex));endDay = str2double (endTime (dayIndex));

katrinaTime =开始时间;

nasaServer = WebMapServer (katrina.ServerURL);iemServer = WebMapServer (nexrad.ServerURL);

katrinarequest = wmsmaprequest（katrina，nasaerver）;nexradrequest = wmsmaprequest（nexrad，iemserver）;

nexradRequest。Latlim = Latlim;nexradRequest。Lonlim = Lonlim;nexradRequest。写成BackgroundColor =黑色;nexradRequest。ImageFormat =“图像/ png”；nexradRequest。ImageHeight = ImageHeight;nexradRequest。ImageWidth = ImageWidth;

videoFilename = fullfile (pwd,“wmsanimated.avi”）;如果存在(videoFilename“文件”删除(videoFilename)结尾作家= VideoWriter (videoFilename);作家。帧速率= 1;作家。质量= 100;打开(作家)

无花果=图;Usamap (latlim,lonlim) hstates = geoshow(states, lonlim)'facecholor'，'没有任何'）;hmap = [];为了k =跑步：终点％更新时间值并为每个服务器分配时间属性。currentday = num2str（k）;KatrinAtime（DayIndex）=截止日期;nexradtime = [katrinatime（1：end-1）'：00z'];katrinaRequest。时间= katrinaTime; nexradRequest.Time = nexradTime;％从服务器（或文件）中检索Katrina的WMS映射这个时间段的％。cachefile = datafile（[“katrina_”num2str（截止日）“.tif”]）;如果useInternet katrinaMap = getMap(nasaServer, katrinaRequest.RequestURL);如果~存在(cachefile“文件”) getiffwrite (cachefile, katrinaMap, katrinaRequest.RasterRef)结尾别的katrinaMap = readgeoraster (cachefile);结尾从服务器(或文件)检索NEXRAD图像的WMS地图这个时间段的％。cachefile = datafile（['nexrad_'num2str（截止日）“.tif”]）;如果Underneternet nexradmap = getMap（IEMServer，NexradRequest.Requesturl）;如果~存在(cachefile“文件”)geotiffwrite（缓存文件、nexradMap、nexradRequest.RasterRef）结尾别的nexradMap = readgeoraster (cachefile);结尾%标识nexradMap图像中不包含%背景色。index = any(nexradMap > threshold, 3);索引= repmat(Index，[1 1 3]);%复合nexradMap与katrinaMap。组合=卡特里娜飓风图；组合（指数）=uint8（下一次飓风图（指数）*255）；%删除旧地图并显示新的合成地图。delete(hmap) hmap = geoshow(combination, katrinaRequest.RasterRef); / /删除hmapuistack (hstates“高级”） 标题（{“12幅卡特里娜飓风图像”，．..“与NEXRAD雷达合成”drawnow, nexradTime})%保存当前帧为RGB图像。currentFrame=getframe（图）；RGB=currentFrame.cdata；%为每个RGB帧创建索引图像，以便显示% GIF动画。如果k = = startDay%第一次通过循环，将RGB图像转换为%创建索引图像并将颜色映射保存到%变量,提出。使用cmap转换以后的帧。[框架,提出]= rgb2ind (RGB, 256,'点亮'）;％使用第一帧的大小和总数用于初始化动画的帧数%大到足以包含所有帧的大小。frameSize=大小（帧）；numFrames=结束日-开始日+1；动画=零（[frameSize 1 numFrames]，“喜欢”、框架);别的%使用颜色图从第一个帧转换和%将此帧转换为索引图像。帧= RGB2IND（RGB，CMAP，'点亮'）;结尾%将帧存储到GIF文件的动画数组中。framecount = k  -  startday + 1;动画（：，：1，frameCount）=帧;%将RGB帧写入AVI文件。WriteVideo（作者，RGB）;结尾

关闭（图）关闭（作者）

filename=fullfile（pwd，“wmsanimated.gif”）;如果存在（文件名，“文件”删除(文件名)结尾滞后时间= 2.0;loopCount =正;imwrite(动画,提出,文件名,．..“延迟时间”，延迟时间，“LoopCount”，loopcount）;