Python locals() 的陷阱

2023-02-22Python39

在工作中, 有时候会遇到一种情况: 动态地进行变量赋值, 不管是局部变量还是全局变量, 在我们绞尽脑汁的时候, Python已经为我们解决了这个问题.

Python的命名空间通过一种字典的形式来体现, 而具体到函数也就是 locals() 和 globals(), 分别对应着局部命名空间和全局命名空间. 于是, 我们也就能通过这些方法去实现我们"动态赋值"的需求.

例如:

很自然, 既然 globals能改变全局命名空间, 那理所当然 locals应该也能修改局部命名空间.修改函数内的局部变量.

但事实真是如此吗? 不是!

输出:

程序运行报错了!

但是在第二次 print locals()很清楚能够看到, 局部空间是已经有那些变量了, 其中也有 变量a并且值也为 1, 但是为什么到了 print a却报出 NameError异常?

再看一个例子:

输出:

上下两段代码, 区别就是, 下面的有显示赋值的代码, 虽然也是同样触发了 NameError异常, 但是局部变量 s的值被打印了出来.

这就让我们觉得很纳闷, 难道通过 locals()改变局部变量, 和直接赋值有不同? 想解决这个问题, 只能去看程序运行的真相了, 又得上大杀器 dis~

根源探讨

直接对第二段代码解析:

在上面的字节码可以看到:

从上面罗列出来的几个关键语句的字节码可以看出, 直接赋值/读取和通过 locals()赋值/读取本质是很大不同的. 那么触发 NameError异常, 是否证明通过 locals()[i] = 1 存储的值, 和真正的局部命名空间是不同的两个位置?

想要回答这个问题, 我们得先确定一个东西, 就是真正的局部命名空间如何获取? 其实这个问题, 在上面的字节码上, 已经给出了标准答案了!

真正的局部命名空间, 其实是存在 STORE_FAST 这个对应的数据结构里面. 这个是什么鬼, 这个需要源码来解答:

看到这里, 应该就能明确了, 函数内部的局部命名空间, 实际是就是帧对象的f的成员 f_localsplus, 这是一个数组, 了解函数创建的童鞋可能会比较清楚, 在 CALL_FUNCTION时, 会对这个数组进行初始化, 将形参赋值什么都会按序塞进去, 在字节码  18 61 LOAD_FAST 0 (s)中, 第四列的0, 就是将 f_localsplus第 0 个成员取出来, 也就是值 "s".

所以 STORE_FAST才是真正的将变量存入局部命名空间, 那 locals()又是什么鬼? 为什么看起来就跟真的一样?

这个就需要分析 locals, 对于这个, 字节码可能起不了作用, 直接去看内置函数如何定义的吧:

从上面 PyFrame_FastToLocals已经看出来, locals() 实际上做了下面几件事:

在这简单介绍下上面几个分别是什么鬼:

通过上面的源码, 我们已经很明确知道 locals() 看到的, 的确是函数的局部命名空间的内容, 但是它本身不能代表局部命名空间, 这就好像一个代理, 它收集了A, B, C的东西, 展示给我看, 但是我却不能简单的通过改变这个代理, 来改变A, B, C真正拥有的东西!

这也就是为什么, 当我们通过 locals()[i] = 1的方式去动态赋值时, print a却触发了 NameError异常, 而相反的, globals()确实真正的全局命名空间, 所以一般会说

locals() 只读, globals() 可读可写

Original: https://blog.51cto.com/u_15722381/5483118
Author: ch3nnn
Title: Python locals() 的陷阱

Title: python调用HEG工具批量处理MODIS数据

下面的代码主要用于使用python语言调用NASA官方的MODIS处理工具HEG进行投影坐标转换与重采样批量处理
主要参考

HEG的用户手册：https://newsroom.gsfc.nasa.gov/sdptoolkit/HEG/HEG215/EED2-TP-030_Rev01_HEG_UsersGuide_2.15.pdf
HEG批处理帮助：https://newsroom.gsfc.nasa.gov/sdptoolkit/HEG/HEG_Batch_job_Help.htm

主要的注意事项如下：

首先按照官方指南安装HEG工具，具体安装步骤可参考我的上篇博客：https://www.cnblogs.com/yhpan/p/11298595.html
根据HEG用户手册批量生成批处理参数文件，可以在HEG工具中生成一个文件，拿来自己改改用用
具体调用哪一个工具，参数文件如何写，请一定仔细阅读用户手册，东西全都在上面。一般常用的是resample.exe和swtif.exe，如果 实在无法判断可以先用HEG的GUI处理一个自己的数据，保存一个prm文件，然后根据这个文件中的参数，对照着用户手册一个一个的找，就可以了。
生成参数文件写入时 一定要注意，设定换行符为'\n',fo=open(prmfilename,'w',newline='\n')， 否则由于在windows系统下默认换行符为'\r\n',程序无法运行成功

下面是源码分享

# -*- coding: utf-8 -*-
"""
Created on Sun Feb 16 11:27:19 2020
调用HEG相关工具批处理MODIS数据，主要完成投影坐标转换与重采样
@author: pan
"""
import os

# 设置HEG相关环境变量
os.environ['MRTDATADIR']='D:/MyApps/HEG/HEG_Win/data'
os.environ['PGSHOME']='D:/MyApps/HEG/HEG_Win/TOOLKIT_MTD'
os.environ['MRTBINDIR']='D:/MyApps/HEG/HEG_Win/bin'

# 设置HEG的bin路径
hegpath = 'D:/MyApps/HEG/HEG_Win/bin'
# 指定处理模块的可执行程序文件路径，此处采用resample.exe，可以根据具体的处理问题设置
hegdo = os.path.join(hegpath, 'resample.exe')
hegdo = hegdo.replace('\\', '/') # 全路径以"/"连接

# 指定输入数据的路径
inpath = r'C:\Users\pan\Desktop\Py_ex\data\hdf'
inpath = inpath.replace('\\', '/')
# 指定输出数据的路径
outpath = r'C:\Users\pan\Desktop\Py_ex\data\hdf\out'
outpath = outpath.replace('\\', '/')
# os.chdir(inpath) #改变当前工作目录到输入数据目录

# 获取当前文件夹下的所有hdf文件
allfiles = os.listdir(inpath)
allhdffiles = []
for eachfile in allfiles:
    if os.path.splitext(eachfile)[1] =='.hdf':
        allhdffiles.append(eachfile)
print('--'*20)
print('文件数量为:', len(allhdffiles),',所有hdf文件如下')
print('  '+'\n  '.join(allhdffiles))
print('--'*20)

# prm文件设置模块，需要首先在HEG工具中生成一个参考的prm文件，示例如下
# 设置prm文件存储路径
prmpath = r"C:\Users\pan\Desktop\Py_ex\data\hdf\prm"
prmpath = prmpath.replace('\\', '/')
for eachhdf in allhdffiles:
    prm=['NUM_RUNS = 1\n',
      'BEGIN\n',
      'INPUT_FILENAME = ' + inpath+'/'+eachhdf+'\n',
      'OBJECT_NAME = MODIS_Grid_8Day_1km_LST|\n',
      'FIELD_NAME = LST_Day_1km\n',
      'BAND_NUMBER = 1\n','SPATIAL_SUBSET_UL_CORNER = ( 90.0 -180.0 )\n',
      'SPATIAL_SUBSET_LR_CORNER = ( -90.0 180 )\n',
      'RESAMPLING_TYPE = BI\n',
      'OUTPUT_PROJECTION_TYPE = ALBERS\n',
      'ELLIPSOID_CODE = WGS84\n',
      'OUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 25.0 47.0 105.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0  )\n',
      'OUTPUT_PIXEL_SIZE = 500.0\n',
      'OUTPUT_FILENAME = ' + outpath+'/'+eachhdf+'_out.tif\n',
      'OUTPUT_TYPE = GEO\n',
      'END\n']
    prmfilename=prmpath +'/'+ eachhdf+'.prm'
    prmfilename=prmfilename.replace('\\', '/')
    #这里一定要注意，设定换行符为'\n',否则由于在windows系统下默认换行符为'\r\n',则无法运行成功
    fo=open(prmfilename,'w',newline='\n')
    fo.writelines(prm)
    fo.close()

for eachhdf in allhdffiles:
    prmfilepath=prmpath +'\\'+ eachhdf + '.prm'
    prmfilepath=prmfilepath.replace('\\', '/')
    try:
        resamplefiles = '{0} -P {1}'.format(hegdo, prmfilepath)
        os.system(resamplefiles)
        print(eachhdf + ' has finished')
    except:
        # 提示错误信息
        print(eachhdf + 'was wrong')

Original: https://www.cnblogs.com/yhpan/p/12324465.html
Author: 岁时
Title: python调用HEG工具批量处理MODIS数据

Title: python之sort与sorted函数的区别

python之sort与sorted函数的区别

原创

Rickyyan2022-08-10 11:38:18博主文章分类：python ©著作权

文章标签 sort() sorted() 排序 python 文章分类 Python 后端开发

©著作权归作者所有：来自51CTO博客作者Rickyyan的原创作品，请联系作者获取转载授权，否则将追究法律责任

sort与sorted方法主要区别：

方法

应用对象

是否修改原对象

返回值

sort()

列表

是

None

sorted()

所有可迭代对象

否

返回新列表

sort(self, key, reverse):

reverse默认为False，升序

降序则，reverse=True

l = [9, 4, 2 ,5, 1, 2, 10]re = l.sort(reverse= True)print(l)print(re)#结果[10, 9, 5, 4, 2, 2, 1]None# 集合使用sort()print(set(l).sort())# 报错AttributeError: 'set' object has no attribute 'sort'

sorted(__iterable, key, reverse)

__iterable：可迭代对象，必须

key：排序字段，传入函数，可选

reverse：升序/降序，默认升序

l = [9, 4, 2 ,5, 1, 2, 10]s = set(l)print(s, '类型：',type(s))new_l = sorted(l)new_s = sorted(s)print('返回新列表：', new_l)print('返回新列表：', new_s)print('集合经过sorted()后，返回的是列表', type(new_s))print('原列表不变：', l)print('原集合不变：', s)# 执行结果{1, 2, 4, 5, 9, 10} 类型： 返回新列表： [1, 2, 2, 4, 5, 9, 10]返回新列表： [1, 2, 4, 5, 9, 10]集合经过sorted()后，返回的是列表 原列表不变： [9, 4, 2, 5, 1, 2, 10]原集合不变： {1, 2, 4, 5, 9, 10}

指定key排序

d = dict(z=3, a=9, f=6, e=1)# 以value字段排序，返回key列表new_d1 = sorted(d.keys(), key= lambda x:d[x])print(new_d1)# 以key排序，降序，返回key_value列表new_d2 = sorted(d.items(), key= lambda x:x[0], reverse= True)print(new_d2)# 执行结果['e', 'z', 'f', 'a'][('z', 3), ('f', 6), ('e', 1), ('a', 9)]