zsmj的博客

原创 全球水库水位数据

最近处理数据时，看到一个比较好的全球水库水位数据DAHITI，主要是用测高卫星测量的。

原创 MERRA-2 土壤湿度 2024-10-22

这两天想用下MERRA-2的土壤水数据，发现资料很少，而且MERRA-2的变量名跟平常见到的不太一样，做个笔记，写下自己的一点收获，也希望以后用到的人少花点时间。特别需要注意的是，变量里面的units我们可以看到是“1”，这指的是无量纲的意思，我们可以将它视为百分比（下表示为WetPercent）。我所涉及到的有着三个变量："GWETPROF"，"GWETROOT"，"GWETTOP"如果你计算第一个，是需要的，因为深度不固定，里面长这样。这里比较推荐 用（2）下载，（1）是我意外翻到的网页。

原创 运行PCR-GLOBWB V2模型_2

最近出于个人需要，跑了PCR-globwb的陆地水和地下水数据，如果有需要可以私信（可能忙时顾不上回）

原创 matlab engine for python与python 的版本对应

matlab engine for python与python 的版本对应

原创 论文发表时，如何在保证图片清晰的同时，改变图片尺寸

保存后，我们使用其他软件打开这个pdf,就可以看到100%时，图片已经很小了，期刊基本是A4纸，此时是可以容得下的。本文用到的文件如下，可以看到这个图片在12.5%大小的时候已经占了大半个屏幕了，可以想象它的尺寸有多大。2.将刚刚存储的pdf拉进ps，此时就可以改尺寸大小了，此时 按住alt ctrl +I，A4纸的大小是20左右，我们可以保存为17cm,这样的话就保存好了。上面时17cm，下面是我设置的8cm的，可以看到图片尺寸已经变化了。我们将8cm的进行放大，可以看到，图片的清晰度是没有变化的。

原创 python STL时序分解错误

把pd.DataFrame改成pd.Series.python进行时序分解时出现错误。最终发现是数据格式的问题。

原创 GRACE/FO数据下载网址

①此页面会慢慢更新；②ftp网址推荐使用8Uftp下载。

原创 运行PCR-GLOBWB V2模型

根据下载好数据后就可以开始运行程序（github网址）了。我运行的python版本是python3.8+Anaconda。个人感觉此时python3.8运行的比较平稳。我在考虑运行模型前看了这篇总结，，大家可以多看看关于PCR的资料，心里就大致明白整个流程了。运行模型下来，我感觉最难的可能不是运行模型（有可能是我的配置和运行过程比较顺利），而是根据自己的应用对数据和参数进行修改。

原创 python selenium批量获取PCR-GLOBWB的输入数据

最近要跑一个PCR-GLOBWB的30min的模型数据，但是里面的数据很多，下载量也比较大，我就根据PCR-GLOBWB github上面给的网址，这个网址是散装文件（有一个zenodo网址给的是一个45G大文件，实在是不得行），于是就写了一个批量下载的脚本，这样就不用一个个单击下载了。

原创 如何创建自己免费同步网盘？（类似坚果云，但是薅羊毛）

坚果云非常好用，但是收费比较高哈哈哈，这里是一个类似坚果云的同步盘设置，非常简单就可以设置一个共享网盘文件夹，用起来也很顺手，用的是tera cloud+RaiDrive

原创 蹭下ChatGPT的热度哈哈

蹭下ChatGPT的热度哈哈

原创 插个眼：python colorbar使用科学计数法

No

原创 MODIS MOD16A2 download using python----0

modis下载

原创 处理MOD16A2 using python and arcgis：from hdf to tif----1

from hdf to tif

原创 python selenium文件下载

python selenium爬虫大批量下载文件,但浏览器默认打开文件

原创 GRACE信号中的GIA用什么文件去除？

GIA（glacial isostatic adjustment）的产品找寻路径，也算是GIA的的文件介绍吧

原创 如何用python做出一个对应shapfile文件的格网

这篇文章的目的是告诉大家怎么用python批量做一个shapefile内的规则格网。其实arcgis已经能够完成这些问题，但为了更省事，我们可以用一些程序来解决他们。中间涉及到一些python包的安装，自行安装吧。（一）不多说，直接上代码吧！！import osimport numpy as npimport mathimport pandas as pdfrom tqdm import tqdmfrom shapely.geometry import Polygon, Point,

原创 GRACE重力卫星水文应用基本概念之trend

原创 Python绘图-怎么获取某个值在色带中对应的颜色

这两天遇到一个问题，感觉还挺有意思。在使用matlplotlib的pyplot进行绘图时，需要根据我的一个具体的值，获取在色带中所对应的颜色。相关的资料很少，并且有的可能不顶用。这里简单的记录一种实现该功能的方式。import matplotlibnorm = matplotlib.colors.Normalize(vmin=-40, vmax=40)#色带的范围是-40至40color0=matplotlib.cm.coolwarm_r(norm(10),bytes=True)这里我以

原创 GRACE重力卫星水文应用基本概念之时变重力场

简单说明一下地球时变重力场，由于个人主要是做应用，可能有些理论的地方写的不是很恰当，欢迎批评指正

原创 GRACE重力卫星水文应用基本概念之Anomaly

What?Anomaly在GRACE水文应用中通常指的是 某月的水储量相对于某距平时间段（如GRACE mascon数据通常的距平时间段位2004-2009，因为这段时间GRACE比较完整）的差值。其表达的形式通常为TWSA（陆地水储量异常）、SMA（土壤水储量异常） 、GWSA（地下水储量异常）等一些理解Anomaly是GRACE在水文领域中经常用到的词语。GRACE卫星获取的其实是某个地区（或者某个格网）的重力信号，其中绝大多数是由水引起的，但也有其他信号比如矿产开采，泥沙运移等。下文为了便

原创 pandas快速绘图

数据长这样：import matplotlib.ticker as tickerax=ts_pd.loc[['wghm', 'Clsm','Mos', 'Noah', 'Vic'],:].T.plot(title='Model TWSA(Yangtze)',figsize=(9,3),rot=45,grid=True,xticks=range(ts_pd.shape[1]-12))ax.xaxis.set_major_locator(ticker.MultipleLocator(12))..

转载 arcgis操作记录，拆分面图层

https://blog.csdn.net/Daisy_jf/article/details/103273840

转载 arcgis将一个面拆分成多个面polygon

https://blog.csdn.net/ucs426/article/details/103183770

原创 2021-10-07

tensorflow训练时出现模型读取错误的解决在model.fit函数运行时，出现的错误是UnicodeDecodeError: ‘utf-8’ codec can’t decode byte 0xbe in position 0: invalid start byte，我的关键代码是cp_callback = tf.keras.callbacks.ModelCheckpoint(filepath=checkpoint_save_path,save_weights_only=True,save

原创 arcgis融合多个polygon,形成一个单独的文件

融合前的shapefile图层，其属性表为：其中obj_status_n列全是0，我们待会需要选中这一列，必须要选一个全部是同一个值的列。接着找到这个工具点开，之后选择为：点击ok,开始运行，结果如下：...

原创 python画图示例：上下为线，平均值是空心的箱体图

先看图再上代码import numpy as npimport pandas as pdimport matplotlib.pyplot as plt df = pd.DataFrame(np.random.rand(10,5),columns=['A','B','C','D','E'])f = df.boxplot(sym = 'o', whis=1.5, meanline = False,showmeans = Tru

原创 python科研绘图：左上角RMSE，右下角R2,并画出双色带

python科研绘图：左上角RMSE，右下角R2,并画出双色带 1.结果图2.直接上代码1.结果图2.直接上代码import numpy as npimport matplotlib.pyplot as pltimport random#自己做的数据，可自行调整label1=[‘a’,‘b’,‘c’,‘d’,‘e’,‘f’,‘g’,‘h’]arr1=np.zeros((len(label1),len(label1)))*np.nanarr2=np.zeros((len(label1),l

原创 python判断某一年的第n天的年月日

import datetimedef get_date(year=2017,days=123): the_date = datetime.datetime(year,1,1) result_date = the_date + datetime.timedelta(days=days) d = result_date.strftime('%Y-%m-%d') return dget_date(year=2017,days=123)

原创 nc文件的经纬度判定

nc文件的经纬度判定与验证关于nc文件的读取，相信很多小伙伴会遇到经纬度的判定问题（别问我怎么知道的…）。在这里，我将根据搜集到一些nc文件进行分析，对它们的经纬度表示方式加以判读，并导入arcgis中进行验证。希望对大家的读取工作有所帮助。1.GLDAS_NOAH025_M.A201705.021.nc4这个是水文模型GLDAS noah的数据，格网是0.25度。图1为该nc文件 中的变量"LWdown_f_tavg"。...