Python3.0科學計算學習之繪圖(二)

来源:https://www.cnblogs.com/chenzhijuan-324/archive/2019/04/08/10673635.html
-Advertisement-
Play Games

(1) np.mashgrid()函數: 生成網路點坐標矩陣,可以是二維網路矩陣,也可以是三維網路矩陣。其中,每個交叉點就是網路點,描述這些網路點的矩陣就是坐標矩陣(橫坐標矩陣X中的每個元素與縱坐標矩陣Y中對應位置元素,共同構成一個點的完整坐標)。 背景示例:網路點與坐標矩陣的解釋如下: impor ...


(1) np.mashgrid()函數:-----生成網路點坐標矩陣,可以是二維網路矩陣,也可以是三維網路矩陣。其中,每個交叉點就是網路點,描述這些網路點的矩陣就是坐標矩陣(橫坐標矩陣X中的每個元素與縱坐標矩陣Y中對應位置元素,共同構成一個點的完整坐標)。

背景示例:網路點與坐標矩陣的解釋如下:

import numpy as np

import matplotlib.pyplot as plt

x=np.array([[0,1,2],[0,1,2]])        #最簡單的方法是,可以把橫縱坐標矩陣X,Y寫出來,生成坐標點

y=np.array([[0,0,0],[1,1,1]])

plt.plot(x,y,color='red',marker='.',markersize=10,linestyle='-.')      #點的形狀為原點,點設置大一些,線性為電劃線。

plt.grid(True)

plt.show()   

註意:按照矩陣給坐標點信息,matplotlib會把橫坐標矩陣中,每一列的點當做同一條線。正如上例中把plotlinestyle=’ ‘改為linestyle=’-.’就會發現A-D,B-E,C-F是連接的。

                                  

對於很多網路點的情況,可用如下meshgrid()函數方法:

   因註意到坐標矩陣其中有大量的重覆---X的每一行都一樣,Y的每一列都一樣。故基於此規律性,numpy提供的np.meshgrid()函數可以快速生成坐標矩陣X,Y

x=np.linspace(-0.5,2.,10)    

y=np.linspace(-1.5,4.,10)

X,Y=np.meshgrid(x,y)     #輸入的x,y是網路點的橫縱坐標列向量(非矩陣),輸出的XY就是坐標矩陣。

plt.plot(X,Y,color='limegreen',marker='.',linestyle='')

plt.grid(True)

plt.show()

                   

(2) Python可視化庫matplotlib.pyplotcontour()contourf()函數

區別contour()counterf() 函數功能相同,都是畫三維等高線圖的,不同點在於contourf會根據不同的水平值用不同的顏色來填充繪圖(即會對等高線間的區域進行填充)。

contour()counterf() 都是繪製三維圖,其中前兩個參數xy為兩個等長一維數組,第三個參數z為二維數組(表示平面點xiyi映射的函數值)。

matplotlib等高線的繪製步驟:

等高線是三維圖像在二維空間的投影

首先準備三維函數及待投影平面的網格坐標。

import numpy as np

import matplotlib.pyplot as plt

def height(x,y):

    return (1-x/2+x**5+y**3)*np.exp(-x**2-y**2)

x=np.linspace(-3,3,300)

y=np.linspace(-3,3,300)

X,Y=np.meshgrid(x,y)

l 繪製等高線

C=plt.contour(X,Y,height(X,Y),10,colors='black')   # 10代表等高線的密集程度。如果是0,則圖像被一分為二

plt.contourf(X,Y,height(X,Y),10,alpha=0.75,cmap=plt.cm.hot)

#為等高線填充顏色 10表示按照高度分成10alpha表示透明度 cmap表示漸變標準

l 添加高度與數字

plt.clabel(C,inline=True,fontsize=10)     #inline=True表示高度寫在等高線上

plt.xticks(())     

plt.yticks(())

#去掉坐標軸刻度

plt.show()

#顯示圖片

l 最後結果圖為:

                                           

實例:尋找某一函數的最小值

import numpy as np

import matplotlib.pyplot as plt

import scipy.optimize as so

rosenbrockfunction=lambda x,y: (1-x)**2+100*(y-x**2)**2  #使用lambda關鍵字所表示的匿名函數來使代碼變得更加簡潔

x=np.linspace(-.5,2.,100)

y=np.linspace(-1.5,4.,100)

X,Y=np.meshgrid(x,y)

Z=rosenbrockfunction(X,Y)

cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10),cmap='gray')   #以第四個參數給出的值為標準來繪製等高線並用gray顏色繪圖   

#plt.contourf(X,Y,Z,10,alpha=0.75,cmap='gray')

rosen=lambda x:rosenbrockfunction(x[0],x[1])

solution,iterates=so.fmin_powell(rosen,x0=np.array([0,-0.7]),retall=True)

#迭代方法so.fmin_powell利用鮑威爾方法來尋找最小值(由一個給定的起始值x0開始,並且當給出選項retall=True時報告所有的迭代,經過16次迭代後得出x=0,y=0

x,y=zip(*iterates)

plt.plot(x,y,'ko')

plt.plot(x,y,'k:',linewidth=1)

plt.title('Steps of powells method to compute a minimum')

plt.clabel(cs)

#clabel使用該contour創建的一個等值線集對象來註釋對應函數值的水平值

補充:np.logspace()常用於創建等比數列,它也有常用的3個參數,第一個參數表示起始點的指數,第二個參數終止點的指數,第三個參數表示數列的個數;最後,通過base參數可以修改底數

                          

(4) 圖像

Imshow()函數用於在指定的視窗中顯示一幅圖像。它負責對圖像進行處理,並顯示其格式,但是不顯示圖片,其後跟著plt.show()才能顯示出來。

imshow圖像顯示函數:imshowf,G) 其中,f是函數,G是顯示該圖像的灰度級數。

                    Imshow(f,[low,high]) 小於或者等於low的值都顯示為黑色,大於或等於high的值都顯示為白色

                    Imshow(f,[ ]) [ ] 自動將low設為函數最小值,將high設為函數最大值。

Colorbar:增加顏色類標的代碼是plt.colorbar()

例如:import numpy as np

import matplotlib.pyplot as plt

X=np.array([[1,2],[3,4],[5,6]])

plt.imshow(X)

plt.colorbar()

plt.show()   

實例:將數組可視化為圖像,

import numpy as np

import matplotlib.pyplot as plt

# 下麵代碼演示瞭如何通過meshgrid生成複數參數值的矩陣:

def mandelbrot(h,w,maxit=20):

    X,Y=np.meshgrid(np.linspace(-2,0.8,w),np.linspace(-1.4,1.4,h))

    c=X+Y*1j

    z=c

    exceeds=np.zeros(z.shape,dtype=bool)

    for iteration in range(maxit):

        z=z**2+c

        exceeded=np.abs(z) > 4

        exceeds_now=exceeded & (np.logical_not(exceeds))

        exceeds[exceeds_now] = True

        z[exceeded]=2

    return (exceeds)   

#命令imshow將矩陣展示為圖像,所選的顏色貼圖顯示了序列出現的白色無界的區域,而其他區域顯示為黑色。

plt.imshow(mandelbrot(400,400),cmap='gray')

plt.axis('off')    #使用plt.axis('off') 關閉坐標軸

plt.show()       

# 使用imshow的無插值的圖像

                                       

註意:imshow命令預設使用插值使圖像看起來更好,這在矩陣較小時可以看的很清楚。

plt.imshow(mandelbrot(400,400),interpolation='nearest',cmap='gray')

Interpolation代表插值運算,nearest表示最近鄰近插值法,只是插值方式中的一種。

cmap表示繪圖時的樣式,這裡選擇的是ocean主題

#使用最近鄰近插值法所得的圖像

                                        

 


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

-Advertisement-
Play Games
更多相關文章
  • 1.可以把七層協議簡化成四層協議鏈路層 網路層 傳輸層 應用層 2.通過路由器連接的兩個網路網路層ip提供的是一個逐跳協議,提供了一種不可靠的服務,中間有可能會丟傳輸層tcp在ip的基礎上提供了可靠的傳輸層 比喻:tcp就是淘寶賣家 ,ip就是包裹,中間如果包裹丟了,賣家會重新發一個包裹,這裡會有一 ...
  • 鏈表:由一系列不必再記憶體中相連的結構組成,每一個結構均含有表元素和指向後繼結構的指針。 與數組、列表的主要區別: 記憶體不連續; 不能通過下標隨機訪問。 優點: 插入、刪除操作效率高,時間複雜度為o(1); 記憶體利用率高,不會浪費記憶體; 大小不固定,擴展靈活; 缺點: 隨機訪問性差,查找效率低,時間復 ...
  • 死磕 java集合之ConcurrentHashMap源碼分析(一) 它的存儲結構是什麼樣的? 它使用了哪些鎖? 它是怎麼擴容的? 它是否是強一致性的? 它不能解決哪些問題? 它的源碼中使用了哪些不常見的技術? ...
  • '''迭代器:兩個基本方法:iter()和next()迭代器是一個可以記住遍歷的位置的對象。 迭代器對象從集合等第一個元素開始訪問,直到所有的元素被訪問結束,迭代器只能往前不會後退。 迭代器有兩個基本的方法:iter()和next() 字元串,列表或元組對象都可以用於創建迭代器。 迭代器的一大優點是 ...
  • 1. 調用鏈Cat 1.1. 調用鏈演進 1.2. 開源產品比較 1.3. 監控場景 1.4. cat的增值作用 1.5. cat典型報表 1.5.1. 應用報錯大盤 1.5.2. 業務大盤 1.5.3. logView 1.5.4. 可視化的logView 1.5.5. 應用報表(APM) 1.5 ...
  • 一、請求鉤子 在客戶端和伺服器交互的過程中,有些準備工作或稍微工作是需要處理的,比如:在請求開始時,建立資料庫連接;在請求開始時,根據需求進行許可權校驗;在請求結束時,指定數據的交互格式等。為了讓每個視圖函數避免編寫重覆功能的代碼,flask提供了通用設施的功能,即請求鉤子。 例子: 運行結果: 第1 ...
  • 概念 公平鎖/非公平鎖 公平鎖是指多個線程按照申請鎖的順序來獲取鎖。 非公平鎖是指多個線程獲取鎖的順序並不是按照申請鎖的順序,有可能後申請的線程比先申請的線程優先獲取鎖。有可能,會造成優先順序反轉或者饑餓現象。 對於 Java ReentrantLock而言,通過構造函數指定該鎖是否是公平鎖,預設是非 ...
  • "概要" "schema" "changeset" "struct" "map" "總結" 概要 Ecto 中, 對資料庫的操作中經常用到 4 個類型: schema changeset struct map 在 Ecto 的 API 中, 基本都是都是對這 4 個類型的操作, 這 4 個類型的關係 ...
一周排行
    -Advertisement-
    Play Games
  • 前言 本文介紹一款使用 C# 與 WPF 開發的音頻播放器,其界面簡潔大方,操作體驗流暢。該播放器支持多種音頻格式(如 MP4、WMA、OGG、FLAC 等),並具備標記、實時歌詞顯示等功能。 另外,還支持換膚及多語言(中英文)切換。核心音頻處理採用 FFmpeg 組件,獲得了廣泛認可,目前 Git ...
  • OAuth2.0授權驗證-gitee授權碼模式 本文主要介紹如何筆者自己是如何使用gitee提供的OAuth2.0協議完成授權驗證並登錄到自己的系統,完整模式如圖 1、創建應用 打開gitee個人中心->第三方應用->創建應用 創建應用後在我的應用界面,查看已創建應用的Client ID和Clien ...
  • 解決了這個問題:《winForm下,fastReport.net 從.net framework 升級到.net5遇到的錯誤“Operation is not supported on this platform.”》 本文內容轉載自:https://www.fcnsoft.com/Home/Sho ...
  • 國內文章 WPF 從裸 Win 32 的 WM_Pointer 消息獲取觸摸點繪製筆跡 https://www.cnblogs.com/lindexi/p/18390983 本文將告訴大家如何在 WPF 裡面,接收裸 Win 32 的 WM_Pointer 消息,從消息裡面獲取觸摸點信息,使用觸摸點 ...
  • 前言 給大家推薦一個專為新零售快消行業打造了一套高效的進銷存管理系統。 系統不僅具備強大的庫存管理功能,還集成了高性能的輕量級 POS 解決方案,確保頁面載入速度極快,提供良好的用戶體驗。 項目介紹 Dorisoy.POS 是一款基於 .NET 7 和 Angular 4 開發的新零售快消進銷存管理 ...
  • ABP CLI常用的代碼分享 一、確保環境配置正確 安裝.NET CLI: ABP CLI是基於.NET Core或.NET 5/6/7等更高版本構建的,因此首先需要在你的開發環境中安裝.NET CLI。這可以通過訪問Microsoft官網下載並安裝相應版本的.NET SDK來實現。 安裝ABP ...
  • 問題 問題是這樣的:第三方的webapi,需要先調用登陸介面獲取Cookie,訪問其它介面時攜帶Cookie信息。 但使用HttpClient類調用登陸介面,返回的Headers中沒有找到Cookie信息。 分析 首先,使用Postman測試該登陸介面,正常返回Cookie信息,說明是HttpCli ...
  • 國內文章 關於.NET在中國為什麼工資低的分析 https://www.cnblogs.com/thinkingmore/p/18406244 .NET在中國開發者的薪資偏低,主要因市場需求、技術棧選擇和企業文化等因素所致。歷史上,.NET曾因微軟的閉源策略發展受限,儘管後來推出了跨平臺的.NET ...
  • 在WPF開發應用中,動畫不僅可以引起用戶的註意與興趣,而且還使軟體更加便於使用。前面幾篇文章講解了畫筆(Brush),形狀(Shape),幾何圖形(Geometry),變換(Transform)等相關內容,今天繼續講解動畫相關內容和知識點,僅供學習分享使用,如有不足之處,還請指正。 ...
  • 什麼是委托? 委托可以說是把一個方法代入另一個方法執行,相當於指向函數的指針;事件就相當於保存委托的數組; 1.實例化委托的方式: 方式1:通過new創建實例: public delegate void ShowDelegate(); 或者 public delegate string ShowDe ...