vani661-CSDN博客

原创使用matlab画风矢量场与填色图（不…

以下操作均以NCEP的风矢量数据作为演示，有nc数据的信息浏览、读出等操作%显示nc数据的信息ncdisp('E:\NCEP\uwnd.mon.mean.nc');%ncinfo我试了下也可以，返回的矩阵中可以分别查看到信息；最笨就是Grads里打开，但是能解决问题都不错***********以下为系统返回数据****************Source: E:\N

2015-11-05 09:56:11 5893 2

原创 [matlab+ArcGIS]在matlab中使用shp…

%将ArcMAP中导出的数据在matlab中画出来。下一步学习想学习怎么将matlab处理好的数据直接叠加在图中%参考http://blog.csdn.net/rumswell/article/details/14648647与%画shp文件zj='E:\\MAP.shp';%文件地址r=shaperead(zj);%文件读入mapshow(r)%画图%axesm 设置投影方式，详询h

2015-11-05 09:56:08 5756 2

原创 [matlab]计算事件连续发生的频数、…

事件定义为持续或短暂的满足某一阈值条件。例如极端高温为定义为日最高气温大于95%分位数的记录，若有一天的气温条件满足上式，则算是一次极端高温事件，若连续10天都满足大于95%分位的条件，则这10天算是一次连续的极端高温事件。程序的目的就是算出一个区域内事件的发生频率、持续时间和起始日期。%数据预处理。满足事件条件的标为1，不满足的标为0,存入列向量a中%%连续事件个数统计，输入量a为列向量，

2015-11-05 09:56:06 1328 2

原创 [matlab]计算月均值：日数据的华丽…

如果你要写一个程序将日值转成月值你要考虑什么：1、平年、闰年；2、大月小月然后就是某年某月的值在这个数组的位置我竟然愚蠢的想了一天这种基础性的问题。%计算月均值，原数据：per-day%1、制作一个日历，把原数据放进去（一般原数据自己会带日期）%2、搜索数据位置，做均值%timemake()function time=timemake(k1,k2)%k1开始年份，k2结束年份；c

2015-11-05 09:56:03 6657

原创自己动手检查网络—ping[DOS]

终于得到网络中心分配的IP地址，可以不用再厚脸皮蹭别人的无线网了~【感动泪目】这两天逐渐对计算机和网络感兴趣起来，在做各方面的尝试，今天清早就遇到个连接不上Internet的尴尬局面，打电话给中心工作人员，人家教我PING的命令~记录下：问题：连接不上Internet，网络诊断结果是DNS没有反应；操作：cmd->ping 10.xx.xxx.xxx(本地局域网的代理服务器（这个是不是DH

2015-11-05 09:56:00 239

原创 Fortran批量输出ctl文件格式

*对于文件名有规律可循且文件格式整齐统一的grd来说，一个一个复制粘贴太丧尽天良地浪费时间了！*虽然我好久没用Fortran，可能编程的时间比粘贴的时间长。。。。学习嘛~熟悉了就好了*....别问我为什么没用matlab，我刚刚意识到自己傻逼了*原帖地址http://bbs.06climate.com/forum.php?mod=viewthread&tid=7267!try to wr

2015-11-05 09:55:58 840

原创提取nc资料的逐月数据[grads]

**原帖地址http://bbs.06climate.com/forum.php?mod=viewthread&tid=20302**我还在尝试搞清楚出来的grd在matlab里怎么排布。。。想不明白0 0!'reinit'*先设置各月的名称mon.1=Jan;mon.2=Feb;mon.3=Mar;mon.4=Apr;mon.5=May;mon.6=Jun;mon.7=Jul;mon

2015-11-05 09:55:55 3309

转载一个Fortran 读文件的模板

留一个~即使不一定会用。。。原文地址：读文件的模板">一个Fortran 读文件的模板作者：欢乐小汪不知道这种文件存哪里比较好每次都要找一段时间干脆存网上好了program GaugeSelectimplicit noneCHARACTER(len=80):: FilenameCHARACTER(len=80):: arg1 ,arg2CHARACTER(len=80):: bu

2015-11-05 09:55:53 327

原创拒绝用FIND!直接找到矩阵中的位置…

%程序第一部分：做时间序列（亮点是直接列出每月日数）%程序第二部分：做日数据的月均值（亮点是不用if和find而直接反出数据位置，最大程度优化程序）%function sum=m_sum(data,y1,y2)size_d=size(data);�ta为没有日期的数据，列为站点或者矩阵点的数据，y1=1961;y2=2013;DATA(size_d(1),size_d(2)+3)=0;D

2015-11-05 09:55:50 368

原创插值！[matlab]

%反距离插值方法，使用ArcGIS很容易就办到，但是具体算法最近才学会%师兄的程序，我要像师兄看齐~增大阅读面，多做数据，并且善于分析，我要多出成果！！！function Z=IDW2015(data,station,lon,lat,e,num)�ta数据 n3*n4（时间*站点数）%station data 对应坐标的经纬度 n4*2 第一列纬度第二列经度%lon lat 需要插值

2015-11-05 09:55:47 4043

原创做空间插值[matlab]

%计算各站点冬、夏两季降水值%对降水值插值，转换到格点数据上，便于做空间相关clear;clc;%1、读取所有站点资料File=uigetfile({'*.txt';'*.*'},'MultiSelect','on');%得到有数据的站点名info=load('E:\info.txt');DATA=[];L=length(File);summer=[];winter=[];for

2015-11-05 09:55:45 3326 1

原创个人逼叨+输出格式+画图设置（乱放…

赶紧写下来，一会又找不到重新写，会哭死的。对于已经处理好的in/out水汽文件（到现在没找到原来处理的程序放哪里了，找到以后补上来），计算月平均值并且输出%计算水汽通量的月平均、趋势、及年变化%clc;clear;path=['E：\ERA'];lat(54,13)=0;lon(54,13)=0;a(3,13)=0;c(12)=0;%不计算经纬净通量就不用lat、lon数组moun

2015-11-05 09:55:42 219

原创月值降水数据格式的提取和转存[mat…

好久没有写新东西了，感觉以前做的数据写不成文章就这么扔掉好可惜。╭(╯^╰)╮，一定要折腾一篇出来。降水值要乘以0.1我都忘记了！clear;clc;%改程序用于处理原始降水值数据，将其转存为站点名对应的月降水值。info.txt中储存了站点名对应的经纬度%starYear为转存的起始年份，endYear为终止年份，本程序一律要求起始年为1960，原始数据不达标则抛弃File=uiget

2015-11-05 09:55:40 851

原创将NC文件转存为GRD，方便Fortran和…

*直接下载的nc文件为28*19的格点数据，1948-2013年月数据*实现的时间起始为1960年*it=145year=1960while(yearmon=1while(monit=(year-1948)*12+mon*while(it'reinit''sdfopen D:\NCEP\vapor\pr_wtr.mon.mean.nc''set dfile 1''set

2015-11-05 09:55:37 4566 2

原创将格点数据格式转化成适用于ArcMap…

%第一次尝试将自己写的小程序总结起来，都是方便查找和自我提高用%原始截断数据u（7021*15），第一行为经度，第一列为纬度，45年月值数据%结果数据矩阵d(year*mon+2,x*y);y=13;x=14;%x、y为经度、纬度格点数function(Data)=ReShape(u,year,x,y)data=u(2:end,2:end);%仅存放数据，不存经纬度信息for n=0:

2015-11-05 09:55:34 691