XinemaChen的博客

极投影画轮廓线或者格网

radiance_scales，和radiance_offset这两项参数代表波段的增益和偏移量，是辐射定标的系数。比如要计算31波段的辐射亮度，读取到scales为0.00084002，offsets为1577.33972168，带入MODIS辐射定标的通用公式：Radiance=scales*(DN-offsets)，即可以得到该波段的辐射亮度。用HEG处理MOD02/MYD02，提取里面的EV_1KM_Emissive波段，band为11和12（其实就是band 31和32）。2. 辐射亮度转亮温。

【代码】gdal转投影。

这两个读取出来的分别就是经度和纬度，但是读取极投影为3413的影像时，读取出来的时投影坐标，因此在程序中多做了一步变换。基本上就是一个经纬度转影像坐标的一个操作。

在import pyproj时候报错FileNotFoundError: [Errno 2] No such file or directory: 'C:\\Users\\Asus\\AppData\\Roaming\\Python\\Python38\\site-packages\\pyproj.libs\\.load-order-pyproj-3.4.1'

只提取了MODIS数据当中的IST一个波段，输出成tif，其他什么都没做。在处理前判断这个文件是否已经被处理出tif。注意出来的tif是没有进行拉伸的，应该乘以0.01才是真正的IST的值。HEG其实可以批处理可以看我另外一篇博文，不需要写代码。但是对于300以上数量的MODIS影像非常容易自动停止，而且越来越慢。还是打算利用python每个调用。处理和直接调用HEG命名是一样的，注意处理前要先在HEG中生成一个参数化文件。

根据csv中读取的经度和纬度在极投影的地图上显示点，点的颜色按照时间显示。

MODIS批量下载

之前基于SNAP处理了Sentinel-1数据，发现地理校正之后，影像外围存在大量的nodata，使得文件大小大了一倍。对处理好的影像进行了裁剪，把白边去掉，并保留地理信息。

调查了目前能够对Sentinel-1进行处理的python库，发现除了SNAP自带的snappy和另一个完全独立于SNAP的python库外，还有一个基于SNAP本身的XML，但是进行了优化，平均处理每一景2min。目前网上代码很少，因此做一个分享。

程序包含把单个shapefile所有的要素都合并在一起，基于shapefile建立缓冲区两个功能

rasterio的矢量裁剪的方法，但是真正运用的时候报错ValueError: No valid geometry objects found for rasterize，但是实际上我的shp和raster是重叠的，发现是没有将矢量数据的投影与栅格的进行统一。另外的一个需求是需要筛选shp的属性，直接shpdata.属性名即可获取属性。

目标是希望一个dataframe按照其中的一列，datatime，按照6个小时为间隔进行均匀采样，即按照0、6、12、18点，进行时间段内的求取均值data=data.resample('6H', on='datetime').mean()data=data.reset_index(drop=False)重新将index去掉，因为后面需要访问data['datatime']报错了...

我的目的是对下面的这个ftp://sidads.colorado.edu/DATASETS/NOAA/G02202_V4/north/daily每年的文件夹进行遍历，下载1-2月的数据因此需要筛选出1-2月的文件利用了-A 这个参数，后面用-A "*198901*","*198902*" 这样就可以筛选出1989年的1-2月文件了，而且注意用逗号进行分割-r这个参数是递归下载-c这个参数是断点续传我需要对每年进行遍历，所以采用python调用cmd进行批量下载因为我的wget并.

最近从NSIDC下载数据，想批量下载每一年的1-2月数据因此对他官方的代码修改如下只需要修改main方法即可其实主要是在每个循环中对filename_filter进行修改并且每次循环都要将url_list=None，否则后面就不会重新搜索文件def main(argv=None): global short_name, version, time_start, time_end, bounding_box, \ polygon, filename_filter,

其实之前网上有很多代码主要是记录一下，首先是对日期的处理，把csv里面的日期，转换成python可以识别的日期，之后又把它转换成月份/日期我的数据是第一列为日期，其他列的列名都为年份sb.heatmap之后报错，后来发现需要将index设为日期以此记录import matplotlib.pyplot as pltimport numpy as npimport pandas as pdimport seaborn as sbfrom matplotlib import colo

首先我的程序分为两个部分，先把PIOMAS的数据提取出来pro read_piomas_monthlynv = envi(/HEADLESS) infilepath = 'H:\mission\PIOMAS\monthly\' outfilename='H:\mission\PIOMAS\monthly\hdr\' CD,infilepath thesefiles = FILE_SEARCH('heff*',count = numFiles) ;解压 FOR fidx=...

之前能够正常运行，但是突然就不行了。但是利用pip安装显示成功后将源码中的import gdal改为以下from osgeo import gdal就可以运行了

#该程序为了绘制散点图，即横坐标为时间，纵坐标有4个import matplotlib.pyplot as pltimport numpy as npimport pandas as pddata=pd.read_csv(r"H:\mission\speed\3_AIS_position_value.csv") #extent_LS.csvdata['datetime'] = pd.to_datetime(data['datetime'])c = data['piomas']c[c==-1.

最近在做一个散点图，横轴为经度，纵轴为速度，即不是纬度，因此不是那种展现空间关系的。但是数据的经度为70到180，之后-180到-170这样，因此用正常的方法出图会出现，-180和180的散点无法进行连接。于是使用python修改标签的办法。首先将横轴的数据即经度的-180~180变为0-360利用以下函数实现def changto360(data): c = data['lon'] c[c<0]=c+360 data['lon']=c return

最近由于项目需求，需要下载AIS数据，其实也就是船舶位置数据。首先通过百度发现目前一些网站下载的AIS数据都是收费的，如船达通等等。船达通是只能查找近一个月的船舶轨迹信息，下载需要支付收费。但是最后我终于找到了一个免费下载AIS数据的网站，并且下面附上了处理说明，可以通过python代码处理。...

本代码是利用读取的csv获取需要新增的ini信息，再读取ini，在原有的ini文件后面新增信息，但是不能在文件的后面插入。中间涉及到字符串与日期时间类型的转换# -*- coding: utf-8 -*-"""Created on Sun Jan 3 13:32:25 2021@author: 10+1"""import ConfigParserimport csvimport datetime#本程序可以负责利用csv表格生成对应ini及csv记录表，修改一共四个路径.

我在anaconda中运行调用sklearn没有任何问题，但是当我利用CMD命令行时调用，就出现了上面的报错。并且我利用的就是anaconda的python.exe进行调用的所以我当时找不到原因后来发现是由于安装包的来源问题，即版本兼容问题。参看此博客https://blog.csdn.net/llkk20082006/article/details/84675279https://www.lfd.uci.edu/~gohlke/pythonlibs/#在上面的地址把scikit-l

data = pd.DataFrame(columns = df1.columns)#需要先建一个为空的dataframe，注意必须有列名data=data.append(df1.iloc[i-1])data=data.append(df2.iloc[i-1])data=data.append(df3.iloc[i-1])data=data.append(df4.iloc[i-1])#然后用append函数和iloc函数即可...

利用pptx批量生成ppt，其实只是使用了模板，然后把图片变成固定格式，塞到固定位置唯一值得说的是prs.slides.add_slide(prs.slide_layouts[3])，这个序号3是指打开的ppt所对应的模板的序号（而不是ppt的页数）# -*- coding: utf-8 -*-"""Created on Wed May 6 14:07:51 2020@author: XX"""from pptx import Presentationfrom pptx.uti.

最近做的实习，带我的导师直接给我的svg格式的图片，我要把它放到ppt里，但是直接复制似乎不能显示出来但是可以用IE浏览器打开，然后再复制，但是我的图片太多了，直接写个批量会比较好用到了一个cairosvg库，安装之后，运行依然会提示缺少某个库。参看了一些博客，需要安装GTK，安装攻略可看https://www.cnblogs.com/yangjin-55/archive/2009/03/20/2786824.html里面有下载连接，按攻略安装就可以运行了还有另外一个问题是cairos.

最近修改论文已经到了最后的关口，结果被质疑分类精度评估的方法。我原来的方法是利用ArcGIS建立格网，格网中心点作为我的目视解译点，然后再从其中收集信息。...

最近在改文章的时候需要重新计算一个属性，不是那种简单的加减所以要使用到python首先尝试了直接在python里面调用arcpy其中的codeblock就是我要实现的给属性赋值的函数，expression就是调用这个函数，然后输入的值是原有的一些属性字段，因为是基于python，所以用!属性名称!。更多详见官网；https://desktop.arcgis.com/zh-cn/arcm...

当你的数据存在一定的问题时，盒图会显示不出来，并且在subplot中会直接使得之后的盒图都无法显示出来，今天就是把时间类型的处理进去便无法显示后面的，并且盒图比较适合连续性的变量分布，不适合分类型0，1这样的显示数量，这种还是留给直方图会合适很多...

思想是：在数组中选定一个数作为基准，这里选择第一个元素，遍历数组，将>=基准的放到数组右边（因为选择第一个元素作为基准如还有相等元素，必定在基准元素右侧），反之左边，之后再对基准左右两边的数组做相同操作，那么递归停止条件是当数组只剩下一个元素或者为空时。def quickSort(arr): if len(arr)<=1: return arr b...

1.剪绳子# -*- coding:utf-8 -*-class Solution: def cutRope(self, number): # write code here num=number if num==2: return 1 if num==3: return...

1.数组中重复的数字# -*- coding:utf-8 -*-class Solution: # 这里要特别注意~找到任意重复的一个值并赋值到duplication[0] # 函数返回True/False def duplicate(self, numbers, duplication): # write code here s=s...

1.和为S的连续正数序列想的方法有点复杂，看了评论区那个滑动窗口的哇瑟= =# -*- coding:utf-8 -*-class Solution: def FindContinuousSequence(self, tsum): # write code here result=[[]] sumr=0 for i...

1.丑数找2*2没乘过的丑数、3*3没乘过的丑数、5*5没乘过的丑数，三者的最小值作为新丑数# -*- coding:utf-8 -*-class Solution: def GetUglyNumber_Solution(self, index): # write code here if index<7:return index ...

1.把数组排成最小的数这个我觉得借用之前的全排列就可以是一种方法，另外我又想到如果可以让这个数组里面的数进行排列，从左数数字小的排在前，如果第一位相同，就看第二位，若一个数没有第二位就拿前一位来比较即可。然后我写了下面的代码，超时了# -*- coding:utf-8 -*-class Solution: def PrintMinNumber(self, numbers):...

1.二叉树中和为某一值的路径最初没有太看懂题目，以为从根节点到下面的子节点的某一部分就OK了，不知道要遍历到根节点= =对自己很无语啊写出了下面的代码，迷惑行为# -*- coding:utf-8 -*-# class TreeNode:# def __init__(self, x):# self.val = x# self.left ...

1.包含min函数的栈一开始自己没看懂思路，去看了讨论区意思是要两个栈，后来想明白了要用多一个栈存储较小值多一个栈的push条件：栈为空或者小于等于栈顶的元素出栈条件：栈不为空且等于栈顶的元素# -*- coding: utf-8 -*-"""Created on Sun Dec 22 14:41:47 2019@author: Shiyi Chen"""# -*...

1.顺时针打印矩阵这是我第一次的思路，说实话没懂自己为啥不能用，说是语法错误和数组越界？大致的思路是一次递归减少外面的一圈，然后把减少完的数组投入下一次递归= =，【【】】和【【1】】我都测试过了是OK的---------------------------------------------------------------我发现了这个代码的错误，在于没有考虑【【1】【2】【3...

1.树的子结构这一题我其实思路不是很明确，是看了讨论区的# -*- coding:utf-8 -*-# class TreeNode:# def __init__(self, x):# self.val = x# self.left = None# self.right = Noneclass Solution: ...

9.二进制中1的个数class Solution: def NumberOf1(self, n): # write code here num=0 if n<0: num=str(bin(2**32+n)) else: num=str(bin(n)) ...