- 博客(47)
- 问答 (2)
- 收藏
- 关注
RSS订阅原创 决定系数与方差解释率
如何评估回归线究竟在以一个变量预测另一个变量的工作中做得怎么样，可以将已解释的方差除以总方差，这个比值就叫做。简言之，“解释”有驱动和预测的两层含义。“解释”在回归模型中对于。也称为“方差解释率”。
2024-03-16 23:15:20
307
1
原创 Python中的多线程和多进程
1. **并发性需求：** 如果你的任务是计算密集型的，即需要大量的 CPU 计算，那么使用 `multiprocessing` 可能更合适，因为它可以利用多个 CPU 核心并行执行任务。综上所述，如果你的任务是计算密集型的，且你希望实现真正的并行计算，那么使用 `multiprocessing` 是一个不错的选择。2. **资源消耗：** `multiprocessing` 创建的进程是独立的，每个进程都有自己的内存空间，因此在资源消耗方面较高。**以下是多线程和多进程的一些关键区别：**
2023-12-01 21:29:05
191
原创 Basemap的投影
1. Cylindrical Equidistant Projection（等经纬度投影）: 'cyl'这只是一些常见的投影方式，还有其他投影方式可供选择。6. Albers Equal Area Projection（阿尔伯斯等面积投影）: 'aea'8. Stereographic Projection（立体投影）: 'stere'2. Plate Carrée Projection（等面积投影）: 'cea'13. Hammer Projection（哈默投影）: 'hammer'
2023-11-22 15:43:38
368
原创 SMOS数据处理，投影变换，‘EPSG:6933‘转为‘EPSG:4326‘
在处理SMOS数据时，遇到了读取nc数据并存为tif后，影像投影无法改变，因此全球数据无法重叠。源数据的投影为'EPSG:6933'，希望转为'EPSG:4326'。
2023-10-11 14:30:57
337
原创 from rasterio. _base import gdal_version, DLL load failed: 找不到指定的模块”
解决办法：直接在pycharm中安装gdal==3.4.3，rasterio==1.3.8, rioxarray==0.15.0, xarray==2023.8.0.问题2：in rasterio.crs.CRS.from_user_input rasterio.errors.CRSError: CRS is invalid。data.rio.write_crs("EPSG:4326", inplace=True) #要点，要加inplace=True。问题3：reproject。
2023-09-28 17:44:29
537
原创 回归方程显著性检验
3.回归分析的目的：根据已知的资料或数据，找出变量之间的关系表达式(找到回归方程)，用自变量的已知值去推测因变量的值或范围(进行预测)，实际上是研究因果关系。2.相关分析是根据统计数据，通过计算分析变量之间关系的方向和紧密程度，而不能说明变量之间相互关系的具体形式，无法从一个变量的变化来推测另一个变量的变化情况。回归分析能够确切说明变量之间相互关系的具体形式，可以通过一个相关的数学表达式，从一个变量的变化来推测另一个变量的变化情况，使估计和预测成为可能。有显著的线性关系，并不意味着每一个解释变量分别对。
2023-06-29 09:35:33
1958
原创 归纳和演绎
归纳是对现状的总结，是对一些现象的高度凝练，并非把所有现状都进行举证，也并非包括所有的现状和现象，特别是对于不完全归纳，尽管前提是真实的，推理过程也是正确的，但结论不一定完全真实，因为可能有黑天鹅现象的出现，可能会出现一些意外。推理的起点是从认清现状开始的，通过认清现状，得到可能的结果，再从得到的结果中找出共有的属性，再由这些共同属性得出一般规则、一般规律。归纳是由个别事物到一般事物的推理，其结论所判定的知识范围，已经完全超出了前提中包含的知识范围，它是这些现象或事物的高度凝练、高度总结。
2023-05-22 14:55:15
409
原创 pyCM3安装过程的问题
删除卸载所有theano， pip install theano-pymc==1.1.0。降级theano-pymc，pip install theano-pymc==1.1.2。安装pyMC3==3.11.0， 一定要装theano-pymc，而不是theano。将numpy降级到：pip install numpy==1.23.1。将numpy降级到：pip install numpy==1.21.6。降级pymc3，pip install pymc3==3.11.2。# 定义观测数据分布。
2023-05-09 11:20:25
1485
2
原创 pyMC3中The “freeze_support()“ line can be omitted if the program
windows系统运行pyMC3的时候出现错误：
2023-05-03 20:44:37
129
原创 arcgis上下标等的标注，制图时可用
表示相对于常规行间距的调整（单位为磅）；0 磅表示没有调整。表示相对于常规字符间距的调整百分比；0% 表示没有调整。表示相对于常规字符宽度的调整百分比；0% 表示没有调整。表示词间的间距百分比；100% 表示常规间距。设置名称、大小和比例，或全部设置。缺失的颜色属性假定为 0。缺失的颜色属性假定为 0。
2023-02-24 15:33:44
294
原创 多重比较中字母标记
在实际应用时，一般以大写凡具不同标记字母的即为差异显著字母 A.B.C……表示α=0.05显著水平。4. 再以标有b的最大平均数为标准，与以下各未标记的平均数比，凡不显著的继续标以字母b，直至某一个与之相差显著的平均数则标以字母c；2. 将该平均数依次和其以下各平均数相比，凡差异不显著的都标字母a，直至某一个与之相差显著的平均数则标以字母b；3. 再以该标有b的平均数为标准，与上方各个比它大的平均数比，凡不显著的也一律标以字母b；这样各平均数间，凡有一个标记相同字母的即为差异不显著。
2023-01-01 22:11:21
3198
原创 关于python numpy中元素个数和值的求和
a3=np.sum(aa[aa>0]) #大于0的元素值求和。a2=np.sum(aa > 0) #大于0的元素个数求和。在数组统计中，有数组元素值的求和，数组元素个数的求和。a1=np.sum(aa) #整个数组值的求和。
2022-12-13 20:19:06
1439
原创 Windows下python安装GDAL出现错误ImportError: DLL load failed
运行后GDAL出现 ImportError: DLL load failed的错误。python版本，3.8.15。解决方法：下载另外版本安装，
2022-12-01 01:40:50
519
原创 基于python与GDAL获取遥感影像
data = band.ReadAsArray( xoff, yoff, xsize, ysize, t_xsize, t_ysize ) //获取影像数据。band_type = fh.GetRasterBand(1).DataType //数据类型。geotransform = fh.GetGeoTransform() //文件几何信息。ysize = ysize = fh.RasterYSize //文件列数。band = fh.GetRasterBand(1) //获取波段。
2022-11-30 02:21:00
1704
原创 生态系统模型
有些过程无法用理论分析方法导出其模型，但可通过试验测定数据，经过数理统计法求得各变量之间的函数关系，称为统计模型。但由于实际的化工过程比较复杂，影响因素很多，纯用理论方法描叙过程往往是不可能，因而纯理论模型的应用是很有限的。生态系统过程模型：根据生态系统的生理生态学特性，结合影响生态系统过程的观测指标，提出的能够反映生态系统过程的机制模型。白箱模型在控制论中指的是不仅反映输入—输出关系，而且也反映过程的状态的一类模型。得到表示过程各有关变量之间的物理数学关系，这种纯粹从基本理论出发，用。
2022-11-30 02:20:19
1176
原创 广义线性混合模型
GLM一般是指 generalized linear model ，也就是广义线性模型；而非 general linear model，也就是一般线性模型；而GLMM （generalized linear mixed model）是广义线性混合模型。广义线性模型GLM很简单，举个例子，药物的疗效和服用药物的剂量有关。这个相关性可能是多种多样的，可能是简单线性关系（发烧时吃一片药退烧0.1度，两片药退烧0.2度，以此类推；
2022-11-30 02:18:52
2732
1
原创 生态系统过程模型
有些过程无法用理论分析方法导出其模型，但可通过试验测定数据，经过数理统计法求得各变量之间的函数关系，称为统计模型。但由于实际的化工过程比较复杂，影响因素很多，纯用理论方法描叙过程往往是不可能，因而纯理论模型的应用是很有限的。生态系统过程模型：根据生态系统的生理生态学特性，结合影响生态系统过程的观测指标，提出的能够反映生态系统过程的机制模型。白箱模型在控制论中指的是不仅反映输入—输出关系，而且也反映过程的状态的一类模型。得到表示过程各有关变量之间的物理数学关系，这种纯粹从基本理论出发，用。
2022-11-30 02:16:37
605
原创 ArcGIS中对矢量进行投影变换发生形变
此时需要对该矢量文件进行增密处理，以达到投影时不发生形变。具体操作为Editing Tools——>Density。
2022-09-29 14:19:52
637
原创 正反馈效应
随着更多的冰山流向海洋，冰盖将会变小，冰面下降，温度更高，融化的速度更快。所谓的正反馈效应就是某人做了符合他人价值观，让他人感到高兴的、兴奋的事情，并受到夸奖、鼓励，进而做事人就会继续努力的把这件事情做好，而且会越做越好。海冰可以反射80%的阳光，当海冰融化后，颜色较深的海洋会吸收更多的太阳光，从而加速海洋的暖化，而海洋的暖化又反过来使海冰融化更快。同样，海洋的暖化也会启动另一个正反馈效应，那就是随着海洋温度的升高，大量的海底甲烷水合物会释放到大气中，从而加快暖化速度。
2022-09-13 12:30:07
436
原创 IDL fix（x）,floor（x）, round(n), ceil(n) 的区别 (
1）fix(n)的意义是取小于n的整数(是向零点舍入的意思是往零的方向上靠)，这是一类应用在整数取值上的函数，就如同以前我们所研究的求整问题，例如，fix(!2）round(n)的意思是纯粹的四舍五入，意思与我们以前数学中的四舍五入是一样的！3）ceil(n)的意思是向正方向舍入，这一点要与floor(n)向负方向舍入一起讲，能够更好的区分，举例说明吧。
2022-09-13 12:27:37
960
原创 将图像DN值定标维热辐射强度之后,可用Planck函数求解出星上亮度温度
对于第 i = 32 波段, 为 K 32 , 1 = 474 . 6 84780 W?式中, K i 1和 K i 2 是常量，对于第 i = 31 波段, 分别为 K 31 , 1 = 729 .541636 W?32波段亮温=1196.978785 /alog(1+474.684780/b32)31波段亮温=1304.413871/alog(1+729.541636/b31)
2022-09-13 12:26:24
981
原创 实验设计-假重复
这个试验设计在基本逻辑上存在漏洞:任何两个试验小区——即使在不进行任何处理的情况下——土壤微生物生物量就一定是不同的，而且只要取样量(n)足够大这种差异就一定能达到统计上显著的水平[2]。认为这是有n个重复的试验而对两个区的调查结果进行统计检验。可做如下表述：假重复是在采用推论统计学方法通过试验数据检验处理效果过程中，所设计的试验要么是处理没有重复（可能有样本重复），要么是重复不具统计学上的独立性。简单伪重复是Hurlbert指出的若干种伪重复情形中的最简单的一种，但是在中国生态学试验中也是最常见的一种。
2022-09-13 12:25:37
1134
1
原创 碳十四断代法
绝对不是某些反对进化论的宗教分子理解的只能推断的范围等于半衰期。如果需要更古老的测定，完全可以用其他半衰期更长的放射性元素。生物体在活着的时候会因呼吸、进食等不断的从外界摄入碳十四，最终体内碳十四与碳十二的比值会达到与环境一致 （该比值基本不变），当生物体死亡时，碳十四的摄入停止，之后因遗体中碳十四的衰变而使遗体中的碳十四与碳十二的比值发生变化，通过测定碳十四与碳十二的比值就可以测定该生物的死亡年代。生物体的每克碳内含有大约500亿个碳14原子,其中每分钟大约有10个碳14原子衰变。
2022-09-13 12:22:13
3496
原创 什么是关键科学问题
7、关键技术是为关键科学问题服务的，是为解决关键科学问题而提出来的具体的方法，关键技术是看得、见摸得着的，申请书里这一节应该写的越具体越好。11、盖一栋房子，打地基是其内容之一，但关键问题可能是由于盖房子的地方下面土质太软，因此要解决的关键问题是打地基时如何解决土质太软的问题。3、凝练出重要的、有意义的、必须要解决的、不得不做的、非做不可的、能开天辟地继往开来造福人类的、但是还没搞定的“关键科学问题”。6、是科学问题，不是具体技术，而是说此技术所代表的科学问题，你拥有了这个技术，可以解决什么关键科学问题。
2022-09-13 12:20:51
5040
原创 pycharm快捷键设置
改变选中代码的背景颜色，File——Settings——Editor——Color Scheme——General——Editor——Guides——Selection background.11. 在pycharm中书写程序的时候，有时候不注意格式，会造成程序中有很多的警告，这时候，pycharm已经为我们开发了一个快速规整程序格式的入口，快捷键为Ctrl + Alt + L。Alt + J选择相同单词，可以同时编辑多行，Alt + Shift + J取消选择。5. 程序选中高亮设置。
2022-09-13 12:17:14
3006
原创 python安装工具包时候提示没有权限
'c:\\program files\\windowsapps\\pythonsoftwarefoundation.python.3.7_3.7.1218.0_x64__qbz5n2kfra8p0\\lib\\site-packages\\pip-19.0.3.dist-info\\entry_points.txt'Could not install packag es due to an EnvironmentError: [WinError 5] 拒绝访问。在install后面加--user即可。
2022-09-13 12:10:04
1443
原创 相关性和因果性
通常计量上的因果指的是格兰杰因果，定义如下：若在包含了变量X、Y的过去信息的条件下，对变量Y的预测效果要优于只单独由Y的过去信息对Y进行的预测效果，即变量X有助于解释变量Y的将来变化，则认为变量X是引致变量Y的格兰杰原因。因果性是用来进行预测的，即宏观世界中，当X为Y的因，那么X在时间线上要早于Y发生，那么加入X变量之后，对Y的预测优于单纯的ARIMA预测，在符合一致性和无偏性时，那么我们认为X格兰杰引起了Y（即X为Y的格兰杰因）。所以，在写论文的时候，一般都写A和B正相关/负相关，置信水平是多少之类的。.
2022-08-12 09:37:05
1664
原创 太阳短波辐射换算
2、w/m2：瓦每平方米，国际单位制导出单位之一，可表示单位时间通过单位土壤截面的热量，是热通量密度的单位；也可表示投射到单位接收面积的辐射照量，是辐射照度的单位。1、MJ/m2：兆焦每平方米，表示单位面积所受的太阳能辐射能量，是太阳能辐射量的单位。所以：1 w/m2=31.536 MJ/m2/year。月辐射乘以30d，年辐射约等于乘以365d....
2022-08-12 09:36:14
11137
1
原创 python3安装pyemd
3. 报错：ModuleNotFoundError: No module named 'PyEMD', 哪里报错，就打开那个文件，将其中的PyEMD全部替换为pyemd, 如果visualisation报错，则改为from pyemd.visualisation import Visualisation，程序改为from pyemd import EMD, Visualisation即可。...
2022-08-12 09:35:15
1122
原创 协方差的定义
如果有X,Y两个变量，每个时刻的“X值与其均值之差”乘以“Y值与其均值之差”得到一个乘积，再对这每时刻的乘积求和并求出均值（其实是求“期望”，但就不引申太多新概念了，简单认为就是求均值了）可以通俗的理解为：两个变量在变化过程中是同方向变化？同向或反向程度如何？...
2022-08-12 09:34:29
370
转载 vscode自定义的模块，在不同目录调用时，一直报错Nomodulenamedxxx
查了好多资料，一直没解决，最后发现，为什么pycharm可以，但vs code不行，主要是因为pycharm创建时会创建__init__.py文件，所以调用不报错，在vscode文件夹下创建__init__.py就不会报错了。从pycharm换到vscode，被一个错误纠结好久，就是自定义的模块，在不同目录相互调用时，一直报错"No module named 'xxx' "...
2022-08-12 08:58:13
442
原创 python多元线性回归
random_state 在需要设置random_state的地方给其赋一个值，当多次运行此段代码能够得到完全一样的结果，别人运行此代码也可以复现你的过程。若不设置此参数则会随机选择一个种子，执行结果也会因此而不同了。虽然可以对random_state进行调参，但是调参后在训练集上表现好的模型未必在陌生训练集上表现好，所以一般会随便选取一个random_state的值作为参数。线性回归模型的评估可以根据RMSE来判断，最小二乘法与梯度下降法的区别：实现方法和结果不同：最小二乘法是直接对求导找出全局最小，是非
2022-06-29 09:57:57
474
转载 Anaconda建立新的环境，出现CondaHTTPError:HTTP000CONNECTIONFAILEDforurl......解决过程
2020.3.7准备scrapy，使用anaconda创建一个新的环境，执行“conda create -n scrapyEnv python=3.6”，结果出现了“CondaHTTPError: HTTP 000 CONNECTION FAILED for url ”。以下是我的解决步骤。一、电脑的环境操作系统：windows10 64位Anaconda版本：conda 4.8.2二、解决步骤1、打开Power Shell，执行“conda config --remove-key channels”命令
2022-06-29 09:56:22
1122
原创 MODIS原网址无法打开
MODIS原网址无法打开1. 登录2. 查找数据：https://ladsweb.modaps.eosdis.nasa.gov/search/3. 提交申请：https://ladsweb.modaps.eosdis.nasa.gov/archive/orders/501609021/4. 看帮助提示，利用wget下载 wget -e robots=off -m -np -R .html,.tmp -nH --cut-dirs=3 "https://ladsweb.modaps.eosdis.nasa.go
2022-06-29 09:54:28
401
原创 python+Module Not Found Error:No module name dmpl_toolkits.basemap
python3.8运行出现错误：ModuleNotFoundError: No module named 'mpl_toolkits.basemap'，但是basemap是安装上的。
2022-06-29 09:52:55
545
原创 ImportError:DLL load failed while importing_datadir:找不到指定的模块。
pycharm中运行文件出现错误：ImportError: DLL load failed while importing _datadir: 找不到指定的模块。解决：在pycharm中，file——setting——add python interpreter——conda environment——existing environment—，分别配置好python3.8的编译器D:\Anaconda3\envs\Python38\python.exe和conda executable: D:\Anaco
2022-06-29 09:51:00
256
原创 以度为单位的像元面积转换为平方千米（matlab实现）
第一步：function cellwt = GetAreaWeight(varargin)if nargin == 1 deg = varargin{1}; rows=180/deg; cols=360/deg; halfSize = deg/2; tmp = ones(rows,cols,'single'); R = makerefmat(-180+halfSize, 90-halfSize, deg, -deg); EP = almanac('earth','w
2022-06-28 17:12:02
569
空空如也
python制图时，plt.show()出现图像显示缺失的情况
2022-06-29
如何多线程（多进程）加速while循环(语言-python)
2022-03-10
TA创建的收藏夹 TA关注的收藏夹
TA关注的人