Merge pull request #52 from nasa/ghg_tutorial

Add EMIT Methane Tutorial Content

.gitignore

@@ -1,14 +1,20 @@
+# python internals
 *.ipynb_checkpoints
+*__pycache__*
+
+# local development
 *.img*
 *.envi*
 *.hdr*
 *.nc*
+*.tif*
+*.pptx*
 /data/example_out.csv
 /data/isla_gaviota_3.geojson
 /data/outputs
 /data/envi
 /python/daac_data_download_python
 /python/emit_utils
 /ignore
-*__pycache__*
 /data/emit_asset_urls.txt
+/.vscode/

CHANGE_LOG.md

@@ -5,6 +5,21 @@ The format is based on [Keep a Changelog](http://keepachangelog.com/)
 and this project adheres to [Semantic Versioning](http://semver.org/).
 _________________________________________________________________________
 
+## 2024-03-13
+
+> ### Added
+>
+> - Added new `Visualizing Methane Plume Timeseries` notebook  
+> - Added new `Generating Methane Spectral Fingerprint` notebook
+> - Added new `tutorial_utils.py` module which has functions specific to the new CH4 notebooks  
+
+> ### Changed
+>
+> - Updated `setup_instructions.md` to provide more up to date Python setup instructions  
+> - Updated `README.md` to include new notebooks  
+> - Added some functions to `emit_tools.py` to support the new methane tutorial notebooks
+> - Minor updates to `How to Extract Area` and `How to Extract Points` notebooks to improve visualizations
+
 ## 2023-12-01
 
 > ### Changed

README.md

@@ -24,6 +24,8 @@ Below are the resources available for EMIT Data.
 |:---|:---|:---|
 |[Getting EMIT Data using EarthData Search](guides/Getting_EMIT_Data_using_EarthData_Search.md)|Markdown Guide|A thorough walkthrough for using [EarthData Search](https://search.earthdata.nasa.gov/search) to find and download EMIT data|
 |[Exploring EMIT L2A Reflectance](python/tutorials/Exploring_EMIT_L2A_Reflectance.ipynb)|Jupyter Notebook|Explore EMIT L2A Reflectance data using interactive plots|
+|[Visualizing Methane Plume Timeseries](python/tutorials/Visualizing_Methane_Plume_Timeseries.ipynb)|Jupyter Notebook|Find EMIT L2B CH4 Plume Data and build a timeseries of CH4 plume complexes|
+|[Generating_Methane_Spectral_Fingerprint](python/tutorials/Generating_Methane_Spectral_Fingerprint.ipynb)|Jupyter Notebook|Extract Radiance Spectra and build an in-plume/out-of-plume ratio to compare with CH4 absorption coefficient|
 |[How to find and access EMIT data](python/how-tos/How_to_find_and_access_EMIT_data.ipynb)|Jupyter Notebook|Use the `earthaccess` Python library to find and download or stream EMIT data|
 |[How to Convert to ENVI Format](python/how-tos/How_to_Convert_to_ENVI.ipynb)|Jupyter Notebook|Convert from downloaded netCDF4 (.nc) format to .envi format|
 |[How to Orthorectify](python/how-tos/How_to_Orthorectify.ipynb)|Jupyter Notebook|Use the geometry lookup table (GLT) included with the EMIT netCDF4 file to project on a geospatial grid (EPSG:4326)|
@@ -62,6 +64,6 @@ Email: <[email protected]>
 Voice: +1-866-573-3222  
 Organization: Land Processes Distributed Active Archive Center (LP DAAC)¹  
 Website: <https://lpdaac.usgs.gov/>  
-Date last modified: 07-07-2023  
+Date last modified: 03-13-2024  
 
 ¹Work performed under USGS contract G15PD00467 for NASA contract NNG14HH33I.  

data/methane_tutorial/emit20220815t042838_ch4_target

@@ -0,0 +1,285 @@
+001    381.006 0.000000000000000000
+002    388.409 0.000000000000000000
+003    395.816 0.000000000000000000
+004    403.225 0.000000000000000000
+005    410.638 0.000000000000000000
+006    418.054 0.000000000000000000
+007    425.472 0.000000000000000000
+008    432.893 0.000000000000000000
+009    440.317 0.000000000000000000
+010    447.743 0.000000000000000000
+011    455.170 0.000000000000000000
+012    462.599 0.000000000000000000
+013    470.030 0.000000000000000000
+014    477.463 0.000000000000000000
+015    484.897 0.000000000000000000
+016    492.333 0.000000000000000000
+017    499.771 0.000000000000000000
+018    507.210 0.000000000000000000
+019    514.650 0.000000000000000000
+020    522.091 0.000000000000000000
+021    529.533 0.000000000000000000
+022    536.977 0.000000000000000000
+023    544.421 0.000000000000000000
+024    551.867 0.000000000000000000
+025    559.314 0.000000000000000000
+026    566.762 0.000000000000000000
+027    574.209 0.000000000000000000
+028    581.659 0.000000000000000000
+029    589.108 0.000000000000000000
+030    596.558 0.000000000000000000
+031    604.010 0.000000000000000000
+032    611.462 0.000000000000000000
+033    618.915 0.000000000000000000
+034    626.368 0.000000000000000000
+035    633.821 0.000000000000000000
+036    641.276 0.000000000000000000
+037    648.730 0.000000000000000000
+038    656.186 0.000000000000000000
+039    663.641 0.000000000000000000
+040    671.098 0.000000000000000000
+041    678.554 0.000000000000000000
+042    686.010 0.000000000000000000
+043    693.468 0.000000000000000000
+044    700.925 0.000000000000000000
+045    708.384 0.000000000000000000
+046    715.841 0.000000000000000000
+047    723.299 0.000000000000000000
+048    730.759 0.000000000000000000
+049    738.217 0.000000000000000000
+050    745.677 0.000000000000000000
+051    753.136 0.000000000000000000
+052    760.596 0.000000000000000000
+053    768.056 0.000000000000000000
+054    775.516 0.000000000000000000
+055    782.978 0.000000000000000000
+056    790.438 0.000000000000000000
+057    797.899 0.000000000000000000
+058    805.362 0.000000000000000000
+059    812.823 0.000000000000000000
+060    820.285 0.000000000000000000
+061    827.746 0.000000000000000000
+062    835.207 0.000000000000000000
+063    842.670 0.000000000000000000
+064    850.131 0.000000000000000000
+065    857.594 0.000000000000000000
+066    865.055 0.000000000000000000
+067    872.518 0.000000000000000000
+068    879.980 0.000000000000000000
+069    887.441 0.000000000000000000
+070    894.904 0.000000000000000000
+071    902.366 0.000000000000000000
+072    909.829 0.000000000000000000
+073    917.291 0.000000000000000000
+074    924.754 0.000000000000000000
+075    932.216 0.000000000000000000
+076    939.679 0.000000000000000000
+077    947.140 0.000000000000000000
+078    954.603 0.000000000000000000
+079    962.064 0.000000000000000000
+080    969.527 0.000000000000000000
+081    976.988 0.000000000000000000
+082    984.450 0.000000000000000000
+083    991.911 0.000000000000000000
+084    999.373 0.000000000000000000
+085   1006.834 0.000000000000000000
+086   1014.295 0.000000000000000000
+087   1021.757 0.000000000000000000
+088   1029.217 0.000000000000000000
+089   1036.678 0.000000000000000000
+090   1044.138 0.000000000000000000
+091   1051.599 0.000000000000000000
+092   1059.060 0.000000000000000000
+093   1066.520 0.000000000000000000
+094   1073.980 0.000000000000000000
+095   1081.440 0.000000000000000000
+096   1088.900 0.000000000000000000
+097   1096.360 0.000000000000000000
+098   1103.818 0.000000000000000000
+099   1111.278 0.000000000000000000
+100   1118.737 0.000000000000000000
+101   1126.196 0.000000000000000000
+102   1133.655 0.000000000000000000
+103   1141.113 0.000000000000000000
+104   1148.572 0.000000000000000000
+105   1156.030 0.000000000000000000
+106   1163.488 0.000000000000000000
+107   1170.946 0.000000000000000000
+108   1178.404 0.000000000000000000
+109   1185.862 0.000000000000000000
+110   1193.318 0.000000000000000000
+111   1200.776 0.000000000000000000
+112   1208.233 0.000000000000000000
+113   1215.690 0.000000000000000000
+114   1223.147 0.000000000000000000
+115   1230.604 0.000000000000000000
+116   1238.060 0.000000000000000000
+117   1245.515 0.000000000000000000
+118   1252.972 0.000000000000000000
+119   1260.428 0.000000000000000000
+120   1267.883 0.000000000000000000
+121   1275.339 -0.000000069654603082
+122   1282.794 -0.000000111146726415
+123   1290.250 -0.000000182101968164
+124   1297.705 -0.000000305636779985
+125   1305.160 -0.000000527336328276
+126   1312.614 -0.000000929683246631
+127   1320.068 -0.000001671465390813
+128   1327.523 -0.000003054706867060
+129   1334.976 -0.000005628539263286
+130   1342.429 -0.000010447838948433
+131   1349.882 -0.000019460938889804
+132   1357.335 -0.000036627215903753
+133   1364.787 -0.000071207584885208
+134   1372.238 -0.000150330219280241
+135   1379.691 -0.000374439666851373
+136   1387.142 -0.001133920033820656
+137   1394.593 -0.003713146872892919
+138   1402.043 -0.010163969848653337
+139   1409.494 -0.017011214017695429
+140   1416.944 -0.017485140828103148
+141   1424.393 -0.020481224988834982
+142   1431.843 -0.015399980587926871
+143   1439.292 -0.007183837386920701
+144   1446.740 -0.005160594956011227
+145   1454.189 -0.003739930072627756
+146   1461.637 -0.002887781246559040
+147   1469.085 -0.003163939162842973
+148   1476.532 -0.003950452229267255
+149   1483.980 -0.002272845531566601
+150   1491.426 -0.001322976560280106
+151   1498.873 -0.000841631995893533
+152   1506.319 -0.000661007257311044
+153   1513.765 -0.000482083726915260
+154   1521.210 -0.000314814341620599
+155   1528.655 -0.000207755538147690
+156   1536.101 -0.000146404330602690
+157   1543.545 -0.000109917824796095
+158   1550.989 -0.000080034561866423
+159   1558.433 -0.000098030483352345
+160   1565.877 -0.000105905085002603
+161   1573.319 -0.000067902940535459
+162   1580.762 -0.000057276544744732
+163   1588.205 -0.000120224595071614
+164   1595.647 -0.000401908209860957
+165   1603.089 -0.001576972276221576
+166   1610.530 -0.006222410401732700
+167   1617.970 -0.021435322927493949
+168   1625.410 -0.059910292461777682
+169   1632.851 -0.112971337794280702
+170   1640.290 -0.120352695021060635
+171   1647.730 -0.129676205161983693
+172   1655.169 -0.106476159401298842
+173   1662.607 -0.144605482127502377
+174   1670.046 -0.165405726174070178
+175   1677.484 -0.080018266251608433
+176   1684.921 -0.091654209149709925
+177   1692.358 -0.116862482051208180
+178   1699.795 -0.159110003641637521
+179   1707.231 -0.188785289744263840
+180   1714.667 -0.095120403815047841
+181   1722.103 -0.121835208746717272
+182   1729.538 -0.155497119666077821
+183   1736.973 -0.104063918135378949
+184   1744.407 -0.071084409123687700
+185   1751.841 -0.050739589759742680
+186   1759.275 -0.047096044049794450
+187   1766.708 -0.064428685570723135
+188   1774.142 -0.093163445946300161
+189   1781.574 -0.060473895624796645
+190   1789.007 -0.091518120707544448
+191   1796.439 -0.070149779655907019
+192   1803.870 -0.028500757144960604
+193   1811.301 -0.031436803808985371
+194   1818.731 -0.028705900563194626
+195   1826.161 -0.002531317794645903
+196   1833.591 -0.000124534361707594
+197   1841.021 -0.000000106797917028
+198   1848.449 -0.000000000050185741
+199   1855.877 -0.000000013870771764
+200   1863.305 -0.000000103434050046
+201   1870.733 -0.000000050168210930
+202   1878.160 -0.000000105384939200
+203   1885.587 -0.000000275024896038
+204   1893.013 -0.000000404027025009
+205   1900.439 -0.000000186151182586
+206   1907.864 -0.000000113526786294
+207   1915.289 -0.000010222163086953
+208   1922.713 -0.000149251023607931
+209   1930.137 -0.000559261915934703
+210   1937.561 -0.001180059937729890
+211   1944.984 -0.001474258832851471
+212   1952.407 -0.000839866126592593
+213   1959.829 -0.000802114280937387
+214   1967.252 -0.001185584956468647
+215   1974.673 -0.000928373925971665
+216   1982.095 -0.000796093377431769
+217   1989.515 -0.000445132732861969
+218   1996.935 -0.000282454611375082
+219   2004.355 -0.000274540349301289
+220   2011.775 -0.000188231187288176
+221   2019.193 -0.000244587775762389
+222   2026.612 -0.000164175820921038
+223   2034.030 -0.000129539761708694
+224   2041.447 -0.000079308880014316
+225   2048.865 -0.000065495516761232
+226   2056.281 -0.000047772383600339
+227   2063.696 -0.000046003731905308
+228   2071.112 -0.000070168273049554
+229   2078.527 -0.000093521950338807
+230   2085.942 -0.000110860062235196
+231   2093.356 -0.000154363119462358
+232   2100.769 -0.000255881943941033
+233   2108.182 -0.000523277565558211
+234   2115.594 -0.001375186007124925
+235   2123.006 -0.003719877579447772
+236   2130.418 -0.009393931436509179
+237   2137.829 -0.020910727897321599
+238   2145.239 -0.039043649935313118
+239   2152.648 -0.068566906252670715
+240   2160.058 -0.097573402237012843
+241   2167.467 -0.113142164077031565
+242   2174.875 -0.099651068695073369
+243   2182.283 -0.067035094053992911
+244   2189.690 -0.053604723388108656
+245   2197.097 -0.329089281850986182
+246   2204.503 -0.456604549736337129
+247   2211.909 -0.142098157650880630
+248   2219.315 -0.221827819797885795
+249   2226.719 -0.322910943255505289
+250   2234.123 -0.473408806585640218
+251   2241.527 -0.572844899580387867
+252   2248.930 -0.668838327700018942
+253   2256.333 -0.730652529686325547
+254   2263.735 -0.757765945995275580
+255   2271.137 -0.669960748226233860
+256   2278.538 -0.683494057865527971
+257   2285.939 -0.804684649946742936
+258   2293.339 -1.070588085146407797
+259   2300.738 -1.161203428970881246
+260   2308.136 -0.608511216055738879
+261   2315.534 -0.621358606897385535
+262   2322.933 -0.933162539142736835
+263   2330.330 -0.736720258289056851
+264   2337.726 -1.091126161372567349
+265   2345.122 -1.453405372049445088
+266   2352.517 -1.231725418739371136
+267   2359.913 -0.657273677359269226
+268   2367.307 -0.770079200631794691
+269   2374.701 -1.199089582930701114
+270   2382.093 -0.681602805457481642
+271   2389.486 -0.674950280268280189
+272   2396.878 -0.662488986727583384
+273   2404.270 -0.513653895860329612
+274   2411.660 -0.298528928634253010
+275   2419.051 -0.303678622373404872
+276   2426.440 -0.348436616387659437
+277   2433.830 -0.212319449968018131
+278   2441.218 -0.166568524018598063
+279   2448.606 -0.154442898524042083
+280   2455.994 -0.114423875787858245
+281   2463.382 -0.077059434305238886
+282   2470.768 -0.060319040795876841
+283   2478.153 -0.037668219548523431
+284   2485.539 -0.029877843993233233
+285   2492.924 -0.034757379003344489

data/methane_tutorial/methane_inout_points.csv

@@ -0,0 +1,4 @@
+ID,in-plume,latitude,longitude
+0,1,39.4628,53.7743
+1,0,39.4697,53.7792
+2,0,39.4726,53.7797

data/methane_tutorial/rad_band_ratio_points.csv

@@ -0,0 +1,9 @@
+ID,in-plume,x,y
+0,0,508,421
+1,1,508,420
+2,1,508,419
+3,0,508,418
+4,0,508,417
+5,0,508,416
+6,0,508,415
+7,0,508,414

python/how-tos/How_to_Extract_Area.ipynb

@@ -177,7 +177,7 @@
    "outputs": [],
    "source": [
     "ds_850 = ds.sel(wavelengths=850,method='nearest')\n",
-    "ds_850.hvplot.image(cmap='viridis', frame_height=600, frame_width=600, geo=True).opts(title=f\"Reflectance at {ds_850.wavelengths.data:.3f} ({ds_850.wavelengths.units})\")*shape.hvplot(color='#d95f02', alpha=0.5)"
+    "ds_850.hvplot.image(cmap='viridis', frame_height=600, frame_width=600, geo=True, crs='EPSG:4326').opts(title=f\"Reflectance at {ds_850.wavelengths.data:.3f} ({ds_850.wavelengths.units})\")*shape.hvplot(color='#d95f02', alpha=0.5, geo=True, crs='EPSG:4326')"
    ]
   },
   {
@@ -256,7 +256,7 @@
     "Voice: +1-866-573-3222  \n",
     "Organization: Land Processes Distributed Active Archive Center (LP DAAC)¹  \n",
     "Website: <https://lpdaac.usgs.gov/>  \n",
-    "Date last modified: 06-30-2023  \n",
+    "Date last modified: 03-13-2024  \n",
     "\n",
     "¹Work performed under USGS contract G15PD00467 for NASA contract NNG14HH33I.  "
    ]
@@ -278,7 +278,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.4"
+   "version": "3.10.13"
   },
   "vscode": {
    "interpreter": {