Protected: 土壤含水率插值-采用时间序列插值
There is no excerpt because this is a protected post.
Category: learn
Protected: 读取CLDAS数据某个站点的太阳短波辐射
There is no excerpt because this is a protected post.
用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’) […]
矩阵有些行需要重复
%% 矩阵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次 […]
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
根据栅格图相同Value，生成矢量边界(ROI)，剪切并保留原边界矢量的属性表
（作者：秋成） 1、作物分布图为栅格img格式，载入ArcGis后可以看到只有3种像元值，其中0为水面，1为村庄，3为农田。想把农田对应的区域生成ROI矢量图。 2、把img转换为tiff格式，采用右键dataExport 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，这样就完成了！
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文件； […]
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 […]
用矩阵散点做等高线图和梯度场图
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); […]
位图格式图片如何转化为Arcgis下的shape文件
本文以jpg格式的图片为例，jpg格式是位图，它没有坐标系，但有时候又非常有必要把它做成shape文件。本文介绍两种方法，一种是通过以下6个步骤来进行的，另外一种简单一些，只有2个步骤。 1.把jpg图片放到有坐标的Arcgis里面，采用Georeferencing工具进行配准。 2.借助矢量化工具比如R2V等进行矢量化，将矢量化后的结果在CAD中对所有的线进行分层。 3.在步骤1里面选取两个特征点A，B并查询它们的坐标。 4.在CAD文件中输入A,B两个参照点的坐标，用AL对齐命令将整个目标图形旋转并缩放到正确的位置。 5.将目标图形加载到Arcgis里面，用export把相应的面或者线层导出并另存为shape文件。 注意：CAD目标文件要干净，没有多余的信息，最好分层做。采用Arcgis10以上的版本，9.3版本的有可能出问题。 6.将shape文件在catalog里面添加坐标和投影系统。就可大功告成。 另外一种简单的方法是，将jpg图片在Arcgis下配准以后，直接新建面对象，采用编辑面对象的工具Editor直接绘制边界。