پایش و برآورد سطح زیر کشت محصولات باغی و زراعی با تلفیق سری‌های زمانی سنجنده سنتینل-۲ و الگوریتم‌های یادگیری ماشین (مطالعه موردی: منطقه مجن)

نوع مقاله : مقاله کامل علمی پژوهشی

نویسندگان

1 دانشجوی دکتری گروه مهندسی آب، دانشگاه علوم کشاورزی و منابع طبیعی گرگان

2 دانشیار گروه مهندسی اب، دانشکده مهندسی اب و خاک، دانشگاه علوم کشاورزی و منابع طبیعی گرگان، ایران.

3 دانشکده مهندسی آب و خاک، دانشگاه علوم کشاورزی و منابع طبیعی گرگان، گرگان. ایران.

4 گروه مهندسی آب، دانشکده مهندسی آب و خاک، دانشگاه علوم کشاورزی و منابع طبیعی گرگان

5 دکتری علوم اقتصادی، کارشناس مطالعات اقتصادی آب و آبفا، عضو دبیر خانه هیئت ساماندهی و راهبری بازارهای محلی آب وزارت نیرو، تهران

چکیده

سابقه و هدف: پایش دقیق الگو و سطح زیر کشت محصولات کشاورزی، به‌ویژه در مناطق کوهپایه‌ای با ساختار خرده‌مالکی و پراکنده، برای مدیریت بهینه منابع آب و برنامه‌ریزی کشاورزی امری ضروری است. با این وجود، فقدان آمار دقیق و به‌روز، برنامه‌ریزی در این مناطق را با چالش مواجه ساخته است. فناوری سنجش از دور با قابلیت پوشش گسترده و تهیه داده‌های به‌هنگام، راهکاری کارآمد برای پایش سطح زیر کشت به‌شمار می‌رود. این پژوهش با هدف پایش و برآورد سطح زیر کشت محصولات باغی و زراعی غالب در منطقه مجن با تلفیق سری‌های زمانی تصاویر سنتینل-۲ و الگوریتم‌های یادگیری ماشین انجام شد.
مواد و روش‌ها: منطقه مورد مطالعه این پژوهش، شهر مجن، با دارا بودن شرایط خاص توپوگرافی واقع در دامنه جنوبی رشته‌کوه البرز، الگوی کشت پیچیده‌ای متشکل از باغات متراکم (مانند زردآلو، سیب و گیلاس) و اراضی زراعی (عمدتاً گندم) را در خود جای داده است. عدم وجود آمار دقیق و به‌روز از سطح زیر کشت هر یک از این محصولات، برنامه‌ریزی برای تخصیص بهینه منابع آب محدود این منطقه را با مشکل مواجه ساخته است. از این‌رو، دسترسی به یک روش خودکار، دقیق و کم‌هزینه مبتنی بر فناوری سنجش از دور برای پایش الگوی کشت در این منطقه، نه تنها یک نیاز علمی، بلکه یک ضرورت عملیاتی برای مدیریت پایدار منابع آب و کشاورزی است. این پژوهش با تمرکز بر این چالش منطقه‌ای، در پی پر کردن خلأ اطلاعاتی موجود است.در این پژوهش با استفاده از تصاویر چند زمانه ماهواره سنتینل 2، اقدام به شناسایی سطح زیر کشت اراضی موجود در بازار آب مجن با بهره‌گیری از دوره فنولوژیکی گردیده است. در ابتدا اقدام به شناسایی نوع محصولات منطقه شامل زردآلو، سیب، گیلاس و گندم و نیز دوره فنولوژیکی آن‌ها گردید. سپس به عنوان اهداف تعلیمی از میان تمامی باغات و مزارع منطقه 91 باغ سیب، 254 باغ زردآلو، 15 باغ گیلاس و 35 مزرعه گندم و نیز 50 محدوده اراضی بایر و فاقد کشت به عنوان جلوگیری از تداخل طیفی انتخاب شدند. همچنین از شاخص‌های گیاهی تفاضلی نرمال شده (NDVI)، پوشش گیاهی بهبود یافته (EVI)، تفاضلی نرمال شده آب (NDWI) به عنوان کاربردی‌ترین شاخص‌های گیاهی جهت شناسایی پوشش‌های گیاهی، نوع کشت و وضعیت سبزینگی منطقه در طی دوره‌های مختلف نیز استفاده گردیده شد. در نهایت به کمک برداشت‌های میدانی و همچنین مدل‌های نظارت شده ماشین بردار پشتیبان (SVM)، شبکه عصبی (Neural Net)، حداقل فاصله از میانگین (MD) و حداکثر تشابه (MLE)، تصاویر ماهواره‌ای طبقه‌بندی و نوع محصولات و نقشه اراضی کشاورزی بدست آمد.
یافته‌ها: جهت بررسی صحت نتایج، نقشه‌های تولید شده با برداشت‌های میدانی تطبیق داده شد. نتایج بررسی نمودارهای شاخص-های گیاهی به وضوح بیانگر تفکیک نمودارها از یکدیگر و مشاهده هر یک از باغات سیب، زردآلو، گیلاس و مزارع گندم بصورت جداگانه با شدت و رنگ متفاوت بودند، لذا روند تغییرات ساختار محصولات غالب منطقه که شامل شاخص سطح برگ، نوع و ساختار تاج پوششی گیاه، مقادیر رطوبت و شناسایی پهنه‌های آبی از گیاهی می‌باشند را در شناسایی، تفکیک و همچنین آغاز روند رشد تا برداشت و خزان به خوبی نمایش داده‌اند. به جهت ارزیابی دقت مدل‌های طبقه‌بندی از روش اعتبارسنجی متقابل استفاده گردید که شاخص‌های آماری ضریب کاپای و دقت کلی در آن‌ها به کار گرفته شد، همچنین مقادیر دو خطای اومیسیون که بیانگر پیکسل‌هایی است که در مساحت اراضی به اشتباه طبقه بندی شده و به آن کلاس تعلق ندارد و نیز خطای کومیسیون که نشان دهنده اراضی یک کلاس که در کلاسهای دیگر واقع شده است به ترتیب برای محاسبه دقت تولیدکننده و دقت کاربر محاسبه گردیدند. برای روش‌های شبکه عصبی، ماشین بردار پشتیبان، حداکثر تشابه و حداقل فاصله از میانگین مقادیر ضریب کاپای به ترتیب 79/0، 78/0، 66/0 و 60/0 و دقت کلی 89، 88، 81 و 79 درصد بدست آمد که بررسی این اعداد به همراه مقادیر دوخطای کمیسیون و اومیسیون بیانگر عملکرد بهتر روش شبکه عصبی و ماشین بردار پشتیبان نسبت به دو روش دیگر بودند. نتایج نشان داد با استفاده از روش شبکه عصبی از مجموع 2000 هکتار اراضی مجن، مساحت تحت پوشش باغات زردآلو، سیب، گیلاس به ترتیب 1442، 333 و 52 هکتار (1/72، 7/16 و 6/2 درصد) و زراعت گندم 173 هکتار (6/8 درصد) می‌باشد، که در مقایسه با آمار سازمان جهاد کشاورزی از دقت 94 درصدی برخودار است. بررسی‌های میدانی نشان داد عواملی همچون کشت درختان هلو، شلیل، آلبالو و گردو به شکلی محدود در برخی از باغات منطقه به صورت ترکیبی با زردآلو و سیب با توجه به وسعت منطقه و تشابه دوره تقویم زراعی موجب بروز خطا در برآورد سطح زیر کشت این محصولات شده است، همچنین از بین نبردن علف‌های هزر در تعداد معدودی از باغات منطقه در ابتدای فصل بهار موجب بروز خطا در تفکیک این باغات با زمین‌های زراعی گندم شده. کمترین میزان خطا را می‌توان در سطح زیر کشت گندم در منطقه مشاهده نمود زیرا کشت آن بصورت پاییزه بوده و در این دوره محصولی همزمان در منطقه کشت نمی‌گردد.
نتیجه گیری: در نهایت می‌توان بیان داشت کاربرد تصاویر ماهواره‌ای چند زمانه و شاخص فنولوژیکی می‌توانند نتایج نسبتاً مقبولی را در شناسایی نوع محصولات کشاورزی و تعیین سطح زیر کشت آن‌ها ارائه نمایند که از دیدگاه‌های مدیریتی و نظارتی با توجه به مزایای آن از جمله قابلیت دسترسی به نقاط دوردست، سطح گسترده پوشش، دسترسی آسان و سریع به تصاویر و نیز دقت مناسب به جهت تهیه نقشه اراضی کشاورزی حائز اهمیت می‌باشد. در واقع می‌توان بیان داشت با بهره‌گیری از این فناوری در تمامی ابعاد مطالعاتی و مدیریتی می‌توان اطلاعات بسیار مفیدی را پیرامون مدیریت آبیاری و نحوه استفاده آن در مزارع و باغات و نیز کنترل تنش-های آبی به جهت برنامه‌ریزی در بخش‌های مختلف بهره جست و دقت کار را با صرف کمترین هزینه و نظارت دائمی افزایش داد.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Monitoring and Estimation of Horticultural and Agricultural Crop Cultivation Area through Integration of Sentinel-2 Sensor Time Series and Machine Learning Algorithms (Case Study: Mojen Region)

نویسندگان [English]

  • mojtaba shaker 1
  • MOUSA HESAM 2
  • Khalil Ghorbani 3
  • Abotaleb Hezarjaribi 4
  • Mohamad Oshani 5
1 Ph.D Student of Water Engineering Sciences, Dept. of Water Engineering, Gorgan University of Agricultural Sciences and Natural Resources
2 Associate Professor, Department of Water Engineering, Faculty of Water and Soil Engineering, Gorgan University of Agricultural Sciences and Natural Resources, Iran.
3 Professor, Dept. of water Engineering, Faculty of Water and Soil Engineering, Gorgan University of Agricultural Sciences and Natural Resources.Gorgan, Iran.
4 Dept. of water Engineering, Faculty of Water and Soil Engineering, Gorgan University of Agricultural Sciences and Natural Resources, Gorgan, Iran.
5 Associated Prof economic sciences, expert in water and water economic studies, member of the secretary of the local water markets organizing and management board of the Ministry of Energy, Tehran.
چکیده [English]

Abstract
Background and objectives: Accurate monitoring of cropping patterns and cultivated area, particularly in mountainous regions with fragmented and small-scale landholding structures, is essential for optimal water resource management and agricultural planning. However, the lack of precise and up-to-date statistics poses significant challenges for planning in these areas. Remote sensing technology, with its extensive coverage capacity and ability to provide timely data, offers an effective solution for monitoring cultivated area. This research was conducted with the objective of monitoring and estimating the cultivated area of dominant horticultural and agricultural crops in the Mojen region by integrating Sentinel-2 satellite image time series and machine learning algorithms.
Materials and methods: The study area of this research is the city of Mojen. Located on the southern slopes of the Alborz mountain range with specific topographic conditions, it features a complex cropping pattern consisting of dense orchards (such as apricot, apple, and cherry) and agricultural lands (primarily wheat). The lack of accurate and up-to-date data on the cultivation area of each of these crops has hindered planning for the optimal allocation of the region's limited water resources. Consequently, access to an automated, accurate, and low-cost method based on remote sensing technology for monitoring the cropping pattern in this region is not only a scientific need but also an operational necessity for sustainable water resource and agricultural management. This research, focusing on this regional challenge, sought to fill the existing information gap. In this study, multi-temporal Sentinel-2 satellite images were used to identify the cultivated area of lands in the Mojen water market by utilizing the phenological cycle. Initially, the types of crops in the region, including apricot, apple, cherry, and wheat, along with their phenological stages, were identified. Subsequently, from all orchards and farms in the region, 91 apple orchards, 254 apricot orchards, 15 cherry orchards, 35 wheat fields, and also 50 areas of barren and uncultivated land were selected as training samples to prevent spectral interference. Furthermore, the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Normalized Difference Water Index (NDWI) were used as the most practical vegetation indices for identifying vegetation covers, crop type, and the greening status of the region during different periods. Finally, with the help of field surveys and supervised models including Support Vector Machine (SVM), Neural Network (Neural Net), Minimum Distance (MD), and Maximum Likelihood Estimation (MLE), the satellite images were classified, and the crop types and agricultural land map were obtained.
Results: To verify the accuracy of the results, the generated maps were cross-referenced with field surveys. The analysis of vegetation index profiles clearly demonstrated the distinct separation of the profiles, with each of the apple, apricot, and cherry orchards, as well as the wheat fields, being discernible as separate entities with different intensities and colors. This effectively illustrated the phenological dynamics of the region's dominant crops—including Leaf Area Index (LAI), canopy type and structure, moisture content, and the identification of vegetated water bodies—throughout their growth cycle from onset to harvest and dormancy. A cross-validation method was employed to assess the accuracy of the classification models, utilizing the statistical metrics of Kappa coefficient and overall accuracy. Additionally, the values of the two errors—omission error (indicating pixels of a class that were incorrectly excluded from that class) and commission error (indicating pixels incorrectly included in a class)—were calculated to determine producer's accuracy and user's accuracy, respectively. For the Neural Network, Support Vector Machine (SVM), Maximum Likelihood, and Minimum Distance methods, the Kappa coefficients were 0.79, 0.78, 0.66, and 0.60, respectively, with corresponding overall accuracies of 89%, 88%, 81%, and 79%. Examination of these values, along with the commission and omission errors, indicated the superior performance of the Neural Network and Support Vector Machine methods compared to the other two. The results obtained using the Neural Network method showed that from the total 2000 hectares of land in Mojen, the areas covered by apricot, apple, and cherry orchards were 1442, 333, and 52 hectares (72.1, 16.7 & 2.6 percent), respectively, and the area under wheat cultivation was 173 hectares (8.6 percent). This represents a 94% accuracy compared to the statistics from the Agricultural Jihad Organization. Field investigations revealed that factors such as the limited intercropping of peach, nectarine, sour cherry, and walnut trees within some apricot and apple orchards—coupled with the region's vastness and similarities in agricultural calendar phases—contributed to estimation errors for these crops. Furthermore, the presence of unremoved weeds in a small number of orchards in early spring led to errors in distinguishing these orchards from wheat farmlands. The lowest error rate was observed in the estimated area under wheat cultivation, attributable to its autumn sowing season, during which no other crops are simultaneously cultivated in the region.
Conclusion: In summary, the integration of multi-temporal satellite imagery and phenological indices offers a reliable approach for identifying agricultural crop types and estimating cultivated areas with acceptable accuracy. From a management and monitoring perspective, this methodology is highly valuable due to its distinct advantages, including access to remote regions, extensive spatial coverage, rapid and easy data acquisition, and sufficient precision for generating accurate agricultural land-use maps. Furthermore, employing this technology across various research and operational contexts can provide critical insights for improving irrigation management, optimizing water use in farms and orchards, and mitigating water stress. By adopting this approach, stakeholders can enhance the accuracy of agricultural planning and monitoring across different sectors, while minimizing costs and enabling continuous, large-scale oversight.

کلیدواژه‌ها [English]

  • Crop Area Monitoring
  • Sentinel-2
  • Vegetation Indices
  • Machine Learning Algorithms