- 博客(104)
- 收藏
- 关注
RSS订阅原创 MATLAB 如何避免复杂shp文件对inpolygon的影响
*任务描述：**当我想用inpolygon函数将属于非洲的pixel选出来时，发现因为周边小岛的影响，pixel选取有问题，如下图。首先将复杂shp文件查分成简单的shp文件，即将不相交的元素分离开手动删掉小岛（左图为删掉小岛前，右图为删小岛后）3. 用inpolygon函数进行识别，得到下图：会发现非洲大陆识别出来了，但马达加斯加岛却没了。
2025-01-22 14:42:40
406
原创 MATLAB 如何在Imagesc绘制的图上显示陆地边界？
任务描述：我有3600*1800的全球降水分布图，用imagesc出图后，如何加上陆地边界？
2025-01-05 18:13:28
256
原创 MATLAB 如何用代码将figure1中的图像背景设置为透明并调整大小后粘贴到figure2中？
如何用代码将figure1中的图像背景设置为透明并调整大小后粘贴到figure2中？
2024-10-11 15:10:23
455
原创 MATLAB 如何将boxplot的上下端横线去掉，并将上下须改为实线？
MATLAB 如何将boxplot的上下端横线去掉，并将上下须改为实线？
2023-05-23 14:36:01
943
原创 MATLAB Boxplot如何在‘compact‘及‘filled‘模式下改变箱子宽度并隐藏箱型图的中位数线？
MATLAB绘制箱型图的一些技巧
2023-04-09 22:08:44
739
翻译 MATLAB 随机森林超参数优化教程（Hyperparameters Tuning of Random Forest）
基于MATLAB的随机森林超参数优化过程
2022-09-01 11:07:48
5615
2
原创 MATLAB 绘制频率分布直方图及累积频率曲线
需要用到函数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
21204
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
2689
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
4239
原创 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
1015
原创 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
12954
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
3691
原创 皮尔逊Ⅲ型曲线的离均系数Φ值表 免费百度云链接
本着分享原则，提供百度云盘免费下载链接如下：链接：https://pan.baidu.com/s/1bVU4rg7WJ-u5SkL3ypG67w 提取码：yzvk
2022-01-03 20:35:42
832
原创 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
1284
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
8846
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
4940
9
原创 MATLAB 笔记整理（三）数据处理
--------------------------------三维矩阵运算A_std = std(A,[],3); %求第三维上的标准差
2021-12-06 13:07:52
336
原创 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
1562
原创 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
7012
5
原创 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
23897
原创 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
34485
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
4286
2
原创 SPSS 定性变量（名义变量、字符串）如何导入做ANOVA方差分析
任务描述：做方差分析时，发现无法在因子中导入名义变量原因：SPSS的规定因素变量只能是数值变量，不能是字符变量。需要使用Transform（转换）将字符变量转换为数值变量。选择该列，点击Transform（转换） > Auto-Recode（自动重新编码） 将该列字符串变量转换为数值变量即可导入Hope this may bring you some inspirations!...
2021-08-23 15:04:21
9893
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
654
原创 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
5076
转载 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
7610
原创 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
10666
1
原创 MATLAB 设置msgbox若干秒后自动关闭
任务描述：在用app designer 时要不断使用msgbox给出反映，但窗口弹多了也很烦 所以设置自动关闭h = msgbox('成功！);start(timer('timerFcn',@(obj,~)close(h),'StartDelay',1,'stopFcn',@(obj,~)delete(obj)))参考阅读：msgbox(’’)之后如何自动关闭该窗口？Hope this would bring you some inspirations!...
2021-06-08 11:03:23
2670
原创 MATLAB 判断点（坐标）是否在.shp边界（不规则多边形）范围内
任务描述：我有流域的.shp边界以及气象站的坐标，想只用MATLAB就筛选出流域内的气象站点clcclearclose allshp_path = '5876_China_62505400_Boundary_Line.shp'; %.shp文件路径 Map = shaperead(shp_path); %读取.shp文件 xb = Map.X; %边界点的经度坐标yb = Map
2021-06-07 22:32:21
5341
1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人