code running

Question

code running

Opened this issue 5 years ago · 7 comments

Dear yannforget,

Could you please give three examples to how to run ledaps, lasrc and cfmask ater these packages are installed as conda-recipes show. I did not find the codes like do_ledaps.py etc.

Answer 1 · 2019-08-05T13:13:25.000Z

You can download do_ledaps.py from upstream repositories. It should work since ledaps and lasrc binaries are now installed in your conda environment.

The manual procedure is documented in upstream repositories (USGS-EROS/espa-surface-reflectance). I don't have the time to update the documentation, but IIRC the procedure was:

Dowload auxiliary data for ledaps or lasrc if needed.
Set required environment variables for ledaps or lasrc (LEDAPS_AUX_DIR, L8_AUX_DIR, ESUN, ESPA_LAND_MASS_POLYGON, ESPA_SCHEMA).

Then for ledaps:

convert_lpgs_to_espa --mtl <mtl_file>
create_angle_bands --xml <xml_file>
mask_per_pixel_angles.py --xml <xml_file>
lndpm --xml <xml_file> --process_sr=true
lndcal --pfile <lndcal_txt_file>
lndsr --pfile <lndsr_txt_file>
cfmask --xml <xml_file>
convert_espa_to_gtif --xml <xml_file> --gtif <arbitrary_prefix>

Or for lasrc:

convert_lpgs_to_espa --mtl <mtl_file>
create_angle_bands --xml <xml_file>
mask_per_pixel_angles.py --xml <xml_file>
lasrc --xml <xml_file> --aux <aux_file> --write-toa
cfmask --xml <xml_file>
convert_espa_to_gtif --xml <xml_file> --gtif

Answer 2 · 2019-08-05T16:33:58.000Z

Thank you for your reply. Can the Lasrc code be used for Sentinel-2 data? Could you please provide it if there are other git code for Sentinel-2? Yann Forget <notifications@github.com> 于2019年8月5日周一下午9:13写道：

…

You can download do_ledaps.py from upstream repositories. It should work since ledaps and lasrc binaries are now installed in your conda environment. The manual procedure is documented in upstream repositories ( USGS-EROS/espa-surface-reflectance <https://github.com/USGS-EROS/espa-surface-reflectance>). I don't have the time to update the documentation, but IIRC the procedure was: 1. Dowload auxiliary data for ledaps or lasrc if needed. 2. Set required environment variables for ledaps or lasrc ( LEDAPS_AUX_DIR, L8_AUX_DIR, ESUN, ESPA_LAND_MASS_POLYGON, ESPA_SCHEMA). Then for ledaps: convert_lpgs_to_espa --mtl <mtl_file> create_angle_bands --xml <xml_file> mask_per_pixel_angles.py --xml <xml_file> lndpm --xml <xml_file> --process_sr=true lndcal --pfile <lndcal_txt_file> lndsr --pfile <lndsr_txt_file> cfmask --xml <xml_file> convert_espa_to_gtif --xml <xml_file> --gtif <arbitrary_prefix> Or for lasrc: convert_lpgs_to_espa --mtl <mtl_file> create_angle_bands --xml <xml_file> mask_per_pixel_angles.py --xml <xml_file> lasrc --xml <xml_file> --aux <aux_file> --write-toa cfmask --xml <xml_file> convert_espa_to_gtif --xml <xml_file> --gtif — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#4?email_source=notifications&email_token=AH6R2LC62ULOEJF2LIIIPATQDARPNA5CNFSM4IHOKXV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD3RYWGY#issuecomment-518228763>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AH6R2LCPRBQBI2HFN4FZ5QDQDARPNANCNFSM4IHOKXVQ> .

-- ---------------------------------------------------------------------------------------------------------------------------------- Xianghong Che Associate professor of Government Geographic Information System Center, Chinese Academy of Surveying & Mapping (CASM) Email: chexh.15b@gmal.com Address: Lihuachi west Road, Haidian District, 100830, Beijing, China.

Answer 3 · 2019-08-21T10:02:17.000Z

Dear friend, What is the fuction of two commands below? I can get the sr data event if I do not run these two commands. create_angle_bands --xml <xml_file> mask_per_pixel_angles.py --xml <xml_file> Best Regards, Xanghong Che xianghong che <chexh.15b@gmail.com> 于2019年8月6日周二上午12:33写道：

…

Thank you for your reply. Can the Lasrc code be used for Sentinel-2 data? Could you please provide it if there are other git code for Sentinel-2? Yann Forget ***@***.***> 于2019年8月5日周一下午9:13写道： > You can download do_ledaps.py from upstream repositories. It should work > since ledaps and lasrc binaries are now installed in your conda > environment. > > The manual procedure is documented in upstream repositories ( > USGS-EROS/espa-surface-reflectance > <https://github.com/USGS-EROS/espa-surface-reflectance>). I don't have > the time to update the documentation, but IIRC the procedure was: > > 1. Dowload auxiliary data for ledaps or lasrc if needed. > 2. Set required environment variables for ledaps or lasrc ( > LEDAPS_AUX_DIR, L8_AUX_DIR, ESUN, ESPA_LAND_MASS_POLYGON, ESPA_SCHEMA > ). > > Then for ledaps: > > convert_lpgs_to_espa --mtl <mtl_file> > create_angle_bands --xml <xml_file> > mask_per_pixel_angles.py --xml <xml_file> > lndpm --xml <xml_file> --process_sr=true > lndcal --pfile <lndcal_txt_file> > lndsr --pfile <lndsr_txt_file> > cfmask --xml <xml_file> > convert_espa_to_gtif --xml <xml_file> --gtif <arbitrary_prefix> > > Or for lasrc: > > convert_lpgs_to_espa --mtl <mtl_file> > create_angle_bands --xml <xml_file> > mask_per_pixel_angles.py --xml <xml_file> > lasrc --xml <xml_file> --aux <aux_file> --write-toa > cfmask --xml <xml_file> > convert_espa_to_gtif --xml <xml_file> --gtif > > — > You are receiving this because you authored the thread. > Reply to this email directly, view it on GitHub > <#4?email_source=notifications&email_token=AH6R2LC62ULOEJF2LIIIPATQDARPNA5CNFSM4IHOKXV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD3RYWGY#issuecomment-518228763>, > or mute the thread > <https://github.com/notifications/unsubscribe-auth/AH6R2LCPRBQBI2HFN4FZ5QDQDARPNANCNFSM4IHOKXVQ> > . > -- ---------------------------------------------------------------------------------------------------------------------------------- Xianghong Che Associate professor of Government Geographic Information System Center, Chinese Academy of Surveying & Mapping (CASM) Email: ***@***.*** Address: Lihuachi west Road, Haidian District, 100830, Beijing, China.

-- ---------------------------------------------------------------------------------------------------------------------------------- Xianghong Che Associate professor of Government Geographic Information System Center, Chinese Academy of Surveying & Mapping (CASM) Email: chexh.15b@gmal.com Address: Lihuachi west Road, Haidian District, 100830, Beijing, China.

Answer 4 · 2019-08-29T09:44:56.000Z

Sorry for the late response,

create_angle_bands is a binary tool included in the recipe that generates per-pixel zenith and azimuth. For more details regarding each command in LEDAPS, please see the official documentation or the upstream do_ledaps.py script.

PS: No, LEDAPS/LaSRC code cannot be used with Sentinel data, only Landsat

Yann

Answer 5 · 2019-08-29T10:01:21.000Z

Dear Yann, I have read the official documentation, but this is no explanation for create_angle_bands and mask_per_pixel_angles.py. These two code are necessary for atmosphere correction? Yann Forget <notifications@github.com> 于2019年8月29日周四下午5:44写道：

…

Sorry for the late response, create_angle_bands is a binary tool included in the recipe that generates per-pixel zenith and azimuth. For more details regarding each command in LEDAPS, please see the official documentation <https://github.com/USGS-EROS/espa-surface-reflectance/blob/master/ledaps/ledapsSrc/doc/LEDAPS_User_Ref_SR.txt> or the upstream do_ledaps.py <https://github.com/USGS-EROS/espa-surface-reflectance/blob/master/ledaps/ledapsSrc/scripts/do_ledaps.py> script. PS: No, LEDAPS/LaSRC code cannot be used with Sentinel data, only Landsat Yann — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#4?email_source=notifications&email_token=AH6R2LC4JPYCT3ZVJVFOGWTQG6LBRA5CNFSM4IHOKXV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD5N5BUY#issuecomment-526110931>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AH6R2LCLS6LDXD2GELTVFETQG6LBRANCNFSM4IHOKXVQ> .

-- ---------------------------------------------------------------------------------------------------------------------------------- Xianghong Che Associate professor of Government Geographic Information System Center, Chinese Academy of Surveying & Mapping (CASM) Email: chexh.15b@gmal.com Address: Lihuachi west Road, Haidian District, 100830, Beijing, China.

Answer 6 · 2019-08-30T09:22:11.000Z

Yes, it is used by the LEDAPS modules. I'm not sure they are necessary though.

Note that, IIRC, create_angle_bands is for Landsat 8 data, and create_landsat_angle_bands is for Landsat TM & ETM.

Yann

Answer 7 · 2019-09-16T08:21:37.000Z

Dear Yannforget, Sorry to disturb you again. When using create_landsat_angle_bands, I had a issue which is showed below: Please do me a favor. Looking forward to hearing from you. [image: image.png] Yann Forget <notifications@github.com> 于2019年8月30日周五下午5:22写道：

…

Yes, it is used by the LEDAPS modules. I'm not sure they are necessary though. Note that, IIRC, create_angle_bands is for Landsat 8 data, and create_landsat_angle_bands is for Landsat TM & ETM. Yann — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#4?email_source=notifications&email_token=AH6R2LFA7EQFCHMDMTYXPO3QHDRELA5CNFSM4IHOKXV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD5RDOYA#issuecomment-526530400>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AH6R2LFHJ3MFXLXFMYHIFMTQHDRELANCNFSM4IHOKXVQ> .

-- ---------------------------------------------------------------------------------------------------------------------------------- Xianghong Che Associate professor of Government Geographic Information System Center, Chinese Academy of Surveying & Mapping (CASM) Email: chexh.15b@gmal.com Address: Lihuachi west Road, Haidian District, 100830, Beijing, China.