天外来客的花田

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

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

任务描述：输出一个带有字符串和数字的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

任务描述：获取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年

在热图中，坐标轴标签用*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

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

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

# -*- 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/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/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

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

#-*- 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

%% 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

找到指定字符串位置：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);...

任务描述：对时间序列进行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

任务描述：有站点坐标，制作.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

任务描述：有站点坐标，制作.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

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

参考：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 (

任务描述：添加单选按钮控件后，如何根据用户不同的选项返回不同的值？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