天外来客Jonty-CSDN博客

需要用到函数randnhisthistogramcumsum%% 绘制频率分布直方图及频率分布曲线%% 修改日期 2022/5/17clcclearclose all%% 生成随机数y = randn(100,100); %生成100*100的正态分布随机数y = y(:);%% 计算频数及其频率x = -3.95:0.1:3.95;

2022-05-17 23:24:48 17223 1

原创 MATLAB error:14090086:SSL routines:ssl3_get_server_certificate:certificate verify failed

用MATLAB的websave下数据时遇到了上述报错，在外网找到了答案：加一句options.CertificateFilename=('');即可！后记：写博客的初衷是分享经验，同时是算是自己对思路和代码的整理，方便日后处理数据，应该可以帮到很多人。我已免费分享我的心得，如果看官还有其他问题的，那么：知识付费，我的时间和经验正好可以解决你的问题。咨询问题请添加QQ：8193693542022年4月20日...

2022-04-20 15:22:55 2571 1

原创 MATLAB 字符串和数字的矩阵输出为txt

任务描述：输出一个带有字符串和数字的txt，方便arcgis读取%%字符串和数字的矩阵输出为txt% 修改日期 2022/4/16clcclearclose alltitle = {'name','lat','lon','sf'}; % 表头txt = [72.1513 -122.0189 0.4456 69.5040 -149.5686 0.3343 69.4690

2022-04-16 11:00:34 3895

原创 MATLAB 快速标注序列日期

任务描述：获取1900年1月1日至2039年12月31日的日期排列（年-月-日）需要用到datenum函数，将日期转化为日期序数（从公元0年1月1日起的天数）a = datenum(0,1,1)b = datenum(1900,1,1)c = datenum(2039,12,31)d = datenum(2039,12,31)-datenum(1900,1,1)+1%% 运行结果% a = 1% b = 693962% c = 745095% d = 511341900-2039年

2022-03-25 13:46:02 908

原创 MATLAB 更改或隐藏热图（heatmap）的坐标轴刻度

在热图中，坐标轴标签用*data表示：figure;my_matrix = rand(3);heatmap(my_matrix, 'Colormap', parula(3), 'ColorbarVisible', 'on', 'XLabel', 'Time', 'YLabel', 'September')ax = gca;ax.XData = ["Hello" "World" "Thursday"]代码和解决思路来自外网，链接如下：How to change XTick Labels in

2022-03-24 15:50:28 10707 1

原创 ArcGIS 坡度分析及坡度集中在80-90间的原因

坡度分析步骤：工具箱–>空间分析–>表面分析–>坡度坡度集中在80-90间的原因：表面Z单位用其他单位而非地面x，y单位表示需要使用z因子对坡度分析进行校正，计算最终输出表面时，将用 z 因子乘以输入表面的 z 值。如果 x,y 单位和 z 单位采用相同的测量单位，则 z 因子为 1。这是默认值。如果 x,y 单位和 z 单位采用不同的测量单位，则必须将 z 因子设置为适当的因子，否则会得到错误的结果。例如，如果 z 单位是英尺而 x,y 单位是米，则应使用 z 因子 0.3048

2022-01-15 10:39:03 2879

原创皮尔逊Ⅲ型曲线的离均系数Φ值表免费百度云链接

本着分享原则，提供百度云盘免费下载链接如下：链接：https://pan.baidu.com/s/1bVU4rg7WJ-u5SkL3ypG67w 提取码：yzvk

2022-01-03 20:35:42 652

原创 Python(arcpy) 批量shp转raster

# -*- coding: utf-8 -*-'''批量实现shp转raster2021/12/24'''import arcpyimport osimport re#检查Arcgis的权限arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput = 1#设置了输出范围及输出坐标系，这三行也可注释掉arcpy.env.overwriteOutput = Truearcpy.env.extent = arcpy.Ext

2021-12-24 14:31:35 1200 2

原创 MATLAB 绘制空间分布图，颜色代表数值（替代ArcGIS）

%% 修改日期 2021/12/9clcclearclose all%% DEM空间分布图figure(1)worldmap([35 75],[-150 -50]); %研究区范围[DEM,R] = geotiffread('F:\DEM_SRTM\North America\Canada_USA_mianland.tif'); %加载带有坐标信息的图DEM = double(DEM); %需要将

2021-12-10 13:53:31 7592 5

原创 MATLAB 绘制点的地理空间分布，并用点的颜色或大小代表数值

%% 修改日期 2021/12/8clcclearclose all%% 测试数据，第一列代表纬度，第二列经度，第三列则是点的值txt = [46.75296619 -69.1022775 0.312909752 46.80735808 -69.06131914 0.391153485 46.68438136 -68.40644856 -0.524917185 47.68233472 -67.83951645 -0.362872481

2021-12-09 21:14:51 3953 9

原创 MATLAB 笔记整理（三）数据处理

--------------------------------三维矩阵运算A_std = std(A,[],3); %求第三维上的标准差

2021-12-06 13:07:52 279

原创 Python(arcpy) 批量计算流域面积

#-*- coding: UTF-8 -*-'''计算子流域面积'''import arcpyimport globimport osimport xlrdimport rearcpy.env.workspace = "F:WATA_polygon_shp\\"outpath = "F:WATA_shp_with_area_field\\"shps = arcpy.ListFiles("*.shp")for shp in shps: oldname = [substr.sta

2021-11-29 22:20:56 1374

转载 MATLAB MK趋势检验

%% MK趋势检验function MKtrend = MKtrend(time_series)n = length(time_series);for i = 1:n s(i) = sum(time_series(i:end)>time_series(i)) - sum(time_series(i:end)<time_series(i));endS = sum(s);%s是正态分布，均值为零，方差：varS = n*(n-1)*(2*n+5)/18;%当n>1

2021-10-26 21:29:54 5969 1

原创 MATLAB 查找指定字符或字符串所在位置

找到指定字符串位置：k = strfind(str,pat)示例：str = 'Find the starting indices of substrings in a character vector';k = strfind(str,'in')k = 1×5 2 15 19 36 41拆解文件名（除去后缀）str = filename;pat = '_';newStr = extractBefore(str,pat);...

2021-09-27 20:23:19 21162

原创 MATLAB Mann-Kendall突变检验 (mk突变检验)

任务描述：对时间序列进行MK突变检验：将MK突变检验的代码封装为函数，直接调用即可，代码如下：%% MK突变检验function [UF,UB] = MKbreak(time_series)n = length(time_series);%% ---------------------------------正序列计算% 定义累计量序列s，长度=n，初始值=0；% s(k)记录的是第i个时刻(其中i<=k)，其数值大于j时刻(其中j<=i)数值个数的累计数s = zeros

2021-09-16 16:57:57 31295 22

原创 MATLAB 根据经纬坐标生成点的shp文件或流域的边界shp文件

任务描述：有站点坐标，制作.shp文件用于arcgis后续处理[txt,~] = xlsread('C:\Users\DELL\Desktop\Works\数据整理.xlsx',12);location = txt(:,1:3);for i = 1:length(txt) %前面的部分为必要的字段 Point.Geometry = 'Point';%元素类型，Point，Line，Polygon等 Point.X = location(i,1);%坐标X 对应经度 P

2021-08-29 14:30:30 3839 2

原创 SPSS 定性变量（名义变量、字符串）如何导入做ANOVA方差分析

任务描述：做方差分析时，发现无法在因子中导入名义变量原因：SPSS的规定因素变量只能是数值变量，不能是字符变量。需要使用Transform（转换）将字符变量转换为数值变量。选择该列，点击Transform（转换） > Auto-Recode（自动重新编码）将该列字符串变量转换为数值变量即可导入Hope this may bring you some inspirations!...

2021-08-23 15:04:21 8614 1

转载 MATLAB Unsupported shape type PolygonM (type code = 25)

参考：https://gis.stackexchange.com/questions/40613/importing-shapefile-in-matlabMatlab使用shaperead读取shp文件时，出现错误Unsupported shape type PolyLineZ (type code = 13). headerTypeCode = readHeaderTypeCode(shpFileId,callingFcn);Unsupported shape type PolygonM (

2021-08-03 15:34:25 550

原创 MATLAB app designer 制作单选按钮的方法

任务描述：添加单选按钮控件后，如何根据用户不同的选项返回不同的值？selectedButton = app.ButtonGroup.SelectedObject;switch selectedButton.Text case '允许' A= 1; case '禁止' A= 0;endassignin('base','A',A); %将值返回到MATLAB工作区Hope this would bring you some inspiration

2021-06-09 11:10:55 4545

转载 MATLAB app designer/GUI清除工作区（workspace）的用法

近日在调试Matlba GUI程序时，因为不想退出程序后手动Clear All来清理，又需要在过程中对WorkSpace进行清理，否则会引用之前的结果导致错误，找了很多资料，国内的论坛什么的都说用Clear，Clear All这种，根本不行，会报错，因为过程中Clear将会导致定义的handles丢失，最后在Matlab的官方论坛中找到了答案：Clear Global Workspace from GUI Pushbutton“global (main)workspace” I assume/guess

2021-06-09 09:36:49 6780

原创 MATLAB app designer 如何设置背景图片

MATLAB app designer 本身是不带设置背景图片功能的，但是我们能够通过使用app.UIAxes或者app.button来实现，灵感来自于国外的网友：首先要注意的是，在app designer中发生控件重叠覆盖时，不同控件之间的图层层级关系是：新生成的控件会显示在旧空间上边，所以！！！要保证坐标区（UIAxes）是最早插入的控件,这样坐标区就能显示在最下一层然后读入图片到app.UIAxes中就行了~~例如：Hope this would bring you some inspira

2021-06-08 23:03:28 9179 1

空空如也

空空如也