Matlab对图片的读写和操作

可以实现对大图片逐像元的操作: 将n张卫星片子看做是第n层数据,构造一个x*y*n的一个三维数组: A=imread(‘TIFF_1986_NDVI.tif’); B=imread(‘TIFF_1987_NDVI.tif’); X=zeros(3984,4186,2); X(:,:,1)=A; X(:,:,2)=B; 对构造的矩阵进行操作,比如,对每层相应的像元进行排序,相减后得到Y_diff: Y=sort(X,3); Y_diff=Y(:,:,2)-Y(:,:,1); 然后将Y_diff重新保存成tif格式的图片; imwrite(Y_diff,’z.tif’); 在整个过程中,由于每张片子是3984*4186的矩阵,因此,在计算的过程中尽量不要在commond windows窗口中显示。

抓住暑假的尾巴玩(三)

1比赛,比赛开始差点没赶上。 儿童赛很顺利地完成。 2找刘叔叔爬山。刘叔叔在某个山窝窝里超级野兽赛赛道上的第10号点位摄影,小勇士翻了两个山头重要找到了。 3寒暄过后 ,要赶路,下午2:00准时返程离开太舞,到中途大滩时,天色尚早,草原骑马,不一样的感受。小勇士从妈妈陪着,别人牵着,到自己找到平衡感,自己骑马,进步还是蛮大的。

抓住暑假的尾巴玩(一20190822)

早晨7:30经过各种准备,我们驱车离开沈阳,在暑假将要结束之时,开始了2019年暑假的斯巴达克之旅。经过六个小时的长途奔袭,到家了今天的目的地承德。今天的时间比较紧张,安排了两个游玩项目。到达承德之后,我们抓紧时间在酒店登记入住,行李都没有放到房间久全速赶往小布达拉宫。 在布达拉宫,我们去了两个地方。第一个是行宫也叫须弥寺,须弥寺是乾隆70大寿的时候建造的,具有四个第一,分别是建造时间最短,所用黄金最多,距离现在最近,参观的达赖喇嘛最多。在参观寺庙的过程中,对一些寺庙当中的一些禁忌和藏传佛教有了一定的了解。 从这个寺庙出来之后,立即赶往主要目的地小布达拉宫。小布达拉宫是乾隆为了庆祝自己的60大寿和他母亲80大寿而修建。寺庙依山而建,气势宏伟。小布达拉宫的主体建筑下三层为实心,上三层为空心,采用下宽上窄的建筑风格,墙体外中轴线上分布有六尊佛像,红墙白窗,窗子叫盲窗,盲窗又有真假之分。盲窗的功能有防风沙,装饰等作用。 站在金顶之可以鸟瞰整个承德市,远处的棒槌山也尽收眼底。 一路参观游览,到达寺庙金顶。金顶之上,整个屋顶都是金光灿灿,好多游客都开玩笑说要是带个铲子铲点黄金回去 就好了。 清风徐来,夕阳西下,屋檐下清脆悦耳的风铃声显得空灵而悠扬。 此情此景让游客们心平气和,忘却了生活中的一切纷纷扰扰。 离开布达拉宫,来到附近的一家火锅店,享用了一段饕餮大餐,回到宾馆做了短暂的休整。 晚上八点半出发观看大型实景演出-康熙大帝。这是我们第一次看实景表演,表演还是非常震撼的,网上的优酷视频网站也是有的,但比实景差远了。

Hydrus日值到小时值的转换,奇妙的0.24倍

Hydrus结果的T_Level.out文件中,在没有降水和灌溉的情况下,地表的潜在通量+根系的潜在吸水速率(rTop+rRoot),应该就是ET0,但是他们之间差了0.24倍。这个0.24是从哪里出来的呢? 查阅技术手册可知p37-38,潜在通量和根系的潜在吸水速率计算的都是瞬时值。 Variations in potential evaporation and transpiration during the day can be generated with HYDRUS-1D using the assumptions that […]

每年有几次降水,几天降水?

考虑了好久的一个问题了,每年有几次降水,几天降水?可以采用diff函数来解决。 % 2014年一共多少次降水呢? 一共45次降水,共66天降水,共362.90mm data(find(data==0))=nan; % 把降水为0的值也赋值为nan a_neq_0=1-isnan(data); T_diff=diff([0;a_neq_0;0]) ind1=find(T_diff==-1) ind2=find(T_diff==1) P_n=size(ind1,1); % 降水次数 P_day_n=sum(a_neq_0); % 降水天数 参考以下链接: https://www.ilovematlab.cn/thread-544465-1-1.html 另外,关于日期和时间序列的有个有趣的事儿 […]

【转】ArcGIS地图文档(mxd)过大的问题(保存未响应)

  今天在作图的时候无意发现了一个问题,mxd文档居然有6MB多,原先以为是符号库引起的,就移除了加载的所有符号库,可文件体积并没有发生大的变化。上网搜了一下esri的论坛,才发现这个问题原来是geoprocessing引起的,只要在arcmap中执行过工具,就会把工具执行的结果保存在mxd里,因此才造成了文件大的离谱。 解决方法:点击geoprocessing菜单——results,将current session和previous session下的东西都右键remove掉(10.0版本);9.3和之前的版本是在arctoolbox窗口的 results标签页下;最后记得另存,文件便会小很多。 原文链接:http://support.esri.com/en/knowledgebase/techarticles/detail/33187 试了一下果然从原来的6MB多变成了200多k,原来居然一直没有发现这样的问题,汗。。。

如何将Google Earth的矢量化导入GIS(转)

1、新建文件夹: 4 、新建面:选择添加多边形按钮如何将Google Earth的矢量化导入GIS,然后画多边形 5、导出为KMZ格式:在参考点文件夹处右键选择“将位置另存为”,保存为参考点.KMZ 6、打开ArcMap,在ArcToolbox内选择Conversion Tools-From KML-KML To Layer,(注:goolge earth 的kmz和kml格式都可以用这个工具转换。) 在Input KML file选择“参考点.KMZ”,输出文件夹选择相应的文件夹,(注:有时候会出现自己选择的文件夹无法转换的情况,这时候可以将输出文件夹设为默认的,等到输出后,进行复制粘贴操作,转移到自己要存放的文件夹里。) http://blog.sina.com.cn/s/blog_64f9a9a60102wjrf.html