毛发浓密的猿工

原创 arcgis归一化栅格方法

原创 Python 风场高度场温度场图

【代码】Python 风场高度场温度场图。

原创 MaskedArray如何填补为nan

在NumPy中，MaskedArray 是一种特殊的数组，它允许你对数组中的某些元素进行掩码（mask），以便在进行计算时忽略这些元素。掩码数组常用于处理缺失值或无效数据。如果你想将MaskedArray中的掩码元素填补为nan（Not a Number），你可以使用filled方法。

原创 pre_min[0:10, 2:3] = pre和pre_min[0:10, 2] = pre区别

这个切片表达式选择了pre_min数组中的第0到第9行（包括第0行但不包括第10行），以及第3列（在Python中，索引是从0开始的，所以2:3表示第3列）。例如，如果您尝试将一个形状为(10, 1)的数组赋值给pre_min[0:10, 2]，或者将一个形状为(10,)的数组赋值给pre_min[0:10, 2:3]，您会得到一个形状不匹配的错误。对于您提到的两个切片表达式，pre_min[0:10, 2:3] 和 pre_min[0:10, 2]，它们有本质的区别，主要体现在所选的维度和形状上。

原创 2*2子图/空间插值图/ax平移

【代码】2*2子图插值图。

原创 并行下载数据

【代码】并行下载数据。

原创 pandas 按相同站号重新整合出一个dataframe

这可以通过left_on和right_on参数来实现，这两个参数分别指定左侧DataFrame和右侧DataFrame中用作合并键的列名。在这个例子中，尽管df1中的ID_num列是整数类型，而df2中的ID_str列是字符串类型，但我们仍然可以通过left_on和right_on参数将它们作为合并键。如果合并后你想要保留原始的列名，可以使用left_index=False和right_index=False参数，并可能需要在合并后重命名列。

原创 python geopandas白化/掩膜提取

注意：在运行这段代码之前，需要确保已经安装了所有必要的库，并且 test.txt 和 dishi.shp 文件存在于指定的路径下。此外，由于代码中使用了中文字符和路径，可能会遇到编码或路径问题，需要根据实际情况进行调整。以下代码主要是利用 Python 中的几个库（geopandas, pandas, scipy, numpy, matplotlib）来进行空间插值并绘制地图。

原创 Python 邮箱发送功能

【代码】Python 邮箱发送功能。

原创 按日将区域站数据插值到格点

【代码】按日将区域站数据插值到格点。

原创 python画图设置中文宋体，数字和英文新罗马

【代码】python画图设置中文宋体，数字和英文新罗马。

原创 bbox_to_anchor调整图例位置

【代码】bbox_to_anchor调整图例位置。

原创 笔记：pycharm当有多个plt.show()时候，只显示第一个plt.show()

【代码】笔记：pycharm当有多个plt.show()时候，只显示第一个plt.show()

原创 笔记：安装geopandas和cartopy时冲突问题

geopandas和cartopy存在相互依赖的底层包，安装的时候容易冲突。

原创 《现代气候统计诊断与预测技术》趋势显著性检验

【代码】趋势显著性检验。

原创 arcgis融合（dissolve）

方法：

原创 笔记：绘图进阶

【代码】笔记：绘图进阶。

原创 osgeo shapefile安装

shapefile 要用conda install pyshp安装。osgeo 要用conda install gdal 安装。关于以下两个库的安装要注意。

原创 Python 画空间插值图代码

【代码】Python 画插值图代码。

原创 python晴雨TS评分

【代码】python晴雨TS评分。

原创 Python 并行计算

在这个示例中，我们定义了一个模拟耗时操作的函数slow_function。然后，我们创建了一个包含数字的列表numbers。通过使用Parallel和delayed，我们并行地对列表中的每个数字执行了slow_function函数。最后，我们打印了结果。请注意，在使用Parallel时，你可以通过设置n_jobs参数来控制并行任务的数量。在上面的示例中，我们设置了n_jobs=-1，这意味着用最大资源进行并行计算。

原创 并行计算出现 pickle.PicklingError

这个问题是由于多进程模块multiprocessing在Windows系统上运行时的限制所导致的。在Windows系统上，multiprocessing模块使用pickle来序列化和反序列化函数，但是，对于在__main__模块中定义的函数，pickle无法正确地序列化和反序列化。要解决这个问题，你可以将my_function定义移到单独的模块中，然后在主程序中导入该模块。这样，当在Windows系统上运行主程序时，就不会再遇到pickle.PicklingError问题。

原创 Python计算日照总辐射量

【代码】Python计算日照总辐射量。

原创 不同行列的数组的合并

【代码】不同行列的数组的合并。

原创 debug报错Connection to Python debugger failed: Socket closed

出现的原因是项目中的code.py文件与调用debug相关的模块时内置的文件重名。重命名自己的code.py文件名。

原创 data[data ＞ 10000] = np.nan报错解决方法

用如下方法转为float类型即可解决。原因是该data为object类型。

原创 Python 获取最临近格点坐标索引

【代码】Python 气象格点数据插值到站点。

原创 pyinstaller PermissionError: [Errno 13] Permission denied:

关闭杀毒软件和windows实时防护。pyinstaller 安装报错。

原创 python 自定义图例可能遇到的问题

要在 Python 中创建自定义的 colormap（cmap），可以使用matplotlib.colors.LinearSegmentedColormap 类。这个类允许您创建由多个颜色组成的线性分段色带。默认情况下，LinearSegmentedColormap 类使用线性插值方法，这可能会导致颜色偏移。即使已经设置了 gamma=1.0，依旧会有颜色偏移。

原创 windows下eccodes安装方法

【代码】windows下eccodes安装方法。

原创 xgboost

【代码】xgboost。

原创 Python micaps数据转nc

【代码】Python micaps数据转nc。python读取micapspython保存成nc

原创 笔记：np查找满足条件的数组下标

代码】笔记：np查找满足条件的数组下标。

原创 python批量读取txt

【代码】python批量读取txt。

原创 Arcgis 计算两个栅格图层相关性

arcgis计算栅格图层相关性

原创 Python 将数字字符串数组转为数字数组的最快方法

将数字字符串转数组最快方法

原创 训练后的随机森林模型导出和加载

运用joblib导出和加载训练后的随机森林模型

原创 Matlab 读取nc文件报错解决方法

读取nc文件报如下错误的解决方法： 把nc文件从移动硬盘盒拷贝到电脑里。Warning: The following error was caught while executing'internal.matlab.imagesci.nc' class destructor:Error using netcdflibThe NetCDF library encountered an error during execution of 'close'function - 'Not a valid I

原创 python pycwt小波分析

安装模块pycwt这个库在官方文档中没有提供conda的安装方法，官网推荐使用pip进行安装pip install pycwt。但是，conda也可以用如下方法安装，推荐conda安装conda install -c conda-forge/label/gcc7 pycwt官网例子https://github.com/regeirk/pycwt"""In this example we will load the NINO3 sea surface temperature anomaly d

原创 python 全球降水四季空间分布图 居中经线问题

# -*- coding: utf-8 -*-"""@Features: 空间分布图 季节平均@Date：2022/4/30"""import matplotlib.pyplot as pltimport numpy as npimport netCDF4 as ncimport cartopy.crs as ccrsfrom cartopy.util import add_cyclic_pointdef sea_mean(start): temp = [] for