个人学习笔记

解决方法是：把区域字段换成FID居然就可以了。如果还不行的话 ，要记得输出表的位置也改成系统默认的，路径都是英文。问题描述：当我按照下图操作时，报了999999的错误。

参考：高分辨率遥感影像在哪免费下载？- Endless的回答 - 知乎 https://www.zhihu.com/question/373272521/answer/2786776699。第三步：提交，出现此弹窗，等待邮件回复。（收到邮件回来更新）第三步一开始我没有弹出该窗口，在浏览器做以下修改后弹窗成功。

一、tif中只有红（band1）绿（band2）蓝（band3）近红外（band4）四个波段

参考：https://www.cnblogs.com/youzi-xuchongyou/p/9241891.html错误：因为路径中有中文mymxd = "E:/myPy/in_data/china400/简称.mxd"解决方法：路径前加中文mymxd = u"E:/myPy/in_data/china400/简称.mxd"

下载地址：FROM-GLC10 - (2017 V0.1.3 with minor improvements) File List下载界面：前面的0是纬度，100是经度，间距都是2举个例子，要下载北京市土地利用类型，先看北京市范围：于是下载了【40_116】【40_114】【38_116】【38_114】这四幅，把它们在arcmap里镶嵌裁剪，设置一下符号系统：每个值代表的土地类型如下：http://data.ess.tsinghua.edu.cn/...

USGS：https://earthexplorer.usgs.gov/

参考：https://www.pythonheidong.com/blog/article/486038/253098ead152c6bec4a7/问题描述：我的研究区只有左图中带颜色的部分，但是最开始用matplotlib着色的时候，那些nodata的全被当0着成了红色。把ndvi.tif（代码中的gray）的ndarray打印出来，发现nodata值全显示为-3.40282347e+38，但是gray[gray==-3.40282347e+38]=np.nan无效。解决方法：gray[(gr

压缩包：影像：landsat08_L1TP_条带号行编号_获取日期_处理日期_01_T1_ANG(投影信息)/MTL(影像元数据)/BQA(影像质量信息)/B1-B11(波段)

二、辐射定标右键影像查看元数据，辐射校正需要的绝对定标系数gains和offsets是ENVI系统自带的，无法自动更新，因此需要手动修改。从中国资源卫星应用中心获取2020年最新定标系数：http://www.cresda.com/CN/Downloads/dbcs/index.shtmlhttp://www.cresda.com/CN/Downloads/dbcs/images/2020/11/11/240A9946ECED64C3CB4D56B8A8764F35.pdf我的影像是GF1WFV3，

1、全球森林变化https://www.globalforestwatch.org/除了森林增减等，还有气候、火情、生物等数据2、生命地图https://mol.org3、全球捕鱼观察https://globalfishingwatch.org/map/

参考：利用百度地图坐标拾取系统获取经纬度，并将其坐标转换保存成SHP文件txt文件长这样（最后一行不是第一行也可以）：代码：import osfrom osgeo import osr, ogrdef coordinateTrans(txtPath): folder = str(os.path.split(txtPath)[:-1][0]) driver = ogr.GetDriverByName("ESRI Shapefile") out_shp_path = os.

1、天天向上的力量def sum(n): a=1 for i in range(365): if i%7 in [0,6]: a=a*(1-0.01) else: a=a*(1+n) return an=0.01while sum(n)<37.78: n=n+0.001print("工作日的努力参数是: {:.3f}".format(n))...

参考：https://www.runoob.com/note/27030import ospath = "E:/网页下载/change/"fileList = os.listdir(path)for i in range(len(fileList)): oldname = path + os.sep + fileList[i] # os.sep添加系统分隔符 newname = path + os.sep + fileList[i][:8] + fileList[i][-4:]

创建虚拟环境：conda create -n myenv python=3.6删除虚拟环境：conda remove -n myenv --allmyenv是想创建/删除的虚拟环境名称

from flask import Flaskapp=Flask(__name__)@app.route("/")def hello_world(): return 'Hello World!'if __name__=='__main__': app.run()

小姜说等她生日的时候给她写个生日快乐歌MIDI To RTTL:http://midi.mathewvp.com/midi2RTTL.php单音轨下载网址参考：https://www.bilibili.com/video/BV1LK4y1L7eR?share_source=copy_web%matplotlib inlinefrom pylab import *from IPython.display import Audiofrom scipy.signal import square..

import fileinputfile = "./in.txt"for line in fileinput.input(file): print(line,end="")fileinput.close()path = "./in.txt"file = open(path, "r", encoding="utf-8")file_text = file.readlines()for line in file_text: print(line, end="")file.clos

1、win+R->cmd2、输入"cd whl存放文件夹路径"，例如cd c:\python37_64，回车3、pip install ______.whl

代码：import osdef ListFilesToTxt(dir, file, wildcard, recursion): exts = wildcard.split(" ") files = os.listdir(dir) for name in files: fullname = os.path.join(dir, name) if (os.path.isdir(fullname) & recursion):

import pandas as pddef readwrite(input_file,output_file): data=pd.read_csv(input_file,sep=',',header=None,engine='python') data.to_csv(output_file,sep=',',columns={1},header=None,index=False)input_file=r"产量.csv"output_file=r"LinearRegression.c

1、蟒蛇import turtleturtle.setup(650, 350, 200, 200)turtle.penup()turtle.fd(-250)turtle.pendown()turtle.pensize(25)turtle.pencolor("purple")turtle.seth(-40)for i in range(4): turtle.circle(40, 80) turtle.circle(-40, 80)turtle.circle(40, 80/2

参考：https://blog.csdn.net/feng98ren/article/details/89968699打包时输入pyinstaller -F xxx.py --hidden-import sklearn.neighbors.typedefs成功解决

二分支结构紧凑形式：guess=eval(input())print("猜{}了".format("对" if guess==99 else "错"))

gray = image(f, 1) img = np.array(gray) for i in img: for j in range(len(i)): if i[j]<=0: i[j]=0 elif i[j]<=0.2: i[j]=0.2 elif i[j]<=0.4: i[j]...

import pandas as pdimport numpy as npfrom sklearn import linear_modeldf=pd.read_csv(r"LinearRegression.csv",engine='python')col=df.shape[1]X=df.iloc[:,2:col]y=df.iloc[:,1:2]X=np.matrix(X.values)y=np.matrix(y.values)regr=linear_model.LinearRegre

%s 字符串%d 十进制整数%i 十进制整数%f 浮点数（默认小数点后6位）默认右对齐，加“-"变左对齐例如：x,y,z="john",20,1.85print("5Name:%8s Age:%4d Height:%5.2f"%(x,y,z))# 5Name: john Age: 20 Height: 1.85print("6Name:{:8s} Age:{:4d} Height:{:5.2f}".format(x, y, z))# 6Name:john Age: 20

代码：txtFile=r"E:\▲▲▲\研一下\Jane Eyre.txt"text=open(txtFile).read()text=text.lower()for m in [';',':',',','!','.','?','-','"',"'",'(',')']: text=text.replace(m," ")#.split()默认分割空格words=text.split()wordStat={}for word in words: #wordStat.get(word,0

参考：https://blog.csdn.net/jacke121/article/details/113059922首先进入我本来创建的python环境（名字叫lianxi）：activate lianxi然后这样就可以了。。conda install python=3.10查看python版本，成功后就变成3.10了...

第一次在自己电脑上搞python，做个记录先下载anaconda，https://repo.anaconda.com/archive/选了这个版本的，没有选最新的：装在这里了：两个选项都选了：pycharm下载的2019专业版，勾选了这三个：

import gdalshp = r".\shapefile\shapefile.shp" # 圈选范围的路径inputImage = r".\image\image.tif" # 遥感影像的路径dataset = gdal.Open(inputImage) # 打开遥感影像outputImage = r".\image\image_mask.tif" # 按照圈选范围提取出的影像所存放的路径gdal.Warp(outputImage, dataset, cutlineDSName=

第1章-1 从键盘输入两个数，求它们的和并输出a=input()b=input()c=eval(a)d=eval(b)print(c+d)

改下目录名称就可# -*- coding: utf-8 -*-# @Time: 2022/4/26 18:40import arcpy# 解决中文乱码import sysreload(sys)sys.setdefaultencoding('utf8')# 要投影的shp所在的目录arcpy.env.workspace = "E:/▲▲▲/2022-项目/220422-肃州玉米/肃州/"arcpy.env.overwriteOutput = True# 投影后的文件要存放的目录ou

1、圆面积的计算交互式：文件式：2、同切圆的绘制交互式：文件式：3、绘制五角星

import ospath = "D:/anaconda3/envs/lianxi/study/220415"file = os.listdir(path)for i in file: print(i)

废话一篇

hhhhh

笔记