چرا NDVI لندست با اعداد خام جواب متفاوتی می‌دهد؟

برای یک پیکسل، یک محاسبه به NDVI حدود ۰٫۲۸۶ می‌رسد و محاسبهٔ دیگری به ۰٫۵۹۵. لازم نیست یکی از تصویرها عوض شده باشد؛ کافی است در محاسبهٔ اول، اعداد ذخیره‌شدهٔ باندها به‌جای بازتاب سطح استفاده شده باشند.

در داده‌های بازتاب سطح Level-2 از Collection 2 لندست، پیش از محاسبهٔ شاخص باید وضعیت اعمال ضریب و آفست روشن باشد. اول بپرسید فایل چه چیزی دارد: باندهای ذخیره‌شده، بازتابی که قبلاً تبدیل شده، یا NDVI آماده؟ هرکدام مسیر متفاوتی دارد و اعمال دوبارهٔ تبدیل هم خطاست.

مثال عددی را از نوع محصول شروع کنیم

این راهنما دربارهٔ باندهای بازتاب سطح Landsat 8–9 Collection 2 Level-2 است. بر اساس راهنمای ضریب مقیاس USGS، عدد ذخیره‌شدهٔ این باندها در ۰٫۰۰۰۰۲۷۵ ضرب می‌شود و سپس ۰٫۲ از آن کم می‌شود. مقدار جمع‌شوندهٔ منفی ۰٫۲ همان آفست است.

برای لندست ۸ و ۹، باند ۵ فروسرخ نزدیک و باند ۴ قرمز است. فرمول NDVI در مستندات USGS، اختلاف بازتاب فروسرخ نزدیک و قرمز را بر مجموع آن‌ها تقسیم می‌کند. عکس رنگیِ نقشه یا اسکرین‌شات، جای این داده‌های عددی را نمی‌گیرد.

قبل از اجرای فرمول، توضیحات محصول و خروجی سرویس را بخوانید. ممکن است سامانه‌ای که فایل را از آن گرفته‌اید، ضریب و آفست را قبلاً اعمال کرده باشد. در این حالت باید از بازتاب حاصل استفاده کنید، نه اینکه تبدیل را بار دیگر انجام دهید.

محاسبهٔ یک پیکسل ساختگی، قدم‌به‌قدم

اعداد این بخش مثال ساختگی برای تمرین محاسبه هستند؛ از یک مزرعه یا تصویر واقعی برداشت نشده‌اند. محاسبه در ۲ نوامبر ۲۰۲۶ با حساب ده‌دهی روی رایانه انجام شد. عدد ذخیره‌شدهٔ باند قرمز را ۱۰٬۰۰۰ و فروسرخ نزدیک را ۱۸٬۰۰۰ در نظر بگیرید.

حالا با بازتاب تبدیل‌شده حساب می‌کنیم:

(0.295 - 0.075) / (0.295 + 0.075) = 0.594594...

ولی اگر همان اعداد ذخیره‌شده را مستقیم در نسبت بگذاریم:

(18000 - 10000) / (18000 + 10000) = 0.285714...

جواب اول تقریباً ۰٫۵۹۵ و جواب دوم تقریباً ۰٫۲۸۶ است. هر دو ممکن است در یک نقشهٔ رنگی عادی به نظر برسند. پس معقول‌بودن ظاهر عدد، صحت مسیر محاسبه را نشان نمی‌دهد.

چرا ضریب حذف می‌شود، اما آفست نه؟

اگر تبدیل فقط ضرب هر دو باند در یک ضریب مشترک بود، آن ضریب از صورت و مخرج کسر حذف می‌شد. اما مقدار جمع‌شونده در تفاضل دو باند حذف می‌شود و در مجموع آن‌ها دو بار باقی می‌ماند. اگر اعداد ذخیره‌شده را با N و R، ضریب را با a و آفست را با b بنویسیم، داریم:

NDVI = a(N - R) / [a(N + R) + 2b]

در محصول مورد بحث، b صفر نیست. نادیده‌گرفتن همین جزء، جواب تمرین را تغییر داد. این دلیل برای همین نوع ورودی معتبر است؛ نباید ضریب لندست را به هر مجموعه‌دادهٔ پوشش گیاهی تعمیم داد.

اگر فایل از قبل NDVI باشد چه؟

محصول آمادهٔ NDVI مسیر جداگانه‌ای دارد. در مشخصات محصول NDVI سازمان USGS ضریب ۰٫۰۰۰۱ آمده است. بنابراین یک مقدار معتبر ذخیره‌شدهٔ ۶٬۰۰۰، شاخص ۰٫۶ را نشان می‌دهد. آفست منفی ۰٫۲ باندهای بازتاب سطح را نباید به این شاخص آماده اعمال کرد.

نوع محصول، ضریب، مقدارهای مخصوص دادهٔ خالی و اشباع را از مستندات همان فایل تأیید کنید؛ اسم فایل به‌تنهایی کافی نیست. اعداد اولیه و مراحل تبدیل را نگه دارید تا مشخص باشد شاخص را از باندها ساخته‌اید یا یک شاخص آماده را خوانده‌اید.

دادهٔ نامعتبر را با صفر پر نکنید

درست‌بودن تبدیل عددی، کیفیت پیکسل را تضمین نمی‌کند. راهنمای باندهای کنترل کیفیت USGS پرچم‌های دادهٔ خالی، ابر، سایهٔ ابر و برف را در QA_PIXEL و اشباع رادیومتری را در QA_RADSAT توضیح می‌دهد. قواعد کنارگذاشتن پیکسل‌ها را متناسب با تحلیل انتخاب و ثبت کنید. عبارت «بدون ابر» جای شرح این قواعد را نمی‌گیرد.

تمرین دوم: در سه موقعیت، مقدارها ۰٫۲، «فاقد داده» و ۰٫۶ هستند. میانگین دو مقدار معتبر ۰٫۴ است. اگر موقعیت خالی را صفر بگذارید، میانگین به حدود ۰٫۲۶۷ می‌رسد و معنای نتیجه تغییر می‌کند. کنار میانگین بنویسید «دو پیکسل معتبر از سه پیکسل». همچنین اگر مجموع بازتاب دو باند صفر شود، کسر NDVI تعریف نشده است و باید برای این حالت هم قاعدهٔ صریح دادهٔ نامعتبر داشته باشید.

از عدد صحیح تا ادعای قابل دفاع

NDVI برابر ۰٫۶ به معنی «۶۰ درصد پوشش گیاهی» یا «۶۰ درصد سلامت» نیست. این عدد به‌تنهایی نوع محصول کشاورزی یا سلامت گیاه را مشخص نمی‌کند. ترکیب اجزای داخل یک پیکسل، تاریخ مشاهده و داده‌های کنارگذاشته‌شده بر مقایسه اثر دارند. برای اینکه ببینید ویژگی موردنظر اصلاً در مقیاس تصویر جداشدنی است یا نه، راهنمای تفکیک‌پذیری تصاویر ماهواره‌ای را ببینید.

هنگام بررسی یک ادعای سبزشدن یا خشک‌شدن، شناسه و تاریخ محصول، نام باندها، وضعیت اعمال ضریب و آفست، قواعد کنترل کیفیت و تعداد پیکسل‌های معتبر را درخواست کنید. این موارد را همراه فرمول در گزارش شواهد بیاورید. بدون این اطلاعات، محاسبه هنوز قابل تأیید نیست؛ حتی اگر نقشه ظاهر قانع‌کننده‌ای داشته باشد.