地理信息
文章平均质量分 65
江北20190411
。
展开
-
Google Earth Engine (GEE) 提取某一位置时间序列值
前言做时间序列相关算法的同学,经常需要下载年际的数据来进行试验。而且算法研究的比较重要的步骤就是利用实测站点数据进行验证。本文讲述如何使用GEE下载指定卫星遥感/再分析数据集指定位置(实测站点位置)的时间序列值。提取指定点年序列值-imageCollection整个过程可以在imageCollection上进行操作,并且中间处理过程的返回对象都是image,这种情景非常简单直接。配合上篇博客,也可以进行QC过滤。首先定义点的经纬度坐标。然后过滤数据集合,这里以MODIS的NDVI产品MOD13A2为原创 2022-03-20 17:55:25 · 6994 阅读 · 5 评论 -
Google Earth Engine (GEE) 实现对MODIS产品批量质量控制
前言最近处在学位论文初稿完成和申博之间的空档期,所以想学点新东西。GEE火了很久了,一直没用真正使用过,还是坚持下载数据到本地,然后用Python处理。主要是怕一旦接触GEE,就会花费很多心思在上面。再者,我本身是做算法的,不仅仅是数据分析,需要使用本地数据来应用算法。现在觉得下载到本地然后处理数据有很多不确定性,不如直接在GEE上完成预处理然后导出,可以省多很多精力给算法研究本身。所以,今后主要关注于使用GEE进行数据预处理和导出的相关实现。言归正传,https://spatialthoughts.c原创 2022-03-13 22:04:48 · 2951 阅读 · 4 评论 -
Python GDAL工具使用及使用VRT格式数据处理
前言VRT及虚拟栅格文件。使用方法先找到gadl工具exe文件的路径D:\Anaconda\envs\geopandas\Library\bin对数据列表构建vrt:首先在工具路径下打开cmd,或者打开cmd输入完整工具路径构建虚拟栅格gdalbuildvrt -input_file_list E:\GeoLearn\dem_vrt_test\data_list.txt E:\GeoLearn\dem_vrt_test\data_vrt.vrt对虚拟栅格进行重投影, 并指定输出分辨率原创 2021-09-28 11:26:46 · 2637 阅读 · 0 评论 -
GWR4 软件输入数据制作
前言GWR软件是实现地理加权回归建模的专业软件。GWR软件下载地址Click here。软件下载后自带了User guidance,也可以参照这篇博客。软件要求输入的数据格式为.csv。数据中要包括投影坐标,因变量,解释变量。数据制作我这里以地表温度(LST)做因变量,NDVI做解释变量为例,使用Arcgis制作.csv输入文件。LST和NDVI的原始数据是两幅投影坐标的栅格影像。首先在Arcgis中打开影像,然后使用栅格转点工具将栅格转换为点状shp文件。得到每个像素点的shp文件。原创 2021-04-26 10:53:05 · 2292 阅读 · 8 评论 -
FY-4A 静止卫星圆盘数据几何校正
前言最近好多网友有FY-4的几何校正需求,就来更新一下。数据准备L1级的FY-4数据是HDF格式的全圆盘数据,数据集里面没有经纬度信息。每个位置对应的经纬度数据需要单独下载,如下图。这里我下载了4 km的经纬度查找表。文件下载解压后,是一个.raw格式的二进制文件。下面用python解析这个文件。import numpy as npfile = r'F:\Py_project\geoFY-4A\FullMask_Grid_4000.raw'raw_image = np.fromfile原创 2021-04-20 15:21:18 · 4601 阅读 · 9 评论 -
Python+Cartopy绘制中国地图
前言之前看了很多cartopy的教程,刚好手上有一份很久前下载的ERA5水汽数据,就练习一下吧。maskout白化的代码和colorbar分割来自气象家园的大佬。代码及数据网盘链接提取码:vtmk代码import cartopy.crs as ccrsimport matplotlib.pyplot as plt import numpy as npimport netCDF4 as ncimport cartopy.feature as cfeatureimport cartopy.io原创 2021-04-18 11:46:13 · 5402 阅读 · 6 评论 -
Python+Basemap 经纬度标签的旋转与间距调整
前言这是个非常个人的需求,原图如下,可以看到经度的标注是水平的,多子图并列的时候占空间很大,也不够美观。所以目的是将它垂直显示。垂直显示经度标签basemap绘制投影影像参看前面的博客。下面两句是分别是绘制经纬线,查看它俩的方法和属性,发现并没有控制标签角度的部分。但是basemap是基于matplotlib库的,所以可以通过操作matplotlib来实现目的。drawparallels = m.drawparallels(parallels,labels=[True,False,False,原创 2021-03-15 21:11:55 · 1708 阅读 · 0 评论 -
SURFRAD (Surface Radiation Budget) Network 数据下载
下载地址Global Monitoring Laboratory点进任意一个站点FTP,很清楚的看到数据的下载列表和READMEREADME对数据质量的说明:文件第138行给出了数据的内容和格式的描述需要注意的是每一种实测数据成列分布,质量控制代码紧跟在实测数据之后,如下所示。计算地表温度所需的F↓和F↑分别对应数据中以下内容。...原创 2020-10-12 11:03:24 · 1665 阅读 · 2 评论 -
Python+Basemap绘制已投影的影像
前言数据一副具有albers投影的地温影像,而Basemap的基本参数中要求输入影像的左下右上的经纬度坐标,所以关键在于如何将投影坐标转化为大地经纬度坐标。在坐标转换过程中,用到了pyproj库,首先要定义转换前后的坐标类型,定义WGS84坐标系可用pyproj.CRS.from_epsg(4326)进行定义,但是由于我的pyproj的数据库路径有些问题,所以首先采用gdal读取一幅WGS84坐标的影像的wkt投影信息,然后pyproj.CRS.from_wkt( )进行定义要转换成的坐标系。关于影像原原创 2020-09-27 09:54:24 · 1933 阅读 · 0 评论 -
HWSD全球土壤数据下载处理
下载地址:http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/原创 2020-04-30 10:33:57 · 27284 阅读 · 40 评论 -
ArcGIS插值可视化气象数据
前言地图底图的制作和气象数据的预处理见我上两篇博客ArcGIS精美中国地图制作(详解)python处理日值气象数据在前面我们已经将日值的地温数据处理成了年均地温,添加了站点经纬度坐标并导出为excel格式。得到如下数据:导入数据在Arcmap中 文件-添加XY数据 将年均地温和站点位置数据导入。注意:excel表格的格式要另存为97-2003工作簿形式,否则会有"arcgi...原创 2020-04-14 10:47:28 · 19944 阅读 · 8 评论 -
ArcGIS精美中国地图制作(详解)
今天看了smileliaohua老师的博客,学到很多,于是进行实验复现一下博客中的中国地图底图。下文中步骤为smileliaohua老师的步骤,但截图是我进行复现时重新取的图片,并且完善了老师博客中有些地方步骤不太详尽的地方。由于我的arcgis10.2为汉化版本,所以步骤里中英文我都进行了标注。原创 2020-04-12 14:08:20 · 28445 阅读 · 18 评论