最近工作需要用到最遠(yuǎn)點(diǎn)采樣fps算法, 找了一些blog學(xué)習(xí)學(xué)習(xí), 寫了個(gè)原型: 測(cè)試一下, 跑的沒問(wèn)題, 但是速度感人 所有點(diǎn)的歐氏距離計(jì)算做個(gè)向量化優(yōu)化: 這一下子快了1...
發(fā)簡(jiǎn)信
IP屬地：云南
-
FarthestPointSampling(FPS)算法學(xué)習(xí)&實(shí)踐
-
記一個(gè)gcc7代編譯器下allocate-source中有關(guān)array bounds的離譜bug
Bug情況 最近寫Fortran代碼的時(shí)候遇到了個(gè)怪bug, 代碼如下: 然后離譜的來(lái)了, 我發(fā)現(xiàn)在gcc750和gcc710里面跑出來(lái)的shape_data的lbound都...
-
cmake - 編寫單元測(cè)試在cmake中構(gòu)建測(cè)試?yán)褪莈nable_testing()+add_test()這一套。本篇文章是我在學(xué)習(xí)cmake構(gòu)建單元測(cè)試的一個(gè)總結(jié)筆記。 工程概覽 如下是我所定義的...
-
WPS ungrib編譯報(bào)錯(cuò) undefined reference to 'jpc_decode'
最近在一個(gè)新機(jī)子上用Spack配置好編譯器, hdf5, netcdf-c, netcdf-fortran, jasper等一系列依賴庫(kù)后編譯了最新版本的WRF-4.4.2和...
-
地理坐標(biāo)轉(zhuǎn)換-proj4
目標(biāo)：WGS-84 轉(zhuǎn)為 CGCS2000坐標(biāo) 使用proj4 測(cè)試數(shù)據(jù) 舉例：天安門在不同類型坐標(biāo)系下的經(jīng)緯度值 WGS-84經(jīng)緯度：116.391349,39.90737...
-
但是這不是netCDF4嗎??
利用xarray寫netcdf文件
三維模式輸出的數(shù)據(jù)很大，傳輸起來(lái)很麻煩。有時(shí)候我們只需要幾個(gè)關(guān)鍵變量，因此把他們抽取出來(lái)寫成nc文件比較方便。 例如，要單獨(dú)保存megan輸出的isoprene變量isop為...
-
EBAS觀測(cè)數(shù)據(jù)處理
記錄一下 EBAS-Thredds[https://ebas.nilu.no/thredds/]數(shù)據(jù)的處理流程, 以全站點(diǎn)的ozone為例 基于前一個(gè)EBAS觀測(cè)數(shù)據(jù)下載[h...
-
m3diff對(duì)比分析ioapi-netcdf文件
最近要在一個(gè)新機(jī)器上跑CMAQ, 先對(duì)機(jī)器和模式做個(gè)校驗(yàn), 以前都是跑出來(lái)自己寫個(gè)ncl腳本比對(duì)的, 現(xiàn)在覺得用ioapi的m3diff也挺好, 不過(guò)還需要基于那個(gè)結(jié)果進(jìn)行后...
-
EBAS觀測(cè)數(shù)據(jù)下載
記錄一下 下載EBAS-Thredds[https://ebas.nilu.no/thredds/]數(shù)據(jù)的流程, 以全站點(diǎn)的ozone為例 1. 下載ebas-thredds...
-
CFSR 6小時(shí)數(shù)據(jù)批量下載腳本-WRF模擬用
WRF官方[https://www2.mmm.ucar.edu/wrf/users/download/free_data.html]提供的CFSR數(shù)據(jù)下載地址是在UCAR[h...
-
FORTRAN調(diào)用C函數(shù)的參數(shù)不一致問(wèn)題-自動(dòng)額外傳參特性
最近在debug的時(shí)候看到CMAQ ioapi的code, envgets.c中和ENVSTR函數(shù)和m3utilio.F中的接口參數(shù)不一致, 但是可以正常跑, 之前沒注意過(guò),...
-
CS61C-2022-su20作業(yè)記錄
記錄一下CS61C 20su的各個(gè)lab和proj的心得, 代碼放在了這里[https://github.com/Roadelse/CS61C-2022-su20] proj...
-
ImageMagick日常使用格式記錄
記錄一些日常用ImageMagick處理圖像的格式 1. 剪切圖的一部分 magick convert src.png -gravity East -crop 200x300...
-
觀測(cè)數(shù)據(jù)處理-國(guó)控環(huán)境監(jiān)測(cè)站點(diǎn)
工作需要用到我國(guó)環(huán)境監(jiān)測(cè)站點(diǎn)的污染物濃度數(shù)據(jù)(感謝大佬的分享), 數(shù)據(jù)很全, 不過(guò)csv格式在分析的時(shí)候尤其大量數(shù)據(jù)分析的時(shí)候并不友好, 所以一般要二次處理一下 上學(xué)的時(shí)候搞...
-
powershell數(shù)組賦值遇到的bug, 開頭的逗號(hào)不是一直都能用的
今天寫Powershell的時(shí)候遇到個(gè)bug, 搜了半天沒搜到(不過(guò)最后解決后感覺以前似乎看到過(guò), 可能我搜的內(nèi)容不對(duì)) 簡(jiǎn)單地說(shuō), powershell弄個(gè)數(shù)組一般就直接逗...
-
numpy slice broadcast的坑, 和ncl或者xarray不同
過(guò)去用ncl作array slice習(xí)慣了, 不同的dimension之間就是直接乘起來(lái), 例如 即可得到一個(gè)2x2的數(shù)組 結(jié)果今天在用numpy處理數(shù)據(jù)的時(shí)候發(fā)現(xiàn) 就只能得...
-
Sublime修改Build Python的默認(rèn)指令
網(wǎng)上查到的配置sublime直接運(yùn)行python腳本的方法大多是新建一個(gè)Build python的設(shè)置, 而直接以默認(rèn)的Python build運(yùn)行的話可能會(huì)報(bào)錯(cuò)提示找不到指...
-
好像目前是隱藏了中國(guó)站點(diǎn)..但是手動(dòng)輸入站點(diǎn)編號(hào)比如北京是54511, 歷史數(shù)據(jù)還是可以看的
不顧20年的似乎不行
另外歷史數(shù)據(jù)我之前都down了,有需要的話可以發(fā)給你下載wyoming大學(xué)的探空數(shù)據(jù)
水平有限,歡迎指正交流,共同進(jìn)步! 俄懷明大學(xué)將全球探空站的數(shù)據(jù)共享使用,并通過(guò)計(jì)算了很多的變量指數(shù).非常有用,近期有這方面的需求,暫時(shí)寫了一段未經(jīng)優(yōu)化的代碼下載.如果要大量...
-
mpi4py 快速上手
在上一篇中我們介紹了如何安裝和使用 mpi4py，下面我們以幾個(gè)簡(jiǎn)單的例子來(lái)展示怎么使用 mpi4py 來(lái)進(jìn)行并行編程，以使讀者能夠快速地上手使用 mpi4py。這些例子來(lái)自...