- 博客(36)
- 资源 (1)
- 收藏
- 关注
原创 Gee中的数字和矩阵混合运算方法
Gee中的数字和矩阵混合运算方法print(Math.cos(40 * Math.PI/180)) // Client-sideprint(ee.Number(40 * Math.PI/180).cos()) // Server-side (For ee.Number and ee.Image)https://gis.stackexchange.com/questions/359702/calculate-sine-and-cosine-in-google-earth-engine...
2021-09-27 14:30:40 622
原创 python栅格数据出图可视化
import numpy as npimport cartopy.crs as ccrsimport cartopy.mpl.ticker as ctickerimport matplotlib.pyplot as pltimport cartopy.feature as cfeaturefrom mpl_toolkits.basemap import Basemapimport cartopy.io.shapereader as shpreaderfrom cartopy.mpl.tick.
2021-06-22 15:27:29 2570 3
原创 (WWF)上的GIS数据集
世界自然基金会(WWF)上的GIS数据集,其数据网址为:https://www.worldwildlife.org/pages/conservation-science-data-and-tools
2021-06-21 15:29:59 302
原创 python文件操作
1.获取特定路径下面特定类型文件(特定后缀名)#使用glob包from glob import glob#文件夹路径input_dir = r''#获得input_dir下面的txt文件列表file_list = glob(input_dir + '*.txt')for each_txt in file_list: print(each_txt)#打印文件路径,此种方法得到的文件名为绝对路径,是全称2.获取特定路径下面的所有文件夹#使用os包import os#文件夹路径inpu
2021-06-10 15:20:48 1099 2
原创 scipy.optimize.curve_fit参数简介
scipy.optimize.curve_fit函数参数列表scipy.optimize.curve_fit(f, xdata, ydata, p0=None, sigma=None, absolute_sigma=False, check_finite=True, bounds=- inf, inf, method=None, jac=None, **kwargs)[source]使用非线性最小二乘拟合函数f。参数列表f:callable可调用模型函数,f(x, …)。它必须将自变量作为
2021-06-03 17:39:43 6349
原创 A atmospheric correction tool大气校正工具
S2 surface reflectance images, cloud free selected, were down loaded from the THEIA-LAND site and orthorectified and corrected for atmosphericeffects using the MAJA processor for cloud detection and atmospheric correction ( https://www.theia-land.fr , acce
2021-05-31 17:35:28 671
原创 获取卫星的过境时间
获取卫星的过境时间https://app.spectator.earth/http://www.satflare.com/home.asphttps://www.n2yo.com/?s=40697%7c39084%7c25682https://www.agi.com/about
2021-05-28 09:47:25 5775
原创 gee选择一年中不同时间的数据
gee选择一年中不同时间的数据//filter collection by range, day of year, and boundsvar collection = ee.ImageCollection('COPERNICUS/S2') .filterDate('2015-05-01', '2017-05-31') .filter(ee.Filter.dayOfYear(121, 151)) .filterBounds(GP);
2021-05-27 11:21:10 1573
原创 Chiaa
Chiahttps://www.chia.net/faq/#faq-1What is Chia?Chia was incorporated in August of 2017 to develop an improved blockchain and smart transaction platform. We are building the Chia Network to improve the global financial and payments systems. Chia is the
2021-05-25 18:37:58 226 1
原创 Compaq Visual FortranV 6.7 正式版
Compaq Visual FortranV 6.7 正式版https://www.greenxf.com/soft/214445.htmlCompaq Visual Fortran(Fortran语言编译工具)是一款专业实用流畅度高的Fortran语言编译软件。不知道怎么进行Fortran语言编译?那就快试试绿色先锋小编推荐的Compaq Visual Fortran正式版下载使用。他可以对矩阵和复数进行运算,而且接近数学公式的自然描述,此版本是Compaq Visual Fortran破解版,内附注
2021-05-21 10:26:09 5988 8
原创 下载哨兵2号数据
下载哨兵2号数据(Export Google Earth Engine RGB Sentinel-2 Imagery to Google Drive using Python API)# This is the cloud masking function provided by GEE but adapted for use in Python.def maskS2clouds(image): qa = image.select('QA60') # Bits 10 and 11 are
2021-05-19 18:39:23 1208
原创 gee去云处理Landsat、Sentinel和Modis影像
gee去云处理Landsat、Sentinel和Modis影像1.Landsat 和 MODIS影像的去云函数function cloudfree_mod09a1(image){ var qa = image.select('StateQA') var cloudState = bitwiseExtract(qa, 0, 1) var cloudShadowState = bitwiseExtract(qa, 2) var cirrusState = bitwiseExtract(qa
2021-05-19 16:23:48 9126 2
原创 Gee获取影像每个像元的经纬度,以及把他们处理成为Numpy矩阵变量
Gee获取影像每个像元的经纬度,以及把他们处理成为Numpy矩阵变量#获取影像的经纬度def getImgData(image): image = image.addBands(ee.Image.pixelLonLat()) roi = image.geometry() data = image.reduceRegion( reducer = ee.Reducer.toList(), geometry = roi, scale = 10, ma
2021-05-19 16:19:24 2777 4
原创 python做方程的参数拟合,计算RMSE和R
python做方程的参数拟合,计算RMSE和R#参数拟合import numpy, scipy, matplotlibimport matplotlib.pyplot as pltfrom scipy.optimize import curve_fitfrom scipy.optimize import differential_evolutionimport warningsxData=numpy.array([0,0,100,200,250,300,400], dtype="float"
2021-05-11 20:55:45 2519
原创 将Pandas DataFrame转换为字典
将Pandas DataFrame转换为字典该to_dict()方法将列名设置为字典键将“ID”列设置为索引然后转置DataFrame是实现此目的的一种方法。to_dict()还接受一个’orient’参数,您需要该参数才能输出每列的值列表。否则,{index: value}将为每列返回表单的字典。可以使用以下行完成这些步骤:>>> df.set_index('ID').T.to_dict('list'){'p': [1, 3, 2], 'q': [4, 3, 2], 'r': [
2021-05-10 18:57:07 18169 3
原创 Cubist software
Cubist software (version 2.05)Free Software DownloadsThe C source code for single-threaded Linux versions of our software is now available under the GNU General Public License. (These are the full versions, not limited demonstration versions.) The source
2021-05-10 10:17:14 1081
原创 International Soil Moisture Network
International Soil Moisture Networkhttps://ismn.geo.tuwien.ac.at/en/
2021-05-08 11:36:57 855
原创 Calculate NDVI from Sentinel-2 with Python API
Calculate NDVI from Sentinel-2 with Python API# get indexesdef getNDVI(image): # Normalized difference vegetation index (NDVI) ndvi = image.normalizedDifference(['nir','red']).rename("ndvi") image = image.addBands(ndvi) return(image
2021-05-05 12:16:08 205
原创 MODIS R包下载数据
MODIS R包下载数据generic source constructorgetAOD <- function(user = 'EarthDataUsername' ,pw = 'EarthDataPassword' ,dir = 'OutputDestinationFolder' ,product = 'MCD19A2' # default product for MAIAC AOD
2021-04-28 23:53:53 793
原创 R语言绘图
R语言绘图https://wklchris.github.io/R-plotting-basic.htmlR语言绘图笔记汇总https://zhuanlan.zhihu.com/p/25074456R三大绘图系统https://blog.csdn.net/qq_34386891/article/details/80447729?utm_medium=distribute.pc_relevant.none-task-blog-baidujs_baidulandingword-0&s
2021-04-28 18:15:17 96
原创 pandas数据类型中的空值处理
pandas包的一些常用处理import pandas as pd#找到空值的索引data['property'].isnull().counts#将空值转为定值data['property'] = data['property'].fillna('99999')#找出空值所在索引data[(data.property == '99999')].index.tolist()#删除固定行,1和2行data = data.drop([1,2])参考文献https://blog.cs
2021-04-28 10:15:41 1092
原创 The Remote Sensing Code Library
The Remote Sensing Code Library (RSCL)The Remote Sensing Code Library (RSCL) is a free online registry of software codes of interest to remote sensing scientists and engineers. The RSCL is indexed by IEEE and is citable by using the DOI assigned to each c
2021-04-26 15:04:01 253
原创 gee 把数据下载到本地
gee 把数据下载到本地#方法一from google.colab import filesdf.to_csv('filename.csv') files.download('filename.csv')#方法二#!pip install -U -q PyDrivefrom pydrive.auth import GoogleAuthfrom pydrive.drive import GoogleDrivefrom google.colab import authfrom google
2021-04-25 16:13:18 5720 1
原创 gee导出DataFrame为csv
gee导出DataFrame为csv# Import Drive API and authenticate.from google.colab import drive# Mount your Drive to the Colab VM.drive.mount('/gdrive')# Write the DataFrame to CSV file.with open('/gdrive/My Drive/foo.csv', 'w') as f: df.to_csv(f)https:/
2021-04-25 15:59:02 595
原创 gee 获取时间序列
gee 获取时间序列import eeee.Authenticate()ee.Initialize()collection = ee.ImageCollection('NASA/FLDAS/NOAH01/C/GL/M/V001')taken = collection.filterDate('2019-01-01', '2019-12-31').sort('system:time_start', False)poi = ee.Geometry.Point([112.621391, -7.983
2021-04-25 14:46:23 1624
原创 gee 获取固定行列的数据
gee 获取固定行列的数据var L5fs1998post = ee.ImageCollection('LANDSAT/LT5_SR') //Landsat 5 Surface reflectance.filterDate('1998-1-9', '1998-2-27') //filter to date bounds .filter(ee.Filter.eq('wrs_path', 91))//filter to path and row.filter(ee.Filter.eq('wrs_row'
2021-04-23 16:05:21 438
原创 MODIS R
MODIS包也在github上,上面有常见问题模块,可以向开发者提问和提供bug.1,在EARTHDATA注册用户。到官方网站注册https://earthdata.nasa.gov/,比如你的用户名和密码分别是:用户名:james;密码:lileihanmeimei1232,安装MODIS包install.packages(“MODIS”)3,下载数据流程3.1 配置用户名和密码第一个是earthdata的账户配置,在MODIS包中 EarthdataLogin函数来配置,假如你的earth
2021-04-22 17:41:38 775
原创 NDVI 300m 10-day
Copernicus Global Land ServiceDelay in Sentinel-1 and Sentinel-3 based products(https://land.copernicus.eu/global/content/delay-sentinel-1-and-sentinel-3-based-products)Published on:2021-03-12Due to a fire at the ESA ground segment infrastructure, h
2021-04-22 16:27:18 168
原创 Water Cloud Model
Water Cloud Modelsoftware packageRhttps://cran.r-project.org/web/packages/WCM/index.htmlpythonhttps://github.com/multiply-org/forward-operators水云模型https://forum.step.esa.int/t/water-cloud-model/16456ATTEMA, E. P. W. & ULABY, F. T. 1978. Vegeta
2021-04-22 15:18:02 2504 4
原创 时间序列滤波之软件TIMESAT
时间序列滤波之软件TIMESAT1.TimeSat安装包2.MATLAB2018b(matlab每年分上下半年发布两个版本,一个是上半年的a,和下半年的b。建议通常用b,因为通常更稳定)(虽然该博文发布的时间,正逢美国对我国进行实体清单限制的当下。尤其是MATLAB授权。只希望国产软件借此东风,扶摇直上)TimeSat是一款不错的时间序列滤波软件,很专业,在提取植被物候参数方面有着很好的应用,本文主要讲解TimeSat 3.3软件在Windows和Matlab平台的安装及TimeSat 3.3软件的使
2021-04-22 09:43:04 2263
原创 gee 日期
gee 日期#GEE中的日期类#定义方式一date = ee.Date("2019-7-1")print(date.getInfo())#格式化输出print(date.format("yyymmdd").getInfo())#定义方式二date = ee.Date.fromYMD(2019,7,1)print(date.format("yyymmdd").getInfo())#定义方式三date = ee.Date.parse(""yyyymmdd,"20190701")pri
2021-04-21 11:10:55 862 1
原创 gee system time 转换为 date
gee system time 转换为 date// Random dummy Imagevar image = ee.Image("COPERNICUS/S2_SR/20200401T001611_20200401T001605_T01WCR");// system:time_start is saved as miliseconds (epoch):var epoch = image.get("system:time_start")print("Epoch system:time_start
2021-04-21 10:24:18 1514
转载 Google Earth Engine(区域统计)
Google Earth Engine(区域统计)提示:在通过GEE计算遥感的地表参量以后,我们就需要对计算出来的数据进行统计。GEE上面画图的功能并不是很齐全,得到的效果也不尽如人意。因此我们就需要将GEE对区域的统计量进行导出,导入到本地以后,再进行绘图。本期的的主要内容就是利用GEE对区域内的地表参量(如植被覆盖度、叶面积指数、蒸散量,本次以NDVI为例)进行统计(如平均值、面积等)。最后把统计的值导出到本地。var table = ee.FeatureCollection("users/b
2021-04-19 15:54:19 1092 1
原创 Google earth engine重采样
Google earth engine重采样代码.代码.# An highlighted blockimport eeee.Initialize()#获取landsat数据image_30m = ee.Image('LANDSAT/LC08/C01/T1/LC08_044034_20170614')#获取其中一个波段band2 = image_30m.select('B2')proj = band2.projection().getInfo()crs = proj['crs'#重采样为
2021-04-19 09:28:59 1127
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人