基於Python進行小波分析

来源:https://www.cnblogs.com/qsgeo/p/18330961
-Advertisement-
Play Games

在氣象學和環境科學的研究中,理解和預測氣象數據的周期性變化至關重要。小波分析作為一種高效的數學工具,近年來在氣象數據的周期性分析中得到了廣泛應用。本文將詳細介紹如何通過Python進行小波分析,以探究氣象數據中的周期性變化。 ...


在氣象學和環境科學的研究中,理解和預測氣象數據的周期性變化至關重要。小波分析作為一種高效的數學工具,近年來在氣象數據的周期性分析中得到了廣泛應用。本文將詳細介紹如何通過Python進行小波分析,以探究氣象數據中的周期性變化。

1 數據來源及下載方式

西北農林科技大學的彭守璋研究員在國家青藏高原科學數據中心公開發佈了氣溫、降水、乾燥度等氣象數據,本文所使用的數據為基於其中的 中國1km解析度逐月降水量數據集 轉換得到的某地區歷年降水量數據,可通過FTP進行下載。

Filezilla 桌面版軟體為例,輸入官網中給定的主機號、用戶名、密碼、埠等信息後即可將數據從遠程站點下載至本地。

image

2 代碼編寫

2.1 導入相關模塊

import openpyxl
import pywt
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from matplotlib import rcParams

2.2 設置地圖字體及字型大小

config = {
    "mathtext.fontset": "stix",
    "font.family": "serif",
    "font.serif": ["Times New Roman"],
    "font.size": 24,
    "axes.unicode_minus": False # 處理負號,即-號
}
rcParams.update(config)

2.3 從Excel表讀取數據

# 通過pandas模塊讀取Excel表,並提取表格中的年份及降水量數據
file_path = "file_path"	# 此處為Excel文件路徑
data = pd.read_excel(file_path)
Year_list = list(data["Year"])
Pre_list = list(data["Pre"])

2.4 進行小波變換,繪製小波等值線圖

scales = np.arange(1, 31)	# 設置小波分析的時間尺度
coef, freqs = pywt.cwt(Pre_list, scales, "morl")	# 對降水量數據進行Morlet小波變換

fig, ax = plt.subplots(figsize=(15, 10))

# 繪製小波繫數圖像
im = plt.imshow(abs(coef), extent=[Year_list[0], Year_list[-1], 30, 1],
                interpolation="bilinear", cmap="gray", aspect="auto",
                vmax=abs(coef).max(), vmin=-abs(coef).max())
plt.colorbar(im, ax=ax)  # 添加色帶圖例

# 繪製小波繫數等值線
contour = plt.contour(Year_list, scales, coef.real, colors="black", linewidths=1)

ax.invert_yaxis()   # 反轉Y軸,使時間尺度從大到小排列
ax.set_xticks(np.arange(1901, 2024, 20))
ax.set_yticks(np.arange(5, 31, 5))
ax.set_xlabel("Year")
ax.set_ylabel("Scale")
fig.savefig("小波等值線圖.jpg", dpi=600)

2.5 計算小波方差,繪製小波方差圖

variance = np.var(coef, axis=1)	# 計算小波方差
variance1 = variance / 10000	# 數值較大,將數值縮小10000倍後顯示在結果圖中


fig, ax = plt.subplots(figsize=(15, 10))

ax.plot(scales, variance1, "k-")	# 繪製方差曲線
ax.set_xlabel("Scale (a)")
# \times為LaTeX語法中的叉乘號
# 在Python中,反斜杠"\"為轉義字元,此處需同時輸入兩個反斜杠
ax.set_ylabel("Variance $\mathrm{(\\times 10^4)}$")
ax.set_xticks(np.arange(0, 31, 5))
ax.set_yticks(np.arange(0, 9, 2))
fig.savefig("小波方差圖.jpg", dpi=600)

3 等值線及方差圖示例

image

image


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 引言 在JavaScript開發中,設計模式是解決特定問題的有效手段。單例模式(Singleton Pattern)是其中一種常見且有用的模式。儘管網上有許多關於單例模式的解釋和實現,本篇將從實際工作中的需求出發,探討如何更好地理解和應用單例模式,以編寫更復用、更高效的代碼。 什麼是單例模式? 單例 ...
  • 手寫 Hibernate 系列 手寫 Hibernate ORM 框架 00-hibernate 簡介 手寫 Hibernate ORM 框架 00-環境準備 手寫 Hibernate ORM 框架 01-註解常量定義 手寫 Hibernate ORM 框架 02-實體 Bean 定義,建表語句自動 ...
  • 本文主要介紹了在使用Python進行面向對象編程時,異常的層級和如何使用繼承關係完成自定義自己項目中異常類,並以感測器數據採集為例進行講解。 ...
  • 拓展閱讀 Java Servlet 教程-20-自己手寫實現 spring mvc 整體思路 Java Servlet 教程-21-自己手寫 spring mvc 簡單實現 Spring Web MVC-00-重學 mvc mvc-01-Model-View-Controller 概覽 mvc-02 ...
  • 題目要求 給定兩個字元串 s 和 t ,編寫一個函數來判斷 t 是否是 s 的字母異位詞。 註意:若 s 和 t 中每個字元出現的次數都相同,則稱 s 和 t 互為字母異位詞。 由於字元在電腦記憶體中是以ASCII碼或Unicode編碼的形式存儲的,我們可以得出'a'在ASCII表中的值是97,'A ...
  • PART1: Java基礎知識概述與Java的下載安裝 1)Java語言概述: ① Java的發展史: 詹姆斯·高斯林(James Gosling) 1977年獲得了加拿大卡爾加里大學電腦科學學士學位,1983年獲得了美國卡內基梅隆大學電腦科學博士學位,畢業後到IBM工作,設計IBM第一代工作站 ...
  • 1、簡述C++中命名空間的作用。 答:避免重覆定義全局變數的問題。 2、定義兩個命名空間A 和 B 分別在A中和B中定義變數value。在main函數中將兩個空間的value列印出來。 #include "iostream" using namespace std; namespace A { in ...
  • 寫在前面 昨天說了一下Java中的數據類型、運算符、選擇語句、迴圈語句部分的基礎知識,今天寫的編程題就是來檢驗這部分的成果,來看看你能寫出來幾題。答案也是僅供參考,如果有更好的解法歡迎在下麵留言! 題目展示 1.輸入自己的名字,年齡和性別,分別用不同的變數接收,並將輸入的信息做輸出。 代碼: pub ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發(一):使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說,近來想嘗試開發一款安卓APP,考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP,畢竟是使用Visual Studio開發工具,使用起來也比較的順手,結合微軟官方的教程進行了安卓 ...
  • 前言 QuestPDF 是一個開源 .NET 庫,用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫,它徹底改變了我們生成 PDF 文檔的方 ...
  • 項目地址 項目後端地址: https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址: ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址: http://tree ...
  • 話不多說,直接開乾 一.下載 1.官方鏈接下載: https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe,運行開始下載SQL server; 二. ...
  • 前言 隨著物聯網(IoT)技術的迅猛發展,MQTT(消息隊列遙測傳輸)協議憑藉其輕量級和高效性,已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫,為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog支持多種接收器用於日誌存儲,增強器用於添加屬性,LogContext管理動態屬性,支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項,適用於不同需求。 ...
  • 目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據,爬取信息時需要模擬瀏覽器行為,否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下: 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • 1.前言 什麼是熱更新 游戲或者軟體更新時,無需重新下載客戶端進行安裝,而是在應用程式啟動的情況下,在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR,Xlua,ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles,Addressable, ...
  • 本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺,其實像阿裡雲,騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後,它會送 ...
  • 通過以下方式可以高效,並保證數據同步的可靠性 1.API設計 使用RESTful設計,確保API端點明確,並使用適當的HTTP方法(如POST用於創建,PUT用於更新)。 設計清晰的請求和響應模型,以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證,確保接收到的數據符合預期格 ...