i.landsat.toargrass - ক্লাউডে অনলাইন

কার্যক্রম:

NAME এর

i.landsat.toar - বায়ুমণ্ডলের উচ্চতা বা প্রতিফলন এবং তাপমাত্রা গণনা করে
Landsat MSS/TM/ETM+/OLI

KEYWORDS

চিত্র, রেডিওমেট্রিক রূপান্তর, উজ্জ্বলতা, প্রতিফলন, উজ্জ্বলতা তাপমাত্রা, ল্যান্ডস্যাট,
বায়ুমণ্ডলীয় সংশোধন

সাইনোপিসিস

i.landsat.toar
i.landsat.toar --help
i.landsat.toar [-rnp] ইনপুট=বেসনাম আউটপুট=বেসনাম [মেটফাইল=নাম] [সেন্সর=স্ট্রিং]
[পদ্ধতি=স্ট্রিং] [তারিখ=YYYY-MM-DD] [সূর্যের_উচ্চতা=ভাসা] [পণ্য_তারিখ=YYYY-MM-DD]
[লাভ করা=স্ট্রিং] [শতাংশ=ভাসা] [পিক্সেল=পূর্ণসংখ্যা] [রেলেগ=ভাসা]
[lsatmet=স্ট্রিং[,স্ট্রিং,...]] [স্কেল=ভাসা] [---ওভাররাইট] [---সাহায্য] [---ভার্বোস]
[---শান্ত] [---ui]

পতাকা:
-r
সমস্ত ব্যান্ডের জন্য প্রতিফলনের পরিবর্তে আউটপুট অ্যাট-সেন্সর উজ্জ্বলতা

-n
ইনপুট রাস্টার মানচিত্রগুলি কোডের পরিবর্তে ব্যান্ডের সংখ্যাকে এক্সটেনশন হিসাবে ব্যবহার করে

-p
প্রিন্ট আউটপুট মেটাডেটা তথ্য

--ওভাররাইট করুন
আউটপুট ফাইলগুলিকে বিদ্যমান ফাইলগুলিকে ওভাররাইট করার অনুমতি দিন

--help
প্রিন্ট ব্যবহারের সারাংশ

-- ভারবোস
ভার্বোস মডিউল আউটপুট

-- শান্ত
শান্ত মডিউল আউটপুট

--ui
জোর করে GUI ডায়ালগ চালু করুন

পরামিতি:
ইনপুট=বেসনাম [প্রয়োজনীয়]
ইনপুট রাস্টার ব্যান্ডের মূল নাম
উদাহরণ: 'বি.' B.1, B.2,... এর জন্য

আউটপুট=বেসনাম [প্রয়োজনীয়]
আউটপুট রাস্টার মানচিত্রের উপসর্গ
উদাহরণ: 'B.toar.' B.toar.1, B.toar.2, ... তৈরি করে

মেটফাইল=নাম
Landsat মেটাডেটা ফাইলের নাম (.met বা MTL.txt)

সেন্সর=স্ট্রিং
মহাকাশযান সেন্সর
শুধুমাত্র 'মেটফাইল' না দিলেই প্রয়োজন (স্যানিটির জন্য প্রস্তাবিত)
বিকল্প: mss1, mss2, mss3, mss4, mss5, tm4, tm5, tm7, oli8
mss1: Landsat-1 MSS
mss2: Landsat-2 MSS
mss3: Landsat-3 MSS
mss4: Landsat-4 MSS
mss5: Landsat-5 MSS
tm4: Landsat-4 TM
tm5: Landsat-5 TM
tm7: Landsat-7 ETM+
oli8: Landsat_8 OLI/TIRS

পদ্ধতি=স্ট্রিং
বায়ুমণ্ডলীয় সংশোধন পদ্ধতি
বায়ুমণ্ডলীয় সংশোধন পদ্ধতি
বিকল্প: অসংশোধিত, dos1, dos2, dos2b, dos3, দুই৩৯
ডিফল্ট: অসংশোধিত

তারিখ=YYYY-MM-DD
ছবি অধিগ্রহণের তারিখ (yyyy-mm-dd)
'মেটফাইল' না দিলেই প্রয়োজন

সূর্যের_উচ্চতা=ভাসা
ডিগ্রিতে সূর্যের উচ্চতা
'মেটফাইল' না দিলেই প্রয়োজন

পণ্য_তারিখ=YYYY-MM-DD
ছবি তৈরির তারিখ (yyyy-mm-dd)
'মেটফাইল' না দিলেই প্রয়োজন

লাভ করা=স্ট্রিং
সমস্ত Landsat ETM+ ব্যান্ডের লাভ (H/L) (1-5,61,62,7,8)
'মেটফাইল' না দিলেই প্রয়োজন

শতাংশ=ভাসা
পথের দীপ্তিতে সৌর দীপ্তির শতাংশ
শুধুমাত্র 'পদ্ধতি' কোনো DOS হলেই প্রয়োজন
ডিফল্ট: 0.01

পিক্সেল=পূর্ণসংখ্যা
ন্যূনতম পিক্সেল ডিজিটাল সংখ্যাকে অন্ধকার বস্তু হিসাবে বিবেচনা করতে হবে
শুধুমাত্র 'পদ্ধতি' কোনো DOS হলেই প্রয়োজন
ডিফল্ট: 1000

রেলেগ=ভাসা
Rayleigh বায়ুমণ্ডল (বিস্তৃত আকাশ বিকিরণ)
শুধুমাত্র 'পদ্ধতি' DOS3 হলেই প্রয়োজনীয়
ডিফল্ট: 0.0

lsatmet=স্ট্রিং [, স্ট্রিং,...]
প্রদত্ত মেটাডেটার জন্য সংরক্ষিত রিটার্ন মান
শুধুমাত্র 'metfile' এবং -p দিলেই প্রয়োজন
বিকল্প: সংখ্যা, সৃষ্টি, তারিখ, সূর্যালোক, সেন্সর, ব্যান্ড সুনাজ, সময়
সংখ্যা: ল্যান্ডস্যাট নম্বর
সৃষ্টি: তৈরির সময় স্ট্যাম্প
তারিখ: তারিখ
sun_elev: সূর্যের উচ্চতা
সেন্সর: সেন্সর
ব্যান্ড: ব্যান্ড গণনা
সুনাজ: সূর্য আজিমুথ কোণ
সময়: সময়

স্কেল=ভাসা
আউটপুট জন্য স্কেল ফ্যাক্টর
ডিফল্ট: 1.0

বর্ণনাঃ

i.landsat.toar ল্যান্ডস্যাট চিত্রের ক্রমাঙ্কিত ডিজিটাল সংখ্যা রূপান্তর করতে ব্যবহৃত হয়
বায়ুমণ্ডলের শীর্ষে থাকা পণ্য বা বায়ুমণ্ডলের শীর্ষে প্রতিফলন এবং তাপমাত্রা
(TM এবং ETM+ সেন্সরের ব্যান্ড 6)। ঐচ্ছিকভাবে, এটি গণনা করতে ব্যবহার করা যেতে পারে
বায়ুমণ্ডলীয় সংশোধন (DOS পদ্ধতি) সহ পৃষ্ঠের উজ্জ্বলতা বা প্রতিফলন।

সাধারণত, এটি করার জন্য উত্পাদন তারিখ, অধিগ্রহণের তারিখ এবং সৌর উচ্চতা হয়
প্রয়োজন তাছাড়া, Landsat-7 ETM+ এর জন্য নয়টির লাভ (উচ্চ বা নিম্ন) প্রয়োজন
নিজ নিজ ব্যান্ড

ঐচ্ছিকভাবে (প্রস্তাবিত), ডেটা মেটাডেটা ফাইল (.met বা MTL.txt) থেকে পড়া যেতে পারে
সমস্ত Landsat MSS, TM, ETM+ এবং OLI/TIRS। তবে সৌর উচ্চতা থাকলে মান দেওয়া হয়
মেটাডেটা ফাইল ওভাররাইট করা হয়. এটি প্রয়োজনীয় যখন .met ফাইলে ডেটা থাকে
ভুল বা সঠিক নয়। এছাড়াও, যদি অধিগ্রহণ বা উত্পাদন তারিখ পাওয়া যায় না
মেটাডেটা ফাইল তারপর কমান্ড লাইন মান ব্যবহার করা হয়.

দৃষ্টি আকর্ষণ করছি: ইনপুট রাস্টারে QCALmin এর চেয়ে ছোট যে কোনো নাল মান নাল ইন সেট করা আছে
আউটপুট রাস্টার এবং এটি সমীকরণের অন্তর্ভুক্ত নয়।

অসংশোধিত এ-সেন্সর মূল্যবোধ (পদ্ধতি=অসংশোধিত, ডিফল্ট)

স্ট্যান্ডার্ড জ্যামিতিক এবং রেডিওমেট্রিক সংশোধনের ফলে একটি ক্রমাঙ্কিত ডিজিটাল সংখ্যা হয়
(QCAL = DN) ছবি। আলোকসজ্জা জ্যামিতির প্রভাবকে আরও মানসম্মত করতে, QCAL
ছবিগুলি প্রথমে অ্যাট-সেন্সর রেডিয়েন্সে এবং তারপরে অ্যাট-সেন্সর প্রতিফলনে রূপান্তরিত হয়।
থার্মাল ব্যান্ডটি প্রথমে QCAL থেকে অ্যাট-সেন্সর রেডিয়েন্সে রূপান্তরিত হয় এবং তারপর কার্যকর হয়
কেলভিন ডিগ্রীতে সেন্সর তাপমাত্রা।

রেডিওমেট্রিক ক্রমাঙ্কন QCAL কে রূপান্তর করে এ-সেন্সর তেজ, একটি রেডিওমেট্রিক পরিমাণ
সমীকরণ ব্যবহার করে W/(m² * sr * µm) এ পরিমাপ করা হয়:

· লাভ = (Lmax - Lmin) / (QCALmax - QCALmin)

পক্ষপাত = Lmin - লাভ * QCALmin

দীপ্তি = লাভ * QCAL + পক্ষপাত
কোথায়, Lmax এবং লিমিন ক্রমাঙ্কন ধ্রুবক, এবং QCALmax এবং QCALmin হয়
QCAL-এ রিস্কেল করা রেডিয়েন্সের সীমার সর্বোচ্চ এবং সর্বনিম্ন পয়েন্ট।

তারপর, গণনা করতে এ-সেন্সর প্রতিফলন সমীকরণগুলি হল:

· সূর্যের আলো = [Esun * sin(e)] / (PI * d^2)

· প্রতিফলন = তেজ / সূর্যের আলো
কোথায়, d জ্যোতির্বিজ্ঞানের এককে পৃথিবী-সূর্য দূরত্ব, e সৌর উচ্চতা কোণ,
এবং এসুন W/(m² * µm) তে গড় সৌর বায়ুমণ্ডলীয় বিকিরণ।

সরলীকৃত এট-সারফেস মূল্যবোধ (পদ্ধতি=ডস[1-4])

বায়ুমণ্ডলীয় সংশোধন এবং প্রতিফলন ক্রমাঙ্কন পথের দীপ্তিকে সরিয়ে দেয়, অর্থাৎ
বায়ুমণ্ডল থেকে বিপথগামী আলো, এবং সৌর আলোকসজ্জার বর্ণালী প্রভাব। আউটপুট করতে
এই সহজ এট-সারফেস তেজ এবং এট-সারফেস প্রতিফলন, সমীকরণগুলি হল (এর জন্য নয়
তাপীয় ব্যান্ড):

· সূর্যের আলো = TAUv * [Esun * sin(e) * TAUz + Esky] / (PI * d^2)

· তেজ_পথ = তেজ_অন্ধকার - শতাংশ * সূর্যের আলো

· তেজ = (এট-সেন্সর_রেডিয়েন্স - দীপ্তি_পথ)

· প্রতিফলন = তেজ / সূর্যের আলো
কোথায়, শতাংশ 0.0 এবং 1.0 এর মধ্যে একটি মান (সাধারণত 0.01), এসকি ছড়িয়ে পড়া আকাশ
বিকিরণ TAUz সূর্য থেকে পথ বরাবর বায়ুমণ্ডলীয় প্রেরণ
স্থল পৃষ্ঠ, এবং TAUv ভূমি থেকে পথ বরাবর বায়ুমণ্ডলীয় প্রেরণ
সেন্সর পৃষ্ঠ. দীপ্তি_অন্ধকার অন্ধকার থেকে গণনা করা অ্যাট-সেন্সর উজ্জ্বলতা
অবজেক্ট, অর্থাৎ সম্পূর্ণ ইমেজের জন্য কমপক্ষে 'ডার্ক_প্যারামিটার' (সাধারণত 1000) পিক্সেল সহ DN।
মানগুলো হলো,

· DOS1: TAUv = 1.0, TAUz = 1.0 এবং Esky = 0.0

· DOS2: TAUv = 1.0, Esky = 0.0, এবং TAUz = sin(e) সর্বাধিক তরঙ্গ সহ সমস্ত ব্যান্ডের জন্য
দৈর্ঘ্য 1 এর কম। (যেমন ব্যান্ড 4-6 MSS, 1-4 TM, এবং 1-4 ETM+) অন্যান্য ব্যান্ড TAUz =
1.0

· DOS3: TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = rayleigh

· DOS4: TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = PI *
দীপ্তি_অন্ধকার
দৃষ্টি আকর্ষণ করছি: আউটপুট রেডিয়েন্স অস্পর্শ থাকে (অর্থাৎ নেতিবাচক হলে 0.0 সেট করা হয় না)
তারা সম্ভাব্য নেতিবাচক মান. যাইহোক, যখন আউটপুট প্রতিফলন 0.0 সেট করা হয়
একটি নেতিবাচক মান প্রাপ্ত.

নোট

আউটপুট রাস্টার সেল মানগুলি এর সাথে পুনরায় স্কেল করা যেতে পারে স্কেল প্যারামিটার (যেমন, 100 ইঞ্চি সহ
মধ্যে প্রতিফলিত আউটপুট ব্যবহার করার ক্ষেত্রে i.gensigset).

On ল্যান্ডস্যাট -8 মেটাডাটা ফাইল
NASA ল্যান্ডস্যাট ডেটার জন্য L1G মেটাডেটা ফাইলের (LDCM-DFCB-004.pdf) একটি কাঠামো রিপোর্ট করেছে
কন্টিনিউটি মিশন (অর্থাৎ ল্যান্ডস্যাট-৮)।

NASA ডিজিটাল রূপান্তর করার জন্য প্রয়োজনীয় তথ্য MIN_MAX_RADIANCE গ্রুপে ধরে রেখেছে৷
দীপ্তি মানগুলিতে সংখ্যা (DN)। তারপর, i.landsat.toar সম্ভাব্য মান প্রতিস্থাপন করে
মেটাডেটা মান সহ মান। ফলাফল মেটাডা দ্বারা রিপোর্ট করা মানের সাথে মেলে
RADIOMETRIC_RESCALING গ্রুপে ফাইল।

এছাড়াও, NASA সব ব্যান্ডের জন্য সর্বোচ্চ-মিনিট মান এবং
লাভ-পক্ষপাত মান। এটি অদ্ভুত যে সমস্ত ব্যান্ডের প্রতিফলনের একই পরিসর রয়েছে। এছাড়াও,
তারা সরাসরি ডিএন থেকে প্রতিফলন গণনা করার জন্য ওয়েব পৃষ্ঠায় লিখেছিল, প্রথমে
RADIOMETRIC_RESCALING মান এবং দ্বিতীয়টি পাপ (সূর্য_উচ্চতা) দ্বারা বিভক্ত।

এটি একটি সাধারণ রিস্কেলিং

· প্রতিফলন = তেজ / সূর্যের তেজ = (DN * RADIANCE_MULT + RADIANCE_ADD) /
সূর্যালোক

· এখন প্রতিফলন = DN * REFLECTANCE_MULT + REFLECTANCE_ADD

· তারপর REFLECTANCE_MULT = RADIANCE_MULT / সূর্যের তেজ

· এবং REFLECTANCE_ADD = RADIANCE_ADD / সূর্যের আলো

সমস্যা দেখা দেয় যখন সূর্যের আলো এবং গণনা করার জন্য আমাদের ESUN মান (প্রদান করা হয়নি) প্রয়োজন
ডস। আমরা অনুমান করি যে REFLECTANCE_MAXIMUM RADIANCE_MAXIMUM এর সাথে মিলে যায়, তারপর

· REFLECTANCE_MAXIMUM / sin(e) = RADIANCE_MAXIMUM / সূর্যালোক

· Esun = (PI * d^2) * RADIANCE_MAXIMUM / REFLECTANCE_MAXIMUM
কোথায় d মেটাডেটা ফাইল দ্বারা প্রদত্ত পৃথিবী-সূর্য দূরত্ব বা এর ভিতরে গণনা করা হয়
প্রোগ্রাম.

সার্জারির i.landsat.toar Lmax, Lmin, এবং Esun ব্যবহার করা চালিয়ে যেতে NASA পুনরায় স্কেলিং ফিরিয়ে দেয়
DN কে রেডিয়েন্সে এবং রেডিয়েন্সকে রিফ্ল্যাল্যান্সে রূপান্তর করতে ধ্রুবক গণনা করার মানগুলি
"প্রথাগত" সমীকরণ এবং সহজ বায়ুমণ্ডলীয় সংশোধন। দৃষ্টি আকর্ষণ করছি: যখন সর্বোচ্চ
মান প্রদান করা হয় না, i.landsat.toar থেকে Lmax, Lmin, এবং Esun গণনা করার চেষ্টা করে
RADIOMETRIC_RESCALING (পরীক্ষায় ফলাফল একই ছিল)।

ক্রমাঙ্কন ধ্রুবক
ভার্বোস মোডে (পতাকা -- ভারবোস), প্রোগ্রাম মৌলিক স্যাটেলাইট তথ্য লিখুন এবং
রূপান্তরে ব্যবহৃত পরামিতি।

উৎপাদন তারিখ একটি সঠিক মান নয় কিন্তু সঠিক ক্রমাঙ্কন প্রয়োগ করা প্রয়োজন
ধ্রুবক, যা তারিখে পরিবর্তিত হয়েছে:

· Landsat-1 MSS: কখনই না

· Landsat-2 MSS: জুলাই 16, 1975

· Landsat-3 MSS: জুন 1, 1978

· Landsat-4 MSS: 26 আগস্ট, 1982 এবং 1 এপ্রিল, 1983

· Landsat-4 TM: আগস্ট 1, 1983 এবং 15 জানুয়ারী, 1984

· Landsat-5 MSS: 6 এপ্রিল, 1984 এবং 9 নভেম্বর, 1984

· Landsat-5 TM: মে 4, 2003 এবং এপ্রিল, 2, 2007

· Landsat-7 ETM+: জুলাই 1, 2000

· Landsat-8 OLI/TIRS: 2013 সালে চালু হয়েছে

উদাহরণ

মেটাডাটা ফাইল উদাহরণ
ব্যান্ড রাস্টার 7_203, 30.1_203 [...]-এ Landsat-30.2 ETM+ এর ডিজিটাল নম্বরগুলিকে রূপান্তর করুন
আউটপুট ফাইল 203_30.1_toar, 203_30.2_toar [...] এবং
আউটপুট ফাইলগুলিতে সেন্সর তাপমাত্রা 293_39.61_toar এবং 293_39.62_toar:
i.landsat.toar ইনপুট=203_30। আউটপুট=_toar \
metfile=p203r030_7x20010620.met
or
i.landsat.toar input=L5121060_06020060714। \
output=L5121060_06020060714_toar \
metfile=L5121060_06020060714_MTL.txt
or
i.landsat.toar input=LC80160352013134LGN03_B output=toar \
metfile=LC80160352013134LGN03_MTL.txt sensor=oli8 date=2013-05-14

DOS1 উদাহরণ
DOS1 ব্যবহার করে প্রতিফলিত করার জন্য DN:
# চ্যানেলগুলির নাম পরিবর্তন করুন বা i.landsat.toar এর ইনপুট স্কিমের সাথে মিল করার জন্য একটি অনুলিপি তৈরি করুন:
g.copy raster=lsat7_2002_10,lsat7_2002.1
g.copy raster=lsat7_2002_20,lsat7_2002.2
g.copy raster=lsat7_2002_30,lsat7_2002.3
g.copy raster=lsat7_2002_40,lsat7_2002.4
g.copy raster=lsat7_2002_50,lsat7_2002.5
g.copy raster=lsat7_2002_61,lsat7_2002.61
g.copy raster=lsat7_2002_62,lsat7_2002.62
g.copy raster=lsat7_2002_70,lsat7_2002.7
g.copy raster=lsat7_2002_80,lsat7_2002.8
DOS1 ব্যবহার করে DN থেকে প্রতিফলিত মানগুলির গণনা (থেকে প্রাপ্ত মেটাডেটা
p016r035_7x20020524.met.gz):
i.landsat.toar ইনপুট=lsat7_2002। output=lsat7_2002_toar. সেন্সর = tm7 \
পদ্ধতি=dos1 তারিখ=2002-05-24 sun_elevation=64.7730999 \
পণ্য_তারিখ=2004-02-12 লাভ=HHHLHLHHL
ফলস্বরূপ ল্যান্ডস্যাট চ্যানেলগুলির নাম হল lsat7_2002_toar.1 .. lsat7_2002_toar.8.

তথ্যসূত্র

· চন্দর জি., বিএল মারখাম এবং ডিএল হেল্ডার, 2009: রিমোট সেন্সিং অফ এনভায়রনমেন্ট,
ফ্লাইট। 113

· চান্দের জিএইচ এবং বি. মার্কহাম, 2003.: জিওসায়েন্স অ্যান্ড রিমোট বিষয়ে IEEE লেনদেন
সেন্সিং, ভলিউম। 41, না। 11.

· শ্যাভেজ পিএস, জুনিয়র 1996. চিত্র-ভিত্তিক বায়ুমণ্ডলীয় সংশোধন - পুনর্বিবেচনা এবং
উন্নত। ফটোগ্রামমেট্রিক ইঞ্জিনিয়ারিং এবং রিমোট সেন্সিং 62(9): 1025-1036

হুয়াং এট আল: অ্যাট-স্যাটেলাইট প্রতিফলন, 2002: একটি প্রথম অর্ডার স্বাভাবিককরণ
Landsat 7 ETM+ ছবি।

· আর. আইরিশ: ল্যান্ডস্যাট 7. বিজ্ঞান ডেটা ব্যবহারকারীদের হ্যান্ডবুক। ফেব্রুয়ারি 17, 2007; 15 মে 2011।

মার্কহাম বিএল এবং জেএল বার্কার, 1986: ল্যান্ডস্যাট এমএসএস এবং টিএম পোস্ট-ক্যালিব্রেশন ডায়নামিক
রেঞ্জ, এক্সোআটমস্ফিয়ারিক রিফ্লেক্টেন্স এবং অ্যাট-স্যাটেলাইট তাপমাত্রা। EOSAT ল্যান্ডস্যাট
প্রযুক্তিগত নোট, নং 1।

· মোরান এমএস, আরডি জ্যাকসন, পিএন স্লেটার এবং পিএম টেইলেট, 1992: রিমোট সেন্সিং অফ
পরিবেশ, ভলিউম। 41.

· গান এট আল, 2001: ল্যান্ডস্যাট টিএম ডেটা ব্যবহার করে শ্রেণিবিন্যাস এবং পরিবর্তন সনাক্তকরণ, কখন
এবং কিভাবে বায়ুমণ্ডলীয় প্রভাব সংশোধন করতে? রিমোট সেন্সিং অফ এনভায়রনমেন্ট, ভলিউম। 75।

onworks.net পরিষেবা ব্যবহার করে অনলাইনে i.landsat.toargrass ব্যবহার করুন