Message Boards

GROUPS:

Staff Picks Data Science Curated Data Wolfram Language Geographic Information System Events & Media Wolfram Function Repository

Arnoud Buzing

Mapping Wildfires

Arnoud Buzing, Wolfram Research

Posted 2 years ago

5857 Views

2 Replies

7 Total Likes

Follow this post

This story originally appeared in Towards Data Science: Mapping Wildfires with the Wolfram Language By Arnoud Buzing Wildfires have been spreading in multiple states recently. California is suffering the most with some of its largest wildfires in history. Oregon has 20 large wildfires covering over 400,000 acres. And even Washington is not being spared this year. To get a sense of where the current wildfires are, I implemented a function called WildfireData which sources data from NASA’s MODIS instruments aboard its Terra and Aqua satellites (code included at the end of this post): https://modis.gsfc.nasa.gov/ This WildfireData function can be accessed through the following global ResourceFunction URI: WildfireData=ResourceFunction["user:arnoudb/DeployedResources/Function/WildfireData"] We can now immediately use this function to get the currently available wildfire data from NASA: data=WildfireData[] The output is a dataset with several columns. The first column holds the geographical location of a wildfire observation and the second column holds its brightness: We can start by plotting all the locations on a map: GeoGraphics[{Red,Point[data[All,"Location"]]}] The output shows the observed wildfires: At this zoom level, it looks like there are quite a few wildfires in the South, but we can zoom into, for example, Mississippi to see up close: ms=Entity["AdministrativeDivision",{"Mississippi","UnitedStates"}]; GeoGraphics[{Blue,Polygon[ms],Red,Point[data[All,"Location"]]},GeoRangems] Up close it becomes clear these wildfires are quite scattered and standalone: We can compare Mississippi to Oregon to see the difference: or=Entity["AdministrativeDivision",{"Oregon","UnitedStates"}]; GeoGraphics[{Blue,Polygon[or],Red,Point[data[All,"Location"]]},GeoRangeor] There are a far greater number of wildfires in Oregon: We can zoom in on any area, for example, the area within 40 miles of the center of Los Angeles: la=Entity["City",{"LosAngeles","California","UnitedStates"}]; GeoGraphics[{Red,Point[data[All,"Location"]]},GeoRangeGeoDisk[la,Quantity[40,"Miles"]]] There are many wildfires in the San Gabriel mountains: These wildfires have been named the Bobcat Fire (#bobcatfire on Twitter): There are also large wildfires east of Salem, Oregon: area=Entity["City",{"Salem","Oregon","UnitedStates"}] GeoGraphics[{Red,Point[data[All,"Location"]]},GeoRangeGeoDisk[area,Quantity[60,"Miles"]]] This shows a very big cluster of wildfires: Unfortunately, the wildfire season in the West is not expected to end until November. So far, in California alone, over 3,000,000 acres have burned with an estimated cost of $800,000,000. Oregon has lost over 600,000 acres and Washington has lost over 500,000 acres. Code for the WildfireData function: WildfireData[]:=Module[{base,csv1,csv2,csv,ds},base="https://firms.modaps.eosdis.nasa.gov/data/active_fire/c6/csv/";csv1=Once[Import[URLBuild[{base,"MODIS_C6_USA_contiguous_and_Hawaii_24h.csv"}]]];csv2=Once[Import[URLBuild[{base,"MODIS_C6_Alaska_24h.csv"}]]];csv=If[Length[csv2]>1,Join[csv1,csv2],csv1];ds=Dataset[Map[<\|"Location"GeoPosition[{First[#],#[[2]]}],"Brightness"#[[3]],"Scan"#[[4]],"Track"#[[5]],"AcquisitionDate"DateObject[{#[[6]]<>" "<>IntegerString[#[[7]],10,4],{"Year","-","Month","-","Day"," ","Hour","Minute"}}],"Satellite"#[[8]],"ConfidenceFraction"#[[9]]/100.,"Version"#[[10]],"BrightnessT31"#[[11]],"FRP"#[[12]],"DayNight"#[[13]]\|>&,Rest[csv]]]]

POSTED BY: Arnoud Buzing

2 Replies

Sort By:

Mads Bahrami

Mads Bahrami, Wolfram Research Inc

Posted 2 years ago

Very interesting post! A question about "AcquisitionDate": is there any info on what time zone the "acq_time" is collected? Maybe there is a convention for this sort of observations? If I evaluate WildfireData, depending on my location, I will get a different time, which may even implies a different day, compared to the actual measurement time/day. Right?

POSTED BY: Mads Bahrami

Moderation Team

Moderation Team, WOLFRAM

Posted 2 years ago

-- you have earned *Featured Contributor Badge* Your exceptional post has been selected for our editorial column *Staff Picks* http://wolfr.am/StaffPicks and Your Profile is now distinguished by a *Featured Contributor Badge* and is displayed on the Featured Contributor Board. Thank you!

POSTED BY: Moderation Team

Reply to this discussion

Reply Preview

Attachments

Remove Add a file to this post

Follow this discussion

or Discard

Group Abstract

Feedback

Enable JavaScript to interact with content and submit forms on Wolfram websites. Learn how »