- 博客(38)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 M1 Mac打开Jupyter notebook
该问题是个路径问题，通过添加PATH环境变量就行了，设置环境变量时需要注意，zshrc和bash_profile中都可以设置，二者的区别是：每打开一个新的终端，zshrc就会被source（也就是设置的PATH被刷新一次），而bash_profile每次都需要手动source。将 if sys.platform[:3] == "win": 改成 elif sys.platform[:3] == "win": 就可以了。这样的操作过后，就不是command找不到的问题了，而是。
2024-10-14 21:02:37
378
原创 解决git@github.com: Permission denied (publickey).fatal: Could not read from remote repository.
git clone过程中遇到的问题git@github.com: Permission denied (publickey).fatal: Could not read from remote repository.
2022-11-01 10:41:37
969
原创 遥感影像叠加在谷歌地球（Google earth）上
谷歌地球可以直接打开tif文件，但是效果很差，会默认转成灰度图，如果想把五颜六色的tif（用了自己设置的colorbar）叠加在谷歌地球上该如何操作呢？ENVI提供了一个工具（SPEAR Google Earth Bridge）可以解决该问题。Add files把自己的文件导入进来，然后一直next就可以导出一个kml了，结束后默认会打开谷歌地球。效果如图：...
2022-03-08 16:57:49
3829
6
原创 PIE engine初体验
PIE engine是GEE的国产版本（由航天宏图开发）。语法二者基本一致又不完全相同，数据集PIE engine正在迅速扩大，其优势是：1、不需要科学上网即可实现遥感数据云处理2、有很多国产数据集（如高分系列）3、在网页coding时，会自动提示 相应的函数和变量劣势是：1、目前数据集时间序列较短，有很多都是从2020年开始，且有的仅限于某些地区。2、很小的任务经常会报一些用户权限不够的错误，但是同样的任务在GEE却可以正常运行。网址：PIE-Engine 遥感与地理信息云
2021-09-19 15:31:30
2159
原创 根据卫星轨道参数（Two Line Element）模拟星下点轨迹
最方便获取的卫星轨道参数是北美防空司令部（NORAD）发布的两行根数（Two Line Element）主要根数记录在第二行：（inclination）轨道倾角，升交点经度，轨道偏心率，近地点角距，平近点角，（mean motion）平运动速率。星下点轨迹绘制clcclose allclear all%%W_EARTH =7.29e-5; %地球自转速度GGC =3.986e5;%地球引力常数%轨道六根数a =6517.3;%6717.3; %地球半径e =0.0031..
2021-09-13 13:43:29
9493
7
原创 javascript打印
这里写自定义目录标题欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入欢迎使用Markdown编辑器你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Mar
2021-07-31 17:51:09
85
原创 GEE区域分析（生成时间序列 time series）
对于单个省生成时间序列//单个省的时间序列var neimeng = ee.FeatureCollection('users/peisipand/cn_province') .filter(ee.Filter.eq('OBJECTID',5)).geometry();Map.addLayer(neimeng)var CH4_col = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_CH4')
2021-05-27 22:05:32
4993
2
原创 GEE 区域统计（求各省的中值并做进一步处理）
目标：对中国的CH4（TROPOMI）提取异常值做法：对每个省的CH4求中值，然后用该省的CH4影像减该中值，每个省都这样做，最终获得一张全国的影像。最帅的代码：var china = ee.FeatureCollection("users/peisipand/cn_province");var image = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_CH4') .select('CH4_column_volume_mixing_ratio_dry
2021-05-27 16:09:46
2574
4
原创 从txt中获取廓线
SFIT提供了温湿压以及各种气体的垂直廓线。clc;clear;filename = 'reference.txt';fid = fopen(filename,'r');all = {'alt','pre','tem','H2O','CO2','O3','N2O','CO','CH4','O2'};if fid<0 warndlg('打开文件失败!'); return;else for i=1:length(all) FormatString=rep
2021-04-24 22:02:51
182
原创 关于GEE导出数据时 含空值的栅格图像
GEE提供的卫星数据大多数是有丢失的，若不加操作直接导出为tif，matlab读出来的矩阵该是空值的都是nan，然而arcgis并不能识别空值为nodata，导致制图非常难受。加上unmask函数可以让空值变为某个数字，这里设为-9999，方便后续处理。selfMask()则可以将图像的负值都判定为空。Export.image.toDrive({ image: xinjiang_CH4_anomaly.unmask(-9999), description: 'xinjiang_ch4_anom
2021-04-20 22:29:05
4274
10
原创 GEE每7天平均影像并输出（以tropomi为例）
var shp = ee.FeatureCollection("users/peisipand/cn_province");// Define an area to perform interpolation over.var yunnan = shp.filter(ee.Filter.eq('OBJECTID', 25));Map.addLayer(ee.Image().paint(yunnan, 0, 2), {}, 'yunnanborder');var first_date = ee.D
2021-04-13 15:22:04
1549
9
原创 GEE栅格插值（以TROPOMI的NO2数据为例）
GEE插值GEE提供了IDW和Kriging两种插值var shp = ee.FeatureCollection("users/peisipand/cn_province");// Define an area to perform interpolation over.var yunnan = shp.filter(ee.Filter.eq('OBJECTID', 25));Map.addLayer(ee.Image().paint(yunnan, 0, 2), {}, 'yunnanbor
2021-04-12 22:03:06
2973
6
原创 wireguard服务器配置
虚拟主机连接好之后，先运行下面这个指令来更新主机的配置软件：apt-get update -y && apt-get install curl -y#!/bin/bashrand(){ min=$1 max=$(($2-$min+1)) num=$(cat /dev/urandom | head -n 10 | cksum | awk -F ' ' '{print $1}') echo $(($num%$max+$min)) }wiregu
2021-02-27 17:48:57
5487
2
原创 OMI数据下载和绘图
下载OMI数据（以对流层NO2为例）：OMI-NO2下载地址选择想要的时间范围，筛选，download all然后大概是这个界面存为txt，上传至超算。wget --user=***** --password=******* -i NO2-download.txt处理成月均数据：clc;clear;datadir='/scratch/zppei/OMI/NO2/'; %指定批量数据所在的文件夹month = string(zeros(12,1));for j = 2015:20
2021-02-27 17:26:40
4958
10
原创 GEE常用语句
注释代码快捷键command + /计算表达式add subtract multiply dividevar step1 = image.expression( 'R75>2.5 && D75>0.3 && B7>0.5',{ 'D75':image.select('B7').subtract(image.select('B5')), 'B7':image.select('B7'), 'R75':imag..
2021-01-22 16:14:40
1269
1
原创 GEE非监督分类，生成图例
// Load a pre-computed Landsat composite for input.var input = ee.Image('LANDSAT/LE7_TOA_1YEAR/2001');// Define a region in which to generate a sample of the input.var region = ee.Geometry.Rectangle(29.7, 30, 32.5, 31.7);// Display the sample region.
2021-01-22 16:12:44
766
原创 批量下载ECWMF-INTERM数据
ECMWF（European Centre for Medium-Range Weather Forecasts 简称ECMWF），下载地址为http://apps.ecmwf.int/datasets/data/interim_full_daily/账号注册完之后，打开https://api.ecmwf.int/v1/key/查看API KEY1.安装ECMWF KEY（三种方法选其一）https://confluence.ecmwf.int/display/WEBAPI/Access+E..
2021-01-05 13:39:00
1052
原创 The NetCDF library encountered an error during execution of ‘open‘ function
>> >> >> >> {Error using netcdflibThe NetCDF library encountered an error during execution of 'open' function -'HDF error (NC_EHDFERR)'.Error in netcdf.open (line 67) [varargout{:}] = netcdflib ( 'open', filename, varar
2020-10-26 21:37:25
1629
原创 利用mmap工具箱显示shp、tif和散点图（m_map）
close all, clear all, clcfpni='E:\MeteoInfo\map\china.shp';%中国地图shp文件所在路径China=shaperead(fpni);boux=[China(:).X];bouy=[China(:).Y];%分别是获取经度X信息和纬度Y信息set(gcf,'position',[0 0 1440 780]);%设置图形窗口位置和大小m_proj('Mercator','lon',[-170 170],'lat',[-80 80]);%设置投影
2020-09-21 14:59:58
2088
2
原创 武汉大学超算中心运行matlab程序
武汉大学超算中心运行matlab程序以前的超算运行matlab比较容易，写好.m文件，复制个matlab_mdcs_slurm.sh改一下，输入指令./matlab_mdcs_slurm.sh就ok了。可参考如何正确优雅地使用武大超算进行MATLAB计算新的改变昨天还能运行的程序，今天就报错。我￥%*&#错误如下：./matlab_mdcs_slurm.sh: line 42: 46789 Killed意思大概是程序没问题，提交上去之后被杀死了。问了超算管理员之后，得到的回复是好
2020-08-28 19:41:35
2654
原创 OCO数据批量下载
https://earthdata.nasa.gov/以OCO2数据为例：这里可以选择筛选时间和空间范围wget --user=**** --password=**** -i7308445156-download.txt (****替换为对应账号密码，后面是含下载链接的文件名)
2020-08-25 16:47:50
2941
6
原创 运行WRF
cd /project/hebin/Build_WRF/WPS./geogrid.exe 在WPS目录生成相应的 geo_em.d0x.nc 文件ncl ./util/plotgrids_new.ncl 绘制研究区./link_grib.csh ../DATA/riyadh/fnl*./ungrib.exe >& ungrib_data.log./metgrid.execd /project/hebin/Build_WRF/WRF/test/em_real.
2020-08-25 07:12:54
1461
原创 WRF气象数据之NECP FNL数据批量下载
https://rda.ucar.edu/datasets/ds083.21.在网站里找到FNL数据2.下载脚本3.linux系统下输入指令./download.csh password（后面跟密码）
2020-08-25 06:11:07
9382
4
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人