پالم روغنی (Elaeis guineensis Jacq.) در دهههای اخیر رونق تولید داشته است و منطقه تولیدی اصلی آن در داخل مناطق استوایی است که بهترین شرایط بیوفیزیکی را برآورده میکنند. سرمایه گذاران و همچنین متخصصان زمین فضایی به طور فزاینده ای به بهترین شرایط رشد و برداشت علاقه مند هستند. بنابراین هدف این مقاله انتخاب بهترین محل کاشت نخل روغنی از طریق ترکیب بهترین روش ها می باشد. منطقه مورد مطالعه ناحیه Njimom واقع در غرب کامرون، انتقالی بین مناطق استوایی و آب و هوایی است. در همان محیط GIS، ترکیب خطی وزنی (WLC) و فرآیند تحلیل سلسله مراتبی فازی (FAHP) به ترتیب تفاوتهای ظریف بین قابلیت و مناسب بودن را برجسته میکنند، در حالی که تابع سودمندی (UF) به ارزیابی جنبههای پایداری کمک میکند. اولین نتایج شامل هشت لایه است که شرایط طبیعی را نشان می دهد، یعنی بارندگی، دما، آفتاب، شیب، ارتفاع، غنای خاک، رطوبت خاک و پوشش جنگل، که در شش طبقه از 5 تا 0 بر اساس مقیاس استاندارد فائو رتبه بندی شده اند. آنها با استفاده از روش ساده WLC برای دادن “لایه قابلیت” متقاطع می شوند. نتایج دوم شامل سه لایه مربوط به محدودیتهای اقتصادی-اجتماعی برای تولید، بهعنوان منطقه ساختهشده، فاصله تا جاده و فاصله تا رودخانه است. این لایه ها به صورت باینری با 1 و 0 کدگذاری مجدد می شوند، آنها با بردار اولویت FAHP و رویکرد عضویت وزن می شوند تا “لایه مناسب” را ارائه دهند. سپس، تعداد وقوع هر یک از جنبه های پایداری در هر یک از دو فرآیند قبلی برای انجام UF شمارش می شود. مقدار حاصل، یعنی 0.
کلید واژه ها
روغن پالم , GIS , WLC , FAHP , UF , قابلیت , مناسب , پایداری
1. مقدمه
• il Palm (Elaeis guineensis Jacq.) یکی از مهم ترین گیاهان درختی استوایی جهان است. عمدتاً برای تولید صنعتی روغن نباتی [ 1 ] کشت می شود. کاربردهای مختلف شامل ساخت صابون، مواد شوینده و سایر لوازم بهداشتی در صنایع شیمیایی اولئو برای نخل روغن هسته [ 2 ] [ 3 ] و همچنین مواد تشکیل دهنده بسیاری از غذاها، در صنایع چرم، کنجاله و شیمیایی، برای نخل روغن خام است. [ 4 ] . با توجه به [ 5 ]، تخمین زده می شود 50 درصد از محصولات بسته بندی شده فروخته شده در سوپرمارکت ها حاوی برخی از روغن های موجود در همه جا هستند. اخیراً در تولید بیودیزل نقش داشته است. علاوه بر این، [ 6 ] [ 7 ] [ 8 ] [9 ] تخمین می زنند که افزایش تقاضا برای روغن های گیاهی در آینده عمدتاً توسط روغن پالم به جای سایر روغن های گیاهی برآورده می شود. این رونق تولید در دهه های گذشته را توضیح می دهد که بیش از 30 درصد از تولید جهانی روغن نباتی را نشان می دهد [ 10 ] و آن را به پرمصرف ترین روغن خوراکی در جهان تبدیل می کند [ 4 ] [ 11 ]. انتظار می رود تقاضا برای روغن پالم تا سال 2020 دو برابر شود [ 12 ].
منطقه تولیدی اصلی در داخل مناطق استوایی بین 10 درجه شمالی و 10 درجه شمالی از استوا واقع شده است، جایی که شرایط بیوفیزیکی و محیطی برای رشد آن ترکیب شده است [ 1 ] [ 11 ]. به این دلایل، 43 کشور اول مولد کشورهای استوایی هستند که در جنوب شرقی آسیا، آفریقا، آمریکای مرکزی و جنوبی قرار دارند [ 3 ] [ 13 ]. در اکثر این کشورها، دولت ها به طور فزاینده ای کشت آن را به عنوان یکی از عوامل اصلی کاهش فقر و همچنین استقلال غذا و انرژی ترویج می کنند [ 14 ].] . سپس درختان خرما به دو صورت پرورش داده می شوند. آنها در باغ های روستایی مشاهده می شوند که روغن را برای مصرف محلی با عملکرد و مقدار روغن نسبتاً کمتر تهیه می کنند. با این حال، برای مقاصد صنعتی، درختان خرما عمدتاً در مناطق یکنواخت بزرگ 3000 تا 5000 هکتار رشد و برداشت میشوند [ 1 ] [ 11 ]. سپس پرورشدهنده یا سرمایهگذار باید از روشهای مدرن برای بهترین برداشت و مزایای اقتصادی استفاده کند [ 11] . او باید همزمان دو بعد دیگر پایداری را که عبارتند از توسعه اجتماعی و حفاظت یا احیای محیط زیست، ادغام کند. تکنیکهای جغرافیایی توانایی رویارویی با این چالش را نشان میدهند و نتایج پردازش به طور گسترده توسط شرکتها، سازمانهای غیردولتی یا افراد دیگری که میخواهند در هر فعالیت کشاورزی سرمایهگذاری کنند، درخواست میشود.
در واقع، بهبود و دقت تکنیکهای مکانی که به نام سیستمهای اطلاعات جغرافیایی (GIS) نیز شناخته میشوند، امروزه در فرآیندهای تصمیمگیری در دنیای واقعی، از جمله برای شیوههای کشاورزی قابل توجه است [ 15 ]. یا برای انتخاب بهترین مکان، مزایای اقتصادی، گسترش یا نظارت فضایی، سنجش از راه دور و تجزیه و تحلیل فضایی پتانسیل قابل توجهی برای حمایت از کشت نخل روغنی مدرن دارد. [ 6 ] یک مدل مناسب جهانی برای کشت نخل روغنی را با ترکیب نقشههای هر عامل بیوفیزیکی مربوطه که آب و هوا، خاک و توپوگرافی است، محاسبه کرده است. پروتکل خاصی توسط [ 16 ] برای زمین مناسب نخل روغنی در جنوب غربی کامرون طراحی و آزمایش شد. [ 17] جنگل زدایی را پیش بینی کرده و مناسب بودن زمین را برای نخل روغنی در پاپوآ گینه نو ارزیابی کرده است. اخیراً، [ 18 ] مدلهایی از مناسب بودن کشاورزی آبی برای دستیابی به امنیت غذایی در یک محیط خشک ارائه کرده است. [ 19 ] شاخص های فازی را برای ارزیابی تناسب زمین کشاورزی ترکیب کرده است.
• پس از ارزیابی تناسب زمین و رشد محصولات، تمرکز بر پایش برای اهداف مختلف مانند کاربری زمین ، جنگل زدایی، سن، شمارش، بیماری ها و غیره است . در مالزی. به همین ترتیب، [ 21 ] گسترش نخل روغنی در اندونزی، مالزی و پاپوآ گینه نو را در بازه زمانی 2010 تا 2050 پیش بینی کرده است . مناطق رشد و ترسیم تفاوت در مراحل نخلستان. [ 28 ] و [ 29] همین کار را با استفاده از داده های LiDAR، PALSAR 50-m و WorldView-3 Imagery، برای گسترش مزارع نخل روغنی، شمارش درختان نخل و تخمین سن انجام داد. [ 30 ] با استفاده از داده های ماهواره ای نوری SPOT 7، نقشه توزیع روغن نخل کوچک را تولید کرده و مقادیر پس پراکندگی ALOS PALSAR-2 را با ارتفاع و سن نخل روغنی در خاک معدنی و تورب زمین بررسی کرده است. [ 27 ] و [ 31 ] یک همجوشی داده های چند حسگر را برای نقشه برداری از نخل آفریقا آزمایش کرده اند. [ 32 ] از سنجش از دور برای تخمین زیست توده استفاده کرده است، در حالی که [ 33 ] استراتژی های کشاورزی نخلستان ها، زیست توده و انرژی تولید شده را برای مدیریت طولانی مدت منابع نخل روغن به هم مرتبط کرده است. [ 34] بر تبدیل مناطق جنگلی به کشت نخل روغنی تمرکز کرده است، [ 35 ] این کار را در منطقه جنوب غربی کامرون نیز انجام داد، با گسترش نخل های روغنی، جنگل زدایی و آسیاب های غیررسمی ارتباط داشت و [ 36 ] نخلستان ها را ارزیابی کرد. [ 22 ] و [ 37 ] ذخایر کربن روغن نخل و اثرات زیست محیطی آن را نشان داده اند. به طور کلی، [ 38 ] اخیراً همه این کاربردهای سنجش از دور را بررسی کرد و هفت گروه بهعنوان طبقهبندی پوشش زمین، تشخیص تغییر، شمارش درختان، تخمین سن، تخمین AGB و تولید کربن، تشخیص آفات و بیماریها، و تخمین عملکرد متمایز شدند.
در این مطالعه، هدف اصلی یک مدل سازی فضایی برای انتخاب بهترین پالمسایت روغنی برای شرایط کاشت و رشد پایدار، از نظر ارزیابی مناسب/قابلیت است [ 39 ] [ 40 ]. مناسب بودن اغلب به وعده های اجتماعی-اقتصادی اشاره دارد و به چارچوب های مورد استفاده برای ارزیابی ظرفیت زمین در حمایت از تولید اولیه اشاره می کند [ 15 ] [ 41 ]. در حالی که قابلیت معمولاً پتانسیل فیزیکی یا طبیعی محیطی را به عنوان آب و هوا، امداد، خاک، هیدرولوژی و پوشش گیاهی برای حمایت از یک کاربری معین نشان میدهد [ 15 ] [ 42 ].] . این تفاوتهای ظریف در این مقاله برای ادغام بهتر سه جنبه توسعه کشاورزی پایدار، یعنی حفظ زمین، آب، منابع ژنتیکی گیاهی و جانوری، از نظر زیستمحیطی غیر تخریبکننده، از لحاظ فنی مناسب، از نظر اقتصادی بادوام و از نظر اجتماعی قابل قبول در نظر گرفته شده است [ 11 ]. .
2. روش شناسی
2.1. چارچوب کلی
طبق [ 1 ] [ 11 ] [ 43 ] [ 44]، کشت نخل روغنی به شرایط خاص طبیعی و اقتصادی نیاز دارد. برای جنبه های طبیعی، آب و هوا، خاک، قرار گرفتن در معرض نور خورشید و محیط فیزیکی عمدتا مورد توجه است. نخل روغنی در جایی که در تمام طول سال گرم است به خوبی رشد می کند و دمای آن بین 25 تا 28 درجه سانتیگراد است و هرگز کمتر از 15 درجه سانتیگراد است. درختان خرما برای فتوسنتز قوی باید دائماً در معرض آفتاب زیاد قرار گیرند، بنابراین برگها بزرگ میشوند، میوهها به خوبی میرسند و روغن بیشتری در میوهها وجود خواهد داشت. شرایط بهینه بارندگی 2000 تا 2500 میلی متر در سال با حداقل 100 میلی متر در ماه، با 1000 میلی متر و دوره خشکی 5 ماهه در سال به عنوان بدترین رطوبت قابل زندگی است. زمین باید مسطح باشد زیرا نخل روغنی در مناطق پست زیر 300 – بهترین رشد را دارد. 400 متر ارتفاع برای جلوگیری یا به حداقل رساندن خطر فرسایش. خاک باید عمیق و قابل نفوذ باشد تا بتواند رشد خوب ریشه ها را فراهم کند و به اندازه کافی غنی از نمک های معدنی باشد تا بتواند خوشه های بزرگ میوه تولید کند. تامین رطوبت برای تنظیم دما معمولاً از طریق وجود جنگل های انبوه اطراف مهم است.
علاوه بر این، جنبه های اقتصادی به نوع مزرعه ای که باید ایجاد شود بستگی دارد. [ 11 ] سه نوع مزرعه نخل روغنی را متمایز می کند. یک مزرعه روغن نخل در مقیاس کوچک ممکن است 7.5 هکتار را برای کشاورزان محلی پوشش دهد. مزرعه در مقیاس متوسط بین 10 تا 500 هکتار است. این نوع مزرعه معمولاً از شیوه های زراعی مدرن مانند فاصله گذاری بوته ها، کشت پوششی، کود دهی، وجین حلقه ای، هرس استفاده می کند. مزارع در مقیاس بزرگ مساحتی بیش از 500 هکتار را پوشش می دهند. اینها شرکت های دولتی هستند که برای تامین نیازهای مصرف داخلی کشور و تامین مازاد صادرات ایجاد شده اند. برای اندازه های متوسط و بزرگ، کارخانه ها و کارخانه ها باید نزدیک به بهره برداری باشند و همچنین زیرساخت های جاده ای برای حمل و نقل آسان. چنین گفت [ 45 ] و [ 46] برخی از چارچوب های اساسی طبقه بندی مناسب شرایط رشد را پیشنهاد کرده اند ( جدول 1 ). رتبه بندی در مقیاس پیشنهادی 5 تا 0 از شرایط بهینه تا بدترین، با ارزیابی زبانی فازی در مقیاس های میانی انجام شده است ( جدول 2 ).
2.2. روش های ریخته گری و سابقه مرتبط
نخل روغنی که در شرایط مختلف رشد می کند به اندازه کافی بیانگر یک مشکل تصمیم گیری چند معیاره (MCDM) است. عمدتاً در ترکیب اطلاعات از چندین معیار برای تشکیل یک شاخص واحد ارزیابی نقش دارد [ 47 ]. بسته به زمینه تجزیه و تحلیل GIS، می توان آن را از تجزیه و تحلیل تخصیص منابع یا تجزیه و تحلیل مناسب/قابلیت واجد شرایط دانست، همانطور که در این مقاله [ 39 ] [ 40 ]] . سپس، انتخاب اصول و روشهای حمایتی که منجر به تجزیه و تحلیل میشوند، بر اساس سه جنبه است که توسط شرایط کاشت، یعنی قابلیت طبیعی، تناسب اجتماعی-اقتصادی و پایداری کلی مطرح میشود. روش های مورد استفاده شامل ترکیب خطی وزنی (WLC)، فرآیند تحلیل سلسله مراتبی (AHP) و تنوع فازی آن و تابع مطلوبیت (UF) می باشد.
・ ترکیب خطی وزنی (WLC)
WLC یا وزندهی افزودنی ساده یکی از روشهای ارزیابی چند معیاره (MCE) به عنوان یک قواعد تصمیمگیری مبتنی بر GIS برای تجزیه و تحلیل تناسب زمین است [ 48 ] [ 49 ]. عملاً، مؤلفهها و شرایط مناسب بر روی نقشهها منتقل میشوند و مطابق با انتخاب مقیاس تحلیلگر استاندارد میشوند. برای هر یک از آنها وزن های با اهمیت نسبی در نظر گرفته شده است. سپس، هر دو وزن و نقشه های استاندارد شده با استفاده از رابطه زیر (1) ترکیب می شوند:
اس= ∑wمنایکسمنS=∑wixi(1)
که در آن S تناسب است، w i وزنی برای عامل i، و x i امتیاز معیار عامل i است.
فرآیند AHP گاهی اوقات از پیاده سازی GIS WLC پشتیبانی می کند، اگر وزن پیش فرض، 1، باید تغییر کند.
・ فرآیند تحلیل سلسله مراتبی فازی (FAHP)
AHP در دهه 70 توسط [ 50 ] به عنوان ابزار کمکی برای تصمیم گیرندگان برای یافتن “بهترین” جایگزین از بین چندین گزینه با در نظر گرفتن تعدادی از معیارهای متناقض توسعه یافت [ 49 ] [ 50 ]. این فرآیند بر اساس سه اصل استوار است، یعنی تجزیه، قضاوت مقایسه ای و ترکیب داده ها [ 50 ] [ 51 ] [ 52 ] [ 53 ]. بسته به مطالعه، مجموعه فازی را می توان با AHP [ 54 ] ادغام کرد. ایده تئوری مجموعه های فازی این است که با ابهام و عدم وجود ویژگی ذاتی در فرمول بندی ترجیحات، محدودیت ها و اهداف انسان بپردازد [ 55 ]] . اطلاعات با مقادیر زبانی یا فازی به جای اعداد بیان می شود و این رویکرد معیارها و ترکیبی از داده های کمی و کیفی را کنترل می کند [ 56 ]. محبوب ترین عدد فازی، عدد فازی مثلثی (TFN) است که با سه پارامتر به صورت l، m، u نشان داده می شود که مربوط به کوچکترین مقدار ممکن، امیدوارکننده ترین مقدار و بزرگترین مقدار ممکن است. برای یک عدد فازی نشان داده شده است آ˜A˜= (l، m، u)، تابع عضویت به صورت زیر تعریف می شود (معادله (2)، شکل 1 ):
μآ˜( x ) =⎧⎩⎨⎪⎪⎪⎪x – lm − l اگر l ≤ x ≤ m u − xu − m اگر m ≤ x ≤ u 0 در غیر این صورت μA˜(x)={x−lm−l if l≤x≤mu−xu−m if m≤x≤u0 otherwise(2)
علاوه بر این، متغیرهای زبانی که امکان ترجمه زبان طبیعی را به گزارههای منطقی یا عددی میدهند، به صورت سلسله مراتبی بهعنوان یک شبکه تصمیمگیری فازی زبانی، بر اساس ترجیحات متخصص یا رابطه علی بین مقادیر [ 57 ] سازماندهی میشوند ( جدول 3 ، شکل 2 ، جدول 4) . ، شکل 3 ).
・ تابع ابزار (UF)
هر هدف چقدر برای پایداری کشت نخل روغنی مفید است؟ پاسخ به این سوال ادغام تابع ابزار را توضیح می دهد. UF نسبت دادن منطقی مقادیر به سبدهای مختلف مصرف است، به گونه ای که مطلوب ترین سبدها مقادیر بالاتری نسبت به سبدهای کمتر دریافت می کنند [ 58 ]. اصل، رتبه بندی ترجیحات مصرف کننده است. بیان شده است که یک سبد (ایکس1،ایکس2)(x1,x2)به سبد ترجیح داده می شود (y1،y2)(y1,y2)، فقط در صورتی که سطح ابزار U از (ایکس1،ایکس2)(x1,x2)بالاتر از یکی از (y1،y2)(y1,y2). معادله (3) فرمول را نشان می دهد.
(ایکس1،ایکس2) > (y1،y2) اگر U (ایکس1،ایکس2) > U(y1،y2)(x1,x2)>(y1,y2) if U(x1,x2)>U(y1,y2)(3)
UF مورد استفاده در این مطالعه به صورت زیر نوشته شده است (معادله (4)):
U(ایکس1،ایکس2) =ایکسc / ( c + d)1ایکسد/ ( ج + د)2U(x1,x2)=x1c/(c+d)x2d/(c+d)(4)
که در آن c و d اعداد مثبت با رابطه هستند:
c > 0 ، d> 0 ، c + d= 1c>0,d>0,c+d=1
شکل 1 . تصویر یک عدد فازی مثلثی.
شکل 2 . توابع عضویت مقادیر زبانی مقیاس اولویت فازی.
شکل 3 . توابع عضویت ارزش های زبانی برای روابط علی.
از مجموع موارد فوق، این مقاله دو آرزوی خاص دارد. اولین مورد نظری است، با در نظر گرفتن مناسب بودن و قابلیت به عنوان اهداف متفاوت بر خلاف اکثر مطالعات، به عنوان مثال [ 6 ] [ 16 ] [ 17 ] [ 18 ] [ 39 ] و [ 40 ]، در حالی که پایداری را ادغام می کند تا دیگری را نشان دهد. رویکرد نسبت به موارد [ 59 ] [ 60 ] . مورد دوم فنی ترین است، با هدف قرار دادن بهترین ترکیب روش ها برای تجسم سه ملاحظات نظری مذکور در همان تجزیه و تحلیل GIS، همانطور که توسط [ 59 ] و [ 60 ] امتحان شده است.
2.3. منطقه مطالعه قبل از اجرای روشها
منطقه مورد مطالعه ناحیه Njimom است که بین عرضهای جغرافیایی شمالی 5˚40′-6˚ و شرقی 10˚48′-11˚20’، در بخش اسمی و منطقه غربی کامرون واقع شده است ( شکل 4 ). مساحت آن 780 کیلومتر مربع است، در قلب ارتفاعات انتقالی بین جنوب و شمال کامرون که از نظر آب و هوا و ویژگی های بیوفیزیکی متضاد هستند. دمای متوسط 25 درجه سانتیگراد السیوس و 1200 تا 2000 میلی متر بارندگی مشخصه آب و هوای سالانه آن است. پوشش گیاهی بین جزایر جنگلی کمکم، مراتع و درختچهها مخلوط است. ارتفاع بین 650 تا 1533 است و نقطه اوج آن کوه یاوو است. خاکها آتشفشانی کمتر تکامل یافتهاند و بسته به ارتفاع، روی شیبها یا خاک رس مرطوب نزدیک رودخانهها قرمز تیره ناپایدار هستند. شبکه مهمی از رودخانه های کوچک هیدروگرافی را تشکیل می دهد و رودخانه اصلی مگنا است. شبکه راه ها عمدتاً از نوع روستایی است و راه ملی شماره 6 بخش شمال کشور را باز می کند. فعالیت های اقتصادی عمدتاً کشاورزی، شکار و تجارت است.
2.4. شرایط-تولید لایه ها، رتبه بندی و گروه بندی
لایه ها بر اساس مقیاس های استاندارد (نگاه کنید به جدول 2 )، و همچنین نیازهای سرمایه گذاران و ویژگی های منطقه مورد مطالعه تولید شدند. داده های اصلی از انواع شطرنجی و برداری هستند. تصویر ماهوارهای Landsat 8 OLI-TIRS، با مختصات مسیر و ردیف 186-56، مربوط به 27 ژانویه 2018 و منطقه مورد مطالعه را پوشش میدهد، از وبسایت رسمی سازمان زمینشناسی ایالات متحده (USGS) به نشانی https://earthexplorer دانلود شده است. usgs.gov/ ( جدول 5 ).
شکل 4 . محل مطالعه.
همانطور که تمام محصولات سطح 1 خام هستند، تصویر نیاز به پیش پردازش دارد. نوارهای آبی (B2) تا مادون قرمز با طول موج کوتاه (B7) مورد استفاده در این مطالعه انباشته شدند. سپس کالیبراسیون رادیومتریک از 16 بیت به 8 بیت و اصلاحات جوی انجام شد. پس از آن، برای افزایش وضوح فضایی، نوار پانکروماتیک (B8) از قبل پردازش شد و به تصویر انباشته چند طیفی ادغام شد.
به همین ترتیب، مدل ارتفاعی دیجیتالی مدل رادار شاتل (SRTM) که در وبسایت USGS نیز موجود است، دانلود و برای برآوردن الزامات آنالیز شرطی شد. با این دو مجموعه داده و بردارهای موجود، می توان تحلیل را انجام داد. نرم افزار مورد استفاده Arcgis 10.4.1، Erdas Imagine 2015 و Super Decision Software (SDS) نسخه 3.0 است.
• دما و بارندگی – برخی از داده های پراکنده موجود در اداره محلی کشاورزی، با ثبت مستقیم دما در مزرعه هر 300 یا 500 متر جمع آوری و تکمیل شد. آنها بین 22 تا 27 درجه سانتیگراد در طول سال 2018 متغیر هستند. برای بارندگی، نقشه منطقه مورد مطالعه چاپ شد و از داده های مناطق اطراف مجهز به ایستگاه های هواشناسی برای تخمین بارش های سالانه محلی استفاده شد. اداره کشاورزی سپس بر اساس داده های فومبان (جنوب)، ماسنگام (جنوب شرقی)، مگبا (شمال شرقی) و بامندا (شمال غربی)، ناحیه نجیموم به چهار ناحیه بارندگی سالانه در فاصله 1206 میلی متر و 2145 تقسیم شده است. میلی متر ( شکل 5 ). هر دو داده با استفاده از عملیات چگالی نقطه شطرنجی شدند.
• ارتفاع، شیب و آفتاب – این سه مجموعه داده با استفاده از DEM-SRTM شرطی شده تولید شدند. ارتفاع مورد نظر با اعمال یک آستانه به دست آمد. شیب بر حسب درجه تولید شد و یک آستانه اعمال شد. برای تابش آفتاب، تابش خورشیدی منطقه انجام شد. در نرم افزار ArcGIS، این شامل استخراج تابش خورشیدی ورودی از یک سطح شطرنجی است. فاصله درخواستی 365 روز برای سال 2018 بود، واحد وات ساعت بر متر مربع (W∙H/m 2 ) است، در حالی که عرض جغرافیایی به دقت شیب خورشید و موقعیت خورشید کمک می کند.
• پوشش جنگلی – این لایه با ترکیب سه شاخص طیفی که معمولاً در مدل تراکم تاج پوشش جنگل (FCD) استفاده می شود، تولید شد [ 61 ].
– شاخص گیاهی پیشرفته (AVI)، که از درجه قدرت پاسخ مادون قرمز برای برجسته کردن تفاوت های ظریف در تراکم تاج پوشش استفاده می کند. (معادله (5))
A Vمن=( ( B 5 + 1 ) ∗ ( 256 − B 4 ) ∗ ( B 5 − B 4 ) )—————————–√3AVI=((B5+1)∗(256−B4)∗(B5−B4))3(5)
– Canopy Shadow Index یا Scaled Shadow Index (SI, SSI)، که مخالف توده های جنگلی و محورهای طیفی ناحیه باز است. (معادله (6))
اسمن=( 256 − B 3 ) ∗ ( 256 − B 4 )——————-√SI=(256−B3)∗(256−B4)(6)
– شاخص خاک لخت (BI)، افزایش مناطق خاک برهنه، زمین های آیش، پوشش گیاهی با پس زمینه مشخص، در مخالفت با AVI. (معادله (7))
B I=( B 5 + B 3 ) – B 4( B 5 + B 3 ) + B 4BI=(B5+B3)−B4(B5+B3)+B4(7)
در ماشین حساب شطرنجی، مجموع AVI و SI برای افزایش حس علامت جنگل محاسبه شد. سپس، BI از آن لایه جدید کم شد تا اثرات خاک به حداقل برسد.
• کیفیت خاک – دو نئو باند تولید شد. برای ارزیابی غنای خاک، سه شاخص برگرفته از [ 62 ] و [ 63 ] به دنبال معادلات (8)، (9) و (10) محاسبه شد.
مواد معدنی خاک رس سیم= B 6 / B 7CM=B6/B7(8)
شاخص ترکیب خاک اسسیمن=( B 6 – B 5 )( B 6 + B 5 )SCI=(B6−B5)(B6+B5)(9)
شاخص مواد آلی ای اممن= 1 / B32OMI=1/B32(10)
سه لایه انباشته شدند، تجزیه و تحلیل مؤلفه اصلی (PCA) برای استخراج بهترین اطلاعات انجام شد، و مؤلفه ای که به ترتیب به CM، SCI و OMI مرتبط بود، به عنوان بهترین بیان غنای خاک حفظ شد
شکل 5 . نمودارهای بارندگی و منحنی های دما 2018 فومبان، ماسنگام، مگبا و بامندا.
در همین منطق، رطوبت خاک با محاسبه شاخص سطح آب زمین (LWSI) در رابطه (11) مورد ارزیابی قرار گرفت.
L Wاسمن= B 5 − B 6 / B 5 + B 6LWSI=B5−B6/B5+B6(11)
سپس شبکه استریم دیجیتالی شد. تجزیه و تحلیل آمار کانونی به دنبال همسایگی شعاع 20 متری انجام شد. برای پایان دادن به این مرحله، تراکم زهکشی منطقه مورد مطالعه طبق رابطه (12) محاسبه شد.
D D = طول شبکه جریان ( کیلومتر ) / منطقه مطالعه (کیلومتر2)DD=Streamnetworklength(km)/StudyArea(km2)(12)
با توجه به اینکه خاک ها باید زهکشی خوبی داشته باشند تا بتوانند ریشه های خوب ثابتی داشته باشند، اما به طور دائم زهکشی نشده باشند تا رشد و برداشت را به خطر بیندازند [ 11 ]، لایه های LWSI و DD در یک ماشین حساب شطرنجی ادغام شدند و به ترتیب 2/3 (67٪) وزن شدند. و 1/3 (33%)، برای به دست آوردن رطوبت خاک.
• ساخته شده – برخی از مطالعات نشان داده اند که شاخص های ساخته شده اطلاعات دقیق تری نسبت به طبقه بندی تحت نظارت یا بدون نظارت می دهند [ 64 ]. به عنوان مثال، شاخص ساخت و ساز خشک (DBI) که اخیراً توسعه یافته است در برجسته کردن ویژگیهای شهرها در آب و هوای خشک مهم است [ 65 ]. فرمول آن با رابطه (13) به دست می آید.
D B I= ( B 2 – B 10 / B 2 + B 10 ) – ND VمنDBI=(B2−B10/B2+B10)−NDVI(13)
با NDVI به معنی شاخص گیاهی تفاوت نرمال شده و به دنبال معادله (14) بیان می شود.
نD Vمن= ( B 5 – B 4 / B 5 + B 4 )NDVI=(B5−B4/B5+B4)(14)
با این حال، منطقه مورد مطالعه دقیقاً منطقهای خشک با خاکهای ماسهای و تپهها و پوشش گیاهی استپی کمکم نیست، بلکه ویژگیهای طبیعی مانند درختچهها و علفزارها، خاکهای قرمز تیره ناپایدار و مرطوب بسته به ارتفاع است [ 65 ]، در حالی که دهها خانه ساخته شده است. با مواد خام مانند آجر و کاه بر آن مسلط است. هنگامی که نتایج با طبقهبندی نظارت شده و شاخصهای پوشش گیاهی بهعنوان AVI و NDVI مقایسه شدند، این شرایط باعث ایجاد برخی سردرگمیها و وارونگی بازتاب پیکسلها شد. سپس، وارونگی عبارت نسبت در فرمول اولیه، همانطور که در رابطه (15) نشان داده شده است، برای به دست آوردن نتایج مورد انتظار، که استخراج ساخته شده است، آزمایش شد:
منD B I= ( B 10 – B 2 / B 10 + B 2 ) – ND VمنIDBI=(B10−B2/B10+B2)−NDVI(15)
بر اساس آستانه مقادیر، نتیجه در دو کلاس مناسبت 1 و نامناسب بودن 0 کدگذاری شد.
• شبکه جاده ها – با فرض اینکه طراحی یک مزرعه نخل روغنی مدرن، نزدیکی یک کارخانه برای تحول مفید را درخواست می کند [ 11 ]، حمل و نقل آسان محصولات نیاز به یک شبکه جاده ای مبله دارد. منطقه مورد مطالعه به لطف مسیرهای روستایی و جاده ملی 6 باز شده است ( شکل 1 را ببینید ) و گفته شد که کارخانه مزرعه باید 1 کیلومتر به جاده نزدیک باشد. شبکه جاده با استفاده از تابع فاصله اقلیدسی در 4000 متر دیجیتالی و شطرنجی شد تا کل منطقه مورد مطالعه را پوشش دهد. سپس یک آستانه 1000 متری اعمال شد و در دو کلاس مناسبت، 1 در داخل این فاصله و نامناسب بودن 0 که بیش از این فاصله است، مجددا کدگذاری شد.
• رودخانه ها – هدف حفظ موقعیت فعلی واحدهای پیش ساخته روغن نخل محلی و سایر کارخانه های صنایع دستی بود. برای این منظور، بر روی یک نقشه چاپی که با نشانگرهای فضایی به عنوان رودخانه ها و تپه ها طراحی شده است، نقشه های شناختی دستی در سرتاسر منطقه با جمعیت هایی با پیروی از رویکرد مشارکتی GIS (PGIS) پیشنهاد شده توسط [ 66 ] و مختصات جغرافیایی یازده واحد ترسیم شده است. با GPS ضبط شده اند. یک تصویر شطرنجی با استفاده از تابع فاصله اقلیدیان در 1500 متر برای پوشش منطقه مورد مطالعه ایجاد شد. یک آستانه 100 متری برای شطرنج حاصل به عنوان محیط محدود به رودخانه ها با توجه به درخواست جمعیت اعمال شد. تصویر باینری دارای 1 به عنوان مناسب برای حفاظت از آسیاب در این فاصله، و 0 به عنوان نامناسب است که برابر با کمتر از 100 متر است.
از نظر فنی، لایه ها در قابلیت طبیعی و مناسب اجتماعی-اقتصادی مونتاژ شدند. از یک طرف، گروهی از لایههای متعلق به قابلیت طبیعی مجدداً طبقهبندی شدند و در مقیاس ارزیابی از 0 تا 5 که نشاندهنده کمترین تا بالاترین قابلیت است، کدگذاری مجدد شدند. و از آنجا که آنها شرایط ثابت هستند یا برای بهره وری کمتر انعطاف پذیر هستند، لایه ها به روش ساده روش مجموع WLC برای به دست آوردن “لایه قابلیت” تلاقی داده شدند. از سوی دیگر، لایههای شرایط اجتماعی-اقتصادی پس از انجام فرآیند AHP تجمیع شدند. ماتریس مقایسه زوجی ساخته شده است که بیان می کند که آسیاب های محلی در اطراف شبکه رودخانه 3 برابر کمتر از کارخانه مدرن اهمیت دارند، اما 5 برابر بیشتر از کارخانه های ساخته شده. در حالی که نزدیکی بین جاده شبکه و اهمیت کارخانه مدرن بالاترین است،
شبکه رودخانه شبکه جاده ای ساخته شده ⎡⎣⎢131/5 _ _1/3 _ _11/7 _ _571⎤⎦⎥RivernetworkRoadnetworkBuilt-up [11/353171/51/71](16)
این وزن های نسبی برای به دست آوردن بردار اولویت نرمال شدند. سازگاری قضاوت ها از طریق نسبت شاخص سازگاری (CI) به شاخص تصادفی (RI) در رابطه (17) ارزیابی شد.
سیR = ( Cمن/ آر آی) < 0.10CR=(CI/RI)<0.10(17)
“لایه مناسب” از مجموع این سه مولفه وزن شده توسط بردار اولویت ثابت آنها در یک ماشین حساب شطرنجی حاصل شده است. الگوهای بصری و نتایج عددی با ارسال تصویر خروجی به عملیات عضویت فازی خطی تایید شد.
• پس از تولید “لایه قابلیت” و “لایه مناسب”، ارزیابی سطح در نظر گرفتن پایداری مهم بود. گام اول شامل شمارش تعداد رخدادهای بعد اجتماعی، اقتصادی یا محیطی در هر پردازش لایه است ( جدول 6 ). سپس پس از آن، تابع مطلوبیت برای دستیابی به مقادیر پایداری فردی آنها اعمال شد.
در هر دو مورد، x 1 ، x 2 و x 3 به ترتیب به جنبه های اجتماعی، اقتصادی و زیست محیطی پایداری اشاره دارد. در داخل هر گروه از لایهها، نسبت بین هر جنبه از تعداد رخدادها و کل ایجاد شده است. بنابراین، ریشه مکعب همه نتایج محاسبه شد و مقادیر جدید برای یافتن ضریب مطلوبیت پایداری (SUC) برای هر گروه ضرب شدند. مقادیر غیر اتفاقی یا صفر به صورت 1 محاسبه شد.
به عنوان آخرین مرحله، در یک ماشین حساب شطرنجی، هر لایه قبل از انجام مجموع وزنی آنها در ضریب مطلوبیت پایداری خود (SUC) ضرب شد.
فلوچارت نمودار ارائه شده در شکل 6 مراحل اصلی روش را ترکیب می کند.
شکل 6 . فلوچارت متدولوژی.
3. نتایج
3.1. شرایط-لایه ها و ارزش های خام و کلاس ها
به طور کلی یازده نقشه لایه شرایط فردی از پردازش بالا حاصل می شود. هشت مورد از آنها متعلق به قابلیت طبیعی و سه مورد به مناسبت اجتماعی-اقتصادی برای کاشت نخل روغنی است. در داخل دو گروه لایه های غیر ترکیبی و ترکیبی یا نئو باند وجود دارد.
برای قابلیت طبیعی، لایه های غیر مرکب عبارتند از بارندگی (1206 – 2145 میلی متر)، دما (20.2 درجه سانتیگراد – 27 درجه سانتیگراد)، آفتاب (13.58 – 2103508)، شیب (0˚ – 5.41 درجه) و ارتفاع (696 -). 1533 متر) ( شکل 7 ). آنها مقادیر واقعی شرایط واقعی را برای نخل روغنی می دهند که تقریباً در اطراف آنهایی که نشان داده شده است رشد می کند [ 11 ]. لایه های مرکب یا نئو باند، به عنوان غنای خاک (-0.0618 – 0.00015)، رطوبت خاک (3 – 12) و پوشش جنگل (0.67 – 258.7) ( شکل 7 ). آنها با ترکیب شبیه سازی های مختلف، شرایط واقعی را تقریب می کنند.
شکل 7 . شرایط قابلیت طبیعی-لایه ها.
با توجه به تناسب اجتماعی-اقتصادی، لایههای غیر مرکب فاصله تا شبکه راه (0-3837 متر) و تا رودخانه (0-1500) است ( شکل 8 ). آنها نیازهای سرمایه گذاران و تصور جمعیت را بیان می کنند. تنها لایه تشکیلشده، نقشه ساختهشده (-1.37 – 1.477) است، که مقادیر آن تقریباً منطقهای است که باید توسط مزرعه اجتناب شود ( شکل 8 ).
در مقایسه با استانداردهای پرورش نخل روغنی [ 11 ]، این شرایط عموماً برای کاشت نخل روغنی در منطقه به طور متوسط برآورده می شوند ( جدول 2 را ببینید). فقط ارتفاع (696 – 1533 متر)، که بهتر است بین 400 تا 500 متر باشد، و همچنین خاک هایی که بهترین محتوای معدنی آنها رس است، برخی محدودیت های دید اول را نشان می دهد ( جدول 2 را ببینید ). سپس، همانطور که در بالا گفته شد، بهترین راه برای ارزیابی آن، رمزگذاری مجدد لایه ها بر اساس مقیاس استاندارد بود.
لایههای شرایط قابلیت طبیعی به شش کلاس جدید مجدداً طبقهبندی میشوند. مقیاس زبانی فازی رو به کاهش است، با قابلیت «بیشترین»، «بالا»، «متوسط»، «حاشیه»، «کم» و «کمترین». بر اساس دادههای پردازششده اولیه، همه لایههای شرایط بالاترین تا پایینترین مقادیر، بهطور کامل مدل را ادغام کردند، به جز ارتفاعی که رتبهبندی دادهها معکوس شد، و بالاترین امتیاز و واجد شرایط را به پایینترین ارتفاعات داد. برای لایههایی که دادههای افراطی وجود نداشت، برای نتایج بهتر، تعداد کلاسها با مقادیر امتیاز تطبیق داده شد ( جدول 7 ، شکل 9 ).
شکل 8 . شرایط مناسب اجتماعی-اقتصادی-لایه ها.
شکل 9 . شرایط قابلیت طبیعی-لایه ها دوباره طبقه بندی شدند.
لایههای شرایط مناسب اقتصادی-اجتماعی دارای دو طبقه «مناسب» و «نامناسب» هستند. بالاترین مقادیر ساخته شده برای مزرعه روغنی نامناسب است. فاصله 4000 متری اطراف جاده برای دسترسی آسان برای کارخانه مناسب است. فاصله بیش از 100 متر از رودخانه برای کشاورزی و کارخانه برای حفظ فعالیت تولید نخل روغن محلی و سپس جلوگیری از درگیری با جمعیت محلی مناسب است ( جدول 8 ، شکل 10 ).
3.2. توانمندی و لایه های ترکیبی اجتماعی-اقتصادی
لایه قابلیت حاصل از WLC روندهایی را با مقادیر 23 تا 45 نشان می دهد. بیشترین قابلیت در قسمت شمالی منطقه مورد مطالعه ثبت شده است. نمرات آن با بالاترین غنای خاک، رطوبت خاک و پوشش جنگلی، متوسط آفتاب، بارندگی و دما، کمترین ارتفاع و شیب مطابقت دارد ( شکل 11 ).
لایه مناسب نتیجه عملکرد AHP است ( جدول 9 ).
سپس، سه لایه رمزگذاری شده مجدد در بردارهای اولویت آنها ضرب می شوند، به ترتیب لایه شبکه رودخانه با 0.2828، لایه شبکه جاده با 0.6434 و لایه ساخته شده با 0.0738. مجموع عضویت فازی خطی آنها در بازه بین 0 و 1، با غلبه قابل مشاهده بر لایه شبکه راه است ( شکل 11 ).
3.3. توزین ضرایب سودمند پایداری (SUC) و بهترین نقشه های سایت نهایی
پس از محاسبات UF، هر لایه یک SUC رتبه بندی شده در مقیاس فازی 0 تا 1 به دست می آورد. هیچ یک از هر دو گروه ابعاد پایداری را در 100% برآورده نمی کنند. پردازش قابلیت طبیعی آن را در 54.2% بیان می کند که مقدار نرمال شده 0.542 است، در حالی که پردازش اجتماعی-اقتصادی آن را در 31.5% بیان می کند که یک مقدار نرمال شده 0.315 است ( جدول 10 ).
این دو لایه در شش کلاس برای مقیاس نقشه نهایی ارزیابی مجدد کدگذاری می شوند. سپس با استفاده از مقادیر نرمال شده SUC، لایه قابلیت در 0.542 و لایه مناسب در 0.315 ضرب می شود. مجموع دو لایه در ماشین حساب شطرنجی آخرین نقشه را داده بود.
شکل 10 . شرایط مناسب اقتصادی-اجتماعی طبقه بندی مجدد لایه ها.
شکل 11. نقشه های قابلیت و تناسب ترکیب و طبقه بندی مجدد شده اند.
نقشه نهایی در چهار طرح مختلف گذاشته شده است ( شکل 12 ). مورد اول روندهای اصلی را از بهترین تا کمترین یا بدترین شرایط نشان می دهد. مقادیر در بهترین حالت بین 5.142 و در بدترین حالت 0.857 است. منطقه اصلی برای ارزش های بالا در شمال منطقه Njimom است. این نقشه در شش کلاس برای تفکیک و اولویت بندی سایت ها در نقشه دوم مجددا کدگذاری شده است. این رتبه بندی از اولویت اول تا ششم برای سرمایه گذار است، با مساحتی بین 35700 هکتار تا 153 هکتار ( شکل 13).). هنگام اعمال رمزگذاری مجدد باینری بر اساس منطق “بله”/”خیر” در این نقشه، با در نظر گرفتن بالاترین کلاس معتبر و سایرین به عنوان نامعتبر، این امکان را فراهم می کند تا بهترین سایت منطقه مورد مطالعه را در نقشه سوم برجسته کنید. منطقه اولویت اول بهترین سایت در مساحت 35700 هکتار است. ده مزرعه بیش از یازده مربوط به آسیاب های محلی مورد استفاده برای تجزیه و تحلیل در محدوده سایت قرار دارند. مجموع حداقل مساحت تخمینی آنها، هنگام تصرف به عنوان 7.5 هکتار طبق [ 11 ] برای کشاورزان کوچک (7.5 × 10 = 750) (به جدول 2 مراجعه کنید )، از این 35700 هکتار کسر می شود، و مساحت کل تقریبی قابل بهره برداری حدود است. 34950 هکتار; که 44.8 درصد از 78000 هکتار (780 کیلومتر مربع ) مساحت کل ناحیه است ( شکل 14)). سپس محل انتخاب شده برای مزارع نخل روغنی با اندازه بزرگ قابل استفاده است ( جدول 2 را ببینید ). در نهایت، حد بهترین سایت ترسیم می شود و مرکزی ترین موقعیت چند ضلعی برای قرار دادن و نصب کارخانه مدرن در داخل منطقه کاشت و مربوط به یازده آسیاب محلی ثبت شده در منطقه انتخاب می شود.
شکل 12 . بهترین منطقه سایت
شکل 13 . شش طبقه از مناطق دارای اولویت در هکتار.
شکل 14 . نسبت بهترین مساحت سایت نسبت به سایرین.
4. بحث ها و چشم انداز
با توجه به هدف اصلی انتخاب بهترین محل کاشت نخل روغنی، این مطالعه با رعایت حداکثر استانداردهای رتبه بندی بین المللی انجام شده است. نتیجه اصلی پردازش مکان بالقوه نشاندادهشدهترین مکان را در مناطق پست شمالی ناحیه Njimom نشان میدهد که در شکل 12 نشان داده شده است. مقایسه بصری آن با لایههای ورودی، تأثیر زیادی از روندهای «لایه قابلیت» را نشان میدهد ( شکل 11 را ببینید ). یک گام به عقب و همچنان بر اساس مقایسه بصری، هفت لایه در هشت شرط، منطقه شمالی را به عنوان مطلوبترین، به شدت برجسته میکنند، به جز لایه پوشش جنگلی که مقادیر متوسط و پایین را در آنجا ثبت میکند ( شکل 7 و شکل 9 را ببینید.) در حالی که لایهها بهعنوان شیب و ارتفاع به وضوح، اما همچنین بارندگی و دما، عمیقاً بر شکل مکان انتخابی تأثیر میگذارند.
علاوه بر این، توضیح اصلی در مورد تسلط روندهای “لایه قابلیت” بر روی “لایه مناسب” در انتخاب بهترین سایت، فرآیند توزین با استفاده از ضرایب مطلوبیت پایداری (SUC) است. در واقع، با بیش از 50٪ (0.542) و فقط 30٪ (0.315)، نتیجه بیشتر نماینده لایه اول نسبت به لایه دوم است. در نتیجه، هنگام برآورد مناطق دارای بالاترین ارزش قابلیت ( شکل 11 ) و بهترین مساحت سایت، آنها به ترتیب 23349 هکتار و 34950 هکتار هستند. تفاوت 11601 هکتار نشان دهنده سهم “لایه مناسب” در انتخاب نهایی است که 33.2٪ یا 1/3 از کل مساحت است، در مقابل 66.8٪ یا 2/3 که توسط مقدار “لایه قابلیت” بالا نشان داده شده است.
• n نکته دیگر، اگرچه سایت با توجه به شرایط موجود محلی مناسب ترین است، اما 100% را برآورده نمی کند، زیرا لایه های شیب و ارتفاع بالاترین رتبه 5 را ثبت نمی کنند. سپس بالاترین رتبه امتیاز نهایی را ثبت می کند. نقشه بر اساس استانداردهای بین المللی حداکثر 4 می باشد.
به همین ترتیب، در مورد هدف فرسایش، استانداردهای گفته شده و تمام مطالعات مطالعه شده در مورد هدف تناسب زمین نخل روغنی بیشتر بر اجتناب از هر گونه خطری تمرکز می کنند تا ادغام آن با فرآیند مدل سازی. برای ادامه این مطالعه، چشم انداز اصلی از یک سو ساختن یک مدل کامل از معادله جهانی تلفات خاک (RUSLE) و ادغام آن با مدل است. از سوی دیگر، ادغام وضعیت تخریب زمین در فرآیند، پیروی از مدل ارائه شده توسط [ 67 ] می تواند اطلاعات بیشتری را برای انتخاب مکان به ارمغان آورد. این دو جزء ممکن است نتیجه را تقویت یا کاهش دهند.
تجزیه و تحلیل همچنین راه را برای اعمال DBI برای استخراج ساخته شده و در نهایت برای آزمایش نسبت معکوس این فرمول برای نتایج بهتر باز کرده است (به معادلات (13) و (15) مراجعه کنید). این نکتهای است که نویسندگان قبلاً روی آن کار میکنند تا درک پارامترهای بازتابی ساختهشده در مواد خام مربوط به باندهای مادون قرمز حرارتی و آبی طیف الکترومغناطیسی را در سنجش از دور شهرهای مختلف با آب و هوای خشک عمیقتر کنند.
5. نتیجه گیری
این مقاله برای ترکیب روشهای WLC، FAHP و UF در یک محیط GIS آزمایش کرده است. هدف انتخاب بهترین محل کاشت نخل روغنی بود. سه جهت قابلیت و تناسب بود که در راس آن ایده پایداری بود. سپس به نظر می رسد که در ناحیه Njimom، با ویژگی های بیوفیزیکی انتقالی بین آب و هوای استوایی و گرمسیری، برخی از الزامات برای کاشت نخل روغنی در ضلع شمالی آن برآورده شده است. بنابراین، استفاده از رویکرد ساده WLC به پیروی از مقیاس استاندارد شده شرایط مناسب مانند بارندگی، دما، آفتاب، شیب، ارتفاع، غنای خاک و رطوبت و پوشش جنگل کمک میکند که برای رشد و برداشت بهتر باید رعایت شود. به همین ترتیب، استفاده از تجزیه و تحلیل انعطاف پذیر به عنوان FAHP به ادغام افکار و انتخاب های انسانی در تجزیه و تحلیل به عنوان ساخته شده کمک می کند. محل کارخانه یا کارخانه های جمعیت محلی. افزایش دیگر سطح در نظر گرفتن پایداری 0.857 است که بیشتر به عنوان ارزیابی کننده ابعاد کلی اجتماعی، اقتصادی و زیست محیطی یک پروژه کشتزار استفاده می شود. در نهایت، ترکیب این سه روش بهترین مکان را برای کاشت نخل روغنی در زمینی به مساحت 34950 هکتار به ارمغان آورده است و منطقه نجیموم را مستعد میزبانی کشاورزی مدرن و پایدار برای توسعه محلی می کند.
بدون دیدگاه