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

首先是我的数据下载方面，从官网上直接下载没有使用代码，因此所有下载的数据都放到一个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

最近在改文章的时候需要重新计算一个属性，不是那种简单的加减所以要使用到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...