- 博客(37)
- 收藏
- 关注
RSS订阅原创 Python学习记录 20250320
定义：（1）函数嵌套（2）内层函数使用外层函数局部变量（3）外层函数的返回值是内层函数的函数名。在不改变原有代码的情况下，添加新的功能，例如：给人一把斧头，去砍树。含义：（1）函数调用自身函数的行为。（2）有一个正确的返回条件。语法：@函数名 切记：：没有小括号，是函数名。（1）标准版：（其实就是标准闭包函数的写法）例子2:计算前100个数字的和。例子3：内外层函数都带参数。例子1：打印斐波那契函数。例子1：不带参数类型的。例子2：外层函数带参数。
2025-03-20 14:23:34
123
原创 Python学习记录 20250319
下载：在cmd窗口输入pip install 文件名，但是需要进入相关的环境，因为我是虚拟环境，所以需要进入虚拟环境中安装。含义:项目结构中的文件夹或者目录，与普通文件夹相比，包含有一个__init__.py文件。一个py文件就是一个模块，里面定义了一些函数和变量，需要的时候可以导入并使用这些模块。1.文件在当前程序执行（即自己执行自己）：__name__ =="__main__"但是需要注意的是，需要在 __init__ 文件中，定义 __all__ 变量，捕获异常是为了检测到异常时，代码还能继续运行。
2025-03-19 14:57:56
378
原创 Python学习记录 20250318
lambda结合if判断，其实用到了三目运算符（代码量大、逻辑复杂时，不建议使用）匿名函数相较于普通函数我认为只是减少了代码量，关键词 lambda。关键字global可以将 局部变量声明为全局变量。去除字典、元组、列表 ，直接获取里面的数据。1.作用域：全局变量、局部变量。首先将内置函数显示出来。
2025-03-18 21:20:08
346
原创 .npy文件介绍
这种格式保留了数组的维度、数据类型（dtype）、形状（shape）等元信息，加载时无需手动解析，非常适合科学计算场景。print(data.shape) # 输出数组维度，例如 (3, 4) 表示 3 行 4 列的矩阵。data = np.load("your_file.npy") # 返回一个 NumPy 数组。加载后可直接操作数组（如 `data.shape`, `data.dtype` 查看属性）保存为 .npy：用 np.save("filename.npy", array)
2025-03-15 20:15:36
439
原创 20250309学习记录
IDM又过期了，气死我了啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊，这个为什么总是网页，各种方法都试过了。第一项，更新公众号，修改稿件。第二项，下载哨兵一号数据。
2025-03-09 22:15:20
112
原创 20250308学习记录
总结来说，这些数据产品的差异主要体现在分辨率、数据类型（如高程、坡向、坡度、坡位）以及具体的应用场景上。高分辨率数据（如30米）通常用于更精确的地形分析，而90米分辨率的数据则适用于较大范围的地形特征分析。某区高程值范围在3-××区间，说明该地区海拔较低（因为数值很小），地势起伏小，较为平坦（因为变化区间小）首先我需要的分辨率是30m，其次GDEMV3 30M用于生成全球数字高程模型，提供更高精度的地形数据，（3）筛选出来一共有15条记录，一个一个下载就可以了，也很快的。第二部分，下载DEM高程数据。
2025-03-08 22:14:01
619
原创 20250307学习记录
编写代码转换为tiff，并且将不同的土地利用类型设置为不同的颜色，一共有多少种类型呢，要查看一下文档吧，按照文档里面的要求来。根本就找不到教程，所以我在公众号上找到了别人处理过的数据，应该也是可以的，也算是有了一点收获，空间分辨率是500米。之后把名称添加上去，官网应该是英文版的，我打算把他翻译成中文的，这样更方便一些。◦ 绥化市辖：安达市、肇东市、兰西县、青冈县、明水县、望奎县（部分）第三部分，扩大研究区，发现之前选择的区域一半是农田，一半是城市。◦ 哈尔滨市辖：双城区、呼兰区、巴彦县、宾县（部分）
2025-03-07 23:02:21
633
原创 20250306学习记录
MCD12数据是利用Terra和Aqua两颗卫星遥感获得的土地覆盖类型数据，MCD12C1是其中分辨率约为0.05°的版本，适合做全球尺度的绘图。MCD12Q1为与之相似的但精度更高的数据。下载smap，生气，这个SMAP总是不能批量下载，难道还让我一张张下载码。第四项，看一下论文，把毕业论文的目录构思。投递完成，好累呀，公众号也更新了，圆满。第三项，下载SMAP和哨兵一号数据。哨兵一号，区域确定不了，没法下载。第二项，处理昨天下载好的数据。第六项，更新明天的公众号。大家早上好呀，我起床了。
2025-03-06 22:53:04
159
原创 20250305学习记录
下好了，中文文献很简单，可以在知网上直接抓取，但是外文文献需要在谷歌索引上下载，并且有的需要交钱才能下，所以我就幅值DOI号，之后在sci-hub上面搜索下载。第六部分：打算开始写毕业论文了，参照别人的硕士论文的格式写一下吧，希望在三月能写完，然后四月五月，实习一个月玩一个月，六月再回学校。目的：下载MODIS land cover数据，参考下面的博客。第一部分：接着安装Zotero中文小组插件，今天必须下载成功。参考上面的博客，安装成功，各种插件也安好了，就是不知道怎么用。之前学的，假装是今天学的。
2025-03-05 21:48:12
261
原创 20250304学习记录
知云论文翻译神器：（马上就下了一个，感觉还不错，比之前自己一句一句在有道上翻译快多了）C刊：如果要以本科、硕士、博士的身份发表，会比较困难，一般要请大牛老师当一作。第二部分，读文献的几个小工具，都看了几篇文献了，竟然到现在才知道，真是笨蛋。第三部分：突然又发现一个很多人推荐的Zotero中文小组，这我必须下载试试。第一部分，先来了解一下各种论文期刊吧，毕竟也是这把岁数了，还什么都不懂呢。小绿鲸英文文献阅读器：据说要收费，还没用。第四部分：下载MODIS数据，这不太会呀。第五部分：回宿舍睡觉，啦啦啦。
2025-03-04 21:41:43
175
原创 20250223学习记录
但是后来用Python处理的时候，直接就是对整个的.hdf5文件处理，当时没有注意这一块，所以就没有这个疑问了。之前HDFview查看.hdf5文件的时候，看到土壤湿度数据是分为AM和PM，当时我有一个这样的疑问。所以又产生了一个疑问，我用Python处理的数据到底是什么东西。（1）如果将上午和下午的数据都用上，最终的数据点会不会更多。（2）这两种数据怎么组合呢，怎么叠置呢。这么看，好像用AM数据也没有什么问题吧。今天突然看到一篇论文里面的一部分。查看自己的Python代码。所以我用的是AM数据。
2025-02-23 23:50:30
330
原创 Python基础题目练习1-5
python中re模块提供了正则表达式的功能，常用的有四个方法(match、search、findall)都可以用于匹配字符串，，正则表达式是对字符串提取的一套规则，我们把这个规则用正则里面的特定语法表达出来，去匹配满足这个规则的字符串。1.编写一个程序，找到2000年至3200年（包括在内）中所有可被7整除但不能被5整除的所有数字，得到的数字按逗号分隔，打印在一行上。列表[]：list 有序且可变的数据集合，可以存储不同类型的元素。print(k) # 打印列表形式的数字。
2025-02-21 22:26:35
994
原创 Day-7 Python学习
float 整型转换为浮点型，自动添一位小数；字符串也可以转浮点型。int 浮点型转换为整型；input 默认输入的是字符串。
2025-02-20 15:35:20
156
原创 Day-6 Python学习
先回顾一下什么是字典，在我们的日常生活中，字典是一个物品，里面有所有的汉字，我们可以按照拼音或者部首查找 我们 所需要的汉字。(1)求字典的长度 print(len(aa))（4）.values() 返回字典里面包含的所有值。学习的时候就只干学习一件事，什么也不能打扰我。（2）.keys() 返回字典中的键名。我爱学习，我爱学习，我爱学习*100。要休息一下，明天再学习，哭哭。for循环只取出键名。
2025-02-20 13:41:47
189
原创 20250217学习记录
因为之前处理的tif文件的空间分辨率都比较高，而这里的空间分辨率只有1千米，所以才会出现这种情况。不管是用Python，还是用Arcgis，裁剪出来的tiff文件并没有严格按照shp的外围进行。问题一：使用shp文件对tiff文件进行裁剪，但是出现一个问题。今天已经是正月二十了，要开始好好学习了。
2025-02-18 23:46:26
109
原创 解决matplotlib与PyCharm中的特定matplotlib backend不兼容问题
前提：Python虚拟环境第一部分：解决方案1.将已经下载的matplotlib库卸载，代码是 conda uninstall matplotlib2.下载3.8版本的matplotlib,代码是 conda install matplotlib==3.8之后再python中就可以正确运行了第二部分：介绍一下 matplotlib库的作用Matplotlib库中的常见函数众多，这些函数涵盖了从基本图形绘制到图形细节定制的各种功能。
2025-02-05 12:34:24
988
原创 Day-5 Python学习
如果元组中只有一个元素，那么类型就是该元素的类型，如果还想让其是一个元组，那么需要再唯一的元素后面加一个逗号；aa.clear() 清空字典aa中的键和值，输出变成了{}，此时可以重新增改。函数不可以被修改的时候；函数的参数和返回值（什么意思，不懂）基本格式：字典名={键1：值1，键2：值2,...}aa=[1,2,"吃饭","睡觉",'打豆豆嘟']bb=(1,2,"吃饭","睡觉",'打豆豆嘟')字典中的 键 具有唯一性，但是 值 可以重复。学习的时候就只干学习一件事，什么也不能打扰我。
2025-01-15 21:33:07
162
原创 Day-4 Python学习
（6）字符串的查找 （7）字符串的判断 （8）字符串的修改 略，感觉这种地方应该不会总用，对于我这种需要速成的人来说，暂时不需要了解这么多啊。（1）字符串拼接 + ，如果用在数字之间表示加法，如果用在字符串之间表示字符串连接。（3）in 成员运算符：检查字符串中是否包含某个子字符串。格式：列表名=【元素1，元素2，元素3，......】格式：print(st[起始位置：末位置：步长])字符串编码的本质就是字符串与二进制数字的一一编码。剩下略，感觉这个地方不是重点，以后也用不到。
2025-01-15 20:47:25
133
原创 Day-3 Python学习
range(start,stop,range)，遵循包前不包后原则，print("1-100的总和是%d"%sum)循环体（也就是条件满足时执行的代码）range用来计算循环次数，相当于一个计数器。continue 结束当前循环，进入下一循环。#利用while循环嵌套打印出9*9乘法表。4.break 和continue。#用for循环计算1-100的值。input内容默认是字符串类型。都是专门在循环中使用的关键字。需要注意——哪里需要空格。
2025-01-14 22:41:11
277
原创 SMAP数据下载
1.首先进入数据下载官网 SMAPThe NSIDC DAAC SMAP data collection includes data products derived from an L-band radiometer and high-resolution L-band radar instrument that make up the orbiting observatory of the Soil Moisture Active Passive (SMAP) satellite mission.
2025-01-13 17:22:10
1103
原创 Day-2 Python学习
r回车 将当前位置移到本行开头。\t 制表符 空四个字符，也称缩进。需要注意的是，//是向下取整，与四舍五入无关。2.比较运算符、逻辑运算符、三目运算。输入完成之后，点击回车键，输入完成。input("请输入姓名： ")\n换行符 换行。input() 输入函数。顺利输出 a和b都在笑。not 表示相反的结果。
2025-01-11 23:36:35
140
原创 Day-1 Python学习
（3）新建.py文件：在新建的项目上点击右键--新建--Python文件。end默认是换行符，如果修改end的值，那么就可以改为空格，或者是逗号。print 字符串需要加引号，print 数字可以不用加引号。（4）编写第一个程序：很简单，直接print。可以通过Debug调试，看到程序执行的顺序。Python是强类型的动态脚本语言。（2）缩进错误： print需要顶格写。（1）关闭项目：文件--关闭项目。（2）新建项目：新建项目--创建。三对单引号或者三对双引号。多看，多思考，多尝试。
2025-01-07 16:46:12
245
原创 哨兵一号数据下载
2.哨兵一号合成孔径雷达可以在四种专有模式下采集数据：SM IW EW WV。1.哨兵一号是一项成像雷达任务，可以在C波段提供连续的全天候、昼夜图像。点击Filter，上传自己的shp文件，剩下的参数按照下图配置。打开IDM--点击任务--从剪切板中批量下载 自动导入。点击Search,依次加入购物车。点击 Copy URLs。
2024-12-27 13:07:03
1086
原创 使用Python批量裁剪tiff格式数据（12.26）
（1）可能是因为，裁剪之后的tiff图像没有空间参考坐标系，但是.shp文件是有坐标系的，所以裁剪才会失败，并且需要给他们添加相同的坐标系才可以。将转换之后的tiff格式的数据导入Arcgis,这个范围值不知道有什么问题，但是在后面Python裁剪的时候，总是报错。真的是焦头烂额，又要写开题报告，又要完成导师给的任务，新手不会用Python，伤心，我要吃好吃的。（2）询问同门，怎么添加坐标系，回复是，它本身是什么坐标系，那么变换之后还是什么坐标系。（3）感觉不太对，同门应该是没有理解我的意思。
2024-12-26 22:09:57
229
原创 使用Python实现.hdf5到.tiff格式的转换（坐标系修改1206）
下载SMAP L3级土壤湿度产品，以.hdf5格式存储，使用HEGTool工具转换格式为.tiff，转换之后的.tiff文件的地理坐标系是GCS_WGS_1984（可以在HEGTool中自己选择坐标系），但是发现数据导入Arcgis之后存在变形（尤其在两极附近非常明显），所以这篇文章记录一下坐标系的选择。“”“”根据导师的指导，好像不用更改坐标系，导师说我添加的WGS 1984坐标系是多此一举，直接设置GDAL库中的GetGeoTransform()函数就可以，根据这个函数获取仿射矩阵信息即可“”“”“”
2024-12-23 21:54:19
1068
原创 用Python GDAL使用.shp文件裁剪.tif文件
裁剪失败: D:/output_1202: Free disk space available is 82176339968 bytes, whereas 27688915463940 are at least necessary. You can disable this check by defining the CHECK_DISK_FREE_SPACE configuration option to FALSE.print(f"Shapefile 文件不存在: {shape_fn}")
2024-12-04 13:58:21
410
1
原创 python实现shp对tif的裁剪1203
不管使用哪一个代码，都是一直报错Attempt to create 0x1 dataset is illegal,sizes must be larger than zero.那就不得不考虑，是原始数据的问题，猜想可能是.shp文件的问题（文件可能太小了，可能？1.解决今天下午遇到的问题，查看了很多博客，可能是由于影像太大，所以才会导致裁剪失败的。2.所以只能用Python批量处理,参考博客。首先，查看坐标系是否一致。求助，这个怎么解决呢。
2024-12-03 13:39:57
360
原创 用Arcgis提取全球影像中的一个区的数据
询问前辈，说最好用Python进行处理，因为本身我要处理的事一个月的数据，而不是一天的数据，所以还是要继续尝试一下python实现shp对tif的裁剪。这对吗，看不出来啊，但是两幅图像地理坐标系都一样的（都是GCS_WGS_1984）又遇到问题了，裁剪成功，但是不显示呀？打开属性表——按属性选择——“县名”=“×××”——根据所选要素创建图层。(之前打开的时候总是出错，今天打开，突然就对了，什么也没有改呀)经验证，不是软件的问题，应该是数据的问题。我选择的是黑龙江省哈尔滨市某区，（2）栅格数据的裁剪。
2024-12-02 18:56:06
381
原创 坐标系WGS_84_UTM_zone_47N与GCS_WGS_1984的区别
UTM 投影是一种等角横轴割圆柱投影，它将地球表面划分为一系列南北纵带（投影带），每个带再划分为纬差 6 度的四边形（虽然通常提及的是 8 度带，但考虑到中央经线两侧各 3 度的范围，实际覆盖的纬度差为 6 度）。每个投影带都有一个特定的中心经线，WGS_84_UTM_zone_47N 的中心经线大致位于北纬 47 度附近的某条经线上。投影坐标系与地理坐标系（也称为经纬度坐标系）的主要区别在于，投影坐标系通过将地球表面的点投影到一个平面上来表示位置，而地理坐标系则直接使用经度和纬度来表示地球上的位置。
2024-12-01 18:39:15
932
原创 1130修改tif影像的地理坐标系
第二种方法的实践，图像导入Arcgis看起来是失败的，但是实际上，查看坐标系是对的，已经转换成GCS_WGS_1984了，但是两种方法得到的图像显然是错误的，第一种方法的实践，不知道什么原因失败了，但是查看图层属性，地理坐标系确实已经转换成GCS_WGS_1984了，不知道为什么图像显示成这样呢。第二，在HEG里面转换，可能是因为hdf5文件转tiff的时候，一些参数设置错误，所以导致的，呢就重新设置一下参数，转换一下。坐标系：WGS_84_UTM_zone_47N。可以看出来，就是坐标系不一致所导致的。
2024-12-01 18:17:16
852
原创 基于Windows11的HEG安装教程
3.HEG的使用可以参照这篇博客，因为我是要将SMAP L3 h5格式转换为tiff，最后转换成功。4.但是现在出现一个问题，就是导入边界线的坐标，与生成的tif文件的坐标是不同的。2.完全参照博客进行下载，注意路径就可以了，最好是默认路径或者给出的参考路径。这很明显是错误的，我再检查一下在下一篇博客记录。所以将文件导入arcgis中之后出现这种情况。我下载的是JDK17，安装在默认路径。1.配置好JAVA环境。
2024-11-30 15:17:03
290
原创 下载.hdf5数据集并提取土壤湿度数据
我选取的是 某省某区（100平方公里左右）2019年七月的数据，一共有8个文件，300MB左右，直接下载大概需要三个小时（不确定是不是自己的网速太慢了^_^）对于小白来说，整个过程真的是艰难，但是大家不要放弃，坚持到底，一点点来。可视化数据：以图形方式呈现数据集，支持图像和数组数据的可视化。首先，需要在Python中安装h5py库，大家可以参考。编辑数据：在某些情况下，可以直接修改数据集中的数据。将数据导入之后，可以看到如下所示。，具体的下载过程可以参考。上下载数据，我下载的是。具体的下载过程，可以看。
2024-11-28 17:54:10
569
原创 在Python虚拟环境中安装h5py包
输入conda env list，查询Anaconda中全部环境，（重点看最后三行，*表示Anaconda当前正在使用的环境）（3）第一步中，Pycharm中的虚拟环境在gdal3.10的下面，所以我们要进入gdal3.10环境，来安装h5py。在网上查到的很多教程是直接使用pip install h5py来安装，但是我试了很久都没有成功。（4）验证是否安装成功，在Pycharm该虚拟环境下面，输入代码。（1）首先，要弄清楚，我是要在自己创建的虚拟环境中安装。（2）点击开始，找到如图所示。
2024-11-28 17:15:25
628
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人