GeoSuper的博客

原创 哨兵3号—sentinel-3

目录1、数据介绍1、数据介绍ref：Clevers J G P W, Gitelson A A. Remote estimation of crop and grass chlorophyll and nitrogen content using red-edge bands on Sentinel-2 and-3[J]. International Journal of Applied Earth Observation and Geoinformation, 2013, 23: 344-351.

转载 MODIS数据知识积累

@[TOPIC] (MODIS数据积累)MODIS数据使用的投影sinusoidal 投影

原创 利用USGS EarthExplorer 批量下载sentinel-2、landsat 8 Level1、2等影像

打开USGSEarthExplorer网站https://earthexplorer.usgs.gov/（要注册登录） 选择自己的研究区 （单击鼠标左键或者放大到自己的研究区use map） ...

原创 哨兵2号（Sentinel-2）介绍、下载、预处理及批处理

哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。6月29日，在轨运行4天的哨兵-2A卫星，传回了第一景数据，幅宽290km，卫星第一次扫描的范围是从瑞典开始，经过中欧和地中.........

原创 Sentinel-2批量大气校正_基于Sen2Cor

将sentinel-2数据放在一个文件夹下1、2、for /D %s in (E:\1\S2A_MSIL1C*) do L2A_process --resolution=10 --refresh %s（分辨率自己根据需要修改）实践证明：for /D %s in (E:\1\S2A_MSIL1C*) do L2A_process %s会比较好些，不容易出错。如果文件名是这种的‘S.........

原创 利用MRT对MODIS数据进行批量重投影+批量波段合成

写在前面：官方渠道已经下载不到MRT了，为什么呢？因为根本不需要使用MRT,可以直接下载处理好的MODIS：详见：https://mp.weixin.qq.com/s/HkKdap1rcpGidftXsuUbRQMRT下载：https://pan.baidu.com/s/1_pDwv2lghAzTpaLGv7IcOw提取码：reygMRT安装：1、解压（其中m.........

原创 Landsat 介绍，大气校正产品申请下载

(Wulder M A, Loveland T R, Roy D P, et al. Current status of Landsat program, science, and applications[J]. Remote Sensing of Environment, 2019: 127-147.) 波段介绍 OLI陆地成像仪 波段名称 ...

翻译 辐射定标、大气校正

辐射校正指在光学遥感数据获取过程中，产生的一切与辐射有关的误差的校正（包括辐射定标和大气校正）。三者关系如图：大气校正的准备过程为辐射定标表观反射率：表观反射率就是指大气层顶的反射率,辐射定标的结果之一，大气层顶表观反射率，简称表观反射率，又称视反射率。英文表示为：apparent reflectance辐射定标：将记录的原始DN值转换为大气外层表面反射率（或称为辐射亮度值）。...

原创 全球制霸生成器，我只能得5分了

全球国家不一样的饿可视化效果

原创 rgee学习笔记

rgee学习笔记01_rgee_hello world02_rgee_simple mapdisplay03_rgee_finding images04_rgee_band math01_rgee_hello worldlibrary(rgee)ee_Initialize()print("Hello World")# [1] "Hello World"#python styleee$String("Hello World from Earth Engine!")$getInfo()#

原创 36_geemap学习笔记 | 质量镶嵌（最大NDVI）

greenest = withNDVI.qualityMosaic('NDVI')

转载 34_geemap学习笔记 | 像素值的交互式提取

转载 32_geemap学习笔记 | 监督分类

断了学习geemap很久了，今天继续今天迎来了吴秋生老师geemap的更新，所以咱们先更新一波

原创 Python&R | nc批量转tif

nc数据：NetCDF(network Common Data Form)网络通用数据格式。NetCDF 文件中的数据以数组形式存储。例如：某个位置处随时间变化的温度以一维数组的形式存储。某个区域内在指定时间的温度以二维数组的形式存储。三维 (3D) 数据（如某个区域内随时间变化的温度）四维 (4D) 数据（如某个区域内随时间和高度变化的温度）以一系列二维数组的形式存储。本次使用的样例数据为2018-2020年的月均温度netCDF数据下面主要介绍如何使用Python和R语言批量将上述nc格

原创 遥感影像预处理之~自动配准

方法一、ArcGIS Pro利用ArcGIS Pro的自动配准查看控制点表可以看到自动打了49个点，还是均匀分布于影像的局部结果如下图所示，然后另存为新的栅格即可。总结：首先进行自动配准，如果不好的话，可再次进行自动配准，然后打开控制点表格，删除误差比较大的点，边删边看下自动配准的怎么样，基本都能达到满意的结果。方法二、Python(arosics包)arosics作者githubarosics使用手册根据使用手册安装arosics及相关依赖库纸上得来终觉浅，实践一下：impor

转载 31_geemap学习笔记 | 非监督分类

原文：31 Unsupervised Classificationhttps://developers.google.com/earth-engine/clusteringThe ee.Clusterer package handles unsupervised classification (or clustering) in Earth Engine. These algorithms are currently based on the algorithms with the same name

转载 30_geemap学习笔记 | 影像属性统计

30 image pops statsimport eeimport geemapgeemap.set_proxy(port=10809)Map = geemap.Map()Mapcentroid = ee.Geometry.Point(117.1439, 30.7538)landsat = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') \ .filterBounds(centroid)\

转载 26_geemap学习笔记 | 修改图层透明度

import eeimport geemapgeemap.set_proxy(port=10809)Map = geemap.Map(center=(40, -100), zoom=4)# Add Earth Engine datasetdem = ee.Image('USGS/SRTMGL1_003')landcover = ee.Image("ESA/GLOBCOVER_L4_200901_200912_V2_3").select('landcover')landsat7 = ee.I

转载 24_geemap学习笔记 | 发布地图

import eeimport geemapimport geemap.foliumap as gfmapgeemap.set_proxy(port=10809)Map = geemap.Map(center=[38,105], zoom = 4)ele = ee.Image('USGS/GMTED2010')Map.setCenter(117.5,32,6)Map.addLayer(ele, {}, 'Elevation')shade = ee.Terrain.hillshade(el

转载 23_geemap学习笔记 | 导入资源

23 import assetsimport osimport eeimport geemapgeemap.set_proxy(10809)Map = geemap.Map()Mapout_dir = os.path.join(os.path.expanduser('~'), 'Downloads')if not os.path.exists(out_dir): os.makedirs(out_dir)out_csv = os.path.join(out_dir, 'ee_

转载 22_geemap学习笔记 | 导入脚本

import eeimport geemapgeemap.set_proxy(port= 10809)Map = geemap.Map()Mapgeemap.ee_search()

转载 21_geemap学习笔记 | 导出地图为html/png/jpg

import eeimport geemapgeemap.set_proxy(port=10809)geemap.show_youtube('h0pz3S6Tvx0')

转载 20_geemap学习笔记 | 时序检查

20_timeseries inspectorNAIP: National Agriculture Imagery ProgramThe National Agriculture Imagery Program (NAIP) acquires aerial imagery during the agricultural growing seasons in the continental U.S.NAIP projects are contracted each year based upon ava

转载 19_geemap学习笔记 | 搜索地点和数据集

19 search places and datasetsimport eeimport geemapgeemap.set_proxy(port=10809)geemap.show_youtube('lwtgzrHrXj8')Map = geemap.Map()MapMap.search_locationsMap.search_loc_geomee.Geometry({# "functionInvocationValue": {# "functionName"

转载 18_geemap学习笔记 | 创建Landsat延时影像

import geemapgeemap.set_proxy(port=10809)geemap.show_youtube('OwjSJnGWKJs')Map = geemap.Map(center=(31,118),zoom =6)QS_shp = "F:/work/20211013/qian.shp"QS_shp = geemap.shp_to_ee(QS_shp)Map.addLayer(QS_shp, {}, '潜山市shp')Map

转载 17_geemap学习笔记 | 在gif上添加颜色条

17 add colorbar to gif下载一个GIFimport osimport geemapgeemap.set_proxy(port=10809)from geemap import *url = 'https://i.imgur.com/MSde1om.gif'out_dir = os.path.join(os.path.expanduser('~'),'Downloads')if not os.path.exists(out_dir): os.makedirs(o

转载 16_geemap学习笔记 | 创建文本动画

16 add animated text# !pip install geemapimport osimport geemapgeemap.show_youtube('fDnDVuM_Ke4')可以下载一个gif图，也可以使用LT-GEE Time Series Animator自己创作一个。in_gif = os.path.abspath('../data/animation.gif')out_dir = os.path.join(os.path.expanduser('~'),'Dow

转载 15_geemap学习笔记 | js转py

15 convert js to pyimport eeimport geemapgeemap.set_proxy(10809)geemap.show_youtube('nAzZjKKd4w0')实例1# 实例1js_snippet = """// Load an image.var image = ee.Image('LANDSAT/LC08/C01/T1_TOA/LC08_044034_20140318');// Define the visualization parame

转载 14_geemap学习笔记 | 图例

14_legendsimport eeimport geemapgeemap.set_proxy(port=10809)geemap.show_youtube('NwnW_qOkNRw')Available builtin legendslegends = geemap.builtin_legendsfor legend in legends: print(legend)# NLCD# NWI# MODIS/051/MCD12Q1# GLOBCOVER# JAXA/P

转载 13_geemap学习笔记 | 按组分区统计

原文：zonal statistics by groupimport eeimport geemapgeemap.set_proxy(10809)geemap.show_youtube('cORcGGH03gg')Map = geemap.Map(center=(38,105))Mapimport requestsimport jsonurl = 'https://geo.datav.aliyun.com/areas_v3/bound/100000_full.json'China

转载 12_geemap学习笔记 | 分区统计

import eeimport geemapimport osgeemap.set_proxy(10809)geemap.show_youtube('ou-Xm3CLitM')Map = geemap.Map(center=(38,105))Map# Add Earth Engine datasetdem = ee.Image('USGS/SRTMGL1_003')#Set visualization parametersdem_vis = { 'min': 0,

转载 11_geemap学习笔记 | 导出影像

导出图像原文入口Download an ee.ImageDownload an ee.ImageCollectionexport pixels as a Numpy array原文入口import eeimport osimport geemapgeemap.set_proxy(port=10809)# geemap.show_youtube('_6JOA-iiEGU')ee.Initialize()Map = geemap.Map()MapDownload an ee.Imag

转载 10_geemap学习笔记 | 矢量数据

10_shapefilesimport geemapgeemap.set_proxy(port=10809)geemap.show_youtube('OlNlqfj4uHo')Map = geemap.Map(center = (38, 105))Map

原创 RGB颜色与十六进制颜色码

英文代码形像颜色HEX格式RGB格式LightPink浅粉色#FFB6C1255,182,193Pink粉红#FFC0CB255,192,203Crimson猩红#DC143C220,20,60LavenderBlush偏红的淡紫色#FFF0F5255,240,245PaleVioletRed苍白的紫罗兰红色#DB7093219,112,147HotPink热情的粉红#FF69B4255,105,180Dee...

转载 09_geemap学习笔记 | 绘图

09_绘图import eeimport geemap# geemap.show_youtube('PDab8mkAFL0')Landsat7 = ee.Image('LE7_TOA_5YEAR/1999_2003') \ .select([0, 1, 2, 3, 4, 6])Landsat_vis = { 'bands': ['B4', 'B3', 'B2'], 'gamma': 1.4}Map.addLayer(La

转载 08_geemap学习笔记 | ee_js_to_ipynb

自动转换谷歌地球引擎js代码至python代码import geemapgeemap.set_proxy(port=10809)geemap.show_youtube('RpIaalFk4H8')Convert Earth Engine JavaScripts to Python scriptsimport osfrom geemap.conversion import *# create a temporary working directorywork_dir = os.path.

转载 07_geemap学习笔记 | geojson

07_geojsonimport geemapimport jsonimport osimport requestsfrom geemap import geojson_to_ee, ee_to_geojsonfrom ipyleaflet import GeoJSONgeemap.set_proxy(port=10809)geemap.show_youtube('DbK_SRgrCH')Map = geemap.Map(center = (38, 105))Mapfile_

转载 06_geemap学习笔记 | 标记聚类

06_marker_clusterimport geemapimport jsonimport osimport requestsfrom geemap import geojson_to_ee, ee_to_geojsonfrom ipyleaflet import GeoJSON, Marker, MarkerClustergeemap.set_proxy(port=10809)geemap.show_youtube('4HycJPrwpuo')Map = geemap.Map(

转载 05_geemap学习笔记 | 绘图工具

05_drawing_toolsshp1.filterBounds(roi)(求roi在shp1内所涉及的矢量部分(shp内的n块))裁剪：image1.clip(roi)(基于roi裁剪image1)import eeimport geemapgeemap.set_proxy(port=10809)# geemap.show_youtube('N7rK2aV1R4c')Map = geemap.Map(center = (38 ,105))Map叠加地图# add earth e

转载 04_geemap学习笔记 | 创建拆分（双面板）地图

04_split_panel_mapsimport geemapgeemap.set_proxy(10809)# geemap.show_youtube('9EUTX8j-YVM')Map = geemap.Map()Map.split_map()MapMap = geemap.Map(center = (38, 105))Map.split_map(left_layer='HYBRID', right_layer = 'ROADMAP')Map查看geemap包里有多少基础地