问题：两张卫片，坐标系不一样，行和列不一样，无法进行波段计算。虽然不用修改坐标系，它们在ArcGIS下目视重合的（因为ArcGIS的投影和坐标是动态的）。但是因为有行和列不同，它们不能在Envi下进行波段计算。解决办法：采用Arcgis重采样。 可以在ArcGIS下将两张片子同时打开动态匹配一下，假如打开PE.tif，然后将两张片子同时进行重采样，在Export Raster Data页面栏进行如下设置：Extent→ Data Frame，Spatial Reference →Data Frame，Cell Size(cx,cy): 30 30. 然后Save. 两张片子都这样操作，保存为新的文件，然后在Envi中打开后，就会发现他们的坐标系、投影、分辨率、行列数都是一样的，这样就可以进行波段计算了。
可以实现对大图片逐像元的操作： 将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窗口中显示。
下载网址： http://print2flash.com/ 教训： 沈阳农业大学的网页要用IE浏览器啊，其他的往往不好使。
1比赛，比赛开始差点没赶上。 儿童赛很顺利地完成。 2找刘叔叔爬山。刘叔叔在某个山窝窝里超级野兽赛赛道上的第10号点位摄影，小勇士翻了两个山头重要找到了。 3寒暄过后 ，要赶路，下午2:00准时返程离开太舞，到中途大滩时，天色尚早，草原骑马，不一样的感受。小勇士从妈妈陪着，别人牵着，到自己找到平衡感，自己骑马，进步还是蛮大的。
有意思的饭和赛场 从崇礼到太舞的路太难走了，因为冬奥会，到处在施工，没多少路，开了1个小时。赛场的景色还是不错的，看看赛道。进行注册，一切都在布置。
早晨7:30经过各种准备，我们驱车离开沈阳，在暑假将要结束之时，开始了2019年暑假的斯巴达克之旅。经过六个小时的长途奔袭，到家了今天的目的地承德。今天的时间比较紧张，安排了两个游玩项目。到达承德之后，我们抓紧时间在酒店登记入住，行李都没有放到房间久全速赶往小布达拉宫。 在布达拉宫，我们去了两个地方。第一个是行宫也叫须弥寺，须弥寺是乾隆70大寿的时候建造的，具有四个第一，分别是建造时间最短，所用黄金最多，距离现在最近，参观的达赖喇嘛最多。在参观寺庙的过程中，对一些寺庙当中的一些禁忌和藏传佛教有了一定的了解。 从这个寺庙出来之后，立即赶往主要目的地小布达拉宫。小布达拉宫是乾隆为了庆祝自己的60大寿和他母亲80大寿而修建。寺庙依山而建，气势宏伟。小布达拉宫的主体建筑下三层为实心，上三层为空心，采用下宽上窄的建筑风格，墙体外中轴线上分布有六尊佛像，红墙白窗，窗子叫盲窗，盲窗又有真假之分。盲窗的功能有防风沙，装饰等作用。 站在金顶之可以鸟瞰整个承德市，远处的棒槌山也尽收眼底。 一路参观游览，到达寺庙金顶。金顶之上，整个屋顶都是金光灿灿，好多游客都开玩笑说要是带个铲子铲点黄金回去 就好了。 清风徐来，夕阳西下，屋檐下清脆悦耳的风铃声显得空灵而悠扬。 此情此景让游客们心平气和，忘却了生活中的一切纷纷扰扰。 离开布达拉宫，来到附近的一家火锅店，享用了一段饕餮大餐，回到宾馆做了短暂的休整。 晚上八点半出发观看大型实景演出-康熙大帝。这是我们第一次看实景表演，表演还是非常震撼的，网上的优酷视频网站也是有的，但比实景差远了。
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 另外，关于日期和时间序列的有个有趣的事儿 […]
今天在作图的时候无意发现了一个问题，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，原来居然一直没有发现这样的问题，汗。。。
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