- 博客(36)
- 收藏
- 关注
RSS订阅原创 安装tf2.1
conda create -n TF2.1 python=3.7>conda activate TF2.1>pip install tensorflow==2.1.0 -i https://pypi.tuna.tsinghua.edu.cn/simple/pythonimport tensorflow as tftf.__version__
2021-06-07 17:54:06
235
原创 下载电离层产品
#!/bin/bashyear=2019sitename=${site[j]}for ((i=1;i<=365;i=i+1))dodoy=$(printf "%03d" "$i")yr=${year:2:4}filename=igsg${doy}0.${yr}i.Zwget --auth-no-challenge https://cddis.nasa.gov/archive/gnss/products/ionex/${year}/${doy}/${filename}gunzip $.
2021-04-21 19:36:41
780
原创 python小技巧
撤消操作：Ctrl+Z恢复操作：Ctrl+Shift+Z批量注释：Ctrl+/取消注释：Ctrl+/对其格式：选择代码按Tab取消对其格式：Shift+Tab
2021-04-15 17:12:39
101
原创 安装与卸载tensorflow2.2
pip install tensorflow-gpu==2.2 -i https://pypi.doubanio.com/simple/pip uninstall tensorflow
2021-04-08 17:27:12
270
原创 gamit10.72处理BDS数据
1/下载数据p文件ftp://igs.gnsswhu.cn/pub/gnss/mgex/daily/rinex3/2020/001/20p/o文件ftp://ftp.geodetic.gov.hk/rinex3/2020/001/hkpc/30s/ftp://gssc.esa.int/gnss/data/daily/2020/001/sp3文件ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/2086/ftp://gssc.esa.int/gnss/
2021-01-17 12:06:38
814
原创 使用sh_metutil生成采样300秒的ztd
#!/bin/bashfor ((i=1;i<=365;i=i+1))dodoy=$(printf "%03d" "$i")cd ${doy}sh_metutil -f o*a.${doy} -z z*.${doy} -i 300cd ..done
2020-12-24 17:05:27
276
2
原创 提取igs产品zpd
%提取zpdclear all;clc;day='缺失天数'; %定义变量名day，记录缺失的文件年积日ZPD = [];for n=1:31 %zpd文件天数 ns=num2str(n,'%03d'); % 输值转字符串，固定三位数字（001,010,100） name=['hksl',ns,'0.19zpd']; % 字符串拼接得到文件名称 ffid=fopen(na.
2020-10-22 22:13:16
719
原创 更新gamit的tables表问题
网站下载说明：1/http://garner.ucsd.edu/pub/gamit/tables/（格网表文件不全）采用迅雷下载不了其他文件除了VMF文件但是可以采用filezilla下载，在虚拟机上和win10上均可以，亲测可用2/http://sopac.ucsd.edu/GAMIT.shtml直接点击目标文件即可下载但是发现有些文件不全，如Solar, Lunar, Nutation仅到2016年的3/ftp://everest.mit.edu//pub/GRIDS/（格网
2020-10-15 21:41:09
1788
1
原创 查找空met
% 查找是否存在空文件clear;clc;Empty=[];site = ['bjfs';'alic';'darw';'daej';'hkoh';'hknp';'hkst';'hkpc';'hksc';'hksl';'hkkt';'hkss';'hklt';'hkws';'t430'];for j=1:15 for year=19 yr = num2str(year); for i=1:9 doy = num2str(i); .
2020-10-12 22:20:05
84
原创 查找空rinex
% 查找是否存在空文件clear;clc;Empty=[];site = ['bjfs';'alic';'darw';'daej';'hkoh';'hknp';'hkst';'hkpc';'hksc';'hksl';'hkkt';'hkss';'hklt';'hkws';'t430'];for j=1:15 for year=19 yr = num2str(year); for i=1:9 doy = num2str(i); .
2020-10-12 22:19:27
80
原创 对流层延迟产品网站
igs产品：https://cddis.nasa.gov/archive/gnss/products/troposphere/zpd/ftp://cddis.gsfc.nasa.gov/gnss/products/troposphere/zpd/（Anonymous ftp service will be discontinued on October 31, 2020.）ftp://gssc.esa.int/gnss/products/troposphere_zpd/%以下内容取自igs官
2020-10-05 20:43:43
1862
原创 提取探空站pwv并存储
% ①数据提取pwv并存储clear;clc;year=2018;month=10;b=[];PWV=[];pwv=[];month_day=7;for day=1:month_day for time=[00,12];%以下将年月日一天的次数数字转为字符串 years=num2str(year); months=num2str(month); days=num2str(day); times=num2str(time).
2020-10-05 15:11:37
1111
10
原创 在gamit中对rinex3的转换
下载的观测源文件：HKKT00HKG_R_20180010000_01D_30S_MO.crx.gz在gamit中采用如下方式1/解压：gunzip HKKT00HKG_R_20180010000_01D_30S_MO.crx.gz结果：HKKT00HKG_R_20180010000_01D_30S_MO.crx2/重命名：把.crx命名成.dsh_rename_rinex3 HKKT00HKG_R_20180010000_01D_30S_MO.crx结果：hkkt0010.18
2020-10-04 11:55:02
1157
1
原创 GAMIT10.61_处理BDS数据
数据下载广播星历：ftp://ftp.cddis.eosdis.nasa.gov/pub/gps/data/campaign/mgex/daily/rinex3/2018/001/18p/ (chmod +w brdm0010.18p)精密星历：http://www.igs.gnsswhu.cn/index.php/Home/DataProduct/mgex.html（文件只读改为可读写：...
2020-08-23 09:11:25
595
1
原创 GAMIT基线解算_分部处理
分步处理2017年001天（chan shao bifs）数据准备工作：1）mkdir demo2）cd demo3）mkdir rinex brdc igs4）cd brdc两种下载方式：第一种方式：采用gamit自带的下载脚本sh_get_nav -archive sopac -yr 2017 -doy 001 -ndays 1第二种方式：采用...
2020-08-23 09:11:06
1757
1
原创 GAMIT基线解算_批处理
分步处理2017年001,002天（shao hkkt hkks，后面两个站点为香港CORS站）数据准备工作：1）mkdir demo2）cd demo3）mkdir rinex brdc igs4）cd brdc采用wget下载wget -c ftp://nfs.kasi.re.kr/gps/data/daily/2017/001/17n/brdc001...
2020-08-23 09:10:31
2037
2
原创 win10 1909安装虚拟机14时在与桌面交互文件时出现问题
VMware虚拟机复制文件卡死的解决以前用的是win10 1803版本，用VMware虚拟机复制文件时没有问题，然后前几天更新了win10到1909，导致VMware版本不适配，需要更新到15.1以上就才能解决这个问题。
2020-06-13 12:13:49
192
原创 计算VMF1和探空Tm的精度
clc;clear;name1=['chan','.txt'];[site,mjd,ah,aw,zhd,zwd,tm]=textread(name1,'%s%s%s%s%s%s%s%*[^\n]','delimiter',' ');%提取前5列数据n=length(tm);for i=1:n%在1-n行数据内判断数据是否为空，如果为空则跳出本次循环 Tm(i,:)=str2num(tm{i}); s=i/4+1; j=floor(s); if mod(i,4)==.
2020-05-09 16:08:57
354
原创 提取VMF1中的Tm
clc;clear;name=['y2018','.vmf1_g'];fid = fopen(name);a=[];a1=[];A1=[];if fid>0 i=0; while feof(fid) == 0 i = i+1; tline{i} = fgetl(fid); a=tline(i); if i>22 a1=cell2mat(a); if (a1(1.
2020-05-09 16:08:24
242
1
原创 初次使用gamit10.71批处理解算小总结
本次结算使用chan shao jfng三个测站，由于使用gamit自带的脚本下载rinex，sp3，brdc等文件时速度太慢（据说cddis给限速了），本文采用了wget（wget在cddis上下载速度还行），以及武汉大学igs数据中心（采用迅雷）下载的，采用批处理。1/刚输入命令结算提示没有sites.defaults，查看了下gamit10.71的安装包后发现，安装包里的tables.1...
2020-05-03 20:18:05
2459
6
原创 东北地区探空站
嫩江,50557,125.23,49.16,243.0伊春,50774,128.90,47.71,232.0哈尔滨,50953,126.76,45.73,143.0通辽,54135,122.26,43.60,180.0长春,54161,125.21,43.90,238.0延吉,54292,129.46,42.88,178.0临江,54374,126.91,41.71,333.0...
2020-04-23 11:00:09
611
原创 GNSS各大网站
http://ggosatm.hg.tuwien.ac.at/DELAY/SITE/GPS/2019/y2019.vmf1_g在测站点初提供ha,wa,zhd,zwd,tm,p,t,e# station name (8 characters)# modified Julian date# hydrostatic coefficient a# wet coefficient a# hyd...
2020-04-22 18:44:43
1749
原创 my_get_chanrinex
#!/bin/bashsitename=chanyear=2018for ((day=274;day<=280;day=day+1))doyr=${year:2:4}filename=${sitename}${day}0.${yr}d.Zwget --http-user=anonymous --http-passwd=zou_yuxue2016@126.com -nd htt...
2020-04-20 21:51:14
126
原创 my_get_gbmclk
#!/bin/bashwget --http-user=anonymous --http-passwd=zou_yuxue2016@126.com -nd http://vmf.geo.tuwien.ac.at/trop_products/GNSS/VMF1/VMF1_OP/yearly/y2018.vmf1_g
2020-04-20 21:50:41
86
原创 my_get_igs
#!/bin/bashfor (( weeks=2021;weeks<=2025;weeks=weeks+1))doday=0filename=${weeks}${day}.sp3.Zwget --http-user=anonymous --http-passwd=zou_yuxue2016@126.com -nd ftp://cddis.gsfc.nasa.gov/pub/gp...
2020-04-20 21:50:08
159
2
原创 my_get_zpd
#!/bin/bashsite=chanfor ((day=274;day<=280;day=day+1))dofilename=${site}${day}0.18zpd.gzwget --http-user=anonymous --http-passwd=zou_yuxue2016@126.com -nd ftp://cddis.gsfc.nasa.gov/gps/produc...
2020-04-20 21:48:47
185
原创 win7安装OpenCV
1、下载OpenCVhttps://opencv.org/releases/2、解压（即安装）双击软件后点击运行3、将解压后的文件夹opencv拷贝到任意目录下，如：D:\Program Files\opencv4、配置环境变量→→→把D:\Program Files\opencv\build\x64\vc15\bin添加到Path后面，依次点击确定后即完成此步...
2020-03-04 17:48:01
1273
原创 arcgis统计面内点的数量
1、arctoolbox→分析工具→叠加分析→空间连接2、双击生成的新图层图层和图层属性对话框1）符号系统→显示（选择任意色带）2）值字段→Join_Count3）色带4）添加所有值3、确定结果：...
2019-12-03 19:33:52
15108
原创 arcgis绘核密度分布图
1、从网站上获取相应区域POI以及1:100万底图数据2、统一坐标系ARCtoolbox→数据管理工具→投影和变换→投影3、统一环境地理处理→环境4、核密度分析spatial analyst 工具→密度分析→核密度分析（此步骤生成的栅格图范围为矩形，需进一步提取到底图的面域内）5、矩形栅格提取到底图面域内spatial analyst 工具→提取分析→按掩...
2019-11-29 16:54:44
10860
原创 m_map绘点
clear;clc;LATLIMS=[36 42];LONLIMS=[120 130];%边界范围m_proj('lambert','lon',LONLIMS,'lat',LATLIMS);%投影模式m_coast('patch',[0.1 0.5,0.5]);%，改变数字，颜色填充m_grid('box','fancy');%没有网格，边框相间[X,Y]=m_ll2xy(125,4...
2019-11-02 16:52:25
1142
原创 MATLAB安装m_map
http://www.eos.ubc.ca/~rich/map.html1、以上网址是m_map工具箱的网站，找到m_map工具箱文件下下来，名字叫gzipped tar-file, orzip archive下其中一个即可。2、解压，得到 m_map文件夹，把他放到你的MATLAB安装目录下，D:\Program Files\MATLAB\R2016a\toolbox3、...
2019-10-30 21:24:46
6005
1
原创 linux解压GNSS相关文件
1.解压.Z文件：gunzip brdm0010.18p.Z2.解压.gz文件：gunzip JFNG00CHN_R_20180010000_01D_30S_MO.crx.gz3.解压.18d文件：sh_crx2rnx -f jfng0010.18d4.文件只读改为可读写：chmod +w wum19821.sp3你可以使用chmod命令，为改文件提供其他的权限。u表示用户、a表示...
2019-07-13 16:25:26
483
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人