learn – zhangxudong

November 15, 2024November 15, 2024learn

根据太阳辐射来计算日照时数DLc 和根据比湿计算相对湿度RHu

% 张旭东更新20240724 % DLc——平均日照时数[hour] % Rn——冠层表面净辐射[MJ/m^2/day] % Rns——净短波辐射[MJ/m^2/day] % Rnl——净长波辐射[MJ/m^2/day] % Rs——入射的太阳辐射[MJ/m^2/day] % Ra24——碧空太阳总辐射[MJ/m^2/day] % latitude——地理纬度转换为[度] % elevation——海拔高度，[m] % RHmean——平均相对湿度，% […]

April 30, 2024April 30, 2024learn

Protected: 土壤含水率插值-采用时间序列插值

There is no excerpt because this is a protected post.

April 26, 2024April 26, 2024learn

Protected: 读取CLDAS数据某个站点的太阳短波辐射

There is no excerpt because this is a protected post.

July 1, 2023July 1, 2023learn

用python去除txt文件中的换行符

将需要去除换行的txt文本文档，命名为 text.txt 运行py文件后，会生成 “去除换行后的文本.txt” def qchh(): f1 = open(‘./text.txt’,’r’,encoding=’utf-8′) #将需要去除换行的txt文本重命名text.txt f2 = open(‘./去除换行后的文本.txt’,’w’,encoding=’utf-8′) #去除换行之后，生成新的txt文本 –“去除换行后的文本.txt” for a in f1.readlines(): data1=a.strip(‘\n’) […]

November 10, 2020November 10, 2020learn

矩阵有些行需要重复

%% 矩阵a每一行重复n次 a为已知矩阵 m=size(a,1); % a有多少行？ b=[] n=5; for i=1:m b0=repmat(a(i,:),n,1); b=[b;b0]; end %% 矩阵a每m行重复n次，a为已知矩阵 b=[]; n=5; m=5; k=size(a,1)/m; %循环k次 […]

April 30, 2020May 3, 2020learn

ArcGIS常用工具：union、merge、append、dissolve区别

转 https://blog.csdn.net/sisiiiiiii/article/details/39178913 回归 https://ww2.mathworks.cn/help/stats/fitlm.html?#opennewwindow https://ww2.mathworks.cn/help/stats/regression-and-anova.html?category=regression-and-anova&s_tid=CRUX_gn_documentation_regression-and-anova

April 17, 2020April 17, 2020learn

根据栅格图相同Value，生成矢量边界(ROI)，剪切并保留原边界矢量的属性表

（作者：秋成） 1、作物分布图为栅格img格式，载入ArcGis后可以看到只有3种像元值，其中0为水面，1为村庄，3为农田。想把农田对应的区域生成ROI矢量图。 2、把img转换为tiff格式，采用右键dataExport data 3、根据剪裁后的tiff格式文件提取农田分布图。打开Editor–>Start Editing；打开tiff格式作物分布图的属性表，选择目标行（本例中选择Value为2的作物像元）；接着打开conversion–>from raster–>raster to polygon把选中的tiff部分转换为shp格式。提取出的农田部分生成shp文件如下： 4、如想把农田按照行政去分开，进行如下操作。依次单击Analysis–>Extract–>clip 框1为要剪切的图形（带有地区信息的边界文件，如可以分割农田的县级行政边界）框2为模板（即步骤3生成的shp图）框3字段根据个人需要进行选择。单击OK，这样就完成了！

April 17, 2020April 17, 2020learn

CAD转shp并直接在Arcgis下配准

（作者：秋成）已知配准点的坐标，目标是将无坐标和无投影的CAD文件转成ArcGIS对象。一、dwg转shp格式 1）先将配准点（已知坐标点/控制点）放入CAD图层，以便后期使用。如下图所示，十字相交的地方就是配准点。 2）使用ArcToolboxàData ManagementàFeaturesàFeatures to xx (要转换的形状点Point、线Line、面Polygon) ，导出位置选择默认（Arcgis默认的数据库） 3）导出成shpfile格式并且选择导出位置和名称。 4）在Arc catalog下面对需要配准的矢量图设置投影方式（右键点击该图，在shapefile属性表的shape字段属性中，定义相应的投影方式。二、在ArcGis下直接配准 1) 通过excel或记事本建立一个文本文件，输入几个配准点的坐标值，x表示经度，y表示纬度，可以是公里网格坐标，和上面的定义的投影方式对应，保存； 2) 单击Editor –>Start Editing，编辑的前面导出的shpfile文件； […]

August 28, 2019learn

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 […]

October 30, 2018October 30, 2018learn

用矩阵散点做等高线图和梯度场图

matlab提供了可以根据矩阵做等高线的命令contour，和求梯度的函数gradient，可以根据梯度用quiver命令将梯度做成箭头图。但要注意一些参数的设置。一般情况下，例如： x=-2:0.2:2; y=-2:0.2:2; [X,Y]=meshgrid(x,y); Z=X.^3+Y.^2; contour(X,Y,Z); [DX,DY]=gradient(Z,0.2,0.2); hold on; quiver(X,Y,DX,DY); 作图如下：设置一下不同的坐标，纵坐标的范围为 [-20 20]，间距为2。 x=-2:0.2:2; y=-20:2:20; [X,Y] = meshgrid(x,y); […]