soft ware for research
1.估计土壤参数ROSETTA Soil hydraulic parameters for the soil in the lysimeter were estimated using the ROSETTA (Schaap et al., 2001) software. […]
Month: November 2019
中国各类基金资助项目英文翻译（中英文对照）（转）
http://muchong.com/html/201801/5648771.html 不完全统计表明，国内期刊多习惯于将基金资助项目的信息作为论文首页的脚注，国外期刊则多将其作为“致谢”的一部分标注。通常情况下只需列出项目的资助号即可，不必标注具体的项目名称。如: Supported by the National Natural Science Foundation of China (Grant No. 60171009), the Major State Basic Research Development […]
怎样让图片的行和列相等（如何让卫片大小相等）
问题：两张卫片，坐标系不一样，行和列不一样，无法进行波段计算。虽然不用修改坐标系，它们在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中打开后，就会发现他们的坐标系、投影、分辨率、行列数都是一样的，这样就可以进行波段计算了。
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窗口中显示。