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

帮同学写的一个程序，目的是：把很多个NC数据合并在一起，并且每个.nc数据的变量名称是相同，我平常读.nc数据比较多，较少涉及写入。我先使用了IDL但是由于原数据的属性信息也很多，因此最后放弃。而MATLAB有基于现有的数据创建一个模板，对数据进行生成。代码如下：1.先对需要的.NC数据进行提取2.提取所需的数据属性等3.对数据进行赋值ticclose all;clear...

本代码是利用读取的csv获取需要新增的ini信息，再读取ini，在原有的ini文件后面新增信息，但是不能在文件的后面插入。中间涉及到字符串与日期时间类型的转换# -*- coding: utf-8 -*-"""Created on Sun Jan3 13:32:25 2021@author: 10+1"""import ConfigParserimport csvimport datetime#本

首先是我的数据下载方面，从官网上直接下载没有使用代码，因此所有下载的数据都放到一个NC文件中了（当然2m温度和气压是分开的）我主要下载的区域是北纬30-90，每日只下载四个时间，0 6 12 18，也只下载6个月的，1-3、10-12月nc文件可以尝试使用panoly打开看看，每个月和每个时间是按顺序排好的我利用以下的代码来处理2m温度和海平面气压数据值得注意的点是，在读取数据时，还要读取校正系数

我想要利用matlab读取geotif后，输入经纬度得到该图像的坐标以及DN值。首先在网上能找到的是利用图像的坐标值得到经纬度坐标，就是这个博客。那么首先是[A,refmat,bbox] = geotiffread('H:\georef.tif');info = geotiffinfo('H:\georef.tif');A获取的是图像的矩阵本身，info好像获取的是各种图像的信...

首先需要打开影像，然后在上面建立ROI利用Basic Tool-Region of Intertrest-Create Class Image from ROI可以将基于ROI的mask建立，具有和影像相同的大小和投影信息。生成的mask，maks值为1，非mask值为0

Sentinel 1文件名命名规则如下图：官网的详细说明：https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1-sar/naming-conventions以下为例子：S1A_EW_GRDM_1SDH_20180112T082556_20180112T082700_020119_0224E3_7F9BS1A...

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

最近需要下载同一景的sentinel数据，最后搜了一下决定使用python使用的python包如该链接https://pypi.org/project/sentinelsat/里面也有详细说明文档因此这里只粘贴代码了from sentinelsat.sentinel import SentinelAPI, read_geojson, geojson_to_wktfrom dat...