ขณะนี้ฉันกำลังทำงานกับการประมาณมวลชีวภาพโดยใช้ภาพถ่ายดาวเทียม ฉันจะกำหนดพื้นหลังของคำถามของฉันอย่างรวดเร็วแล้วอธิบายคำถามเชิงสถิติที่ฉันกำลังทำอยู่
พื้นหลัง
ปัญหา
ฉันพยายามประเมินมวลชีวภาพของพื้นที่ในฝรั่งเศส คำตอบของฉันคือความหนาแน่นปริมาตรไม้ไอน้ำ (เป็น ) ซึ่งมากหรือน้อยตามสัดส่วนของชีวมวล (ขึ้นอยู่กับความหนาแน่นของไม้ ... )
ตัวแปรอิสระที่ฉันมีคือดัชนีพืชพรรณที่ได้จากการสะท้อนกลับที่วัดได้ในพื้นที่นี้ (ดาวเทียมที่ใช้ในการศึกษาคือ MODIS สำหรับผู้ที่รู้ว่ามัน) ดัชนีเหล่านี้เป็นตัวอย่าง NDVI, EVI ฯลฯ ฉันมีแผนที่ของดัชนีและความละเอียดของแผนที่คือ 250m
มีความสัมพันธ์ที่แข็งแกร่งระหว่างดัชนีเหล่านี้กับปริมาณในป่าชนิดเดียวกัน (ชีวนิเวศและภูมิอากาศ) ดังนั้นฉันจึงพยายามที่จะลดความหนาแน่นของปริมาณเทียบกับตัวบ่งชี้เหล่านี้ (อันที่จริงเวลาของพวกเขา) ในแปลงสินค้าคงคลังที่ฉันรู้ปริมาณ
สินค้าคงเหลือป่าไม้
ปริมาณในแปลงเหล่านี้ประมาณด้วยวิธีการสุ่มตัวอย่างต่อไปนี้:
- โหนดสินค้าคงคลังจะถูกวางในตารางปกติที่ครอบคลุมพื้นที่
- พล็อตเชื่อมต่อกับแต่ละโหนดและกระบวนการสินค้าคงคลัง (ชนิดต้นไม้ปริมาณความสูงของหลังคาเป็นต้น) เกิดขึ้นในพล็อตนี้ แน่นอนฉันสนใจเฉพาะพล็อตสินค้าคงคลังและค่าดัชนีพืชพรรณของฉันคือค่าของพิกเซลที่มีพล็อต
กระบวนการสินค้าคงคลังในพล็อตมีดังต่อไปนี้:
- วัดต้นไม้ที่มีเส้นผ่านศูนย์กลาง> 37.5 ซม. ในวงกลมรัศมี 15 ม
- วัดต้นไม้ที่มีเส้นผ่านศูนย์กลาง> 22.5 ซม. ในวงกลมรัศมี 9 ม
- วัดต้นไม้ที่มีเส้นผ่านศูนย์กลาง> 7.5 ซม. ในวงกลมรัศมี 6 ม
จากนั้นคำนวณความหนาแน่นของปริมาตรโดยใช้ปัจจัยการขยาย
สำหรับพล็อตแต่ละครั้งฉันสามารถเข้าถึงข้อมูลสำหรับต้นไม้ที่วัดได้ทั้งหมด
ยิ่งกว่านั้นสำหรับต้นไม้แต่ละต้นฉันมีความไม่แน่นอนเกี่ยวกับปริมาณเนื่องจากการใช้สมการอัลโลเมตริก (ให้เราบอกว่า 10%)
สถิติมีความสำคัญ ...
เพื่อให้การถดถอยของฉันมีความแม่นยำยิ่งขึ้นฉันต้องการความแปรปรวน / CI ของการวัดนี้แต่ละครั้ง ขึ้นอยู่กับจำนวนต้นไม้ที่สุ่มตัวอย่างและปริมาณความหนาแน่นที่พบ
ดังนั้นฉันมีสองปัญหา:
วิธีการบัญชีสำหรับข้อเท็จจริงที่ว่าดัชนีพืชของฉันมีการวัดมากกว่าพิกเซล 250m?
ฉันสามารถสันนิษฐานได้ว่าความหนาแน่นของปริมาณคงที่มากกว่าหนึ่งพิกเซลและฉันลองพิกเซลนี้ด้วยพล็อตโฆษณาหนึ่งรายการ
จะประเมินความแปรปรวนของความหนาแน่นของโวลุ่มได้อย่างไร
ฉันคิดว่าฉันสามารถใช้ bootstrapping กับประชากรต้นไม้ แต่จำนวนต้นไม้ที่วัดได้ทั้งหมดของฉันนั้นค่อนข้างเล็ก (จาก 7 ถึง 20 ... ) ยิ่งกว่านั้นฉันจะคำนึงถึงความจริงที่ว่าฉันกำลังวัดต้นไม้ในวงกลมต่าง ๆ ขึ้นอยู่กับขนาดของมันอย่างไร และความแปรปรวนควรเปลี่ยนไปอย่างไรถ้าฉันดูพิกเซลทั้งหมด
ฉันยังคิดว่าฉันสามารถใช้การจำลองแบบมอนติคาร์โลเพื่อจำลองป่าแล้วสุ่มตัวอย่างป่านี้ด้วยแผนการเพื่อดูว่าเกิดอะไรขึ้น ...
ฉันไม่มีภูมิหลังทางสถิติที่แข็งแกร่งดังนั้นฉันจึงแพ้บ้างเล็กน้อย!