คำถามในหัวข้อ 'scipy'

วิธีลบคอลัมน์ใน numpy.array
ฉันต้องการลบคอลัมน์ที่เลือกใน numpy.array นี่คือสิ่งที่ฉันทำ: n [397]: a = array([[ NaN, 2., 3., NaN], .....: [ 1., 2., 3., 9]]) In [398]: print a [[ NaN 2. 3. NaN] [ 1. 2. 3. 9.]] In [399]: z = any(isnan(a),...
217586 มุมมอง
schedule 21.12.2023

การดำเนินการแบบ pointwise บนเมทริกซ์ scipy.sparse
เป็นไปได้ไหมที่จะใช้ตัวอย่าง numpy.exp หรือตัวดำเนินการ pointwise ที่คล้ายกันกับองค์ประกอบทั้งหมดใน scipy.sparse.lil_matrix หรือรูปแบบเมทริกซ์กระจัดกระจายอื่น import numpy from scipy.sparse import lil_matrix x = numpy.ones((10,10)) y =...
1260 มุมมอง
schedule 07.12.2023

ปรับแต่งเครื่องหมายขีดเพียงด้านเดียวใน matplotlib โดยใช้สัน
ฉันมีแถบแนวนอนของ matplotlib ที่วาดดังนี้: import matplotlib.pyplot as plt from numpy import * from scipy import * bars = arange(5) + 0.1 vals = rand(5) print bars, vals plt.figure(figsize=(5,5), dpi=100) spines = ["bottom"] ax = plt.subplot(1, 1,...
11451 มุมมอง
schedule 30.11.2023

สร้างฟังก์ชันสำหรับการคำนวณความแปรปรวนร่วม
มีวิธีใดในหลามในการรับเมทริกซ์ความแปรปรวนร่วมโดยพิจารณาจากค่าเฉลี่ยและจุดข้อมูลตัวอย่าง ตัวอย่าง: mean = [3 3.6] data = [[1 2] [2 3] [3 3] [4 5] [5 5]]...
9615 มุมมอง
schedule 25.12.2023

การลดอาร์เรย์จำนวนมากสำหรับการวาดแผนภูมิ
ฉันต้องการวาดแผนภูมิในแอปพลิเคชันหลามของฉัน แต่อาร์เรย์ numpy ของแหล่งที่มามีขนาดใหญ่เกินไปสำหรับการทำเช่นนี้ (ประมาณ 1'000'000+) ฉันต้องการหาค่าเฉลี่ยขององค์ประกอบข้างเคียง แนวคิดแรกคือการทำในรูปแบบ C++: step = 19000 # every 19 seconds (for...
784 มุมมอง
schedule 24.11.2023

เหตุใด ndimage.label(image) จึงติดป้ายกำกับรูปร่างพิเศษในกลุ่มที่พิจารณา
เมื่อฉันใช้ ndimage.label(img) นำเข้าจากแพ็คเกจ scipy เพื่อติดป้ายกำกับรูปภาพ PNG ระดับสีเทา มันจะทำงานเช่นนี้ ฉันมีสองภาพที่มีรูปร่างบางอย่างที่สร้างโดย Photoshop : รูปภาพแรก : http://imageshack.us/a/img140/8669/onehx.png...
5681 มุมมอง
schedule 06.11.2023

การวนซ้ำแบบ scipy kmeans หมายความว่าอย่างไร
ฉันใช้อัลกอริธึม kmeans2 ตั้งแต่สีพิกเซลแบบ scipy ไปจนถึงสีคลัสเตอร์ในรูปภาพ เพื่อให้ได้สีโดยเฉลี่ยอันดับต้นๆ ของรูปภาพ http://docs.scipy.org/doc/scipy/reference/generated/scipy.cluster.vq.kmeans2.html#scipy.cluster.vq.kmeans2...
3411 มุมมอง
schedule 01.01.2024

การค้นหาจุดสูงสุดและเซนทรอยด์อย่างรวดเร็วใน python
ฉันกำลังพยายามพัฒนาอัลกอริธึมที่รวดเร็วใน python เพื่อค้นหาพีคในรูปภาพ จากนั้นค้นหาเซนทรอยด์ของพีคเหล่านั้น ฉันได้เขียนโค้ดต่อไปนี้โดยใช้ scipy.ndimage.label และ ndimage.find_objects เพื่อค้นหาวัตถุ นี่ดูเหมือนจะเป็นจุดคอขวดในโค้ด และใช้เวลาประมาณ 7...
12631 มุมมอง

การซิงโครไนซ์ข้อมูลใน Python (หรือ Excel)
ฉันใช้ Python บ่อยครั้ง (และ Excel เป็นครั้งคราว) เพื่อประมวลผลและเปรียบเทียบข้อมูลการทดสอบระหว่างการทดลองหลายครั้ง ในบางกรณีข้อมูลอาจไม่ตรงกันซึ่งทำให้การเปรียบเทียบโดยตรงทำได้ยาก ตัวอย่างเช่น ข้อกำหนดการทดสอบทั่วไปจะเป็นดังนี้: 1) Stabilize...
1332 มุมมอง
schedule 04.01.2024

การใช้ scipy.optimize.curve_fit ภายในคลาส
ฉันมีชั้นเรียนที่อธิบายฟังก์ชันทางคณิตศาสตร์ คลาสจะต้องมีกำลังสองน้อยที่สุดที่พอดีกับตัวเองจึงจะส่งข้อมูลได้ นั่นคือคุณสามารถเรียกวิธีการดังนี้: classinstance.Fit(x,y) และจะปรับตัวแปรภายในให้เหมาะสมกับข้อมูลมากที่สุด ฉันกำลังพยายามใช้...
1967 มุมมอง

คำสั่งนำเข้า Python ในลูป: การนำเข้ารันการวนซ้ำทุกครั้งหรือไม่
สำหรับโค้ดที่ฉันกำลังเขียน ฉันกำลังรัน scipy.curve_fit() หลายหมื่นครั้ง ฉันสังเกตเห็นใน ซอร์สโค้ด curve_fit() ที่เกี่ยวข้อง โดยเฉพาะใน บรรทัด 430 และ 431 ใน แหล่งที่มา (ในฟังก์ชัน leastsq() ) มีคำสั่ง import สองคำสั่ง: from numpy.dual...
961 มุมมอง
schedule 10.01.2024

ใน scipy อะไรคือจุดสำคัญของฟังก์ชันระยะทางสองฟังก์ชันที่แตกต่างกันที่ใช้ในการจัดกลุ่มแบบลำดับชั้น?
มีฟังก์ชันระยะทางหนึ่งที่ฉันสามารถส่งผ่านไปยัง pdist ใช้เพื่อสร้างเมทริกซ์ระยะทางที่กำหนดให้กับ linkage มีฟังก์ชันระยะทางที่สองที่ฉันสามารถส่งผ่านไปยัง linkage เป็นหน่วยเมตริกได้ เหตุใดจึงมีฟังก์ชันระยะทางที่เป็นไปได้สองฟังก์ชัน...
155 มุมมอง
schedule 01.11.2023

ค้นหาดัชนีของจุดที่ใกล้ที่สุดบนตาราง 2D ที่ไม่ใช่สี่เหลี่ยมได้อย่างมีประสิทธิภาพ
ฉันมีตาราง lon/lat ที่ผิดปกติ (ไม่ใช่สี่เหลี่ยม) และมีจุดจำนวนมากในพิกัด lon/lat ซึ่งควรสอดคล้องกับจุดบนตาราง (แม้ว่ามันอาจจะผิดไปเล็กน้อยด้วยเหตุผลเชิงตัวเลข) ตอนนี้ฉันต้องการดัชนีของจุด lon/lat ที่สอดคล้องกัน ฉันได้เขียนฟังก์ชันที่ทำสิ่งนี้...
1953 มุมมอง

TypeError: function() รับ 2 อาร์กิวเมนต์ (ให้ 1 รายการ) (python)
import numpy as np import scipy.optimize as spo def function(x,y): return (np.sin(x*y+y)*np.exp(-1*(x**2+y**2)))**-1 xi=[0,0] answer=spo.fmin(function,xi) print 'the answer is', answer ฉันกำลังพยายามย่อฟังก์ชั่นนี้ให้เล็กสุด...
5053 มุมมอง
schedule 29.10.2023

ใช้ Odeint เพื่อแก้สมการเชิงอนุพันธ์อันดับหนึ่งสองอันด้วยเวกเตอร์เป็นอินพุต [ซ้ำกัน]
ฉันมีระบบสมการเชิงอนุพันธ์สองแบบ: dri/dt=tan(al) dal/dt=(vz-C_alz)/C_aln โดยที่ vz เรียกว่านิพจน์ซึ่งขึ้นอยู่กับ l , C_alz และ C_aln เท่านั้น ซึ่งเป็นนิพจน์สองรายการของตัวแปร ri และ al...
129 มุมมอง
schedule 15.01.2024

วิธีการประมาณค่า scipy interp1d fill_value = tuple ไม่ทำงาน
ฉันต้องการคาดการณ์ฟังก์ชันที่เหมาะสม scipy.interpolate.interp1d ควรจะสามารถทำได้ (ดูตัวอย่างเอกสาร) แต่ฉันได้รับ "ValueError: ค่าใน x_new อยู่ต่ำกว่าช่วงการแก้ไข" การใช้: python 2.7.12, numpy 1.13.3, scipy 0.19.1 fill_value : array-like...
5753 มุมมอง
schedule 25.10.2023

การกรองพิกัดตามระยะทางจากจุดหนึ่ง
ฉันมีสองอาร์เรย์พูดว่า: A = np.array([[ 1. , 1. , 0.5 ], [ 2. , 2. , 0.7 ], [ 3. , 4. , 1.2 ], [ 4. , 3. , 2.33], [ 1. , 2. , 0.5 ], [ 6. , 5. , 0.3 ],...
694 มุมมอง
schedule 31.12.2023

ฉันจะคำนวณค่าเฉลี่ยโดยปล่อย NaN และค่าผิดปกติจาก dataframe ในรูปแบบนี้ได้อย่างไร
ฉันมี dataframe ในรูปแบบด้านล่าง: Original Dataframe | x | value1 | value2 | value3 | value4 ---|-----|----------|----------|----------|----------- 0 | 1 | 1 | NaN | 3 | 1 1 | 2 | 4 | NaN |...
285 มุมมอง
schedule 12.01.2024

วิธีสุ่มตัวอย่างจากการแจกแจงที่กำหนด CDF ใน Python
ฉันต้องการดึงตัวอย่างจากการแจกแจงความน่าจะเป็นด้วย CDF 1 - e^(-x^2) มีวิธีใน python/scipy/etc หรือไม่ เพื่อให้คุณสามารถสุ่มตัวอย่างจากการแจกแจงความน่าจะเป็นโดยให้เฉพาะ CDF เท่านั้น
433 มุมมอง
schedule 22.12.2023

ส่งกลับระยะทางไปยังจุดที่ใกล้ที่สุดโดยกลุ่ม - python
ฉันหวังว่าจะคำนวณระยะทางไปยังจุดที่ใกล้ที่สุดโดยจัดกลุ่มตามรายการเฉพาะ โดยเฉพาะการใช้ด้านล่างนี้ calculate_distances จะวัดระยะทางระหว่าง id แต่ละรายการกับจุดที่เหลือ ฉันหวังว่าจะคืนระยะทางไปยังจุดที่ใกล้ที่สุด แต่สำหรับแต่ละรายการใน Group...
79 มุมมอง
schedule 22.11.2023