CHEN_BR/BLOG

转载 机器学习常用方法网址（第一批汇集）【转】

（1）2016-09-26，30分钟学会用scikit-learn的基本回归方法（线性、决策树、SVM、KNN）和集成方法（随机森林，Adaboost和GBRT）https://blog.youkuaiyun.com/u010900574/article/details/52666291 竟然KNN这个计算效能最差的算法效果最好（2）2018-12-26，机器学习的几种方法（knn，逻辑回归，S...

原创 《白话统计》笔记 v.1

《白话统计》笔记 v.11. 变异，统计学的基础2. 概率分布2. 1 用weibull分布寻找生存规律2. 2 用logistic分布探索疾病流行规律2. 3 “普通”的正态分布2. 4 标准正态分布2. 5 小样本分布——t 分布3. 关于统计资料类型的思考3. 1 寻找cut-off值的多种方法3.1.1 根据专业和经验3.1.2 利用ROC曲线找到cut-off值（二分类结局）3.1...

原创 RFsp — Random Forest for spatial data (R tutorial) （翻译）

原文 ：https://github.com/thengl/GeoMLA摘要：本教程介绍了如何使用随机森林生成空间和时空预测(即使用随机森林从点观测生成地图)。空间自相关，特别是在交叉验证残差中仍然存在时，表明预测可能存在偏差，这是次优的。为了解释这一点，我们使用随机森林(如ranger包中实现的那样)结合到抽样位置的地理距离来拟合模型和预测值。我们描述八个典型情况下感兴趣的空间预测的应用:...

原创 关于ERA5再分析逐小时数据中的变量expver

来源 http://bbs.06climate.com/forum.php?mod=viewthread&tid=99863最近又开始使用ERA5的该数据集下载到了2021 发现出现expver变量 不是很理解它选1或者5的意思 用python尝试选了下读出来都有问题 我是想把位势和温度逐小时数据处理成日均数据 想请问一下大家在出现这个变量是怎么处理的ERA5 hourly data on pressure levels from 1979 to present2021年的数据有expver

原创 Matlab与Google Earth一起自定义M_Map岸线

来源：https://my.oschina.net/u/4579695/blog/4487173https://www.eoas.ubc.ca/~rich/map.htmlhttp://www.ngdc.noaa.gov/mgg/shorelines/data/gshhshttps://ww2.mathworks.cn/matlabcentral/fileexchange/35642-kml2structhttps://ww2.mathworks.cn/matlabcentral/fi.

原创 安装nc4

./configure --prefix=/home/theo/netcdf4/zlib-1.2.11makemake checkmake install./configure --prefix=/home/theo/netcdf4/hdf-1.8.20 --with-zlib=/home/theo/netcdf4/zlib-1.2.11 CC=icc CXX=icpc --enable-cxx CFLAGS=-O3 CXXFLAGS=-O3 FC=ifort F77=ifort F90=ifort

原创 bash: ./config.guess: /bin/sh^M: 坏的解释器: 没有那个文件或目录

[theo@localhost work]$ ./config.guessbash: ./config.guess: /bin/sh^M: 坏的解释器: 没有那个文件或目录[theo@localhost work]$ sed -i 's/\r//′urfilesed：无法读取urfile：没有那个文件或目录[theo@localhostwork]//' urfilesed：无法读取 urfile：没有那个文件或目录[theo@localhost work]//′urfilesed：无法读取urfil

转载 vmware 虚拟机共享文件夹后找不到共享目录

执行一下这个命令vmhgfs-fuse /mnt/hgfs/ 再看 mnt 目录

原创 查看sklearn中所有的模型评估指标

import sklearnsorted(sklearn.metrics.SCORERS.keys())['accuracy','adjusted_mutual_info_score','adjusted_rand_score','average_precision','balanced_accuracy','completeness_score','explained_variance','f1','f1_macro','f1_micro','f1_sample...

原创 圆形断面临界水深

% -- 计算 抛物线 形渠道的 临界水深计算clc; clear allbeta=0.001:0.01:10; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度Mk=8*(1/g*beta.^2).^(1/3); hk025=0.207*(Mk).^0.7704; % 单位 m；plot(beta,hk025,'b-','linewidth',2); hold on axis([0.001,3,0.01,1]) xlabel('无量纲.

原创 梯形断面临界水深莫洛图

% -- 计算梯形渠道的 临界水深计算clc; clear allbeta=0.01:0.01:1000; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度 % plot(lam,entak0,'k-','linewidth',2); hold on% 迭代一次% plot(lam,entak1,'r-','linewidth',2); hold on % 当 m=0 时； 矩形hk00=(1/g*beta.^2).^(1/3);plot(beta.

原创 三角形临界水深

% -- 计算 三角形 形渠道的 临界水深计算clc; clear allbeta=0.01:0.01:1000; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度 Q=beta*g^0.5; % 当 m=0.5 时； m=0.5;hk05=(2*beta.^2/m/m).^0.2; % 单位 m；plot(Q,hk05,'k-','linewidth',2); hold on % 当 m=1 时； m=1;hk10=(2*beta.^2/m/.

原创 抛物线断面临界水深莫洛图

% -- 计算 抛物线 形渠道的 临界水深计算clc; clear allbeta=0.01:0.01:1000; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度 Q=beta*g^0.5; % 当 m=0.5 时； m=0.25;hk025=(27/64*beta.^2/m).^0.25; % 单位 m；plot(Q,hk025,'b-','linewidth',2); hold on % 当 m=0.5 时； m=0.5;hk05=(.

原创 圆形断面正常水深莫洛图

% -- 计算 抛物线 形渠道的 正常水深计算clc; clear allM0=0.001:0.01:8; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度h0=(0.0104+0.3356*M0) ./(1+0.3732*M0-0.0962*M0.^2); beta=(M0/(2^2.6)).^(1/0.6); plot(beta,h0,'b-','linewidth',2); hold on axis([0.001,0.4,0.01,1]).

原创 梯形断面正常水深莫洛图

% -- 计算梯形渠道的 正常水深计算clc; clear allh0=0.001:0.01:20; % 无量纲参数的取值范围 β=αQ/b^2.5;g=9.81; % 重力加速度 % plot(lam,entak0,'k-','linewidth',2); hold on% 迭代一次% plot(lam,entak1,'r-','linewidth',2); hold on % 当 m=0 时； 矩形beta=h0.^(5/3)./(1+2*h0).^(2/3) ;plot(.

原创 通用Linux目录名称

目录 用法 / 虚拟目录的根目录。此处通常没有文件 /bin 二进制目录，存储许多GNU用户级别的实用程序 /boot 引导目录，储存引导文件 /dev 设备目录，Linux在该目录中创建设备节点 /etc 系统配置文件目录 /home 主目录，Linux在该目录中创建用户目录 /lib 库目录，存储系统和应用程序库文件 /media 媒体目录，可移动媒体设备常用的挂载点 /mnt 挂载目录，另一个可移动媒体设备

原创 libappindicator-gtk3 被 sunloginclient-10.1.1.38139-1.x86_64 需要

yum install libappindicator-gtk3

原创 /var/run/yum.pid 已被锁定

网上的解决办法：直接在终端运行 rm -f /var/run/yum.pid 将该文件删除，然后再次运行yum

原创 核密度估计

具体原理及定义：传送门https://en.wikipedia.org/wiki/Density_estimationMATLAB 代码实现如下： 1 % Kernel Density Estimation 2 % 只能处理正半轴密度 3 function [t, y_true, tt, y_KDE] = KernelDensityEstimation(x) 4 % clear 5 6 % x = px_last; 7 % x = px_last_tu; 8 %% 9 %参数初.

原创 QT5

Qt5官方在线文档：https://doc.qt.io/qt-5/Qt官方下载链接：http://download.qt.io/archive/qt/Qt官方下载链接：http://download.qt.io/new_archive/qt/Qt君整理的下载地址：http://qthub.com/download/Qt镜像列表：http://download.qt.io/static/mirrorlist/Qt清华镜像：https://mirrors.tuna.tsinghua.edu....

原创 查看所有可用的QGIS中的算法功能

# 查看可用的所有QGIS功能from processing.core.Processing import Processingfrom qgis.analysis import QgsNativeAlgorithmsProcessing.initialize()QgsApplication.processingRegistry().addProvider(QgsNativeAlgorithms())for alg in QgsApplication.processingRegistry()..

原创 matlab_cheatsheet

% make the cheatsheet for matlab% by : Peijinclear;clc;close allvis_ax = 'on';ftsz=0.85;fig_size = 800;fig_0 = figure('color','w','position',[0, 0, fig_size*1.414,fig_size]);set(fig_0,'renderer','Painters')% mainax_header = axes('position',[0...

原创 water_level_plot.m

Plotting program x = 1:1:10000; x1 = CumHOH_jaar_rot_orig; y1 = Levels_orig; x2 = PDF; y2 = v; hold on hl1 = scatter(x1,y1,'.r'); hold on; ax1 = gca; set(ax1,'XColor','k','YColor','k') ax2 = axes('Position',get(ax1,'Position'),'XAxisLocatio.

原创 water_level.m

clear all; clc; figure(1), clf; data = load('waterstand_HOH.txt'); % loading water levels afmeting = size(data); som = afmeting(:,1); % size of the dataset replace = true; % Bootstrap with replacing jaar = 98; % Numbe.

原创 横坐标

ax1.XAxis.MinorTick='on';ax1.XAxis.MinorTickValues=[256789];

原创 观测风如何变成10m风

原创 SUB-DIRECTORY STRUCTURE

SUB-DIRECTORY STRUCTURE: YYYY - Year of archived data containing archived storm data in comma delimited ATCF format. files: aBBCCYYYY.dat - guidance information bBBCCYYYY.dat - best track information fBBCCY..

原创 Mapping Toolbox 1

worldmap worldload coastlineswhos[latcells, loncells] = polysplit(coastlat, coastlon);numel(latcells)plotm(coastlat, coastlon)

原创 Mapping Toolbox 2

h = worldmap('Europe');getm(h,'MapProjection')geoshow('landareas.shp', 'FaceColor', [0.15 0.5 0.15])geoshow('worldlakes.shp', 'FaceColor', 'cyan')geoshow('worldrivers.shp', 'Color', 'blue')geoshow('worldcities.shp', 'Marker', '.',... .

原创 什么是KML？

“KML 是一种文件格式，用于在地球浏览器（例如 Google 地球、Google 地图和 Google 地图移动版）中显示地理数据。”Keyhole 标记语言 (KML) 是一种基于 XML 的文件格式，可用于表示应用程序（如ArcGIS Explorer和 Google Earth）中的地理要素。KML 允许您在地图与 globe 上绘制点、线和面，并与他人共享这些信息。您也可使用 KML 来指定文本、图片、电影或者用户单击要素后出现其他的 GIS 服务的链接。许多 KML 客户端应用程序都是免..

原创 utm with python

原创 matlab里如何实现每循环一次就把结果保存到一个新矩阵里？

每次循环产生一个一维的矩阵，希望把每次循环产生的结果依次写入一个新矩阵中，一次结果占一行。怎样创建新矩阵？怎么实现新结果不会覆盖原来的？我是初学者，求大神指教~ a = zeros(10,5);for i = 1:10 ... a(i,:) = ...endhttp://www.ilovematlab.com/thread-269562-1-1.html...

原创 Basemap画shp

# -*- coding: utf-8 -*-"""Created on Mon Nov 16 19:37:49 2020@author: chenb"""# shape格式文件的处理# Basemap处理矢量图的方法和其他库文件有很大的不同。# 先从简单的入手，用最简单的办法显示一个矢量图（山东半岛）from mpl_toolkits.basemap import Basemapimport matplotlib.pyplot as pltmap = Basemap(pro.

原创 Basemap

Unable to open boundary dataset file. Only the 'crude' and 'low',resolution datasets are installed by default.If you are requesting an, 'intermediate', 'high' or 'full'resolution dataset, you may need to download and install thosefiles separately with.

原创 【python】Basemap 报错KeyError: PROJ_LIB的解决办法

https://blog.youkuaiyun.com/weixin_45656790/article/details/108932364https://blog.youkuaiyun.com/weixin_45589713/article/details/104111271

原创 三维条形图和折线图

close all;clc;clear allf=gcf;f.Color='w';%=====================这里是绘图，数据是随意的Z=6+4*rand(10,1);Y=2+8*rand(10,1);BarHandle=bar3([1:10],Z,0.7);hold onLineHandel=plot(Y,[1:10],'Color',[0 0 0],'LineWidth',1);PointHandel=plot(Y,[1:10],'o','MarkerFaceCo...

原创 Installation instructions for VisIt

This document contains instructions for installing VisIt on Windows, MacOS X,and Unix systems. The Windows and MacOS X installation instructions are first,followed by the Unix installation instructions.Windows:========1. Download the visit3_1_3.w..

原创 Netcdf安装

Netcdf4.4的安装过程（附netcdf4.1.3的安装过程）https://blog.youkuaiyun.com/schumacher2016/article/details/80850240

原创 安装ncl

所有的NCL二进制源文件都在 Earth System Grid https://www.earthsystemgrid.org/dataset/ncl.html中列出。目前最新的版本是 NCL 6.4.0 binaries到你需要的版本，下载吧~ 如何安装NCL二进制文件？ 首先需要确定的是，你想把文件安装在哪里。 比如你想安装在/usr/local/ncl-6.4.0这个文件夹中: 那么首先你需要创建这个文件夹: mkdir /usr/local/ncl-6.4.0

原创 mt参数